diff --git a/README.md b/README.md
index d3a29c5025cda8e46565f42649eb821ad077f21e..f4bd0d7ae2e18fa3e4cf6a1b9b30454130777d1c 100644
--- a/README.md
+++ b/README.md
@@ -24,6 +24,7 @@
   * [Validity ranges](#validity-ranges)
 * [License](#license)
 * [Contact](#contact)
+* [How to cite this work](#cite)
 * [Acknowledgements](#acknowledgements)
 
 
@@ -98,7 +99,13 @@ Due to biases at the boundary of the RMF modeling, induced by features in the in
 Valentina Fioretti (INAF OAS Bologna) - valentina.fioretti<at>inaf.it
 Project Link: [https://ahead.iaps.inaf.it/](https://ahead.iaps.inaf.it/)
 
+<!-- CITE -->
+## How to cite this work
 
+Please make sure to cite the following papers if you use the proton response files:
+
+- [V. Fioretti et al., "Unveiling the origin of XMM-Newton soft proton flares I. Design and validation of a response matrix for proton spectral analysis", A&A, 691, A229 (2024)](https://www.aanda.org/articles/aa/full_html/2024/11/aa51231-24/aa51231-24.html)
+- [T. Mineo et al., "Unveiling the origin of XMM-Newton soft proton flares II. Systematics in the proton spectral analysis", A&A, 691, A230 (2024)](https://www.aanda.org/articles/aa/full_html/2024/11/aa51232-24/aa51232-24.html)
 
 <!-- ACKNOWLEDGEMENTS -->
 ## Acknowledgements