Are Hetero-metallapentalenes Aromatic or Not? A DFT Investigation

Aromaticity is one of the most basic concepts in organic chemistry. The planar Möbius aromatic metallapentalynes and metallapentalenes have attracted considerable attention in the past few years. However, the aromaticity of metallapentalenes containing heteroatoms (such as B, N, and O), termed as hetero-metallapentalenes, is rarely studied. Herein, the stability and aromaticity of a series of hetero-metallapentalenes are theoretically investigated. The results reveal lower aromaticity in metallaborapentalene, comparable aromaticity in metallazapentalene, and nonaromaticity in metalloxapentalene relative to that of metallapentalene. Moreover, the effect of Lewis bases on the aromaticity and stability of metallaborapentalene is discussed. These results provide useful information for experimental chemists to realize more hetero-metallapentalenes.     

Sodium Carbazolide and Derivatives as Solid-State Electrolytes for Sodium-Ion Batteries

Replacing widely used organic liquid electrolytes with solid-state electrolytes (SSEs) could effectively solve the safety issues in sodium-ion batteries. Efforts on seeking novel solid-state electrolytes have been continued for decades. However, issues about SSEs still exist, such as low ionic conductivity at ambient temperature, difficulty in manufacturing, low electrochemical stability, poor compatibility with electrodes, etc. Here, sodium carbazolide (Na-CZ) and its THF-coordinated derivatives are rationally fabricated as Na⁺ conductors, and two of their crystal structures are successfully solved. Among these materials, THF-coordinated complexes exhibit fast Na⁺ conductivities, i.e., 1.20×10⁻⁴ S cm⁻¹ and 1.95×10⁻³ S cm⁻¹ at 90 °C for Na-CZ-1THF and Na-CZ-2THF, respectively, which are among the top Na⁺ conductors under the same condition. Furthermore, stable Na plating/stripping is observed even over 400 h cycling, showing outstanding interfacial stability and compatibility against Na electrode. More advantages such as ease of synthesis, low-cost, and cold pressing for molding can be obtained. In situ NMR results revealed that the evaporation of THF may play an essential role in the Na⁺ migration, where the movement of THF creates defects/vacancies and facilitates the migration of Na⁺.     

Catalytic Mechanisms of Transfer Hydrogenation of Azobenzene with Ammonia Borane by Pincer Bismuth Complex: Crucial Role of C=N Functional Group on the Pincer Ligand

Transfer hydrogenation of azobenzene with ammonia borane mediated by pincer bismuth complex 1 was systematically investigated through density functional theory calculations. An unusual metal-ligand cooperation mechanism was disclosed, in which the saturation/regeneration of the C=N functional group on the pincer ligand plays an essential role. The reaction is initiated by the hydrogenation of the C=N bond (saturation) with ammonia borane to afford 3^CN , which is the rate-determining step with Gibbs energy barrier (ΔG^≠) and Gibbs reaction energy (ΔG) of 25.6 and −7.3 kcal/mol, respectively. 3^CN is then converted to a Bi−H intermediate through a water-bridged pathway, which is followed up with the transfer hydrogenation of azobenzene to produce the final product N,N′-diphenylhydrazine and regenerate the catalyst. Finally, the catalyst could be improved by substituting the phenyl group for the tert-butyl group on the pincer ligand, where the ΔG^≠ value (rate-determining step) decreases to 24.0 kcal/mol.     

Implementation of molecular symmetry in valence bond calculation

A novel strategy for the construction of many-electron symmetry-adapted wave function is proposed for ab initio valence bond (VB) calculations and is implemented for valence bond self-consistent filed (VBSCF) and breathing orbital valence bond (BOVB) methods with various orbital optimization algorithms. Symmetry-adapted VB functions are constructed by the projection operator of symmetry group. The many-electron symmetry-adapted wave function is expressed in terms of symmetry-adapted VB functions, and thus the VB calculations can be performed with the molecular symmetry restriction. Test results show that molecular symmetry reduces the computational cost of both the iteration numbers and CPU time. Furthermore, excited states with specific symmetry can be conveniently obtained in VB calculations by using symmetry-adapted VB functions.     

基于氮杂Mn咔咯二维纳米材料的单原子Mn中心催化氧化烷烃制醇

咔咯是由四个吡咯共轭相连而形成的具有芳香性的新型卟啉类大环化合物,但咔咯分子中存在一个直接连结两个吡咯环的C?C键,与卟啉相比,仅仅是少了一个“meso”位置的C原子.因此,在结构上,咔咯含有三个“吡咯型”氮原子和一个“吡啶型”氮原子,当咔咯失去三个内氢原子后变成了三价阴离子,易与金属形成高价态的稳定配合物.氮杂咔咯是一种咔咯的meso位上的C被取代为N的咔咯衍生物.与正常的咔咯相比,它更易于与过渡金属形成稳定配合物.正是由于这些独特的结构特点,使其在金属催化、染料敏化太阳能电池、光敏剂、金属传感器、甚至在医学上都有很好的应用前景.金属有机大环均相催化剂的非均相化,是改进反应产物分离和实现催化剂循环使用的最简单有效方法之一.环境友好的Mn氮杂咔咯催化剂,在温和条件下可以利用氧气直接将有机底物氧化.本文选用Mn氮杂咔咯催化剂作为基本构建单元,通过理论计算,构建了一种新型的具有高催化活性的含Mn氮杂咔咯环结构单元的二维纳米催化材料.我们分别使用高斯软件(Gaussian 09)和维也纳从头算模拟软件包(VASP)对孤立分子和周期性体系进行结构优化以及性质的计算.在这种二维材料中,每一个Mn原子作为相对独立的金属单原子中心(SAC),保留了单环中Mn金属中心的高催化活性.在温和的光照条件下,Mn金属中心可以直接活化氧气生成类自由基[Mn]-O-O中心,随后[Mn]-O-O中心可以有效地通过夺取有机底物中的H和紧接着新生自由基的偶合反应,选择性氧化C?H键为C?OH键.另外,通过沿[Mn]-O-O反应轴施加不同强度的外电场,可对此二维纳米材料的催化反应活性进行精细调控.本文为实验上制备基于Mn氮杂咔咯的非均相催化剂以及单原子Mn基催化剂提供了理论依据.     

