Collisionless N-body Simulations Research Articles

Dissipationless collapse and the dynamical mass–ellipticity relation of elliptical galaxies in Newtonian gravity and MOND

Context. Recent observational studies proposed an empirical relation between the dark-to-total mass ratio and ellipticity in elliptical galaxies based on their observed total dynamical mass-to-light ratio data M/L = (14.1 ± 5.4)ϵ. In other words, the larger the content of dark matter in the galaxy, the more the stellar component will be flattened. If true, this observation appears to be in stark contrast with the commonly accepted galaxy formation scenario, whereby this process takes place inside dark halos with reasonably spherical symmetry. Aims. Comparing the processes of dissipationless galaxy formation in different theories of gravity and the emergence of the galaxy scaling relations therein provides an important framework within which, in principle, one can discriminate between these processes. Methods. By means of collisionless N-body simulations in modified Newtonian dynamics (MOND) and Newtonian gravity with and without active dark matter halos, with both spherical and clumpy initial structure, I study the trends of intrinsic and projected ellipticities, Sérsic index, and anisotropy with the total dynamical-to-stellar mass ratio. Results. I show that the end products of both cold spherical collapses and mergers of smaller clumps show an increasing departure from spherical symmetry for increasing values of the total dynamical-to-stellar mass ratio, at least in a range of halo masses. The equivalent Newtonian systems of the end products of MOND collapses show a similar behaviour. However, the M/L relation obtained from the numerical experiments in both gravities is rather different from that reported by Deur and coauthors.

A mechanism of bar formation in disc galaxies: Synchronization of apsidal precession

ABSTRACT We discuss the mechanism(s) of bar formation in isolated and tidally interacting disc galaxies using the results of idealized collision-less N-body simulations of the galaxies. In order to better understand the mechanism, we investigate orbital eccentricities (e), epochs of apocentre passages (ta), azimuthal angles at ta (φa), precession rates (Ωpre), for individual stars, as well as bar strengths represented by relative m = 2 Fourier amplitude (A2) and bar pattern speeds (Ωbar). The main results are as follows. A significant fraction of stars with initially different φa and Ωpre in an isolated disc galaxy can have similar values within several dynamical time-scales. This synchronization of φa and Ωpre, which is referred to as apsidal precession synchronization (‘APS’) in this study, is caused by the enhanced strength of the tangential component of gravitational force. A weak seed bar (A2 < 0.1) is first formed through APS in local regions of a disc, then the bar grows due to APS. In the bar growth phase (0.1 < A2 < 0.4), APS can proceed more efficiently due to stronger tangential force from the bar so that it can enhance the bar strength further. This positive feedback loop in APS is the key physical mechanism of bar growth in isolated stellar discs. Bar formation can be severely suppressed in discs with lower disc mass fractions and/or higher Q parameters due to much less efficient APS. APS proceeds more rapidly and more efficiently due to strong tidal perturbation in the formation of tidal bars compared to spontaneous bar formation.

Migration and heating in the galactic disc from encounters between Sagittarius and the Milky Way

ABSTRACT Stars born on near-circular orbits in spiral galaxies can subsequently migrate to different orbits due to interactions with non-axisymmetric disturbances within the disc such as bars or spiral arms. This paper extends the study of migration to examine the role of external influences using the example of the interaction of the Sagittarius dwarf galaxy (Sgr) with the Milky Way (MW). We first make impulse approximation estimates to characterize the influence of Sgr disc passages. The tidal forcing from Sgr can produce changes in both guiding radius ΔRg and orbital eccentricity, as quantified by the maximum radial excursion ΔRmax. These changes follow a quadrupole-like pattern across the face of the disc, with amplitude increasing with Galactocentric radius. We next examine a collisionless N-body simulation of a Sgr-like satellite interacting with an MW-like galaxy and find that Sgr’s influence in the outer disc dominates the secular evolution of orbits between disc passages. Finally, we use the same simulation to explore possible observable signatures of Sgr-induced migration by painting the simulation with different age stellar populations. We find that following Sgr disc passages, the migration it induces manifests within an annulus as an approximate quadrupole in azimuthal metallicity variations (δ[Fe/H]), along with systematic variations in orbital eccentricity, ΔRmax. These systematic variations can persist for several rotational periods. We conclude that this combination of signatures may be used to distinguish between the different migration mechanisms shaping the chemical abundance patterns of the MW’s thin disc.

