In[]:=
n=4;Ψ=ConstantArray[0,n];Ψ[[1]]= 0.03596901902499999|0,0〉|30,30〉+0.056822724205|0,1〉|30,29〉;
In[]:=
SetFockSpaceSize[20];aL=TensorProductOperator[AnnihilationOperator[],IdentityOperator[]];aR=TensorProductOperator[IdentityOperator[],AnnihilationOperator[]];ApplyOperator[operator_,state_]:=operator.state;ApplyOperator[aL,Ψ[[1]]]
Out[]=
TensorProductOperator[AnnihilationOperator[],IdentityOperator[]].(0.035969|0,0〉|30,30〉+0.0568227|0,1〉|30,29〉)