MathGroup Archive 2000

[Date Index] [Thread Index] [Author Index]

Search the Archive

want to modify NDSolve--molecular dynamics with mathematica


Hello,

I was wondering whether I can apply NDSolve to molecular dynamics, on the
standpoint of periodical boundary condition.

The eqation I want to solve is the simple Newton equation (Let me simplify
the problem).
  m D[x, {x,2}] == F
However, periodical boundary condition makes it difficult to apply NDSolve.
That is, once
  0<x[t]<a0
is not satisfied (a0 is a constant), I have to move the atom so that this
condition is satisfied.
In concrete,
  if x[t]<0  --->  x[t]=x[t]+a0
  if x[t]>a0  --->  x[t]=x[t]-a0
and then restart to solve the equation.

Is there any way to realize this by directly applying NDSolve?
If not, how can I modify NDSolve?

-Toshi



  • Prev by Date: Re: Mathlink probs under Unix
  • Next by Date: Re: Mathlink probs under Unix
  • Previous by thread: RE: Relative Font Size in Plots / Graphics`
  • Next by thread: Re: want to modify NDSolve--molecular dynamics with mathematica