From 1286ad2014f743b556dc596d103beb7b20789051 Mon Sep 17 00:00:00 2001 From: Karthik Chandrashekara Date: Wed, 18 Sep 2024 17:44:15 +0200 Subject: [PATCH] Added plot of expected Rabi oscillations during K-D scattering for a given trap depth. --- EstimatesForAccordionLattice.m | 26 ++++++++++++++++++++++++++ 1 file changed, 26 insertions(+) diff --git a/EstimatesForAccordionLattice.m b/EstimatesForAccordionLattice.m index 19cac75..631272a 100644 --- a/EstimatesForAccordionLattice.m +++ b/EstimatesForAccordionLattice.m @@ -160,4 +160,30 @@ xlim([0 4]); xlabel('Trap depths (E_r)', FontSize=16) ylabel('Rabi oscillation frequency (kHz)', 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 + +%% 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)); + +PulseDurations = linspace(1E-6, 10E-6, 1000); + +alpha = TwoPhotonRecoilEnergy / PlanckConstantReduced; +beta = TrapDepth / PlanckConstantReduced; +C = beta^2 / ((2*beta^2) + (4*alpha^2)); + +RabiOscillations = C .* (sin(0.5 .* PulseDurations .* (RabiOscillationFrequency))); + +figure(6); +set(gcf,'Position',[100 100 950 750]) +plot(PulseDurations .* 1E6, RabiOscillations, LineWidth=2.0) +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