In [2]:
%pylab inline
Populating the interactive namespace from numpy and matplotlib
WARNING: pylab import has clobbered these variables: ['f']
`%matplotlib` prevents importing * from pylab and numpy
In [3]:
h = 6.62606957e-34 # Planck Contant
hbar = h/(2*pi) # Planck Constant over 2 Pi
heV = 4.136E-15 # h in eV
qe = 1.60217657e-19 # Electric Charge
hbar = h/2/pi
DelAl = 0.180e-3 * qe # V
kB = 1.3806488e-23 # Boltzmann Constant
e0 = 8.854e-12 # Permittivity vacuum, F/m
In [4]:
# 1 : ground
# 2 : botton Cap
# 3 : Top Cap
# 4 : Gate
C12 = 113*1e-15
C13 = 84*1e-15
C14 =58*1e-15
C23 =50.134*1e-15
C24 = 0.607*1e-15
C34 =9.826*1e-15
Cj = 2 * 1e-15 # Junctions
Cg = (C34+C13)*(C24+C12)/(C34+C13+C12+C24) # Gate capacitance
Ct = Cg + C23 + Cj # Total Capacitance
print("Ct = ", Ct/1e-15, 'fF')
Ec = qe**2/(2*Ct)/h
print('Ec =',Ec/1e9,'GHz')
R = 20e3 # Ohms
T = 20e-3 # K
Ej =2 * DelAl*tanh(DelAl/(2*kB*T))/(8*qe**2*R)
#Ej = 2*DelAl/(8*qe*R)
print('Ej =', Ej/1e9,'GHz')
E01 = sqrt(8*Ec*Ej)
print('w01 = ', E01/1e9, 'GHz')
beta = Cg/Ct
Vrms = sqrt(hbar*2*pi*5e9/2/beta)
mat = 1/2*(Ej/8/Ec)**(1/4)/sqrt(2)
g = 2*beta*qe*Vrms*mat/h
print('g =', g,'MHz')
Ct = 103.5206664513361 fF
Ec = 0.1871146133104413 GHz
Ej = 14.0433959785 GHz
w01 = 4.58495331107 GHz
g = 271.34241887 MHz
In [5]:
Vrms,mat
Out[5]:
(1.8267802478099438e-12, 0.61876197725847726)
In [6]:
# Cavity kappa
Q = 20e3
freqCav = 5e9
kappa = freqCav/Q
print('Kappa =', kappa/1e9, ' GHz')
Kappa = 0.00025 GHz
In [7]:
# qubit
Qqubit = 4.3e9/1e6
freqQubit = 4.3e9
gamma = freqQubit/Qqubit
print('Gamma =', gamma/1e9, ' GHz')
Gamma = 0.001 GHz
In [8]:
ls
Device Estimation .ipynb
Program Recover IQ from pulsed data.ipynb
Rabi Fit.ipynb
Simulation cavity qubit filter.ipynb
data/
In [8]:
In [8]:
In [8]:
In [8]:
In [8]:
In [8]:
In [9]:
3.14158*432
Out[9]:
1357.16256
In [9]:
In [17]:
f = 5.05e9
T=30e-3
n = 1/(exp(h*f/(kB*T))-1)
neq=0.5*(1/tanh(beta*h*f/2))-.5
n,neq
Out[17]:
(0.00031016188426027783, 6.0204540298315753e+23)
In [2]:
import qutip.testing as qt
In [3]:
qt.run()
Transform 2-level to eigenbasis and back ... ok
Transform 10-level real-values to eigenbasis and back ... ok
Transform 10-level to eigenbasis and back ... ok
Transform 10-level imag to eigenbasis and back ... ok
Consistency between transformations of kets and density matrices ... ok
Check diagonalization via eigenbasis transformation ... ok
brmesolve: simple qubit ... ok
brmesolve: c_ops alone ... ok
brmesolve: c_ops with a_ops ... ok
brmesolve: harmonic oscillator, zero temperature ... ok
brmesolve: harmonic oscillator, finite temperature ... ok
brmesolve: harmonic oscillator, finite temperature, states ... ok
brmesolve: Jaynes-Cummings model, zero temperature ... ok
correlation: legacy me and es for oscillator in coherent initial state ... ok
correlation: comparing me and es for oscillator in coherent initial state ... ok
correlation: comparing me and mc for driven oscillator in ground state ... ok
correlation: legacy me and es for oscillator in steady-state ... ok
correlation: comparing me and es for oscillator in steady-state ... ok
correlation: legacy spectrum from es and pi methods ... ok
correlation: comparing spectrum from es and fft methods ... ok
correlation: comparing spectrum from es and pi methods ... ok
correlation: comparing TLS emission correlations (str-list td format) ... ok
correlation: comparing TLS emission correlations (fn-list td format) ... ok
correlation: comparing TLS emission correlations (fn td format) ... ok
Counting statistics: current and current noise in a DQD model ... ok
qutip.tests.test_cqed.Testcqed.test_dispersivecqed_ISWAP ... ok
qutip.tests.test_cqed.Testcqed.test_dispersivecqed_combination ... ok
Diagonalization of random two-level system ... ok
Diagonalization of composite systems ... ok
Excitation-number-restricted state-space: full state space ... ok
Excitation-number-restricted state space: single excitations ... ok
Excitation-number-restricted state space: two excitations ... ok
Excitation-number-restricted state space: fock states ... ok
Excitation-number-restricted state space: identity operator ... ok
Excitation-number-restricted state space: thermal density operator (I) ... ok
Excitation-number-restricted state space: thermal density operator (II) ... ok
Entropy: von-Neumann entropy ... ok
Entropy: Linear entropy ... ok
Entropy: Concurrence ... ok
Entropy: Mutual information ... ok
Entropy: Conditional entropy ... ok
Entropy: Entangling power ... ok
qutip.tests.test_expect.TestExpect.testExpectSolverCompatibility ... ok
qutip.tests.test_expect.TestExpect.testOperatorDensityMatrix ... ok
qutip.tests.test_expect.TestExpect.testOperatorKet ... ok
qutip.tests.test_expect.TestExpect.testOperatorListState ... ok
qutip.tests.test_expect.TestExpect.testOperatorListStateList ... ok
qutip.tests.test_expect.TestExpect.testOperatorStateList ... ok
qutip.tests.test_fileio.TestFileIO.testRWComplexDecimal ... ok
qutip.tests.test_fileio.TestFileIO.testRWComplexDefault ... ok
qutip.tests.test_fileio.TestFileIO.testRWComplexExp ... ok
qutip.tests.test_fileio.TestFileIO.testRWRealDecimal ... ok
qutip.tests.test_fileio.TestFileIO.testRWRealDefault ... ok
qutip.tests.test_fileio.TestFileIO.testRWRealExp ... ok
qutip.tests.test_fileio.TestFileIO.testRWSeparatorDetection ... ok
qutip.tests.test_floquet.TestFloquet.testFloquetUnitary ... ok
qutip.tests.test_gates.TestGates.testExpandGate1toN ... ok
qutip.tests.test_gates.TestGates.testExpandGate2toN ... ok
qutip.tests.test_gates.TestGates.testExpandGate2toNSwap ... ok
qutip.tests.test_gates.TestGates.testExpandGate3toN ... ok
qutip.tests.test_gates.TestGates.testExpandGate3toN_permutation ... ok
qutip.tests.test_gates.TestGates.testSwapGate ... ok
Graph: Graph Degree ... ok
Graph: Breadth-First Search ... ok
Graph: Reverse Cuthill-McKee Ordering (simple) ... ok
Graph: Reverse Cuthill-McKee Ordering (Bucky) ... ok
Graph: Reverse Cuthill-McKee Ordering (boost) ... ok
Graph: Reverse Cuthill-McKee Ordering (qutip) ... ok
Graph: Maximum Bipartite Matching ... ok
Graph: Weighted Bipartite Matching ... ok
Monte-carlo: Constant H with no collapse ops (expect) ... ok
Monte-carlo: Constant H with no collapse ops (states) ... ok
Monte-carlo: Constant H with no collapse ops (expect and states) ... ok
Monte-carlo: Constant H (str format) with no collapse ops (expect) ... ok
Monte-carlo: Constant H (func format) with no collapse ops (expect) ... ok
Monte-carlo: Constant H (str format) with no collapse ops (states) ... ok
Monte-carlo: Constant H (func format) with no collapse ops (states) ... ok
Monte-carlo: Check for stored collapse operators and times ... ok
Monte-carlo: Constant H with constant collapse ... ok
Monte-carlo: Constant H with constant collapse (states) ... ok
Monte-carlo: Constant H with constant collapse (states) ... ok
Monte-carlo: Constant H with single collapse operator ... ok
Monte-carlo: Constant H with single expect operator ... ok
Monte-carlo: Collapse terms constant (func format) ... ok
Monte-carlo: Collapse terms constant (str format) ... ok
Monte-carlo: Time-dependent H (func format) ... ok
Monte-carlo: Time-dependent H (str format) ... ok
Monte-carlo: check for correct dtypes (average_states=True) ... ok
Monte-carlo: check for correct dtypes (average_states=False) ... ok
Monte-carlo: check reusing seeds ... ok
Monte-carlo: check not reusing seeds ... ok
Monte-carlo: list of trajectories ... ok
mcsolve_f90: Constant H with no collapse ops (expect) ... ok
mcsolve_f90: Constant H with no collapse ops (states) ... ok
mcsolve_f90: Constant H with constant collapse ... ok
mcsolve_f90: Constant H with single collapse operator ... ok
mcsolve_f90: Constant H with single expect operator ... ok
mcsolve_f90: check for correct dtypes (average_states=True) ... ok
mcsolve_f90: check for correct dtypes (average_states=False) ... ok
qutip.tests.test_mesolve.TestJCModelEvolution.testCase1 ... ok
qutip.tests.test_mesolve.TestJCModelEvolution.testCase2 ... ok
qutip.tests.test_mesolve.TestJCModelEvolution.testCase3 ... ok
qutip.tests.test_mesolve.TestJCModelEvolution.testQubitDynamics1 ... ok
qutip.tests.test_mesolve.TestJCModelEvolution.testQubitDynamics2 ... ok
qutip.tests.test_mesolve.TestMESolveTDDecay.testMETDDecayAsArray ... ok
qutip.tests.test_mesolve.TestMESolveTDDecay.testMETDDecayAsFuncList ... ok
qutip.tests.test_mesolve.TestMESolveTDDecay.testMETDDecayAsPartFuncList ... ok
qutip.tests.test_mesolve.TestMESolveTDDecay.testMETDDecayAsStrList ... ok
qutip.tests.test_mesolve.TestMESolverConstDecay.testMEDecay ... ok
qutip.tests.test_mesolve.TestMESolverConstDecay.testMEDecayAsFuncList ... ok
qutip.tests.test_mesolve.TestMESolverConstDecay.testMEDecayAsStrList ... ok
qutip.tests.test_mesolve.TestMESolverConstDecay.testMEDecaySingleCollapse ... ok
qutip.tests.test_mesolve.TestMESolverConstDecay.testMEDecaySingleExpect ... ok
Metrics: Fidelity / trace distance limiting cases ... ok
Metrics: Fidelity, mixed state inequality ... ok
Metrics: Fidelity, invariance under unitary trans. ... ok
Metrics: Trace dist., invariance under unitary trans. ... ok
Metrics: Trace dist. & Fidelity mixed/mixed inequality ... ok
Metrics: Trace dist. & Fidelity mixed/pure inequality ... ok
Metrics: Check avg gate fidelities for random ... ok
Metrics: Hilbert distance. ... ok
Operator CSR Type: commutator ... ok
Operator CSR Type: create ... ok
Operator CSR Type: destroy ... ok
Operator CSR Type: displace ... ok
Operator CSR Type: spin ops ... ok
Operator CSR Type: momentum ... ok
Operator CSR Type: num ... ok
Operator CSR Type: phase ... ok
Operator CSR Type: position ... ok
Operator CSR Type: qdiags ... ok
Operator CSR Type: qeye/identity ... ok
Operator CSR Type: qutrit ops ... ok
Operator CSR Type: squeeze ... ok
Operator CSR Type: squeezing ... ok
Operator CSR Type: zero_oper ... ok
Spin 1/2 operators ... ok
Spin 3/2 operators ... ok
Spin 2 operators ... ok
Spin 5/2 operators ... ok
Destruction operator ... ok
Creation operator ... ok
Identity operator ... ok
Identity operator ... ok
Number operator ... ok
Squeezing operator ... ok
Displacement operator ... ok
parfor ... ok
parallel_map ... ok
serial_map ... ok
parfor ... ok
partial transpose of bipartite systems ... ok
partial transpose: comparing sparse and dense implementations ... ok
partial transpose: randomized tests on tripartite system ... ok
qutip.tests.test_qft.TestQFT.testQFTComparison ... ok
Qobj data ... ok
Qobj type ... ok
Qobj Hermicity ... ok
Qobj shape ... ok
Qobj addition ... ok
Qobj subtraction ... ok
Qobj multiplication ... ok
Qobj division ... ok
Qobj power ... ok
Qobj negation ... ok
Qobj equals ... ok
Qobj getitem ... ok
Qobj multiplication type ... ok
Qobj conjugate ... ok
Qobj adjoint (dagger) ... ok
Qobj diagonals ... ok
Qobj eigenenergies ... ok
Qobj eigenstates ... ok
Qobj expm ... ok
Qobj expm (explicit sparse) ... ok
Qobj expm (explicit dense) ... ok
Qobj sqrtm ... ok
Qobj full ... ok
Qobj norm ... ok
Qobj permute ... ok
Qobj ket type ... ok
Qobj bra type ... ok
Qobj operator type ... ok
Qobj superoperator type ... ok
Checks that binary ops preserve 'superrep'. ... ok
Checks that binary ops preserve 'superrep'. ... ok
Checks that binary ops preserve 'superrep'. ... ok
Checks that binary ops preserve 'superrep'. ... ok
Checks that binary ops preserve 'superrep'. ... ok
Checks that binary ops preserve 'superrep'. ... ok
Checks that binary ops preserve 'superrep'. ... ok
Checks that binary ops preserve 'superrep'. ... ok
Checks that binary ops preserve 'superrep'. ... ok
Checks that binary ops preserve 'superrep'. ... ok
Checks that binary ops preserve 'superrep'. ... ok
Checks that binary ops preserve 'superrep'. ... ok
Checks that binary ops preserve 'superrep'. ... ok
Checks that binary ops preserve 'superrep'. ... ok
Checks that binary ops preserve 'superrep'. ... ok
mul and tensor of skew-Hermitian operators report ``isherm = True``. ... ok
Tensor: Checks that super_tensor respects states. ... ok
Tensor: Super_tensor correctly tensors on underlying spaces. ... ok
Composite: Tests compositing unitaries and superoperators. ... ok
Composite: Tests compositing states and density operators. ... ok
quantum process tomography for snot gate ... ok
quantum process tomography for cnot gate ... ok
Test mesolve qubit, with dissipation ... ok
Test mesolve qubit, no dissipation ... ok
Test essolve qubit, with dissipation ... ok
Test mcsolve qubit, with dissipation ... ok
Test mcsolve qubit, no dissipation ... ok
qutip.tests.test_qubitcircuit.TestQubitCircuit.testCNOTtoCSIGN ... ok
qutip.tests.test_qubitcircuit.TestQubitCircuit.testCNOTtoISWAP ... ok
qutip.tests.test_qubitcircuit.TestQubitCircuit.testCNOTtoSQRTISWAP ... ok
qutip.tests.test_qubitcircuit.TestQubitCircuit.testCNOTtoSQRTSWAP ... ok
qutip.tests.test_qubitcircuit.TestQubitCircuit.testCSIGNtoCNOT ... ok
qutip.tests.test_qubitcircuit.TestQubitCircuit.testISWAPtoCNOT ... ok
qutip.tests.test_qubitcircuit.TestQubitCircuit.testSWAPtoCNOT ... ok
qutip.tests.test_qubitcircuit.TestQubitCircuit.testadjacentgates ... ok
qutip.tests.test_rand.TestRand.testRandUnitary ... ok
qutip.tests.test_rand.TestRand.testRanddm ... ok
qutip.tests.test_rand.TestRand.testRandherm ... ok
qutip.tests.test_rand.TestRand.testRandket ... ok
Sparse eigs Hermitian ... ok
Sparse eigs non-Hermitian ... ok
Sparse eigvals only Hermitian. ... ok
Dense eigs Hermitian. ... ok
Dense eigs non-Hermitian ... ok
Dense eigvals only Hermitian ... ok
Sparse: Symmetric Permute ... ok
Sparse: Nonsymmetric Permute ... ok
Sparse: Symmetric Reverse Permute ... ok
Sparse: Nonsymmetric Reverse Permute ... ok
Sparse: Bandwidth ... ok
Sparse: Profile ... ok
qutip.tests.test_spinchain.TestSpinChain.test_circular_ISWAP ... ok
qutip.tests.test_spinchain.TestSpinChain.test_circular_SQRTISWAP ... ok
qutip.tests.test_spinchain.TestSpinChain.test_circular_combination ... ok
qutip.tests.test_spinchain.TestSpinChain.test_linear_ISWAP ... ok
qutip.tests.test_spinchain.TestSpinChain.test_linear_SQRTISWAP ... ok
qutip.tests.test_spinchain.TestSpinChain.test_linear_combination ... ok
State CSR Type: basis ... ok
State CSR Type: bell_state ... ok
State CSR Type: bra ... ok
State CSR Type: coherent ... ok
State CSR Type: coherent_dm ... ok
State CSR Type: fock ... ok
State CSR Type: fock_dm ... ok
State CSR Type: ghz_state ... ok
State CSR Type: ket ... ok
State CSR Type: ket2dm ... ok
State CSR Type: maximall_mixed_dm ... ok
State CSR Type: spin_coherent ... ok
State CSR Type: spin_state ... ok
State CSR Type: thermal_dm ... ok
State CSR Type: w_state ... ok
qutip.tests.test_states.TestStates.testCoherentDensityMatrix ... ok
qutip.tests.test_states.TestStates.testFockDensityMatrix ... ok
qutip.tests.test_states.TestStates.testThermalDensityMatrix ... ok
Steady state: Thermal qubit - direct solver ... ok
Steady state: Thermal qubit - eigen solver ... ok
Steady state: Thermal qubit - power solver ... ok
Steady state: Thermal qubit - iterative-gmres solver ... ok
Steady state: Thermal qubit - iterative-bicgstab solver ... ok
Steady state: Thermal qubit - iterative-lgmres solver ... ok
Steady state: Thermal HO - direct solver ... ok
Steady state: Thermal HO - eigen solver ... ok
Steady state: Thermal HO - power solver ... ok
Steady state: Thermal HO - iterative-gmres solver ... ok
Steady state: Thermal HO - iterative-bicgstab solver ... ok
Steady state: Thermal HO - iterative-lgmres solver ... ok
Steady state: Driven cavity - direct solver ... ok
Steady state: Driven cavity - eigen solver ... ok
Steady state: Driven cavity - power solver ... ok
Steady state: Driven cavity - iterative-gmres solver ... ok
Steady state: Driven cavity - iterative-bicgstab solver ... ok
Steady state: Driven cavity - iterative-lgmres solver ... ok
Stochastic: smesolve: photo-current ... ok
Stochastic: smesolve: homodyne ... ok
Stochastic: smesolve: heterodyne ... ok
Stochastic: ssesolve: photo-current ... ok
Stochastic: ssesolve: homodyne ... ok
Stochastic: ssesolve: heterodyne ... ok
qutip.tests.test_subsys_apply.TestSubsysApply.test_ComplexSingleApply ... ok
qutip.tests.test_subsys_apply.TestSubsysApply.test_ComplexSuperApply ... ok
qutip.tests.test_subsys_apply.TestSubsysApply.test_SimpleSingleApply ... ok
qutip.tests.test_subsys_apply.TestSubsysApply.test_SimpleSuperApply ... ok
qutip.tests.test_subsystem_apply.TestSubsystemApply.test_ComplexSingleApply ... ok
qutip.tests.test_subsystem_apply.TestSubsystemApply.test_ComplexSuperApply ... ok
qutip.tests.test_subsystem_apply.TestSubsystemApply.test_SimpleSingleApply ... ok
qutip.tests.test_subsystem_apply.TestSubsystemApply.test_SimpleSuperApply ... ok
qutip.tests.test_superop_reps.TestSuperopReps.test_ChoiKrausChoi ... ok
qutip.tests.test_superop_reps.TestSuperopReps.test_ChoiPreservesSelf ... ok
qutip.tests.test_superop_reps.TestSuperopReps.test_SuperChoiChiSuper ... ok
qutip.tests.test_superop_reps.TestSuperopReps.test_SuperChoiSuper ... ok
qutip.tests.test_superop_reps.TestSuperopReps.test_SuperPreservesSelf ... ok
qutip.tests.test_superop_reps.TestSuperopReps.test_choi_tr ... ok
qutip.tests.test_superop_reps.TestSuperopReps.test_known_iscptp ... ok
qutip.tests.test_superop_reps.TestSuperopReps.test_random_iscptp ... ok
qutip.tests.test_superoper.TestMatVec.testLiouvillianImplem ... ok
qutip.tests.test_superoper.TestMatVec.testMatrixVecMat ... ok
qutip.tests.test_superoper.TestMatVec.testOperatorSpostAppl ... ok
qutip.tests.test_superoper.TestMatVec.testOperatorSpreAppl ... ok
qutip.tests.test_superoper.TestMatVec.testOperatorUnitaryTransform ... ok
qutip.tests.test_superoper.TestMatVec.testOperatorVector ... ok
qutip.tests.test_superoper.TestMatVec.testVecMatIndexCompability ... ok
qutip.tests.test_superoper.TestMatVec.testVecMatIndexConversion ... ok
qutip.tests.test_superoper.TestMatVec.testVecMatVec ... ok
qutip.tests.test_superoper.TestMatVec.test_reshuffle ... ok
qutip.tests.test_superoper.TestMatVec.test_sprepost ... ok
qutip.tests.test_superoperator.TestMatrixVector.testLiouvillianImplementations ... ok
qutip.tests.test_superoperator.TestMatrixVector.testMatrixVectorMatrix ... ok
qutip.tests.test_superoperator.TestMatrixVector.testOperatorSpostApplication ... ok
qutip.tests.test_superoperator.TestMatrixVector.testOperatorSpreApplication ... ok
qutip.tests.test_superoperator.TestMatrixVector.testOperatorUnitaryTransformation ... ok
qutip.tests.test_superoperator.TestMatrixVector.testOperatorVector ... ok
qutip.tests.test_superoperator.TestMatrixVector.testVectorMatrixIndexCompability ... ok
qutip.tests.test_superoperator.TestMatrixVector.testVectorMatrixIndexConversion ... ok
qutip.tests.test_superoperator.TestMatrixVector.testVectorMatrixVector ... ok
qutip.tests.test_superoperator.TestMatrixVector.test_reshuffle ... ok
qutip.tests.test_superoperator.TestMatrixVector.test_sprepost ... ok
td_format_check: monte-carlo ... ok
qutip.tests.test_tensor.test_flatten ... ok
qutip.tests.test_tensor.test_enumerate_flat ... ok
qutip.tests.test_tensor.test_deep_remove ... ok
qutip.tests.test_tensor.test_unflatten ... ok
qutip.tests.test_tensor.test_tensor_contract ... ok
Three-level atom: States ... ok
Three-level atom: Operators ... ok
utilities: energy unit conversions ... ok
wigner: test wigner function calculation for coherent states ... ok
wigner: test wigner function calculation for Fock states ... ok
wigner: compare wigner methods for random density matrices ... ok
wigner: compare wigner methods for random state vectors ... ok
Wigner: Compare Wigner fft and iterative for rand. ket ... ok
Wigner: Compare Wigner fft and iterative for rand. dm ... ok
----------------------------------------------------------------------
Ran 336 tests in 518.423s
OK
In [1]:
from qutip.ipynbtools import version_table
version_table()
Out[1]:
Software Version matplotlib 1.4.2 IPython 2.3.1 Cython 0.21.2 Numpy 1.9.1 QuTiP 3.1.0 SciPy 0.15.1 Python 3.4.2 (default, Oct 19 2014, 13:34:42)
[GCC 4.2.1 Compatible Apple LLVM 6.0 (clang-600.0.51)] OS posix [darwin] Tue Feb 03 15:45:41 2015 EST
In [ ]:
Content source: rouxinol/HybridSystem
Similar notebooks: