Changing aromatic properties through stacking: the face-to-face dimer of Ni(II) bis(pentafluorophenyl)norcorrole.

Nuclear magnetic resonance (NMR) shielding constants have been calculated for Ni(II) bis(pentafluorophenyl)norcorrole and its face-to-face stacked dimer at the Hartree-Fock (HF), second-order Møller-Plesset perturbation theory (MP2), complete-active-space self-consistent-field (CASSCF) levels as well as at density functional theory (DFT) levels using several functionals. The calculated 1H NMR shielding constants agree rather well with the experimental ones. The shielding constants of N and Ni calculated at DFT, HF, and MP2 levels differ from those obtained in the CASSCF calculations due to near-degeneracy effects at the Ni atom. The calculated magnetically induced current densities show that the monomer is antiaromatic, sustaining a strong global paratropic ring current, and the dimer is aromatic, sustaining a strong diatropic ring current. Qualitatively the same current density is obtained at the employed levels of theory. The most accurate ring-current strengths are probably obtained at the MP2 level. The aromatic dimer has a short intermolecular distance of less than 3 Å. The intermolecular interaction changes the nature of the frontier orbitals leading to a formal double bond between the norcorrole macrocycles.

Molecular phosphorescence enhancement by the plasmon field of metal nanoparticles.

A theoretical model is proposed that allows the estimation of the quantum yield of phosphorescence of dye molecules in the vicinity of plasmonic nanoparticles. For this purpose, the rate constants of the radiative and nonradiative intramolecular transitions for rhodamine 123 (Rh123) and brominated rhodamine (Rh123-2Br) dyes have been calculated. The plasmon effect of Ag nanoparticles on various types of luminescence processes has been studied both theoretically and experimentally. We show that in the presence of a plasmonic nanoparticle, the efficiency of the immediate and delayed fluorescence increases significantly. The phosphorescence rate of the rhodamine dyes also increases near plasmonic nanoparticles. The long-lived luminescence i.e., delayed fluorescence and phosphorescence is more enhanced for Rh123-2Br than for Rh123. The largest phosphorescence quantum yield is obtained when the dye molecule is at a distance of 4-6 nm from the nanoparticle surface. Our results can be used in the design of plasmon-enhancing nanostructures for light-emitting media, organic light-emitting diodes, photovoltaic devices, and catalysts for activation of molecular oxygen.

Organometallic copper(II) complex of meso-meso N-methyl N-confused pyrrole-bridged doubly N-methyl N-confused hexaphyrin.

Synthesis, spectroscopic and theoretical characterization of a hitherto unknown meso-meso N-confused N-methylpyrrole-bridged doubly N-confused hexaphyrin (molecule 5) and its organometallic copper(II) complex (molecule 6) are reported herein. The absence of Q-type bands in the UV-Vis spectrum and the high chemical shifts of the inner proton signals of 5 suggest its globally non-aromaticity. The spectroscopic evidence of non-aromaticity for 5 and the paramagnetic nature of 6, are fully supported by density functional theory (DFT) calculations of the UV-Vis spectra, electron paramagnetic resonance (EPR) g-tensor parameters, and the magnetically induced current density strengths obtained with the gauge-including magnetically induced currents (GIMIC) method.

Close Stacking of Antiaromatic Ni(II) Norcorrole Originating from a Four-Electron Multicentered Bonding Interaction.

A π-conjugated molecule with one electronic spin often forms a π-stacked dimer through molecular orbital interactions between two unpaired electrons. The bonding is recognized as a multicentered two-electron interaction between the two π-conjugated molecules. Here, we disclose a multicentered bonding interaction between two antiaromatic molecules involving four electrons. We have synthesized an antiaromatic porphyrin analogue, Ni(II) bis(pentafluorophenyl)norcorrole. Its dimer adopts a face-to-face stacked structure with an extremely short stacking distance of 2.97 Å. The close stacking originates from a multicenter four-electron bonding interaction between the two molecules. The bonding electrons were experimentally observed via synchrotron X-ray diffraction analysis and corroborated by theoretical calculations. The intermolecular interaction of the molecular orbitals imparts the stacked dimer with aromatic character that is distinctly different from that of its monomer.

Experimental and computational studies of the optical properties of 2,5,8-tris(phenylthiolato)heptazine with an inverted singlet-triplet gap.

