root / ase / test / maxwellboltzmann.py @ 13
Historique | Voir | Annoter | Télécharger (395 octet)
| 1 | 1 | tkerber | from ase.md.velocitydistribution import MaxwellBoltzmannDistribution |
|---|---|---|---|
| 2 | 1 | tkerber | from ase.lattice.cubic import FaceCenteredCubic |
| 3 | 1 | tkerber | |
| 4 | 1 | tkerber | atoms = FaceCenteredCubic(size=(50,50,50), symbol="Cu", pbc=False) |
| 5 | 1 | tkerber | print "Number of atoms:", len(atoms) |
| 6 | 1 | tkerber | MaxwellBoltzmannDistribution(atoms, 0.1)
|
| 7 | 1 | tkerber | temp = atoms.get_kinetic_energy() / (1.5 * len(atoms)) |
| 8 | 1 | tkerber | |
| 9 | 1 | tkerber | print "Temperature", temp, " (should be 0.1)" |
| 10 | 1 | tkerber | assert abs(temp - 0.1) < 1e-3 |