diff --git a/doc/src/Howto_rheo.rst b/doc/src/Howto_rheo.rst
index 4d5734271a..2407fd484d 100644
--- a/doc/src/Howto_rheo.rst
+++ b/doc/src/Howto_rheo.rst
@@ -15,8 +15,9 @@ details of the system, or develop new capabilities. For instance, the numerics
 associated with calculating gradients, reproducing kernels, etc. are separated
 into distinct classes to simplify the development of new integration schemes
 which can call these calculations. Additional numerical details can be found in
-:ref:`(Clemmer) <howto_rheo_clemmer>`. Example movies illustrating some of these
-capabilities are found at https://www.lammps.org/movies.html#rheopackage.
+:ref:`(Palermo) <howto_rheo_palermo>` and :ref:`(Clemmer) <howto_rheo_clemmer>`.
+Example movies illustrating some of these capabilities are found at
+https://www.lammps.org/movies.html#rheopackage.
 
 Note, if you simply want to run a traditional SPH simulation, the :ref:`SPH package
 <PKG-SPH>` package is likely better suited for your application. It has fewer advanced
@@ -107,6 +108,10 @@ criteria for creating/deleting a bond or altering force calculations).
 
 ----------
 
+.. _howto_rheo_palermo:
+
+**(Palermo)** Palermo, Wolf, Clemmer, O'Connor, Phys. Fluids, 36, 113337 (2024).
+
 .. _howto_rheo_clemmer:
 
 **(Clemmer)** Clemmer, Pierce, O'Connor, Nevins, Jones, Lechman, Tencer, Appl. Math. Model., 130, 310-326 (2024).
diff --git a/doc/src/fix_rheo.rst b/doc/src/fix_rheo.rst
index eb88ef0536..617aa0eac5 100644
--- a/doc/src/fix_rheo.rst
+++ b/doc/src/fix_rheo.rst
@@ -46,8 +46,10 @@ Description
 .. versionadded:: 29Aug2024
 
 Perform time integration for RHEO particles, updating positions, velocities,
-and densities. For an overview of other features available in the RHEO package,
-see :doc:`the RHEO howto <Howto_rheo>`.
+and densities. For a detailed breakdown of the integration timestep and
+numerical details, see :ref:`(Palermo) <rheo_palermo>`. For an overview
+and list of other features available in the RHEO package, see
+:doc:`the RHEO howto <Howto_rheo>`.
 
 The type of kernel is specified using *kstyle* and the cutoff is *cut*. Four
 kernels are currently available. The *quintic* kernel is a standard quintic
@@ -70,7 +72,12 @@ and velocity of solid particles are alternatively reconstructed for every
 fluid-solid interaction to ensure no-slip and pressure-balanced boundaries.
 This is done by estimating the location of the fluid-solid interface and
 extrapolating fluid particle properties across the interface to calculate a
-temporary apparent density and velocity for a solid particle.
+temporary apparent density and velocity for a solid particle. The numerical
+details are the same as those described in
+:ref:`(Palermo) <howto_rheo_palermo>` except there is an additional
+restriction that the reconstructed solid density cannot be less than the
+equilibrium density. This prevents fluid particles from sticking to solid
+surfaces.
 
 A modified form of Fickian particle shifting can be enabled with the
 *shift* keyword. This effectively shifts particle positions to generate a
@@ -163,6 +170,10 @@ Default
 
 ----------
 
+.. _rheo_palermo:
+
+**(Palermo)** Palermo, Wolf, Clemmer, O'Connor, Phys. Fluids, 36, 113337 (2024).
+
 .. _fix_rheo_hu:
 
 **(Hu)** Hu, and Adams J. Comp. Physics, 213, 844-861 (2006).