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

Minor corrections.

Browse Source
This commit is contained in:
Karthik 2024-09-18 22:24:49 +02:00
parent 2afd1adec2
commit b62203046e
1 changed files with 3 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
EstimatesForAccordionLattice.m
View File

@ -173,20 +173,18 @@ TwoPhotonRecoilEnergy = (2*PlanckConstantReduced*2*pi/Wavelength)^2 /
TrapDepthsInUnitsOfRecoilEnergy = TrapDepth ./ TwoPhotonRecoilEnergy;
RabiOscillationFrequency = (1/PlanckConstantReduced) .* (sqrt(TrapDepth.^2/2 + TwoPhotonRecoilEnergy^2));
PulseDurations = linspace(1E-6, 150E-6, 1000);
alpha = TwoPhotonRecoilEnergy / PlanckConstantReduced;
beta = TrapDepth / PlanckConstantReduced;
C = beta^2 / ((2*beta^2) + (4*alpha^2));
RabiOscillations = C .* (sin(0.5 .* PulseDurations .* (RabiOscillationFrequency)));
PulseDurations = linspace(1E-6, 150E-6, 1000);
PopulationInFirstOrders = C .* (sin(0.5 .* PulseDurations .* (RabiOscillationFrequency)));
figure(6);
set(gcf,'Position',[100 100 950 750])
plot(PulseDurations .* 1E6, RabiOscillations, LineWidth=2.0, DisplayName=['\bf Power = ' num2str(Power) ' W / Trap depth = ' num2str(round(TrapDepthsInUnitsOfRecoilEnergy, 1)) ' E_r'])
plot(PulseDurations .* 1E6, PopulationInFirstOrders, 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 Expected Rabi oscillation', FontSize=16)
grid on
legend(FontSize=16)