diff --git a/examples/gjf/molecular_dynamics_results/README.md b/examples/gjf/molecular_dynamics_results/README.md
new file mode 100644
index 0000000000..cd3fe81b53
--- /dev/null
+++ b/examples/gjf/molecular_dynamics_results/README.md
@@ -0,0 +1,5 @@
+# LAMMPS GJF-2GJ THERMOSTAT EXAMPLE
+
+## GJF-2GJ THERMOSTAT
+
+This directory contain a series of graphs, which contain the results from numermous molecular dynamics simulations. All simulations are run in the NVT ensemble. Two systems are reported, guaiacol and argon. The damping parameter and the timestep are varied. Also the temperature is varied for argon. GJF U is the half-step velocity and GJF V is the onsite velocity. GJF U and GJF V represent exactly the same configurational statistics.