File: ea.py

package info (click to toggle)
python-ase 3.26.0-3
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid
  • size: 15,484 kB
  • sloc: python: 148,112; xml: 2,728; makefile: 110; javascript: 47
file content (13 lines) | stat: -rw-r--r-- 332 bytes parent folder | download 
1
2
3
4
5
6
7
8
9
10
11
12
13
 
from ase.db import connect

refs = connect('refs.db')
db = connect('ads.db')

for row in db.select():
    ea = (
        row.energy
        - refs.get(formula=row.ads).energy
        - refs.get(layers=row.layers, surf=row.surf).energy
    )
    h = row.positions[-1, 2] - row.positions[-2, 2]
    db.update(row.id, height=h, ea=ea)