---
title: Non-Linear Linear Alignment in Cosmology
url: https://www.emergentmind.com/topics/non-linear-linear-alignment-nla
type: topic
---

# Non-Linear Linear Alignment in Cosmology

Searching arXiv for the cited NLA intrinsic-alignment papers and closely related work.
Non-Linear Linear Alignment (NLA) is a phenomenological intrinsic-alignment model used in weak-lensing cosmology to describe how galaxy shapes correlate with the large-scale tidal field while replacing the linear matter power spectrum of the original linear-alignment prescription with a non-linear matter power spectrum. In this sense, the phrase “non-linear” refers to the use of $P_{\rm NL}$ or $P_\delta^{\rm nl}(k,z)$ rather than to a non-linear shape response in the baseline ansatz itself. NLA has become a standard description of intrinsic alignments in two-point cosmic-shear analyses, has been directly calibrated on spectroscopic galaxy samples, and has also been the subject of perturbative critiques showing where it breaks down and how it can be extended [1811.09598] [1708.09247] [2509.25166].

## 1. Definition and theoretical origin

The starting point is the “instantaneous tidal-alignment” or linear-alignment (LA) model, in which galaxy shapes—primarily luminous red or elliptical galaxies—respond linearly to the large-scale tidal field
$$
s_{ij}(x)\equiv\left(\frac{\partial_i\partial_j}{\nabla^2}-\frac{\delta_{ij}}{3}\right)\delta(x).
$$
In the notation summarized by Blazek et al., the LA model writes
$$
P_{II}^{\rm LA}(k,z)=F^2(z)\,P_{\rm lin}(k,z),\qquad
P_{GI}^{\rm LA}(k,z)=F(z)\,P_{\rm lin}(k,z),
$$
with
$$
F(z)\equiv -\,C_1(z)\,\rho_{\rm crit,0}\,\Omega_{m,0}\,G^{-1}(z),
$$
where $G(z)$ is the linear growth factor normalized as $G(0)=1$ [1708.09247].

The NLA prescription keeps the same tidal-alignment ansatz but replaces the linear matter power spectrum by a non-linear one:
$$
P_{II}^{\rm NLA}(k,z)=F^2(z)\,P_{\rm NL}(k,z),\qquad
P_{GI}^{\rm NLA}(k,z)=F(z)\,P_{\rm NL}(k,z).
$$
This replacement is explicitly described as a phenomenological fix. The 2025 simulation study states that observationally “the alignments appear stronger than can be captured by purely linear theory,” and that most analyses find a better fit when $P_\delta^{\rm lin}\to P_\delta^{\rm NL}$, for example via Halofit. That same study describes the model as yielding percent-level agreement with data on scales $5$–$100'$ in the contexts tested there [2509.25166].

A common misconception is that the name implies a non-linear alignment law. The baseline NLA construction does not do that: it preserves linear tidal coupling and changes the matter clustering input. This distinction matters because later perturbative work identifies genuinely higher-order intrinsic-alignment operators beyond NLA.

## 2. Mathematical formulation and standard parameterizations

In the empirical calibration carried out with KiDS, GAMA, and SDSS, the NLA model is written directly in terms of the intrinsic–intrinsic and matter–intrinsic three-dimensional power spectra:
$$
P_{II}(k,z)=A_{\rm IA}^2
\biggl[C_1\,\frac{\bar\rho(z)\,a^2(z)}{D(z)}\biggr]^2
P_{\delta}^{\rm nl}(k,z),
$$
$$
P_{\delta I}(k,z)=-\,A_{\rm IA}
\biggl[C_1\,\frac{\bar\rho(z)\,a^2(z)}{D(z)}\biggr]
P_{\delta}^{\rm nl}(k,z).
$$
Here $A_{\rm IA}$ is a dimensionless alignment amplitude; $C_1=5\times10^{-14}\,h^{-2}M_\odot^{-1}\,\mathrm{Mpc}^3$ normalizes the original “SuperCOSMOS” measurement; $\bar\rho(z)$ is the mean matter density; $a(z)=(1+z)^{-1}$ is the scale factor; $D(z)$ is the linear growth factor with $D(0)=1$; and $P_\delta^{\rm nl}(k,z)$ is the non-linear matter power spectrum [1811.09598].

