diff --git a/nipype/interfaces/afni/model.py b/nipype/interfaces/afni/model.py
index 544a4afa4f..7c7cef4ded 100644
--- a/nipype/interfaces/afni/model.py
+++ b/nipype/interfaces/afni/model.py
@@ -191,6 +191,8 @@ class DeconvolveInputSpec(AFNICommandInputSpec):
 class DeconvolveOutputSpec(AFNICommandOutputSpec):
     out_file = File(desc='output statistics file',
                     exists=True)
+    reml_script = File(desc='Autogenerated script for 3dREML')
+    x1D = File(desc='save out X matrix')
 
 
 class Deconvolve(AFNICommand):
@@ -230,7 +232,7 @@ class Deconvolve(AFNICommand):
 
     _cmd = '3dDeconvolve'
     input_spec = DeconvolveInputSpec
-    output_spec = AFNICommandOutputSpec
+    output_spec = DeconvolveOutputSpec
 
     def _parse_inputs(self, skip=None):
         if skip is None:
@@ -249,5 +251,10 @@ def _list_outputs(self):
             else:
                 outputs['x1D'] = self.inputs.x1D
 
+        _gen_fname_opts = {}
+        _gen_fname_opts['basename'] = self.inputs.out_file
+        _gen_fname_opts['cwd'] = os.getcwd()
+        outputs['reml_script'] = self._gen_fname(suffix='.REML_cmd', **_gen_fname_opts))
+
         outputs['out_file'] = self.inputs.out_file
         return outputs