Bar-driven leading spiral arms in a counter-rotating dark matter halo

ABSTRACT An overwhelming majority of galactic spiral arms trail with respect to the rotation of the galaxy, though a small sample of leading spiral arms has been observed. The formation of these leading spirals is not well understood. Here we show, using collisionless N-body simulations, that a barred disc galaxy in a retrograde dark matter halo can produce long-lived (∼3 Gyr) leading spiral arms. Due to the strong resonant coupling of the disc to the halo, the bar slows rapidly and spiral perturbations are forced ahead of the bar. We predict that such a system, if observed, will also host a dark matter wake oriented perpendicular to the stellar bar. More generally, we propose that any mechanism that rapidly decelerates the stellar bar will allow leading spiral arms to flourish.

Effective N-body models of composite collisionless stellar systems

ABSTRACT Gas-poor galaxies can be modelled as composite collisionless stellar systems, with a dark matter halo and one or more stellar components, representing different stellar populations. The dynamical evolution of such composite systems is often studied with numerical N-body simulations, whose initial conditions typically require realizations with particles of stationary galaxy models. We present a novel method to conceive these N-body realizations, which allows one to exploit at best a collisionless N-body simulation that follows their evolution. The method is based on the use of an effective N-body model of a composite system, which is in fact realized as a one-component system of particles that is interpreted a posteriori as a multicomponent system, by assigning in post-processing fractions of each particle’s mass to different components. Examples of astrophysical applications are N-body simulations that aim to reproduce the observed properties of interacting galaxies, satellite galaxies, and stellar streams. As a case study we apply our method to an N-body simulation of tidal stripping of a two-component (dark matter and stars) satellite dwarf galaxy orbiting in the gravitational potential of the Milky Way.

Higher order initial conditions for mixed baryon–CDM simulations

ABSTRACT We present a novel approach to generate higher order initial conditions (ICs) for cosmological simulations that take into account the distinct evolution of baryons and dark matter. We focus on the numerical implementation and the validation of its performance, based on both collisionless N-body simulations and full hydrodynamic Eulerian and Lagrangian simulations. We improve in various ways over previous approaches that were limited to first-order Lagrangian perturbation theory (LPT). Specifically, we (1) generalize nth-order LPT to multifluid systems, allowing 2LPT or 3LPT ICs for two-fluid simulations, (2) employ a novel propagator perturbation theory to set up ICs for Eulerian codes that are fully consistent with 1LPT or 2LPT, (3) demonstrate that our ICs resolve previous problems of two-fluid simulations by using variations in particle masses that eliminate spurious deviations from expected perturbative results, (4) show that the improvements achieved by going to higher order PT are comparable to those seen for single-fluid ICs, and (5) demonstrate the excellent (i.e. few per cent level) agreement between Eulerian and Lagrangian simulations, once high-quality initial conditions are used. The rigorous development of the underlying perturbation theory is presented in a companion paper. All presented algorithms are implemented in the monofonic music-2 package that we make publicly available.

Double X/Peanut structures in barred galaxies – insights from an N-body simulation

ABSTRACT Boxy, peanut-, or X-shaped ‘bulges’ are observed in a large fraction of barred galaxies viewed in, or close to, edge-on projection, as well as in the Milky Way. They are the product of dynamical instabilities occurring in stellar bars, which cause the latter to buckle and thicken vertically. Recent studies have found nearby galaxies that harbour two such features arising at different radial scales, in a nested configuration. In this paper, we explore the formation of such double peanuts, using a collisionless N-body simulation of a pure disc evolving in isolation within a live dark matter halo, which we analyse in a completely analogous way to observations of real galaxies. In the simulation, we find a stable double configuration consisting of two X/peanut structures associated with the same galactic bar – rotating with the same pattern speed – but with different morphology, formation time, and evolution. The inner, conventional peanut-shaped structure forms early via the buckling of the bar, and experiences little evolution once it stabilizes. This feature is consistent in terms of size, strength, and morphology, with peanut structures observed in nearby galaxies. The outer structure, however, displays a strong X, or ‘bow-tie’, morphology. It forms just after the inner peanut, and gradually extends in time (within 1–1.5 Gyr) to almost the end of the bar, a radial scale where ansae occur. We conclude that, although both structures form, and are dynamically coupled to, the same bar, they are supported by inherently different mechanisms.

