We discuss an extension of the massively parallel cosmological simulation code gadget-2, which now enables investigation of multiple and distinct gravitational force laws, provided they are dominated by a constant scaling of the Newtonian force. In addition to simplifying investigations of a universally modified force law, the ngravs extension provides a foundation for state-of-the-art collisionless cosmological simulations of exotic gravitational scenarios. We briefly review the algorithms used by gadget-2, and present our extension to multiple gravities, highlighting additional features that facilitate consideration of exotic force laws. We discuss the accuracy and performance of the ngravs extension, both internally and with an unaltered gadget-2, in the relevant operational modes. The ngravs extension is publicly released to the research community. Program summaryProgram title: Gadget-2.0.7-ngravsCatalogue identifier: AFAV_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AFAV_v1_0.htmlProgram obtainable from: CPC Program Library, Queen’s University, Belfast, N. IrelandLicensing provisions: GNU GPL v2No. of lines in distributed program, including test data, etc.: 36504No. of bytes in distributed program, including test data, etc.: 3919508Distribution format: tar.gzProgramming language: C.Computer: Commodity.Operating system: Unix.RAM: 256MB+Classification: 1.9, 4.12.External routines: GSL, FFTW3, HDF5 (optional)Nature of problem:N-body cosmological codes are traditionally designed to investigate a single gravitating species interacting via the Newtonian force law. There exist viable extensions to General Relativity [1], however, which predict weakfield, slow-motion limits featuring distinct gravitational force laws between distinct particle species. To enable investigation and constraint of these theories with available astrophysical data, a necessary first step is to extend an N-body simulator to handle distinct gravitating species.Solution method:The massively parallel Barnes–Hut tree, Fast Fourier Transform, and sorting routines of the versatile and well-vetted N-body simulator [2] gadget-2 were extended to support D distinct gravitationally interacting species. The tree implementation now vectorizes over each species’ monopole masses and positions, the Fourier routines now handle active and passive gravitational masses separately, and the sorting routines now group all particle data by type. An additional file was introduced allowing the user to specify all D2 gravitational interactions: real space, Fourier space, and lattice summation corrections. To improve monopole approximations in scenarios where the scale of the gravitational interaction depends on the mass itself, an optional tracking of the number of bodies contributing to any particular monopole approximation has been written.Restrictions:Force laws with mass-dependent scale lengths are not amenable to the implemented Fourier methods (or even the traditional [3] Fourier approach). Nodes containing highly heterogeneous collections of particles with different mass-dependent scale lengths may not be well-approximated, even with the additional tracking introduced. The collisional “gas” species can only interact via a single gravitational force law.Unusual features:The extension allows consideration of quite general and exotic interactions between bodies, and will serve as a common platform for model-dependent adjustments to the cosmological background evolution.Additional comments:Data file format is identical to that of gadget-2. Configuration file format is unchanged, save for the addition of required bindings between particle species and gravitational type.Running time:Typical running times are ≤2Dx those of gadget-2, where D is an integer between 1 and 6.