From 4096e685a1facfdca0125b121d1fdfba34821422 Mon Sep 17 00:00:00 2001 From: Karthik Chandrashekara Date: Sun, 11 Jul 2021 06:13:49 +0200 Subject: [PATCH] Corrected errors in the formula for the acceleration due to the push beam. --- .../@MOTSimulator/accelerationDueToPushBeam.m | 12 ++++++------ 1 file changed, 6 insertions(+), 6 deletions(-) diff --git a/MOT Capture Process Simulation/@MOTSimulator/accelerationDueToPushBeam.m b/MOT Capture Process Simulation/@MOTSimulator/accelerationDueToPushBeam.m index 0472c67..907d702 100644 --- a/MOT Capture Process Simulation/@MOTSimulator/accelerationDueToPushBeam.m +++ b/MOT Capture Process Simulation/@MOTSimulator/accelerationDueToPushBeam.m @@ -4,17 +4,17 @@ function ret = accelerationDueToPushBeam(this, PositionVector, VelocityVector) WaveVectorEndPoint = WaveVectorEndPoint./norm(WaveVectorEndPoint); Origin=[0,0,0]; - SaturationIntensity = this.calculateLocalSaturationIntensity(this.PushBeamSaturationIntensity, PositionVector, Origin, WaveVectorEndPoint, this.PushBeamRadius, this.PushBeamWaist); + SaturationIntensity = this.calculateLocalSaturationIntensity(this.PushBeamSaturationParameter, PositionVector, Origin, WaveVectorEndPoint, this.PushBeamRadius, this.PushBeamWaist); - DopplerShift = (VelocityVector * WaveVectorEndPoint(1:3)') * this.PushBeamWaveVector; + DopplerShift = (VelocityVector * WaveVectorEndPoint(1:3)') * this.PushBeamWaveNumber; Detuning = this.PushBeamDetuning - DopplerShift; - - a_push = (Helper.PhysicsConstants.PlanckConstantReduced * this.PushBeamWaveVector * WaveVectorEndPoint(1:3)/Helper.PhysicsConstants.Dy164Mass).*(this.PushBeamLinewidth * 0.5) .* ... - (SaturationIntensity/(1 + SaturationIntensity + (4 * (Detuning./this.PushBeamLinewidth).^2))); + + s_push = SaturationIntensity/(1 + SaturationIntensity + (4 * (Detuning./this.PushBeamLinewidth).^2)); + a_push = (Helper.PhysicsConstants.PlanckConstantReduced * this.PushBeamWaveNumber * WaveVectorEndPoint(1:3)/Helper.PhysicsConstants.Dy164Mass).*(this.PushBeamLinewidth * 0.5) .* (s_push/(1+s_push)); if this.SpontaneousEmission - a_scatter = this.accelerationDueToSpontaneousEmissionProcess(SaturationIntensity, SaturationIntensity, Detuning, this.PushBeamLinewidth, this.PushBeamWaveVector); + a_scatter = this.accelerationDueToSpontaneousEmissionProcess(s_push, s_push, Detuning, this.PushBeamLinewidth, this.PushBeamWaveNumber); else a_scatter = [0,0,0]; end