diff --git a/openmc/mgxs/library.py b/openmc/mgxs/library.py
index d9518e78d..671b0b8f0 100644
--- a/openmc/mgxs/library.py
+++ b/openmc/mgxs/library.py
@@ -1094,6 +1094,11 @@ class Library:
             xsdata.set_beta_mgxs(mymgxs, xs_type=xs_type, nuclide=[nuclide],
                                  subdomain=subdomain)
 
+        if 'decay-rate' in self.mgxs_types:
+            mymgxs = self.get_mgxs(domain, 'decay-rate')
+            xsdata.set_decay_rate_mgxs(mymgxs, xs_types=xs_type, nuclide=[nuclide],
+                                subdomain=subdomain)
+
         # If multiplicity matrix is available, prefer that
         if 'multiplicity matrix' in self.mgxs_types:
             mymgxs = self.get_mgxs(domain, 'multiplicity matrix')
diff --git a/openmc/mgxs_library.py b/openmc/mgxs_library.py
index c2a4a9db7..e3c830c23 100644
--- a/openmc/mgxs_library.py
+++ b/openmc/mgxs_library.py
@@ -2044,7 +2044,7 @@ class XSdata:
                     xs_grp.create_dataset("beta", data=self._beta[i])
 
                 if self._decay_rate[i] is not None:
-                    xs_grp.create_dataset("decay rate",
+                    xs_grp.create_dataset("decay-rate",
                                           data=self._decay_rate[i])
 
             if self._scatter_matrix[i] is None:
@@ -2210,7 +2210,7 @@ class XSdata:
             xs_types = ['total', 'absorption', 'fission', 'kappa-fission',
                         'chi', 'chi-prompt', 'chi-delayed', 'nu-fission',
                         'prompt-nu-fission', 'delayed-nu-fission', 'beta',
-                        'decay rate', 'inverse-velocity']
+                        'decay-rate', 'inverse-velocity']
 
             temperature_group = group[temp]
 
diff --git a/openmc/plotter.py b/openmc/plotter.py
index 5da0e81e5..dfa6a13ba 100644
--- a/openmc/plotter.py
+++ b/openmc/plotter.py
@@ -16,7 +16,7 @@ PLOT_TYPES = ['total', 'scatter', 'elastic', 'inelastic', 'fission',
 PLOT_TYPES_MGXS = ['total', 'absorption', 'scatter', 'fission',
                    'kappa-fission', 'nu-fission', 'prompt-nu-fission',
                    'deleyed-nu-fission', 'chi', 'chi-prompt', 'chi-delayed',
-                   'inverse-velocity', 'beta', 'decay rate', 'unity']
+                   'inverse-velocity', 'beta', 'decay-rate', 'unity']
 # Create a dictionary which can be used to convert PLOT_TYPES_MGXS to the
 # openmc.XSdata attribute name needed to access the data
 _PLOT_MGXS_ATTR = {line: line.replace(' ', '_').replace('-', '_')