运用NMR和DFT研究Michael加成反应的动力学行为及反应机理：(CO)5W=C(OEt)C2Ph作为底物

利用核磁共振方法研究了取代吡唑对炔基Fischer卡宾化合物的Michael加成的动力学行为，该反应为典型的二级反应。当吡唑的3，5-位由较大基团取代时，反应速率常数变小，而活化焓和活化熵明显增大。利用密度泛函理论研究了炔基钨卡宾为底物的Michael加成反应机理，发现吡唑上取代基团的增大可以导致第三步反应的活化能大于第一步，从而使反应的决速步骤由原来的第一步转变为第三步。     

The cycloadditions of various substituted carbenes, silylenes, and germylenes onto the diamond (100) surface: a theoretical exploration

The cycloadditions of 21 singlet substituted carbenes, silylenes, and germylenes onto the diamond (100) surface have been theoretically studied by means of density functional theory coupled with effective cluster models. The calculated reaction energies and reaction pathways have disclosed that the substituents play an important effect on the reaction profiles for the additions of carbenes, silylenes, and germylenes onto the diamond (100) surface. Our theoretical investigations illustrate that, irrespective of carbenes, silylenes, and germylenes, the cycloadditions of those with electropositive substituents (such as H and CH(3)) onto diamond (100) are much more favorable than those with electronegative and pi-donating substituents (such as F and NH(2)) both thermodynamically and kinetically. In broad perspective, we believe that a similar reactivity trend can also be extended to that of Si (100), Ge (100), fullerene, single-walled carbon nanotube, disilenes, digermenes, silenes, and germenes because all of these materials feature an analogous bonding motif.     

Reactivity and mechanism of stable heterocyclic silylenes with carbon tetrachloride

The potential energy surfaces for the reactions of stable silylenes with carbon tetrachloride have been characterized in detail using density functional theory [B3LYP/6-311G(d)], including zero-point corrections. Five stable silylene species (1-5) have been chosen in this work as model reactants. The activation barriers and enthalpies of the reactions are compared to determine the relative reactivity of the stable silylenes on the reaction potential energy surface. Our theoretical findings suggest that stable silylene 5, which has two carbon atoms bonded to the silicon center and does not contain a resonance structure, is relatively unstable with respect to the reaction with haloalkanes, in comparison with the other stable silylenes (1-4). Of the three possible reaction paths, Cl abstraction (path 1), CCl3 abstraction (path 2), and CCl4 insertion (path 3), path 1 is found to be most favorable, with a very low activation energy and a large exothermicity. In short, electronic as well as steric factors play a dominant role in determining the chemical reactivity of the stable silylene species kinetically as well as thermodynamically. Furthermore, a configuration mixing model based on the work of Pross and Shaik is used to rationalize the computational results. The results obtained allow a number of predictions to be made.     

29Si nuclear magnetic resonance of dialkylpolysilanes

The silicon-29 nuclear magnetic resonance spectra of a number of dialky- and alkylmethylpoly-silanes are reported. Polysilanes composed of asymmetrically substituted silylenes (i.e., alkylmethylsilylenes) exhibited very broad resonance lines attributed to diastereomeric chemical shifts of stereogenic silylenes alpha and beta to the observed nucles Symmetrically substituted polysilanes showed a single narrow peak. The ²⁹Si chemical shifts for these polysilanes decrease with increasing steric bulk of the substituents, varying inversely with the electronic excitation energy.     

The nuclear-spin-rotation constants of HCY, HSiY, and SiY(2) (Y=F, Cl): an ab initio study

The nuclear-spin-rotation constants of fluoro- (HCF) and chloro- (HCCl) carbene, of the corresponding silylenes (HSiF and HSiCl), and of difluoro- and dichlorosilylene (SiF(2) and SiCl(2)) are quantum-chemically investigated employing the coupled-cluster singles and doubles model augmented by a perturbative treatment of triple excitations together with various sequences of correlation-consistent basis sets. Theoretical best estimates are obtained through consideration of corrections for core correlation and of zero-point vibrational contributions. In addition, nuclear quadrupole coupling constants for the chlorine containing species are determined. A thorough comparison with experiment is made.     

The role of multifunctional kinetics during early-stage silicon hydride pyrolysis: reactivity of Si2H2 isomers with SiH4 and Si2H6

Kinetic parameters for the dominant pathways during the addition of the four Si(2)H(2) isomers, i.e., trans-HSiSiH, SiSiH(2), Si(H)SiH, and Si(H(2))Si, to monosilane, SiH(4), and disilane, Si(2)H(6), have been calculated using G3//B3LYP, statistical thermodynamics, conventional and variational transition state theory, and internal rotation corrections. The direct addition products of the multifunctional Si(2)H(2) isomers were monofunctional substituted silylenes, hydrogen-bridged species, and silenes. During addition to monosilane and disilane, the SiSiH(2) isomer was found to be most reactive over the temperature range of 800 to 1200 K. Revised parameters for the Evans-Polanyi correlation and a representative pre-exponential factor for multifunctional silicon hydride addition and elimination reaction families under pyrolysis conditions were regressed from the reactions in this study. This revised kinetic correlation was found to capture the activation energies and rate coefficients better than the current literature methods.     

