From 2568b8d678c0624dcdd9c057b9a48fc23c09a350 Mon Sep 17 00:00:00 2001 From: Karthik Chandrashekara Date: Mon, 18 Nov 2024 22:42:52 +0100 Subject: [PATCH] Latest working version - tested on cluster. --- .../+BdGSolver2D/solveBogoliubovdeGennesIn2D.m | 4 ++-- Dipolar-Gas-Simulator/+Scripts/run_locally.m | 8 +++++--- Dipolar-Gas-Simulator/+Scripts/run_on_cluster.m | 6 +++--- 3 files changed, 10 insertions(+), 8 deletions(-) diff --git a/Dipolar-Gas-Simulator/+BdGSolver2D/solveBogoliubovdeGennesIn2D.m b/Dipolar-Gas-Simulator/+BdGSolver2D/solveBogoliubovdeGennesIn2D.m index f09c0c6..44a1430 100644 --- a/Dipolar-Gas-Simulator/+BdGSolver2D/solveBogoliubovdeGennesIn2D.m +++ b/Dipolar-Gas-Simulator/+BdGSolver2D/solveBogoliubovdeGennesIn2D.m @@ -5,11 +5,11 @@ gdd = Params.gdd; gammaQF = Params.gammaQF; KEop = 0.5*(Transf.KX.^2+Transf.KY.^2); -Ez = (0.25*VParams.ell^2) + (0.25*Params.gz*VParams.ell^2); -muchem_tilde = muchem - Ez; g_pf_2D = 1/(sqrt(2*pi)*VParams.ell); gQF_pf_2D = sqrt(2/5)/(pi^(3/4)*VParams.ell^(3/2)); +Ez = (0.25/VParams.ell^2) + (0.25*Params.gz*VParams.ell^2); +muchem_tilde = muchem - Ez; % eigs only works with column vectors psi = psi.'; diff --git a/Dipolar-Gas-Simulator/+Scripts/run_locally.m b/Dipolar-Gas-Simulator/+Scripts/run_locally.m index 66f584d..3c0ffbf 100644 --- a/Dipolar-Gas-Simulator/+Scripts/run_locally.m +++ b/Dipolar-Gas-Simulator/+Scripts/run_locally.m @@ -52,7 +52,7 @@ Plotter.visualizeGSWavefunction(Params.njob) OptionsStruct = struct; -OptionsStruct.NumberOfAtoms = 1.6E8; +OptionsStruct.NumberOfAtoms = 1.634E8; OptionsStruct.DipolarPolarAngle = 0; OptionsStruct.DipolarAzimuthAngle = 0; OptionsStruct.ScatteringLength = 102.2518; % Critical point: 102.515; Triangular phase: 98.0676; Stripe phase: 102.2518; Honeycomb phase: 102.6441 @@ -62,12 +62,12 @@ OptionsStruct.TrapPotentialType = 'None'; OptionsStruct.NumberOfGridPoints = [256, 256]; OptionsStruct.Dimensions = [100, 100]; % Critical point: 6.996; Triangular phase: 7.5; Stripe phase: 6.972; Honeycomb phase: 6.239 for both for Atom Number fixed to 1E5 -OptionsStruct.TimeStepSize = 500E-6; % in s +OptionsStruct.TimeStepSize = 100E-6; % in s OptionsStruct.MinimumTimeStepSize = 1E-5; % in s OptionsStruct.TimeCutOff = 2E6; % in s OptionsStruct.EnergyTolerance = 5E-10; OptionsStruct.ResidualTolerance = 1E-04; -OptionsStruct.NoiseScaleFactor = 8; +OptionsStruct.NoiseScaleFactor = 4; OptionsStruct.MaxIterations = 20; OptionsStruct.VariationalWidth = 4; @@ -99,3 +99,5 @@ Plotter.visualizeWavefunction2D(psi,Params,Transf) %% - Plot GS wavefunction Plotter.visualizeGSWavefunction2D(Params.njob) +%% +Plotter.visualizeGSWavefunction2D(1) diff --git a/Dipolar-Gas-Simulator/+Scripts/run_on_cluster.m b/Dipolar-Gas-Simulator/+Scripts/run_on_cluster.m index a74a7c3..8e8badf 100644 --- a/Dipolar-Gas-Simulator/+Scripts/run_on_cluster.m +++ b/Dipolar-Gas-Simulator/+Scripts/run_on_cluster.m @@ -2,7 +2,7 @@ OptionsStruct = struct; -OptionsStruct.NumberOfAtoms = 1.6E8; +OptionsStruct.NumberOfAtoms = 1.634E8; OptionsStruct.DipolarPolarAngle = 0; OptionsStruct.DipolarAzimuthAngle = 0; OptionsStruct.ScatteringLength = 102.2518; % Critical point: 102.515; Triangular phase: 98.0676; Stripe phase: 102.2518; Honeycomb phase: 102.6441 @@ -12,12 +12,12 @@ OptionsStruct.TrapPotentialType = 'None'; OptionsStruct.NumberOfGridPoints = [256, 256]; OptionsStruct.Dimensions = [100, 100]; % Critical point: 6.996; Triangular phase: 7.5; Stripe phase: 6.972; Honeycomb phase: 6.239 for both for Atom Number fixed to 1E5 -OptionsStruct.TimeStepSize = 500E-6; % in s +OptionsStruct.TimeStepSize = 100E-6; % in s OptionsStruct.MinimumTimeStepSize = 1E-5; % in s OptionsStruct.TimeCutOff = 2E6; % in s OptionsStruct.EnergyTolerance = 5E-10; OptionsStruct.ResidualTolerance = 1E-04; -OptionsStruct.NoiseScaleFactor = 8; +OptionsStruct.NoiseScaleFactor = 4; OptionsStruct.MaxIterations = 20; OptionsStruct.VariationalWidth = 4;