Calculating Force on a Particle using the Time Dependent Schrödinger Equation

AwesomeTrains
Messages
115
Reaction score
3
Hey!
1. Homework Statement

I've been given the time dependent schrödinger equation in momentum space and have to calculate the force, as a function of the position, acting on a particle with mass m.
\vec{F}(\vec{r})=-\nabla{V(\vec{r})}

(\frac{\vec{p}^2}{2m}-a\nabla_\vec{p}^2)\psi(\vec{p},t)=i\frac{\partial}{\partial t}\psi(\vec{p},t) and a∈ℝ

2. Homework Equations

Fourier transform:
\psi(\vec{p},t)=\frac{1}{\sqrt{2\pi}}\int \psi(\vec{r},t)e^{i\vec{p}\vec{r}} \, dr^3
\psi(\vec{r},t)=\frac{1}{\sqrt{2\pi}}\int \psi(\vec{p},t)e^{-i\vec{p}\vec{r}} \, dp^3

Dirac delta function:
\delta(\vec{p}-\vec{p}´)\frac{1}{2\pi}\int e^{i\vec{r}(\vec{p}-\vec{p}´)} \, dr

ħ=1
3. The Attempt at a Solution

I tried plugging in the expression for \psi(\vec{p},t) into the given equation, to get it transformed into position space.
I got this far:
\int (\frac{\vec{p}^2}{2m}-i\frac{\partial}{\partial t}) \int \psi(\vec{r}^{'},t) e^{i\vec{p}(\vec{r}^{'}-\vec{r})} \, dr´^3 \, dp^3=-a\psi(\vec{r},t)

I wanted to use the property of the delta function again but there's that \vec{p}^2 which is messing things up. How do I continue from here? Or is this a wrong approach/ has I made a mistake?
Any hints are appreciated :smile:
Kind regards
Alex
 
Physics news on Phys.org
AwesomeTrains said:
Hey!
1. Homework Statement

I've been given the time dependent schrödinger equation in momentum space and have to calculate the force, as a function of the position, acting on a particle with mass m.
\vec{F}(\vec{r})=-\nabla{V(\vec{r})}

(\frac{\vec{p}^2}{2m}-a\nabla_\vec{p}^2)\psi(\vec{p},t)=i\frac{\partial}{\partial t}\psi(\vec{p},t) and a∈ℝ

2. Homework Equations

Fourier transform:
\psi(\vec{p},t)=\frac{1}{\sqrt{2\pi}}\int \psi(\vec{r},t)e^{i\vec{p}\vec{r}} \, dr^3
\psi(\vec{r},t)=\frac{1}{\sqrt{2\pi}}\int \psi(\vec{p},t)e^{-i\vec{p}\vec{r}} \, dp^3

Dirac delta function:
\delta(\vec{p}-\vec{p}´)\frac{1}{2\pi}\int e^{i\vec{r}(\vec{p}-\vec{p}´)} \, dr

Hi, I am not sure how you are expected to use the Dirac Delta function, but my best guess at finding an answer would be to solve for the wave function in momentum space, then do the furier transform, and then solve for the potential in the Scroedinger Equation for position space.
 
Brage said:
Hi, I am not sure how you are expected to use the Dirac Delta function, but my best guess at finding an answer would be to solve for the wave function in momentum space, then do the furier transform, and then solve for the potential in the Scroedinger Equation for position space.

Thanks for the reply.
I tried doing that but got stuck already in the first step. I couldn't find the solution in momentumspace therefore I tried transforming the whole equation.
Can you explain how you solve this equation? Do I have to guess a solution?
 
Hi AwesomeTrains, I see what you mean I started the solution and ran into the same problem. I tried another approach that I am not 100% sure is correct although it does give me an answer for the force. Basically try substititing in the equation for the furier transform instead of you wave function. You can then express your formula as
\int(\frac{\vec{p}^2}{2m}-a\nabla_\vec{p}^2)[\psi(\vec{r},t)e^{i\vec{p}\cdotp\vec{r}}]d\vec{r}^3=\int i\frac{\partial}{\partial t}\psi(\vec{r},t)e^{i\vec{p}\cdotp\vec{r}}d\vec{r}^3
then substitute in the hamiltonian for ##i\frac{\partial}{\partial t}\psi(\vec{r},t)## and compare like therms in your integral.
 
