skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on May 7, 2026

Title: Optical spectra of complex aggregates and crystals: Vibronic band structure and Davydov splitting
A thorough understanding of the nature and pattern of intermolecular interactions in molecular aggregates and crystals is a prerequisite for the design of the next generation of functional materials. In systems with multiple, symmetrically equivalent molecules per unit cell, each excited state of the isolated molecule splits into several Davydov components that appear in the absorption spectra in up to three orthogonally polarized transitions. In this work, a Frenkel–Holstein Hamiltonian is adopted to simulate the vibronic structure of the Davydov components in aggregates and crystals with up to four molecules per unit cell, where electrostatic intermolecular interactions define either 1D or 2D structures. Analysis shows that vibronic signatures report directly on the electronic couplings that contribute to the Davydov splitting and the exciton band shapes. Specifically, the vibronic signature of a given Davydov component is solely determined by its free excitonic shift. For crystals with two molecules per unit cell, the lower and upper Davydov components can each exhibit J-like or H-like behavior, resulting in JJ, JH, and HH components in order of increasing energy, all with unique vibronic signatures. Under certain conditions, null points can exist in either band, leading to a monomer-like absorption spectrum for the corresponding Davydov component. In crystals with four symmetrically equivalent molecules per unit cell, the J- or H-nature of the three orthogonally polarized Davydov components results in four possible combinations, JJJ, JJH, JHH, and HHH in order of increasing energy, all readily identified through vibronic signatures.  more » « less
Award ID(s):
2221923
PAR ID:
10588381
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
American Institute of Physics
Date Published:
Journal Name:
The Journal of Chemical Physics
Volume:
162
Issue:
17
ISSN:
0021-9606
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. In this work we study a low-cost two-dimensional ferromagnetic semiconductor with possible applications in biomedicine, solar cells, spintronics, and energy and hydrogen storage. From first principle calculations we describe the unique electronic, transport, optical, and magnetic properties of a π-conjugated micropore polymer (CMP) with three iron atoms placed in the middle of an isolated pore locally resembling heme complexes. This material exhibits strong Fe-localized d z2 bands. The bandgap is direct and equal to 0.28 eV. The valence band is doubly degenerate at the Γ -point and for larger k -wavevectors the HOMO band becomes flat with low contribution to charge mobility. The absorption coefficient is roughly isotropic. The conductivity is also isotropic with the nonzero contribution in the energy range 0.3–8 eV. The xy -component of the imaginary part of the dielectric tensor determines the magneto-optical Faraday and Kerr rotation. Nonvanishing rotation is observed in the interval of 0.5–5.0 eV. This material is found to be a ferromagnet of an Ising type with long-range exchange interactions with a very high magnetic moment per unit cell, m = 6 μ B . The exchange integral is calculated by two independent methods: (a) from the energy difference between ferromagnetic and antiferromagnetic states and (b) from a magnon dispersion curve. In the former case J nn = 27 μeV. In the latter case the magnon dispersion is fitted by the Ising model with the nearest and next-nearest neighbor spin interactions. From these estimations we find that J nn = 19.5 μeV and J nnn = −3 μeV. Despite the different nature of the calculations, the exchange integrals are only within 28% difference. 
    more » « less
  2. Abstract The excited-state properties of molecular crystals are important for applications in organic electronic devices. TheGWapproximation and Bethe-Salpeter equation (GW+BSE) is the state-of-the-art method for calculating the excited-state properties of crystalline solids with periodic boundary conditions. We present the PAH101 dataset ofGW+BSE calculations for 101 molecular crystals of polycyclic aromatic hydrocarbons (PAHs) with up to  ~500 atoms in the unit cell. To the best of our knowledge, this is the firstGW+BSE dataset for molecular crystals. The data records include theGWquasiparticle band structure, the fundamental band gap, the static dielectric constant, the first singlet exciton energy (optical gap), the first triplet exciton energy, the dielectric function, and optical absorption spectra for light polarized along the three lattice vectors. The dataset can be used to (i) discover materials with desired electronic/optical properties, (ii) identify correlations between DFT andGW+BSE quantities, and (iii) train machine learned models to help in materials discovery efforts. 
    more » « less
  3. The excited-state properties of molecular crystals are important for applications in organic electronic devices. The GW approximation and Bethe-Salpeter equation (GW+BSE) is the state-of-the-art method for calculating the excited-state properties of crystalline solids with periodic boundary conditions. We present the PAH101 dataset of GW +BSE calculations for 101 molecular crystals of polycyclic aromatic hydrocarbons (PAHs) with up to ∼500 atoms in the unit cell. The data records include the GW quasiparticle band structure, the fundamental band gap, the static dielectric constant, the first singlet exciton energy (optical gap), the first triplet exciton energy, the dielectric function, and optical absorption spectra for light polarized along the three lattice vectors. In addition, the dataset includes the density functional theory (DFT) single-molecule and crystal features used in Liu et al. [npj Computational Materials, 8, 70 (2022)]. We envision the dataset being used to (i) identify correlations between DFT and GW +BSE quantities, (ii) discover materials with desired electronic/ optical properties in the dataset itself, and (iii) train machine-learned models to help in materials discovery efforts. We provide examples to illustrate these three use cases. 
    more » « less
  4. The excited-state properties of molecular crystals are important for applications in organic electronic devices. The GW approximation and Bethe-Salpeter equation (GW +BSE) is the state-of-the-art method for calculating the excited-state properties of crystalline solids with periodic boundary conditions. We present the PAH101 dataset of GW +BSE calculations for 101 molecular crystals of polycyclic aromatic hydrocarbons (PAHs) with up to ∼500 atoms in the unit cell. The data records include the GW quasiparticle band structure, the fundamental band gap, the static dielectric constant, the first singlet exciton energy (optical gap), the first triplet exciton energy, the dielectric function, and optical absorption spectra for light polarized along the three lattice vectors. In addition, the dataset includes the density functional theory (DFT) single-molecule and crystal features used in Liu et al. [npj Computational Materials, 8, 70 (2022)]. We envision the dataset being used to (i) identify correlations between DFT and GW +BSE quantities, (ii) discover materials with desired electronic/ optical properties in the dataset itself, and (iii) train machine-learned models to help in materials discovery efforts. We provide examples to illustrate these three use cases. 
    more » « less
  5. Molecular aggregates with long-range excitonic couplings have drastically different photophysical properties compared to their monomer counterparts. From Kasha's model for one-dimensional systems, positive or negative excitonic couplings lead to blue or red-shifted optical spectra with respect to the monomers, labeled H-and J-aggregates, respectively. The overall excitonic couplings in higher dimensional systems are much more complicated and cannot be simply classified from their spectral shifts alone. Here, we provide a unified classification for extended 2D aggregates using temperature dependent peak shifts, thermal broadening, and quantum yields. We discuss the examples of six 2D aggregates with J-like absorption spectra but quite drastic changes in quantum yields and superradiance. We find the origin of the differences is, in fact, a different excitonic band structure where the bright state is lower energy than the monomer but still away from the band edge. We call this an “I-aggregate.” Our results provide a description of the complex excitonic behaviors that cannot be explained solely on Kasha's model. Furthermore, such properties can be tuned with the packing geometries within the aggregates providing supramolecular pathways for controlling them. This will allow for precise optimizations of aggregate properties in their applications across the areas of optoelectronics, photonics, excitonic energy transfer, and shortwave infrared technologies. 
    more » « less