Theoretical study of silylene substituent effects on the abstraction reactions with oxirane and thiirane

The potential energy surfaces for the abstraction reactions of silylenes with oxirane and thiirane have been characterized in detail using density functional theory (B3LYP) as well as the ab initio method (QCISD), including zero-point corrections. Five silylene species including SiH(2), Si(CH(3))(2), Si(NH(2))(2), Si(OH)(2), and SiF(2) have been chosen in this work as model reactants. All the interactions involve the initial formation of a donor-acceptor ylide-like complex followed by a heteroatom shift via a two-center transition state. The complexation energies, activation barriers, and enthalpies of the reactions were used comparatively to determine the relative silylenic reactivity, as well as the influence of substituents on the reaction potential energy surface. As a result, our theoretical investigations suggest that, irrespective of deoxygenation and desulfurization, the alkyl-substituted silylene abstractions are much more favorable than those of the pi donor-substituted silylenes. Moreover, for a given silylene, while both deoxygenation and desulfurization are facile processes, the deoxygenation reaction is more exothermic as well as more kinetically favorable. Furthermore, a configuration mixing model based on the work of Pross and Shaik is used to rationalize the computational results. The results obtained allow a number of predictions to be made.     

Ab initio molecular dynamics study of water at constant pressure using converged basis sets and empirical dispersion corrections

It is generally believed that studies of liquid water using the generalized gradient approximation to density functional theory require dispersion corrections in order to obtain reasonably accurate structural and dynamical properties. Here, we report on an ab initio molecular dynamics study of water in the isothermal-isobaric ensemble using a converged discrete variable representation basis set and an empirical dispersion correction due to Grimme [J. Comp. Chem. 27, 1787 (2006)]. At 300 K and an applied pressure of 1 bar, the density obtained without dispersion corrections is approximately 0.92?g/cm(3) while that obtained with dispersion corrections is 1.07 g/cm(3), indicating that the empirical dispersion correction overestimates the density by almost as much as it is underestimated without the correction for this converged basis. Radial distribution functions exhibit a loss of structure in the second solvation shell. Comparison of our results with other studies using the same empirical correction suggests the cause of the discrepancy: the Grimme dispersion correction is parameterized for use with a particular basis set; this parameterization is sensitive to this choice and, therefore, is not transferable to other basis sets.     

Electron deficient bridges involving silylenes: a theoretical study

Ab initio and density functional studies show that silylenes can form complexes with BH(3) and the resultant complexes possess 3c-2e bridges. The complexation energy for the formation of the these H-bridged structures is in the range of 18-46 kcal/mol. The characteristics of the electron deficient bridges depend on the substituents attached to the silylenes. With an increase in the pi-donating capacity of the substituents, the exothermicity of complex formation gets reduced but the kinetic stability of the H-bridged structures increase. The natural bond orbital analysis shows that all the H-bridged structures are associated with sigma(B-H)-->ppi(Si) second-order delocalization, which is responsible for the origin of the 3c-2e bonds. The complexation energies of the silylene-BH(3) complexes have been shown to have a correlation to the singlet-triplet energy gaps of silylenes.     

Electron binding capabilities of some silylenes having small singlet-triplet splittings or triplet ground states

Several silyl and alkaline metal substituted silylenes have been investigated using the CAS-ACPF method in conjunction with the aug-cc-pVTZ basis sets. Silylsilylene and disilylsilylene are found to have singlet ground states with DeltaEST(-) values of 0.676 and 0.319 eV, respectively. The adiabatic ground state electron affinities are found to be 1.572 and 2.361 eV for HSiSiH(3) and Si(SiH(3))(2). respectively. Both silylenes possesses a stable 2A1 excited negative ion state, with respective adiabatic EA values of 0.037 and 1.000 eV. In contrast, all silylenes with at least one akaline metal substituent exhibit triplet neutral ground states. The metalated silylenes HSiLi, HSiNa, LiSiSiH(3), NaSiLi, SiLi(2), and SiNa(2) have adiabatic ground state EAs somewhat below 1 eV, but each of these negatively charged system possesses up to three bound excited negative ion states, some of which are dipole-bound states.     

Exploring the oxidative cyclization of substituted N-aryl enamines: Pd-catalyzed formation of indoles from anilines

The direct Pd-catalyzed oxidative coupling of two C-H-bonds within N-aryl-enamines 1 allows the efficient formation of differently substituted indoles 2. In this cross-dehydrogenative coupling, many different functional groups are tolerated and the starting material N-aryl-enamines 1 can be easily prepared in one step from commercially available anilines. In addition, the whole sequence can also be run in a one-pot fashion. Optimization data, mechanistic insight, substrate scope, and applications are reported in this full paper.     

Reversible Silylene Insertion Reactions into Si−H and P−H σ‐Bonds at Room Temperature

Phosphine‐stabilized silylenes react with silanes and a phosphine by silylene insertion into E?H σ‐bonds (E=Si,P) at room temperature to give the corresponding silanes. Of special interest, the process occurs reversibly at room temperature. These results demonstrate that both the oxidative addition (typical reaction for transient silylenes) and the reductive elimination processes can proceed at the silicon center under mild reaction conditions. DFT calculations provide insight into the importance of the coordination of the silicon center to achieve the reductive elimination step.     

