We study the imprints of boxy/peanut structures on the 2D line-of-sight kinematics of simulated disk galaxies. The models under study belong to a family with varying initial gas fraction and halo triaxiality, plus few other control runs with different structural parameters; the kinematic information was extracted using the Voronoi-binning technique and parametrised up to the fourth order of a Gauss-Hermite series. Building on a previous work for the long-slit case, we investigate the 2D kinematic behaviour in the edge-on projection as a function of the boxy/peanut strength and position angle; we find that for the strongest structures the highest moments show characteristic features away from the midplane in a range of position angles. We also discuss the masking effect of a classical bulge and the ambiguity in discriminating kinematically this spherically-symmetric component from a boxy/peanut bulge seen end-on. Regarding the face-on case, we extend existing results to encompass the effect of a second buckling and find that this phenomenon spurs an additional set of even deeper minima in the fourth moment. Finally, we show how the results evolve when inclining the disk away from perfectly edge-on and face-on. The behaviour of stars born during the course of the simulations is discussed and confronted to that of the pre-existing disk. The general aim of our study is providing a handle to identify boxy/peanut structure and their properties in latest generation IFU observations of nearby disk galaxies.

We examine the effects that the modelling of a Boxy/Peanut (B/P) bulge will have on the estimates of the stellar gravitational potential, forces, orbital structure and bar strength of barred galaxies. We present a method for obtaining the potential of disc galaxies from surface density images, assuming a vertical density distribution (height function), which is let to vary with position, thus enabling it to represent the geometry of a B/P. We construct a B/P height function after the results from a high-resolution, N-body+SPH simulation of an isolated galaxy and compare the resulting dynamical model to those obtained with the commonly used, position-independent "flat" height functions. We show that methods that do not allow for a B/P can induce errors in the forces in the bar region of up to 40% and demonstrate that this has a significant impact on the orbital structure of the model, which in turn determines its kinematics and morphology. Furthermore, we show that the bar strength is reduced in the presence of a B/P. We conclude that neglecting the vertical extent of a B/P can introduce considerable errors in the dynamical modelling. We also examine the errors introduced in the model due to uncertainties in the parameters of the B/P and show that even for generous but realistic values of the uncertainties, the error will be noticeably less than that of not modelling a B/P bulge at all.

