The nonlinear optical susceptibilities are analyzed in their optical propagation length considering the incoherent-coherent contributions with the dependence on the relaxation times. In addition, some symmetry properties of the Four-Wave mixing (FWM) signal concerning the frequency detuning exchange are used in the study, without restricting ourselves to the maxima of population oscillations. We analyze the behaviors of these optical responses in a vibronic coupling scheme using a molecule model consisting of two coupled harmonic curves of electronic energies with displaced minima in position and nuclear energies. The solvent effect in our model is treated through the natural Bohr frequency shift to a time-dependent function, with explicit manifestations in its comparison as if the upper state were broadened. Our results for both FWM signal propagation and nonlinear optical susceptibilities are sensitive to the intramolecular coupling parameters, solvent stochasticity, and perturbation order for the treatment of the weak probe beam. The consideration of a second-order probe beam in FWM would not represent a significant contribution to the predicted results, since its percentage contribution is negligible within the model development framework; however, this contribution becomes important when increasing the ratio between longitudinal and transverse relaxation times.