Clar valence bond representation of pi-bonding in carbon nanotubes

The application of the Clar aromatic sextet valence bond (VB) model to extended, defect-free single-walled carbon nanotubes (CNTs) with roll-up vectors (m, n) provides a real space model of their electronic structure. If m - n = 3k, where k is an integer, then all pi-electrons can be represented by aromatic sextets, and the CNT is fully benzenoid; the converse is also true. Since m - n = 3k is known to be a necessary criterion for conductivity in CNTs, only fully benzenoid CNTs are metallic, and only potentially metallic CNTs are fully benzenoid. This behavior contrasts with that of planar polycyclic aromatic hydrocarbons, in which the fully benzenoid structures are known to have large HOMO-LUMO gaps. For CNTs that are not fully benzenoid, e.g., m - n = 3k + l, where l = 1 or 2 and k is an integer, a seam of double bonds wraps about an otherwise benzenoid CNT at the chiral angle - 60 degrees or the chiral angle, respectively. Nucleus-independent chemical shift calculations on hydrogen-terminated CNT segments support this, and show that the magnetic manifestation of aromatic sextets is not due to electron correlation. The resonance hybrid of the Clar VB structures corresponds to patterns occasionally observed in scanning tunneling microscopy images of CNTs.     

Optical Activity of Spin-Forbidden Electronic Transitions in Metal Complexes from Time-Dependent Density Functional Theory with Spin-Orbit Coupling

The calculation of magnetic transition dipole moments and rotatory strengths was implemented at the zeroth-order regular approximation (ZORA) two-component relativistic time-dependent density functional theory (TDDFT) level. The circular dichroism of the spin-forbidden ligand-field transitions of tris(ethylenediamine)cobalt(III) computed in this way agrees very well with available measurements. Phosphorescence dissymmetry factors and the corresponding lifetimes are evaluated for three N-heterocyclic-carbene-based iridium complexes, two of which contain helicene moieties, and for two platinahelicenes. The agreement with experimental data is satisfactory. The calculations reproduce the signs and order of magnitude of , and the large variations of phosphorescence lifetimes among the systems. The electron spin contribution to the magnetic transition dipole moment is shown to be important in all of the computations.     

Molar Entropy and Enthalpy of CO Adsorbed in Zeolites as Derived from VTIR Data: Role of Intermolecular Modes

Detailed analysis of recently reported variable-temperature IR (VTIR) spectra of carbon monoxide adsorbed in alkaline zeolites shows how, not only the corresponding values of standard adsorption enthalpy ( ) and entropy ( ) can be obtained, but also the thermodynamic values of molar entropy and enthalpy which characterize the adsorbed gas phase. In addition, it is shown that the so obtained molar entropy data can lead to new insights into soft molecular modes, which would be hardly accessible by conventional IR spectroscopic techniques.     

Cover Feature: How Different are the Diamagnetic and Paramagnetic Contributions to Off-Nucleus Shielding in Aromatic and Antiaromatic Rings? (ChemPhysChem 9/2023)

The spatial variations in the diamagnetic and paramagnetic contributions to the off-nucleus isotropic shielding, , and to the zz component of the off-nucleus shielding tensor, , around benzene (C₆H₆) and cyclobutadiene (C₄H₄) are investigated using complete-active-space self-consistent field wavefunctions. Despite the substantial differences between and around the aromatic C₆H₆ and the antiaromatic C₄H₄, the diamagnetic and paramagnetic contributions to these quantities, and , and and , are found to behave similarly in the two molecules, shielding and deshielding, respectively, each ring and its surroundings. The different signs of the most popular aromaticity criterion, the nucleus-independent chemical shift (NICS), in C₆H₆ and C₄H₄ are shown to follow from a change in the balance between the respective diamagnetic and paramagnetic contributions. Thus, the different NICS values for antiaromatic and antiaromatic molecules cannot be attributed to differences in the ease of access to excited states only; differences in the electron density, which determines the overall bonding picture, also play an important role.     

Canonical orbital contributions to the magnetic fields induced by global and local diatropic and paratropic ring currents

