Source code for gnrs.energy.vasp

"""
This module computes the energy using VASP DFT code with ASE.

This source code is licensed under the BSD-3-Clause license found in the
LICENSE file in the root directory of this source tree.
"""
from __future__ import annotations

__author__ = ["Yi Yang", "Rithwik Tom"]
__email__ = "yiy5@andrew.cmu.edu"
__group__ = "https://www.noamarom.com/"

import subprocess
from ase import Atoms
from ase.calculators.vasp import Vasp

from gnrs.core.energy import EnergyCalculatorABC




class VASP(EnergyCalculatorABC):
    """
    Computes the energy using VASP DFT code with ASE
    https://wiki.fysik.dtu.dk/ase/ase/calculators/vasp.html
    """



    def __init__(self, *args) -> None:
        super().__init__(*args)
        if self.tsk_set.get("use_slurm", False):
            cmd = "scontrol show hostname $SLURM_JOB_NODELIST"
            all_hosts = subprocess.check_output(cmd, shell=True).decode().strip().split('\n')
            host = all_hosts[self.rank+1]
            command = f'mpirun -host {host} -np {self.tsk_set["num_cores"]} {self.tsk_set["command"]}'
        else:
            command = self.tsk_set["command"]
        self.calc = Vasp(command=command, **self.tsk_set['energy_settings'])




    def initialize(self) -> None:
        """
        Initialize the energy calculator.
        """
        pass




    def compute(self, xtal: Atoms) -> None:
        """
        Compute the energy of the crystal.
        """
        xtal.calc = self.calc
        try:
            energy = xtal.get_potential_energy()
        except Exception:
            energy = 0
        xtal.info[self.energy_name] = energy




    def finalize(self) -> None:
        """
        Finalize the energy calculator.
        """
        pass