On the ridges, undulations & streams in Gaia DR2: Linking the topography of phase-space to the orbital structure of an N-body bar

Abstract We explore the origin of phase-space substructures revealed by the second Gaia data release in the disc of the Milky Way, such as the ridges in the Vφ-r plane, the undulations in the Vφ-r-Vr space and the streams in the Vφ-Vr plane. We use a collisionless N-body simulation with co-spatial thin and thick discs, along with orbit integration, to study the orbital structure close to the Outer Lindblad Resonance (OLR) of the bar. We find that a prominent, long-lived ridge is formed in the Vφ-r plane due to the OLR which translates to streams in the Vφ-Vr plane and examine which closed periodic and trapped librating orbits are responsible for these features. We find that orbits which carry out small librations around the x1(1) family are preferentially found at negative Vr, giving rise to a ‘horn’-like feature, while orbits with larger libration amplitudes, trapped around the x1(2) and x1(1) families, constitute the positive Vr substructure, i.e. the Hercules-like feature. This changing libration amplitude of orbits will translate to a changing ratio of thin/thick disc stars, which could have implications on the metallicity distribution in this plane. We find that a scenario in which the Sun is placed close to the OLR gives rise to a strong asymmetry in Vr in the Vφ-Vr plane (i.e. Hercules vs. ‘the horn’) and subsequently to undulations in the Vφ-r-Vr space. We also explore a scenario in which the Sun is placed closer to the bar corotation and find that the bar perturbation alone cannot give rise to the these features.

Formation of LISA Black Hole Binaries in Merging Dwarf Galaxies: The Imprint of Dark Matter

Abstract Theoretical models for the expected merger rates of intermediate-mass black holes (IMBHs) are vital for planned gravitational-wave detection experiments such as the Laser Interferometer Space Antenna (LISA). Using collisionless N-body simulations of dwarf galaxy (DG) mergers, we examine how the orbital decay of IMBHs and the efficiency of IMBH binary formation depend on the central dark matter (DM) density profile of the merging DGs. Specifically, we explore various asymptotic inner slopes γ of the DG’s DM density distribution, ranging from steep cusps (γ = 1) to shallower density profiles (γ < 1), motivated by well-known baryonic-feedback effects as well as by DM models that differ from cold DM at the scales of DGs. We find that the inner DM slope is crucial for the formation (or lack thereof) of an IMBH binary; only mergers between DGs with cuspy DM profiles (γ = 1) are favorable to forming a hard IMBH binary, whereas when γ < 1 the IMBHs stall at a separation of 50–100 pc. Consequently, the rate of LISA signals from IMBH coalescence will be determined by the fraction of DGs with a cuspy DM profile. Conversely, the LISA event rates at IMBH mass scales offer in principle a novel way to place constraints on the inner structure of DM halos in DGs and address the core–cusp controversy. We also show that, with spatial resolutions of ∼0.1 kpc, as often adopted in cosmological simulations, all IMBHs stall, independent of γ. This suggests caution should be taken when employing cosmological simulations of galaxy formation to study BH dynamics in DGs.

Formation and structure of ultralight bosonic dark matter halos

We simulate the formation and evolution of ultralight bosonic dark matter halos from cosmological initial conditions. Using zoom-in techniques we are able to resolve the detailed interior structure of the halos. We observe the formation of solitonic cores and confirm the core-halo mass relation previously found by Schive et al. The cores exhibit strong quasi-normal oscillations that remain largely undamped on evolutionary timescales. On the other hand, no conclusive growth of the core mass by condensation or relaxation can be detected. In the incoherent halo surrounding the cores, the scalar field density profiles and velocity distributions show no significant deviation from collisionless N-body simulations on scales larger than the coherence length. Our results are consistent with the core properties being determined mainly by the coherence length at the time of virialization, whereas the Schr\"odinger-Vlasov correspondence explains the halo properties when averaged on scales greater than the coherence length.

Around the Way: Testing ΛCDM with Milky Way Stellar Stream Constraints

