03 / Research project · Summer 2025

Organic
cocrystals

First-principles calculations at the molecular and crystal scales to investigate the electronic properties of new donor–acceptor materials.

ProgramNYU Undergraduate Summer Research Program

Research groupProf. Rakesh K. Behera · NYU Tandon

FocusComputational materials science

Tools used
  • Gaussian 16
  • MedeA 3.10
  • VASP
  • NYU HPC
  • Density functional theory
Ball-and-stick view of the SS49 organic cocrystal in its crystalline unit cell
5Cocrystal candidates shown
1.1–1.5 eVDiffuse-basis molecular gaps
0.066 eVCalculated SS49 crystal gap
01

Research question

Organic cocrystals combine electron donors and acceptors in a repeating structure. Small changes in the molecules or how they stack can change their electronic properties. Our study examined newly synthesized cocrystals based on a benzimidazole-hydride (BIH) donor and TCNQ derivatives as acceptors.

I investigated their electronic structure with computational calculations on NYU's high-performance computing resources. The work used crystal structures determined by X-ray diffraction by collaborators and compared isolated donor–acceptor pairs with their periodic crystal structures.

Chemical structures of five BIH and TCNQ-based organic cocrystal candidates, labeled SS46, SS47, SS49, SS50, and SS51
The five candidate structures considered in the study. Molecular and periodic calculations answer different questions about their electronic behavior.
02

Computational approach

For isolated donor–acceptor pairs, we used Gaussian 16 with three B3LYP basis sets, checked optimized structures with frequency calculations, and collected the energies and locations of the frontier orbitals. We also calculated TTF–TCNQ under the same conditions as a reference.

At the crystal scale, we used the X-ray-derived structural files in MedeA and VASP for periodic calculations, including electronic band structures and density of states. Keeping the two scales separate matters: a molecular HOMO–LUMO gap is not the same quantity as a crystal band gap.

From measured structure to calculated properties

X-ray diffraction supplied cell parameters and atomic coordinates. The computational workflow then evaluated electronic properties; the study did not measure electrical conductivity.

03

Molecular results

With the diffuse basis sets highlighted in the final poster, calculated HOMO–LUMO gaps for the new donor–acceptor pairs were broadly consistent at 1.1–1.5 eV. The TTF–TCNQ reference was approximately 0.4 eV under the same calculation conditions.

Orbital visualizations placed the frontier orbitals primarily on the TCNQ acceptor. This suggests limited donor–acceptor coupling for the isolated molecular pairs and motivates looking at the complete crystal.

Calculated HOMO and LUMO energy levels for SS46, SS47, SS49, SS50, SS51, and reference TTF–TCNQ
Molecular energy-level comparison across three B3LYP basis sets, with TTF–TCNQ at right as a reference.
04

Crystal results

Initial periodic calculations put the band gaps of most candidates below 2 eV. SS49 had the narrowest calculated gap, 0.066 eV, making it a candidate for further investigation of charge transport. A narrow calculated gap suggests a possibility of higher conductivity; conductivity was not experimentally established here.

Molecular arrangement and unit cell of the SS49 cocrystal
SS49 crystalline unit cell.
Calculated electronic band structure of SS49 plotted against energy
Calculated SS49 band structure. The periodic result is distinct from the molecular gap above.
05

Next steps & resources

The calculations provided an initial comparison of the cocrystals and identified SS49 for closer study. Further analysis of the band structure, density of states, and charge transport was planned in the final poster. These results are computational predictions, not a measured device performance.

Three-minute UGSRP research presentation. Read the transcript.

Poster authors: Yupeng Li, Kripa Shankar SR, Swagat K. Mohapatra, and Rakesh K. Behera. The study was carried out through NYU's Undergraduate Summer Research Program with support from NYU high-performance computing resources.