自然轨道福井函数和成键活性描述符应用于解释苯硫醌和1, 3-二烯的[2+4]和[4+2]环加成反应中的成键机理

苯硫醌与脂肪族烯烃可以发生[2+4]和[4+2]环加成反应。为了解释这些环加成反应中的成键过程,本文使用了自然轨道福井函数(NOFF)与成键活性描述符。自然轨道福井函数揭示了苯硫醌和脂肪族烯烃的键或轨道的亲电性,表明电子供体的成键轨道和电子受体的反键/成键轨道之间发生了电子转移,然后成环,在这一过程中有两个共价键形成,得到了环状产物。成键活性描述符表明共价键比较容易在一个分子中具有较大f_k1⁺值的k1原子与另一个分子中具有较大f_k2^-值的k2原子之间形成。自然轨道福井函数与成键活性描述符都可以有效解释苯硫醌与1, 3-二烯之间的[2+4]与[4+2]环加成反应的机理。     

Removing electrons can increase the electron density: a computational study of negative Fukui functions

Ab initio and density-functional theory calculations for a family of substituted acetylenes show that removing electrons from these molecules causes the electron density along the C-C bond to increase. This result contradicts the predictions of simple frontier molecular orbital theory, but it is easily explained using the nucleophilic Fukui function-provided that one is willing to allow for the Fukui function to be negative. Negative Fukui functions emerge as key indicators of redox-induced electron rearrangements, where oxidation of an entire molecule (acetylene) leads to reduction of a specific region of the molecule (along the bond axis, between the carbon atoms). Remarkably, further oxidization of these substituted acetylenes (one can remove as many as four electrons!) causes the electron density along the C-C bond to increase even more. This work provides substantial evidence that the molecular Fukui function is sometimes negative and reveals that this is due to orbital relaxation.     

An example where orbital relaxation is an important contribution to the Fukui function

Density-functional electronic structure calculations are performed on the molecules Cr2(hpp)4, Mo2(hpp)4, and W2(hpp)4, where the bridging ligand, hpp, is the anion of 1,3,4,6,7,8-hexahydro-2H-pyrimido[1,2-a]pyrimidine. The calculated electronic densities are used to determine the Fukui functions. These molecules are unique not only in their ability as electron donors but also because orbital relaxation plays a decisive role in their reactivity. Unlike other examples in the literature, the reactivity of these compounds cannot be expressed solely in terms of the highest occupied and lowest unoccupied Kohn-Sham orbitals but only using the Fukui function, which includes the effects of orbital relaxation.     

Bicyclo[2.2.1]hepta-2-ene-5-yl-7-ylium radical cation: a theoretical validation of a bishomoaromatic radical cation intermediate

[structure: see text] A natural bond orbital analysis of the distonic bicyclo[2.2.1]hepta-2-ene-5-yl-7-ylium radical cation interprets its structure and radical character by a three-center two-electron bond between C2, C3, and C7 (a bishomoaromatic stabilization) and a singly occupied orbital on C5, n(5). Moreover, B3LYP/6-311+G(d,p) ESR parameters, which agree excellently with experiment, are interpreted in terms of spin polarization in the natural hybrids of sigma(C5-H5), and a dual hyperconjugative effect involving n(5), sigma(C1-H1a), sigma(C1-H1b), and antibonding counterparts.     

Inner and outer phases of antibonding orbitals of transition metal bonds: olefin metathesis

The frontier orbital theory is applied to ruthenium olefin metathesis. The formal [2+2] cycloaddition step, that is, the key step involved in the catalytic cycle of the reaction, is found to be favored by the phases of the HOMO and LUMO, in sharp contrast to [2+2] cycloaddition reactions between olefins. In the LUMO of transition metal part, a d-orbital overlaps out of phase with the vacant p-orbital of the carbene in the inner space of the metal–carbon π bond as is expected, but the remote lobe of the d orbital overlaps in phase in the outer space of the bond. This is a characteristic feature of the antibonding orbitals of transition metal bonds. The outer orbital phase plays more important role in the interaction.     

Organocatalytic, enantioselective intramolecular [6+2] cycloaddition reaction for the formation of tricyclopentanoids and insight on its mechanism from a computational study

