diff --git a/src/global.F90 b/src/global.F90
index f00a1d368f..a5720464de 100644
--- a/src/global.F90
+++ b/src/global.F90
@@ -352,6 +352,9 @@ module global
   ! Write out initial source
   logical :: write_initial_source = .false.
 
+  ! Whether create fission neutrons or not. Only applied for MODE_FIXEDSOURCE
+  logical :: create_fission_neutrons = .true.
+
   ! ============================================================================
   ! CMFD VARIABLES
 
diff --git a/src/input_xml.F90 b/src/input_xml.F90
index db08eedaab..c92ceee8c5 100644
--- a/src/input_xml.F90
+++ b/src/input_xml.F90
@@ -1102,6 +1102,19 @@ contains
       end select
     end if
 
+    ! Check whether create fission sites
+    if (run_mode == MODE_FIXEDSOURCE) then
+      if (check_for_node(doc, "create_fission_neutrons")) then
+        call get_node_value(doc, "create_fission_neutrons", temp_str)
+        temp_str = to_lower(temp_str)
+        if (trim(temp_str) == 'true' .or. trim(temp_str) == '1') then
+          create_fission_neutrons = .true.
+        else if (trim(temp_str) == 'false' .or. trim(temp_str) == '0') then
+          create_fission_neutrons = .false.
+        end if
+      end if
+    end if
+
     ! Close settings XML file
     call close_xmldoc(doc)
 
diff --git a/src/physics.F90 b/src/physics.F90
index 6a6a9787d5..232effebf3 100644
--- a/src/physics.F90
+++ b/src/physics.F90
@@ -96,10 +96,11 @@ contains
     ! absorption (including fission)
 
     if (nuc % fissionable) then
-      call sample_fission(i_nuclide, i_reaction)
       if (run_mode == MODE_EIGENVALUE) then
+        call sample_fission(i_nuclide, i_reaction)
         call create_fission_sites(p, i_nuclide, i_reaction, fission_bank, n_bank)
-      elseif (run_mode == MODE_FIXEDSOURCE) then
+      elseif (run_mode == MODE_FIXEDSOURCE .and. create_fission_neutrons) then
+        call sample_fission(i_nuclide, i_reaction)
         call create_fission_sites(p, i_nuclide, i_reaction, &
              p % secondary_bank, p % n_secondary)
       end if
diff --git a/src/physics_mg.F90 b/src/physics_mg.F90
index 2e5e467c14..c60fc6e2e1 100644
--- a/src/physics_mg.F90
+++ b/src/physics_mg.F90
@@ -73,7 +73,7 @@ contains
     if (mat % fissionable) then
       if (run_mode == MODE_EIGENVALUE) then
         call create_fission_sites(p, fission_bank, n_bank)
-      elseif (run_mode == MODE_FIXEDSOURCE) then
+      elseif (run_mode == MODE_FIXEDSOURCE .and. create_fission_neutrons) then
         call create_fission_sites(p, p % secondary_bank, p % n_secondary)
       end if
     end if