VASP plugins: linking the Viennaab-initioSimulation Package with Python

Implementing novel features and experimental algorithms into widely adopted density functional theory (DFT) codes is frequently hindered by complex legacy architectures and the use of compiled languages such as Fortran. These production codes, while optimized for high-performance computing clusters, present significant hurdles for software development and rapid prototyping, often requiring deep expertise in the code's internal structure to modify. To address this challenge, we present a Python plugin infrastructure for the Viennaab-initioSimulation Package (VASP) that combines computational efficiency with the flexibility of high-level scripting. Our architecture uses a C++ intermediate layer andto expose VASP data asarrays via shared memory buffers, ensuring high performance without data duplication. We implement two categories of plugins: those that modify quantities at the end of each converged self-consistent field (SCF) cycle, such as,and, and those that operate during the SCF cycle, such as. We demonstrate the utility of our implementation through three applications, structure relaxation using thelibrary, implementing an implicit solvent model, and adding thedispersion corrections. For a representative post-SCF plugin (), this coupling introduces a callback overhead of well underper ionic step relative to a natively linked implementation. This infrastructure effectively bridges the gap between high-performance electronic structure routines and the widespread scientific Python ecosystem.