AbstractIn this study, the and the two equivalent / polarized Raman spectra of a LiNbO3 single crystal have been recorded and used as a benchmark test for the density functional theory (DFT) calculation of the longitudinal modes and of their Raman activity. The theoretical approach, based on periodic boundary conditions and a linear combination of atomic orbitals (LCAO), provides excellent predictions of phonon wavenumbers and relative bands intensities for both A1 and E Longitudinal Optical (LO) modes and complements a previous paper limited to the study of Transverse Optical (TO) modes. Overall, the present investigation demonstrates that the LCAO approach, as implemented in the CRYSTAL software, gives results of similar accuracy for the TO and the LO phonons features. By means of a band deconvolution scheme applied to the experimental spectra, we present, for the first time, a quantitative comparison between experimental and theoretically predicted polarized Raman band intensities of LiNbO3 LO modes. This analysis highlights the role of the suitable determination of the static and high‐frequency dielectric matrices that are needed for the prediction of the TO/LO frequency split but also the first nonlinear electric susceptibility tensor for an accurate description of the Raman intensity pattern.