diff --git a/conftest.py b/conftest.py
index f877d016af..c39b1d52b7 100644
--- a/conftest.py
+++ b/conftest.py
@@ -9,9 +9,9 @@
 import argparse
 import os
 import re
+import subprocess
 from pathlib import Path
 
-import procrunner
 import pytest
 
 
@@ -50,7 +50,7 @@ def ccp4():
         pytest.skip("CCP4 installation required for this test")
 
     try:
-        result = procrunner.run(["refmac5", "-i"], print_stdout=False)
+        result = subprocess.run(["refmac5", "-i"], capture_output=True)
     except OSError:
         pytest.skip(
             "CCP4 installation required for this test - Could not find CCP4 executable"
@@ -74,7 +74,7 @@ def xds():
     Skip the test if XDS is not installed.
     """
     try:
-        result = procrunner.run(["xds"], print_stdout=False)
+        result = subprocess.run(["xds"], capture_output=True)
     except OSError:
         pytest.skip("XDS installation required for this test")
     version = re.search(rb"BUILT=([0-9]+)\)", result.stdout)
@@ -100,7 +100,7 @@ def ensure_repository_is_clean():
     if not os.getenv("CHECK_CLEAN_WORKDIR"):
         return
     print("Working directory:", _repository)
-    status = procrunner.run(("git", "status", "-s"), working_directory=_repository)
+    status = subprocess.run(("git", "status", "-s"), cwd=_repository)
     if status.stdout:
         assert False, "Working directory is not clean"
 
diff --git a/newsfragments/684.misc b/newsfragments/684.misc
new file mode 100644
index 0000000000..496e026215
--- /dev/null
+++ b/newsfragments/684.misc
@@ -0,0 +1 @@
+Replace ``procrunner.run()`` usage with standard Python library routine ``subprocess.run()``
diff --git a/setup.cfg b/setup.cfg
index d520e15605..4b97dcc990 100644
--- a/setup.cfg
+++ b/setup.cfg
@@ -29,7 +29,6 @@ include_package_data = True
 install_requires =
     dials-data>=2.0
     Jinja2
-    procrunner
     pyyaml
     tabulate
 packages = find:
diff --git a/src/xia2/Handlers/Versions.py b/src/xia2/Handlers/Versions.py
index 4023998466..5c84597f09 100644
--- a/src/xia2/Handlers/Versions.py
+++ b/src/xia2/Handlers/Versions.py
@@ -5,8 +5,6 @@
 import subprocess
 from collections.abc import Mapping
 
-import procrunner
-
 import dials.util.version
 
 import xia2.XIA2Version
@@ -45,8 +43,10 @@ def get_xia2_version():
 @functools.lru_cache(maxsize=None)
 def get_xds_version():
     try:
-        result = procrunner.run(
-            ["xds"], print_stdout=False, print_stderr=False, stdin=subprocess.DEVNULL
+        result = subprocess.run(
+            ["xds"],
+            stdin=subprocess.DEVNULL,
+            capture_output=True,
         )
     except OSError:
         pass
@@ -58,11 +58,10 @@ def get_xds_version():
 
 @functools.lru_cache(maxsize=None)
 def get_aimless_version():
-    result = procrunner.run(
+    result = subprocess.run(
         ["aimless", "--no-input"],
-        print_stdout=False,
-        print_stderr=False,
         stdin=subprocess.DEVNULL,
+        capture_output=True,
     )
     version = re.search(rb"version\s\d+\.\d+\.\d+", result.stdout)
     if version:
@@ -72,8 +71,10 @@ def get_aimless_version():
 
 @functools.lru_cache(maxsize=None)
 def get_pointless_version():
-    result = procrunner.run(
-        ["pointless"], print_stdout=False, print_stderr=False, stdin=subprocess.DEVNULL
+    result = subprocess.run(
+        ["pointless"],
+        stdin=subprocess.DEVNULL,
+        capture_output=True,
     )
     version = re.search(rb"version\s\d+\.\d+\.\d+", result.stdout)
     if version:
diff --git a/src/xia2/Modules/SSX/data_integration_standard.py b/src/xia2/Modules/SSX/data_integration_standard.py
index 20a697a55e..073d6be581 100644
--- a/src/xia2/Modules/SSX/data_integration_standard.py
+++ b/src/xia2/Modules/SSX/data_integration_standard.py
@@ -6,11 +6,11 @@
 import math
 import os
 import pathlib
+import subprocess
 from dataclasses import asdict, dataclass, field
 from typing import List, Optional, Tuple
 
 import numpy as np
-import procrunner
 
 from dials.util.mp import multi_node_parallel_map
 from dxtbx.serialize import load
@@ -278,7 +278,7 @@ def run_import(working_directory: pathlib.Path, file_input: FileInput) -> None:
     else:
         xia2_logger.notice(banner("Importing"))  # type: ignore
     with record_step("dials.import"):
-        result = procrunner.run(import_command, working_directory=working_directory)
+        result = subprocess.run(import_command, cwd=working_directory)
         if result.returncode or result.stderr:
             raise ValueError(
                 "dials.import returned error status:\n" + str(result.stderr)
diff --git a/src/xia2/cli/make_sphinx_html.py b/src/xia2/cli/make_sphinx_html.py
index 8bf211d183..2ebf0badb7 100644
--- a/src/xia2/cli/make_sphinx_html.py
+++ b/src/xia2/cli/make_sphinx_html.py
@@ -4,8 +4,6 @@
 import shutil
 import subprocess
 
-import procrunner
-
 import xia2
 
 
@@ -18,15 +16,11 @@ def run():
     sphinx_dir = xia2_dir / "doc" / "sphinx"
     if dest_dir.is_dir():
         shutil.rmtree(dest_dir)
-    result = procrunner.run(
-        ["make", "clean"], working_directory=sphinx_dir, stdin=subprocess.DEVNULL
-    )
+    result = subprocess.run(["make", "clean"], cwd=sphinx_dir, stdin=subprocess.DEVNULL)
     if result.returncode:
         exit(f"make clean failed with exit code {result.returncode}")
 
-    result = procrunner.run(
-        ["make", "html"], working_directory=sphinx_dir, stdin=subprocess.DEVNULL
-    )
+    result = subprocess.run(["make", "html"], cwd=sphinx_dir, stdin=subprocess.DEVNULL)
     if result.returncode:
         exit(f"make html failed with exit code {result.returncode}")
 
diff --git a/tests/Applications/test_xia2setup.py b/tests/Applications/test_xia2setup.py
index 2a923f661c..5565054932 100644
--- a/tests/Applications/test_xia2setup.py
+++ b/tests/Applications/test_xia2setup.py
@@ -1,6 +1,8 @@
 from __future__ import annotations
 
-import procrunner
+import os
+import subprocess
+
 import pytest
 
 from xia2.Handlers.XInfo import XInfo
@@ -18,13 +20,13 @@ def insulin_with_missing_image(dials_data, tmp_path):
 
 
 def test_write_xinfo_insulin_with_missing_image(insulin_with_missing_image, tmp_path):
-    result = procrunner.run(
+    result = subprocess.run(
         [
             "xia2.setup",
             f"image={insulin_with_missing_image.parent.joinpath('insulin_1_001.img')}",
         ],
-        environment_override={"CCP4": tmp_path},
-        working_directory=tmp_path,
+        env={"CCP4": tmp_path, **os.environ},
+        cwd=tmp_path,
     )
     assert not result.returncode
     assert not result.stderr
@@ -37,14 +39,14 @@ def test_write_xinfo_insulin_with_missing_image(insulin_with_missing_image, tmp_
 
 
 def test_write_xinfo_template_missing_images(insulin_with_missing_image, tmp_path):
-    result = procrunner.run(
+    result = subprocess.run(
         [
             "xia2.setup",
             f"image={insulin_with_missing_image.parent.joinpath('insulin_1_001.img:1:22')}",
             "read_all_image_headers=False",
         ],
-        environment_override={"CCP4": tmp_path},
-        working_directory=tmp_path,
+        env={"CCP4": tmp_path, **os.environ},
+        cwd=tmp_path,
     )
     assert not result.returncode
     assert not result.stderr
@@ -56,15 +58,15 @@ def test_write_xinfo_template_missing_images(insulin_with_missing_image, tmp_pat
 
 
 def test_write_xinfo_split_sweep(dials_data, tmp_path):
-    result = procrunner.run(
+    result = subprocess.run(
         [
             "xia2.setup",
             f"image={dials_data('insulin', pathlib=True) / 'insulin_1_001.img:1:22'}",
             f"image={dials_data('insulin', pathlib=True) / 'insulin_1_001.img:23:45'}",
             "read_all_image_headers=False",
         ],
-        environment_override={"CCP4": tmp_path},
-        working_directory=tmp_path,
+        env={"CCP4": tmp_path, **os.environ},
+        cwd=tmp_path,
     )
     assert not result.returncode
     assert not result.stderr
@@ -78,14 +80,14 @@ def test_write_xinfo_split_sweep(dials_data, tmp_path):
 
 def test_write_xinfo_unroll(dials_data, tmp_path):
     # This test partially exercises the fix to https://github.com/xia2/xia2/issues/498 with a different syntax
-    result = procrunner.run(
+    result = subprocess.run(
         [
             "xia2.setup",
             f"image={dials_data('insulin', pathlib=True) / 'insulin_1_001.img:1:45:15'}",
             "read_all_image_headers=False",
         ],
-        environment_override={"CCP4": tmp_path},
-        working_directory=tmp_path,
+        env={"CCP4": tmp_path, **os.environ},
+        cwd=tmp_path,
     )
     assert not result.returncode
     assert not result.stderr
diff --git a/tests/regression/test_X4_wide.py b/tests/regression/test_X4_wide.py
index 208bd93dfd..cadccaf962 100644
--- a/tests/regression/test_X4_wide.py
+++ b/tests/regression/test_X4_wide.py
@@ -1,6 +1,7 @@
 from __future__ import annotations
 
-import procrunner
+import subprocess
+
 import pytest
 
 import iotbx.mtz
@@ -46,9 +47,10 @@ def _split_xinfo(data_dir, tmp_path) -> str:
 
 @pytest.mark.parametrize("pipeline,scaler", (("dials", "xdsa"), ("3dii", "dials")))
 def test_incompatible_pipeline_scaler(pipeline, scaler, tmp_path, ccp4):
-    result = procrunner.run(
+    result = subprocess.run(
         ["xia2", f"pipeline={pipeline}", "nproc=1", f"scaler={scaler}"],
-        working_directory=tmp_path,
+        cwd=tmp_path,
+        capture_output=True,
     )
     assert result.returncode
     assert (
@@ -67,7 +69,7 @@ def test_dials_aimless(regression_test, dials_data, tmp_path, ccp4):
         "truncate=cctbx",
         dials_data("x4wide", pathlib=True),
     ]
-    result = procrunner.run(command_line, working_directory=tmp_path)
+    result = subprocess.run(command_line, cwd=tmp_path, capture_output=True)
     success, issues = xia2.Test.regression.check_result(
         "X4_wide.dials-aimless", result, tmp_path, ccp4, expected_space_group="P41212"
     )
@@ -86,7 +88,7 @@ def test_dials_aimless_with_dials_pipeline(regression_test, dials_data, tmp_path
         "truncate=cctbx",
         dials_data("x4wide", pathlib=True),
     ]
-    result = procrunner.run(command_line, working_directory=tmp_path)
+    result = subprocess.run(command_line, cwd=tmp_path, capture_output=True)
     success, issues = xia2.Test.regression.check_result(
         "X4_wide.dials-aimless", result, tmp_path, ccp4
     )
@@ -106,7 +108,7 @@ def test_dials(regression_test, dials_data, tmp_path, ccp4):
         "crystal=bar",
         dials_data("x4wide", pathlib=True),
     ]
-    result = procrunner.run(command_line, working_directory=tmp_path)
+    result = subprocess.run(command_line, cwd=tmp_path, capture_output=True)
     scaled_expt_file = tmp_path / "DataFiles" / "foo_bar_scaled.expt"
     assert scaled_expt_file.is_file()
     scaled_expt = load.experiment_list(scaled_expt_file)
@@ -153,7 +155,7 @@ def test_dials_aimless_split(regression_test, dials_data, tmp_path, ccp4):
         "trust_beam_centre=True",
         "xinfo=%s" % _split_xinfo(dials_data("x4wide", pathlib=True), tmp_path),
     ]
-    result = procrunner.run(command_line, working_directory=tmp_path)
+    result = subprocess.run(command_line, cwd=tmp_path, capture_output=True)
     success, issues = xia2.Test.regression.check_result(
         "X4_wide_split.dials-aimless", result, tmp_path, ccp4
     )
@@ -170,7 +172,7 @@ def test_dials_split(regression_test, dials_data, tmp_path, ccp4):
         "xinfo=%s" % _split_xinfo(dials_data("x4wide", pathlib=True), tmp_path),
         "mode=parallel",
     ]
-    result = procrunner.run(command_line, working_directory=tmp_path)
+    result = subprocess.run(command_line, cwd=tmp_path, capture_output=True)
     success, issues = xia2.Test.regression.check_result(
         "X4_wide_split.dials",
         result,
@@ -193,7 +195,7 @@ def test_xds(regression_test, dials_data, tmp_path, ccp4, xds):
         "read_all_image_headers=False",
         dials_data("x4wide", pathlib=True),
     ]
-    result = procrunner.run(command_line, working_directory=tmp_path)
+    result = subprocess.run(command_line, cwd=tmp_path, capture_output=True)
     success, issues = xia2.Test.regression.check_result(
         "X4_wide.xds", result, tmp_path, ccp4, xds, expected_space_group="P41212"
     )
@@ -210,7 +212,7 @@ def test_xds_split(regression_test, dials_data, tmp_path, ccp4, xds):
         "trust_beam_centre=True",
         "xinfo=%s" % _split_xinfo(dials_data("x4wide", pathlib=True), tmp_path),
     ]
-    result = procrunner.run(command_line, working_directory=tmp_path)
+    result = subprocess.run(command_line, cwd=tmp_path, capture_output=True)
     success, issues = xia2.Test.regression.check_result(
         "X4_wide_split.xds", result, tmp_path, ccp4, xds
     )
@@ -226,7 +228,7 @@ def test_xds_ccp4a(regression_test, dials_data, tmp_path, ccp4, xds):
         "trust_beam_centre=True",
         dials_data("x4wide", pathlib=True),
     ]
-    result = procrunner.run(command_line, working_directory=tmp_path)
+    result = subprocess.run(command_line, cwd=tmp_path, capture_output=True)
     success, issues = xia2.Test.regression.check_result(
         "X4_wide.ccp4a", result, tmp_path, ccp4, xds
     )
@@ -245,7 +247,7 @@ def test_xds_ccp4a_split(regression_test, dials_data, tmp_path, ccp4, xds):
         "mode=parallel",
         "xinfo=%s" % _split_xinfo(dials_data("x4wide", pathlib=True), tmp_path),
     ]
-    result = procrunner.run(command_line, working_directory=tmp_path)
+    result = subprocess.run(command_line, cwd=tmp_path, capture_output=True)
     success, issues = xia2.Test.regression.check_result(
         "X4_wide_split.ccp4a", result, tmp_path, ccp4, xds
     )
@@ -269,7 +271,7 @@ def test_space_group_dials(
         "image=%s"
         % dials_data("x4wide", pathlib=True).joinpath("X4_wide_M1S4_2_0001.cbf:20:30"),
     ]
-    result = procrunner.run(command_line, working_directory=tmp_path)
+    result = subprocess.run(command_line, cwd=tmp_path, capture_output=True)
     success, issues = xia2.Test.regression.check_result(
         "X4_wide.space_group.%s" % pipeline,
         result,
@@ -296,7 +298,7 @@ def test_space_group_3dii(
         "image=%s"
         % dials_data("x4wide", pathlib=True).joinpath("X4_wide_M1S4_2_0001.cbf:20:30"),
     ]
-    result = procrunner.run(command_line, working_directory=tmp_path)
+    result = subprocess.run(command_line, cwd=tmp_path, capture_output=True)
     success, issues = xia2.Test.regression.check_result(
         "X4_wide.space_group.3dii",
         result,
diff --git a/tests/regression/test_mad_example.py b/tests/regression/test_mad_example.py
index 014cae8ff0..cfe7de3683 100644
--- a/tests/regression/test_mad_example.py
+++ b/tests/regression/test_mad_example.py
@@ -1,6 +1,6 @@
 from __future__ import annotations
 
-import procrunner
+import subprocess
 
 import xia2.Test.regression
 
@@ -25,7 +25,7 @@ def test_dials(regression_test, dials_data, tmp_path, ccp4):
         "trust_beam_centre=True",
         dials_data("fumarase", pathlib=True),
     ]
-    result = procrunner.run(command_line, working_directory=tmp_path)
+    result = subprocess.run(command_line, cwd=tmp_path)
     success, issues = xia2.Test.regression.check_result(
         "mad_example.dials",
         result,
@@ -46,7 +46,7 @@ def test_dials_aimless(regression_test, dials_data, tmp_path, ccp4):
         "trust_beam_centre=True",
         dials_data("fumarase", pathlib=True),
     ]
-    result = procrunner.run(command_line, working_directory=tmp_path)
+    result = subprocess.run(command_line, cwd=tmp_path)
     success, issues = xia2.Test.regression.check_result(
         "mad_example.dials-aimless",
         result,
@@ -67,7 +67,7 @@ def test_xds(regression_test, dials_data, tmp_path, ccp4, xds):
         "trust_beam_centre=True",
         dials_data("fumarase", pathlib=True),
     ]
-    result = procrunner.run(command_line, working_directory=tmp_path)
+    result = subprocess.run(command_line, cwd=tmp_path)
     success, issues = xia2.Test.regression.check_result(
         "mad_example.xds",
         result,
@@ -90,7 +90,7 @@ def test_xds_ccp4a(regression_test, dials_data, tmp_path, ccp4, xds):
         "scaler=ccp4a",
         dials_data("fumarase", pathlib=True),
     ]
-    result = procrunner.run(command_line, working_directory=tmp_path)
+    result = subprocess.run(command_line, cwd=tmp_path)
     success, issues = xia2.Test.regression.check_result(
         "mad_example.ccp4a",
         result,
diff --git a/tests/regression/test_multiple_sweeps.py b/tests/regression/test_multiple_sweeps.py
index f507a11447..31bb279a26 100644
--- a/tests/regression/test_multiple_sweeps.py
+++ b/tests/regression/test_multiple_sweeps.py
@@ -10,7 +10,8 @@
 
 from __future__ import annotations
 
-import procrunner
+import subprocess
+
 import pytest
 
 
@@ -22,7 +23,7 @@ def test_multiple_sweeps(multi_sweep_type, ccp4, dials_data, tmp_path):
     """
     Run xia2 with various different multiple-sweep options.
 
-    Run xia2 using procrunner and look for errors or timeouts.  Run it with a reduced
+    Run xia2 using subprocess.run() and look for errors or timeouts.  Run it with a reduced
     number of reflections per degree required for profile modelling and turn off the
     Xtriage analysis, since we won't have enough reflections for the default settings
     of either to be successful.
@@ -48,9 +49,10 @@ def test_multiple_sweeps(multi_sweep_type, ccp4, dials_data, tmp_path):
         # Don't run the Xtriage analysis — we don't have enough reflections overall.
         "xtriage_analysis=False",
     ]
