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

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.