diff --git a/Dipolar-Gas-Simulator/+Scripts/run_locally.m b/Dipolar-Gas-Simulator/+Scripts/run_locally.m
index adbbde8..19b5c9c 100644
--- a/Dipolar-Gas-Simulator/+Scripts/run_locally.m
+++ b/Dipolar-Gas-Simulator/+Scripts/run_locally.m
@@ -549,18 +549,54 @@ SaveDirectory = './Results/Data_3D/AnisotropicTrap/TiltedDipoles45';
 JobNumber     = 0;
 Plotter.visualizeGSWavefunction(SaveDirectory, JobNumber)
 %%
+SaveDirectory = './Results/Data_3D/GradientDescent/aS_079_theta_000_phi_000_N_90000';
+JobNumber     = 0;
+Plotter.visualizeGSWavefunction(SaveDirectory, JobNumber)
+%%
 SaveDirectory = './Results/Data_3D/GradientDescent/aS_080_theta_000_phi_000_N_90000';
 JobNumber     = 0;
 Plotter.visualizeGSWavefunction(SaveDirectory, JobNumber)
 %%
+SaveDirectory = './Results/Data_3D/GradientDescent/aS_081_theta_000_phi_000_N_90000';
+JobNumber     = 0;
+Plotter.visualizeGSWavefunction(SaveDirectory, JobNumber)
+%%
+SaveDirectory = './Results/Data_3D/GradientDescent/aS_082_theta_000_phi_000_N_90000';
+JobNumber     = 0;
+Plotter.visualizeGSWavefunction(SaveDirectory, JobNumber)
+%%
+SaveDirectory = './Results/Data_3D/GradientDescent/aS_083_theta_000_phi_000_N_90000';
+JobNumber     = 0;
+Plotter.visualizeGSWavefunction(SaveDirectory, JobNumber)
+%%
+SaveDirectory = './Results/Data_3D/GradientDescent/aS_084_theta_000_phi_000_N_90000';
+JobNumber     = 0;
+Plotter.visualizeGSWavefunction(SaveDirectory, JobNumber)
+%%
 SaveDirectory = './Results/Data_3D/GradientDescent/aS_085_theta_000_phi_000_N_90000';
 JobNumber     = 0;
 Plotter.visualizeGSWavefunction(SaveDirectory, JobNumber)
 %%
+SaveDirectory = './Results/Data_3D/GradientDescent/aS_086_theta_000_phi_000_N_90000';
+JobNumber     = 0;
+Plotter.visualizeGSWavefunction(SaveDirectory, JobNumber)
+%%
+SaveDirectory = './Results/Data_3D/GradientDescent/aS_087_theta_000_phi_000_N_90000';
+JobNumber     = 0;
+Plotter.visualizeGSWavefunction(SaveDirectory, JobNumber)
+%%
+SaveDirectory = './Results/Data_3D/GradientDescent/aS_088_theta_000_phi_000_N_90000';
+JobNumber     = 0;
+Plotter.visualizeGSWavefunction(SaveDirectory, JobNumber)
+%%
+SaveDirectory = './Results/Data_3D/GradientDescent/aS_089_theta_000_phi_000_N_90000';
+JobNumber     = 0;
+Plotter.visualizeGSWavefunction(SaveDirectory, JobNumber)
+%%
 SaveDirectory = './Results/Data_3D/GradientDescent/aS_090_theta_000_phi_000_N_90000';
 JobNumber     = 0;
 Plotter.visualizeGSWavefunction(SaveDirectory, JobNumber)
 %%
 SaveDirectory = './Results/Data_3D/GradientDescent/';
-JobNumber     = 0;
+JobNumber     = 1;
 Plotter.visualizeGSWavefunction(SaveDirectory, JobNumber)
\ No newline at end of file
diff --git a/Dipolar-Gas-Simulator/+Scripts/run_on_cluster.m b/Dipolar-Gas-Simulator/+Scripts/run_on_cluster.m
index 4927b21..a245b07 100644
--- a/Dipolar-Gas-Simulator/+Scripts/run_on_cluster.m
+++ b/Dipolar-Gas-Simulator/+Scripts/run_on_cluster.m
@@ -16,6 +16,11 @@ OptionsStruct.CutoffType              = 'Cylindrical';
 OptionsStruct.SimulationMode          = 'EnergyMinimization';     % 'ImaginaryTimeEvolution' | 'RealTimeEvolution' | 'EnergyMinimization'
 OptionsStruct.GradientDescentMethod   = 'NonLinearCGD';           % 'HeavyBall' | 'NonLinearCGD' 
 OptionsStruct.MaxIterationsForGD      = 15000;
