1. Limited time only! Sign up for a free 30min personal tutor trial with Chegg Tutors
    Dismiss Notice
Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Homework Help: Quantum particle in a magnetic field

  1. May 22, 2014 #1
    1. The problem statement, all variables and given/known data
    A particle with electrical charge [itex]q[/itex] and mass [itex]m[/itex] is in a electromagnetic field described by [itex]\phi (\vec{r}, t)[/itex] and [itex]A(\vec{r}, t)[/itex]. Its Hamiltonian is as follows:
    [itex]H = \frac{1}{2m} \left ( \frac{\hbar}{i}\vec{\nabla}-\frac{q}{c} \vec{A} (\vec{r}, t) \right ) ^2 +q\phi (\vec{r}, t)[/itex]

    The conservation of charge guarantees the continuity equation is fulfilled:
    [itex]\frac{\partial}{\partial t} \rho (\vec{r},t)+\vec{\nabla}·\vec{j}(\vec{r},t) = 0[/itex],
    where [itex]\rho = q\left|\Psi(x)\right|^2[/itex] is the charge density.

    Find the current density [itex]\vec{j}(\vec{r},t)[/itex]

    2. Relevant equations

    The Hamiltonian of a particle in an electromagnetic field.
    [itex]H = \frac{1}{2m} \left ( \frac{\hbar}{i}\vec{\nabla}-\frac{q}{c} \vec{A} (\vec{r}, t) \right ) ^2 +q\phi (\vec{r}, t)[/itex]

    The continuity equation
    [itex]\frac{\partial}{\partial t} \rho (\vec{r},t)+\vec{\nabla}·\vec{j}(\vec{r},t) = 0[/itex]

    3. The attempt at a solution
    I've tried stating the time-dependent Shcrodinger equation [itex]\hat{H}\Psi (x) = i\hbar\frac{\partial \Psi(x)}{\partial t}[/itex] and solve for the time-derivative of the wavefunction, which gives:
    [itex]\frac{\partial \Psi}{\partial t} = \left [ \frac{\vec{\nabla}^2}{2m}+\frac{q}{\hbar imc}\vec{A}\vec{\nabla}-\frac{1}{2\hbar^2 m} \left( \frac{q}{c} \vec{A}\right )^2 -\frac{q}{2 \hbar m}\phi\right]\Psi[/itex]


    [itex]\frac{\partial\rho}{\partial t} = q2|\Psi(x)|\frac{\partial |\Psi(x)|}{\partial t} = \left [ \frac{q\vec{\nabla}^2}{m}+\frac{2q^2}{\hbar imc}\vec{A}\vec{\nabla}-\frac{q}{\hbar^2 m} \left( \frac{q}{c} \vec{A}\right )^2 -\frac{q^2}{ \hbar m}\phi\right]\Psi^2[/itex]

    And then, from the continuity equation we get
    [itex]\vec{\nabla}·\vec{j}(\vec{r},t) = \frac{\partial}{\partial t} \rho (\vec{r},t) = \left [ \frac{q\vec{\nabla}^2}{m}+\frac{2q^2}{\hbar imc}\vec{A}\vec{\nabla}-\frac{q}{\hbar^2 m} \left( \frac{q}{c} \vec{A}\right )^2 -\frac{q^2}{ \hbar m}\phi\right]\Psi^2[/itex]

    But I'm stuck there. How do I "remove" the nabla operators?
    Last edited: May 22, 2014
  2. jcsd
  3. May 22, 2014 #2
    Question edited with more information.
  4. May 22, 2014 #3


    User Avatar
    Gold Member

    You must be careful at your steps and how you deal with complex functions and operators.
    Problematic points in your approach
    No1: your derivative of [itex]\rho[/itex] is not a useful thing... it's better to try to write [itex]|\psi|^{2}= \psi^{\dagger} \psi[/itex] and then use the derivative on each... [itex](\frac{\partial \psi}{\partial t}) \psi^{\dagger} + \psi \frac{\partial \psi^{\dagger}}{\partial t}[/itex]
    No2: you had [itex]|\psi|[/itex] on the left and then you moved it on the right, passing it through operators like nabla to write [itex]\psi^{2}[/itex]
    No3: You don't know the derivative of [itex]|\psi|[/itex]... as you've written it, it must be a real number, yet you have real or complex operators [itex] \frac{q^{2}}{hm} \phi[/itex] or [itex] \frac{2q^{2}}{himc} A∇[/itex] acting on a complex wavefunction...
  5. May 23, 2014 #4
    Yeah, that didn't feel right at all. Should I have left the wavefunction at the left of the operators after expanding the squared term?

    I don't get what's the problem with my derivative. Can't I use the chain rule?
  6. May 23, 2014 #5


    User Avatar
    Gold Member

    your time derivative gives you quantities that are not easy to work on... Schroedinger's eq. gives you the time evolution of the wavefunction, not its module. Both yours and mine use the chain rule, but yours is not useful- you cannot work easily on it, if you can at all... Doing the otherway, you can always take the conjugate of the Schr eq and insert the expression for [itex] \frac{\partial \psi^{*}}{\partial t}[/itex]
    So even if you leave your wf at the left, things are not right because you put wrong expression for the time derivative.
  7. May 23, 2014 #6
    Oh, right, what I did looks stupid now. Thank you!

    Now: how would the nabla operator work on the wavefunction? Since [itex]\Psi(x)[/itex] is an scalar I'm guessing:
    [itex]\vec{\nabla}\Psi(x) = \left (\frac{\partial \Psi(x)}{\partial x}, \frac{\partial \Psi(x)}{\partial y}, \frac{\partial \Psi(x)}{\partial z} \right ) = \frac{\partial \Psi(x)}{\partial x} · \vec{\hat{x}}[/itex], is that right?

    And then [itex]\vec{\nabla}^2\Psi(x) = \frac{\partial ^2\Psi(x)}{\partial x^2}[/itex]
    Last edited: May 23, 2014
  8. May 23, 2014 #7


    User Avatar
    Gold Member

    No... The [itex]x[/itex] in the wavefunction doesn't really mean the [itex]x[/itex] variable alone- that would be true only if you had 1D problem, and then you nabla would just be the partial of x alone. The [itex]x[/itex] in the wavefunction is a vector, that means:
    [itex] \psi(x)= \psi(\vec{x})=\psi(x,y,z)[/itex]

    why do you want to act with the nabla on the wavefunction in the first place?
  9. May 23, 2014 #8
    I don't really need that for this problem, I just wanted to know how would that work. Maybe I shouldn't have put it here, sorry.

    Using what you said I've managed to find a very simple solution. I'd appreciate if anyone gave a look at it:
    [itex]i\hbar \frac{\partial \psi}{\partial t}= \hat{H} \psi \rightarrow \frac{\partial \psi}{\partial t} = \frac{1}{i\hbar} \hat{H} \psi =\frac{1}{2mi\hbar} \left ( \frac{\hbar}{i}\vec{\nabla}-\frac{q}{c} \vec{A} (\vec{r}, t) \right ) ^2 \psi+ \frac{q\phi}{i\hbar} \psi[/itex]
    Expanding the square:
    [itex]\frac{\partial \psi}{\partial t} = \frac{1}{2i\hbar m} \left( -\hbar ^2\vec{\nabla}^2 - \frac{\hbar q}{ic}\vec{\nabla}\vec{A} +\left ( \frac{q}{c} \vec{A}\right )^2\right )\psi + \frac{q\phi}{i\hbar} \psi[/itex]
    [itex]\frac{\partial \psi}{\partial t} = \left ( \frac{i\hbar}{2m} \vec{\nabla}^2 + \frac{q}{2mc} \vec{\nabla}\vec{A} - \frac{i}{2m\hbar} \left( \frac{q}{c} \vec{A}\right)^2 - \frac{iq\phi}{\hbar} \right )\psi[/itex]
    And therefore:
    [itex]\frac{\partial \psi ^\dagger}{\partial t} = \left ( -\frac{i\hbar}{2m} \vec{\nabla}^2 + \frac{q}{2mc} \vec{\nabla}\vec{A} + \frac{i}{2m\hbar} \left( \frac{q}{c} \vec{A}\right)^2 + \frac{iq\phi}{\hbar} \right )\psi[/itex]
    [itex]\frac{\partial \psi}{\partial t}\psi ^\dagger = \left ( \frac{i\hbar}{2m} \vec{\nabla}^2 + \frac{q}{2mc} \vec{\nabla}\vec{A} - \frac{i}{2m\hbar} \left( \frac{q}{c} \vec{A}\right)^2 - \frac{iq\phi}{\hbar} \right )|\psi|^2[/itex] and similarly for [itex]\psi\frac{\partial \psi^\dagger}{\partial t}[/itex]
    And from there is easy:
    [itex]\frac{\partial \rho}{\partial t} = q\frac{\partial }{\partial t}|\psi|^2 = q \frac{\partial}{\partial t} \psi ^\dagger \psi = q \left(\frac{\partial \psi^\dagger}{\partial t} \psi + \psi ^\dagger \frac{\partial \psi}{\partial t} \right) = \left (\frac{q^2}{mc}\vec{\nabla}\vec{A}\right) |\psi|^2[/itex] (because most terms in the parentheses cancel each other).
    And then from the continuity equation I finally obtain for [itex]\vec{\nabla}\vec{j}[/itex]:
    [itex]\vec{\nabla}\vec{j} = -\frac{\partial\rho}{\partial t} = \frac{q^2}{mc}\vec{\nabla}\vec{A} |\psi|^2[/itex]

    Now my question is: can I just take those nablas out?
  10. May 23, 2014 #9


    User Avatar
    Gold Member

    You should be more careful when you expand the square of operators....
    [itex]∇A \ne A∇[/itex]
    For this kind of terms you need to do something like:
    [itex]∇(A f)= A∇f + (∇A)f [/itex]
    and then you can use the first term [itex]A∇[/itex] with the same term coming from the square...eg:
    [itex](c ∇+ d A)(c ∇+dA)= c^{2} ∇^{2} + d^{2}A^{2} + cd A∇ + cd ∇A[/itex]
    in the 3rd term, nabla acts on the wavefunction alone, in the 4th term nabla acts on both A and the wavefunction:
    [itex](c ∇+ d A)^{2}= c^{2} ∇^{2} + d^{2}A^{2} + cd A∇+ cd A∇+ cd (∇A)[/itex]
    [itex](c ∇+ d A)^{2}= c^{2} ∇^{2} + d^{2}A^{2} + 2cd A∇+ cd (∇A)[/itex]

    Also the derivative of the conjugate, you should also take the complex for [itex]\psi[/itex] at the RHS.
    Also avoid writing things like [itex](\hat{O} \psi) \psi^{*}=\hat{O} |\psi|^{2}[/itex]
    The reason is that in the first case the operator acts on [itex]\psi[/itex] and you suddenly make it act on [itex]|\psi|^{2}[/itex].
    What you should do is try to move nablas around to create a needed expression (as the one I give in the end of this post)

    In order to move nablas out, you have to do what I did above for the square (work with [itex]∇(AB)=(∇A)B + A(∇B) [/itex])
    and don't be in haste to reach a result, think carefully for each point.

    Afterwards, if you have an equation:
    [itex]∇J= ∇(Something)[/itex]
    you can deduce that
    [itex] J= Something[/itex]
    Last edited: May 23, 2014
  11. May 23, 2014 #10


    User Avatar
    Gold Member

  12. May 25, 2014 #11
    Would this be the correct Hamiltonian with the expanded square?
    [itex]\hat{H}=\frac{1}{2m}\left(-\hbar^2\vec{\nabla}^2 -\frac{\hbar q}{ic}\left ( \vec{\nabla}·\vec{A} + \vec{A}·\vec{\nabla}\right ) + \left(\frac{q}{c}\vec{A}\right)^2 \right) + q\phi[/itex]
    Last edited: May 25, 2014
  13. Jun 13, 2014 #12
    After solving a couple of formalism problems I had I´ve come to this:
    [itex] \hat{H}\psi = -\frac{\hbar^2}{2m}\nabla^2\psi-\frac{\hbar q}{2imc}(\vec{\nabla}\cdot\vec{A})\psi+\frac{\hbar q}{imc}\vec{A}\cdot(\vec{\nabla}\psi)+\frac{q^2\vec{A^2}}{2mc^2}\psi+ \frac{q \phi}{2m}\psi[/itex]

    Can you check if its right?
  14. Jun 27, 2014 #13
    I understand that in this step the [itex] cd ∇A[/itex] term at the end has been developed into [itex]cd A∇+ cd (∇A)[/itex].

    But isn´t [itex]c^{2} ∇^{2} = c^{2}∇∇[/itex] a particular case of [itex] cd ∇A[/itex] with [itex]d = c[/itex] and [itex] A = ∇[/itex]?

    Why doesn´t the same thing happen with [itex]c^{2} ∇^{2}[/itex]? That is, why don´t we have [itex]c^{2} ∇^{2} =c^{2} ∇∇+ c^{2} (∇∇)[/itex]?
  15. Jun 27, 2014 #14
    I've done it again from scratch. Can you tell me if I got it right this time before I go on?
    [itex]-\frac{\hbar}{i}\frac{\partial \Psi}{\partial t} = \hat{H}\Psi[/itex]
    [itex]-\frac{\hbar}{i}\frac{\partial \Psi}{\partial t} = \left [\frac{1}{2m} \left ( \frac{\hbar}{i}\vec{\nabla}-\frac{q}{c} \vec{A} (\vec{r}, t) \right ) ^2 +q\phi (\vec{r}, t)\right ]\Psi[/itex]
    [itex]-\frac{\hbar}{i}\frac{\partial \Psi}{\partial t} = \left [ \frac{1}{2m}\left ( -\hbar ^2 \nabla ^2-\frac{\hbar q}{ic}\left ( (\nabla A) + (A \nabla)\right ) +\frac{q^2}{c^2}A^2\right ) +q\phi\right ]\Psi[/itex]
  16. Jun 27, 2014 #15


    User Avatar
    Gold Member

    you must understand that it's not a matter of "definition", but a simple fact that you are working with operators. You have an operator [itex]\hat{O}[/itex]
    which in general doesn't commute with what is right to it. In this case you have partial derivatives of a function which depends on position [itex]\vec{A}[/itex].
    And you have that operator acting on the wavefunction after that...so in lines:
    [itex] \hat{O}_{x} A(x) \psi(x) [/itex]
    If the operator [itex]\hat{O}_{x} \equiv \frac{\partial}{\partial x}[/itex]
    what will you have?

    [itex] \frac{\partial}{\partial x} A(x) \psi(x) [/itex]

    [itex] \frac{\partial A(x)}{\partial x} \psi(x) + A(x) \frac{\partial \psi(x) }{\partial x} = (\frac{\partial A(x)}{\partial x}+ A(x) \frac{\partial}{\partial x}) \psi(x) [/itex]

    If you write [itex]\hat{G} \psi[/itex] , then the operator [itex]\hat{G}=\frac{\partial A(x)}{\partial x}+ A(x) \frac{\partial}{\partial x}[/itex]
    And this relation holds in general, even if [itex]O,A[/itex] commutes...in that case the first term in G would be zero. That term corresponds to the commutator of [itex][O,A][/itex] of O as I wrote above. [itex]O A \psi = A O \psi + [O,A] \psi = (A O + [O,A]) \psi [/itex]
    In the case you set [itex]A= ∇[/itex] you have [itex]∇∇[/itex] and the two operators commute- and no extra term has to appear.

    (eg another more complicated would be to try to find the [itex] [ ∇^{2}, A(x) ] [/itex]

    [itex] ∇^{2} Ag = ∇∇(Ag) = ∇ [ (∇A) g + A (∇g) ] [/itex]
    [itex]= (∇^{2}A) g + (∇A) (∇g) + (∇A)(∇g) + A (∇^{2}g)[/itex]

    so as an operator, taking g out:

    [itex] ∇^{2}A = (∇^{2}A) + 2 (∇A) ∇ + A ∇^{2} [/itex]
    or that

    [itex] [ ∇^{2}, A(x) ] = (∇^{2}A) + 2 (∇A) ∇ [/itex]
    Last edited: Jun 27, 2014
  17. Jun 27, 2014 #16


    User Avatar
    Gold Member

    I know that this might seem confusing, or tedious work, but once you do some scholastic calculations with it,, being aware and careful at each step, it's easy to "grab the sense" and it'll be faster for you next time. But you have to do them, because it's a matter of understanding operators and how they act or how you can treat them...as you may already know the result of multiplying two operators is not as it may seem at first... and the confusion is because of the commutation relations.
    I guess that since you are dealing with the Electromagnetism Hamiltonian you have already seen the Hamiltonian of the Harmonic Oscillator in the form of the ladder operators and how you find not [itex]H= \alpha^{\dagger} \alpha [/itex] but you also have an additional [itex]+\frac{1}{2}[/itex] term... Although [itex]\alpha= a x + i c p [/itex] and [itex] \alpha^{\dagger}= ax - i c p [/itex] with [itex]a,c[/itex] chosen appropriately...
    The same calculations apply in showing the Hamiltonian form from them...

    [itex]\alpha^{\dagger} \alpha = (ax - i c p )(a x + i c p) \ne a^{2} x^{2} + c^{2} p^{2}=H [/itex]
    [itex] \alpha^{\dagger} \alpha = a^{2} x^{2} + c^{2} p^{2} + iac xp - iac px = H + iac [x,p][/itex]

    NEVERMIND, for your question:

    [itex](a ∇+b A) ( a ∇ + b A) \psi[/itex]

    To make more obvious what the 1st parenthesis operator acts on:

    [itex](a ∇+b A) ( a ∇\psi + b A \psi) [/itex]

    multiplying each one with the other:

    [itex] a^{2} ∇^{2} \psi + b^2 A^2 \psi + ab ∇(A \psi) + ba A (∇\psi) [/itex]

    Using Leibniz Rule for the derivative in the 3rd term:

    [itex] a^{2} ∇^{2} \psi + b^2 A^2 \psi + ab (∇A) \psi + ab A (∇ \psi) + ba A (∇ \psi) [/itex]

    adding the same terms:

    [itex] a^{2} ∇^{2} \psi + b^2 A^2 \psi + ab (∇A) \psi + 2 ab A (∇ \psi) [/itex]

    taking out [itex]\psi[/itex] since it's all in the right:

    [itex] (a^{2} ∇^{2} + b^2 A^2 + ab (∇A) + 2 ab A ∇ ) \psi [/itex]

    Identifying the operator I had to act on [itex]\psi[/itex] at first with the one I extracted in the last step:

    [itex](a ∇+b A) ( a ∇ + b A)= a^{2} ∇^{2} + b^2 A^2 + ab (∇A) + 2 ab A ∇[/itex]

    In your case [itex] a= \frac{h}{i} [/itex] and [itex]b= -\frac{q}{c}[/itex]

    [itex](\frac{h}{i} ∇- \frac{q}{c} A)^{2} = (\frac{h}{i})^{2} ∇^{2} + (\frac{q}{c})^2 A^2 - \frac{hq}{ic} (∇A) - \frac{2hq}{ci} A ∇[/itex]

    [itex](\frac{h}{i} ∇- \frac{q}{c} A)^{2} = -h^{2} ∇^{2} + \frac{q^{2}}{c^{2}} A^2 - \frac{hq}{ic} [ (∇A) + 2 A ∇ ] [/itex]
    Last edited: Jun 27, 2014
  18. Jun 28, 2014 #17
    I think my problem is that I lack a proper way of denoting what is an operator and what is a vector or a derivative. Let me ask a general question:

    If we were to denote "chained" operators like [itex]\hat{A_1}(\hat{A_2}(...(\hat{A_n}\Psi)...) ...)) [/itex] when is it appropriate to take out the hats anI would say the correct way of solving that would be to apply the transformations of the first operator, then the second, then the third... until we have just the wavefunction or derivatives of it along with some other terms. Is this right?
    d apply the transformations corresponding to the operators (derivative, multiplication...)?

    I would say the correct way of solving that would be to apply the transformations of the first operator, then the second, then the third... until we have just the wavefunction or derivatives of it along with some other terms. Is this right?

    Also, for the case we are on: if we write ∇∇ in component form [itex]\partial q_i(\partial q_j \Psi)[/itex] (with summation over i and j) we will have [itex]\partial q_i(\partial q_j \Psi) = (\partial q_i\partial q_j)\Psi + \partial q_j \partial q_i(\Psi) = 2\partial q_j \partial q_i(\Psi) = 2∇²[/itex]. I guess I did something wrong somewhere, but I don't know where.

    I can actually move forward without understanding this, but I'd like to get my head around it.

    PD: Thank you very much for taking the time to write those answers.
  19. Jun 28, 2014 #18


    User Avatar
    Gold Member

    The operators with the wavefunctions as that act on what is right to them.
    When you have something of the form :
    [itex]A_{1} A_{2} f(x) [/itex]
    Then it means that [itex]A_{1}[/itex] acts on the result you have from [itex]A_{2}f(x)[/itex]

    [itex]A_{1} (A_{2} f(x)) [/itex]

    To generalize it if you have:
    [itex] A_{n} A_{n-1} .... A_{2} A_{1} \psi(x) = A_{n} (A_{n-1} (.... (A_{2} (A_{1} \psi(x) )))...)) [/itex]
    So [itex]A_{2}[/itex] will have to act on your result of [itex]A_{1} \psi(x)[/itex], then the result is acted on by [itex]A_{3}[/itex] then the result is acted on blah blah blah... until you'll get the result of [itex]A_{n-1} G(x)[/itex] (G(x) the general result you achieved from the previous actions) on which you'll act with [itex]A_{n}[/itex].

    In case [itex]A_{1}=A_{2}[/itex] (or a more general statement is that if [itex]A_{1},A_{2}[/itex] commute), then:

    [itex] A_{1} A_{2} f(x) = A_{2} A_{1} f(x) [/itex]

    In your case the main thing is that you don't have a commutor:

    [itex] \partial_{i} \partial_{j} \psi [/itex]
    you can interchange the two partial derivatives without any extra term appearing because [itex][ \partial_{i}, \partial_{j}]=0[/itex].
    [itex] \partial_{j} \partial_{i} \psi [/itex]

    As I said, the extra terms is not just a definition so that you can apply them any where- they are a result of the operators not commuting. In the EM case you have [itex]\partial[/itex] acting on a function of [itex]x,~~A(x)[/itex]... So they don't really have to commute (as the mommentum operator p -partial derivative wrt to x, in position representation- does not commute with position x- function of x, A(x)=x , in position repr).

    Now going back to my [itex]A_{1,2}[/itex] operators. If they didn't commute:
    [itex] A_{1} A_{2} \psi = A_{2} A_{1} \psi + [A_{1},A_{2}] \psi [/itex]
    you can see that immediately by replacing the commutator.
    [itex]A_{2} A_{1} \psi + [A_{1},A_{2}] \psi = A_{2} A_{1} \psi + A_{1}A_{2} \psi - A_{2}A_{1} \psi = A_{1}A_{2} \psi[/itex]

    This is general, when [itex]A_{1},A_{2}[/itex] commute, then the extra term [itex][A_{1},A_{2}]=0 [/itex] and so it doesn't appear... that's why it doesn't with the partial derivatives you gave before...Because:
    [itex]\frac{\partial^{2} f(x)}{\partial x \partial x} =\frac{\partial^{2} f(x)}{\partial x \partial x}[/itex]
    or to make ti more obvious:
    [itex]\frac{\partial^{2} f(x,y)}{\partial x \partial y} =\frac{\partial^{2} f(x,y)}{\partial y \partial x}[/itex]
    Last edited: Jun 28, 2014
  20. Jun 28, 2014 #19


    User Avatar
    Gold Member

    Suppose you have a vector:

    On which you act on with the two matrices.....
    [itex] A_{1}=
    1 & 0\\
    0 & -1 \\

    [itex] A_{2}=

    0 & 1\\
    1 & 0 \\

    What's the result of:
    [itex]A_{1} A_{2}[/itex]

    What's the result of:
    [itex] A_{2} A_{1}[/itex]

    Do [itex]A_{1,2}[/itex] commute?
    The fact that they don't commute means that their action on [itex]x[/itex] depends on the order they act on... [itex]A_{1}A_{2} x \ne A_{2} A_{1} x [/itex]
    The extra term that appears is their commutator...
    [itex]A_{1}A_{2} x = A_{2} A_{1} x + [A_{1},A_{2}] x[/itex]
    you can put in the matrices and reconfirm the above result.

    Suppose now, instead of [itex] A_{1} [/itex] as given above, you write 1 in the place of -1, then what's
    [itex]A_{1} A_{2}[/itex]
    [itex] A_{2} A_{1}[/itex]
    Do they commute?
    You can then use:
    [itex]A_{1} A_{2} x = A_{2} A_{1} x + [A_{1},A_{2}]x [/itex]
    and show
    [itex] A_{1} A_{2}x= A_{2} A_{1} x [/itex]
  21. Jun 28, 2014 #20


    User Avatar
    Staff Emeritus
    Science Advisor
    Gold Member

  22. Jun 28, 2014 #21
    I know all that, I just can't make it fit with the ∇∇ case when I write it explicitly [itex]\partial q_i(\partial q_j \Psi) = (\partial q_i\partial q_j)\Psi + \partial q_j \partial q_i(\Psi) = 2\partial q_j \partial q_i(\Psi) = 2∇²[/itex].

    But since it must be that I'm just missing some stupid detail I'm gonna leave it for now and I will tackle this another time.

    Thank you very much, ChrisVer, for being so patient with me. I really appreciate it.

    Going back to the problem. If I state the Schrödinger equation with the expanded Hamiltonian I get:

    [itex]i\hbar\frac{\partial \Psi}{\partial t} = \frac{1}{2m}( -h^{2} ∇^{2} + \frac{q^{2}}{c^{2}} A^2 - \frac{hq}{ic} [ (∇A) + 2 A ∇ ] )\Psi + q \phi \Psi[/itex]
    [itex]\frac{\partial \Psi}{\partial t} = \frac{1}{2i\hbar m}( -h^{2} ∇^{2} + \frac{q^{2}}{c^{2}} A^2 - \frac{hq}{ic} [ (∇A) + 2 A ∇ ] )\Psi + \frac{q \phi}{i\hbar} \Psi[/itex]

    And for the conjugate
    [itex]\frac{\partial \Psi ^*}{\partial t} = \frac{-1}{2i\hbar m}( -h^{2} ∇^{2} + \frac{q^{2}}{c^{2}} A^2 + \frac{hq}{ic} [ (∇A) + 2 A ∇ ] )\Psi ^*- \frac{q \phi}{i\hbar} \Psi^*[/itex]
    [itex]\frac{\partial \Psi ^*}{\partial t} = \frac{1}{2i\hbar m}( h^{2} ∇^{2} - \frac{q^{2}}{c^{2}} A^2 - \frac{hq}{ic} [ (∇A) + 2 A ∇ ] )\Psi ^*- \frac{q \phi}{i\hbar} \Psi^*[/itex]

    So now the charge density derivative looks like
    [itex]\frac{\partial \rho}{\partial t} = q \frac{\partial}{\partial t}\vert \Psi ^2 \vert =q \frac{\partial}{\partial t}\Psi ^* \Psi = q\left ( \frac{\partial\Psi ^*}{\partial t} \Psi + \frac{\partial\Psi}{\partial t}\Psi^* \right)[/itex]

    Now I have some doubts. If I put in the above expressions I will get terms like ## \Psi\frac{-q^2}{c^2}A^2\Psi^*\Psi## and## \Psi\frac{q^2}{c^2}A^2\Psi\Psi^*## which cancel, but also terms like ##-\hbar^2 \nabla ^2\Psi ^* \Psi ## and ##\hbar^2 \nabla ^2\Psi \Psi^* ##, which I'm not sure how to work with. Can I conmute the wavefunctions on the second one so that both terms look the same and I can cancel them, or does the ∇² operator apply before and therefore I have ##-\Psi \nabla ^2\Psi^* ## and ## \Psi^* \nabla ^2\Psi ##?
  23. Jun 29, 2014 #22


    User Avatar
    Gold Member

    you get terms like:
    [itex] \psi^* [\operators] \psi [/itex]
    In general the current is something which has:
    [itex] \psi^* \hat{p}' \psi - \psi \hat{p}' \psi^*[/itex]
    to this form you are supposed to reach. This form has also some physical meaning. This thing happens to give you [itex] \partial_{t} \rho = ∇ J [/itex]
    where [itex]\hat{p}'[/itex] is the canonical momentum... in case of your equation, it's not only the derivative, but also has to be modified in order to contain the [itex]A[/itex]. By checking the Schroedinger's equation as you gave, you can immediately see that [itex]\hat{p}'= \frac{h}{i} \partial + \frac{q}{c} A [/itex]

    again you must see what the operator you are using on acts on... [itex]\partial_{t} \psi [/itex] is a function, not operator.

    Have you ever worked in finding the current of the normal shroedinger equation-with no potential.
    Last edited: Jun 29, 2014
  24. Jun 29, 2014 #23


    User Avatar
    Gold Member

    Why do you write the 1st step in the equation?

    the [itex] \partial q_j \psi [/itex] is a function which you then take the derivative of [itex]\partial q_i [/itex]
    So what's the meaning of the 1st step?
    Instead of writing it as such why don't you write it like:

    [itex] \frac{\partial}{\partial q_i} \frac{\partial}{\partial q_j} \psi[/itex]

    which is by calculus:

    [itex] \frac{\partial^{2}}{\partial q_i \partial q_j} \psi[/itex]
    Last edited: Jun 29, 2014
  25. Jun 29, 2014 #24
    Because it acts on both the derivative operator and the wavefunction, so I use the pdoruct rule.
    No, I actually haven't. I'll look up some examples of that.

    Looking at that form for the current: should I have put the wavefunctions like this?

    [itex]\frac{\partial \rho}{\partial t} = q\left ( \frac{\partial\Psi ^*}{\partial t} \Psi + \Psi^*\frac{\partial\Psi}{\partial t} \right)[/itex]

    [itex]\frac{\partial \Psi ^*}{\partial t} = \Psi ^*\frac{1}{2i\hbar m}( h^{2} ∇^{2} - \frac{q^{2}}{c^{2}} A^2 - \frac{hq}{ic} [ (∇A) + 2 A ∇ ] )- \Psi^*\frac{q \phi}{i\hbar}[/itex]

    (Rather than how I wrote them here, when I took the conjugate of the wavefunction derivative and derived the product in the charge density)
    From this I think I can reach an expression like yours.
  26. Jun 29, 2014 #25


    User Avatar
    Gold Member

    have you ever used in calculus that:

    [itex] \frac{\partial}{\partial x} \frac{\partial}{\partial y} F(x,y) = (\frac{\partial}{\partial x}\frac{\partial}{\partial y}) F(x,y) + \frac{\partial}{\partial y} [\frac{\partial}{\partial x} F(x,y)] [/itex]


    If you haven't tried that for the easiest thing, I'd hardly recommend that you try it before going in more general stuff...

    From eq (3) to (4) you have to take out one ∇, by the same product rule...
    [itex] f^* (x) ∇^{2} f(x) = ∇ ( f^* (x) ∇f(x) ) - (∇f^* (x) ) (∇f(x))[/itex]
    the 2nd term in the above, cancels out by a similar appearing from the other term, and the 1st one together with one more coming from the other term again will give you something like:
    [itex]∇ (f^* (x) ∇f(x) - f(x)∇f^* (x) )[/itex]
    And you get the probability current once you identify the above with:
    [itex]∇J [/itex]
    Last edited: Jun 29, 2014
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook

Have something to add?
Draft saved Draft deleted