The induced magnetic field (IMF) of naphthalene, biphenyl, biphenylene, benzocyclobutadiene, and pentalene is dissected to contributions from the total π system, canonical π‐molecular orbitals (CMO), and HOMO→π* excitations, to evaluate and interpret relative global and local diatropicity and paratropicity. Maps of the IMF of the total π system reveal its relative strength and topology that corresponds to global and local diatropic and paratropic ring currents. The total π magnetic response is determined by this of canonical HOMOs and particularly by paratropic contributions of rotational excitations from HOMOs to unoccupied π * orbitals. Low energy excitations and similar nodal structure of HOMO and π * induce strong paratropic fields that dominate on antiaromatic rings. High energy excitations and different nodal structures lead to weak paratropic contributions of canonical HOMOs, which are overwhelmed by diatropic response of lower energy canonical orbitals in aromatic rings. CMO‐IMF analysis is found in agreement with ring current analysis. © 2017 Wiley Periodicals, Inc.     

On the Mechanical Properties of Popgraphene-Based Nanotubes: a Reactive Molecular Dynamics Study

Carbon-based tubular materials have sparked a great interest in future electronics and optoelectronics device applications. In this work, we computationally studied the mechanical properties of nanotubes generated from popgraphene (PopNTs). Popgraphene is a 2D carbon allotrope composed of 5-8-5 rings. We carried out fully atomistic reactive (ReaxFF) molecular dynamics for PopNTs of different chiralities ( and ) and/or diameters and at different temperatures (from 300 up to 1200 K). Results showed that the tubes are thermally stable (at least up to 1200 K). All tubes presented stress/strain curves with a quasi-linear behavior followed by an abrupt drop of stress values. Interestingly, armchair-like PopNTs ( ) can stand a higher strain load before fracturing when contrasted to the zigzag-like ones ( ). Moreover, it was obtained that Young's modulus (Y_Mod) (750–900 GPa) and ultimate strength (σ_US) (120–150 GPa) values are similar to the ones reported for conventional armchair and zigzag carbon nanotubes. Y_Mod values obtained for PopNTs are not significantly temperature-dependent. While the σ_US values for the showed a quasi-linear dependence with the temperature, the exhibited no clear trends.     

Tunnelling Dynamics of Kicked Systems: A Route to Tunnelling Control

The effects of discontinuously time-varying perturbations on the dynamics of a particle moving in harmonic, symmetric double well and symmetric triple well potentials, are investigated both classically and quantum mechanically. The quantum dynamics is followed using the time-dependent Fourier grid Hamiltonian (TDFGH) method while the classical dynamics is analyzed within the framework of classical Hamiltonian mechanics. Depending on the spatial symmetry of the perturbation and the characteristic features of the reversal time , different types of ‘phase space’ structures are observed in each of the potentials. For symmetric double and triple well potentials, quantum dynamics reveals that complete destruction of tunnelling (CDT) can be achieved in the presence of a time-dependent spatially asymmetric perturbing field that is continuous in time. Any discontinuity in time-variation of the perturbation may induce over the barrier transition. The relevance of these results in the context of (i) tunnelling control and (ii) quantum computing with 3-state or 2-state quantum registers is briefly discussed.     

Double-Exchange Magnetic Interactions in High-Temperature Ferromagnetic Iron Chalcogenide Monolayers

Smythite ( ) is an iron-based chalcogenide with a lamellar structure, different from the compositionally identical mineral greigite. Owing to their natural abundance, such transition metal chalcogenides are promising materials for low-cost spintronic-based devices. Herein, we discuss the charge transfer processes and complex magnetic ordering in a two-dimensional (2D) smythite lattice. We find that redox couple and complex magnetic ordering are governing factors in the charge transfer processes. A very strong ferromagnetic in-lattice coupling is also observed, which is attributed to the presence of three Fe-centres. To describe the magnetic behaviour molecular and periodic approaches have been considered. We found a substantial increase in Curie temperature with applied mechanical stress due to opening of the double exchange interaction angle. We also observe an in-plane Jahn−Teller distortion, which is further confirmed by the spin−orbit counter plot. Our study thus provides an insight into the double exchange mechanism favoured by the redox couple and results in a strong ferromagnetic ordering.     