We describe a NOVel form of Adaptive softening (NovA) for collisionless $N$-body simulations, implemented in the Ramses adaptive mesh refinement code. We introduce a refinement criterion that the particle distribution within each cell be sufficiently isotropic, as measured by its moment of inertia tensor. In this way, collapse is only refined if it occurs along all three axes, ensuring that the softening $\epsilon$ is always of order twice the largest inter-particle spacing in a cell. This more conservative force softening criterion is designed to minimise spurious two-body effects, while maintaining high force resolution in collapsed regions of the flow. We test NovA using an antisymmetric perturbed plane wave collapse (`Valinia' test) before applying it to warm dark matter (WDM) simulations. For the Valinia test, we show that -- unlike the standard $N$-body method -- NovA produces no numerical fragmentation while still being able to correctly capture fine caustics and shells around the collapsing regions. For the WDM simulations, we find that NovA converges significantly more rapidly than standard $N$-body, producing little or no spurious halos on small scales. We show, however, that determining whether or not halos exist below the free streaming mass $M_{\rm fs}$ is complicated by the fact that our halo finder (AHF) likely incorrectly labels some caustics and criss-crossing filaments as halos, while one or two particularly massive filaments appear to fragment in any version of NovA where refinement is allowed. Such massive filaments may be physically unstable to collapse, as is the case for infinite, static, self-gravitating cylinders. We will use NovA in forthcoming papers to study the issue of halo formation below $M_{\rm fs}$; filament stability; and to obtain new constraints on the temperature of dark matter.

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.

Dynamical simulations are a fundamental tool for studying the secular evolution of disc galaxies. Even at their maximum resolution, they still follow a limited number of particles and typically resolve scales of the order of a few tens of parsecs. Generally, the spatial resolution is defined by (some multiple of) the softening length, whose value is set as a compromise between the desired resolution and the need for limiting small-scale noise. Several works have studied the question whether a softening scale fixed in space and time provides a good enough modelling of an astrophysical system. Here we address this question within the context of dynamical simulations and disc instabilities. We first follow the evolution of a galaxy-like object in isolation and then set up a simulation of an idealised merger event. Alongside a run using the standard fixed-softening approach, we performed simulations where the softening lengths were let to vary from particle to particle according to the evolution of the local density field in space and time. Even though during the most violent phases of the merging the fixed-softening simulation tends to underestimate the resulting matter densities, as far as the evolution of the disc component is concerned we found no significant differences among the runs. We conclude that using an appropriate fixed softening scale is a safe approach to the problem of modelling an N-body, non-cosmological disc galaxy system.

Galaxies inhabiting a cluster environment experience significant evolution in their orbital motions throughout time; this is accompanied by changes in the anisotropy parameter, measuring the relative importance of radial and tangential motions for a given class of objects. Along with orbital changes, galaxies in clusters are well known to undergo severe alteration in their hot/cold gas content and star formation properties. Understanding the link between the changes in the internal properties of galaxies and their orbital motion is of crucial importance in the study of galaxy evolution, as it could unveil the primary mechanism responsible for its environmental dependence. Do the changes in the internal properties happen in parallel with those in the orbital motion? Or are the orbital features at the time of infall what determines the fate of the member galaxies? Alternatively: are the properties of galaxies at a given time related to the coeval orbital anisotropy or are they better related to the anisotropy at infall? In order to answer these questions, we studied the orbital evolution of different galaxy populations in the semi-analytic models of Guo et al. (2011) applied on to the Millennium Simulation. For each class of objects, characterised by different internal properties (such as age, star formation rate and colour), we studied the anisotropy profile at redshift zero and its evolution by tracing the progenitors back in time. We conclude that the orbital properties at infall strongly influence the subsequent evolution of the internal features of galaxies and that the overall anisotropy of the galaxy population tends to increase with time.

Cosmological simulations of structure formation follow the collisionless evolution of dark matter starting from a nearly homogeneous field at early times down to the highly clustered configuration at redshift zero. The density field is sampled by a number of particles in number infinitely smaller than those believed to be its actual components and this limits the mass and spatial scales over which we can trust the results of a simulation. Softening of the gravitational force is introduced in collisionless simulations to limit the importance of close encounters between these particles. The scale of softening is generally fixed and chosen as a compromise between the need for high spatial resolution and the need to limit the particle noise. In the scenario of cosmological simulations, where the density field evolves to a highly inhomogeneous state, this compromise results in an appropriate choice only for a certain class of objects, the others being subject to either a biased or a noisy dynamical description. We have implemented adaptive gravitational softening lengths in the cosmological simulation code GADGET; the formalism allows the softening scale to vary in space and time according to the density of the environment, at the price of modifying the equation of motion for the particles in order to be consistent with the new dependencies introduced in the system's Lagrangian. We have applied the technique to a number of test cases and to a set of cosmological simulations of structure formation. We conclude that the use of adaptive softening enhances the clustering of particles at small scales, a result visible in the amplitude of the correlation function and in the inner profile of massive objects, thereby anticipating the results expected from much higher resolution simulations.

We present the first numerical, N-body, hydrodynamical, chemical simulations of cosmic structure formation in the framework of non-Gaussian models. We study the impact of primordial non-Gaussianities on early chemistry (e, H, H+, H-, He, He+, He++, H2, H2+, D, D+, HD, HeH+), molecular and atomic gas cooling, star formation, metal (C, O, Si, Fe, Mg, S) enrichment, population III (popIII) and population II-I (popII) transition, and on the evolution of "visible" objects. We find that non-Gaussianities can have some consequences on baryonic structure formation at very early epochs, but the subsequent evolution at later times washes out any difference among the various models. When assuming reasonable values for primordial non-Gaussian perturbations, it turns out that they are responsible for: (i) altering early molecular fractions in the cold, dense gas phase of ~10 per cent; (ii) inducing small temperature fluctuations of <~10 per cent during the cosmic evolution of primordial objects; (iii) influencing the onset of the first star formation events, at z>~15, and of the popIII/popII transition of up to some 10^7yr; (iv) determining variations of <~10 per cent in the gas cloud and stellar mass distributions after the formation of the first structures; (v) causing only mild variations in the chemical history of the Universe. We stress, though, that purely non-Gaussian effects might be difficult to address, since they are strictly twisted with additional physical phenomena (e.g. primordial gas bulk flows, unknown primordial popIII stellar mass function, etc.) that have similar or stronger impact on the behaviour of the baryons.

[abridged] We present a detailed comparison of fundamental dark matter halo properties retrieved by a substantial number of different halo finders. These codes span a wide range of techniques including friends-of-friends (FOF), spherical-overdensity (SO) and phase-space based algorithms. We further introduce a robust (and publicly available) suite of test scenarios that allows halo finder developers to compare the performance of their codes against those presented here. This set includes mock haloes containing various levels and distributions of substructure at a range of resolutions as well as a cosmological simulation of the large-scale structure of the universe. All the halo finding codes tested could successfully recover the spatial location of our mock haloes. They further returned lists of particles (potentially) belonging to the object that led to coinciding values for the maximum of the circular velocity profile and the radius where it is reached. All the finders based in configuration space struggled to recover substructure that was located close to the centre of the host halo and the radial dependence of the mass recovered varies from finder to finder. Those finders based in phase space could resolve central substructure although they found difficulties in accurately recovering its properties. Via a resolution study we found that most of the finders could not reliably recover substructure containing fewer than 30-40 particles. However, also here the phase space finders excelled by resolving substructure down to 10-20 particles. By comparing the halo finders using a high resolution cosmological volume we found that they agree remarkably well on fundamental properties of astrophysical significance (e.g. mass, position, velocity, and peak of the rotation curve).

Several inflationary models predict the possibility that the primordial perturbations of the density field may contain a degree of non-Gaussianity which would influence the subsequent evolution of cosmic structures at large scales. In order to study their impact, we use a set of three cosmological DM-only simulations starting from initial conditions with different levels of non-Gaussianity: f_NL=0,+/-100. More specifically, we focus on the distribution of galaxy clusters at different redshifts and, using suitable scaling relations, we determine their X-ray and SZ signals. Our analysis allows us to estimate the differences in the logN-logS and logN-logY due to the different initial conditions and to predict the cluster counts at different redshifts expected for future surveys (eROSITA and SPT). We also use a second set of simulations assuming a different cosmological scenario to estimate how the dependence on f_NL is degenerate with respect to other parameters. Our results indicate that the effects introduced by a realistic amount of primordial non-Gaussianity are small when compared to the ones connected with current uncertainties in cosmological parameters, particularly with sigma_8. However, if future surveys will be associated with optical follow-up campaigns to determine the cluster redshift, an analysis of the samples at z>1 can provide significant constraints on f_NL. In particular we predict that the SPT cluster survey will be able to detect ~1000 clusters at z>1 for the Gaussian case, with a difference of 15-20 per cent associated to f_NL=+/-100.

The description of the abundance and clustering of halos for non-Gaussian initial conditions has recently received renewed interest, motivated by the forthcoming large galaxy and cluster surveys, which can potentially yield constraints of order unity on the non-Gaussianity parameter f_{NL}. We present tests on N-body simulations of analytical formulae describing the halo abundance and clustering for non-Gaussian initial conditions. We calibrate the analytic non-Gaussian mass function of Matarrese et al.(2000) and LoVerde et al.(2008) and the analytic description of clustering of halos for non-Gaussian initial conditions on N-body simulations. We find excellent agreement between the simulations and the analytic predictions if we make the corrections delta_c --> delta_c X sqrt{q} and delta_c --> \delta_c X q where q ~ 0.75, in the density threshold for gravitational collapse and in the non-Gaussian fractional correction to the halo bias, respectively. We discuss the implications of this correction on present and forecasted primordial non-Gaussianity constraints. We confirm that the non-Gaussian halo bias offers a robust and highly competitive test of primordial non-Gaussianity.