Photophysical properties of the three-fold symmetric 2,5,8-tris(phenylthiolato)heptazine molecule (1) are studied from combined experimental and computational viewpoints. The intense blue photoemission of 1 in the solid state and in toluene solution is proposed to have a fluorescent origin on the basis of a relatively short emission lifetime and no detectable triplet decay. Calculations at correlated ab initio levels of theory also show that 1 has a large inverted singlet-triplet (IST) gap, a non-vanishing spin-orbit coupling matrix element between the first excited singlet and triplet states, and a fast intersystem crossing rate constant that leads to singlet population from the higher-lying triplet state. The IST gap implies that the first excited singlet state is the lowest excited one, agreeing with the measured fluorescent behaviour of 1. IST gaps are also obtained for the oxygen-containing (2) and selenium-containing (3) analogues of 1 at the ADC(2) level of theory, but not for the tellurium one (4). Calculations of the magnetically induced current density demonstrate that the heptazine core of 1 is globally non-aromatic due to the alternation of carbon and nitrogen atoms along its external rim.

Understanding Aromaticity in [5]Helicene-Bridged Cyclophanes: A Comprehensive Study.

This study explores the aromaticity of doubly [5]helicene-bridged (1,4)cyclophane and triply [5]helicene-bridged (1,3,5)cyclophane via calculations of the magnetic response and of electronic aromaticity indices. The primary objective is to assess the π-electron delocalization to determine whether they sustain global ring currents associated with π aromaticity. The molecules show local ring currents in the presence of an external magnetic field. The ring currents flow diatropically in the stacked six-membered rings and in the helicene arms. However, these π currents are not interconnected due to the discontinuity of the π delocalization at the C-C single bonds connecting the central six-membered rings to the helicene arms. Electronic indices suggest that the helicene-arm systems have significantly smaller electron delocalization than benzene. The reduction in the delocalization does not compromise their ability to exhibit ring currents in the presence of an external magnetic field. The analysis provides further evidence that the magnetic criteria yield a different degree of aromaticity for the helicene arms than obtained in the calculation of the electronic aromaticity indices. However, both approaches confirm that the studied molecules are not globally aromatic.

Metallaantiaromaticity of 10-Platinacorrole Complexes.

The aromaticity of cyclic π-conjugated organometallic compounds is known as metallaaromaticity. π-Conjugated metallacycles, such as metallabenzenes and metallapentalenes, have been investigated in order to understand the involvement of the d electrons from the metal center in the π-conjugated systems of the organic ligands. Here, we report the synthesis of Pd(II) 10-platinacorrole complexes with cyclooctadiene (COD) and norbornadiene (NBD) ligands. While the Pd(II) 10-platinacorrole COD complex adopts a distorted structure without showing appreciable antiaromaticity, the corresponding NBD complex exhibits a distinct antiaromatic character due to its highly planar conformation. Detailed density functional theory (DFT) calculations revealed that two d orbitals are involved in macrocyclic π-conjugation. We furthermore demonstrated that Craig-Möbius antiaromaticity is not present in the studied system. The synthesis of 10-platinacorroles enables a systematic comparison of the antiaromaticity and aromaticity of closely related porphyrin analogues, providing a better understanding of π-conjugation that involves d orbitals.

Revisiting Gauge-Independent Kinetic Energy Densities in Meta-GGAs and Local Hybrid Calculations of Magnetizabilities.

In a recent study [J. Chem. Theory Comput. 2021, 17, 1457-1468], some of us examined the accuracy of magnetizabilities calculated with density functionals representing the local density approximation (LDA), generalized gradient approximation (GGA), meta-GGA (mGGA), as well as global hybrid (GH) and range-separated (RS) hybrid functionals by assessment against accurate reference values obtained with coupled-cluster theory with singles, doubles, and perturbative triples [CCSD(T)]. Our study was later extended to local hybrid (LH) functionals by Holzer et al. [J. Chem. Theory Comput. 2021, 17, 2928-2947]; in this work, we examine a larger selection of LH functionals, also including range-separated LH (RSLH) functionals and strong-correlation LH (scLH) functionals. Holzer et al. also studied the importance of the physically correct handling of the magnetic gauge dependence of the kinetic energy density (τ) in mGGA calculations by comparing the Maximoff-Scuseria formulation of τ used in our aforementioned study to the more physical current-density extension derived by Dobson. In this work, we also revisit this comparison with a larger selection of mGGA functionals. We find that the newly tested LH, RSLH, and scLH functionals outperform all of the functionals considered in the previous studies. The various LH functionals afford the seven lowest mean absolute errors while also showing remarkably small standard deviations and mean errors. Most strikingly, the best two functionals are scLHs that also perform remarkably well in cases with significant multiconfigurational character, such as the ozone molecule, which is traditionally excluded from statistical error evaluations due to its large errors with common density functionals.

