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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Magnetic response properties of carbon nano-onions.

The magnetic response of a number of double- and triple-layer carbon nano-onions (CNOs) is analyzed by calculating the magnetically induced current density and the induced magnetic field using the pseudo-π model. Qualitatively the same magnetic response was obtained in calculations at the all-electron level. The calculations show that the CNOs exhibit strong net diatropic (paratropic) ring currents when the external magnetic field points perpendicularly to one of the six-membered (five-membered) rings. They are deshielded inside and shielded outside the CNO; the latter dominates for larger CNOs. The magnetic response originates from a combination of spherical paratropic current densities on the inside of each carbon layer and diatropic current densities on the outside of them. The quantitative differences in the aromaticity of the CNOs as compared to single fullerenes are discussed in terms of ring-current strengths. The magnetic response of some of the CNOs is approximately the sum of the magnetic response of the individual layers, whereas deviations are significant for CNOs containing C80. For the largest CNOs, the deviation from the sum of the fullerene contributions is larger, especially when the external magnetic field is perpendicular to a six-membered ring.

Non-intersecting ring currents in [12]infinitene.

The aromaticity of the newly synthesized [12]infinitene is addressed via analysis of the magnetically induced current density and the induced magnetic field. Our calculations reveal that [12]infinitene responds to an external magnetic field by creating two current-density pathways that flow diatropically along the edges of the molecule. The current-density pathways do not intersect. The entire structure is completely shielded suggesting that the infinitene molecule is aromatic, contrary to what the Möbius rule for twisted annulene structures predicts. We also studied the dication of [12]infinitene, which sustains two paratropic ring currents flowing along the edges. The space between the stacked rings at the crossing point is shorter for the dication as compared to the neutral molecule. Hence, a strong through-space current density appears at the crossing point of π-π stacked rings.

On the origin of the inverted singlet-triplet gap of the 5th generation light-emitting molecules.

Excitation energies of the lowest singlet and triplet state of molecules whose first excited singlet state lies energetically below the first triplet state have been studied computationally at (time-dependent) density functional theory, coupled-cluster, and second-order multiconfiguration perturbation theory levels. The calculations at the ab initio levels show that the singlet-triplet gap is inverted as compared to the one expected from Hund's rule, whereas all density functionals yield the triplet state as the lowest excited state. Double excitations responsible for the inverted singlet-triplet gap have been identified. Employing the spin-flip and ΔSCF methods, singlet-triplet inversion was obtained at the density functional theory level for some of the studied molecules.

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.

Core-electron contributions to the molecular magnetic response.

Orbital contributions to the magnetic response depend on the method used to compute them. Here, we show that dissecting nuclear magnetic shielding tensors using natural localized molecular orbitals (NLMOs) leads to anomalous core contributions. The arbitrariness of the assignment might significantly affect the interpretation of the magnetic response of nonplanar molecules such as C60 or [14]helicene and the assessment of their aromatic character. We solve this problem by computing the core- and σ-components of the induced magnetic field (and NICS) and the magnetically induced current density by removing the valence electrons (RVE). We estimate the core contributions to the magnetic response by performing calculations on the corresponding highly charged molecules, such as C6H630+ for benzene, using gauge-including atomic orbitals and canonical molecular orbitals (CMOs). The orbital contributions to nuclear magnetic shielding tensors are usually estimated by employing a natural chemical shielding (NCS) analysis in NLMO or CMO bases. The RVE approach shows that the core contribution to the magnetic response is small and localized at the nuclei, contrary to what NCS calculations suggest. This may lead to a completely incorrect interpretation of the magnetic σ-orbital response of nonplanar structures, which may play a major role in the overall magnetic shielding of the system. The RVE approach is thus a simple and inexpensive way to determine the magnetic response of the core- and σ-electrons.