diff --git a/README.md b/README.md
index 1afa1fa..dee95dd 100644
--- a/README.md
+++ b/README.md
@@ -1,7 +1,9 @@
 #sensBVP
 
-This is the public repository for the code sensBVP. Its use is in
-calculating sensitivity coefficients of ignition delay time to rate
+This is the public repository for the code sensBVP as cited in the paper
+[Direct Sensitivity Analysis for ignition delay
+times](https://doi.org/10.1016/j.combustflame.2019.08.007). Its use is
+in calculating the sensitivity of ignition delay time to rate
 coefficients. In sensBVP, this is carried out by transforming the
 initial value problem to a boundary value problem, reducing the time to
 calculate the sensitivity coefficients by at least an order of
@@ -31,3 +33,22 @@ species mass fractions, temperature, and pressure can be found in
 "ignitionSensitivities.dat".
 
 Try the cases in the examples directory.
+
+If you use this code in your work, please cite: 
+
+```
+@article{GURURAJAN2019478,
+title = "Direct sensitivity analysis for ignition delay times",
+journal = "Combustion and Flame",
+volume = "209",
+pages = "478 - 480",
+year = "2019",
+issn = "0010-2180",
+doi = "https://doi.org/10.1016/j.combustflame.2019.08.007",
+url = "http://www.sciencedirect.com/science/article/pii/S0010218019303645",
+author = "Vyaas Gururajan and Fokion N. Egolfopoulos",
+keywords = "Ignition, Sensitivity, Jacobian, Chemical kinetics",
+}
+
+
+```