Current-Density Calculations on Zn-Porphyrin40 Nanorings.

Two porphyrinoid nanorings have been studied computationally. They were built by linking 40 Zn-porphyrin units with butadiyne bridges. The molecular structures belonging to the D40h point group were fully optimized with the Turbomole program at the density functional theory (DFT) level using the B3LYP functional and the def2-SVP basis sets. The aromatic character was studied at the DFT level by calculating the magnetically induced current-density (MICD) susceptibility using the GIMIC program. The neutral molecules are globally non-aromatic with aromatic Zn-porphyrin units. Charged nanorings could not be studied because almost degenerate frontier orbitals led to vanishing optical gaps for the cations. Since DFT calculations of the MICD are computationally expensive, we also calculated the MICD using three pseudo-π models. Appropriate pseudo-π models were constructed by removing the outer hydrogen atoms and replacing all carbon and nitrogen atoms with hydrogen atoms. The central Zn atom was either replaced with a beryllium atom or with two inner hydrogen atoms. Calculations with the computationally inexpensive pseudo-π models yielded qualitatively the same magnetic response as obtained in the all-electron calculations.

Influence of plasmons on the luminescence properties of solvatochromic merocyanine dyes with different solvatochromism.

The effect of localized surface plasmon resonance (LSPR) of a system consisting of a highly dipolar merocyanine dye and a silver nanoparticle (NP) was studied experimentally and theoretically. A theoretical model for estimating the fluorescence quantum yield (φfl) using quantum chemical calculations of intramolecular and intermolecular electronic transition rate constants was developed. Calculations show that the main deactivation channels of the lowest excited singlet state of the studied merocyanines are internal conversion (kIC(S1 → S0)) and fluorescence (kr(S1 → S0)). The intersystem-crossing transition has a low probability due to the large energy difference between the singlet and triplet levels. In the presence of plasmonic NPs, the fluorescence quantum yield is increased by a factor of two according to both experiment and computations. The calculated values of φfl, when considering changes in kr(S1 → S0) and the energy-transfer rate constant (ktransfer) from the dye to the NP was also twice as large at distances of 6-8 nm between the NP and the dye molecule. We also found that the LSPR effect can be increased or decreased depending on the value of the dielectric constant (εm) of the environment.

The magnetically induced current density of the [12]infinitene dianion.

The magnetically induced current-density susceptibility (MICD) of the [12]infinitene dianion and the induced magnetic field around it have been calculated at the density functional theory level. Separation of the MICD into diatropic and paratropic contributions shows that the MICD is dominated by the diatropic contribution in contrast to the notion that it is antiaromatic, which was reported in a recently published article. The MICD of the [12]infinitene dianion exhibits several through-space MICD pathways, whereas it sustains only weak local paratropic current-density contributions. We identified four main current-density pathways of which two are similar to the ones for neutral [12]infinitene. It is difficult to judge from calculations of the nucleus independent shielding constants and the induced magnetic field around the [12]infinitene dianion whether it sustains diatropic or paratropic ring currents.

Aromaticity: Quo Vadis.

Aromaticity is one of the most deeply rooted concepts in chemistry. But why, if two-thirds of existing compounds can be classified as aromatic, is there no consensus on what aromaticity is? σ-, π-, δ-, spherical, Möbius, or all-metal aromaticity… why are so many attributes needed to specify a property? Is aromaticity a dubious concept? This perspective aims to reflect where the aromaticity community is and where it is going.

On the antiaromatic-aromatic-antiaromatic transition of the stacked cyclobutadiene dimer.

