From 551c16dd848779d161a00c83c122bb7e7dfe0d21 Mon Sep 17 00:00:00 2001 From: Karthik Chandrashekara Date: Mon, 18 Nov 2024 13:43:12 +0100 Subject: [PATCH] Cosmetic changes. --- Dipolar-Gas-Simulator/+Scripts/run_locally.m | 8 +- .../@DipolarGas/setupParameters.m | 74 +++++++++---------- 2 files changed, 41 insertions(+), 41 deletions(-) diff --git a/Dipolar-Gas-Simulator/+Scripts/run_locally.m b/Dipolar-Gas-Simulator/+Scripts/run_locally.m index 7eb96d0..2d18a99 100644 --- a/Dipolar-Gas-Simulator/+Scripts/run_locally.m +++ b/Dipolar-Gas-Simulator/+Scripts/run_locally.m @@ -55,15 +55,15 @@ OptionsStruct = struct; OptionsStruct.NumberOfAtoms = 1E5; OptionsStruct.DipolarPolarAngle = 0; OptionsStruct.DipolarAzimuthAngle = 0; -OptionsStruct.ScatteringLength = 98.0676; % Critical point: 102.515; Triangular phase: 98.0676; Stripe phase: 102.2518; Honeycomb phase: 102.6441 +OptionsStruct.ScatteringLength = 102.2518; % Critical point: 102.515; Triangular phase: 98.0676; Stripe phase: 102.2518; Honeycomb phase: 102.6441 OptionsStruct.TrapFrequencies = [10, 10, 72.4]; OptionsStruct.TrapPotentialType = 'None'; OptionsStruct.NumberOfGridPoints = [128, 128]; -OptionsStruct.Dimensions = [7.5, 7.5]; % 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 = 5E-6; % in s -OptionsStruct.TimeCutOff = 1E5; % in s +OptionsStruct.Dimensions = [6.972, 6.972]; % 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 = 5E-6; % in s +OptionsStruct.TimeCutOff = 1E5; % in s OptionsStruct.EnergyTolerance = 5E-10; OptionsStruct.ResidualTolerance = 1E-05; diff --git a/Dipolar-Gas-Simulator/+VariationalSolver2D/@DipolarGas/setupParameters.m b/Dipolar-Gas-Simulator/+VariationalSolver2D/@DipolarGas/setupParameters.m index e7fd853..4d6226f 100644 --- a/Dipolar-Gas-Simulator/+VariationalSolver2D/@DipolarGas/setupParameters.m +++ b/Dipolar-Gas-Simulator/+VariationalSolver2D/@DipolarGas/setupParameters.m @@ -1,43 +1,43 @@ function [Params] = setupParameters(this) - CONSTANTS = Helper.PhysicsConstants; - hbar = CONSTANTS.PlanckConstantReduced; % [J.s] - kbol = CONSTANTS.BoltzmannConstant; % [J/K] - mu0 = CONSTANTS.VacuumPermeability; % [N/A^2] - muB = CONSTANTS.BohrMagneton; % [J/T] - a0 = CONSTANTS.BohrRadius; % [m] - m0 = CONSTANTS.AtomicMassUnit; % [kg] - w0 = 2*pi*61.6582; % Angular frequency unit [s^-1] - mu0factor = 0.3049584233607396; % =(m0/me)*pi*alpha^2 -- me=mass of electron, alpha=fine struct. const. + CONSTANTS = Helper.PhysicsConstants; + hbar = CONSTANTS.PlanckConstantReduced; % [J.s] + kbol = CONSTANTS.BoltzmannConstant; % [J/K] + mu0 = CONSTANTS.VacuumPermeability; % [N/A^2] + muB = CONSTANTS.BohrMagneton; % [J/T] + a0 = CONSTANTS.BohrRadius; % [m] + m0 = CONSTANTS.AtomicMassUnit; % [kg] + w0 = 2*pi*61.6582; % Angular frequency unit [s^-1] + mu0factor = 0.3049584233607396; % =(m0/me)*pi*alpha^2 -- me=mass of electron, alpha=fine struct. const. % mu0=mu0factor *hbar^2*a0/(m0*muB^2) % Number of points in each direction - Params.Nx = this.NumberOfGridPoints(1); - Params.Ny = this.NumberOfGridPoints(2); + Params.Nx = this.NumberOfGridPoints(1); + Params.Ny = this.NumberOfGridPoints(2); % Dimensions (in units of l0) - Params.Lx = this.Dimensions(1); - Params.Ly = this.Dimensions(2); + Params.Lx = this.Dimensions(1); + Params.Ly = this.Dimensions(2); % Mass, length scale - Params.m = CONSTANTS.Dy164Mass; - l0 = sqrt(hbar/(Params.m*w0)); % Defining a harmonic oscillator length + Params.m = CONSTANTS.Dy164Mass; + l0 = sqrt(hbar/(Params.m*w0)); % Defining a harmonic oscillator length % Atom numbers - Params.N = this.NumberOfAtoms; + Params.N = this.NumberOfAtoms; % Dipole angle - Params.theta = this.DipolarPolarAngle; % pi/2 dipoles along x, theta=0 dipoles along z - Params.phi = this.DipolarAzimuthAngle; + Params.theta = this.DipolarPolarAngle; % pi/2 dipoles along x, theta=0 dipoles along z + Params.phi = this.DipolarAzimuthAngle; % Dipole lengths (units of muB) - Params.mu = CONSTANTS.DyMagneticMoment; + Params.mu = CONSTANTS.DyMagneticMoment; % Scattering lengths - Params.as = this.ScatteringLength*a0; + Params.as = this.ScatteringLength*a0; % Trapping frequencies - Params.wx = 2*pi*this.TrapFrequencies(1); - Params.wy = 2*pi*this.TrapFrequencies(2); - Params.wz = 2*pi*this.TrapFrequencies(3); + Params.wx = 2*pi*this.TrapFrequencies(1); + Params.wy = 2*pi*this.TrapFrequencies(2); + Params.wz = 2*pi*this.TrapFrequencies(3); % Tolerances Params.Etol = this.EnergyTolerance; @@ -64,18 +64,18 @@ function [Params] = setupParameters(this) % ================ Parameters defined by those above ================ % % Contact interaction strength (units of l0/m) - Params.gs = 4*pi*Params.as/l0; + Params.gs = 4*pi*Params.as/l0; % Dipole lengths - Params.add = mu0*Params.mu^2*Params.m/(12*pi*hbar^2); + Params.add = mu0*Params.mu^2*Params.m/(12*pi*hbar^2); % DDI strength - Params.gdd = 12*pi*Params.add/l0; %sometimes the 12 is a 4 --> depends on how Vdk (DDI) is defined + Params.gdd = 12*pi*Params.add/l0; %sometimes the 12 is a 4 --> depends on how Vdk (DDI) is defined % Trap gamma - Params.gx = (Params.wx/w0)^2; - Params.gy = (Params.wy/w0)^2; - Params.gz = (Params.wz/w0)^2; + Params.gx = (Params.wx/w0)^2; + Params.gy = (Params.wy/w0)^2; + Params.gz = (Params.wz/w0)^2; % == Calculate LHY correction to account for quantum fluctuations == % @@ -90,14 +90,14 @@ function [Params] = setupParameters(this) Q5 = real(Q5); end - Params.gammaQF = 128/3*sqrt(pi*(Params.as/l0)^5)*Q5; + Params.gammaQF = 128/3*sqrt(pi*(Params.as/l0)^5)*Q5; % Loading the rest into Params - Params.hbar = hbar; - Params.kbol = kbol; - Params.mu0 = mu0; - Params.muB = muB; - Params.a0 = a0; - Params.w0 = w0; - Params.l0 = l0; + Params.hbar = hbar; + Params.kbol = kbol; + Params.mu0 = mu0; + Params.muB = muB; + Params.a0 = a0; + Params.w0 = w0; + Params.l0 = l0; end \ No newline at end of file