A luminosity-dependent generalization, denoted “NLA-$\beta$,” replaces the amplitude by
$$
A_{\rm IA}\longrightarrow A_{\beta}\,\bigl\langle L/L_{\rm piv}\bigr\rangle^{\beta},
$$
where $L_{\rm piv}\simeq4.6\times10^{10}\,L_\odot$ corresponds to $M_r=-22$, and $\langle L/L_{\rm piv}\rangle$ is the sample mean. Some analyses also introduce a redshift-evolution factor
$$
A_{\rm IA}\,\bigl[(1+z)/(1+z_0)\bigr]^{\eta}.
$$
In the KiDS+GAMA calibration no $\eta$-term was fitted because the redshift baseline $z\lesssim0.5$ was short; the summary states that effectively $\eta\approx0$ [1811.09598].

In tomographic cosmic-shear notation, the same model can be projected into angular statistics as
$$
C^{ij,II}_\ell = \int \frac{n_i(\chi)n_j(\chi)}{\chi^2}\,
P_{II}\bigl(k=(\ell+1/2)/\chi,z(\chi)\bigr)\,d\chi,
$$
$$
C^{ij,GI}_\ell = \int \frac{q_i(\chi)n_j(\chi)+q_j(\chi)n_i(\chi)}{\chi^2}\,
P_{GI}(k,z)\,d\chi,
$$
followed by
$$
\xi_{+/-}^{ij}(\theta)=\frac{1}{2\pi}\int_0^\infty
\left[C_\ell^{ij,GG}+C_\ell^{ij,GI}+C_\ell^{ij,II}\right]
J_{0/4}(\ell\theta)\,\ell\,d\ell.
$$
The 2025 implementation notes that many cosmic-shear codes equivalently write $P_{\delta I}(k,z)=A_{\rm IA}(z)\,P^{NL}_{\delta\delta}(k,z)$, with $A_{\rm IA}(z)=-\,C_1(z)a^2(z)$ absorbing the normalization and redshift dependence [2509.25166].

## 3. Empirical calibration in KiDS, GAMA, and SDSS

A direct calibration of NLA amplitudes was obtained using r-band shape measurements from KiDS-450 imaging over the three equatorial GAMA fields, totaling $180\,\mathrm{deg}^2$, with spectroscopic redshifts from GAMA to $r<19.8$ at $\ge98\%$ completeness, together with an SDSS Main sample with $r<17.77$ over $\simeq3340\,\mathrm{deg}^2$. Shapes were measured with the DEIMOS moment-based method for KiDS+GAMA and with the “re-Gaussianisation” pipeline for SDSS Main. Galaxies were separated into red and blue populations with the rest-frame colour cut $g-r=0.66$; GAMA was further split at $z=0.26$ into the two redshift bins Z1 and Z2, yielding Z1B, Z1R, Z2B, and Z2R subsamples [1811.09598].

The analysis estimated the galaxy–intrinsic shear correlation $\hat\xi_{g+}(r_p,\Pi)$ via the M06 estimator, subtracting random–shape pairs to remove additive systematics, and projected it as
$$
w_{g+}(r_p)=\int_{-\Pi_{\max}}^{\Pi_{\max}}\xi_{g+}(r_p,\Pi)\,d\Pi,
\qquad \Pi_{\max}=60\,h^{-1}\mathrm{Mpc}.
$$
Galaxy clustering was estimated with the Landy–Szalay estimator and projected into $w_{gg}(r_p)$. The two observables were then fitted jointly, with $w_{gg}$ constraining the linear galaxy bias $b_g$, on scales $r_p>6\,h^{-1}\,\mathrm{Mpc}$ [1811.09598].

The main parameter constraints were the following.

| Sample/model | Constraint | Interpretation |
|---|---:|---|
| Blue galaxies | $A_{\rm IA}^B = 0.21^{+0.37}_{-0.36}$ | Consistent with zero |
| Red galaxies | $A_{\rm IA}^R = 3.18^{+0.47}_{-0.46}$ | $9\sigma$ detection of radial alignments |
| Combined sample | $A_{\rm IA}^{\rm all} = 1.06^{+0.47}_{-0.46}$ | No colour split |
| Red NLA-$\beta$ | $\beta_R = 0.18^{+0.20}_{-0.22}$ | Consistent with zero |
| Blue NLA-$\beta$ | $\beta_B = 2.06^{+2.20}_{-2.82}$ | Consistent with zero |

