Practical Implementation: Continuity and The Minimum Image Convention
To implement periodic boundary conditions in practice, at least two steps are needed.
The first is to make an object which leaves the simulation cell on one side enter back on the other. This is of course a simple operation, and could in code be e.g. (for the x dimension, assuming an orthogonal unit cell centered on the origin):
if (periodicx) then if (x < -xsize*0.5) x=x+xsize if (x >= xsize*0.5) x=x-xsize endifThe second is to make sure that every distance between atoms, or other vector calculated from one atom to another, has a length and direction which corresponds to the minimum image criterion. This can be achieved as follows to calculate e.g. the x direction distance component from atom i to atom j:
if (periodicx) then dx = x(j) - x(i) if (abs(dx) > xsize*0.5) dx = dx - sign(xsize,dx) endifNaturally both operations should be repeated in all 3 dimensions.
These operations can be written in much more compact form for orthorhombic cells if the origin is shifted to a corner of the box. Then we have, in one dimension, for positions and distances respectively:
! After x(i) update without regard to PBC: x(i)=x(i)-floor(x(i)/xsize)*xsize !For a box with the origin at the lower left vertex ! Works for xs lying in any image. dx=x(j)-x(i) dx=dx-nint(dx/(0.5*xsize))*xsizeFor non-orthorhombic cells the situation can be considerably more complicated.
In simulations of ionic systems considerably more complicated operations may be needed to handle the long-range Coulomb interactions.
Read more about this topic: Periodic Boundary Conditions
Famous quotes containing the words practical, continuity, minimum, image and/or convention:
“The best way to teach a child restraint and generosity is to be a model of those qualities yourself. If your child sees that you want a particular item but refrain from buying it, either because it isn’t practical or because you can’t afford it, he will begin to understand restraint. Likewise, if you donate books or clothing to charity, take him with you to distribute the items to teach him about generosity.”
—Lawrence Balter (20th century)
“Continuous eloquence wearies.... Grandeur must be abandoned to be appreciated. Continuity in everything is unpleasant. Cold is agreeable, that we may get warm.”
—Blaise Pascal (1623–1662)
“After decades of unappreciated drudgery, American women just don’t do housework any more—that is, beyond the minimum that is required in order to clear a path from the bedroom to the front door so they can get off to work in the mourning.”
—Barbara Ehrenreich (20th century)
“We have no participation in Being, because all human nature is ever midway between being born and dying, giving off only a vague image and shadow of itself, and a weak and uncertain opinion. And if you chance to fix your thoughts on trying to grasp its essence, it would be neither more nor less than if your tried to clutch water.”
—Michel de Montaigne (1533–1592)
“The metaphor of the king as the shepherd of his people goes back to ancient Egypt. Perhaps the use of this particular convention is due to the fact that, being stupid, affectionate, gregarious, and easily stampeded, the societies formed by sheep are most like human ones.”
—Northrop Frye (b. 1912)