We have studied the changes in the aromatic nature of two cyclobutadiene (C4H4) molecules on decreasing the intermolecular distance and approaching the cubane structure in a face-to-face fashion. The analysis based on the calculations of the magnetically induced current density and the induced magnetic field shows that the aromaticity of the two C4H4 rings changes from a strongly antiaromatic character at long distances to an aromatic transition state of stacked C4H4 rings at intermediate internuclear distances when approaching the antiaromatic state of cubane.

Gold(i)-containing light-emitting molecules with an inverted singlet-triplet gap.

Delayed fluorescence from molecules with an inverted singlet-triplet gap (DFIST) is the consequence of the unusual reverse order of the lowest excited singlet (S1) and triplet (T1) states of thermally activated delayed fluorescence (TADF) emitters. Heptazine (1,3,4,6,7,9,9b-heptaazaphenalene) derivatives have an inverted singlet-triplet gap thanks to the combination of multiple resonance (MR) effects and a significant double excitation character. Here, we study computationally the effect of gold(i) metalation and coordination on the optical properties of heptazine (molecule 4) and the phosphine-functionalized 2,5,8-tris(dimethylphosphino)heptazine derivatives (molecules 1-3). Ab initio calculations at the approximate second-order coupled cluster (CC2) and extended multiconfigurational quasi degenerate perturbation theory at the second order (XMC-QDPT2) levels show that molecules 1-4 have an inverted singlet-triplet gap due to the alternating spatial localization of the electron and hole of the exciton in the heptazine core. A non-vanishing one-electron spin-orbit coupling operator matrix element between T1 and and a fast S1 ← T1 intersystem crossing rate constant (k ISC) calculated at the XMC-QDPT2(12,12) level of theory for molecule 4 suggest that this new family of complexes may be the first organometallic DFIST emitters reported.

Current-density pathways in figure-eight-shaped octaphyrins.

We have calculated the current density induced by an external magnetic field in a set of figure-eight-shaped expanded porphyrinoids. The studied octaphyrins can be divided into three classes (N2, N4, and N6) based on the number of the inner hydrogen atoms of the pyrrole rings. Using the Runge-Kutta method, the current density is split into diatropic and paratropic contributions that are analyzed separately. The calculations show that one common ring current consists of two rather independent pathways. Each of them follows the outer side of the molecular frame of one half of the molecule and passes to the inner side of the frame on the other half. The ring-current pathways are similar to the ones for [12]infinitene. However, the current density of the octaphyrins is more complex having many branching points and pathways. Vertical through-space current-density pathways pass in the middle of the molecules through a plane that is parallel to the figure-eight-shaped view of the molecules when the magnetic field is perpendicular to the plane. The isolectronic N2 and the N4 dication sustain a weak paratropic ring current inside the molecule, which is also observed in the 1H NMR magnetic shielding constant of the inner hydrogen atoms. The diatropic current-density contribution dominates in the studied molecules. For the N4 and N6 molecules, the global current-density pathways are only diatropic and N6 sustains the strongest global diatropic current-density flux of 13.2 nA T-1.

A "gold standard" computational proof for the existence of gold(III) aurophilicity.

The existence of aurophilic gold(III)⋯gold(III) interactions has for a long time been neglected due to structural arguments and comparison with the aurophilicity of gold(I) compounds. We show with calculations at the CCSD(T) level of theory that the [AuIII(CH3)3(NH3)]2 dimer has a metallophilic dispersion interaction between the gold(III) atoms of 10.5 kJ mol-1. The aurophilic interaction is illustrated by topological QTAIM calculations and IRI analysis.

A theoretical study of M-M' polar-covalent bonding in heterobimetallic multinuclear organometallic complexes of monovalent group 11 metal centres.

Complexes with closed-shell (d10-d10) interactions have been studied for their interesting luminescence properties in organic light-emitting diode (OLED) devices. The present computational study aims at understanding the chemical bonding/interactions in a series of molecules with unusually short metal-metal bond distances between monovalent coinage-metal (d10-d10) centres. The investigated molecules include pentanuclear complexes with M or M' = Cu(I), Ag(I), or Au(I) and Mes = 2,4,6-Me3C6H2. In such complexes, the M-M' distances are up to 50-100 pm shorter than typical metallophilic bonds in homometallic analogues. Characterization and analysis of the chemical bond strength was performed using ab initio methods, density functional theory methods including a semi-empirical treatment of dispersion interactions (DFT-D3) and semi-empirical calculations at the extended Hückel theory (EHT) level. Population analysis suggests that hybridization occurs by mixing the (n + 1)s and (n + 1)p orbitals of M with the (nd) orbitals of M'. The orbital mixing plays a pivotal role in the polydentated polar-covalency/dative M-M' bonds that distinguish this bonding from the weaker metallophilic interactions.

