From 2afd1adec2c3fc09566224b0a0e6cd4113ac120a Mon Sep 17 00:00:00 2001 From: Karthik Chandrashekara Date: Wed, 18 Sep 2024 18:34:27 +0200 Subject: [PATCH] Added power/trap depth label to the Rabi oscillation plot. --- EstimatesForAccordionLattice.m | 25 ++++++++++++++----------- 1 file changed, 14 insertions(+), 11 deletions(-) diff --git a/EstimatesForAccordionLattice.m b/EstimatesForAccordionLattice.m index 631272a..33dfbaa 100644 --- a/EstimatesForAccordionLattice.m +++ b/EstimatesForAccordionLattice.m @@ -164,15 +164,16 @@ grid on %% Rabi oscillations in the first order in Kapitza-Dirac scattering -a = 180 * (AtomicUnitOfPolarizability / (2 * SpeedOfLight * VacuumPermittivity)); -waist_y = 250E-6; -waist_z = 50E-6; -Power = 4; -TrapDepth = ((8 * a .* Power) ./ (pi * waist_y * waist_z)); -TwoPhotonRecoilEnergy = (2*PlanckConstantReduced*2*pi/Wavelength)^2 / (2 * Dy164Mass); -RabiOscillationFrequency = (1/PlanckConstantReduced) .* (sqrt(TrapDepth.^2/2 + TwoPhotonRecoilEnergy^2)); +a = 180 * (AtomicUnitOfPolarizability / (2 * SpeedOfLight * VacuumPermittivity)); +waist_y = 250E-6; +waist_z = 50E-6; +Power = 0.2; +TrapDepth = ((8 * a .* Power) ./ (pi * waist_y * waist_z)); +TwoPhotonRecoilEnergy = (2*PlanckConstantReduced*2*pi/Wavelength)^2 / (2 * Dy164Mass); +TrapDepthsInUnitsOfRecoilEnergy = TrapDepth ./ TwoPhotonRecoilEnergy; +RabiOscillationFrequency = (1/PlanckConstantReduced) .* (sqrt(TrapDepth.^2/2 + TwoPhotonRecoilEnergy^2)); -PulseDurations = linspace(1E-6, 10E-6, 1000); +PulseDurations = linspace(1E-6, 150E-6, 1000); alpha = TwoPhotonRecoilEnergy / PlanckConstantReduced; beta = TrapDepth / PlanckConstantReduced; @@ -182,8 +183,10 @@ RabiOscillations = C .* (sin(0.5 .* PulseDurations .* (RabiOscillationF figure(6); set(gcf,'Position',[100 100 950 750]) -plot(PulseDurations .* 1E6, RabiOscillations, LineWidth=2.0) +plot(PulseDurations .* 1E6, RabiOscillations, LineWidth=2.0, DisplayName=['\bf Power = ' num2str(Power) ' W / Trap depth = ' num2str(round(TrapDepthsInUnitsOfRecoilEnergy, 1)) ' E_r']) xlabel('Pulse duration (µs)', FontSize=16) ylabel('Fraction of atoms in first order', FontSize=16) -% title(['\bf Upper bound = ' num2str(round(max(RabiOscillationFrequency .* 1E-3),1)) ' kHz; Lower bound = ' num2str(round(min(RabiOscillationFrequency .* 1E-3),1)) ' kHz'], FontSize=16) -grid on \ No newline at end of file +title('\bf Expected Rabi oscillation', FontSize=16) +grid on +legend(FontSize=16) +