How to relax only neighboring atoms  [SOLVED]

Total energy, geometry optimization, DFT+U, spin....

Moderator: bguster

Locked
mboaks
Posts: 10
Joined: Thu Feb 16, 2017 6:55 pm

How to relax only neighboring atoms

Post by mboaks » Wed Apr 05, 2017 8:58 pm

Hi,

Could anyone let me know if it is possible to relax only part of a structure? For example, I have a 4x4x4 Si supercell with 1/4th of a pore (diameter=3 nm) in a corner. I need to relax only the Si atoms near the pore. Is it possible in abinit? I have all the atom positions in Angstrom unit.

Thanks in advance.

Regards,
Mawla Boaks

ilukacevic
Posts: 271
Joined: Sat Jan 16, 2010 12:05 pm
Location: Dept. of Physics, University J. J. Strossmayer, Osijek, Croatia
Contact:

Re: How to relax only neighboring atoms  [SOLVED]

Post by ilukacevic » Thu Apr 06, 2017 8:36 am

Dear Mawla,

You can use 'natfix' and 'iatfix' variables to specify exactly which atoms you want to fix. The ones you do not specify will be relaxed.

Best regards,

Igor

mboaks
Posts: 10
Joined: Thu Feb 16, 2017 6:55 pm

Re: How to relax only neighboring atoms

Post by mboaks » Sat Apr 08, 2017 10:40 pm

Thank you so much. It totally works!

Locked