Traditional computational methods for linear spectroscopy simulation have predominantly focused on incorporating vibronic effects through Franck-Condon (FC) spectral simulations 1-17 using either an adiabatic Hessian or a vertical gradient approach for efficiency. In practice, harmonic potential energy surfaces are generally assumed within the FC approach as the harmonic oscillator nuclear wave functions are readily available upon geometry optimization and normal mode analysis. The FC methods often model solvent effects implicitly, utilizing polarizable continuum models to approximate the influence of the surrounding medium. This harmonic FC approach works a) cisborn@ucmerced.edu well for rigid molecules or systems with weak solvent interactions. However, the absence of direct solute-solvent interactions and the harmonic treatment of the potential energy surfaces can lead to more narrow spectra than observed experimentally. Explicit solvent effects are traditionally incorporated into spectral simulations by employing a nuclear ensemble approach. This method involves sampling chromophore-solvent configurations, generally through a molecular dynamics (MD) trajectory, and naturally captures vibrational degrees of freedom of the chromophore. This approach does not incorporate vibronic effects, often leading to poor spectral shape. The ensemble approach has seen limited application in simulating fluorescence spectra, likely due to the computational demands of performing excited state molecular dynamics to obtain chromophore-solvent configurations on the excited state potential energy surface. There are multiple advantages of the nuclear ensemble approach, such as being able to sample anharmonic regions of the potential, the ease of including effects from higher-lying excited states, capturing inhomogeneous broadening due to explicit interactions with solvent, and accounting for temperature effects of both the chromophore and solvent via finite temperature MD. Some recent approaches to simulating absorption spectra have made progress in capturing both specific environmental and vibronic effects through the use of MD sampling of chromophore-environment configurations. For example, employing a cumulant expansion for the ensemble average of the time-ordered exponential of the linear response function treats chromophore vibrational and environmental degrees of freedom on equal footing, allowing vibronic absorption spectra to be generated from energy gap correlation function computed along an MD trajectory. Molecular dynamics, along with snapshot clustering for unique features, has also been exploited within a variational and perturbative approach to spectroscopy. An alternative approach based on the separation of classical (soft) and quantum (stiff) nuclear degrees of freedom expresses the absorption spectrum as a conformational integral of vibronic spectra of the stiff coordinates. The approach developed by some of the authors of the present work is similar in spirit, where vibronic spectra of the chromophore are generated within a frozen environment in order to combine features of both the ensemble and Franck-Condon approaches, creating hybrid ensemble-Franck-Condon (E-FC) approaches that capture both explicit environment and vibronic effects. In previous work, we have outlined variations of the E-FC approach with differing levels of approximation and corresponding computational cost, and shown good agreement with experiment for the simulated absorption spectral shapes of chromophores in various explicit solvent environments.