Diphenylprolinol silyl ether was found to be an effective organocatalyst for promoting the asymmetric, catalytic, intramolecular [6 + 2] cycloaddition reactions of fulvenes substituted at the exocyclic 6-position with a δ-formylalkyl group to afford synthetically useful linear triquinane derivatives in good yields and excellent enantioselectivities. The cis-fused triquinane derivatives were obtained exclusively; the trans-fused isomers were not detected among the reaction products. The intramolecular [6 + 2] cycloaddition occurs between the fulvene functionality (6π) and the enamine double bond (2π) generated from the formyl group in the substrates and the diphenylprolinol silyl ether. The absolute configuration of the reaction products was determined by vibrational circular dichroism. The reaction mechanism was investigated using molecular orbital calculations, B3LYP and MP2 geometry optimizations, and subsequent single-point energy evaluations on model reaction sequences. These calculations revealed the following: (i) The intermolecular [6 + 2] cycloaddition of a fulvene and an enamine double bond proceeds in a stepwise mechanism via a zwitterionic intermediate. (ii) On the other hand, the intramolecular [6 + 2] cycloaddition leading to the cis-fused triquinane skeleton proceeds in a concerted mechanism via a highly asynchronous transition state. (iii) The fulvene functionality and the enamine double bond adopt the gauche-syn conformation during the C-C bond formation processes in the [6 + 2] cycloaddition. (iv) The energy profiles calculated for the intramolecular reaction explain the observed exclusive formation of the cis-fused triquinane derivatives in the [6 + 2] cycloaddition reactions. The reasons for the enantioselectivity seen in these [6 + 2] cycloaddition reactions are also discussed.     

1-乙胺基-3-甲基咪唑四氟硼酸盐吸收CO_2的理论研究

采用密度泛函理论(DFT)对离子液体1-乙胺基-3-甲基咪唑四氟硼酸盐([NH2e-mim][BF4])吸收CO_2的反应机理进行了研究.在B3LYP/6-311++G(d,p)计算水平下,对离子液体[NH2e-mim][BF4]的结构及与CO_2反应的中间体、过渡态和产物进行了全优化,获得了优化结构的构型参数、振动频率和热力学数据.利用自然键轨道(NBO)分析了离子液体[NH2e-mim][BF4]和CO_2的自然电荷布居.计算结果表明,通过阳离子[NH2e-mim]+自偶解离产生的阳离子[NH3e-mim]2+能与阴离子[BF4]-结合形成更强的离子键.根据反应吉布斯自由能变(ΔG0—)和焓变(ΔH0—)的计算结果,判断离子液体[NH2e-mim][BF4]吸收CO_2按理论摩尔比2∶1分步进行反应,吸收过程中质子的转移需克服52.51 k J/mol的能垒.     

Synthesis,electrochemistry, and photophysical properties of pyrazolino[60]fullerene–1,8-naphthalimide fluorescent derivatives

A series of pyrazolino[60]fullerene–1,8-naphthalimide (Pz[60]–NI) fluorescent derivatives were synthesized in one pot by a [3+2] dipolar cycloaddition between C₆₀ and functionalized hydrazones in good yield. In contrast with 4-aziridino[60]fullerene–1,8-naphtalimide dyads, Pz[60]–NI derivatives present stronger fluorescence intensity. Electrochemical study revealed that Pz[60]–NI presents better electron accepting character than the parent C₆₀. The natural bond orbital of the dyads were calculated using density functional theory method and found that the sp3 nitrogen atom in the pyrazoline ring plays a key role in the charge transfer process.     

Quantum computational,spectroscopic and molecular docking studies on 6-amino-3-bromo-2-methylpyridine

In order to confirm the vibrational assignments, the density functional hypothesis has been used for 6-amino-3-bromo-2-methylpyridine (6A3B2MP). The entire energy distribution is used to materialize the several vibrational modes of 6A3B2MP (TED). The DFT/B3LYP approach is used to examine the molecular optimum limits and electronic characteristics of 6A3B2MP. Investigations have been made into the molecular orbital, Fukui function analysis, natural bond orbital (NBO), and molecular electrostatic potential (MEP) characteristics of 6A3B2MP. Theoretical studies have focused on the UV–vis spectra with different solvents. According to the outcomes of the molecular docking technique, the 6A3B2MP is docked with two target proteins, which is crucial for the emergence of cancer. Therefore, this research paves the door for the development of medicines that are specially formulated.     

