We computed the FCFs using the vertical gradient (VG) method as implemented in the ezFCF package . At the optimized EOM-EA-CCSD/cc-pVDZ geometry for the ground state, we calculated the harmonic frequencies and modes for the ground state. At the same geometry with the same method we then calculated the gradient for the excited state. VG approximation finds the equilibrium geometries for the excited state. We use that excited state geometry for the parallel approximation harmonic simulation of the FCFs.