+OptionsStruct.TimeStepSize            = 1e-3;            % in s
+OptionsStruct.MinimumTimeStepSize     = 1E-6;            % in s
+OptionsStruct.TimeCutOff              = 1E5;             % in s
+OptionsStruct.EnergyTolerance         = 5E-10;
+OptionsStruct.ResidualTolerance       = 1E-05;
 OptionsStruct.NoiseScaleFactor        = 0.010;
 
 OptionsStruct.PlotLive                = false;
diff --git a/Dipolar-Gas-Simulator/bwhpc_matlab_gpe_sim_gpu.slurm b/Dipolar-Gas-Simulator/bwhpc_matlab_gpe_sim_gpu.slurm
index 5ada7c5..b2b3b90 100644
--- a/Dipolar-Gas-Simulator/bwhpc_matlab_gpe_sim_gpu.slurm
+++ b/Dipolar-Gas-Simulator/bwhpc_matlab_gpe_sim_gpu.slurm
@@ -9,7 +9,7 @@
 #SBATCH --gres=gpu:1
 #SBATCH --mem=12G 
 # Estimated wallclock time for job
-#SBATCH --time=03:00:00
+#SBATCH --time=00:40:00
 #SBATCH --job-name=simulation 
 #SBATCH --error=simulation.err
 #SBATCH --output=simulation.out 
@@ -27,7 +27,6 @@ echo "Number of GPUs allocated to job:      $SLURM_GPUS"
 # Load module 
 module load math/matlab/R2023a
 
-echo Directory is `pwd`
 echo "Initiating Job..."
 
 # Start a Matlab program
diff --git a/Dipolar-Gas-Simulator/submit_jobs.sh b/Dipolar-Gas-Simulator/submit_jobs.sh
index bf39401..e086fa2 100644
--- a/Dipolar-Gas-Simulator/submit_jobs.sh
+++ b/Dipolar-Gas-Simulator/submit_jobs.sh
@@ -1,11 +1,11 @@
 # ----------- Define scan ranges -----------
 
-# Use space-separated or comma-separated floats/integers (can be bracketed)
-SCATTERING_LENGTH_RANGE="[85.0]"
-POLAR_ANGLE_RANGE="[0]"
+# Use space-separated floating-point/integer values
+SCATTERING_LENGTH_RANGE="[79.0 80.0 81.0 82.0 83.0 84.0 85.0 86.0 87.0 88.0 89.0 90.0]"
+POLAR_ANGLE_RANGE="[0.0]"
 AZIMUTHAL_ANGLE_RANGE="[0.0]"
-NUM_ATOMS_LIST="[90000]"
-CHUNK_SIZE=1
+NUM_ATOMS_LIST="[50000 60000 70000 80000]"
+CHUNK_SIZE=4
 
 # ----------- Count total combinations for SLURM array -----------
 
@@ -32,24 +32,15 @@ sbatch --export=SCATTERING_LENGTH_RANGE="$SCATTERING_LENGTH_RANGE",POLAR_ANGLE_R
 #SBATCH --cpus-per-task=8
 #SBATCH --gres=gpu:1
 #SBATCH --mem=16G
-#SBATCH --time=00:40:00
+#SBATCH --time=03:30:00
 #SBATCH --job-name=simulation
 #SBATCH --error=simulation_%A_%a.err
 #SBATCH --output=simulation_%A_%a.out
 #SBATCH --array=1-${numArrayJobs}
 
-echo "Working Directory:                    $PWD"
-echo "Running on host                       $HOSTNAME"
-echo "Job id:                               $SLURM_JOB_ID"
-echo "Job name:                             $SLURM_JOB_NAME"
-echo "Number of nodes allocated to job:     $SLURM_JOB_NUM_NODES"
-echo "Number of GPUs allocated to job:      $SLURM_GPUS"
-echo "SLURM task index:                     \$SLURM_ARRAY_TASK_ID"
-
 # Load MATLAB
 module load math/matlab/R2023a
 
-echo Directory is pwd
 echo "Initiating Job..."
 
 # Run MATLAB wrapper with this batch index