Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

A Monte Carlo Wavefunction Methods

  1. Mar 26, 2016 #1
    From the Theory of Open Quantum Systems; the Euler scheme is given by:

    ##\psi_{k+1} = \psi_{k} + D_1(\psi_k)\Delta t + D_2(\psi_k) \Delta W_k##

    and is a scheme of order 1. What does the order of convergence mean? From my understanding higher order schemes require fewer interations to give a decent approximation. Is there anything more than that?
    Last edited: Mar 26, 2016
  2. jcsd
  3. Mar 31, 2016 #2
    Thanks for the post! This is an automated courtesy bump. Sorry you aren't generating responses at the moment. Do you have any further information, come to any new conclusions or is it possible to reword the post?
  4. Apr 1, 2016 #3

    A. Neumaier

    User Avatar
    Science Advisor

    Order ##p## means that the error per step is ##O(\Delta t^{p+1})##.

    As a consequence, for a nonstiff differential equation, the total error after intergation with a single-step method over a time of order 1 is ##O(\Delta t^{p})##. If ##\Delta t## is sufficiently small this makes higher order methods far more useful, in the sense that to get the same global accuracy, far bigger (and hence far fewer) steps can be taken.

    For multistep methods and in the stiff case, the total error may be much larger, and additional stability conditions beyond the order are needed to get the same conclusion.

    The possible conclusions about the global error are again more complicated (and not easy to summarize) in case of using these methods for stochastic differential equations, as in your case. (But all this has nothing to do with physics.)
    Last edited: Apr 1, 2016
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook

Have something to add?
Draft saved Draft deleted