ax.set_color_cycle([cm.jet(k) for k in linspace(0,1,N_DET)]) pilier = OmR()#0.304) pilier.set_params(m=25.e-9, f_mech=4.e6, q=2e6, lambda_nm = 1064., losses=15.e-6, transmission_input=30e-6, transmission_output=1.e-6, length=150e-6, i_incident=1000e-6, temp=1., delta_hz=10, detection_losses=0.0) pilier.freq = linspace(3.99e6,4.01e6, 1000) """ DETS = linspace(0,200e3, N_DET) extent = [-90, 90, min(pilier.freq), max(pilier.freq)] for index, delta_hz in enumerate(DETS): subplot(N_DET,1,index+1) pilier.delta_hz = delta_hz img = array([10*log10(abs(pilier.left.output.spectrum_sym(phi))) for phi in linspace(-pi/2,pi/2)]) img = img.T imshow(img, aspect='auto',origin='lower', extent=extent) cb = colorbar() contour(img, levels=[0.], linewidths=3, colors='w', extent=extent) ylabel('$\Delta/2 \pi$') cb.set_label('Noise (dB)') #plot(pilier.freq, len(pilier.freq)*[1.], ':k')
lambda_nm = 1064., losses=1.e-6, transmission_input=30e-6, transmission_output=10.e-6, length=0.00764, i_incident=1e-3, temp=3.8e-4, delta_hz=1., detection_losses=0.304) o.gamma_m_hz = 2560. tau = 1./(2*pi*1.7e6) o.tau_input = tau/0.4 o.tau_output = tau/0.6 o.n_photons = 1.1e8 o.freq = linspace(1.5e6, 1.56e6, 1000) #o.right.output._get_losses = lambda: if True: figname = 'PRX Regal fig1' close(figname) figure(figname, figsize=(15,10)) dets = (3e3, 6e3, 13e3, 20e3) for index, delta_hz in enumerate(dets): subplot(len(dets),2,2*index+1) o.delta_hz = delta_hz o.n_photons = 1.1e8
o = OmR() o.set_params(m=6.75e-12, f_mech=1e6, q=1000., lambda_nm = 1064., losses=1.e-6, transmission_input=300e-6, length=0.00764, i_incident=1e-3, delta_hz=1.) tau = 1./(2*pi*1.7e6) o.tau_in = tau/0.4 o.tau_ex = tau/0.6 o.n_th = 0 o.freq = linspace(0.99e6, 1.01e6, 400) o.delta_hz = 0 figure('SB asymmetry on resonance', figsize=(10,13)) subplot(211) o.n_th = 0 o.cooperativity = 0.1 title('Cooperativity = 0.1, n_th = 0') plot(o.freq, o.right.output.spectrum_ss('upper'), '-r', label='Upper sideband') plot(o.freq, o.right.output.spectrum_ss('lower'), '-b', label='Lower sideband') xlabel('Frequency (Hz)')
lambda_nm = 1064., losses=1.e-6, transmission_input=30e-6, transmission_output=10.e-6, length=0.00764, i_incident=1e-3, temp=3.8e-4, delta_hz=1., detection_losses=0.304) cav.gamma_m_hz = 2560. tau = 1./(2*pi*1.7e6) cav.tau_input = tau/0.4 cav.tau_output = tau/0.6 cav.n_photons = 1.1e8 cav.freq = linspace(1.5e6, 1.56e6, 1000) figure('mean_fields_omr', figsize=(12.9, 10.9)) last=None for index, ((input_trans, out_trans, loss), title_) in enumerate((((1e-6,5e-6, 0), 'undercoupled'), ((5e-6, 5e-6,0), 'critically coupled'), ((5e-6, 0.5e-6,0), 'over coupled'), ((1e-6, 1e-6, 5e-6), 'loss dominated'))): last = subplot(2, 2 ,index + 1, sharex=last) cav.transmission_input = input_trans cav.transmission_output = out_trans cav.losses = loss title(title_) reflected = [] transmitted = [] deltas = linspace(-25*cav.bandwidth_hz, 25*cav.bandwidth_hz, 300) for delta_hz in deltas:
if current_dir=='': current_dir = '.' ### Arcizet 2006 arcizet = OmR() arcizet.set_params(m=190e-9, f_mech=814e3, q=1.e4, lambda_nm = 1064., losses=1.e-6, transmission_input=2*pi/3e4, transmission_output=1e-6, length=0.0024, i_incident=2e-4, temp=300, delta_hz=1.) arcizet.freq=linspace(813.e3, 816.e3, 1000) kappa_hz = arcizet.kappa_hz name = "Arcizet2006" figure(name) for power in array((1.4,2.5,4.5,6.5,9.5)): shifts = [] cooling = [] phis = linspace(-4,4,100) arcizet.delta_hz = 0 arcizet.i_intra = power for index, ph in enumerate(phis): #(0.03, 0.06, 0.09, 0.11, 0.13)):# arcizet.delta_hz = ph*kappa_hz/2 y = arcizet.right.output.spectrum_sym(pi/2)#dummy_spec() c = curve.Curve()