---
title: Tidal Adaptive Softening in Cosmological Simulations
url: https://www.emergentmind.com/topics/tidal-adaptive-softening
type: topic
---

# Tidal Adaptive Softening in Cosmological Simulations

Searching arXiv for recent and foundational papers on tidal adaptive softening and closely related adaptive force-softening methods.
Searching arXiv for tidal adaptive softening, NovA, and adaptive gravitational softening.
Tidal adaptive softening denotes a class of adaptive gravitational force-resolution schemes for collisionless $N$-body and related cosmological simulations in which the effective softening length, or the refinement level that determines force resolution, is regulated by the local anisotropy or strength of the gravitational tidal field rather than by density alone. In the contemporary literature, the term covers two closely related lines of development. One is an explicitly tidal-tensor-based prescription, in which each particle’s softening is set from the norm of the local tidal tensor, as in the GIZMO implementation of conservative adaptive softening [2212.06851] and its cosmological assessment in warm dark matter simulations [2507.16930]. The other is an anisotropy-gated adaptive mesh refinement strategy, exemplified by NovA in RAMSES, which withholds refinement in sheets and filaments and allows it only when collapse is sufficiently isotropic, thereby making the effective softening conservative in regions of one- or two-dimensional collapse [1503.02689]. Across these variants, the shared objective is to prevent the force resolution from becoming spuriously fine in anisotropic structures, where discreteness noise and two-body scattering seed artificial fragmentation.

## 1. Conceptual definition and scope

Tidal adaptive softening is motivated by a specific failure mode of conventional collisionless simulations: discretization of a collisionless fluid with finite-mass macroparticles generates localized force perturbations that can dominate the physical small-scale signal when the matter power spectrum is suppressed, especially in warm dark matter cosmologies [2507.16930]. In standard formulations with fixed or density-only adaptive softenings, force resolution often tightens as soon as local density increases, even if the physical collapse is highly anisotropic. In sheets and filaments this causes the softening to track the shortest spacing while remaining too small along the long axis, amplifying two-body noise and producing the familiar “beads-on-a-string” fragmentation [1503.02689].

The central idea of tidal adaptive softening is that force resolution should respond to the geometry of collapse set by the tidal field. If matter is compressing along all principal directions, finer resolution is justified; if compression occurs along only one or two directions, refinement should remain conservative [1503.02689]. This logic can be implemented directly through the tidal tensor, using the local Hessian of the gravitational potential to set a particle softening scale [2212.06851, 2507.16930], or indirectly through an isotropy test on the particle configuration inside an AMR cell, using the inertia tensor as a proxy for local tidal anisotropy [1503.02689]. In both cases, the purpose is not merely adaptive resolution in the generic sense, but adaptive resolution conditioned on anisotropic collapse.

This distinguishes tidal adaptive softening from density-based adaptive softening, including earlier conservative schemes in GADGET that adapt the softening length as $h \propto \rho^{-1/3}$ via a target neighbor number [1107.2942]. Those methods improve spatial resolution in dense regions and can sharpen inner profiles, but they are not tidal-field adaptations and do not explicitly suppress refinement in filamentary or sheet-like collapse [1107.2942].

## 2. Mathematical formulations

Two mathematical formulations dominate the literature.

The explicitly tidal formulation defines the tidal tensor as
$$
T_{ij} \equiv \frac{\partial^2 \Phi}{\partial x_i \partial x_j},
$$
or, in the GIZMO-based formulation,
$$
T = - \nabla \otimes \nabla \hat{\Phi},
$$
where $\hat{\Phi}$ is a local potential estimator [2212.06851, 2507.16930]. A scalar tidal strength is then constructed from the Frobenius norm,
$$
||T||^2 = \sum_i \lambda_i^2,
$$
with $\lambda_i$ the eigenvalues of $T$ [2507.16930]. The softening rule is
$$
\epsilon_a = \xi \left(\frac{G m_a}{||T||_a}\right)^{1/3},
$$
with $\xi \sim 2$ in the GIZMO implementation [2507.16930]. The 2022 conservative formulation further introduces smooth floor or ceiling bounds through
$$
\epsilon_i = [ \epsilon_b^{n} + \epsilon_{tid,i}^{n} ]^{1/n},
$$
with $\epsilon_{tid,i} = \xi ( G m_i / ||T||_i )^{1/3}$ and $n=\pm 1$ typical [2212.06851]. This construction is intended to preserve differentiability when imposing minimum or maximum softenings.