-    result = procrunner.run(
+    result = subprocess.run(
         command + [f"image={str(image)}" for image in images],
-        working_directory=tmp_path,
+        capture_output=True,
+        cwd=tmp_path,
     )
 
     assert not result.returncode
diff --git a/tests/regression/test_overload.py b/tests/regression/test_overload.py
index af91084cf8..9b7ba9dbe0 100644
--- a/tests/regression/test_overload.py
+++ b/tests/regression/test_overload.py
@@ -1,11 +1,13 @@
 from __future__ import annotations
 
-import procrunner
+import subprocess
 
 
 def test(dials_data, tmp_path):
     images = list(dials_data("centroid_test_data", pathlib=True).glob("centroid*.cbf"))
 
-    result = procrunner.run(["xia2.overload"] + images, working_directory=tmp_path)
+    result = subprocess.run(
+        ["xia2.overload"] + images, cwd=tmp_path, capture_output=True
+    )
     assert not result.returncode and not result.stderr
     assert (tmp_path / "overload.json").is_file()
diff --git a/tests/regression/test_small_molecule.py b/tests/regression/test_small_molecule.py
index 53685f68b5..ebe04098f8 100644
--- a/tests/regression/test_small_molecule.py
+++ b/tests/regression/test_small_molecule.py
@@ -1,6 +1,6 @@
 from __future__ import annotations
 
-import procrunner
+import subprocess
 
 import xia2.Test.regression
 
@@ -22,7 +22,7 @@ def test_dials_aimless(regression_test, dials_data, tmp_path, ccp4):
         "trust_beam_centre=True",
         dials_data("small_molecule_example", pathlib=True),
     ]