Interpreting Aromaticity and Antiaromaticity through Bifurcation Analysis of the Induced Magnetic Field

In all molecules, a current density is induced when the molecule is subjected to an external magnetic field. In turn, this current density creates a particular magnetic field. In this work, the bifurcation value of the induced magnetic field is analyzed in a representative set of aromatic, non-aromatic and antiaromatic monocycles, as well as a set of polycyclic hydrocarbons. The results show that the bifurcation value of the ring-shaped domain adequately classifies the studied molecules according to their aromatic character. For aromatic and nonaromatic molecules, it is possible to analyze two ring-shaped domains, one diatropic (inside the molecular ring) and one paratropic (outside the molecular ring). Meanwhile, for antiaromatic rings, only a diatropic ring-shaped domain (outside the molecular ring) is possible to analyze, since the paratropic domain (inside the molecular ring) is irreducible with the maximum value (attractor) at the center of the molecular ring. In some of the studied cases, i. e., in heteroatomic species, bifurcation values do not follow aromaticity trends and present some inconsistencies in comparison to ring currents strengths, showing that this approximation provides only a qualitative estimation about (anti)aromaticity.     

The phenalenyl motif: a magnetic chameleon

The 12pi cation (3) and 14pi anion (4) derived from the phenalenyl radical (2) support diatropic ("aromatic") perimeter ring currents, but isoelectronic replacement of the central atom by either boron (5) or nitrogen (6) leads to paratropic ("antiaromatic") current; the ipsocentric approach to molecular magnetic response accounts for all four patterns in terms of competition between translationally and rotationally allowed virtual pi-pi* excitations.     

Controlling Mirror Symmetry Breaking and Network Formation in Liquid Crystalline Cubic,Isotropic Liquid and Crystalline Phases of Benzil-Based Polycatenars

Spontaneous development of chirality in systems composed of achiral molecules is important for new routes to asymmetric synthesis, chiral superstructures and materials, as well as for the understanding of the mechanisms of emergence of prebiotic chirality. Herein, it is shown that the 4,4′-diphenylbenzil unit is a universal transiently chiral bent building block for the design of multi-chained (polycatenar) rod-like molecules capable of forming a wide variety of helically twisted network structures in the liquid, the liquid crystalline (LC) and the crystalline state. Single polar substituents at the apex of tricatenar molecules support the formation of the achiral (racemic) cubic double network phase with Ia d symmetry and relatively small twist along the networks. The combination of an alkyl chain with fluorine substitution leads to the homogeneously chiral triple network phase with I23 space group, and in addition, provides a mirror symmetry broken liquid. Replacing F by Cl or Br further increases the twist, leading to a short pitch double gyroid Ia d phase, which is achiral again. The effects of the structural variations on the network structures, either leading to achiral phases or chiral conglomerates are analyzed.     

1H NMR patterns of antiaromatic systems: Paratropic displacement dependence upon LUMO-HOMO energy gap

¹H NMR chemical shifts of antiaromatic species reveal an enhanced paratropic displacement. The high-field shifts exhibited by doubly charged benzenoid polycycles were shown to be strongly related to the LUMO-HOMO energy gaps in these antiaromatic systems. As the gap decreases, a larger paratropic shift was observed.     

Visualizing electron delocalization in contorted polycyclic aromatic hydrocarbons