Atom and Bond Fukui Functions and Matrices: A Hirshfeld‐I Atoms‐in‐Molecule Approach

The Fukui function is often used in its atom‐condensed form by isolating it from the molecular Fukui function using a chosen weight function for the atom in the molecule. Recently, Fukui functions and matrices for both atoms and bonds separately were introduced for semiempirical and ab initio levels of theory using Hückel and Mulliken atoms‐in‐molecule models. In this work, a double partitioning method of the Fukui matrix is proposed within the Hirshfeld‐I atoms‐in‐molecule framework. Diagonalizing the resulting atomic and bond matrices gives eigenvalues and eigenvectors (Fukui orbitals) describing the reactivity of atoms and bonds. The Fukui function is the diagonal element of the Fukui matrix and may be resolved in atom and bond contributions. The extra information contained in the atom and bond resolution of the Fukui matrices and functions is highlighted. The effect of the choice of weight function arising from the Hirshfeld‐I approach to obtain atom‐ and bond‐condensed Fukui functions is studied. A comparison of the results with those generated by using the Mulliken atoms‐in‐molecule approach shows low correlation between the two partitioning schemes.     

Differentiating between trivalent lanthanides and actinides

The reactions of LnCl(3) with molten boric acid result in the formation of Ln[B(4)O(6)(OH)(2)Cl] (Ln = La-Nd), Ln(4)[B(18)O(25)(OH)(13)Cl(3)] (Ln = Sm, Eu), or Ln[B(6)O(9)(OH)(3)] (Ln = Y, Eu-Lu). The reactions of AnCl(3) (An = Pu, Am, Cm) with molten boric acid under the same conditions yield Pu[B(4)O(6)(OH)(2)Cl] and Pu(2)[B(13)O(19)(OH)(5)Cl(2)(H(2)O)(3)], Am[B(9)O(13)(OH)(4)]·H(2)O, or Cm(2)[B(14)O(20)(OH)(7)(H(2)O)(2)Cl]. These compounds possess three-dimensional network structures where rare earth borate layers are joined together by BO(3) and/or BO(4) groups. There is a shift from 10-coordinate Ln(3+) and An(3+) cations with capped triangular cupola geometries for the early members of both series to 9-coordinate hula-hoop geometries for the later elements. Cm(3+) is anomalous in that it contains both 9- and 10-coordinate metal ions. Despite these materials being synthesized under identical conditions, the two series do not parallel one another. Electronic structure calculations with multireference, CASSCF, and density functional theory (DFT) methods reveal the An 5f orbitals to be localized and predominately uninvolved in bonding. For the Pu(III) borates, a Pu 6p orbital is observed with delocalized electron density on basal oxygen atoms contrasting the Am(III) and Cm(III) borates, where a basal O 2p orbital delocalizes to the An 6d orbital. The electronic structure of the Ce(III) borate is similar to the Pu(III) complexes in that the Ce 4f orbital is localized and noninteracting, but the Ce 5p orbital shows no interaction with the coordinating ligands. Natural bond orbital and natural population analyses at the DFT level illustrate distinctive larger Pu 5f atomic occupancy relative to Am and Cm 5f, as well as unique involvement and occupancy of the An 6d orbitals.     

A combined theoretical and experimental study of simple terminal group 6 nitride and phosphide N[triple bond]MX3 and P[triple bond]MX3 molecules

Organometallic complexes containing terminal metal nitrides and phosphides are important synthetic reagents. Laser-ablated group 6 metal atoms react with NF 3, PF 3, and PCl 3 to form the simple lowest energy N[triple bond]MF 3, and P[triple bond]MX 3 products following insertion and halogen transfer, with the exception of P[triple bond]CrF3, which is a higher energy species and is not observed. The E[triple bond]MX3 pnictide metal trihalide molecules are identified from both argon and neon matrix infrared spectra and frequencies calculated by density functional theory and multiconfigurational second-order perturbation theory (CASSCF/CASPT2). These simple terminal nitrides involve strong triple bonds, which range from 2.80 to 2.77 to 2.59 natural bond order for M = W, Mo, and Cr, respectively, as computed by CASSCF/CASPT2, and the M[triple bond]N stretching frequencies also follow this order. The terminal phosphides are weaker with bond orders 2.74, 2.67, and 2.18, respectively, as the more diffuse 3p orbitals are less effective for bonding to the more compact metal valence d orbitals.     

