We present a new simulation technique in which any chosen mode $\mathbf{k}$ of the density contrast field can be amplified by an amplitude $\Delta$. These amplified-mode simulations allow us to study the response of the halo density field to a long-wavelength mode other than the DC mode. In this sense they are a generalization of the separate-universe simulations to finite-wavelength modes. In particular, we use these simulations to obtain robust measurements of the first higher-derivative bias of dark matter halos $b_{\nabla^2 \delta}$. We find a negative bias at all mass considered, roughly following the $-R_L^2(M)$ relation, the Lagrangian radius of halos squared, as naively expected. We compare our results with those obtained from a fit to the 1-loop halo-matter power spectrum, as well as with the recent results from Abidi and Baldauf (2018), and to the prediction from the peak theory.

We present direct measurements of cubic bias parameters of dark matter halos from the halo-matter-matter-matter trispectrum. We measure this statistic efficiently by cross-correlating the halo field measured in N-body simulations with specific third-order nonlocal transformations of the initial density field in the same simulation. Together with the recent Abidi & Baldauf (2018), these are the first measurements of halo bias using the four-point function that have been reported to date. We also obtain constraints on the quadratic bias parameters. For all individual cubic parameters involving the tidal field $\mathcal{K}_{ij}$, we find broad consistency with the prediction of the Lagrangian local-in-matter-density ansatz, with some indications of a positive Lagrangian coefficient $b_{\rm td}^L$ multiplying the time derivative of $\mathcal{K}_{ij}$. For the quadratic tidal bias ($b_{K^2}$), we obtain a significant detection of a negative Lagrangian tidal bias.

We measure the large-scale bias of dark matter halos in simulations with non-Gaussian initial conditions of the local type, and compare this bias to the response of the mass function to a change in the primordial amplitude of fluctuations. The two are found to be consistent, as expected from physical arguments, for three halo-finder algorithms which use different Spherical Overdensity (SO) and Friends-of-Friends (FoF) methods. On the other hand, we find that the commonly used prediction for universal mass functions, that the scale-dependent bias is proportional to the first-order Gaussian Lagrangian bias, does not yield a good agreement with the measurements. For all halo finders, high-mass halos show a non-Gaussian bias suppressed by 10-15% relative to the universal mass function prediction. For SO halos, this deviation changes sign at low masses, where the non-Gaussian bias becomes larger than the universal prediction.

We present precise measurements of the assembly bias of dark matter halos, i.e. the dependence of halo bias on other properties than the mass, using curved "separate universe" N-body simulations which effectively incorporate an infinite-wavelength matter overdensity into the background density. This method measures the LIMD (local-in-matter-density) bias parameters $b_n$ in the large-scale limit. We focus on the dependence of the first two Eulerian biases $b^E_1$ and $b^E_2$ on four halo properties: the concentration, spin, mass accretion rate, and ellipticity. We quantitatively compare our results with previous works in which assembly bias was measured on fairly small scales. Despite this difference, our findings are in good agreement with previous results. We also look at the joint dependence of bias on two halo properties in addition to the mass. Finally, using the excursion set peaks model, we attempt to shed new insights on how assembly bias arises in this analytical model.

We present accurate measurements of the linear, quadratic, and cubic local bias of dark matter halos, using curved "separate universe" N-body simulations which effectively incorporate an infinite-wavelength overdensity. This can be seen as an exact implementation of the peak-background split argument. We compare the results with the linear and quadratic bias measured from the halo-matter power spectrum and bispectrum, and find good agreement. On the other hand, the standard peak-background split applied to the Sheth & Tormen (1999) and Tinker et al. (2008) halo mass functions matches the measured linear bias parameter only at the level of 10%. The prediction from the excursion set-peaks approach performs much better, which can be attributed to the stochastic moving barrier employed in the excursion set-peaks prediction. We also provide convenient fitting formulas for the nonlinear bias parameters $b_2(b_1)$ and $b_3(b_1)$, which work well over a range of redshifts.

The dark matter halos that host galaxies and clusters form out of initial high-density patches, providing a biased tracer of the linear matter density field. In the simplest local bias approximation, the halo field is treated as a perturbative series in the average overdensity of the Lagrangian patch. In more realistic models, however, additional quantities will affect the clustering of halo-patches, and this expansion becomes a function of several stochastic variables. In this paper, we present a general multivariate expansion scheme that can parametrize the clustering of any biased Lagrangian tracer, given only the variables involved and their symmetry (in our case rotational invariance). This approach is based on an expansion in the orthonormal polynomials associated with the relevant variables, so that no renormalization of the coefficients ever occurs. We provide explicit expression for the series coefficients, or Lagrangian bias parameters, in the case of peaks of the linear density field. As an application of our formalism, we present a simple derivation of the original BBKS formula, and compute the non-Gaussian bias in the presence of a primordial trispectrum of the local shape.