Abstract Recent analyses of the Pal 5 and GD-1 tidal streams suggest that the inner dark matter halo of the Milky Way is close to spherical, in tension with predictions from collisionless N-body simulations of cosmological structure formation. We use the Eris simulation to test whether the combination of dissipative physics and hierarchical structure formation can produce Milky Way–like galaxies whose dark matter halos match the tidal stream constraints from the GD-1 and Pal 5 clusters. We use a dynamical model of the simulated Eris galaxy to generate many realizations of the GD-1 and Pal 5 tidal streams, marginalize over observational uncertainties in the cluster galactocentric positions and velocities, and compare with the observational constraints. We find that the total density and potential of Eris contributed by baryons and dark matter satisfies constraints from the existing Milky Way stellar stream data, as the baryons both round and redistribute the dark matter during the dissipative formation of the galaxy, and provide a centrally concentrated mass distribution that rounds the inner potential. The Eris dark matter halo or a spherical Navarro–Frenk–White dark matter work comparably well in modeling the stream data. In contrast, the equivalent dark matter–only ErisDark simulation produces a prolate halo that cannot reproduce the observed stream data. The ongoing Gaia mission will provide decisive tests of the consistency between and Milky Way streams, and should distinguish between models like Eris and more spherical halos.

Dark-ages Reionization and Galaxy Formation Simulation – XIV. Gas accretion, cooling, and star formation in dwarf galaxies at high redshift

We study dwarf galaxy formation at high redshift ($z\ge5$) using a suite of high- resolution, cosmological hydrodynamic simulations and a semi-analytic model (SAM). We focus on gas accretion, cooling and star formation in this work by isolating the relevant process from reionization and supernova feedback, which will be further discussed in a companion paper. We apply the SAM to halo merger trees constructed from a collisionless N-body simulation sharing identical initial conditions to the hydrodynamic suite, and calibrate the free parameters against the stellar mass function predicted by the hydrodynamic simulations at z = 5. By making comparisons of the star formation history and gas components calculated by the two modelling techniques, we find that semi-analytic prescriptions that are commonly adopted in the literature of low-redshift galaxy formation do not accurately represent dwarf galaxy properties in the hydrodynamic simulation at earlier times. We propose 3 modifications to SAMs that will provide more accurate high-redshift simulations. These include 1) the halo mass and baryon fraction which are overestimated by collisionless N-body simulations; 2) the star formation efficiency which follows a different cosmic evolutionary path from the hydrodynamic simulation; and 3) the cooling rate which is not well defined for dwarf galaxies at high redshift. Accurate semi-analytic modelling of dwarf galaxy formation informed by detailed hydrodynamical modelling will facilitate reliable semi-analytic predictions over the large volumes needed for the study of reionization.

On the cavity of a debris disc carved by a giant planet

One possible explanation of the cavity in debris discs is the gravitational perturbation of an embedded giant planet. Planetesimals passing close to a massive body are dynamically stirred resulting in a cleared region known as the chaotic zone. Theory of overlapping mean-motion resonances predicts the width of this cavity. To test whether this cavity is identical to the chaotic zone, we investigate the formation of cavities by means of collisionless N-body simulations assuming a 1.25-10 Jupiter mass planet with eccentricities of 0-0.9. Synthetic images at millimetre wavelengths are calculated to determine the cavity properties by fitting an ellipse to 14 percent contour level. Depending on the planetary eccentricity, e_pl, the elliptic cavity wall rotates as the planet orbits with the same (e_pl<0.2) or half (e_pl>0.2) period that of the planet. The cavity centre is offset from the star along the semi-major axis of the planet with a distance of d=0.1q^-0.17e_pl^0.5 in units of cavity size towards the planet's orbital apocentre, where q is the planet-to-star mass ratio. Pericentre (apocentre) glow develops for e_pl<0.05 (e_pl>0.1), while both are present for 0.05<=e_pl<=0.1. Empirical formulae are derived for the sizes of the cavities: da_cav=2.35q^0.36 and da_cav=7.87q^0.37e_pl^0.38 for e_pl<=0.05 and e_pl>0.05, respectively. The cavity eccentricity, e_cav, equals to that of the planet only for 0.3<=e_pl<=0.6. A new method based on ALMA observations for estimating the orbital parameters and mass of the planet carving the cavity is also given.

Deadly Dark Matter Cusps versus Faint and Extended Star Clusters: Eridanus II and Andromeda XXV