Gold(I)···Lanthanide(III) Bonds in Discrete Heterobimetallic Compounds: A Combined Computational and Topological Study.

The chemical nature of the ligand-unsupported gold(I)-lanthanide(III) bond in the proposed [LnIII(η5-Cp)2][AuIPh2] (Ln-Au; LnIII = LaIII, EuIII, or LuIII; Cp = cyclopentadienide; Ph = phenyl) models is examined from a theoretical viewpoint. The covalent bond-like Au-Ln distances (Au-La, 2.95 Å; Au-Eu, 2.85 Å; Au-Lu, 2.78 Å) result from a strong interaction between the oppositely charged fragments (ΔEintMP2 > 600 kJ mol-1), including the aforementioned metal-metal bond and additional LnIII-Cipso and C-H···π interactions. The Au-Ln bond has been characterized as a chemical bond rather than a strong metallophilic interaction with the aid of energy decomposition analysis, interaction region indicator, and quantum theory of atoms in molecules topological tools. The chemical nature of the Au-Ln bond cannot be fully ascribed to a covalent or an ionic model; an intermediate situation or a charge shift bond is proposed. The [AuIPh2]- anion has also been identified as a suitable lanthanide(III) emission sensitizer for La-Au and Lu-Au.

B3Al4+: A Three-Dimensional Molecular Reuleaux Triangle.

We systematically explore the potential energy surface of the B3Al4+ combination of atoms. The putative global minimum corresponds to a structure formed by an Al4 square facing a B3 triangle. Interestingly, the dynamical behavior can be described as a Reuleaux molecular triangle since it involves the rotation of the B3 triangle at the top of the Al4 square. The molecular dynamics simulations, corroborating with the very small rotational barriers of the B3 triangle, show its nearly free rotation on the Al4 ring, confirming the fluxional character of the cluster. Moreover, while the chemical bonding analysis suggests that the multicenter interaction between the two fragments determines its fluxionality, the magnetic response analysis reveals this cluster as a true and fully three-dimensional aromatic system.

Time-dependent density functional theory studies of the optical and electronic properties of the [M25(MPA)18]- (M = Au, Ag, MPA = SCH2CH2COOH) clusters.

Electronic absorption and emission spectra of [M25(MPA)18]- (M = Au, Ag; MPA = SCH2CH2COOH) clusters have been recorded for energies below 5 eV at the time-dependent density functional theory (TDDFT) level using B3LYP and TPSSh functionals and compared to the calculated ones using the computationally inexpensive simplified TDDFT (sTDDFT) approach. The results show a qualitative agreement between the TDDFT and sTDDFT approaches used here, which were also in line with the experimental and theoretical spectra previously reported. However, the sTDDFT calculations were several orders of magnitude faster than those obtained by TDDFT. Our results support that sTDDFT provides an excellent balance between accuracy and low computational cost, becoming a suitable approach for studying the absorption and emission spectra of noble-metal clusters of sizes that would be unaffordable by the traditional TDDFT methods. The main peaks of the experimental absorption spectrum of [Au25(MPA)18]- have been previously assigned, whereas [Ag25(MPA)18]- has not been synthesised. However, its absorption spectrum resulted in having similar features to the experimental spectrum of [Ag25(GSH)18]- (GSH = glutathione), used to validate our results. The emission spectra, which to date have not been reported either from experimental or theoretical means, were simulated by using the molecular structure of the first excited triplet state (T1). The emission spectra were obtained by comparing them to those of [M25(GSH)18]- since no experimental luminescence spectra have been reported for [M25(MPA)18]-. The calculations suggest that [Ag25(SR)18]- (SR = thiolate) clusters have a weak luminescence band in the NIR region. Finally, solvent shifts were found to be minor, whereas the absorption bands seem to be significantly redshifted in solid-state materials.