The NovA formulation in RAMSES does not define $\epsilon$ directly from $T_{ij}$. Instead, it computes the moment of inertia tensor of the particle distribution within each cell about the cell’s centroid,
$$
I_{ij} = \sum_a m_a (x_a - x_c)_i (x_a - x_c)_j,
$$
diagonalizes it to obtain eigenvalues $\lambda_1 \ge \lambda_2 \ge \lambda_3$, and forms the anisotropy measure
$$
q_i \equiv \lambda_1 / \lambda_3.
$$
A cell is refined only if
$$
N_{\rm cell} > N_c \quad \text{and} \quad q_i < \eta,
$$
where $\eta$ is a tolerance above unity calibrated from shot-noise realizations [1503.02689]. For $N_c = 32$, the reported values are $\sigma_q \approx 0.169$ and $\eta \approx 1.084$ [1503.02689]. This refinement gate embodies the rule “refine only if collapse is occurring along all three axes” [1503.02689].

The physical interpretation given for both approaches is similar. In the Zel’dovich picture, sheets correspond to collapse along one principal direction, filaments to collapse along two, and halos to collapse along three [1503.02689]. A tidal-aware method therefore aims to reduce $\epsilon$ only when the local state resembles bona fide three-dimensional collapse.

## 3. Conservative dynamics and implementation frameworks

A major issue with any adaptive softening law is conservation. If the softening varies in space and time, the equations of motion must be modified to account for the induced dependence of the gravitational Lagrangian on particle positions through the softening field. Density-based conservative softening in GADGET derived this through a symmetrized Lagrangian and explicit $\nabla h$ correction terms [1107.2942]. The same conservative logic was generalized to arbitrary softening rules, including tidal softening, in the GIZMO formulation [2212.06851].

For adaptive softening in GADGET, the Lagrangian is written as
$$
L = \sum_{i=0}^N \frac{1}{2}m_i\mathbf{v}_i^2 - \frac{G}{2}\sum_{i=0}^N \sum_{j=0}^N m_i m_j \left[ \frac{\phi_{ij}(h_i) + \phi_{ij}(h_j)}{2}\right],
$$
and the resulting equations of motion contain the standard symmetric softened force plus correction terms involving
$$
\zeta_i \equiv \frac{\partial h_i}{\partial \rho_i}\sum_{k=0}^N m_k \frac{\partial \phi_{ik}(h_i)}{\partial h_i},
$$
and
$$
\Omega_i \equiv 1 - \frac{\partial h_i}{\partial \rho_i}\sum_{k=0}^N m_k \frac{\partial W_{ik}(h_i)}{\partial h_i}
$$
[1107.2942]. The paper stresses that omitting the correction term causes loss of energy conservation and biases the low-mass halo population [1107.2942].

The 2022 GIZMO treatment extends this beyond density-based rules. Starting from a softened pairwise potential $\tilde{\phi}_{ab}(r_{ab}; \epsilon_a, \epsilon_b)$, the acceleration is
$$
\frac{d p_a}{dt} = - \sum_b m_a m_b [ \nabla_a \tilde{\phi}_{ab} + \Upsilon_{ab} - \Upsilon_{ba} ],
$$
where $\Upsilon_{ab}$ and $\Upsilon_{ba}$ are the energy- and momentum-conserving “grad-$\epsilon$” terms required for arbitrary adaptive softenings [2212.06851]. For the tidal rule, these corrections depend on derivatives of the local tidal estimator and on
$$
\frac{\partial \epsilon_a}{\partial G_{a,\alpha\beta}} = - \epsilon_a \left[ \frac{T_{a,\alpha\beta}}{3 ||T_a||^2} \right]
$$
[2212.06851].

