Search

Complexes with closed-shell (d10–d10) interactions have been studied for their interesting luminescence properties in organic light-emitting diode (OLED) devices.

Substitution with gold(i)-containing moieties results in non-vanishing oscillator strengths and spin–orbit coupling leading to fast intersystem crossing in light-emitting heptazine derivates with an inverted singlet–triplet gap.

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 though 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. 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 H-1 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.

The molecular structures of porphyrinoid cages were obtained by constructing small polyhedral graphs whose vertices have degree-4. The initial structures were then fully optimized at the density functional theory (DFT) level using the generalized gradient approximation. Some of polyhedral vertices were replaced with Zn-porphyrin units and their edges were replaced with ethyne or butadiyne bridges or connected by fusing two neighboring Zn-porphyrin units. Molecule 1 is an ethynebridge porphyrinoid nanotube, whose ends are sealed with a Zn-porphyrin. Molecule 2 is the corresponding open porphyrinoid nanotube. Molecule 3 is a clam-like porphyrinoid cage, whose shells consist of fused Zn-porphyrins, and the two halves are connected via butadiyne bridges. Molecule 4 is a cross-belt of fused Zn-porphyrins, and molecule 5 is a cross-belt of Zn-porphyrins connected with butadiyne bridges. The magnetically induced current density of the optimized porphyrinoid cages was calculated for determining the aromatic character, the degree of aromaticity and the current-density pathways. The current-density calculations were performed at the DFT level with the gauge-including magnetically induced currents (GIMIC) method using the B3LYP hybrid functional and def2-SVP basis sets. Calculations of the current densities show that molecule 2 sustains a paratropic ring current around the nanotube, whereas sealing the ends as in molecule 1 leads to an almost nonaromatic nanotube. Fusing porphyrinoids as in molecules 3 and 4 results in complicated current density pathways that differ from the ones usually appearing in porphyrinoids. The aromatic character of molecules 4 and 5 changes upon oxidation. The neutral molecule 4 is antiaromatic, whereas the dication is nonaromatic. Molecule 5 is nonaromatic, and its dication is aromatic.

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.

Magnetically induced ring currents are calculated from the magnetic shielding tensor by employing the Ampere-Maxwell law. The feasibility of the method is demonstrated by integrating the zz component of the shielding tensor along the symmetry axis of highly symmetric ring-shaped aromatic, antiaromatic and nonaromatic molecules. The calculated ring-current strengths agree perfectly with the ones obtained by integrating the current-density flux passing through a plane cutting half the molecular ring. The method can be used in combination with all electronic structure codes capable of calculating nuclear magnetic resonance (NMR) shielding tensors in general points in space. We also show that nucleus independent chemical shifts (NICS) along the symmetry axis are related to the spatial derivative of the strength of the global ring-current along the z axis.

[M25(MPA)18]− (M = Au, Ag; MPA = SCH2CH2COOH) clusters are modeled and their electronic and optical properties are described.

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.