Abstract The recent detection of two faint and extended star clusters in the central regions of two Local Group dwarf galaxies, Eridanus II and Andromeda XXV, raises the question of whether clusters with such low densities can survive the tidal field of cold dark matter halos with central density cusps. Using both analytic arguments and a suite of collisionless N-body simulations, I show that these clusters are extremely fragile and quickly disrupted in the presence of central cusps with . Furthermore, the scenario in which the clusters were originally more massive and sank to the center of the halo requires extreme fine tuning and does not naturally reproduce the observed systems. In turn, these clusters are long lived in cored halos, whose central regions are safe shelters for . The only viable scenario for hosts that have preserved their primordial cusp to the present time is that the clusters formed at rest at the bottom of the potential, which is easily tested by measurement of the clusters proper velocity within the host. This offers means to readily probe the central density profile of two dwarf galaxies as faint as and , in which stellar feedback is unlikely to be effective.

Angular momentum properties of haloes and their baryon content in the Illustris simulation

The angular momentum properties of virialised dark matter haloes have been measured with good statistics in collisionless N-body simulations, but an equally accurate analysis of the baryonic spin is still missing. We employ the Illustris simulation suite, one of the first simulations of galaxy formation with full hydrodynamics that produces a realistic galaxy population in a sizeable volume, to quantify the baryonic spin properties for more than $\sim$ 320,000 haloes. We first compare the systematic differences between different spin parameter and halo definitions, and the impact of sample selection criteria on the derived properties. We confirm that dark matter only haloes exhibit a close to self-similar spin distribution in mass and redshift of lognormal form. However, the physics of galaxy formation radically changes the baryonic spin distribution. While the dark matter component remains largely unaffected, strong trends with mass and redshift appear for the spin of diffuse gas and the formed stellar component. With time the baryons staying bound to the halo develop a misalignment of their spin vector with respect to dark matter, and increase their specific angular momentum by a factor of $\sim$ 1.3 in the non-radiative case and $\sim$ 1.8 in the full physics setup at z = 0. We show that this enhancement in baryonic spin can be explained by the combined effect of specific angular momentum transfer from dark matter onto gas during mergers and from feedback expelling low specific angular momentum gas from the halo. Our results challenge certain models for spin evolution and underline the significant changes induced by baryonic physics in the structure of haloes.

GOTHIC: Gravitational oct-tree code accelerated by hierarchical time step controlling

The tree method is a widely implemented algorithm for collisionless N-body simulations in astrophysics well suited for GPU(s). Adopting hierarchical time stepping can accelerate N-body simulations; however, it is infrequently implemented and its potential remains untested in GPU implementations. We have developed a Gravitational Oct-Tree code accelerated by HIerarchical time step Controlling named GOTHIC, which adopts both the tree method and the hierarchical time step. The code adopts some adaptive optimizations by monitoring the execution time of each function on-the-fly and minimizes the time-to-solution by balancing the measured time of multiple functions. Results of performance measurements with realistic particle distribution performed on NVIDIA Tesla M2090, K20X, and GeForce GTX TITAN X, which are representative GPUs of the Fermi, Kepler, and Maxwell generation of GPUs, show that the hierarchical time step achieves a speedup by a factor of around 3–5 times compared to the shared time step. The measured elapsed time per step of GOTHIC is 0.30 s or 0.44 s on GTX TITAN X when the particle distribution represents the Andromeda galaxy or the NFW sphere, respectively, with 224= 16,777,216 particles. The averaged performance of the code corresponds to 10–30% of the theoretical single precision peak performance of the GPU.

DAMPING OF THE MILKY WAY BAR BY MANIFOLD-DRIVEN SPIRALS

ABSTRACT We describe a new phenomenon of “bar damping” that may have played an important role in shaping the Milky Way bar and bulge as well as its spiral structure. We use a collisionless N-body simulation of a Milky Way–like galaxy initially composed of a dark matter halo and an exponential disk with a Toomre parameter slightly above unity. In this configuration, dominated by the disk in the center, a bar forms relatively quickly, after 1 Gyr of evolution. This is immediately followed by the formation of two manifold-driven spiral arms and the outflow of stars that modifies the potential in the vicinity of the bar, apparently shifting the position of the L 1/L 2 Lagrange points. This modification leads to the shortening of the bar and the creation of a next generation of manifold-driven spiral arms at a smaller radius. The process repeats itself a few times over the next 0.5 Gyr resulting in further substantial weakening and shortening of the bar. The time when the damping comes to an end coincides with the first buckling episode in the bar that rebuilds the orbital structure so that no more new spiral arms are formed. The morphology of the bar and the spiral structure at this time show remarkable similarity to the present properties of the Milky Way. Later on, the bar starts to grow rather steadily again, weakened only by subsequent buckling episodes occurring at more distant parts of the disk.