The same work also emphasizes symmetrization of pairwise softenings. It argues that force-averaging and simple softening-averaging can produce large errors or over-scattering when $\epsilon_a \neq \epsilon_b$, and recommends a “maximum-softening” symmetrization,
$$
\tilde{\phi}_{ab} = \phi^0(r; \max\{\epsilon_a,\epsilon_b\}),
$$
or smooth approximations that converge to this limit [2212.06851]. The implementation is public in GIZMO [2212.06851].

## 4. Relation to local inter-particle spacing and anisotropic collapse

A recurring theoretical theme is the relation between softening and inter-particle spacing. In adaptive schemes tied to density, a common rule is
$$
\epsilon \sim 2 d,
$$
with $d$ the local mean inter-particle spacing [1503.02689]. For isotropic collapse this is natural, because all principal spacings are comparable. In anisotropic collapse, however, the relevant spacings differ by direction. If the principal axis extents are associated with the eigenvalues of the inertia tensor, then a density-only criterion effectively ties $\epsilon$ to the shortest spacing, which is numerically hazardous in filaments and sheets [1503.02689].

NovA’s anisotropy gate is designed to prevent this. By withholding refinement until $q_i < \eta$, it keeps the effective softening of order twice the largest inter-particle spacing in the cell,
$$
\epsilon \gtrsim 2 \max[d_i] \equiv 2 a,
$$
rather than allowing it to undershoot along the long axis [1503.02689]. This is why NovA is described as a “tidal-aware” adaptive softening method even though the algorithm is expressed through $I_{ij}$ instead of $T_{ij}$ [1503.02689].

The explicitly tidal rule reaches a similar endpoint from a different direction. In a smooth Keplerian background with enclosed mass $M_{\rm enc}$ at scale $\ell$, the tidal norm is estimated as
$$
||T|| \approx \sqrt{6} \, G M_{\rm enc}/\ell^3,
$$
leading to a Hill-radius scaling $r_{\rm tidal}\approx (m/M_{\rm enc})^{1/3}\ell$ [2212.06851]. For a cubic spline kernel, the reported relation is $\epsilon/r_{\rm tidal} \approx 0.74\,\xi$ [2212.06851]. In isotropic cusps with local density $\rho$, the tidal norm tends to
$$
||T|| \to (4\pi/\sqrt{3})G\rho,
$$
so that
$$
\epsilon \approx 0.52\,\xi \,(m/\rho)^{1/3},
$$
that is, of order the local mean inter-particle spacing [2212.06851]. This suggests that the tidal rule interpolates between a Hill-radius criterion in external tides and a density-like spacing criterion when local self-gravity dominates.

## 5. Numerical performance in idealized and cosmological tests

The numerical record is mixed and depends strongly on the problem class.

In controlled anisotropic-collapse tests, tidal-aware methods show clear gains. NovA was tested on the asymmetric perturbed plane wave or “Valinia” collapse and produced no numerical fragmentation while still capturing fine caustics and shells around the collapsing regions [1503.02689]. In the same class of idealized filamentary-collapse tests, the 2025 cosmological study reports that tidal adaptive softening in GIZMO significantly improves force accuracy over fixed softening [2507.16930]. Using an anti-symmetric perturbed plane-wave collapse with $L = 10$, $N = 256^3$, $z_{\rm ini} = 99$, a $512^3$ PM mesh, and $\epsilon_{\rm fixed} = (\text{mean inter-particle spacing})/40$, the authors find that TAS yields approximately $40\%$ more particles with zero spurious $z$-shift than fixed softening, and suppresses large artificial displacements along the unperturbed axis [2507.16930].

In warm dark matter cosmological runs, NovA also improves over standard RAMSES. For a thermal relic with $m_\chi = 0.2\,{\rm keV}$ in cosmologies with $\Omega_0 = 0.27$, $\Omega_\Lambda = 0.73$, $h = 0.705$, and $\sigma_8 = 0.81$, NovA converges more rapidly than standard mass-only refinement and produces little or no spurious halos on small scales, with nearly an order-of-magnitude fewer low-mass halos below $M_{\rm fs} \approx 2.26 \times 10^9\,M_\odot$ at $N = 1024^3$ compared to standard runs [1503.02689]. However, the same work notes that AHF likely mislabels some caustics and criss-crossing filaments as halos and that one or two particularly massive filaments appear to fragment in any NovA version where refinement is allowed [1503.02689].

