Vertical and adiabatic electronic excitations in biphenylene: A theoretical study

Cyano-Polycyclic Aromatic Hydrocarbon Interstellar Candidates: Laboratory Identification, Equilibrium Structure and Astronomical Search of Cyanobiphenylene

The interplay between laboratory rotational spectroscopy and radio astronomical observations provides the most effective procedure for identifying molecules in the interstellar medium (ISM). Following the recent interstellar detections of several Polycyclic Aromatic Hydrocarbons (PAHs) and cyano derivatives in the dense molecular cloud TMC-1, it is reasonable to consider searching for other cyano-PAHs in this astronomical source. We present a rotational spectroscopy investigation of the two cyano derivatives of the PAH biphenylene, a plausible reaction product of interstellar benzyne. The rotational spectrum provided molecular parameters for the parent species and 14 monosubstituted isotopologues for each isomer. An accurate equilibrium structure was determined for both isomers using Watson's mass-dependence method (rm(2)), offering information on its uncommon ring union. Astronomical searches for the cyanobiphenylene isomers have been undertaken in TMC-1, using the QUIJOTE line survey. No lines of any isomer were found in this astronomical source, but the experimental data will serve to enable future searches for these species in the ISM.

A New Series of Sandwich-Type 5,5'-Biterphenylenes: Synthetic Challenge, Structural Uniqueness and Photodynamics

A new series of biaryls, bi-linear-terphenylenes (BLTPs), were prepared using the tert-butyllithium-mediated cyclization as the key synthetic step. The three-dimensional structures of the studied compounds were verified using X-ray crystallography and DFT calculations. Tetraaryl(ethynyl)-substituted BLTPs are highly crowded molecules, and the internal rotation around the central C-C bond is restricted due to a high barrier (>50 kcal/mol). These structures contain several aryl/terphenylenyl/aryl sandwiches, where the through-space π-π (TSPP) interactions are strongly reflected in the shielding of 1 H NMR chemical shifts, reduction of oxidation potentials, increasing aromaticity of the central six-membered ring and decreasing antiaromaticity of the four-membered rings in a terphenylenyl moiety based on NICS(0) and iso-chemical shielding surfaces. Despite the restricted C-C bond associated intramolecular TSPP interactions for BLTPs in the ground state, to our surprise, the electronic coupling between two linear terphenylenes (LTPs) in BLTPs in the excited state is weak, so that the excited-state behavior is dominated by the corresponding monomeric LTPs. In other words, all BLTPs undergo ultrafast relaxation dynamics via strong exciton-vibration coupling, acting as a blue-light absorber with essentially no emission.

Antiaromatic Molecules as Magnetic Couplers: A Computational Quest

In this study, we investigate a set of organic diradical structures in which two oxo-verdazyl radicals are selected as radical spin centers that are connected (coupled) via six coupler molecules (CM), resulting in various magnetic (ferromagnetic (FM) or antiferromagnetic (AFM)) characteristics, as reflected by their exchange coupling constants (J). We have designed 12 diradicals with 6-antiaromatic couplers coupled with bis-oxo-verdazyl diradicals with meta-meta (m-m) and para-meta (p-m) positional connectivities. The nature of the magnetic coupling (ferromagnetic, nonmagnetic, or antiferromagnetic) and the magnitude of the exchange constant J depend on the type of coupler, the connecting point between each radical center and CM, the degree of aromaticity of the coupler, and the length of the through-bond distance between radical centers. The computed magnetic exchange coupling constants J for these diradicals at the B3LYP/6-311++G(d,p) and MN12SX/6-311++G(d,p) levels of theory are large for many of these structures, indicating strong ferromagnetic coupling (with positive J values). In some cases, magnetic couplings are observed with J > 1000 cm-1 (B3LYP/6-311++G(d,p)) and strong antiferromagnetic coupling (with negative J values) with J < -1000 cm-1 (B3LYP/6-311++G(d,p)). Similarly, in some cases, magnetic couplings are observed with J > 289 cm-1 (MN12SX/6-311++G(d,p)) and strong antiferromagnetic coupling (with negative J values) with J < -568 cm-1 (MN12SX/6-311++G(d,p)). Furthermore, while numerous studies have reported that the degree of aromaticity of molecular couplers often favors strong ferromagnetic coupling, displaying the high-spin character of diradicals in their ground states, the couplers chosen in this study are characterized as antiaromatic or nonaromatic. The current investigation provides evidence that, remarkably, antiaromatic couplers are able to enhance stability by favoring electronic diradical structures with very strong ferromagnetic coupling when the length of the through-bond distance and connectivity pattern between radical centers are selected in such a way that the FM coupling is optimized. The findings in this study offer new strategies in the design of novel organic materials with interesting magnetic properties for practical applications.