General relativistic screening in cosmological simulations

We revisit the issue of interpreting the results of large volume cosmological simulations in the context of large scale general relativistic effects. We look for simple modifications to the nonlinear evolution of the gravitational potential $\psi$ that lead on large scales to the correct, fully relativistic description of density perturbations in the Newtonian gauge. We note that the relativistic constraint equation for $\psi$ can be cast as a diffusion equation, with a diffusion length scale determined by the expansion of the Universe. Exploiting the weak time evolution of $\psi$ in all regimes of interest, this equation can be further accurately approximated as a Helmholtz equation, with an effective relativistic 'screening' scale $\ell$ related to the Hubble radius. We demonstrate that it is thus possible to carry out N-body simulations in the Newtonian gauge by replacing Poisson's equation with this Helmholtz equation, involving a trivial change in the Green's function kernel. Our results also motivate a simple, approximate (but very accurate) gauge transformation - $\delta_{\rm N}(\mathbf{k}) \approx \delta_{\rm sim}(\mathbf{k})\times (k^2+\ell^{-2})/k^2$ - to convert the density field $\delta_{\rm sim}$ of standard collisionless N-body simulations (initialised in the comoving synchronous gauge) into the Newtonian gauge density $\delta_{\rm N}$ at arbitrary times. A similar conversion can also be written in terms of particle positions. Our results can be interpreted in terms of a Jeans stability criterion induced by the expansion of the Universe. The appearance of the screening scale $\ell$ in the evolution of $\psi$, in particular, leads to a natural resolution of the 'Jeans swindle' in the presence of super-horizon modes.

Contributions to the accreted stellar halo: an atlas of stellar deposition

The accreted component of stellar halos is composed of the contributions of several satellites, falling onto their host with their different masses, at different times, on different orbits. This work uses a suite of idealised, collisionless N-body simulations of minor mergers and a particle tagging technique to understand how these different ingredients shape each contribution to the accreted halo, in both density and kinematics. I find that more massive satellites deposit their stars deeper into the gravitational potential of the host, with a clear segregation enforced by dynamical friction. Earlier accretion events contribute more to the inner regions of the halo; more concentrated subhaloes sink deeper through increased dynamical friction. The orbital circularity of the progenitor at infall is only important for low-mass satellites: dynamical friction efficiently radialises the most massive minor mergers erasing the imprint of the infall orbit for satellite-to-host virial mass ratios $\gtrsim1/20$. The kinematics of the stars contributed by each satellite is also ordered with satellite mass: low-mass satellites contribute fast-moving populations, in both ordered rotation and radial velocity dispersion. In turn, contributions by massive satellites have lower velocity dispersion and lose their angular momentum to dynamical friction, resulting in a strong radial anisotropy.

MODELING THE Lyα FOREST IN COLLISIONLESS SIMULATIONS

ABSTRACT Cosmological hydrodynamic simulations can accurately predict the properties of the intergalactic medium (IGM), but only under the condition of retaining the high spatial resolution necessary to resolve density fluctuations in the IGM. This resolution constraint prohibits simulating large volumes, such as those probed by BOSS and future surveys, like DESI and 4MOST. To overcome this limitation, we present “Iteratively Matched Statistics” (IMS), a novel method to accurately model the Lyα forest with collisionless N-body simulations, where the relevant density fluctuations are unresolved. We use a small-box, high-resolution hydrodynamic simulation to obtain the probability distribution function (PDF) and the power spectrum of the real-space Lyα forest flux. These two statistics are iteratively mapped onto a pseudo-flux field of an N-body simulation, which we construct from the matter density. We demonstrate that our method can reproduce the PDF, line of sight and 3D power spectra of the Lyα forest with good accuracy (7%, 4%, and 7% respectively). We quantify the performance of the commonly used Gaussian smoothing technique and show that it has significantly lower accuracy (20%–80%), especially for N-body simulations with achievable mean inter-particle separations in large-volume simulations. In addition, we show that IMS produces reasonable and smooth spectra, making it a powerful tool for modeling the IGM in large cosmological volumes and for producing realistic “mock” skies for Lyα forest surveys.