Reaction mechanism and chemoselectivity of intermolecular cycloaddition reactions between phenyl-substituted cyclopropenone ketal and methyl vinyl ketone

In this paper, the mechanisms of the intermolecular [3+2] and [1+2] cycloaddition reactions of 1,1/1,3-dipolar π-delocalized singlet vinylcarbenes, which is obtained from cyclopropenone, with an electron-deficient C═O or C═C dipolarophile, to generate five-membered ring products are first disclosed by the density functional theory (DFT). Four reaction pathways, including two concerted [3+2] cycloaddition reaction pathways and two stepwise reaction pathways (an initial [1+2] cycloaddition and then a rearrangement from the [1+2] cycloadducts to the final [3+2] cycloadducts), are investigated at the B3LYP/6-31G(d,p) level of theory. The calculated results reveal that, in contrast to the concerted C═O [3+2] cycloaddition reaction pathway, which is 7.1 kcal/mol more energetically preferred compared with its stepwise reaction pathway, the C═C dipolarophile favors undergoing [1+2] cycloaddition rather than concerted [3+2] cycloaddition (difference of 5.3 kcal/mol). The lowest free energy barrier of the C═O concerted [3+2] cycloaddition reaction pathway shows that it predominates all other reaction pathways. This observation is consistent with the finding that the C═O [3 + 2] cycloadduct is the main product under experimental conditions. In addition, natural bond orbital second-order perturbation charge analyses are carried out to explain the preferred chemoselectivity of C═O to the C═C dipolarophile and the origins of cis-stereoselectivity for C═C [1+2] cycloaddition. Solvent effects are further considered at the B3LYP/6-31G(d,p) level in the solvents CH(3)CN, DMF, THF, CH(2)Cl(2), toluene, and benzene using the PCM model. The results indicate that the relative reaction trends and the main products are insensitive to the polarity of the reaction solvent.     

Theoretical investigation of some O-nitrosyl carboxylate biologic molecules �� A natural bond orbital study

Theoretical study of several O-nitrosyl carboxylate compounds have been performed using quantum computational ab initio RHF and density functional B3LYP and B3PW91 methods with 6-31G** basis set. Geometries obtained from DFT calculations were used to perform the natural bond orbital (NBO) analysis. It is noted that weakness in the O₃-N₂ bond is due to $ n_{O_1 } \to \sigma _{O_3 - N_2 }^* $ n_{O_1 } \to \sigma _{O_3 - N_2 }^* delocalization and is responsible for the longer O₃-N₂ bond lengths in O-nitrosyl carboxylate compounds. It is also noted that decreased occupancy of the localized $ \sigma _{O_3 - N_2 } $ \sigma _{O_3 - N_2 } orbital in the idealized Lewis structure, or increased occupancy of $ \sigma _{O_3 - N_2 }^* $ \sigma _{O_3 - N_2 }^* of the non-Lewis orbital, and their subsequent impact on molecular stability and geometry (bond lengths) are related with the resulting p character of the corresponding sulfur natural hybrid orbital (NHO) of $ \sigma _{O_3 - N_2 } $ \sigma _{O_3 - N_2 } bond orbital. In addition, the charge transfer energy decreases with the increase of the Hammett constants of subsitutent groups.     

The pressure tunning Raman and IR spectral studies on the multinuclear metal carbyne complexes

The Raman and infrared (IR) spectra of four tungsten metal carbyne complexes I, II, IV and V [Cl(CO)2(L)W[triple bond]CC6H4[triple bond](C[triple bond]CC6H4)n[triple bond]N[triple bond]C[triple bond]]2M (L = TMEDA, n = 0, M = PdI2 or ReCl(CO)3; L = DPPE, n = 1, M = PdI2 or ReCl(CO)3) were studied at high external pressure. Their pressure-induced phase transitions were observed near 20kbar (complexes I), 15 kbar (complexes II), 25 kbar (complex IV) and 30 kbar (complex V). The pressure-induced phase transition likely is first order in complex I and the pressure-induced phase transitions of complexes II, IV and V are mostly second order. The pressure sensitivities d nu/dp of nu(W[triple bond]C) are high in the low-pressure phase area and very low in the high-pressure phase area due to the pressure strengthening pi back-bonding from metal W to pi* orbital of C[triple bond]O in fragment Cl(CO)2(L)W[triple bond]C. The pressure strengthening metal pi back-bonding from metal Re or Pd to pi* orbital of C[triple bond]O or C[triple bond]N also happened to both of central metal centers of NCPd(I2)CN in complex I and NCReCl(CO)3CN in complex II.     

