The ability to accurately predict the evolution of tidally stripped haloes is important for understanding galaxy formation and testing the properties of dark matter. Most studies of substructure evolution make predictions based on empirical models of tidal mass loss that are calibrated using numerical simulations. This approach can be accurate in the cases considered, but lacks generality and does not provide a physical understanding of the processes involved. Recently, we demonstrated that truncating NFW distribution functions sharply in energy results in density profiles that resemble those of tidally stripped systems, offering a path to constructing physically motivated models of tidal mass loss. In this work, we review calculations of mass loss based on energy truncation alone, and then consider what secondary effects may modulate mass loss beyond this. We find that a combination of dependence on additional orbital parameters and variations in individual particle energies over an orbit results in a less abrupt truncation in energy space as a subhalo loses mass. Combining the energy truncation approach with a simple prediction for the mass-loss rate, we construct a full model of mass loss that can accurately predict the evolution of a subhalo in terms of a single parameter $\eta_{\rm eff}$. This parameter can be fully determined from the initial orbital and halo properties, and does not require calibration with numerical simulations.

The structural properties of individual dark matter haloes, including shape, spin, concentration, and substructure, are linked to the halo's growth history, but the exact connection between the two is unclear. One open question, in particular, is the effect of major mergers on halo structure. We have performed a large set of simulations of binary equal-mass mergers between isolated haloes with various density profiles, to map out the relationship between the initial conditions and merger parameters and the structure of the final remnant. In this paper we describe our initial set-up and analysis methods, and report on the results for the size, shape, and spin of the merger remnant. The outcomes of mergers are most easily understood in terms of a scaled dimensionless energy parameter $\kappa$ and an angular momentum (or spin) parameter $\lambda$. We find that the axis ratio $c/a$ scales roughly linearly with energy $\kappa$ while the axis ratio $c/b$ scales linearly with spin $\lambda$. Qualitatively, mergers on radial orbits produce prolate remnants, while mergers on tangential orbits produce oblate remnants. The spin of the remnant can be predicted from angular momentum conservation, while the overall size changes as $\sim \kappa^{-5}$, as expected from self-similar scaling at constant mean density. We discuss potential cosmological applications for these simple patterns.

Several lines of evidence suggest that as dark matter haloes grow their scale radius increases, and that the density in their central region drops. Major mergers seem an obvious mechanism to explain both these phenomena, and the resulting patterns in the concentration--mass--redshift relation. To test this possibility, we have simulated equal-mass mergers between haloes with a variety of cosmological density profiles, placed on various different orbits. The remnants typically have higher densities than the initial conditions, but differ only slightly from self-similar scaling predictions. They are reasonably well fit by Einasto profiles, but have parameters distinct from those of the initial conditions. The net internal energy available to the merger remnant, relative to the internal energy of the initial conditions, $\kappa$, has the greatest influence on the properties of the final mass distribution. As expected, energetic encounters produce more extended remnants while mergers of strongly bound systems produce compact remnants. Surprisingly, however, the scale radius of the density profile shows the opposite trend, {\it increasing} in the remnants of low-energy encounters relative to energetic ones. Also even in the most energetic encounters, the density within the scale radius decreases only slightly (by 10--20\%), while for very low-energy systems it increases significantly after the merger. We conclude that while major mergers can produce remnants that are more diffuse at large radii, they are relatively ineffective at changing the central densities of haloes, and seem unlikely to explain the mean trends in the concentration--mass--redshift relation.

We propose a new method for generating equilibrium models of spherical systems of collisionless particles that are finite in extent, but whose central regions resemble dark matter halos from cosmological simulations. This method involves iteratively removing unbound particles from a Navarro-Frenk-White profile truncated sharply at some radius. The resulting models are extremely stable, and thus provide a good starting point for N-body simulations of isolated halos. We provide a code to generate such models for NFW and a variety of other common density profiles. We then develop an analytic approximation to this truncated distribution function. Our method proceeds by analogy with the King model, truncating and shifting the original distribution function of an infinitely extended Navarro-Frenk-White profile in energy space. We show that the density profiles of our models closely resemble the tidally truncated density profiles seen previously in studies of satellite evolution. Pursuing this analogy further with a series of simulations of tidal mass loss, we find that our models provide a good approximation to the full distribution function of tidally stripped systems, thus allowing theoretically motivated phase-space calculations for such systems.