diff --git a/NEWS.md b/NEWS.md
index ccecb6f7..0406c3e5 100644
--- a/NEWS.md
+++ b/NEWS.md
@@ -1,5 +1,19 @@
 # Molly.jl release notes
 
+## v0.9.0 - Mar 2022
+
+- The arguments to `forces` and `accelerations` are made consistent across implementations.
+- Centre of mass motion is removed by default during simulation using `remove_CM_motion!`.
+- Coordinates are centred in the simulation box by default during setup.
+- The `Langevin` integrator and `Verlet` integrator are added.
+- The `MorseBond` potential is added.
+- The GB-Neck2 implicit solvent model is added via `ImplicitSolventGBN2`.
+- The `CubicSplineCutoff` is added.
+- The `rmsd` function is added.
+- The AtomsBase.jl interface is made more complete.
+- The progress bar is removed from simulations.
+- The out-of-place neighbor list type `NeighborListVec` is changed.
+
 ## v0.8.0 - Feb 2022
 
 - General interactions are renamed to pairwise interactions throughout to better reflect their nature. The abstract type is now `PairwiseInteraction` and the keyword argument to `System` is now `pairwise_inters`. General interaction now refers to a new type of interaction that takes in the whole system and returns forces for all atoms, allowing interactions such as neural network potentials acting on the whole system. This is available via the keyword argument `general_inters` to `System`.
diff --git a/Project.toml b/Project.toml
index 6f22f3e4..8780a40b 100644
--- a/Project.toml
+++ b/Project.toml
@@ -1,7 +1,7 @@
 name = "Molly"
 uuid = "aa0f7f06-fcc0-5ec4-a7f3-a573f33f9c4c"
 authors = ["Joe G Greener <jgreener@hotmail.co.uk>"]
-version = "0.8.0"
+version = "0.9.0"
 
 [deps]
 AtomsBase = "a963bdd2-2df7-4f54-a1ee-49d51e6be12a"
diff --git a/docs/src/docs.md b/docs/src/docs.md
index cc0601f7..494893f3 100644
--- a/docs/src/docs.md
+++ b/docs/src/docs.md
@@ -626,6 +626,9 @@ function Molly.simulate!(sys,
         # Apply coupling like this
         apply_coupling!(sys, sim, sim.coupling)
 
+        # Remove centre of mass motion like this
+        remove_CM_motion!(sys)
+
         # Find new neighbors like this
         neighbors = find_neighbors(sys, sys.neighbor_finder, neighbors, step_n;
                                    parallel=parallel)