-    result = procrunner.run(command_line, working_directory=tmp_path)
+    result = subprocess.run(command_line, cwd=tmp_path, capture_output=True)
     success, issues = xia2.Test.regression.check_result(
         "small_molecule.dials-aimless",
         result,
@@ -43,7 +43,7 @@ def test_dials(regression_test, dials_data, tmp_path, ccp4):
         "trust_beam_centre=True",
         dials_data("small_molecule_example", pathlib=True),
     ]
-    result = procrunner.run(command_line, working_directory=tmp_path)
+    result = subprocess.run(command_line, cwd=tmp_path, capture_output=True)
     success, issues = xia2.Test.regression.check_result(
         "small_molecule.dials",
         result,
@@ -64,7 +64,7 @@ def test_xds(regression_test, dials_data, tmp_path, ccp4, xds):
         "trust_beam_centre=True",
         dials_data("small_molecule_example", pathlib=True),
     ]
-    result = procrunner.run(command_line, working_directory=tmp_path)
+    result = subprocess.run(command_line, cwd=tmp_path, capture_output=True)
     success, issues = xia2.Test.regression.check_result(
         "small_molecule.xds",
         result,
@@ -87,7 +87,7 @@ def test_xds_ccp4a(regression_test, dials_data, tmp_path, ccp4, xds):
         "scaler=ccp4a",
         dials_data("small_molecule_example", pathlib=True),
     ]
