root / ase / test / dimer.py @ 4
Historique | Voir | Annoter | Télécharger (548 octet)
| 1 |
from ase import Atom, Atoms |
|---|---|
| 2 |
from ase.calculators.lj import LennardJones |
| 3 |
from ase.constraints import FixBondLength |
| 4 |
|
| 5 |
dimer = Atoms([Atom('X', (0, 0, 0)), |
| 6 |
Atom('X', (0, 0, 1))], |
| 7 |
calculator=LennardJones(), |
| 8 |
constraint=FixBondLength(0, 1)) |
| 9 |
print dimer.get_forces()
|
| 10 |
print dimer.positions
|
| 11 |
dimer.positions[:] += 0.1
|
| 12 |
print dimer.positions
|
| 13 |
dimer.positions[:, 2] += 5.1 |
| 14 |
print dimer.positions
|
| 15 |
dimer.positions[:] = [(1,2,3),(4,5,6)] |
| 16 |
print dimer.positions
|
| 17 |
dimer.set_positions([(1,2,3),(4,5,6.2)]) |
| 18 |
print dimer.positions
|
| 19 |
|