10.26165/JUELICH-DATA/8WA9I7 Wittberg, Clemens(PGI-7, Forschungszentrum Jülich GmbH) Jülich DATA 2026 The corresponding git repository can be found at: https://iffgit.fz-juelich.de/pgi-7/band-alignment-calculations This repository contains the numerical model used to calculate the equilibrium and bias-dependent electrostatic potential, electric field, charge density, and conduction/valence band profiles across arbitrarily complex, multi-layer (metal–insulator–semiconductor) heterostructures. The model self-consistently solves the Poisson equation as a two-point boundary value problem (BVP), including: Position-dependent doping/trap concentrations, permittivities, effective masses, and electron affinities across an arbitrary number of layers. A Gaussian-distributed trap density of states (DOS) per material, with configurable mean depth, standard deviation, and relative weight for multiple trap species. Fermi–Dirac occupation statistics (with spin degeneracy for donor-like traps) for both trapped charge and free carriers in the conduction band (parabolic band approximation). Externally applied bias, including additional voltage-drop terms to capture non-ideal voltage division across the stack (e.g., due to interface dipoles or series resistances). The output of the model is the spatial profile of the conduction and valence band edges, which can be used to interpret band alignment, barrier formation, and accumulation/depletion behavior at oxide interfaces relevant to resistive-switching (memristive) devices. Wittberg, Clemens(PGI-7, Forschungszentrum Jülich GmbH)Deckers, Malte(PGI-7, Forschungszentrum Jülich GmbH)Dittmann, Regina(PGI-7, Forschungszentrum Jülich GmbH)