Stuck on a complex ODE (conjugate issues )

Simplexed
Messages
2
Reaction score
0
Hi everyone.

I'm stuck on a problem, banging of head etc. Basically, I've got a time-evolution problem restricted to the positive quadrant of complex space, where the state of the system \psi \in \mathbb{C} is described by the following type of ODE:

<br /> \frac{{d\psi \left( t \right)}}{{dt}} = i\overline {\psi \left( t \right)} <br />

i.e, the time evolution of the state depends on the conjugate, not the state itself (this is determined by the first principles of my system).

Numerically, I can easily generate the quiver diagrams of this first difference, and then use it to simulate the time evolution of the system. The norm seems to grow roughly exponentially, while the argument exponentially approaches \pi /4 , which is what I expect.

My problem is that I am at a loss of ideas for how to actually integrate the expression, in order to obtain the expression for the time evolution of \psi \left( 0 \right) for any given starting point in the space.

Any assistance on this would be very much appreciated.
 
Last edited:
Physics news on Phys.org
This is just a pair of coupled DE's. Two complex numbers are equal if and only if their real parts are equal and their imaginary parts are equal. Your equation can be written

\frac{d}{dt}(\psi_R +i\psi_I) = -\psi_I + i\psi_R

which splits up into

\frac{d \psi_R}{dt}= -\psi_I

and

\frac{d \psi_I}{dt} = \psi_R

which you should hopefully be able to solve, yes?
 
:redface:

Well, that figures... I was trying to work it while keeping the complex number whole, but completely overlooked treating it as a system of ODEs (and not a complicated one at that).

Thanks a lot for that wake-up call!
 
In this case, this is likely the easier way to solve it: You could also note that if you take the conjugate of your original DE you get

\frac{d \overline{\psi}}{dt} = -i\psi

so taking another derivative of your original equation gives

\frac{d^2\psi}{dt^2} = i\frac{d \overline{\psi}}{dt} = \psi
 
Last edited:
There is the following linear Volterra equation of the second kind $$ y(x)+\int_{0}^{x} K(x-s) y(s)\,{\rm d}s = 1 $$ with kernel $$ K(x-s) = 1 - 4 \sum_{n=1}^{\infty} \dfrac{1}{\lambda_n^2} e^{-\beta \lambda_n^2 (x-s)} $$ where $y(0)=1$, $\beta>0$ and $\lambda_n$ is the $n$-th positive root of the equation $J_0(x)=0$ (here $n$ is a natural number that numbers these positive roots in the order of increasing their values), $J_0(x)$ is the Bessel function of the first kind of zero order. I...
Are there any good visualization tutorials, written or video, that show graphically how separation of variables works? I particularly have the time-independent Schrodinger Equation in mind. There are hundreds of demonstrations out there which essentially distill to copies of one another. However I am trying to visualize in my mind how this process looks graphically - for example plotting t on one axis and x on the other for f(x,t). I have seen other good visual representations of...
Back
Top