By contrast, the 2025 study of Mostoghiu Paun et al. reports that explicit tidal adaptive softening does not substantially reduce spurious halos in full cosmological WDM simulations when paired with optimized initial conditions [2507.16930]. Their production runs use monofonIC with 3LPT, anti-aliasing via Orszag’s $3/2$-rule padding, optional PLT corrections, a late start at $z_{\rm ini}=39$, $L = 25$, $N=512^3$ dark matter particles, a $1024^3$ PM mesh, and $\epsilon_{\rm fixed}=1/40$ of the mean inter-particle spacing at the initial conditions [2507.16930]. In these runs, both fixed and TAS cases show the characteristic low-mass upturn in the cumulative halo mass function near $M_{200} \approx 10^9$, indicating spurious halos [2507.16930]. The ratio TAS/fixed peaks around that scale, with TAS yielding up to $1.14\times$ more halos than fixed, corresponding to approximately $+14\%$ near the 100-particle mass [2507.16930]. Table-based counts obtained with an empirical proto-halo sphericity cut give 10,452 total WDM objects with 5,973 flagged as spurious in the fixed-softening run, versus 13,060 total with 8,004 flagged as spurious in the TAS run [2507.16930].

The same study also reports a measurable shift in WDM halo formation times under TAS. Using the half-maximum mass scale factor $a_{\rm hmm}$, the mean values in WDM are $\bar{a}_{\rm hmm,fixed}\approx 0.44$ and $\bar{a}_{\rm hmm,adapt}\approx 0.41$, a difference $\Delta a \approx 0.03 \approx 0.5\,{\rm Gyr}$; no analogous shift is found in CDM, where the means are both approximately $0.34$ [2507.16930]. This indicates that the softening prescription can bias inferred formation histories even when it does not cure fragmentation.

## 6. Tidal adaptive softening for subhalo tidal evolution

A more recent and distinct use of the term concerns the numerical convergence of tidally stripped subhalos rather than artificial filament fragmentation. In this context, tidal adaptive softening means adapting force resolution to the instantaneous tidal radius of each subhalo [2510.26901]. The key result is a universal force-resolution criterion: in AMR simulations, the minimum instantaneous cell size must satisfy
$$
\Delta x_{\rm min}^{\rm ins} \leq \frac{r_{\rm t,min}}{20},
$$
and in tree codes the Plummer-equivalent softening must satisfy
$$
\varepsilon_{\rm P} \leq \frac{r_{\rm t,min}}{20}
$$
[2510.26901]. Here $r_{\rm t,min}$ is the minimum tidal radius experienced since infall.

The tidal radius is defined through
$$
r_{\rm t} \equiv \left[\frac{G\, m(r_{\rm t})}{\Omega^2 - \frac{d^2 \Phi_{\rm h}}{dr^2}}\right]^{1/3},
$$
with $\Omega \equiv v/r$, $m(r)$ the instantaneous enclosed bound mass profile of the remnant, and $\Phi_{\rm h}$ the host-halo potential [2510.26901]. The minimum tidal radius is
$$
r_{\rm t,min} \equiv \min[r_{\rm t}(t)], \quad t \in [t_{\rm infall}, t_{\rm now}] .
$$
Across AMR simulations spanning subhalo concentrations $4 \le c \le 40$, constant anisotropies $-0.5 \le \beta \le +0.5$, infall mass ratios $-2.5 \le \log(m_{\rm s,0}/M_{\rm h}) \le -4$, circular and eccentric orbits with $0.2 \le r_{\rm circ}(E_{\rm orb})/r_{\rm vir,h} \le 0.8$, and profiles including NFW, Plummer, and Hernquist, the $20$-element rule is reported to be universal [2510.26901].

