Replication Data for: Structure-Preserving Real-Space Quantum Simulation of Maxwell Fields in Inhomogeneous Photonic Structures

Quantum-optical simulations of inhomogeneous media typically rely on a global eigenmode basis that must be recomputed for each geometry. We show that, for lossless, non-dispersive media, the dynamics can instead be formulated without such a basis. From the source-free Maxwell equations we identify a self-adjoint real-space generator of single-photon dynamics and directly construct bosonic field operators, a single-particle Hamiltonian, and a Liouville-von Neumann equation for the single-photon density matrix. A finite integration technique (FIT) discretization transfers this operator structure to the numerical grid while preserving self-adjointness. This ensures unitary evolution and exact photon-number conservation, up to projection accuracy and machine precision. Numerical examples demonstrate single-photon propagation in a dielectric waveguide and reproduce the two-photon Hong–Ou–Mandel dip in a 50/50 coupler. The few-photon dynamics in this formulation is fully determined by the bosonic symmetry of the state space and the linear single-particle propagation, without explicit photon-photon interaction terms. The approach therefore enables direct quantum-optical simulation of inhomogeneous and geometrically complex photonic structures. This dataset includes the program files along with the results generated with them.