Predicting Dinitrogen Activation via Transition-Metal-Involved [4+2] Cycloaddition Reaction

As the strongest triple bond in nature, the N≡N triple bond activation has always been a challenging project in chemistry. On the other hand, since the award of the Nobel Prize in Chemistry in 1950, the Diels-Alder reaction has served as a powerful and widely applied tool in the synthesis of natural products and new materials. However, the application of the Diels-Alder reaction to dinitrogen activation remains less developed. Here we first demonstrate that a transition-metal-involved [4+2] Diels-Alder cycloaddition reaction could be used to activate dinitrogen without an additional reductant by density functional theory calculations. Further study reveals that such a dinitrogen activation by 1-metalla-1,3-dienes screened out from a series of transition metal complexes (38 species) according to the effects of metal center, ligand, and substituents can become favorable both thermodynamically (with an exergonicity of 28.2 kcal mol⁻¹) and kinetically (with an activation energy as low as 13.8 kcal mol⁻¹). Our findings highlight an important application of the Diels-Alder reaction in dinitrogen activation, inviting experimental chemists’ verification.     

随机薛定谔级联方程：理论和应用

随机薛定谔级联方程是一种基于波函数的严格量子动力学方法，它可用于研究耦合到玻色子热库的复杂体系中的量子动力学过程. 本综述从开放量子体系费曼路径积分的影响泛函出发，概述了随机薛定谔级联方程的一般理论框架和各种具体表述形式，并通过对复杂体系中超快激发能量转移过程的模拟来展示方法的应用范例和计算效率.     

Planar Tetracoordinate Silicon in Organic Molecules As Carbenoid-Type Amphoteric Centers: A Computational Study

Designing and synthesizing a stable compound with a planar tetracoordinate silicon (ptSi) center is a challenging goal for chemists. Here, a series of potential aromatic ptSi compounds composed of four conjugated rings shared by a centrally embedded Si atom are theoretically designed and computationally verified. Both Born–Oppenheimer molecular dynamics (BOMD) simulations and potential energy surface scannings verify the high stability and likely existence of these compounds, particularly Si-16-5555 (SiN₄C₈H₈) with 16 π electrons, under standard ambient temperature and pressure. Notably, the Hückel aromaticity rule, which works well for single rings, is inconsistent with the high stability of Si-16-5555 where the 16 p electrons are spread over four five-membered rings fused together. Bonding analyses show that the strong electron donation from the peripheral 12-membered conjugated ring with 16 π electrons to the vacant central atomic orbital Si 3p_z leads to the stabilization for both the ptSi coordination and planar aromaticity. The partial occupation of Si 3p_z results in the peculiar carbenoid-type behaviors for the amphoteric center. By modulating the electron density on the ring with substituent groups, we can regulate the nucleophilic and electrophilic properties of the central Si.     

Constructing spin-adiabatic states for the modeling of spin-crossing reactions. I. A shared-orbital implementation

In the modeling of spin-crossing reactions, it has become popular to directly explore the spin-adiabatic surfaces. Specifically, through constructing spin-adiabatic states from a two-state Hamiltonian (with spin-orbit coupling matrix elements) at each geometry, one can readily employ advanced geometry optimization algorithms to acquire a “transition state” structure, where the spin crossing occurs. In this work, we report the implementation of a fully-variational spin-adiabatic approach based on Kohn-Sham density functional theory spin states (sharing the same set of molecular orbitals) and the Breit-Pauli one-electron spin-orbit operator. For three model spin-crossing reactions (predissociation of N₂O, singlet-triplet conversion in CH₂, and CO addition to Fe(CO)₄), the spin-crossing points were obtained. Our results also indicated the Breit-Pauli one-electron spin-orbit coupling can vary significantly along the reaction pathway on the spin-adiabatic energy surface. On the other hand, due to the restriction that low-spin and high-spin states share the same set of molecular orbitals, the acquired spin-adiabatic energy surface shows a cusp (ie, a first-order discontinuity) at the crossing point, which prevents the use of standard geometry optimization algorithms to pinpoint the crossing point. An extension with this restriction removed is being developed to achieve the smoothness of spin-adiabatic surfaces.     

Parity Alternation of Silicon-doped Ternary Cationic Clusters HC_nSi_2~+(n=1～9)

Systematic study on the electronic/geometrical structures and the parity alternation effect of silicon-doped ternary cationic clusters HCnSi2＋（n = 1 ~9） have been carried out at the coupled cluster level. The ground-state （G-S） isomers of the clusters have been defined. The C, chains of the G-S isomers display polyacetylene-like structures. The even-n cations are more stable than the odd-n ones. Such a trend of even/odd alternation has been elaborated based on concepts of the bond character, atomic charge, incremental binding energy, ionization potential, proton affinity and fragmentation energies of the systems. The findings accord with the relative intensities of HC,,Si2＋ species recorded in the related mass spectrometric experiments.