Weird ground state for my hamiltonian #2433
Replies: 2 comments
-
By the way, the plot of the wells are scaled to be relevant on the graph |
Beta Was this translation helpful? Give feedback.
0 replies
-
Okay guys I found my error ! To define my gaussian well, I used .expm(), it is better if you do |
Beta Was this translation helpful? Give feedback.
0 replies
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
-
Hello everyone,
This is my first post on a forum in general, so I'll do my best to share and answer your questions ! (Sorry for the English, it is not my mother tongue). I first give the idea.
Here is my problem :
At the beginning of a run, I ask Qutip to give me the ground state at T=0, and to plot the wave function and the probability in space. At T=0, my well is not moving.
According to some approximations of a gaussian well by an harmonic one, I should have a density probability with a standard deviation given by a formula. But Qutip gives me a bigger one (and I can see it on the plot). The well is not moving at this moment. The result of Qutip is at the center of my well, everything is perfect, just is it too big according to theory.
I have not found any mistake for the moment in my code...
I gave the total Hamiltonian divided by h_barre in sesolve function.
If I ask to plot other eigenstate, I can see harmonic-oscillator-eigenstates-like curves, they are just a bit to much spread according to theory.
I gave some pictures of the plot, one with the spread gaussian well ground state, and the other if I put an harmonic well.
Feel free to ask more precision !
I want to know why Qutip cannot give me a more correct ground state. Furthermore, if I change the waist of my gaussian well, the ground state is still the same... (probably my gaussian well is defined with a constant term in the exponential, but it is not the case)
Beta Was this translation helpful? Give feedback.
All reactions