When the rule is violated, the bound mass fraction is systematically underestimated. The reported correction fit is
$$
\frac{f_{\rm bound}}{f_{\rm bound}^{\rm truth}} = 0.5\left\{\operatorname{erf}\left[\,2.1\,\log\left(0.24\,\frac{r_{\rm t,min}}{\Delta x_{\rm min}^{\rm ins}}\right)\right]+1\right\},
$$
where $f_{\rm bound}^{\rm truth}$ is measured from a very-high-resolution reference run [2510.26901]. The same work also gives a universal scatter formula for discreteness-noise-driven variance in the bound mass fraction,
$$
\sigma_{\log f} = 1.8\,N_{\rm par}^{-0.5}\,f_{\rm bound}^{-0.6}\,\left(1 - f_{\rm bound}^{0.5}\right)^{0.8},
$$
and notes that as many as $50\%$ of first-order subhalos in the Bolshoi AMR simulation are likely to be force and/or mass unresolved [2510.26901].

This subhalo-focused literature broadens the meaning of tidal adaptive softening. Instead of using the tidal tensor to choose a particle-scale smoothing directly, it uses a physically relevant tidal scale, the instantaneous tidal radius, as the target quantity that must be resolved. A plausible implication is that the term now encompasses both a local field-based prescription and an object-centered adaptive-resolution criterion.

## 7. Debates, limitations, and current interpretation

Three issues dominate current discussion.

The first is whether tidal adaptive softening actually solves artificial fragmentation in cosmological WDM runs. The answer depends on the variant and setup. NovA substantially suppresses spurious refinement in anisotropic structures and strongly reduces low-mass spurious halos relative to standard AMR [1503.02689]. However, explicit tidal-tensor softening in GIZMO, when paired with late-start 3LPT initial conditions, improves force accuracy in idealized filament collapse but does not remove the low-mass WDM mass-function upturn in full cosmological boxes and can slightly increase low-mass counts near the 100-particle scale [2507.16930]. The 2025 study interprets this to mean that, in such runs, the dominant source of spurious fragmentation lies deeper in the $N$-body discretization than what TAS alone addresses [2507.16930].

The second issue is the dependence of spurious-halo identification on post-processing choices. The empirical sphericity-based discriminator used in [2507.16930] is calibrated from proto-halo inertia tensors,
$$
I_{ij} = \sum_{\rm particles} m (\delta_{ij}|x|^2 - x_i x_j),
$$
with sphericity $s = \lambda_c/\lambda_a$, and an empirical mean relation
$$
\bar{s} = \alpha \log_{10}(M_{200,{\rm hmm}} / ?) + \beta,
$$
with $\alpha = 0.0624$ and $\beta = -0.0988$ below $M_{200,{\rm hmm}} = 10^{10}$ [2507.16930]. But that same paper shows that moving from $z_{\rm ini}=39$ to $z_{\rm ini}=99$ lowers mean proto-halo sphericities by approximately $0.03$, or about $6\%$, in both CDM and WDM, so the same empirical cut no longer cleanly separates the two [2507.16930]. This suggests that shape-based cleaning is sensitive to initial-condition choices.

The third issue concerns what should count as the most physically justified adaptive variable. Density-based softening in GADGET is conservative and improves small-scale clustering and inner-profile convergence, but it is not tidal-aware [1107.2942]. Tidal-tensor softening in GIZMO is explicitly translation- and Galilean-invariant and equivalence-principle respecting, avoids the gauge dependence of potential-based rules and the pathological behavior of acceleration-based rules, and imposes negligible timestep penalties compared to neighbor-based schemes [2212.06851]. NovA, although not expressed directly through $T_{ij}$, encodes the same collapse geometry through an isotropy criterion on the particle distribution and enforces a conservative softening in anisotropic regions [1503.02689]. These methods therefore occupy different points in the design space rather than constituting a single algorithm.

Taken together, the literature supports a restrained interpretation. Tidal adaptive softening is best understood as a family of tidal-aware force-resolution strategies that align numerical resolution with the geometry or strength of gravitational collapse. Its most robust successes are in controlling anisotropy-driven force errors, preserving conservative dynamics, and improving convergence for tidally stripped subhalos when the instantaneous tidal scale is explicitly resolved [2212.06851, 2510.26901]. Its efficacy as a standalone cure for cosmological warm-dark-matter fragmentation is more limited and remains contingent on the broader discretization, initialization, and halo-identification pipeline [1503.02689, 2507.16930].

Source: https://www.emergentmind.com/topics/tidal-adaptive-softening