Chaos in phase conjugation: physical vs numerical instabilities

Four-wave mixing equations in photorefractive media are approximated by different dynamical models and treated by different numerical methods. It is shown that the onset of instabilities and irregular behaviour in the same crystal, with a single wave mixing region, may be dependent both on the model used and the numerical method applied. Long-time irregular dynamics following from any finite-order difference schemes should be viewed with caution.