Synthesis of unsymmetrical substituted 1,4-dihydropyridines through thermal and microwave assisted [4+2] cycloadditions of 1-azadienes and allenic esters

Thermal and microwave assisted [4+2] cycloadditions of 1,4-diaryl-1-aza-1,3-butadienes with allenic esters lead to cycloadducts, which after a 1,3-H shift afford variedly substituted unsymmetrical 2-alkyl-1,4-diaryl-3-ethoxycarbonyl-1,4-dihydropyridines in high yields. Reactions carried out under microwave irradiation are cleaner and give higher yields with much shortened reaction times. Density functional theory (DFT) at the B3LYP/6-31G* level has been used to calculate geometric features of the reactants, barrier for s-trans to s-cis and reverse isomerization of azadienes (5a-d, 10a-e), dihedral angles between N(1), C(2), C(3), and C(4) atoms of azadienes along with various indices such as chemical hardness (eta), chemical potential (micro), global electrophilicity (omega), and the difference in global electrophilicity (Deltaomega) between the reacting pairs and Fukui functions (f (+) and f(-)). The results revealed that s-trans is the predominant conformation of azadienes at ambient temperature and the barrier for conversion of the s-trans rotamer of 1-azadienes to s-cis may be the major factor influencing the chemoselectivity, i.e., [4+2] verses [2+2] cycloaddition. The regiochemistry of the observed cycloadditions is collated with the obtained local electrophilicity indices (Fukui functions). Transition states for the formation of both [4+2] and [2+2] cycloadducts as located at the PM3 level indicate that the transition state for the formation of [4+2] cycloadducts has lower energy, again supporting the earlier conclusion that preferred formation of [4+2] cycloaaducts at higher temperature may be a consequence of barrier for s-trans to s-cis transformation of 1-azadienes.     

Natural bond orbital analysis of some S‐nitrosothiols biological molecules

Theoretical study of several S‐nitrosothiols biological molecules has been performed using quantum computational ab initio RHF and density functional B3LYP and B3PW91 methods with 6‐31G(d,p) basis set. Geometries obtained from DFT calculations were used to perform natural bond orbital (NBO) analysis. It is noted that the weakness in the S? N sigma bond is due to n_O1→σ delocalization and is responsible for the longer S? N bond length in S‐nitrosothiols. It is also noted that decreased occupancy of the localized σ_SN orbital in the idealized Lewis structure or increased occupancy of σ of the non‐Lewis orbital, and their subsequent impact on molecular stability and geometry (bond lengths) are related with the resulting p character of the corresponding sulfur natural hybrid orbital of σ_SN bond orbital. In addition, the charge transfer energy decreases with the increasing of the Hammett constants of substituent groups, and the partial charge distribution on the skeletal atoms shows that the electrostatic repulsion or attraction between atoms can give a significant contribution to the intramolecular and intermolecular interaction. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

The Fukui matrix: a simple approach to the analysis of the Fukui function and its positive character

The Fukui matrix is introduced as the derivative of the one-electron reduced density matrix with respect to a change in the number of electrons under constant external potential. The Fukui matrix extends the Fukui function concept: the diagonal of the Fukui matrix is the Fukui function. Diagonalizing the Fukui matrix gives a set of eigenvectors, the Fukui orbitals, and accompanying eigenvalues. At the level of theory used, there is always one dominant eigenvector, with an eigenvalue equal to 1. The remaining eigenvalues are either zero or come in pairs with eigenvalues of the same magnitude but opposite sign. Analysis of the frontier molecular orbital coefficient in the eigenvector with eigenvalue 1 gives information on the quality of the frontier molecular orbital picture. The occurrence of negative Fukui functions can be easily interpreted in terms of the nodal character of the dominant eigenvector versus the characteristics of the remaining eigenvectors and eigenvalues.