Reuse & Permissions

It is not necessary to obtain permission to reuse this article or its components as it is available under the terms of the Creative Commons Attribution 4.0 International license. This license permits unrestricted use, distribution, and reproduction in any medium, provided attribution to the author(s) and the published article's title, journal citation, and DOI are maintained. Please note that some figures may have been included with permission from other third parties. It is your responsibility to obtain the proper permission from the rights holder directly for these figures.

Export citation

Export citation

Choose format for download:

Download Citation
  • Editors' Suggestion
  • Open Access

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

Paolo Lazzaroni, Shubham Sharma, and Mariana Rossi*

  • *Contact author: mariana.rossi@mpsd.mpg.de

Phys. Rev. B 113, 054309 – Published 20 February, 2026

DOI: https://doi.org/10.1103/t4s2-45t8

Abstract

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.

View figure in article

Physics Subject Headings (PhySH)

Article Text

Supplemental Material

References (84)

  1. N. Raimbault, V. Athavale, and M. Rossi, Anharmonic effects in the low-frequency vibrational modes of aspirin and paracetamol crystals, Phys. Rev. Mater. 3, 053605 (2019).
  2. M. Rossi, P. Gasparotto, and M. Ceriotti, Anharmonic and quantum fluctuations in molecular crystals: A first-principles study of the stability of paracetamol, Phys. Rev. Lett. 117, 115702 (2016).
  3. M. Krynski and M. Rossi, Efficient Gaussian process regression for prediction of molecular crystals harmonic free energies, npj Comput. Mater. 7, 169 (2021).
  4. V. A. Konstantinov, Heat transfer by low-frequency phonons and “diffusive” modes in molecular crystals, Low Temp. Phys. 29, 422 (2003).
  5. L. Legenstein, L. Reicht, S. Wieser, M. Simoncelli, and E. Zojer, Heat transport in crystalline organic semiconductors: Coexistence of phonon propagation and tunneling, npj Comput. Mater. 11, 29 (2025).
  6. S. Fratini, S. Ciuchi, D. Mayou, G. T. de Laissardière, and A. Troisi, A map of high-mobility molecular semiconductors, Nat. Mater. 16, 998 (2017).
  7. V. Coropceanu, J. Cornil, D. A. da Silva Filho, Y. Olivier, R. Silbey, and J.-L. Brédas, Charge transport in organic semiconductors, Chem. Rev. 107, 926 (2007).
  8. M. Asher, D. Angerer, R. Korobko, Y. Diskin-Posner, D. A. Egger, and O. Yaffe, Anharmonic lattice vibrations in small-molecule organic semiconductors, Adv. Mater. 32, 1908028 (2020).
  9. M. Asher, R. Jouclas, M. Bardini, Y. Diskin-Posner, N. Kahn, R. Korobko, A. R. Kennedy, L. Silva de Moraes, G. Schweicher, J. Liu, D. Beljonne, Y. Geerts, and O. Yaffe, Chemical modifications suppress anharmonic effects in the lattice dynamics of organic semiconductors, ACS Materials Au 2, 699 (2022).
  10. N. Benshalom, M. Asher, R. Jouclas, R. Korobko, G. Schweicher, J. Liu, Y. Geerts, O. Hellman, and O. Yaffe, Phonon–phonon interactions in the polarization dependence of Raman scattering, J. Phys. Chem. C 127, 18099 (2023).
  11. N. Raimbault, A. Grisafi, M. Ceriotti, and M. Rossi, Using Gaussian process regression to simulate the vibrational Raman spectra of molecular crystals, New J. Phys. 21, 105001 (2019).
  12. E. Berger and H.-P. Komsa, Polarizability models for simulations of finite temperature Raman spectra from machine learning molecular dynamics, Phys. Rev. Mater. 8, 043802 (2024).
  13. B. Gurlek, S. Sharma, P. Lazzaroni, A. Rubio, and M. Rossi, Accurate machine learning interatomic potentials for polyacene molecular crystals: Application to single molecule host-guest systems, npj Comput. Mater. 11, 318 (2025).
  14. M. Grumet, C. Von Scarpatetti, T. Bučko, and D. A. Egger, Delta machine learning for predicting dielectric properties and Raman spectra, J. Phys. Chem. C 128, 6464 (2024).
  15. A. Hashemi, A. V. Krasheninnikov, M. Puska, and H.-P. Komsa, Efficient method for calculating Raman spectra of solids with impurities and alloys and its application to two-dimensional transition metal dichalcogenides, Phys. Rev. Mater. 3, 023806 (2019).
  16. E. Berger, Z.-P. Lv, and H.-P. Komsa, Raman spectra of 2D titanium carbide MXene from machine-learning force field molecular dynamics, J. Mater. Chem. C 11, 1311 (2023).
  17. C. V. Raman and K. S. Krishnan, A new type of secondary radiation, Nature (London) 121, 501 (1928).
  18. J. Kim, J.-U. Lee, J. Lee, H. Ju Park, Z. Lee, C. Lee, and H. Cheong, Anomalous polarization dependence of Raman scattering and crystallographic orientation of black phosphorus, Nanoscale 7, 18708 (2015).
  19. B. Xu, N. Mao, Y. Zhao, L. Tong, and J. Zhang, Polarized Raman spectroscopy for determining crystallographic orientation of low-dimensional materials, J. Phys. Chem. Lett. 12, 7442 (2021).
  20. D. A. McQuarrie, Statistical Mechanics (University Scince Books, Mill Valley, CA, 2000).
  21. D. A. Long, The Raman Effect: A Unified Treatment of the Theory of Raman Scattering by Molecules (Wiley, New York, NY, 2002).
  22. M. Walter and M. Moseler, Ab initio wavelength-dependent Raman spectra: Placzek approximation and beyond, J. Chem. Theory Comput. 16, 576 (2020).
  23. H. A. Kramers and W. Heisenberg, Über die Streuung von Strahlung durch Atome, Z. Phys. 31, 681 (1925).
  24. P. A. M. Dirac and R. H. Fowler, The quantum theory of dispersion, Proc. R. Soc. London, Ser. A 114, 710 (1927).
  25. F. Knoop, N. Benshalom, M. Menahem, P. Gartner, T. Salzillo, O. Yaffe, and O. Hellman, Ab initio theory of the non-resonant Raman effect in crystals at finite temperature in comparison to experiment: The examples of GaN and BaZrS3, arXiv:2412.17711.
  26. See Supplemental Material at http://link.aps.org/supplemental/10.1103/t4s2-45t8 for details on normal modes projections, supercell convergence, perpendicular configuration PO-maps, Raman tensors benchmarking, and full harmonic PO-Raman results.
  27. B. A. DeAngelis, R. E. Newnham, and W. B. White, Factor group analysis of the vibrational spectra of crystals: A review and consolidation, Am. Mineral. 57, 255 (1972).
  28. M. H. Beck, A. Jäckle, G. A. Worth, and H. D. Meyer, The multiconfiguration time-dependent Hartree (MCTDH) method: A highly efficient algorithm for propagating wavepackets, Phys. Rep. 324, 1 (2000).
  29. S. Chatterjee and N. Makri, Real-time path integral methods, quantum master equations, and classical vs quantum memory, J. Phys. Chem. B 123, 10470 (2019).
  30. I. R. Craig and D. E. Manolopoulos, Quantum statistics and classical mechanics: Real time correlation functions from ring polymer molecular dynamics, J. Chem. Phys. 121, 3368 (2004).
  31. B. J. Braams and D. E. Manolopoulos, On the short-time limit of ring polymer molecular dynamics, J. Chem. Phys. 125, 124105 (2006).
  32. S. Habershon, D. E. Manolopoulos, T. E. Markland, and T. F. Miller, Ring-polymer molecular dynamics: Quantum effects in chemical dynamics from classical trajectories in an extended phase space, Annu. Rev. Phys. Chem. 64, 387 (2013).
  33. M. Rossi, M. Ceriotti, and D. E. Manolopoulos, How to remove the spurious resonances from ring polymer molecular dynamics, J. Chem. Phys. 140, 234116 (2014).
  34. S. C. Althorpe, Path integral simulations of condensed-phase vibrational spectroscopy, Annu. Rev. Phys. Chem. 75, 397 (2024).
  35. S. Jang and G. A. Voth, A derivation of centroid molecular dynamics and other approximate time evolution methods for path integral centroid variables, J. Chem. Phys. 111, 2371 (1999).
  36. G. Trenins, M. J. Willatt, and S. C. Althorpe, Path-integral dynamics of water using curvilinear centroids, J. Chem. Phys. 151, 054109 (2019).
  37. J. E. Lawrence, A. Z. Lieberherr, T. Fletcher, and D. E. Manolopoulos, Fast quasi-centroid molecular dynamics for water and ice, J. Phys. Chem. B 127, 9172 (2023).
  38. F. Musil, I. Zaporozhets, F. Noé, C. Clementi, and V. Kapil, Quantum dynamics using path integral coarse-graining, J. Chem. Phys. 157, 181102 (2022).
  39. T. E. Markland and M. Ceriotti, Nuclear quantum effects enter the mainstream, Nat. Rev. Chem. 2, 0109 (2018).
  40. J. Hoja, A. M. Reilly, and A. Tkatchenko, First-principles modeling of molecular crystals: Structures and stabilities, temperature and pressure, WIREs Comput. Mol. Sci. 7, e1294 (2017).
  41. A. Putrino and M. Parrinello, Anharmonic Raman spectra in high-pressure ice from ab initio simulations, Phys. Rev. Lett. 88, 176401 (2002).
  42. M. Pagliai, C. Cavazzoni, G. Cardini, G. Erbacci, M. Parrinello, and V. Schettino, Anharmonic infrared and Raman spectra in Car–Parrinello molecular dynamics simulations, J. Chem. Phys. 128, 224514 (2008).
  43. I. Batatia, D. P. Kovacs, G. Simm, C. Ortner, and G. Csanyi, MACE: Higher order equivariant message passing neural networks for fast and accurate force fields, in Advances in Neural Information Processing Systems, edited by S. Koyejo, S. Mohamed, A. Agarwal, D. Belgrave, K. Cho, and A. Oh (MIT Press, Cambridge, MA, 2022), Vol. 35, pp. 11423–11436.
  44. I. Batatia, S. Batzner, D. P. Kovács, A. Musaelian, G. N. C. Simm, R. Drautz, C. Ortner, B. Kozinsky, and G. Csányi, The design space of E(3)-equivariant atom-centred interatomic potentials, Nat. Mach. Intell. 7, 56 (2025).
  45. V. Kapil, D. P. Kovács, G. Csányi, and A. Michaelides, First-principles spectroscopy of aqueous interfaces using machine-learned electronic and quantum nuclear effects, Faraday Discuss. 249, 50 (2024).
  46. B. A. A. Martin, A. M. Ganose, V. Kapil, and K. T. Butler, General learning of the electric response of inorganic materials, arXiv:2508.17870.
  47. A. Grisafi, D. M. Wilkins, G. Csányi, and M. Ceriotti, Symmetry-adapted machine learning for tensorial properties of atomistic systems, Phys. Rev. Lett. 120, 036002 (2018).
  48. V. L. Deringer, A. P. Bartók, N. Bernstein, D. M. Wilkins, M. Ceriotti, and G. Csányi, Gaussian process regression for materials and molecules, Chem. Rev. 121, 10073 (2021).
  49. M. J. Willatt, F. Musil, and M. Ceriotti, Atom-density representations for machine learning, J. Chem. Phys. 150, 154110 (2019).
  50. A. Jana, S. Shepherd, Y. Litman, and D. M. Wilkins, Learning electronic polarizations in aqueous systems, J. Chem. Inf. Model. 64, 4426 (2024).
  51. J. Behler and M. Parrinello, Generalized neural-network representation of high-dimensional potential-energy surfaces, Phys. Rev. Lett. 98, 146401 (2007).
  52. M. Suzuki, T. Yokoyama, and M. Ito, Polarized Raman spectra of naphthalene and anthracene single crystals, Spectrochim. Acta A: Mol. Spectrosc. 24, 1091 (1968).
  53. G. N. Zhizhin and E. Mukhtarov, Vibrational spectra and structure, in Optical Spectra and Lattice Dynamics of Molecular Crystals (Elsevier Science, Amsterdam, 1995), Vol. 21, pp. 245–267.
  54. H. Seiler, M. Krynski, D. Zahn, S. Hammer, Y. W. Windsor, T. Vasileiadis, J. Pflaum, R. Ernstorfer, M. Rossi, and H. Schwoerer, Nuclear dynamics of singlet exciton fission in pentacene single crystals, Sci. Adv. 7, eabg0869 (2021).
  55. H.-H. Kowalski, Theory of thermal conductivity, Ph.D. thesis, TU Berlin, 2023.
  56. Y. Litman, J. Behler, and M. Rossi, Temperature dependence of the vibrational spectrum of porphycene: A qualitative failure of classical-nuclei molecular dynamics, Faraday Discuss. 221, 526 (2020).
  57. M. Rossi, Progress and challenges in ab initio simulations of quantum nuclei in weakly bonded systems, J. Chem. Phys. 154, 170902 (2021).
  58. J. Wang, S. Olsson, C. Wehmeyer, A. Pérez, N. E. Charron, G. de Fabritiis, F. Noé, and C. Clementi, Machine learning of coarse-grained molecular dynamics force fields, ACS Cent. Sci. 5, 755 (2019).
  59. J. Castro, G. Trenins, V. Kapil, and M. Rossi, Vibrational spectra of materials and molecules from partially-adiabatic elevated-temperature centroid molecular dynamics, J. Chem. Phys. 163, 204102 (2025).
  60. C. Kranert, C. Sturm, R. Schmidt-Grund, and M. Grundmann, Raman tensor formalism for optically anisotropic crystals, Phys. Rev. Lett. 116, 127401 (2016).
  61. J. Behler, Four generations of high-dimensional neural network potentials, Chem. Rev. 121, 10037 (2021).
  62. A. Omranpour and J. Behler, Insights into the structure and dynamics of water at Co3O4 (001) using a high-dimensional neural network potential, Chem. Cat. Chem. 18, e01321 (2026).
  63. N. Artrith, T. Morawietz, and J. Behler, High-dimensional neural-network potentials for multicomponent systems: Applications to zinc oxide, Phys. Rev. B 83, 153101 (2011).
  64. N. Artrith and J. Behler, High-dimensional neural network potentials for metal surfaces: A prototype study for copper, Phys. Rev. B 85, 045439 (2012).
  65. G. C. Sosso, G. Miceli, S. Caravati, J. Behler, and M. Bernasconi, Neural network interatomic potential for the phase change material gete, Phys. Rev. B 85, 174103 (2012).
  66. S. K. Natarajan and J. Behler, Neural network molecular dynamics simulations of solid–liquid interfaces: Water at low-index copper surfaces, Phys. Chem. Chem. Phys. 18, 28704 (2016).
  67. S. Batzner, A. Musaelian, L. Sun, M. Geiger, J. P. Mailoa, M. Kornbluth, N. Molinari, T. E. Smidt, and B. Kozinsky, E(3)-equivariant graph neural networks for data-efficient and accurate interatomic potentials, Nat. Commun. 13, 2453 (2022).
  68. R. Drautz, Atomic cluster expansion for accurate and transferable interatomic potentials, Phys. Rev. B 99, 014104 (2019).
  69. D. P. Kovács, J. H. Moore, N. J. Browning, I. Batatia, J. T. Horton, Y. Pu, V. Kapil, W. C. Witt, I.-B. Magdău, D. J. Cole, and G. Csányi, MACE-OFF: Short-range transferable machine learning force fields for organic molecules, J. Am. Chem. Soc. 147, 17598 (2025).
  70. J. P. Perdew, K. Burke, and M. Ernzerhof, Generalized gradient approximation made simple, Phys. Rev. Lett. 77, 3865 (1996).
  71. A. Tkatchenko and M. Scheffler, Accurate molecular Van der Waals interactions from ground-state electron density and free-atom reference data, Phys. Rev. Lett. 102, 073005 (2009).
  72. V. Blum, R. Gehrke, F. Hanke, P. Havu, V. Havu, X. Ren, K. Reuter, and M. Scheffler, Ab initio molecular simulations with numeric atom-centered orbitals, Comput. Phys. Commun. 180, 2175 (2009).
  73. L. M. Hunnisett, N. Francia, J. Nyman, N. S. Abraham, S. Aitipamula, T. Alkhidir, M. Almehairbi, A. Anelli, D. M. Anstine, J. E. Anthony, J. E. Arnold, F. Bahrami, M. A. Bellucci, G. J. O. Beran, R. M. Bhardwaj, R. Bianco, J. A. Bis, A. D. Boese, J. Bramley, D. E. Braun, et al., The seventh blind test of crystal structure prediction: Structure ranking methods, Acta Crystallogr. B: Struct. Sci. Cryst. Eng. Mater 80, 548 (2024).
  74. S. C. Capelli, A. Albinati, S. A. Mason, and B. T. M. Willis, Molecular motion in crystalline naphthalene: Analysis of multi-temperature X-ray and neutron diffraction data, J. Phys. Chem. A 110, 11695 (2006).
  75. H. Shang, N. Raimbault, P. Rinke, M. Scheffler, M. Rossi, and C. Carbogno, All-electron, real-space perturbation theory for homogeneous electric fields: Theory, implementation, and application within DFT, New J. Phys. 20, 073040 (2018).
  76. GitHub—ACEsuit/mace: MACE—Fast and accurate machine learning interatomic potentials with higher order equivariant message passing—github.com, https://github.com/ACEsuit/mace.
  77. GitHub—alanmlewis/TENSOAP: Public repository for symmetry-adapted Gaussian process regression (SA-GPR)—github.com, https://github.com/alanmlewis/TENSOAP.
  78. A. Togo, First-principles phonon calculations with Phonopy and Phono3py, J. Phys. Soc. Jpn. 92, 012001 (2023).
  79. A. Togo, L. Chaput, T. Tadano, and I. Tanaka, Implementation strategies in phonopy and phono3py, J. Phys.: Condens. Matter 35, 353001 (2023).
  80. Y. Litman, V. Kapil, Y. M. Y. Feldman, D. Tisi, T. Begušić, K. Fidanyan, G. Fraux, J. Higer, M. Kellner, T. E. Li, E. S. Pós, E. Stocco, G. Trenins, B. Hirshberg, M. Rossi, and M. Ceriotti, I-PI 3.0: A flexible and efficient framework for advanced atomistic simulations, J. Chem. Phys. 161, 062504 (2024).
  81. A. Togo, K. Shinohara, and I. Tanaka, Spglib: A software library for crystal symmetry search, Sci. Technol. Adv. Mater.: Methods 4, 2384822 (2024).
  82. G. Bussi, D. Donadio, and M. Parrinello, Canonical sampling through velocity rescaling, J. Chem. Phys. 126, 014101 (2007).
  83. A. Carreras, A. Togo, and I. Tanaka, DynaPhoPy: A code for extracting phonon quasiparticles from molecular dynamics simulations, Comput. Phys. Commun. 221, 221 (2017).
  84. https://doi.org/10.5281/zenodo.18445141.

Outline

Information

Sign In to Your Journals Account

Filter

Filter

Article Lookup

Enter a citation