These measurements were summarized as showing no significant detection of blue-galaxy alignments, robust detections for red galaxies, and no evidence for any scaling of alignments with galaxy luminosity. Additional checks on central versus satellite galaxies revealed that large-scale GI correlations above $6\,h^{-1}\,\mathrm{Mpc}$ arise almost entirely from red central galaxies, while red satellites, blue centrals, and blue satellites show no significant large-scale alignment. This suggests that colour alone does not exhaust the relevant astrophysical degrees of freedom.

## 4. Use in weak-lensing inference

The KiDS+GAMA calibration was explicitly translated into informative priors for cosmic-shear analyses:
$$
A_{\rm IA}^R=3.18\pm0.47,\qquad A_{\rm IA}^B=0.21\pm0.37.
$$
These were proposed as Gaussian priors on the red and blue NLA amplitudes in weak-lensing likelihoods [1811.09598].

Fisher forecasts for a completed KiDS-like $1{,}350\,\mathrm{deg}^2$ survey, tomographically colour-split into five bins each for red and blue galaxies, showed that applying these intrinsic-alignment priors can reduce the marginalized uncertainties on
$$
S_8\equiv\sigma_8\sqrt{\Omega_m/0.3}
$$
and on the dark-energy equation-of-state parameter $w_0$ by up to $62\%$ and $51\%$, respectively. The summary specifies that these are the maximal gains when also marginalizing over realistic photo-$z$ biases [1811.09598].

The same work argues that these priors improve upon de facto wide priors that allow unrealistic levels of intrinsic-alignment contamination. At the same time, it warns that modelling alignments only through a red/blue split may be insufficient when the fraction of central and satellite galaxies varies with depth or survey footprint. A plausible implication is that informative priors on $A_{\rm IA}$ are most reliable when the sample definition is close to that used in the calibration.

## 5. Perturbative generalizations and the breakdown of NLA

Blazek et al. formulate intrinsic alignments as a perturbative expansion analogous to galaxy biasing. Up to second order, the intrinsic shear field is expanded in all spin-2, traceless combinations of the density, tidal, and velocity-shear fields:
$$
\gamma^I_{ij}(x)=C_1\,s_{ij}(x)
+ C_2\!\left[s_{ik}(x)s_{kj}(x)-\frac13\delta_{ij}s_{kl}s_{kl}\right]
+ C_{1\delta}\,[\delta(x)s_{ij}(x)]
+ C_t\,t_{ij}(x)+\cdots.
$$
In this framework, the NLA model is the subset retaining the linear tidal term while replacing the linear matter spectrum by a non-linear one [1708.09247].

The physical assumptions of LA and NLA are stated explicitly: galaxy shapes respond linearly to the large-scale tidal field, higher-order responses are neglected, and the time delay between structure formation and galaxy alignment is ignored. The summary further states that the model is valid only on large, quasi-linear scales where non-Gaussian, nonlinear intrinsic-alignment mechanisms such as tidal torquing are subdominant; empirically, one typically restricts NLA to $k\lesssim0.1$–$0.2\,h\,\mathrm{Mpc}^{-1}$ and $z\lesssim1$, beyond which omission of one-loop and higher corrections leads to order-unity errors [1708.09247].

The one-loop corrections computed in that work include quadratic alignment and density-weighting terms and are reported to contribute $O(1)$ corrections to the total intrinsic-alignment signal at $k\sim0.1\,h^{-1}{\rm Mpc}$. In an LSST-like forecast with five tomographic bins, $30\,\mathrm{gal}\,\mathrm{arcmin}^{-2}$, area $18\,000\,\mathrm{deg}^2$, and $\ell\in[100,1000]$, a data vector generated with the full one-loop intrinsic-alignment model but analyzed with NLA yields
$$
w(a_{\rm piv})|_{\rm true}=-1,\qquad
w(a_{\rm piv})|_{\rm NLA\,fit}=-2.13\pm0.13,
$$
described as a $>8\sigma$ excursion. Allowing the NLA amplitude to scale as $(1+z)^\alpha$ reduces the bias to
$$
w(a_{\rm piv})=-1.13\pm0.13,
$$
but does not eliminate it. By contrast, fitting the full perturbative intrinsic-alignment model gives
$$
w(a_{\rm piv})=-1.03\pm0.08,
$$
with bias $\lesssim0.4\sigma$, while marginalizing over $C_1$ and $C_2$ increases the uncertainty by only $\sim2\%$ relative to the no-intrinsic-alignment analysis [1708.09247].

