diff --git a/docs/src/comparison.md b/docs/src/comparison.md
index df79c6c7..011ff298 100644
--- a/docs/src/comparison.md
+++ b/docs/src/comparison.md
@@ -14,7 +14,7 @@ using DataInterpolations
 using Unitful, UnitfulAstro
 using CairoMakie
 
-lmax = 6
+lmax = 10
 reionization = true
 Hswitch = 1
 Heswitch = 6
diff --git a/docs/src/components.md b/docs/src/components.md
index b6116417..4a8d9c25 100644
--- a/docs/src/components.md
+++ b/docs/src/components.md
@@ -143,7 +143,7 @@ initialization_equations(r)
 SymBoltz.photons
 ```
 ```@example components
-γ = SymBoltz.photons(g; lmax = 6)
+γ = SymBoltz.photons(g; lmax = 10)
 nothing # hide
 ```
 ```@raw html
@@ -251,7 +251,7 @@ initialization_equations(b)
 SymBoltz.massless_neutrinos
 ```
 ```@example components
-ν = SymBoltz.massless_neutrinos(g; lmax = 6)
+ν = SymBoltz.massless_neutrinos(g; lmax = 10)
 nothing # hide
 ```
 ```@raw html
@@ -278,7 +278,7 @@ SymBoltz.momentum_quadrature
 SymBoltz.massive_neutrinos
 ```
 ```@example components
-h = SymBoltz.massive_neutrinos(g; nx=4, lmax = 5)
+h = SymBoltz.massive_neutrinos(g; nx=4, lmax = 10)
 nothing # hide
 ```
 ```@raw html
diff --git a/src/models/cosmologies.jl b/src/models/cosmologies.jl
index d47d227a..cae9997e 100644
--- a/src/models/cosmologies.jl
+++ b/src/models/cosmologies.jl
@@ -1,6 +1,6 @@
 """
     ΛCDM(;
-        lmax = 6,
+        lmax = 10,
         recombination = true,
         reionization = true,
         Hswitch = 1,
@@ -23,7 +23,7 @@
 Create a ΛCDM model.
 """
 function ΛCDM(;
-    lmax = 6,
+    lmax = 10,
     recombination = true,
     reionization = true,
     Hswitch = 1,
diff --git a/src/models/neutrinos.jl b/src/models/neutrinos.jl
index f3cdd36c..0d015d23 100644
--- a/src/models/neutrinos.jl
+++ b/src/models/neutrinos.jl
@@ -1,9 +1,9 @@
 """
-    massless_neutrinos(g; lmax = 6, name = :ν, kwargs...)
+    massless_neutrinos(g; lmax = 10, name = :ν, kwargs...)
 
 Create a particle species for massless neutrinos in the spacetime with metric `g`.
 """
-function massless_neutrinos(g; lmax = 6, name = :ν, kwargs...)
+function massless_neutrinos(g; lmax = 10, name = :ν, kwargs...)
     description = "Massless neutrinos"
     ν = radiation(g; adiabatic = true, name, description, kwargs...) |> background |> complete
 
@@ -58,11 +58,11 @@ end
 
 # TODO: use vector equations and simplify loops
 """
-    massive_neutrinos(g; nx = 4, lmax = 4, name = :h, kwargs...)
+    massive_neutrinos(g; nx = 4, lmax = 10, name = :h, kwargs...)
 
 Create a particle species for massive neutrinos in the spacetime with metric `g`.
 """
-function massive_neutrinos(g; nx = 4, lmax = 4, name = :h, kwargs...)
+function massive_neutrinos(g; nx = 4, lmax = 10, name = :h, kwargs...)
     # compute numerical reduced momenta x = q*c / (kB*T) and Gaussian quadrature weights
     # for approximating integrals ∫dx x² f₀(x) g(x) for any g(x) over the infinite domain (0, ∞),
     # but change variables to transform it into a finite domain (0, 1)
diff --git a/src/models/photons.jl b/src/models/photons.jl
index 45dfc6a4..fd989737 100644
--- a/src/models/photons.jl
+++ b/src/models/photons.jl
@@ -1,9 +1,9 @@
 """
-    photons(g; polarization = true, lmax = 6, name = :γ, kwargs...)
+    photons(g; polarization = true, lmax = 10, name = :γ, kwargs...)
 
 Create a particle species for photons in the spacetime with metric `g`.
 """
-function photons(g; polarization = true, lmax = 6, name = :γ, kwargs...)
+function photons(g; polarization = true, lmax = 10, name = :γ, kwargs...)
     lmax >= 3 || error("Need lmax >= 3")
     description = "Photons"
     γ = radiation(g; adiabatic = true, name, description, kwargs...) |> background |> complete # prevent namespacing in extension below