diff --git a/VERSION b/VERSION
index 7c66fca579..af9764a59e 100644
--- a/VERSION
+++ b/VERSION
@@ -1 +1 @@
-4.7.1
+4.7.2
diff --git a/doc/changelog.rst b/doc/changelog.rst
index 10914320d5..0158514607 100644
--- a/doc/changelog.rst
+++ b/doc/changelog.rst
@@ -6,6 +6,25 @@ Change Log
 
 .. towncrier release notes start
 
+QuTiP 4.7.2 (2023-06-28)
+========================
+
+This is a bugfix release for QuTiP 4.7.X. It adds support for
+numpy 1.25 and scipy 1.11.
+
+Bug Fixes
+---------
+- Fix setting of sso.m_ops in heterodyne smesolver and passing through of sc_ops to photocurrent solver. (#2081 by Bogdan Reznychenko and Simon Cross)
+- Update calls to SciPy eigvalsh and eigsh to pass the range of eigenvalues to return using ``subset_by_index=``. (#2081 by Simon Cross)
+- Fixed bug where some matrices were wrongly found to be hermitian. (#2082 by AGaliciaMartinez)
+
+Miscellaneous
+-------------
+- Fixed typo in stochastic.py (#2049, by  eltociear)
+- `ptrace` always return density matrix (#2185, issue by udevd)
+- `mesolve` can support mixed callable and Qobj for `e_ops` (#2184 issue by balopat)
+
+
 Version 4.7.1 (December 11, 2022)
 +++++++++++++++++++++++++++++++++
 
diff --git a/doc/changes/2049.misc b/doc/changes/2049.misc
deleted file mode 100644
index bea6f72568..0000000000
--- a/doc/changes/2049.misc
+++ /dev/null
@@ -1 +0,0 @@
-Fixed typo in stochastic.py
diff --git a/doc/changes/2081.bugfix b/doc/changes/2081.bugfix
deleted file mode 100644
index 82501946d9..0000000000
--- a/doc/changes/2081.bugfix
+++ /dev/null
@@ -1,2 +0,0 @@
-Fix setting of sso.m_ops in heterodyne smesolver and passing through of sc_ops to photocurrent solver. (`#2081 <https://github.com/qutip/qutip/pull/2081>`_ by Bogdan Reznychenko and Simon Cross)
-Update calls to SciPy eigvalsh and eigsh to pass the range of eigenvalues to return using ``subset_by_index=`` instead of ``eigvals=``. ``eigvals=`` was deprecated in SciPy 1.11 and will be removed in SciPy 1.12. (`#2081 <https://github.com/qutip/qutip/pull/2081>`_ by Simon Cross)
diff --git a/doc/changes/2082.bugfix b/doc/changes/2082.bugfix
deleted file mode 100644
index deb92dfb7d..0000000000
--- a/doc/changes/2082.bugfix
+++ /dev/null
@@ -1 +0,0 @@
-Fixed bug were some matrices were wrongly found to be hermitian.
diff --git a/doc/changes/2083.misc b/doc/changes/2083.misc
deleted file mode 100644
index 4031e62e82..0000000000
--- a/doc/changes/2083.misc
+++ /dev/null
@@ -1 +0,0 @@
-Fix for scipy 1.11.0
diff --git a/doc/changes/2118.bugfix b/doc/changes/2118.bugfix
deleted file mode 100644
index cbb30fe251..0000000000
--- a/doc/changes/2118.bugfix
+++ /dev/null
@@ -1 +0,0 @@
-mesolve can support mixed callable and Qobj e_ops.
diff --git a/doc/changes/2129.misc b/doc/changes/2129.misc
deleted file mode 100644
index 5102b534f9..0000000000
--- a/doc/changes/2129.misc
+++ /dev/null
@@ -1 +0,0 @@
-ptrace ket always return density matrix
diff --git a/doc/changes/2178.misc b/doc/changes/2178.misc
deleted file mode 100644
index e7489e5c5b..0000000000
--- a/doc/changes/2178.misc
+++ /dev/null
@@ -1 +0,0 @@
-Update for numpy 1.25