Magnetically induced current densities (MICDs) of Zn-porphyrinoid nanostructures have been studied at the density functional theory level using the B3LYP functional and the def2-SVP basis set. Six of the studied Zn-porphyrinoid nanostructures consist of two crossing porphyrinoid belts, and one is a porphyrinoid nanoball belonging to the octahedral (

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

Mildly thermal air or HNO3 oxidized activated carbons catalyse oxidative dehydrogenative couplings of benzo[b]fused heteroaryl 2,2’‐dimers, e.g., 2‐(benzofuran‐2‐yl)‐1H‐indole, to chiral 3,3’‐coupled cyclooctatetraenes or carbazole‐type migrative products under O2 atmosphere. DFT calculations show that the radical cation and the Scholl‐type arenium cation mechanisms lead to different products with 2‐(benzofuran‐2‐yl)‐1H‐indole, being in accord with experimental product distributions., Oxidative coupling of benzofused heteroaryl dimers proceeded via divergent mechanistic routes depending on the used substrate and carbocatalyst. DFT level calculations using molecular quinones to compare the thermodynamics between radical and arenium cationic pathways rationalized the experimental behaviour of 2‐(benzofuran‐2‐yl)‐1H‐indole: HNO3 oxidized active carbon (AC) catalyst yielded carbazole derivative while air oxidized AC yielded COTs.

An efficient method for estimating non-adiabatic coupling matrix elements (NACME) and rate constants for internal conversion (k(IC)) is presented. The method, based on Plotnikov's theory, requires only calculations of the electronic wave functions and the corresponding electronic excitation energies. Computationally expensive calculations of the derivatives of the electronic wave function with respect to the nuclear coordinates are avoided. When the main accepting modes of the electronic excitation energy are X-H vibrations, the present method can be used for estimating the efficiency of the energy transfer between donor and acceptor molecules. It can also be used in studies of the influence of hydrogen bonding or solvent effect on fluorescence quenching, in studies of vibronic effects of TADF (thermally activated delayed fluorescence) emitters, and for calculating k(IC). Here, k(IC) and NACME are calculated for free-base porhyrin, magnesium porphyrin, azulene, naphthalene, pyrene and fluorenone interacting with a solvent molecule. Reverse k(IC) and NACME are further calculated for the T-1 -> T-2 transition of dibenzothiophene-S,S-dioxide (PTZ-DBTO2), which is used in TADF applications. Finally, we estimate the efficiency of the energy transfer between two large porphyrinoid dimers.

We give an overview of the molecular response to an external magnetic field perturbing quantum mechanical systems. We present state-of-the-art methods for calculating magnetically-induced current-density susceptibilities. We discuss the essence and properties of current-density susceptibilities and how molecular magnetic properties can be calculated from them. We also review the theory of spin-current densities, how relativity affects current densities and magnetic properties. An overview of the magnetic ring-current criterion for aromaticity is given, which has implications on theoretical and experimental research. The recently reported theory of antiaromaticity and how molecular symmetry affects the magnetic response are discussed and applied to closed-shell paramagnetic molecules. The topology of magnetically induced current densities and its consequences for molecular magnetic properties are also presented with twisted and toroidal molecules as examples.

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 C-60 or [14]helicene and the assessment of their aromatic character. We solve this problem by computing the core- and sigma-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 sigma-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 sigma-electrons.

Cyclo[

Magnetizabilities and magnetically induced current densities have been calculated and analyzed for a series of antiaromatic cyclo[4k]carbons (k = 2-11), iso[n]phlorins (n = 4-8), expanded porphyrinoids, and meso-meso, beta-beta,beta-beta triple-linked porphyrin and isophlorin arrays. The cyclo[4k]carbons with k = 2-6 are predicted to be closed-shell paramagnetic molecules due to the very strong paratropic ring current combined with its large radius. Larger cyclo[4k]carbons with k = 6-11 are diamagnetic because they sustain a paratropic ring current whose strength is weaker than -20 nA T-1, which seems to be the lower threshold value for closed-shell paramagnetism. This holds not only for cyclo[4k]carbons but also for other organic molecules like expanded porphyrinoids and oligomers of porphyrinoids. The present study shows that meso-meso, beta-beta, beta-beta triple-linked linear porphyrin and isophlorin arrays have a domainlike distribution of alternating diatropic and paratropic ring currents. The strength of their local paratropic ring currents is weaker than -20 nA T-1 in each domain. Therefore, linear porphyrin and isophlorin arrays become more diamagnetic with increasing length of the ribbon. For the same reason, square-shaped meso-meso, beta-beta, beta-beta triple-linked free-base porphyrin and isophlorin tetramers as well as the Zn(II) complex of the porphyrin tetramer are diamagnetic. We show that closed-shell molecules with large positive magnetizabilities can be designed by following the principle that a strong paratropic current ring combined with a large ring-current radius leads to closed-shell paramagnetism.

Rate constants for radiative and non-radiative transitions of the [Au(HNCOH)]3 complex and its dimer were calculated within the Herzberg–Teller approximation based on quantum mechanical principles. A high triplet quantum yield was estimated for the monomer. Internal conversion (IC) was found to be the major competing process to the intersystem crossing (ISC) from the lowest excited singlet state (S1) to the lowest triplet state (T1). ISC and IC from the spin-mixed 1 state also dominate the triplet relaxation process resulting in a negligible phosphorescence quantum yield for the monomer. The IC and ISC rate constants of the dimer are considerably smaller due to much lower Franck–Condon factors. For the dimer a triplet quantum yield of 0.71 was estimated using the extended multi-configuration quasi-degenerate second-order perturbation theory (XMCQDPT2) method to calculate the transition energies. Fluorescence is the major competing process to the ISC relaxation of the S1 state of the dimer. The ISC and IC processes are insignificant for the relaxation of the T1 state, resulting in unity phosphorescence quantum yield. The high triplet and phosphorescence quantum yields of the [Au(HNCOH)]3 dimer make it and its higher oligomers potential candidates as dopants for phosphorescent organic light emitting diodes and as down-converters in solid-state lighting systems.

A generating function method was used to simulate the vibrationally resolved absorption and emission spectra of perylene, terrylene and quaterrylene. This method operates on the basis of adiabatic excitation energies and electronic ground and excited state vibrational frequencies. These parameters were calculated using density functional theory with the PBE0 functional for perylene and terrylene and with the BH-LYP functional for quaterrylene. The vertical excitation energies of the lower excited states were calculated using functionals with differing amounts of Hartree-Fock exchange. The optimal functional for each molecule was chosen by comparing these energies to literature excitation energies. Using this technique the calculated absorption spectra and the calculated emission spectrum of perylene were found to be in excellent agreement with the literature experimental spectra after introducing a shift and a scaling factor. The most prominent bands of the absorption spectra were assigned to their respective vibronic transitions.

The magnetically induced current-density susceptibility tensor (CDT) of the lowest singlet and triplet states of the metallocenothiaporphyrins, where the metal is V, Cr, Mn, Fe, Co, Ni, Mo, Tc, Ru, or Rh, have been studied with the gauge-including magnetically induced currents (GIMIC) method. The compounds containing V, Mn, Co, Tc or Rh were studied as cations because the neutral molecules have an odd number of electrons. The calculations show that the aromatic nature of most of the studied molecules follows the Huckel and Baird rules of aromaticity. CDT calculations on the high-spin states of the neutral metallocenothiaporphyrins with V, Mn, Co, Tc or Rh also shows that these molecules follow a unified extended Huckel and Baird aromaticity orbital-count rule stating that molecules with an odd number of occupied conjugated valence orbitals are aromatic, whereas molecules with an even number of occupied conjugated orbitals are antiaromatic.

Cyclo[n]carbons (n = 5, 7, 9,..., 29) composed from an odd number of carbon atoms are studied computationally at density functional theory (DFT) and ab initio complete active space self-consistent field (CASSCF) levels of theory to get insight into their electronic structure and aromaticity. DFT calculations predict a strongly delocalized carbene structure of the cyclo[n]carbons and an aromatic character for all of them. In contrast, calculations at the CASSCF level yield geometrically bent and electronically localized carbene structures leading to an alternating double aromaticity of the odd-number cyclo[n]carbons. CASSCF calculations yield a singlet electronic ground state for the studied cyclo[n]carbons except for C25, whereas at the DFT level the energy difference between the lowest singlet and triplet states depends on the employed functional. The BHandHLYP functional predicts a triplet ground state of the larger odd-number cyclo[n]carbons starting from n = 13. Current-density calculations at the BHandHLYP level using the CASSCFoptimized molecular structures show that there is a through-space delocalization in the cyclo[n]carbons. The current density avoids the carbene carbon atom, leading to an alternating double aromaticity of the oddnumber cyclo[n]carbons satisfying the antiaromatic [4k+1] and aromatic [4k+3] rules. C11, C15, and C19 are aromatic and can be prioritized in future synthesis. We predict a bond-shift phenomenon for the triplet state of the cyclo[n]carbons leading to resonance structures that have different reactivity toward dimerization. Funding Agencies|Swedish Research Council [2020-04600]; Academy of Finland [325369, 314821]; Danish Council for Independent Research [DFF4181-00206, 9040-00265]; VILLUM FONDEN [40871]; Russian Scientific Foundation [18-19-00268-Pi]

We have assessed the accuracy of the magnetic properties of a set of 51 density functional approximations, including both recently published and already established functionals. The accuracy assessment considers a series of 27 small molecules and is based on comparing the predicted magnetizabilities to literature reference values calculated using coupled-cluster theory with full singles and doubles and perturbative triples [CCSD(T)] employing large basis sets. The most accurate magnetizabilities, defined as the smallest mean absolute error, are obtained with the BHandHLYP functional. Three of the six studied Berkeley functionals and the three range-separated Florida functionals also yield accurate magnetizabilities. Also, some older functionals like CAM-B3LYP, KT1, BHLYP (BHandH), B3LYP, and PBE0 perform rather well. In contrast, unsatisfactory performance is generally obtained with Minnesota functionals, which are therefore not recommended for calculations of magnetically induced current density susceptibilities and related magnetic properties such as magnetizabilities and nuclear magnetic shieldings. We also demonstrate that magnetizabilities can be calculated by numerical integration of magnetizability density; we have implemented this approach as a new feature in the gauge-including magnetically induced current (GIMIC) method. Magnetizabilities can be calculated from magnetically induced current density susceptibilities within this approach even when analytical approaches for magnetizabilities as the second derivative of the energy have not been implemented. The magnetizability density can also be visualized, providing additional information that is not otherwise easily accessible on the spatial origin of magnetizabilities.

The magnetically induced current density of an intriguing naphthalene-fused heteroporphyrin has been studied, using the quantum-chemical, gauge-including magnetically induced currents (GIMIC) method. The ring-current strengths and current-density pathways for the heteroporphyrin, its Pd complex, and the analogous quinoline-fused heteroporphyrin provide detailed information about their aromatic properties. The three porphyrinoids have similar current-density pathways and are almost as aromatic as free-base porphyrin. Notably, we show that the global ring current makes a branch at three specific points. Thus, the global current is composed of a total of eight pathways that include 22 pi-electrons, with no contributions from 18-electron pathways.

The π-contribution to the magnetically induced current densities, ring-current strengths, and induced magnetic fields of large planar molecules (as kekulene) and three-dimensional molecules (as [10]cyclophenacene and chiral toroidal nanotubes C2016 and C2196) have been computed using the pseudo-π model with the gauge-including magnetically induced currents method. The magnetic response analysis shows that π-electrons are the main actors of the electron delocalization in carbon systems regardless of their size, suggesting that the π- component of the ring-current strengths can be used for assessing the aromatic character of this kind of molecules. Computations using the pseudo-π model yield current densities and induced magnetic fields that are not contaminated by contributions from core and σ-electrons allowing investigations of large molecular structures as polycyclic aromatic hydrocarbons and cylindrical or toroidal carbon nanotubes.

Cyclo[18]carbon (C18) is studied computationally at density functional theory (DFT) and ab initio levels to obtain insights into its electronic structure, aromaticity, and adsorption properties on the NaCl surface. DFT functionals with a small amount of Hartree-Fock exchange as well as XMC-QDPT2 calculations fail to determine the experimentally observed polyyne molecular structure revealing a cumulene-type geometry. Exchange-correlation functionals with a large amount of Hartree-Fock ex-change as well as ab initio CASSCF calculations yield the polyyne structure is a ground state while the cumulene structure corresponds to the transition state between the two inverted polyyne structures through the Kekule distortion. The polyyne and the cumulene structures are found to be doubly Hückel aromatic. The calculated adsorption energy of cyclo[18]carbon on the NaCl surface is small (37 meV/C) and almost the same for both structures implying that the surface does not stabilize a particular geometry.

We have assessed the accuracy for magnetic properties of a set of 51 density functional approximations, including both recently published as well as already established functionals. The accuracy assessment considers a series of 27 small molecules and is based on comparing the predicted magnetizabilities to literature reference values calculated using coupled cluster theory with full singles and doubles and perturbative triples [CCSD(T)] employing large basis sets. The most accurate magnetizabilities, defined as the smallest mean absolute error, were obtained with the BHandHLYP functional. Three of the six studied Berkeley functionals and the three range-separated Florida functionals also yield accurate magnetizabilities. Also some older functionals like CAM-B3LYP, KT1, BHLYP (BHandH), B3LYP and PBE0 perform rather well. In contrast, unsatisfactory performance was generally obtained with Minnesota functionals, which are therefore not recommended for calculations of magnetically induced current density susceptibilities, and related magnetic properties such as magnetizabilities and nuclear magnetic shieldings. We also demonstrate that magnetizabilities can be calculated by numerical integration of the magnetizability density; we have implemented this approach as a new feature in the gauge-including magnetically induced current method (GIMIC). Magnetizabilities can be calculated from magnetically induced current density susceptibilities within this approach even when analytical approaches for magnetizabilities as the second derivative of the energy have not been implemented. The magnetizability density can also be visualized, providing additional information that is not otherwise easily accessible on the spatial origin of the magnetizabilities., 22 pages, 4 figures + 25 pages of supporting information

We develop a methodology for calculating, analyzing and visualizing nuclear magnetic shielding densities, which are calculated from the current density via the Biot-Savart relation. Atomic contributions to nuclear magnetic shielding constants can be estimated within our framework with a Becke partitioning scheme. The new features have been implemented in the GIMIC program and are applied in this work to the study of the $^1$H and $^{13}$C nuclear magnetic shieldings in benzene (C$_6$H$_6$) and cyclobutadiene (C$_4$H$_4$). The new methodology allows a visual inspection of the spatial origins of the positive (shielding) and negative (deshielding) contributions to the nuclear magnetic shielding constant of a single nucleus, something which has not been hitherto easily accomplished. Analysis of the shielding densities shows that diatropic and paratropic current-density fluxes yield both shielding as well as deshielding contributions, as the shielding or deshielding is determined by the direction of the current-density flux with respect to the studied nucleus instead of the tropicity. Becke partitioning of the magnetic shieldings shows that the magnetic shielding contributions mainly originate from the studied atom and its nearest neighbors, confirming the localized character of nuclear magnetic shieldings., 18 pages, 8 figures, includes supporting information

Current densities are induced in the electronic structure of molecules when they are exposed to external magnetic fields. Aromatic molecular rings sustain net diatropic ring currents, whereas the net ring current in antiaromatic molecular rings is paratropic and flows in the opposite, nonclassical direction. We present computational methods and protocols to calculate, analyze, and visualize magnetically induced current densities in molecules. Calculated ring-current strengths are used for quantifying the degree of aromaticity. The methods have been demonstrated by investigating ring-current strengths and the degree of aromaticity of aromatic, antiaromatic, and nonaromatic six-membered hydrocarbon rings. Current-density pathways and ring-current strengths of aromatic and antiaromatic porphyrinoids and other polycyclic molecules have been studied. The aromaticity and current density of Mobius-twisted molecules have been investigated to find the dependence on the twist and the spatial deformation of the molecular ring. Current densities of fullerene, gaudiene, and toroidal carbon nanotubes have also been studied.

Current densities are induced in the electronic structure of molecules when they are exposed to external magnetic fields. Aromatic molecular rings sustain net diatropic ring currents, whereas the net ring current in antiaromatic molecular rings is paratropic and flows in the opposite, non-classical direction. We present computational methods and protocols to calculate, analyse and visualise magnetically induced current densities in molecules. Calculated ring-current strengths are used for quantifying the degree of aromaticity. The methods have been demonstrated by investigating ring-current strengths and the degree of aromaticity of aromatic, antiaromatic and non-aromatic six-membered hydrocarbon rings. Current-density pathways and ring-current strengths of aromatic and antiaromatic porphyrinoids and other polycyclic molecules have been studied. The aromaticity and current density of Mobius-twisted molecules has been investigated to find the dependence on the twist and the spatial deformation of the molecular ring. Current densities of fullerene, gaudiene and toroidal carbon nanotubes have also been studied.

The recently synthesized cyclo[18]carbon molecule has been characterized in a number of studies by calculating electronic, spectroscopic, and mechanical properties. However, cyclo[18] carbon is only one member of the class of cyclo[n]carbons-standalone carbon allotrope representatives. Many of the larger members of this class of molecules have not been thoroughly investigated. In this work, we calculate the magnetically induced current density of cyclo[n]carbons in order to elucidate how electron delocalization and aromatic properties change with the size of the molecular ring (n), where n is an even number between 6 and 100. We find that the Hiickel rules for aromaticity (4k + 2) and antiaromaticity (4k) become degenerate for large C-n rings (n > 50), which can be understood as a transition from a delocalized electronic structure to a nonaromatic structure with localized current density fluxes in the triple bonds. Actually, the calculations suggest that cyclo[n]carbons with n > 50 are nonaromatic cyclic polyalkynes. The influence of the amount of nonlocal exchange and the asymptotic behavior of the exchange-correlation potential of the employed density functionals on the strength of the magnetically induced ring current and the aromatic character of the large cyclo[n]carbons is also discussed.

The effect of anion complexation on magnetically induced current densities and excitation energies of antiaromatic molecular rings has been investigated by calculations on expanded antiaromatic porphyrinoids including orangarin, rosarin, amethyrin and on a theoretically predicted strongly antiaromatic hydrocarbon ring. Magnetically induced current densities and the lowest vertical excitation energies have been calculated at the density functional theory (DFT) and time-dependent DFT (TDDFT) levels using the M06-2X functional. Similar calculations have been performed on sapphyrin, cyclo[6]carbon and rubyrin, which are aromatic expanded porphyrinoids. The calculations show that anion complexation weakens the strength of the ring currents and the degree of (anti)aromaticity of the studied porphyrinoids and the antiaromatic hydrocarbon ring, because electronic charge is transferred from the anion to the molecular ring. The anion complexation weakens the calculated ring-current strength susceptibility of the antiaromatic porphyrinoids by 5-7 nA T-1 (25-30%), by 6-16 nA T-1 (21-48%) for the aromatic porphyrinoids, and by 8 nA T-1 (27%) for the antiaromatic hydrocarbon ring, whereas the current-density pathways remain the same for most molecules. Calculations on the lowest excited states show that the electronic excitation transfers electron density from the anion to the molecular ring. An antiaromatic heterocyclic molecular ring with five inner NH moieties that was constructed from the antiaromatic hydrocarbon ring was found to have a large Cl- complexation energy of 67.6 kcal mol(-1).

A new method for calculating internal conversion rate constants (k[combining low line]IC), including anharmonic effects and using the Lagrangian multiplier technique, is proposed. The deuteration effect on k[combining low line]IC is investigated for naphthalene, anthracene, free-base porphyrin (H2P) and tetraphenylporphyrin (H2TPP). The results show that anharmonic effects are important when calculating k[combining low line]IC for transitions between electronic states that are energetically separated (ΔE) by more than 20 000-25 000 cm-1. Anharmonic effects are also important when ΔE20 000-25 000 cm-1 and when the accepting modes are X-H stretching vibrations with a frequency larger than 2000 cm-1. The calculations show that there is mixing between the S1 and S2 states of naphthalene induced by non-adiabatic interactions. The non-adiabatic interaction matrix element between the S1 and S2 states is 250 cm-1 and 50 cm-1 for the normal and fully deuterated naphthalene structure and this difference significantly affects the estimated fluorescence quantum yield. Besides aromatic hydrocarbons H2P and H2TPP, the k[combining low line]IC rate constant is also calculated for pyrometene (PM567) and tetraoxa[8]circulene (4B) with a detailed analysis of the effect of the vibrational anharmonicity.

We report the synthesis of novel perhalophenyl three- coordinated gold( I) complexes using 1,2- bis-(diphenylphosphino)benzene (dppBz) as the chelating ligand and [AuR(tht)] (R = C6F5, C6Cl2F3, C6Cl5) as the perhalophenyl-gold(I) source, leading to [AuR(dppBz)] (R = C6F5 (1), C6Cl2F3 (2), C6Cl5 (3)) complexes. The solid-state structures of compounds 2 and 3 consist of discrete three-coordinated Au(I) complexes, which show a distorted trigonal planar geometry for the gold center with dissimilar Au-P distances. The distorted structural arrangement is closely related to its photophysical properties. The studied complexes display very intense emissions at room temperature (RT) and at 77 K in the solid state. Studies of the emissive properties of the complexes at different temperatures suggest that the emissions are phosphorescent at 77 K and exhibit thermally activated delayed fluorescence (TADF) at RT. First-principle calculations of the photophysical processes yielded rate constants for intersystem crossing and reverse intersystem crossing that are in excellent agreement with experimental data.

A large series of dipeptides containing sulfur groups and antimony SbV were modeled to understand their inhibitor activity against Leishmaniasis. The trypanothione reductase (TR), which acts as a reducing agent in several vital processes, is responsible for maintaining the parasite’s cellular thiol redox balance. The antimonic SbV acid (Sb2O5·nH2O) is being evaluated as a drug with inhibitory activity against Leishmaniasis. In the present work, we investigated the inhibitory effect of antimony oxide on (TR) activity modeled as a substrate by probing two model clusters in gas phase and continuum water medium: A [(Sb2O10H8)]−2 coordinated to cysteine, and B [Sb7O28H21] coordinated to trypanothione, including glucose adduct. We report here density functional theory (DFT and DFT-D3) using (B3LYP/LANL2DZ and (TPSS/def2TZVP) results on the binding energy of cysteine and trypanthione complexed to these clusters as possible sites promoting the inhibition process. Upon viewing the results of the computational studies of cluster models and theoretical thermochemistry data for receptor-substrate interactions, identification of ligand-cluster interactions helps to unravel the mechanism of inhibition. The acidity of (Sb2O5,nH2O) leads to great cluster-dipeptide passivation. The electrostatic forces between cluster interface and dipeptide interaction present relevant inhibition effects through proton transfer or mobility from the different amine and ketone groups. The reactivity differences come from the unoccupied lone pairs 5pSb which lie at higher energy but remain available to make a good interaction with the lowest orbital p nitrogen in NH2, and in the CO groups substrate fragment in the zone HOMO-LUMO. Further cluster stability comparisons show a lower Gibbs free energy in B3 (B/trypanthione/glucose) (18 – 30 kcal/mol) at both used level in this study and gives good accurate intramolecular interactions, confirmed by the use of the dispersion-corrected density functional (DFT-D3). Given the dipeptide H-mobility and the (Sb2O5,nH2O) cluster acidity, (donor-acceptor duality), the system is predicted to be potent cluster of the inhibitors by endothermic and spontaneous reaction requiring 3.10 kcal/mol in aqueous medium.

We propose a novel class of gold-containing molecules, which have been designed using conjugated carbon structures as templates. The sp-hybridized carbons of C-2 moieties are replaced with a gold atom and one of the adjacent carbons is replaced by nitrogen. Applying the procedure to hexadehydro[12]annulene yields the well-known cyclic trinuclear gold(i) carbeniate complex. Planar, tubular and cage-shaped complexes can be obtained by taking similar sp-hybridized carbon structures as the starting point.

Rate constants for radiative and non-radiative transitions of the [Au(HN[double bond, length as m-dash]COH)]3 complex and its dimer were calculated within the Herzberg-Teller approximation based on quantum mechanical principles. A high triplet quantum yield was estimated for the monomer. Internal conversion (IC) was found to be the major competing process to the intersystem crossing (ISC) from the lowest excited singlet state (S1) to the lowest triplet state (T1). ISC and IC from the spin-mixed T[combining tilde]1 state also dominate the triplet relaxation process resulting in a negligible phosphorescence quantum yield for the monomer. The IC and ISC rate constants of the dimer are considerably smaller due to much lower Franck-Condon factors. For the dimer a triplet quantum yield of 0.71 was estimated using the extended multi-configuration quasi-degenerate second-order perturbation theory (XMCQDPT2) method to calculate the transition energies. Fluorescence is the major competing process to the ISC relaxation of the S1 state of the dimer. The ISC and IC processes are insignificant for the relaxation of the T1 state, resulting in unity phosphorescence quantum yield. The high triplet and phosphorescence quantum yields of the [Au(HN[double bond, length as m-dash]COH)]3 dimer make it and its higher oligomers potential candidates as dopants for phosphorescent organic light emitting diodes and as down-converters in solid-state lighting systems.

Magnetic response properties of 1,2,3,5-tetrazine derivatives including the newly synthesized 4,6-diphenyl-1,2,3,5-tetrazine have been studied computationally at the density functional theory (DFT) level. Calculations of magnetically induced current densities and induced magnetic fields show that the unsubstituted 1,2,3,5-tetrazine is almost as aromatic as benzene. Separating the magnetic shielding functions into molecular orbital components provided additional insights into the magnetic response. The aromatic character estimated from magnetically induced current densities and induced magnetic fields shows that NICS pi zz(0) values and ring-current strengths yield about the same degree of aromaticity, whereas NICSzz(0) and NICSzz(1) values are contaminated by sigma electron contributions. The studied 1,2,3,5-tetrazine derivatives are less aromatic than the unsubstituted one. Calculations of magnetic response properties of 4,6-diphenyl-1,2,3,5-tetrazine showed that it is the least aromatic among the studied molecules according to the ring-current criterion, while 4,6-[1,2,3,5]-ditetrazinyl-1,2,3,5-tetrazine is as aromatic as 4,6-dimethyl-1,2,3,5-tetrazine and slightly less aromatic than the unsubstituted 1,2,3,5-tetrazine.

The aromaticity of three nonplanar, fully conjugated aza-nanographenes built around a pyrrolo[3,2-b]pyrrole core is assessed through the application of two different computational procedures—GIMIC and NICS. We examine the calculated magnetically induced current densities (GIMIC) and nucleus-independent chemical shifts (NICS). The structural differences between these three apparently similar molecules lead to significantly different aromatic properties. GIMIC analysis indicates that the peripheral diatropic ring current of 3.9 nA/T for the studied bowl-shaped diaza-nanographene is the strongest, followed by the double [6]helicene which lacks seven-membered rings, and is practically nonexistent for the double [5]helicene possessing seven-membered rings. The biggest difference however is that in the two not-fully-fused molecules, the central pyrrole rings possess a significant diatropic current of about 4.1 nA/T, whereas there is no such current in the diaza-nanographene. Moreover, the antiaromaticity of the seven-membered rings is increasing while moving from double [5]helicene to diaza-nanographene (from −2.4 to −6.0 nA/T). The induced currents derived from NICSπ,zz-XY-scan analysis for all of the studied systems are in qualitative agreement with the GIMIC results. Subtle differences may originate from σ-electron currents in GIMIC or inaccuracy of NICSπ,zz values due to the nonplanarity of the systems, but the general picture is similar.

We present fully numerical electronic structure calculations on diatomic molecules exposed to an external magnetic field at the unrestricted Hartree-Fock limit, using a modified version of a recently developed finite element program, HelFEM. We have performed benchmark calculations on a few low-lying states of H2, HeH+, LiH, BeH+, BH, and CH+ as a function of the strength of an external magnetic field parallel to the molecular axis. The employed magnetic fields are in the range of $B=[0,10]~B_0$ atomic units, where $B_0 \approx 2.35 \times 10^5$ T. We have compared the results of the fully numerical calculations to ones obtained with the LONDON code using a large uncontracted gauge-including Cartesian Gaussian (GICG) basis set with exponents adopted from the Dunning aug-cc-pVTZ basis set. By comparison to the fully numerical results, we find that the basis set truncation error in the gauge-including Gaussian basis set is of the order of 1 kcal/mol at zero field, that the truncation error grows rapidly when the strength of the magnetic field increases, and that the largest basis set truncation error at $B=10~B_0$ exceeds 1000 kcal/mol. Studies in larger Gaussian basis sets suggest that reliable results can be obtained in GICG basis sets at fields stronger than $B=B_0$, provided that a sufficient coverage of higher-angular-momentum functions is included in the basis set., Comment: 21 pages and 2 figures, and a further 18 pages and 16 figures of supporting information

International audience; We present a detailed non-relativistic study of the atoms H, He, C and K and the molecule CH4 in the center of a spherical soft confinement potential of the form VN (r) = (r/r0) N with stiffness parameter N and confinement radius r0. The soft confinement potential approaches the hard-wall limit as N → ∞, giving a more detailed picture of spherical confinement. The confined hydrogen atom is considered as a base model: it is treated numerically to obtain ground-and excited-state energies and nodal positions of the eigenstates to study the convergence towards the hard-wall limit. We also derive some important analytical relations. The use of Gaussian basis sets is analyzed. We find that, for increasing stiffness parameter N , the convergence towards the basis-set limit becomes problematic. As an application, we report dipole polarizabilities for different values of N and r0 of hydrogen. For helium, we determine electron correlation effects with varying N and r0, and discuss the virial theorem for both soft and hard confinements in the limit r0 → 0. For carbon, a change in the orbital population from 2s 2 2p 2 to 2s 0 2p 4 is observed with decreasing r0, while, for potassium, we observe a change from the 2 S to 2 D ground state at small r0 values. For CH4, we show that the one-particle density becomes more spherical with increasing confinement. A possible application of soft confinement to atoms and molecules under high pressure is discussed. Dedicated to Prof. Jürgen Gauss on the occasion of his 60th birthday

Comprehensive spectroscopic kinetic studies illustrate an alternative mechanism for the traditional free-carbene intermediated H/D exchange reaction of 1,3-dialkylimidazolium salts under neutral (D2O) and acidic conditions (DCl/D2O 35 wt % solution). The deuteration of high purity [bmim]Cl in D2O is studied at different temperatures, in absence of catalyst or impurities, to yield an activation energy. DFT transition-state modelling, of a small water cluster and [bmim] cation, also yields an activation energy which strongly supports the proposed mechanism. The presence of basic impurities are shown to significantly enhance the exchange reaction, which brings into question the need for further analysis of technical purities of ionic liquids and the implications for a wide range of chemical reactions in such media.

The use of ionic liquids (Cn C1 Im)[BF4 ] with long alkyl chains (n=10, 12) in the ionothermal treatment of Na2 [HgTe2 ] led to lamellar crystal structures with molecular macrocyclic anions [Hg8 Te16 ]8- (1), the heaviest known topological relative of porphyrin. [Hg8 Te16 ]8- differs from porphyrin by the absence of an electronic π-system, which prevents a "global" aromaticity. Quantum chemical studies reveal instead small ring currents in the pyrrole-type five-membered rings that indicate weak local (σ) aromaticity. As a result of their lamellar nature, the compounds are promising candidates for the formation of sheets containing chalcogenidometalate anions.

The aromaticity of fourteen 3-oxo-verdazyl (1–8, 15) and Kuhn verdazyl (9–14) radicals with different substituents has been investigated computationally using the gauge-including magnetically induced current-density (GIMIC) method. The strength of the local and global diatropic and paratropic ring current has been obtained by performing numerical integration of the current density flow. Detailed studies of the current densities show that the aromatic character of the verdazyl ring obeys Huckel's aromaticity rule. The cations with formally 6π electrons in the verdazyl ring are weakly aromatic and the anions with 8π electrons are antiaromatic. The verdazyl ring of the neutral radicals is found to be weakly antiaromatic. The studied molecules sustain a significant global diatropic ring current that flows around the verdazyl and the phenyl rings suggesting that there is a strong coupling between the rings. For the 3-oxo-verdazyl cations, the global ring current is weaker. The spin-current density of the radicals is paratropic forming a ring current in the studied 3-oxo-verdazyls, whereas in the Kuhn verdazyls we found that the spin current density is mainly localized on the nitrogens. The calculated ionization potentials are also compared to the values deduced from cyclic voltammetry data.