-    result = procrunner.run(command_line, working_directory=tmp_path)
+    result = subprocess.run(command_line, cwd=tmp_path, capture_output=True)
     success, issues = xia2.Test.regression.check_result(
         "small_molecule.ccp4a",
         result,
diff --git a/tests/regression/test_ssx.py b/tests/regression/test_ssx.py
index 9aedad9049..8bfb45d097 100644
--- a/tests/regression/test_ssx.py
+++ b/tests/regression/test_ssx.py
@@ -2,8 +2,7 @@
 
 import json
 import os
-
-import procrunner
+import subprocess
 
 from dxtbx.serialize import load
 
@@ -28,7 +27,7 @@ def test_import_without_reference_or_crystal_info(dials_data, tmp_path):
     args = ["dev.xia2.ssx"]
     args.append("image=" + os.fspath(ssx / "merlin0047_1700*.cbf"))
 
-    result = procrunner.run(args, working_directory=tmp_path)
+    result = subprocess.run(args, cwd=tmp_path, capture_output=True)
     assert not result.returncode and not result.stderr
 
     # want to test
@@ -61,7 +60,7 @@ def test_geometry_refinement_and_run_with_reference(dials_data, tmp_path):
     ]
     args.append("image=" + os.fspath(ssx / "merlin0047_1700*.cbf"))
 