Can Baird's and Clar's Rules Combined Explain Triplet State Energies of Polycyclic Conjugated Hydrocarbons with Fused 4nπ- and (4n + 2)π-Rings?

Compounds that can be labeled as "aromatic chameleons" are π-conjugated compounds that are able to adjust their π-electron distributions so as to comply with the different rules of aromaticity in different electronic states. We used quantum chemical calculations to explore how the fusion of benzene rings onto aromatic chameleonic units represented by biphenylene, dibenzocyclooctatetraene, and dibenzo[a,e]pentalene modifies the first triplet excited states (T₁) of the compounds. Decreases in T₁ energies are observed when going from isomers with linear connectivity of the fused benzene rings to those with cis- or trans-bent connectivities. The T₁ energies decreased down to those of the parent (isolated) 4nπ-electron units. Simultaneously, we observe an increased influence of triplet state aromaticity of the central 4n ring as given by Baird's rule and evidenced by geometric, magnetic, and electron density based aromaticity indices (HOMA, NICS-XY, ACID, and FLU). Because of an influence of triplet state aromaticity in the central 4nπ-electron units, the most stabilized compounds retain the triplet excitation in Baird π-quartets or octets, enabling the outer benzene rings to adapt closed-shell singlet Clar π-sextet character. Interestingly, the T₁ energies go down as the total number of aromatic cycles within a molecule in the T₁ state increases.

Electronic structure investigation of biphenylene films

Photoelectron Spectroscopy (PS) and Near-Edge X-ray Absorption Fine Structure (NEXAFS) spectroscopy have been used to investigate the occupied and empty density of states of biphenylene films of different thicknesses, deposited onto a Cu(111) crystal. The obtained results have been compared to previous gas phase spectra and single molecule Density Functional Theory (DFT) calculations to get insights into the possible modification of the molecular electronic structure in the film induced by the adsorption on a surface. Furthermore, NEXAFS measurements allowed characterizing the variation of the molecular arrangement with the film thickness and helped to clarify the substrate-molecule interaction.

The IPEA dilemma in CASPT2

Multi-configurational second order perturbation theory (CASPT2) has become a very popular method for describing excited-state properties since its development in 1990. To account for systematic errors found in the calculation of dissociation energies, an empirical correction applied to the zeroth-order Hamiltonian, called the IPEA shift, was introduced in 2004. The errors were attributed to an unbalanced description of open-shell versus closed-shell electronic states and is believed to also lead to an underestimation of excitation energies. Here we show that the use of the IPEA shift is not justified and the IPEA should not be used to calculate excited states, at least for organic chromophores. This conclusion is the result of three extensive analyses. Firstly, we survey the literature for excitation energies of organic molecules that have been calculated with the unmodified CASPT2 method. We find that the excitation energies of 356 reference values are negligibly underestimated by 0.02 eV. This value is an order of magnitude smaller than the expected error based on the calculation of dissociation energies. Secondly, we perform benchmark full configuration interaction calculations on 137 states of 13 di- and triatomic molecules and compare the results with CASPT2. Also in this case, the excited states are underestimated by only 0.05 eV. Finally, we perform CASPT2 calculations with different IPEA shift values on 309 excited states of 28 organic small and medium-sized organic chromophores. We demonstrate that the size of the IPEA correction scales with the amount of dynamical correlation energy (and thus with the size of the system), and gets immoderate already for the molecules considered here, leading to an overestimation of the excitation energies. It is also found that the IPEA correction strongly depends on the size of the basis set. The dependency on both the size of the system and of the basis set, contradicts the idea of a universal IPEA shift which is able to compensate for systematic CASPT2 errors in the calculation of excited states.

Benchmarking the performance of spin-component scaled CC2 in ground and electronically excited states

A generalization of the spin-component scaling and scaled opposite-spin modifications of second-order Møller-Plesset perturbation theory to the approximate coupled-cluster singles-and-doubles model CC2 (termed SCS-CC2 and SOS-CC2) is discussed and a preliminary implementation of ground and excited state energies and analytic gradients is reported. The computational results for bond distances, harmonic frequencies, adiabatic and 0-0 excitation energies are compared with experimental results to benchmark their performance. It is found that both variants of the spin-scaling increase the robustness of CC2 against strong correlation effects and lead for this method even to somewhat larger improvements than those observed for second-order Møller-Plesset perturbation theory. The spin-component scaling also enhances systematically the accuracy of CC2 for 0-0 excitation energies for pi --> pi* and n --> pi* transitions, if geometries are determined at the same level.