[ase-users] how to get the pressure

andreasm at fysik.dtu.dk andreasm at fysik.dtu.dk
Wed Apr 20 17:41:43 CEST 2011


Hi,

Your missing an a.

pressure=get_isotropic_pressure()
should be
pressure=a.get_isotropic_pressure()


Best wishes
Andreas


> hi ,I have some  questions about ASE.how can i get the pressure of system
> ,such as h2.
> my scribe is following:
>
> from ase import Atoms
> from gpaw import GPAW
> a = Atoms('H2', positions=[(0, 0, 0), (0, 0, 1.1)],
>            calculator=GPAW())
> stress=a.get_stress()
> pressure=get_isotropic_pressure()
> print stress ,pressure
> but it does not work.
> what should i do?
> print h2.get_forces()
>  _______________________________________________
> ase-users mailing list
> ase-users at listserv.fysik.dtu.dk
> https://listserv.fysik.dtu.dk/mailman/listinfo/ase-users





More information about the ase-users mailing list