Electron delocalization in contorted polycyclic aromatic hydrocarbon (PAH) molecules was examined through 3D isotropic magnetic shielding (IMS) contour maps built around the molecules using pseudo-van der Waals surfaces. The resulting maps of electron delocalization provided an intuitive, yet detailed and quantitative evaluation of the aromatic, non aromatic, and antiaromatic character of the local and global conjugated cyclic circuits distributed over the molecules. An attractive pictural feature of the 3D IMS contour maps is that they are reminiscent of the Clar π-sextet model of aromaticity. The difference in delocalization patterns between the two faces of the electron circuits in contorted PAHs was clearly visualized. For π-extended contorted PAHs, some splits of the π system resulted in recognizable patterns typical of smaller PAHs. The differences between the delocalization patterns of diastereomeric chiral PAHs could also be visualized. Mapping IMS on pseudo-van der Waals surfaces around contorted PAHs allowed visualization of their superimposed preferred circuits for electron delocalization and hence their local and global aromaticity patterns.

Electron delocalization in contorted polycyclic aromatic hydrocarbon (PAH) molecules was examined through 3D isotropic magnetic shielding (IMS) contour maps built around the molecules using pseudo-van der Waals surfaces.     

Manipulating Slow Magnetic Relaxation by Light in a Charge Transfer {Fe2Co} Complex

Some cyanide-bridged complexes are known for exhibiting slow magnetic relaxation behavior in a light-induced metastable state. Herein, an unexpected reverse effect is observed for the first time in the S= {Fe^II_LS-Co^III_LS-Fe^III_LS} (HS=high spin, LS=low spin) ground state of a novel V-shaped trinuclear cyanide-bridged {Fe₂Co} complex. In this complex, light-switchable iron-cobalt charge transfer with repeatable off/on switching of slow magnetic relaxation is discovered upon alternating laser irradiation at 785 and 560 nm. An important characteristic of the present compound is that the S= ground state exhibits slow magnetic relaxation before irradiation, whereas this is accelerated after irradiation. This is different from the typical behavior, where the light-induced metastable state exhibits slow magnetic relaxation.     

Temperature Dependence of 1H Paramagnetic Chemical Shifts in Actinide Complexes,Beyond Bleaney's Theory: The AnVIO22+–Dipicolinic Acid Complexes (An=Np,Pu) as an Example

Actinide +VI complexes ( = , and ) with dipicolinic acid derivatives were synthesized and characterized by powder XRD, SQUID magnetometry and NMR spectroscopy. In addition, and complexes were described by first principles CAS based and two-component spin-restricted DFT methods. The analysis of the ¹H paramagnetic NMR chemical shifts for all protons of the ligands according to the X-rays structures shows that the Fermi contact contribution is negligible in agreement with spin density determined by unrestricted DFT. The magnetic susceptibility tensor is determined by combining SQUID, pNMR shifts and Evans’ method. The SO-RASPT2 results fit well the experimental magnetic susceptibility and pNMR chemical shifts. The role of the counterions in the solid phase is pointed out; their presence impacts the magnetic properties of the complex. The temperature dependence of the pNMR chemical shifts has a strong contribution, contrarily to Bleaney's theory for lanthanide complexes. The fitting of the temperature dependence of the pNMR chemical shifts and SQUID magnetic susceptibility by a two-Kramers-doublet model for the complex and a non-Kramers-doublet model for the complex allows for the experimental evaluation of energy gaps and magnetic moments of the paramagnetic center.     

Metal Organic Framework Cubosomes

We demonstrate a general strategy for the synthesis of ordered bicontinuous-structured metal organic frameworks (MOFs) by using polymer cubosomes (PCs) with a double primitive structure (Im m symmetry) as the template. The filling of MOF precursors in the open channel of PCs, followed by their coordination and removal of the template, generates MOF cubosomes with a single primitive topology (Pm m) and average mesopore diameters of 60–65 nm. Mechanism study reveals that the formation of ZIF-8 cubosomes undergoes a new MOF growth process, which involves the formation of individual MOF seeds in the template, their growth and eventual fusion into the cubosomes. Their growth kinetics follows the Avrami equation with an Avrami exponent of n=3 and a growth rate of k=1.33×10⁻⁴, indicating their fast 3D heterogeneous growth mode. Serving as a bioreactor, the ZIF-8 cubosomes show high loading of trypsin enzyme, leading to a high catalytic activity in the proteolysis of bovine serum albumin.     

Electronic Structure of the Low-Lying States of the Triatomic MoS2 Molecule: The Building Block of 2D MoS2

