diff --git a/crates/ekore/refs.bib b/crates/ekore/refs.bib
index 8d155ed48..480619d12 100644
--- a/crates/ekore/refs.bib
+++ b/crates/ekore/refs.bib
@@ -119,3 +119,13 @@ @article{Gluck1995yr
     pages = "4775--4786",
     year = "1996"
 }
+@phdthesis{Carrazza2015dea,
+    author = "Carrazza, Stefano",
+    title = "{Parton distribution functions with QED corrections}",
+    eprint = "1509.00209",
+    archivePrefix = "arXiv",
+    primaryClass = "hep-ph",
+    doi = "10.13130/carrazza-stefano_phd2015-07-06",
+    school = "Milan U.",
+    year = "2015"
+}
diff --git a/crates/ekore/src/anomalous_dimensions/unpolarized/spacelike/aem1.rs b/crates/ekore/src/anomalous_dimensions/unpolarized/spacelike/aem1.rs
index 0a5feaa6b..1a9e317ba 100644
--- a/crates/ekore/src/anomalous_dimensions/unpolarized/spacelike/aem1.rs
+++ b/crates/ekore/src/anomalous_dimensions/unpolarized/spacelike/aem1.rs
@@ -9,36 +9,36 @@ use crate::anomalous_dimensions::unpolarized::spacelike::as1;
 
 /// Compute the leading-order photon-quark anomalous dimension.
 ///
-/// Implements Eq. (2.5) of
+/// Implements Eq. (2.5) of [\[Carrazza:2015dea\]][crate::bib::Carrazza2015dea].
 pub fn gamma_phq(c: &mut Cache, nf: u8) -> Complex<f64> {
     as1::gamma_gq(c, nf) / CF
 }
 
 /// Compute the leading-order quark-photon anomalous dimension.
 ///
-/// Implements Eq. (2.5) of
+/// Implements Eq. (2.5) of [\[Carrazza:2015dea\]][crate::bib::Carrazza2015dea].
+/// However, we are adding the $N_C$ and the $2n_f$ factors from $\theta$ inside the
+/// definition of $\gamma_{q \gamma}^{(0)}(N)$.
 pub fn gamma_qph(c: &mut Cache, nf: u8) -> Complex<f64> {
     as1::gamma_qg(c, nf) / TR * (NC as f64)
 }
 
 /// Compute the leading-order photon-photon anomalous dimension.
 ///
-/// Implements Eq. (2.5) of
+/// Implements Eq. (2.5) of [\[Carrazza:2015dea\]][crate::bib::Carrazza2015dea].
 pub fn gamma_phph(_c: &mut Cache, nf: u8) -> Complex<f64> {
     let cc = ChargeCombinations { nf };
     (4.0 / 3.0 * (NC as f64) * ((cc.nu() as f64) * EU2 + (cc.nd() as f64) * ED2)).into()
 }
 
-/// Compute the leading-order non-singlet QED anomalous dimension
+/// Compute the leading-order non-singlet |QED| anomalous dimension.
 ///
-/// Implements Eq. (2.5) of
+/// Implements Eq. (2.5) of [\[Carrazza:2015dea\]][crate::bib::Carrazza2015dea].
 pub fn gamma_ns(c: &mut Cache, nf: u8) -> Complex<f64> {
     as1::gamma_ns(c, nf) / CF
 }
 
-/// Compute the leading-order singlet QED anomalous dimension matrix
-///
-/// Implements Eq. (2.5) of
+/// Compute the leading-order singlet |QED| anomalous dimension matrix.
 pub fn gamma_singlet(c: &mut Cache, nf: u8) -> [[Complex<f64>; 4]; 4] {
     let cc = ChargeCombinations { nf };
 
@@ -74,9 +74,7 @@ pub fn gamma_singlet(c: &mut Cache, nf: u8) -> [[Complex<f64>; 4]; 4] {
     ]
 }
 
-/// Compute the leading-order valence QED anomalous dimension matrix
-///
-/// Implements Eq. (2.5) of
+/// Compute the leading-order valence |QED| anomalous dimension matrix.
 pub fn gamma_valence(c: &mut Cache, nf: u8) -> [[Complex<f64>; 2]; 2] {
     let cc = ChargeCombinations { nf };
 
diff --git a/crates/ekore/src/bib.rs b/crates/ekore/src/bib.rs
index 175efe194..06533f6e2 100644
--- a/crates/ekore/src/bib.rs
+++ b/crates/ekore/src/bib.rs
@@ -1,4 +1,4 @@
-//! List of References (autogenerated on 2025-01-13T14:38:04.013494).
+//! List of References (autogenerated on 2025-01-13T14:46:34.886936).
 
 #[allow(non_snake_case)]
 /// The Three loop splitting functions in QCD: The Nonsinglet case
@@ -119,3 +119,15 @@ pub fn Bierenbaum2009zt() {}
 ///
 /// DOI: [10.1103/PhysRevD.53.4775](https:dx.doi.org/10.1103/PhysRevD.53.4775)
 pub fn Gluck1995yr() {}
+
+#[allow(non_snake_case)]
+/// Parton distribution functions with QED corrections
+///
+/// Carrazza, Stefano
+///
+/// Published as PhD thesis at Milan U. (2015)
+///
+/// e-Print: [1509.00209](https://arxiv.org/abs/1509.00209)
+///
+/// DOI: [10.13130/carrazza-stefano_phd2015-07-06](https:dx.doi.org/10.13130/carrazza-stefano_phd2015-07-06)
+pub fn Carrazza2015dea() {}