Density functional theory study of pyrophyllite and M-montmorillonites (M = Li, Na, K, Mg, and Ca): role of dispersion interactions

The stacking parameters, lattice constants, bond lengths, and bulk moduli of pyrophyllite and montmorillonites (MMTs), with alkali and alkali earth metal ions, are investigated using density functional theory with and without dispersion corrections. For pyrophyllite, it is found that the inclusion of the dispersion corrections significantly improves the agreement of the calculated values of the lattice parameters and bulk modulus with the experimental values. For the MMTs, the calculations predict that the interlayer spacing varies approximately linearly with the cation radius. The inclusion of dispersion corrections leads to sizable shifts of the interlayer spacings to shorter values. In Li-MMT, compaction of the interlayer distance triggers migration of the Li ion into the tetrahedral sheet and close coordination with basal oxygen atoms. Analysis of electron density distributions shows that the isomorphic octahedral Al(3+)/Mg(2+) substitution in MMT causes an increase of electron density on the basal oxygen atoms of the tetrahedral sheets.     

Silylenes and germylenes: The activation of H-H bond in hydrogen molecule

Possible mechanisms of activation reactions of H₂ with a variety of acyclic and cyclic silylenes and germylenes have been investigated by using the density functional theory (DFT), the second order Møller-Plesset perturbation theory (MP2), and the complete active space self-consistent field (CASSCF) method. Calculation results demonstrate the facile occurrence of the H₂ activation reaction through a concerted mechanism. The relative reactivity of H₂ splitting is closely related to the HOMO−LUMO or the singlet-triplet gaps of silylenes and germylenes. The activation energies of H₂ split by silylenes are smaller than those by germylenes. For N-heterocyclic silylenes and germylenes with the larger singlet-triplet energy gaps, the higher activation barriers are required to reach the transition states. The cyclopenta-2,4-dienylidene silylenes and germylenes are better candidates for activation reaction of H₂ with lower activation barriers. It is also shown that the halogen (F, Cl, Br) substitutions on different ring positions of the cyclopenta-2,4-dienylidene silylenes and germylenes have little influence on the activation energies and the exothermic energies of the insertion reactions with H₂.     

Self-consistent implementation of a nonlocal van der Waals density functional with a Gaussian basis set

Nearly all common density functional approximations fail to properly describe dispersion interactions responsible for binding in van der Waals complexes. Empirical corrections can fix some of the failures but cannot fully grasp the complex physics and may not be reliable for systems dissimilar to the fitting set. In contrast, the recently proposed nonlocal van der Waals density functional (vdW-DF) was derived from first principles, describes dispersion interactions in a seamless fashion, and yields the correct asymptotics. Implementation of this functional is somewhat cumbersome: Nonlocal dependence on the electron density requires numerical double integration over the space variables and functional derivatives are nontrivial. This paper shows how vdW-DF can be implemented self-consistently with Gaussian basis functions. The gradients of the energy with respect to nuclear displacements have also been derived and coded, enabling efficient geometry optimizations. We test the vdW-DF correlation functional in combination with several exchange approximations. We also study the sensitivity of the method to the basis set size and to the quality of the numerical quadrature grid. For weakly interacting systems, acceptable accuracy in semilocal exchange is achieved only with fine grids, whereas for nonlocal vdW-DF correlation even rather coarse grids are sufficient. The current version of vdW-DF is not well suited for pairing with Hartree-Fock exchange, leading to considerable overbinding.     

Re-examining the properties of the aqueous vapor-liquid interface using dispersion corrected density functional theory

First-principles molecular dynamics simulations, in which the forces are computed from electronic structure calculations, have great potential to provide unique insight into structure, dynamics, electronic properties, and chemistry of interfacial systems that is not available from empirical force fields. The majority of current first-principles simulations are driven by forces derived from density functional theory with generalized gradient approximations to the exchange-correlation energy, which do not capture dispersion interactions. We have carried out first-principles molecular dynamics simulations of air-water interfaces employing a particular generalized gradient approximation to the exchange-correlation functional (BLYP), with and without empirical dispersion corrections. We assess the utility of the dispersion corrections by comparison of a variety of structural, dynamic, and thermodynamic properties of bulk and interfacial water with experimental data, as well as other first-principles and force field-based simulations.     

Effect of dispersion corrections on covalent and non-covalent interactions in DFTB calculations

Dispersion corrections in quantum mechanical methods with the focus on non-covalent interactions have been extensively investigated in the past decade. In this paper, we elucidate the role of dispersion corrections in both non-covalent and covalent interactions within the density functional tight binding (DFTB) method. Our results suggest that two dispersion correction models, D3(BJ) and D3(CSO), generally improve different properties including barrier heights, isomerization energies, bond dissociation energies, and non-covalent binding energies. The D3(CSO) model, with fewer dispersion coefficients and DFTB-dependent parameters, was shown to perform as well as the D3(BJ) model.