Resonantly Scattered CIV Emission in Local Star-forming Galaxies: Radiative Transfer Constraints on High-ionization Gas in Reionization Analogs

Rest-frame ultraviolet emission lines are diagnostics of hard ionizing spectra and highly ionized gas in galaxies resembling the sources of cosmic reionization. Nebular CIV emission traces both intrinsic line production and resonant scattering through its velocity-resolved profile. Previous studies have treated CIV empirically, using integrated fluxes, equivalent widths, and line ratios rather than modeling the radiative transfer (RT) shaping the emergent profiles. We present the first systematic RT modeling of resonantly scattered CIV emission profiles in 18 local star-forming galaxies observed with HST/COS. Using the clumpy RT framework PEACOCK, idealized experiments show that profile morphology is governed by column density, gas kinematics, intrinsic equivalent width, and aperture-dependent recovery of scattered emission. We fit continuum-normalized spectra corrected for stellar contributions with a neural-network-accelerated Bayesian pipeline, reproducing P-Cygni-like, double-peaked, double-peaked with two absorption troughs, and blue-bump plus red-peak morphologies. Emission infilling connects these profile classes and complicates kinematic inferences, while aperture losses decouple observed net equivalent widths from intrinsic line production. Combining RT-inferred intrinsic equivalent widths with rest-UV diagnostics, we find that most galaxies are consistent with stellar photoionization, although some may require harder spectra or more extreme ionization conditions. Comparison with CLASSY reveals stronger associations of stellar mass and star formation rate with gas kinematics than with CIV column density. Our results establish CIV as both a tracer of hard ionizing radiation and a resonant-line probe of gas structure and kinematics, with velocity-resolved profiles providing information beyond integrated line strengths in local reionization analogs.