Investigating anharmonicities in polarization-orientation Raman spectra of acene crystals with machine learning

We present a first-principles machine-learning computational framework to investigate anharmonic effects in polarization-orientation (PO) Raman spectra of molecular crystals, focusing on anthracene and naphthalene. By combining machine learning models for interatomic potentials and polarizability tensors, we enable efficient, large-scale simulations that capture temperature-dependent vibrational dynamics beyond the harmonic approximation. Our approach reproduces key qualitative features observed experimentally. We show, systematically, what are the fingerprints of anharmonic lattice dynamics, thermal expansion, and Raman tensor symmetries on PO-Raman intensities. However, we find that the simulated polarization dependence of Raman intensities shows only subtle deviations from quasiharmonic predictions, failing to capture the pronounced temperature-dependent changes that have been reported experimentally in anthracene. We propose that part of these inconsistencies stem from the impossibility to deconvolute certain vibrational peaks when only experimental data is available. This work therefore provides a foundation to improve the interpretation of PO-Raman experiments in complex molecular crystals with the aid of theoretical simulations.