Motivated by the outstanding properties of bis(BF2) core complexes as fluorophore probes, we present a systematic computational study of their vibrationally resolved one- and two-photon absorption spectra in vacuum and in solution. Electronic and vibrational structure calculations were performed using the coupled cluster CC2 method and the Kohn-Sham formulation of density functional theory (DFT). A nonempirical estimation of the inhomogeneous broadening, accomplished using the polarizable embedding (PE) approaches combined with time-dependent DFT and CC2 methods, is used as a key ingredient of the computational protocol employed for simulations of the spectral features in solution. The inhomogeneous broadening is also determined based on the Marcus theory employing linear response and state-specific polarizable continuum model (PCM) methods. It is found that the polarizable embedding CC2 and the state-specific PCM are the most successful approaches for description of environmental broadening. For the 1(1)A(g) -> 1(1)B(u) transition, the non-Condon effects can be safely neglected and a fair agreement between the simulated and experimental band shapes is found. In contrast, the shape of the vibrationally resolved band corresponding to the two-photon allowed 1(1)A(g) -> 2(1)A(g) transition is largely dominated by non-Condon effects. A generalized few-level model was also employed to analyze the mechanism of the electronic two-photon 1(1)A(g) -> 2(1)A(g) excitation. It was found that the most important optical channel involves the 1(1)B(u) excited state. Ramifications of the findings for general band shape modeling are briefly discussed.