The Chemical Space of B, N-substituted Polycyclic Aromatic Hydrocarbons: Combinatorial Enumeration and High-Throughput First-Principles Modeling

Abstract

Combinatorial introduction of heteroatoms in the two-dimensional framework of aromatic hydrocarbons opens up possibilities to design compound libraries exhibiting desirable photovoltaic and photochemical properties. Exhaustive enumeration and first-principles characterization of this chemical space provide indispensable insights for rational compound design strategies. Here, for the smallest seventy-seven Kekulean-benzenoid polycyclic systems, we reveal combinatorial substitution of C atom pairs with the isosteric and isoelectronic B, N pairs to result in 7,453,041,547,842 (7.4 tera) unique molecules. We present comprehensive frequency distributions of this chemical space, analyze trends and discuss a symmetry-controlled selectivity manifestable in synthesis product-yield. Furthermore, by performing high-throughput ab initio density functional theory calculations of over thirty-three thousand (33k) representative molecules, we discuss quantitative trends in the structural stability and inter-property relationships across heteroarenes. Our results indicate a significant fraction of the 33k molecules to be electronically active in the 1.5-2.5 eV region, encompassing the most intense region of the solar spectrum, indicating their suitability as potential light-harvesting molecular components in photo-catalyzed solar cells.

References in corpus (4)

• Boron nitride substrates for high-quality graphene electronics
• Time-dependent density functional study of the electronic spectra of oligoacenes in the charge states -1, 0, +1, and +2
• Water on BN doped benzene: A hard test for exchange-correlation functionals and the impact of exact exchange on weak binding
• A parallel algorithm for the enumeration of benzenoid hydrocarbons

Cited by in corpus (7)

• Deep learning for molecular design - a review of the state of the art
• High-Throughput Design of Peierls and Charge Density Wave Phases in Q1D Organometallic Materials
• Unlocking Inverted Singlet-Triplet Gap in Alternant Hydrocarbons with Heteroatoms
• Band gaps of long-period polytypes of IV, IV-IV, and III-V semiconductors estimated with an Ising-type additivity model
• Data-Driven Modeling of S0 -> S1 Excitation Energy in the BODIPY Chemical Space: High-Throughput Computation, Quantum Machine Learning, and Inverse Design
• Antisymmetry rules of response properties in certain chemical spaces
• Simplifying inverse material design problems for fixed lattices with alchemical chirality