From 5e7940a1b1d061ac5c2379f212dcd70a8101fb34 Mon Sep 17 00:00:00 2001 From: Karthik Chandrashekara Date: Thu, 23 Jan 2025 01:19:31 +0100 Subject: [PATCH] Latest fully functional code - successfully converged to expected ground states. --- Dipolar-Gas-Simulator/+Plotter/plotLive2D.m | 5 +++-- .../+Plotter/visualizeGSWavefunction2D.m | 4 ++-- Dipolar-Gas-Simulator/+Scripts/run_locally.m | 2 +- .../+Scripts/run_on_cluster.m | 18 +++++++++--------- 4 files changed, 15 insertions(+), 14 deletions(-) diff --git a/Dipolar-Gas-Simulator/+Plotter/plotLive2D.m b/Dipolar-Gas-Simulator/+Plotter/plotLive2D.m index f5ebd1f..6be5442 100644 --- a/Dipolar-Gas-Simulator/+Plotter/plotLive2D.m +++ b/Dipolar-Gas-Simulator/+Plotter/plotLive2D.m @@ -25,7 +25,8 @@ function plotLive2D(psi, Params, Transf, Observ, vrun) set(plotxy, 'EdgeColor', 'none'); cbar1 = colorbar; cbar1.Label.Interpreter = 'latex'; - colormap('turbo') + cbar1.Ticks = []; % Disable the ticks + colormap(gca, 'parula') % clim(ax1,[0.00,0.3]); ylabel(cbar1,'$na_{dd}^2$','FontSize',16,'Rotation',270) xlabel('$x$ ($\mu$m)', 'Interpreter', 'latex', 'FontSize', 14) @@ -38,7 +39,7 @@ function plotLive2D(psi, Params, Transf, Observ, vrun) set(plotxy, 'EdgeColor', 'none'); cbar2 = colorbar; cbar2.Label.Interpreter = 'latex'; - colormap('hsv') + colormap(gca, 'hsv') clim([-pi,pi]) ylabel(cbar2,'$\phi$','FontSize',16,'Rotation',270) xlabel('$x$ ($\mu$m)', 'Interpreter', 'latex', 'FontSize', 14) diff --git a/Dipolar-Gas-Simulator/+Plotter/visualizeGSWavefunction2D.m b/Dipolar-Gas-Simulator/+Plotter/visualizeGSWavefunction2D.m index 629eb4e..3d345b9 100644 --- a/Dipolar-Gas-Simulator/+Plotter/visualizeGSWavefunction2D.m +++ b/Dipolar-Gas-Simulator/+Plotter/visualizeGSWavefunction2D.m @@ -45,7 +45,7 @@ function visualizeGSWavefunction2D(folder_path, run_index) set(plotxy, 'EdgeColor', 'none'); cbar1 = colorbar; cbar1.Label.Interpreter = 'latex'; - colormap('turbo') + colormap(gca, 'parula') % clim(ax1,[0.00,0.3]); ylabel(cbar1,'$na_{dd}^2$','FontSize',16,'Rotation',270) xlabel('$x$ ($\mu$m)', 'Interpreter', 'latex', 'FontSize', 14) @@ -58,7 +58,7 @@ function visualizeGSWavefunction2D(folder_path, run_index) set(plotxy, 'EdgeColor', 'none'); cbar2 = colorbar; cbar2.Label.Interpreter = 'latex'; - colormap('hsv') + colormap(gca, 'hsv') clim([-pi,pi]) ylabel(cbar2,'$\phi$','FontSize',16,'Rotation',270) xlabel('$x$ ($\mu$m)', 'Interpreter', 'latex', 'FontSize', 14) diff --git a/Dipolar-Gas-Simulator/+Scripts/run_locally.m b/Dipolar-Gas-Simulator/+Scripts/run_locally.m index 2344383..9cede11 100644 --- a/Dipolar-Gas-Simulator/+Scripts/run_locally.m +++ b/Dipolar-Gas-Simulator/+Scripts/run_locally.m @@ -182,7 +182,7 @@ OptionsStruct.ScatteringLength = 101.35; OptionsStruct.TrapFrequencies = [0, 0, 72.4]; OptionsStruct.TrapPotentialType = 'None'; -OptionsStruct.NumberOfGridPoints = [256, 256]; +OptionsStruct.NumberOfGridPoints = [128, 128]; OptionsStruct.Dimensions = [100, 100]; OptionsStruct.TimeStepSize = 0.005; % in s OptionsStruct.MinimumTimeStepSize = 1E-5; % in s diff --git a/Dipolar-Gas-Simulator/+Scripts/run_on_cluster.m b/Dipolar-Gas-Simulator/+Scripts/run_on_cluster.m index 85fb5f2..dd4140c 100644 --- a/Dipolar-Gas-Simulator/+Scripts/run_on_cluster.m +++ b/Dipolar-Gas-Simulator/+Scripts/run_on_cluster.m @@ -9,7 +9,7 @@ % as = ((as/add)*Params.add)/Params.a0 % Critical point: 102.5133; Triangular phase: 98.0676; Stripe phase: 100.0289; Honeycomb phase: 101.9903 -%{ + %% - Create Variational2D and Calculator object with specified options OptionsStruct = struct; @@ -22,9 +22,9 @@ OptionsStruct.ScatteringLength = 98.0676; OptionsStruct.TrapFrequencies = [0, 0, 72.4]; OptionsStruct.TrapPotentialType = 'None'; -OptionsStruct.NumberOfGridPoints = [256, 256]; +OptionsStruct.NumberOfGridPoints = [128, 128]; OptionsStruct.Dimensions = [100, 100]; -OptionsStruct.TimeStepSize = 500E-6; % in s +OptionsStruct.TimeStepSize = 0.005; % in s OptionsStruct.MinimumTimeStepSize = 1E-5; % in s OptionsStruct.TimeCutOff = 2E6; % in s OptionsStruct.EnergyTolerance = 5E-10; @@ -64,9 +64,9 @@ OptionsStruct.ScatteringLength = 100.0289; OptionsStruct.TrapFrequencies = [0, 0, 72.4]; OptionsStruct.TrapPotentialType = 'None'; -OptionsStruct.NumberOfGridPoints = [256, 256]; +OptionsStruct.NumberOfGridPoints = [128, 128]; OptionsStruct.Dimensions = [100, 100]; -OptionsStruct.TimeStepSize = 500E-6; % in s +OptionsStruct.TimeStepSize = 0.005; % in s OptionsStruct.MinimumTimeStepSize = 1E-5; % in s OptionsStruct.TimeCutOff = 2E6; % in s OptionsStruct.EnergyTolerance = 5E-10; @@ -93,20 +93,20 @@ solver.Potential = pot.trap(); %-% Run Solver %-% [Params, Transf, psi, V, VDk] = solver.run(); -%} + %% - Create Variational2D and Calculator object with specified options OptionsStruct = struct; -OptionsStruct.NumberOfAtoms = 4.148e+07; +OptionsStruct.NumberOfAtoms = 4.0102e+07; OptionsStruct.DipolarPolarAngle = 0; OptionsStruct.DipolarAzimuthAngle = 0; -OptionsStruct.ScatteringLength = 101.35; +OptionsStruct.ScatteringLength = 101.9903; OptionsStruct.TrapFrequencies = [0, 0, 72.4]; OptionsStruct.TrapPotentialType = 'None'; -OptionsStruct.NumberOfGridPoints = [256, 256]; +OptionsStruct.NumberOfGridPoints = [128, 128]; OptionsStruct.Dimensions = [100, 100]; OptionsStruct.TimeStepSize = 0.005; % in s OptionsStruct.MinimumTimeStepSize = 1E-5; % in s