Enabling Ab-Initio Molecular Dynamics under Bias: The CP2K+SMEAGOL Interface for Integrating Density Functional Theory and Non-Equilibrium Green Functions
arXiv:2406.11494 · doi:10.1021/acs.jctc.4c00371
Abstract
Density functional theory (DFT) combined with non-equilibrium Greens functions (NEGF) is a powerful approach to model quantum transport under external bias potentials, at reasonable computational cost. In this work we present a new interface between the popular mixed Gaussian/plane wave electronic structure package CP2K and the NEGF code SMEAGOL, the most feature-rich implementation of DFT-NEGF available for CP2K to-date. The CP2K+SMEAGOL interface includes the implementation of current induced forces. We verify this implementation for a variety of systems: an infinite 1D Au wire, a parallel-plate capacitor and a Au-H2-Au junction. We find good agreement with SMEAGOL calculations performed with SIESTA for the same systems, and with the example of a solvated Au wire demonstrate for the first time that DFT-NEGF can be used to perform molecular dynamics simulations under bias of large-scale condensed phase systems under realistic operating conditions.
References in corpus (22)
- The SIESTA method for ab initio order-N materials simulation
- Efficient index handling of multidimensional periodic boundary conditions
- Density functional method for nonequilibrium electron transport
- CP2K: An Electronic Structure and Molecular Dynamics Software Package -- Quickstep: Efficient and Accurate Electronic Structure Calculations
- QuantumATK: An integrated platform of electronic and atomic-scale modelling tools
- Formation and manipulation of a metallic wire of single gold atoms
- First-Principles Based Matrix-Green's Function Approach to Molecular Electronic Devices: General Formalism
- SIESTA: recent developments and applications
- Improvements on non-equilibrium and transport Green function techniques: the next-generation transiesta
- Advances and challenges in single-molecule electron transport
- Stiff monatomic gold wires with a spinning zigzag geometry
- Efficient implementation of the nonequilibrium Green function method for electronic transport calculations
- Calibration of the length of a chain of single gold atoms
- Effects of self-interaction corrections on the transport properties of phenyl-based molecular junctions
- The Puzzling Stability of Monatomic Gold Wires
- H dissociation over Au-nanowires and the fractional conductance quantum
- Current-induced forces: a simple derivation
- Ab-initio non-equilibrium quantum transport and forces with the real space projector augmented wave method
- Pushing Back the Limit of Ab-initio Quantum Transport Simulations on Hybrid Supercomputers
- Current-induced energy barrier suppression for electromigration from first principles
- First-principles methodology for quantum transport in multiterminal junctions
- Inelastic fingerprints of hydrogen contamination in atomic gold wire systems