<i>Ab initio</i> simulations of EMI-BF4 neutral-surface interactions in electrospray thrusters

Electrospray thrusters promise compact, high specific impulse propulsion for small spacecraft, yet ground characterization remains confounded by secondary species emission and incomplete diagnostics of neutral products. To address these limitations, we perform energy-resolved mixed quantum/classical ab initio molecular dynamics of neutral 1-ethyl-3-methylimidazolium tetrafluoroborate colliding with Au extractor surfaces with impact energies from 10 to 100 eV to resolve fragment species spectra, charge states, kinetic energy partitioning, and scattering geometry. The simulations reveal a three-stage sequence with impact energy: the low energy regime, 10–20 eV, which favors ionic dissociation; intermediate energy regime, between 30 and 40 eV, opens a neutralization window; and high energy regime, greater than 50 eV, drives covalent fragmentation into many light products with mixed charge states. Fractional energy distributions show a transition from few-body, energy-concentrated outcomes in the low energy regime to many-body, energy-dispersed outcomes in the high energy regime. Deflection angle distributions exhibit a strong mass-to-angle anti-correlation such that heavier fragments favor small deflection, while lighter fragments are observed across the full deflection range but dominate larger deflection angles. The fraction of metastables peaks near 50 eV, coinciding with abundant neutral fragment production. Importantly, neutral bombardment still produces charged secondaries at the target even when the upstream ion plume is fully suppressed by a decelerating electrode. These findings provide a basis for de-biasing facility measurements by pairing tandem time-of-flight secondary ion mass spectrometry and a residual gas analyzer with suppression-bias corrections to inform the design of electrospray thrusters that reduce interception on extractor surfaces.