search

qojulia / QuantumOptics.jl

Library for the numerical simulation of closed as well as open quantum systems.
http://qojulia.org
Other
518 stars 101 forks source link

Can't see occupations for composite many-body basis #383

Open Umut-Can-Physics opened 4 months ago

Umut-Can-Physics commented 4 months ago

Example code is following:

b = NLevelBasis(3)
s = bosonstates(b, [0,1,2])
mb = ManyBodyBasis(b,s)
A = destroy(mb,1) ⊗ identityoperator(mb) # composite operator
A.basis_l # print: [ManyBody(onebodybasis=NLevel(N=3), states:10) ⊗ ManyBody(onebodybasis=NLevel(N=3), states:10)]
mb.occupations # one species basis # I can reach it
MB = tensor(mb,mb)
# But there is no occupations tag for MB object
MB.bases # (ManyBody(onebodybasis=NLevel(N=3), states:10), ManyBody(onebodybasis=NLevel(N=3), states:10))
MB.shape # 2-element Vector{Int64}: 10 10

So, ı want to see our occupations or Fock basis for composite mixture. I'm wondering how QoJulia construct basis for my A operator?