These results are often interpreted as evidence that NLA is adequate as a large-scale effective model for current two-point analyses but not as a complete description for Stage IV weak-lensing inference.

## 6. Simulation-level implementations and extended model families

A later simulation study embeds NLA directly into weak-lensing simulations constructed from the Outer-Rim N-body light cone, using $10{,}240^3$ particles in a $4.2\,\mathrm{Gpc}$ box and $57$ curved-sky mass shells up to $z=3$ on Healpix $\mathrm{nside}=8192$. For each shell, the projected tidal field is obtained from the density map, the B-mode is set to zero, and $Q,U$ maps are reconstructed. The resulting tidal-field maps are smoothed with a Gaussian beam of width $\sigma_G$, with tests at $\sigma_G=0.1$ and $0.5\,h^{-1}\,\mathrm{Mpc}$, and then downgraded to $\mathrm{nside}=4096$ [2509.25166].

Mock source galaxies are either placed randomly on the octant, corresponding to the NLA assumption of no density weighting, or are Poisson-sampled from the smoothed shells with mean proportional to $[1+b\,\delta]$ for linear-bias tracers, or from halo-occupation-distribution catalogues for non-linear-bias tracers used in $\delta$-NLA. Intrinsic ellipticities are assigned by interpolating the tidal field at each galaxy position and then combined with simulated reduced shear $g$ through
$$
\epsilon^{\rm obs}=\frac{\epsilon^{\rm int}+g}{1+g\,\epsilon^{\rm int*}}.
$$
The NLA implementation is validated against CosmoSIS predictions using fully non-linear $P(k)$: agreement is reported at the $1$–$2\%$ level for $\theta\gtrsim5'$ in $\xi_+$ and $\theta\gtrsim20'$ in $\xi_-$ when $\sigma_G=0.5\,h^{-1}\,\mathrm{Mpc}$ [2509.25166].

Within this simulation framework, NLA is compared to several related models. The extended NLA model, also called e-NLA or $\delta$-NLA, multiplies the NLA ellipticity by $(1+b_{TA}\delta)$:
$$
\epsilon^{\delta\text{-NLA}}=\epsilon^{\rm NLA}[1+b_{TA}\delta].
$$
The tidal-torque (TT) model instead uses a quadratic coupling,
$$
\epsilon\sim s_{ik}s_{kj}-\tfrac13\delta_{ij}s^2,
$$
with $C_2\propto A_2 D^{-2}(z)$. The study also considers $\delta$-TT and HOD-TATT variants. It reports that the $\delta$-NLA model has by far the largest impact on most non-Gaussian probes, at times more than twice the strength of the NLA; that minima, void profiles, and the lensing PDF are the best probes for model rejection because large differences between models appear in under-dense regions; and that the third-order aperture mass statistic $M^3_{ap}$ and the integrated three-point functions are particularly sensitive to the source-clustering term when low-redshift data are included, often exceeding a $20\%$ impact on the data vector [2509.25166].

In the same work, nested-sampling likelihood analyses varying $\{\Omega_m,S_8,A_{IA}\}$ recover $A_{IA}\approx1$ and values of $S_8$ and $\Omega_m$ consistent with the true Outer-Rim cosmology, within sampling variance, when the simulated data are generated with NLA and analyzed with the corresponding NLA theory. Opening the full TATT parameter space introduces the expected degeneracy between $A_{IA}$ and the quadratic amplitude $C_2$, but $S_8$ remains stable. This suggests that NLA remains a robust phenomenological description on the scales used by current two-point lensing analyses, even as higher-order statistics and next-generation surveys motivate more general intrinsic-alignment models.

Source: https://www.emergentmind.com/topics/non-linear-linear-alignment-nla