-    result = procrunner.run(args, working_directory=tmp_path)
+    result = subprocess.run(args, cwd=tmp_path, capture_output=True)
     assert not result.returncode and not result.stderr
 
     # Check that the processing went straight to geometry refinement and then
@@ -95,7 +94,7 @@ def test_geometry_refinement_and_run_with_reference(dials_data, tmp_path):
     args.append(f"reference_geometry={os.fspath(reference)}")
     args.append("enable_live_reporting=True")
 
-    result = procrunner.run(args, working_directory=tmp_path)
+    result = subprocess.run(args, cwd=tmp_path, capture_output=True)
     assert not result.returncode and not result.stderr
 
     assert without_reference_identifiers != with_reference_identifiers
@@ -126,7 +125,7 @@ def test_full_run_without_reference(dials_data, tmp_path):
     ]
     args.append("image=" + os.fspath(ssx / "merlin0047_1700*.cbf"))
 
-    result = procrunner.run(args, working_directory=tmp_path)
+    result = subprocess.run(args, cwd=tmp_path, capture_output=True)
     assert not result.returncode and not result.stderr
 
     # Now check that the processing went straight to geometry refinement
@@ -172,7 +171,7 @@ def test_full_run_without_reference(dials_data, tmp_path):
         "directory=batch_1/",
         "d_min=2.0",
     ]
