diff --git a/README.md b/README.md
index 023c05e747cee0978d98490ad52abeafeeffa60f..7be3f592c55edc2089a145e0c079864d41a7b5f0 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)
 
 
@@ -89,6 +90,12 @@ 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)
 
 
 <!-- ACKNOWLEDGEMENTS -->