Measurements at z = 0 of the linear bias parameter b₁ (upper panel) and the PNG-response parameter p (lower panel) for different codes. All simulations have been initialised at z = 49 using 2LPT and f_NL = 100 (see Section 3.2). All the points from the same code or configuration have a different marker and are grouped, from left to right, MonofonIC, 2LPTPNG,LPICOLA, and FastPM, with k_max = 1, k_max = 2/3, and k_max = 1/4. the colours indicate different halo mass bins (blue: 1.6 · 10¹² < M(h⁻¹ M_⨀) < 4.0 · 10¹², green: 4.0 · 10¹² < M(h⁻¹M_⨀) < 8.0 · 10¹², orange: 8.0 · 10¹² < M(h⁻¹M_⨀) < 2.4 · 10¹³, and red: 2.4 · 10¹³ < M(h⁻¹M_⨀) < 1.3 · 10¹⁵). The horizontal dashed lines in the upper panel show the values of b₁ measured in our reference simulation, using MonofonIC (Section 3.3). In the lower panel, the horizontal dashed line shows the universality relation for b_ϕ/p. The value of p measured from the simulation run with FastPM assuming k_max = 1/4 is ~2.7σ away from the reference for the first three mass bins (Section 4.2). The values of b₁ and p measured from all other codes and configurations agree within 1σ.