diff --git a/scripts/group/features/multi_gaussian_expansion/simulator.py b/scripts/group/features/multi_gaussian_expansion/simulator.py
index 57d38f8f3..b827a0933 100644
--- a/scripts/group/features/multi_gaussian_expansion/simulator.py
+++ b/scripts/group/features/multi_gaussian_expansion/simulator.py
@@ -226,10 +226,6 @@
 
     Solve for the lensed positions of the source galaxy.
     """
-    import os
-
-    small_datasets = os.environ.pop("PYAUTO_SMALL_DATASETS", None)
-
     solver = al.PointSolver.for_grid(
         grid=al.Grid2D.uniform(shape_native=(500, 500), pixel_scales=0.1),
         pixel_scale_precision=0.001,
@@ -240,9 +236,6 @@
         tracer=tracer, source_plane_coordinate=source_galaxy.bulge.centre
     )
 
-    if small_datasets is not None:
-        os.environ["PYAUTO_SMALL_DATASETS"] = small_datasets
-
     al.output_to_json(
         obj=positions,
         file_path=dataset_path / "positions.json",
diff --git a/scripts/group/features/no_lens_light/simulator.py b/scripts/group/features/no_lens_light/simulator.py
index 2d34099aa..183a8ef41 100644
--- a/scripts/group/features/no_lens_light/simulator.py
+++ b/scripts/group/features/no_lens_light/simulator.py
@@ -254,10 +254,6 @@
 Solve for the lensed positions of the source galaxy, which are used as input for the group
 modeling scripts (e.g. SLaM pipeline) to help the non-linear search converge.
 """
-import os
-
-small_datasets = os.environ.pop("PYAUTO_SMALL_DATASETS", None)
-
 solver = al.PointSolver.for_grid(
     grid=al.Grid2D.uniform(shape_native=(500, 500), pixel_scales=0.1),
     pixel_scale_precision=0.001,
@@ -268,9 +264,6 @@
     tracer=tracer, source_plane_coordinate=source_galaxy.bulge.centre
 )
 
-if small_datasets is not None:
-    os.environ["PYAUTO_SMALL_DATASETS"] = small_datasets
-
 al.output_to_json(
     obj=positions,
     file_path=dataset_path / "positions.json",
diff --git a/scripts/group/simulator.py b/scripts/group/simulator.py
index 0f8c22275..3df992b2e 100644
--- a/scripts/group/simulator.py
+++ b/scripts/group/simulator.py
@@ -291,10 +291,6 @@
 Solve for the lensed positions of the source galaxy, which are used as input for the group
 modeling scripts (e.g. SLaM pipeline) to help the non-linear search converge.
 """
-import os
-
-small_datasets = os.environ.pop("PYAUTO_SMALL_DATASETS", None)
-
 solver = al.PointSolver.for_grid(
     grid=al.Grid2D.uniform(shape_native=(500, 500), pixel_scales=0.1),
     pixel_scale_precision=0.001,
@@ -305,9 +301,6 @@
     tracer=tracer, source_plane_coordinate=source_galaxy.bulge.centre
 )
 
-if small_datasets is not None:
-    os.environ["PYAUTO_SMALL_DATASETS"] = small_datasets
-
 al.output_to_json(
     obj=positions,
     file_path=dataset_path / "positions.json",