Brage said:
Hi AwesomeTrains, I see what you mean I started the solution and ran into the same problem. I tried another approach that I am not 100% sure is correct although it does give me an answer for the force. Basically try substititing in the equation for the furier transform instead of you wave function. You can then express your formula as
\int(\frac{\vec{p}^2}{2m}-a\nabla_\vec{p}^2)[\psi(\vec{r},t)e^{i\vec{p}\cdotp\vec{r}}]d\vec{r}^3=\int i\frac{\partial}{\partial t}\psi(\vec{r},t)e^{i\vec{p}\cdotp\vec{r}}d\vec{r}^3
then substitute in the hamiltonian for ##i\frac{\partial}{\partial t}\psi(\vec{r},t)## and compare like therms in your integral.
Thanks a lot for helping me out :biggrin:
Okay I'll try and continue from where you stopped:
\int(\frac{\vec{p}^2}{2m}-a\nabla_\vec{p}^2)[\psi(\vec{r},t)e^{i\vec{p}\cdotp\vec{r}}]d\vec{r}^3=\int i\frac{\partial}{\partial t}\psi(\vec{r},t)e^{i\vec{p}\cdotp\vec{r}}d\vec{r}^3
Then from the first equation we have that (Dividing by \psi(\vec{r},t)):
(\frac{\vec{p}^2}{2m}-a\nabla_\vec{p}^2)=i\frac{\partial}{\partial t}
Plugging it in:
\int(\frac{\vec{p}^2}{2m}-a\nabla_\vec{p}^2)[\psi(\vec{r},t)e^{i\vec{p}\cdotp\vec{r}}]d\vec{r}^3=\int (\frac{\vec{p}^2}{2m}-a\nabla_\vec{p}^2)\psi(\vec{r},t)e^{i\vec{p}\cdotp\vec{r}}d\vec{r}^3

How do I continue from here? Did I plug in the wrong hamiltonian? :-p
If I would plug in the general one H=T+V I would by comparison, like you said, get that V(\vec{p})=a\nabla_\vec{p}^2 since the first term looks like the kinetic energy.

Do I then transform the potential into position space and then do the derivative to get the force? Is this how you got your result?

Kind regards
Alex
 
$$\int(\frac{\vec{p}^2}{2m}-a\nabla_\vec{p}^2)[\psi(\vec{r},t)e^{i\vec{p}\cdotp\vec{r}}]d\vec{r}^3=\int (\frac{\vec{p}^2}{2m}+a\vec{r}^2)\psi(\vec{r},t)e^{i\vec{p}\cdotp\vec{r}}d\vec{r}^3$$
This is because the laplacian of our momentum coordiantes acts on ##\psi(\vec{r},t)e^{i\vec{p}\cdotp\vec{r}}##, where the only term with ##\vec{p}## dependence is the exponential function. This then means
$$\int \left(\frac{\vec{p}^2}{2m}+a\vec{r}^2\right)\psi(\vec{r},t)e^{i\vec{p}\cdotp\vec{r}}d\vec{r}^3=\int i\frac{\partial}{\partial t}\psi(\vec{r},t)e^{i\vec{p}\cdotp\vec{r}}d\vec{r}^3$$
Also, note that in position space, the hamiltonian becomes
$$\left(\frac{\vec{p}^2}{2m}+V(\vec{r})\right)\psi(\vec{r},t)=i\frac{\partial}{\partial t}\psi(\vec{r},t)$$
This is the nesessary equation to use here, as we cannot use the hamiltonian for momentum space when the operators act on the wavefunction for position space.
 
  • Like
Likes AwesomeTrains
Brage said:
$$\int(\frac{\vec{p}^2}{2m}-a\nabla_\vec{p}^2)[\psi(\vec{r},t)e^{i\vec{p}\cdotp\vec{r}}]d\vec{r}^3=\int (\frac{\vec{p}^2}{2m}+a\vec{r}^2)\psi(\vec{r},t)e^{i\vec{p}\cdotp\vec{r}}d\vec{r}^3$$
This is because the laplacian of our momentum coordiantes acts on ##\psi(\vec{r},t)e^{i\vec{p}\cdotp\vec{r}}##, where the only term with ##\vec{p}## dependence is the exponential function. This then means
$$\int \left(\frac{\vec{p}^2}{2m}+a\vec{r}^2\right)\psi(\vec{r},t)e^{i\vec{p}\cdotp\vec{r}}d\vec{r}^3=\int i\frac{\partial}{\partial t}\psi(\vec{r},t)e^{i\vec{p}\cdotp\vec{r}}d\vec{r}^3$$
Also, note that in position space, the hamiltonian becomes
$$\left(\frac{\vec{p}^2}{2m}+V(\vec{r})\right)\psi(\vec{r},t)=i\frac{\partial}{\partial t}\psi(\vec{r},t)$$
This is the nesessary equation to use here, as we cannot use the hamiltonian for momentum space when the operators act on the wavefunction for position space.

Ah I see, then V(\vec{r})=a\vec{r}^2 and it follows that \vec{F}(\vec{r})=-2a\vec{r} is that what you got aswell? If it is then thanks for the help :)
 
Last edited:
Yeah that is what I got as well, glad I could help.
 
Back
Top