root / ase / optimize / test / C5H12.py @ 14
Historique | Voir | Annoter | Télécharger (2,05 ko)
| 1 |
#/usr/bin/env python
|
|---|---|
| 2 |
#PBS -l nodes=4:ppn=8
|
| 3 |
#PBS -l walltime=13:00:00
|
| 4 |
from ase import Atoms |
| 5 |
from ase.optimize.test import run_test |
| 6 |
from gpaw import GPAW |
| 7 |
from gpaw import Mixer |
| 8 |
from gpaw.poisson import PoissonSolver |
| 9 |
|
| 10 |
name = 'C5H12'
|
| 11 |
|
| 12 |
def get_atoms(): |
| 13 |
atoms = Atoms(symbols='C5H12',
|
| 14 |
pbc=[False, False, False], |
| 15 |
cell=[ |
| 16 |
[ 16.83752497, 0. , 0. ], |
| 17 |
[ 0. , 12.18645905, 0. ], |
| 18 |
[ 0. , 0. , 11.83462179] |
| 19 |
], |
| 20 |
positions=[ |
| 21 |
[ 5.90380523, 5.65545388, 5.91569796], |
| 22 |
[ 7.15617518, 6.52907738, 5.91569796], |
| 23 |
[ 8.41815022, 5.66384716, 5.92196554], |
| 24 |
[ 9.68108996, 6.52891016, 5.91022362], |
| 25 |
[ 10.93006206, 5.65545388, 5.91569796], |
| 26 |
[ 5.00000011, 6.30002353, 5.9163716 ], |
| 27 |
[ 5.88571848, 5.0122839 , 6.82246859], |
| 28 |
[ 5.88625613, 5.01308931, 5.01214155], |
| 29 |
[ 7.14329342, 7.18115393, 6.81640316], |
| 30 |
[ 7.14551332, 7.17200869, 5.00879027], |
| 31 |
[ 8.41609966, 5.00661165, 5.02355167], |
| 32 |
[ 8.41971183, 5.0251482 , 6.83462168], |
| 33 |
[ 9.69568096, 7.18645894, 6.8078633 ], |
| 34 |
[ 9.68914668, 7.16663649, 5.00000011], |
| 35 |
[ 10.95518898, 5.02163182, 6.8289018 ], |
| 36 |
[ 11.83752486, 6.29836826, 5.90274952], |
| 37 |
[ 10.94464142, 5.00000011, 5.01802495] |
| 38 |
]) |
| 39 |
return atoms
|
| 40 |
|
| 41 |
def get_calculator(): |
| 42 |
calc = GPAW(h=0.2,
|
| 43 |
mode = 'lcao',
|
| 44 |
basis = 'szp(dzp)',
|
| 45 |
mixer=Mixer(beta=0.1, nmaxold=5, weight=50.0), |
| 46 |
poissonsolver=PoissonSolver(nn='M', relax='GS'), |
| 47 |
txt='C5H12.txt')
|
| 48 |
return calc
|
| 49 |
|
| 50 |
run_test(get_atoms, get_calculator, name + '-gpaw')
|