Molybdenum disulfide (MoS₂) is the building component of 1D-monolayer, 2D-layered nanosheets and nanotubes having many applications in industry, and it is detected in various molecular systems observed in nature. Here, the electronic structure and the chemical bonding of sixteen low-lying states of the triatomic MoS₂ molecule are investigated, while the connection of the chemical bonding of the isolated MoS₂ molecule to the relevant 2D-MoS₂, is emphasized. The MoS₂ molecule is studied via DFT and multireference methodologies, i. e., MRCISD(+Q)/aug-cc-pVQZ(−PP)_Mo. The ground state, ³B₁, is bent (Mo−S=2.133 Å and ϕ(SMoS)=115.9°) with a dissociation energy to atomic products of 194.7 kcal/mol at MRCISD+Q. In the ground and in the first excited state a double bond is formed between Mo and each S atom, i. e., . These two states differ in which d electrons of Mo are unpaired. The Mo−S bond distances of the calculated states range from 2.108 to 2.505 Å, the SMoS angles range from 104.1 to 180.0°, and the Mo−S bonds are single or double. Potential energy curves and surfaces have been plotted for the ³B₁, ⁵A₁ and ⁵B₁ states. Finally, the low-lying septet states of the triatomic molecule are involved in the material as a building block, explaining the variety of its morphologies.     

Predicting polycyclic aromatic hydrocarbon formation with an automatically generated mechanism for acetylene pyrolysis

Using Reaction Mechanism Generator (RMG), we have automatically constructed a detailed mechanism for acetylene pyrolysis, which predicts formation of polycyclic aromatic hydrocarbons (PAHs) up to pyrene. To improve the data available for formation pathways from naphthalene to pyrene, new high‐pressure limit reaction rate coefficients and species thermochemistry were calculated using a combination of electronic structure data from the literature and new quantum calculations. Pressure‐dependent kinetics for the CH potential energy surface calculated by Zádor et al. were incorporated to ensure accurate pathways for acetylene initiation reactions. After adding these new data into the RMG database, a pressure‐dependent mechanism was generated in a single RMG simulation which captures chemistry from C to C. In general, the RMG‐generated model accurately predicts major species profiles in comparison to plug‐flow reactor data from the literature. The primary shortcoming of the model is that formation of anthracene, phenanthrene, and pyrene are underpredicted, and PAHs beyond pyrene are not captured. Reaction path analysis was performed for the RMG model to identify key pathways. Notable conclusions include the importance of accounting for the acetone impurity in acetylene in accurately predicting formation of odd‐carbon species, the remarkably low contribution of acetylene dimerization to vinylacetylene or diacetylene, and the dominance of the hydrogen abstraction CH addition (HACA) mechanism in the formation pathways to all PAH species in the model. This work demonstrates the improved ability of RMG to model PAH formation, while highlighting the need for more kinetics data for elementary reaction pathways to larger PAHs.     

Analysis of the Binding of Cytokines to Highly Charged Polymer Networks

A model describing the binding of biological signaling proteins to highly charged polymer networks is presented. The networks are formed by polyelectrolyte chains for which the distance between two charges at the chain is smaller than the Bjerrum length. Counterion condensation on such highly charged chains immobilizes a part of the counterions. The Donnan-equilibrium between the polymer network and the aqueous solution with salt concentration $c_s^b$ is used to calculate the salt concentration of the co- and counterions $c_s^g$ entering the network. Two factors are decisive: i) The electrostatic interaction between the network and the protein is given by the Donnan-potential of the network and the net charge of the protein. In addition to this leading term, a second term describes the change in the Born-energy of the proteins when entering the network. ii) The interaction of the protein with the highly charged chains within the network is governed by counterion release: Patches of positive charge at the protein become multivalent counterions of the polyelectrolyte chains thus releasing a concomitant number of condensed counterions. The model compares favorably to experimental data obtained on a set of biohybrid polymer networks composed of crosslinked glycosaminoglycan chains that interact with a mixture of key signaling proteins.