diff --git a/HEU_MET_FAST_003/HEU_MET_FAST_003.py b/HEU_MET_FAST_003/HEU_MET_FAST_003.py
index 79efe11..85f6a6c 100644
--- a/HEU_MET_FAST_003/HEU_MET_FAST_003.py
+++ b/HEU_MET_FAST_003/HEU_MET_FAST_003.py
@@ -104,7 +104,6 @@
             {"block_ids": [1], "xs": xs_oralloy},
             {"block_ids": [2], "xs": xs_tuballoy},
         ],
-        scattering_order=3,
         options={
             "use_precursors": False,
             "verbose_inner_iterations": True,
diff --git a/OpenSn_Logo_CAD/opensn.py b/OpenSn_Logo_CAD/opensn.py
index 99e79da..fab8d33 100644
--- a/OpenSn_Logo_CAD/opensn.py
+++ b/OpenSn_Logo_CAD/opensn.py
@@ -59,7 +59,6 @@
             {"block_ids": [1], "xs": xs_source},
             {"block_ids": [2], "xs": xs_block},
         ],
-        scattering_order=0,
         volumetric_sources=[src],
     )
     ss_solver = SteadyStateSourceSolver(problem=phys)
diff --git a/Six_1g_spherical_benchmarks/Problem_1.py b/Six_1g_spherical_benchmarks/Problem_1.py
index 6b3ff26..dd01803 100644
--- a/Six_1g_spherical_benchmarks/Problem_1.py
+++ b/Six_1g_spherical_benchmarks/Problem_1.py
@@ -102,7 +102,6 @@
         xs_map=[
             {"block_ids": [1], "xs": xs_mat},
         ],
-        scattering_order=0,
         boundary_conditions=[
             {"name": "sphere_surface", "type": "isotropic", "group_strength": bsrc},
         ],
diff --git a/Six_1g_spherical_benchmarks/Problem_2.py b/Six_1g_spherical_benchmarks/Problem_2.py
index da96c6e..f50bd4e 100644
--- a/Six_1g_spherical_benchmarks/Problem_2.py
+++ b/Six_1g_spherical_benchmarks/Problem_2.py
@@ -66,7 +66,6 @@
         xs_map=[
             {"block_ids": [1, 2], "xs": xs_void},
         ],
-        scattering_order=0,
         volumetric_sources=[mg_src],
     )
 
diff --git a/Six_1g_spherical_benchmarks/Problem_3.py b/Six_1g_spherical_benchmarks/Problem_3.py
index 313927d..4bb9e88 100644
--- a/Six_1g_spherical_benchmarks/Problem_3.py
+++ b/Six_1g_spherical_benchmarks/Problem_3.py
@@ -63,7 +63,6 @@
         xs_map=[
             {"block_ids": [1, 2, 3], "xs": xs_void},
         ],
-        scattering_order=0,
         volumetric_sources=[mg_src],
     )
 
diff --git a/Six_1g_spherical_benchmarks/Problem_4.py b/Six_1g_spherical_benchmarks/Problem_4.py
index 8b1afb6..2c3625d 100644
--- a/Six_1g_spherical_benchmarks/Problem_4.py
+++ b/Six_1g_spherical_benchmarks/Problem_4.py
@@ -66,7 +66,6 @@
             {"block_ids": [1], "xs": xs_mat},
             {"block_ids": [2], "xs": xs_void},
         ],
-        scattering_order=0,
         boundary_conditions=[
             {"name": "sphere_surface", "type": "isotropic", "group_strength": bsrc},
         ],
diff --git a/Six_1g_spherical_benchmarks/Problem_5.py b/Six_1g_spherical_benchmarks/Problem_5.py
index 1d5c18c..842df74 100644
--- a/Six_1g_spherical_benchmarks/Problem_5.py
+++ b/Six_1g_spherical_benchmarks/Problem_5.py
@@ -66,7 +66,6 @@
             {"block_ids": [2], "xs": xs_mat},
             {"block_ids": [1, 3], "xs": xs_void},
         ],
-        scattering_order=0,
         boundary_conditions=[
             {"name": "sphere_surface", "type": "isotropic", "group_strength": bsrc},
         ],
diff --git a/Six_1g_spherical_benchmarks/Problem_6.py b/Six_1g_spherical_benchmarks/Problem_6.py
index 3ef913d..0f43506 100644
--- a/Six_1g_spherical_benchmarks/Problem_6.py
+++ b/Six_1g_spherical_benchmarks/Problem_6.py
@@ -66,7 +66,6 @@
             {"block_ids": [1], "xs": xs_mat},
             {"block_ids": [2], "xs": xs_void},
         ],
-        scattering_order=0,
         volumetric_sources=[mg_src],
     )
 
diff --git a/Urban_Source/urban_source.py b/Urban_Source/urban_source.py
index 2e2afe7..da4239d 100644
--- a/Urban_Source/urban_source.py
+++ b/Urban_Source/urban_source.py
@@ -62,7 +62,6 @@
             {"block_ids": [2], "xs": xs_source},
             {"block_ids": [3], "xs": xs_air},
         ],
-        scattering_order=0,
         volumetric_sources=[src],
     )
     ss_solver = SteadyStateSourceSolver(problem=phys)