diff --git a/doc/src/Commands_compute.txt b/doc/src/Commands_compute.txt
index f566702609..0817f44adf 100644
--- a/doc/src/Commands_compute.txt
+++ b/doc/src/Commands_compute.txt
@@ -80,6 +80,7 @@ KOKKOS, o = USER-OMP, t = OPT.
 "meso/e/atom"_compute_meso_e_atom.html,
 "meso/rho/atom"_compute_meso_rho_atom.html,
 "meso/t/atom"_compute_meso_t_atom.html,
+"momentum"_compute_momentum.html,
 "msd"_compute_msd.html,
 "msd/chunk"_compute_msd_chunk.html,
 "msd/nongauss"_compute_msd_nongauss.html,
diff --git a/doc/src/compute.txt b/doc/src/compute.txt
index 87dbee57d6..12c6b141b7 100644
--- a/doc/src/compute.txt
+++ b/doc/src/compute.txt
@@ -229,6 +229,7 @@ compute"_Commands_compute.html doc page are followed by one or more of
 "meso/e/atom"_compute_meso_e_atom.html -
 "meso/rho/atom"_compute_meso_rho_atom.html -
 "meso/t/atom"_compute_meso_t_atom.html -
+"momentum"_compute_momentum.html - translational momentum
 "msd"_compute_msd.html - mean-squared displacement of group of atoms
 "msd/chunk"_compute_msd_chunk.html - mean-squared displacement for each chunk
 "msd/nongauss"_compute_msd_nongauss.html - MSD and non-Gaussian parameter of group of atoms
diff --git a/doc/src/compute_momentum.txt b/doc/src/compute_momentum.txt
new file mode 100644
index 0000000000..4db171d42d
--- /dev/null
+++ b/doc/src/compute_momentum.txt
@@ -0,0 +1,45 @@
+"LAMMPS WWW Site"_lws - "LAMMPS Documentation"_ld - "LAMMPS Commands"_lc :c
+
+:link(lws,http://lammps.sandia.gov)
+:link(ld,Manual.html)
+:link(lc,Commands_all.html)
+
+:line
+
+compute momentum command :h3
+
+[Syntax:]
+
+compute ID group-ID momentum :pre
+
+ID, group-ID are documented in "compute"_compute.html command
+momentum = style name of this compute command :ul
+
+[Examples:]
+
+compute 1 all momentum :pre
+
+[Description:]
+
+Define a computation that calculates the translational momentum
+of a group of particles.
+
+The momentum of each particles is computed as m v, where m and v are
+the mass and velocity of the particle.
+
+[Output info:]
+
+This compute calculates a global vector (the summed momentum) on
+length 3. This value can be used by any command that uses a global
+vector value from a compute as input. See the "Howto
+output"_Howto_output.html doc page for an overview of LAMMPS output
+options.
+
+The vector value calculated by this compute is "extensive". The vector
+value will be in momentum "units"_units.html.
+
+[Restrictions:] none
+
+[Related commands:]
+
+[Default:] none
diff --git a/doc/src/computes.txt b/doc/src/computes.txt
index 926b8da222..1c914365b1 100644
--- a/doc/src/computes.txt
+++ b/doc/src/computes.txt
@@ -57,6 +57,7 @@ Computes :h1
    compute_meso_e_atom
    compute_meso_rho_atom
    compute_meso_t_atom
+   compute_momentum
    compute_msd
    compute_msd_chunk
    compute_msd_nongauss
diff --git a/doc/src/lammps.book b/doc/src/lammps.book
index 2738c9b051..816acd3521 100644
--- a/doc/src/lammps.book
+++ b/doc/src/lammps.book
@@ -468,6 +468,7 @@ compute_ke_rigid.html
 compute_meso_e_atom.html
 compute_meso_rho_atom.html
 compute_meso_t_atom.html
+compute_momentum.html
 compute_msd.html
 compute_msd_chunk.html
 compute_msd_nongauss.html