A Consistent Implementation of Cluster Strong Lensing in Cosmological Simulation Light Cones

Galaxy cluster strong gravitational lensing plays a central role in precision cosmology, yet robust theoretical predictions have lagged behind an abundance of high-quality strong lensing observations. This shortfall reflects both a mismatch between the geometry of the strong-lensing problem and standard cubic simulation boxes, and the fundamental tension between simulation volume and resolution. Consequently, many current forecasts adopt hybrid approaches that extract individual lenses from simulations and combine them with analytic or observed source populations positioned near caustics. These methods often omit correlated and/or uncorrelated line-of-sight (LoS) structure, or include it in ways that do not preserve correlations across redshift. Here we present a fully simulation-based procedure that generates strong-lensing images directly from particle data, drawing the lens, source, and all intervening resolved objects self-consistently from the simulated large-scale structure. Our approach combines a structure-preserving remapping of the simulation volume into a lensing-appropriate geometry with multi-plane ray tracing, enabling the use of uniform simulation boxes that resolve both cluster-scale primary lenses and high-redshift source galaxies. We demonstrate the method by generating example light cones and images using IllustrisTNG data, then use these results to conservatively quantify the impact of LoS structure on image configurations and critical-curve morphology. We find that uncorrelated LoS structure can shift the relative positions of lensed images by several arcseconds, introduces a $\sim 6\%$ scatter in the area of a cluster's primary critical curve, and changes the total critical area within 100$^{\prime\prime}$ of the cluster potential minimum by $16^{+20\%}_{-14\%}$ at a source plane redshift of $z_s=4$.