-    result = procrunner.run(args, working_directory=tmp_path)
+    result = subprocess.run(args, cwd=tmp_path capture_output=True)
     assert not result.returncode and not result.stderr"""
 
 
@@ -189,7 +188,7 @@ def test_stepwise_run_without_reference(dials_data, tmp_path):
     args.append("image=" + os.fspath(ssx / "merlin0047_1700*.cbf"))
     args.append("steps=find_spots")
 
-    result = procrunner.run(args, working_directory=tmp_path)
+    result = subprocess.run(args, cwd=tmp_path, capture_output=True)
     assert not result.returncode and not result.stderr
 
     # Now check that the processing went straight to geometry refinement
@@ -222,18 +221,18 @@ def test_stepwise_run_without_reference(dials_data, tmp_path):
 
     # Now rerun, check nothing further done if forgetting to specify geometry
     args[-1] = "steps=index"
-    result = procrunner.run(args, working_directory=tmp_path)
+    result = subprocess.run(args, cwd=tmp_path, capture_output=True)
     assert not result.returncode and not result.stderr
     check_output(tmp_path, find_spots=True)
 
     # Now specify geometry
     args.append("reference_geometry=geometry_refinement/refined.expt")
-    result = procrunner.run(args, working_directory=tmp_path)
+    result = subprocess.run(args, cwd=tmp_path, capture_output=True)
     assert not result.returncode and not result.stderr
     check_output(tmp_path, find_spots=True, index=True)
 
     # Now specify geometry
     args[-2] = "steps=integrate"
-    result = procrunner.run(args, working_directory=tmp_path)
+    result = subprocess.run(args, cwd=tmp_path, capture_output=True)
     assert not result.returncode and not result.stderr
     check_output(tmp_path, find_spots=True, index=True, integrate=True)
diff --git a/tests/regression/test_vmxi_thaumatin.py b/tests/regression/test_vmxi_thaumatin.py
index 181075ff7e..8b5f59de1d 100644
--- a/tests/regression/test_vmxi_thaumatin.py
+++ b/tests/regression/test_vmxi_thaumatin.py
@@ -1,6 +1,7 @@
 from __future__ import annotations
 
-import procrunner
+import subprocess
+
 import pytest
 
 import xia2.Test.regression
@@ -20,7 +21,7 @@ def test_xia2(pipeline, regression_test, dials_data, tmp_path, ccp4):
         "space_group=P41212",
         f"image={master_h5}",
     ]
-    result = procrunner.run(command_line, working_directory=tmp_path)
+    result = subprocess.run(command_line, cwd=tmp_path, capture_output=True)
     success, issues = xia2.Test.regression.check_result(
         f"vmxi_thaumatin.{pipeline}",
         result,
diff --git a/tests/test_run_xia2.py b/tests/test_run_xia2.py
index 590fb0720b..1a71c2cef1 100644
--- a/tests/test_run_xia2.py
+++ b/tests/test_run_xia2.py
@@ -3,9 +3,9 @@
 
 from __future__ import annotations
 
-import procrunner
+import subprocess
 
 
 def test_start_xia2():
-    result = procrunner.run(["xia2"])
+    result = subprocess.run(["xia2"])
     assert result.returncode == 0