karthik/Calculations
1
0
Fork 0
Code Issues Pull Requests repo.project_board Releases Wiki Activity

Added plot of expected Rabi oscillations during K-D scattering for a given trap depth.

Browse Source
This commit is contained in:
Karthik 2024-09-18 17:44:15 +02:00
parent 85f4af3e60
commit 1286ad2014
1 changed files with 26 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

26
EstimatesForAccordionLattice.m
View File

@ -160,4 +160,30 @@ xlim([0 4]);
xlabel('Trap depths (E_r)', FontSize=16) xlabel('Trap depths (E_r)', FontSize=16)
ylabel('Rabi oscillation frequency (kHz)', 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) 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 grid on