The Dowd-Beckwith ring expansion: a theoretical study

UB3LYP/6-311++G(d,p) and ROMP2/6-311++G(d,p)//UB3LYP/6-311++G(d,p) calculations including the effect of benzene solvent through the PCM-UAHF method render a concerted mechanism without fragmentation as the most favourable one for the Dowd-Beckwith radical ring expansion of the bromomethyl adduct of methyl cyclopentanone-2-carboxylate to yield methyl cyclohexanone-3-carboxylate. The corresponding concerted TS is a bicyclic alcoxy radical.     

NMR spectra, GIAO and charge density calculations of five-membered aromatic heterocycles

The B3LYP/6-31+G(d) molecular geometry optimized structures of 17 five-membered heterocycles were employed together with the gauge including atomic orbitals (GIAO) density functional theory (DFT) method at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p) and B3LYP/6-311+G(2d,p) levels of theory for the calculation of proton and carbon chemicals shifts and coupling constants. The method of geometry optimization for pyrrole (1), N-methylpyrrole (2) and thiophene (7) using the larger 6-311++G(d,p) basis sets at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p), B3LYP/6-31+G(2d,p) and B3LYP/cc-pVTZ levels of theory gave little difference between calculated and experimental values of coupling constants. In general, the (1)H and 13C chemical shifts for all compounds are in good agreement with theoretical calculations using the smaller 6-31 basis set. The values of nJHH(n=3, 4, 5) and rmnJ(CH)(n=1, 2, 3, 4) were predicted well using the larger 6-31+G(d,p) and 6-311++G(d,p) basis sets and at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p), B3LYP/6-31+G(2d,2p) levels of theory. The computed atomic charges [Mülliken; Natural Bond Orbital Analysis (NBO); Merz-Kollman (MK); CHELP and CHELPG] for the B3LYP/6-311++G(d,p) geometry optimized structures of 1-17 were used to explore correlations with the experimental proton and carbon chemical shifts.     

On the flexibility of beta-peptides

The full conformational space was explored for an achiral and two chiral beta-peptide models: namely For-beta-Ala-NH2, For-beta-Abu-NH2, and For-beta-Aib-NH2. Stability and conformational properties of all three model systems were computed at different levels of theory: RHF/3-21G, B3LYP/6-311++G(d,p)//RHF/3-21G, B3LYP/6-311++G(d,p), MP2//B3LYP/6-311++G(d,p), CCSD//B3LYP/6-311++G(d,p), and CCSD(T)//B3LYP/6-311++G(d,p). In addition, ab initio E = E(phi, micro, psi) potential energy hypersurfaces of all three models were determined, and their topologies were analyzed to determine the inherent flexibility properties of these beta-peptide models. Fewer points were found and assigned than expected on the basis of Multidimensional Conformational Analysis (MDCA). Furthermore, it has been demonstrated, that the four-dimensional surface, E = E(phi, mu, psi), can be reduced into a three-dimensional one: E = E[phi, f(phi), psi]. This reduction of dimensionality of freedom of motion suggests that beta-peptides are less flexible than one would have thought. This agrees with experimental data published on the conformational properties of peptides composed of beta-amino acid residues.     

Structural elucidation of nitro-substituted five-membered aromatic heterocycles utilizing GIAO DFT calculations

The GIAO (Gauge Including Atomic Orbitals) DFT (Density Functional Theory) method is applied at the B3LYP/6-31+G(d,p)//B3LYP/6-31+G(d), B3LYP/6-311++G(d,p)//B3LYP/6-31+G(d), B3LYP/6-311+G (2d,p)//B3LYP/6-31+G(d) and B3LYP/6-311++G(d,p)//B3LYP/6-311++G(d,p) levels of theory for the calculation of proton and carbon chemicals shifts and coupling constants for 25 nitro-substituted five-membered heterocycles. Difference (1D NOE) spectra in combination with long-range gHMBC experiments were used as tools for the structural elucidation of nitro-substituted five-membered heterocycles. The assigned NMR data (chemical shifts and coupling constants) for all compounds were found to be in good agreement with theoretical calculations using the GIAO DFT method. The magnitudes of one-bond (1JCH) and long-range (nJCH, n>1) coupling constants were utilized for unambiguous differentiation between regioisomers of nitro-substituted five-membered heterocycles.     

A new scheme for determining the intramolecular seven-membered ring N-H...O=C hydrogen-bonding energies of glycine and alanine peptides

In this paper a new scheme was proposed to calculate the intramolecular hydrogen-bonding energies in peptides and was applied to calculate the intramolecular seven-membered ring N-H...O=C hydrogen-bonding energies of the glycine and alanine peptides. The density-functional theory B3LYP6-31G(d) and B3LYP6-311G(d,p) methods and the second-order Moller-Plesset perturbation theory MP26-31G(d) method were used to calculate the optimal geometries and frequencies of glycine and alanine peptides and related structures. MP26-311++G(d,p), MP26-311++G(3df,2p), and MP2/aug-cc-pVTZ methods were then used to evaluate the single-point energies. It was found that the B3LYP6-31G(d), MP26-31G(d), and B3LYP6-311G(d,p) methods yield almost similar structural parameters for the conformers of the glycine and alanine dipeptides. MP2/aug-cc-pVTZ predicts that the intramolecular seven-membered ring N-H...O=C hydrogen-bonding strength has a value of 5.54 kcal/mol in glycine dipeptide and 5.73 and 5.19 kcal/mol in alanine dipeptides, while the steric repulsive interactions of the seven-membered ring conformers are 4.13 kcal/mol in glycine dipeptide and 6.62 and 3.71 kcal/mol in alanine dipeptides. It was also found that MP26-311++G(3df,2p) gives as accurate intramolecular N-H...O=C hydrogen-bonding energies and steric repulsive interactions as the much more costly MP2/aug-cc-pVTZ does.     

Theoretical study of the conformational energy hypersurface of cyclotrisarcosyl

The multidimensional Conformational Potential Energy Hypersurface (PEHS) of cyclotrisarcosyl was comprehensively investigated at the DFT (B3LYP/6-31G(d), B3LYP/6-31G(d,p) and B3LYP/6-311++G(d,p)), levels of theory. The equilibrium structures, their relative stability, and the Transition State (TS) structures involved in the conformational interconversion pathways were analyzed. Aug-cc-pVTZ//B3LYP/6-311++G(d,p) and MP2/6-31G(d)//B3LYP/6-311++G(d,p) single point calculations predict a symmetric cis-cis-cis crown conformation as the energetically preferred form for this compound, which is in agreement with the experimental data. The conformational interconversion between the global minimum and the twist form requires 20.88 kcal mol-1 at the MP2/6-31G(d)//B3LYP/6-311++G(d,p) level of theory. Our results allow us to form a concise idea about the internal intricacies of the PEHSs of this cyclic tripeptide, describing the conformations as well as the conformational interconversion processes in this hypersurface. In addition, a comparative analysis between the conformational behaviors of cyclotrisarcosyl with that previously reported for cyclotriglycine was carried out     

Theoretical studies on the intramolecular hydrogen bond and tautomerism of 8-mercaptoquinoline in the gaseous phase and in solution using modern DFT methods

A theoretical quantum chemical study of the intramolecular hydrogen bonding interactions in 8-mercaptoquinoline has been carried out. Special attention has been paid to the rotation of S-H bond and intramolecular proton-transfer reactions. Therewith, the B3LYP/6-311++G(d,p), B3LYP/6-31+G(2d,2p), MPW1K/6-311++G(d,p), MPW1K/6-31+G(2d,2p), BH&HLYP/6-311++G(d,p), and G96LYP/6-311++G(d,p) methods have been used. By means of the Onsager and PCM reaction field methods, the effects of solvent on hydrogen-bond energies, conformational equilibria, rotational barriers, and tautomerism in aqueous solution have been studied. These simulations were done at the MPW1K/6-311++G(d,p) and B3LYP/6-311++G(d,p) levels. Natural-bond orbital analysis has been performed to study the intramolecular hydrogen bond (IHB) in the gaseous phase and in aqueous medium. The stability of forms under consideration in solution does not coincide with that in the gaseous phase, underlining a great importance of the electrostatic influence of solvent. Double-proton transfer in the prototropic tautomerization of 8-mercaptoquinoline, one water molecule complex in the gaseous phase and in solution, has been systematically studied. The double-proton transfer occurs concertedly and synchronously. The water-assisted tautomerization is kinetically less, but thermodynamically more favorable, compared to that of the single-proton transfer. As in the case with single-proton transfer, for water-assisted reaction, the tautomerization energies and barrier heights decrease with the increase in dielectric constant, which implies faster and more complete tautomerization of 8-mercaptoquinoline in a polar solvent.     

Inductive effects in radicals calculated from DFT energies; substituted bicyclo[2.2.2]octan-1-yloxy radicals

Energies of a series of 4-substituted 1-oxybicyclo[2.2.2]octan-1-yloxy radicals with 18 various substituents were calculated within the framework of the DFT theory at the levels UB3LYP/6-311+G(d,p)//UB3LYP/6-311+G(d,p) and UB3LYP/6-311++G(2df,p)//UB3LYP/6-311+G(d,p) and compared with similar series of the parent alcohols, their deprotonated and protonated forms calculated at the levels B3LYP/6-311+G(d,p)//B3LYP/6-311+G(d,p) and B3LYP/6-311++G(2df,p)//B3LYP/6-311+G(d,p). The two levels are of the same performance and both are sufficient for molecules of this type according to comparison with scarce experimental gas-phase acidities and basicities. The substituent effects were analyzed in terms of isodesmic equations. In addition to strong dependence on the substituent inductive effect, a slight dependence on the electronegativity of the first atom of the substituent was proven in certain cases. In all aspects, there is no qualitative difference between the effects on radicals and on similar closed shell species. Radicals behave as slightly electron deficient; the substituent effect is weaker than that on the ions but stronger than on neutral molecules.     

Theoretical Studies on the Aminolysis of Ester: Effects of the Catalysis and Substituent to the Reaction of Methyl Indole-3-acetate with Ammonia

A comprehensive exploration of the aminolysis mechanism for methyl indole-3-acetate with ammonia is carried out by employing the B3 LYP/6-311++G(d,p), M06-2 X/6-311++G(d,p) and MP2/6-311++G(d,p)//M06-2 X/6-311++G(d,p) levels. Two alterative reaction channels of the concerted and addition/elimination stepwise processes including the uncatalyzed, base-catalyzed reactions are taken into consideration. Subsequently, the substituent effects and solvent effects in methanol are also evaluated at the M06-2 X/6-311++G(d,p) level. The calculated results indicate that the calculated values of M06-2 X level are quite close to those of MP2, the stepwise pathway has more advantages to the concerted one for all of the reaction processes and the catalyst facilitates the proton migration and decreases the energy barriers as well. It is shown that the most preferred mechanism is the based-catalyzed stepwise process, the substituent of NH2 group slightly accelerates all the aminolysis reaction processes, and the solvent effect does not remarkably change the mechanism of the reaction.     

类硅烯H2C＝SiLiBr与RH (R＝F, OH, NH2)的插入反应

采用DFT B3LYP和QCISD方法研究了类硅烯H2C＝SiLiBr与RH (R＝F, OH, NH2)的插入反应. 在B3LYP/6- 311＋G(d,p)水平上优化了反应势能面上的驻点构型. 结果表明, H2C＝SiLiBr与HF, H2O或NH3发生插入反应的机理相同. QCISD/6-311＋＋G(d,p)//B3LYP/6-311＋G(d,p)计算的三个反应的势垒分别为148.62, 164.42和165.07 kJ•mol－1, 反应热分别为－69.63, －43.02和－28.27 kJ•mol－1. 相同条件下发生插入反应时, 反应活性都是H—F＞H—OH＞H—NH2.     

Theoretical investigation on the GaCl3-catalyzed ring-closing metathesis reaction of N-2,3-butadienyl-2-propynyl-1-amine: three-membered ring versus four-membered ring mechanism

The gallium chloride (GaCl(3))-catalyzed ring-closing metathesis reaction mechanism of N-2,3-butadienyl-2-propynyl-1-amine has been studied at the Becke three-parameter hybrid functional combined with Lee-Yang-Parr correlation functional (B3LYP)/6-31G(d), B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p)//B3LYP/ 6-31G(d) and the second-order M?ller-Plesset perturbation (MP2)/6-311++G(d,p)//B3LYP/6-31+G(d,p) levels. It was found that the final metathesis product can be yielded via a three-membered or four-membered ring mechanism. The three-membered ring pathway is favorable due to its low energy barrier at the rate determining step. The whole reaction is stepwise and strongly exothermic.     

Theoretical study of the N-H…O red-shifted and blue-shifted hydrogen bonds

Theoretical calculations are performed to study the nature of the hydrogen bonds in complexes HCHO…HNO, HCOOH…HNO, HCHO…NH3, HCOOH…NH3, HCHO…NH2F and HCOOH…NH2F. The geometric structures and vibrational frequencies of these six complexes at the MP2/6-31+G(d,p), MP2/6-311++G(d,p), B3LYP/6-31+G(d,p) and B3LYP/6-311++G(d,p) levels are calculated by standard and counterpoise-corrected methods, respectively. The results indicate that in complexes HCHO…HNO and HCOOH…HNO the N-H bond is strongly contracted and N-H…O blue-shifted hydrogen bonds are observed. While in complexes HCHO…NH3, HCOOH…NH3, HCHO…NH2F and HCOOH…NH2F, the N-H bond is elongated and N-H…O red-shifted hydrogen bonds are found. From the natural bond orbital analysis it can be seen that the X-H bond length in the X-H…Y hydrogen bond is controlled by a balance of four main factors in the opposite directions: hyperconjugation, electron density redistribution, rehybridization and structural reorganization. Among them hyperconjugation has the effect of elongating the X-H bond, and the other three factors belong to the bond shortening effects. In complexes HCHO…HNO and HCOOH…HNO, the shortening effects dominate which lead to the blue shift of the N-H stretching frequencies. In complexes HCHO…NH3, HCOOH…NH3, HCHO…NH2F and HCOOH…NH2F where elongating effects are dominant, the N-H…O hydrogen bonds are red-shifted.     

Electronic properties of multifurcated bent hydrogen bonds CH3...Y and CH2...Y

H-bonding angle angleYHX has an important effect on the electronic properties of the H-bond Y...HX, such as intra- and intermolecular hyperconjugations and rehybridization, and topological properties of electron density. We studied the multifurcated bent H-bonds of the proton donors H3CZ (Z = F, Cl, Br), H2CO and H2CF2 with the proton acceptors Cl(-) and Br(-) at the four high levels of theory: MP2/6-311++G(d,p), MP2/6-311++G(2df,2p), MP2/6-311++G(3df,3pd) and QCISD/6-311++G(d,p), and found that they are all blue-shifted. These complexes have large interaction energies, 7-12 kcal mol(-1), and large blue shifts, delta r(HC) = -0.0025 --0.006 A and delta v(HC) = 30-90 cm(-1). The natural bond orbital analysis shows that the blue shifts of these H-bonds Y...HnCZ are mainly caused by three factors: rehybridization; indirect intermolecular hyperconjugation n(Y) -->sigma*(CZ), in that the electron density from n(Y) of the proton acceptor is transferred not to sigma*(CH), but to sigma*(CZ) of the donor; intramolecular hyperconjugation n(Z) -->sigma*(CH), in that the electron density in sigma*(CH) comes back to n(Z) of the donor such that the occupancy in sigma*(CH) decreases. The topological properties of the electron density of the bifurcated H-bonds Y...H2CZ are similar to those of the usual linear H-bonds, there is a bond critical point between Y and each hydrogen, and a ring critical point inside the tetragon YHCH. However, the topological properties of electron density of the trifurcated H-bonds Y...H3CZ are essentially different from those of linear H-bonds, in that the intermolecular bond critical point, which represents a closed-shell interaction, is not between Y and hydrogen, but between Y and carbon.     

Theoretical and experimental study of the absolute configuration of helical structure of (2R,3S)-Rubiginone A2 analog

A novel analogue of (2R,3S)-Rubiginone A₂ was synthesized as a chiral helical model compound via an eight-step procedure (2.7% overall yield). Quantum methods, such as density functional theory (DFT) at different basis sets of 6-311+(d), 6-311++G(2d,p), were used to compute its optical rotation and electronic circular dichroism at the B3LYP/6-311++G(2d,p) level in the gas phase and in solution using PCM model, respectively. UV corrections were performed in electronic circular dichroism (ECD) simulations to match the experimental ECD well. The suitable computational methods, e.g., B3LYP/6-311++G(2d,p)//B3LYP/6-311++G(2d,p) in the gas phase using zero-point energy in Boltzmann statistics, were found and suggested for optical rotation and circular dichroism computations that can be used for absolute configuration determination of chiral helical compounds.     

Mechanistic insights into triterpene synthesis from quantum mechanical calculations. Detection of systematic errors in B3LYP cyclization energies

Most quantum mechanical studies of triterpene synthesis have been done on small models. We calculated mPW1PW91/6-311+G(2d,p)//B3LYP/6-31G* energies for many C30H51O+ intermediates to establish the first comprehensive energy profiles for the cationic cyclization of oxidosqualene to lanosterol, lupeol, and hopen-3beta-ol. Differences among these 3 profiles were attributed to ring strain, steric effects, and proton affinity. Modest activation energy barriers and the ample exothermicity of most annulations indicated that the cationic intermediates rarely need enzymatic stabilization. The course of reaction is guided by hyperconjugation of the carbocationic 2p orbital with parallel C-C and C-H bonds. Hyperconjugation for cations with a horizontal 2p orbital (in the plane of the ABCD ring system) leads to annulation and ring expansion. If the 2p orbital becomes vertical, hyperconjugation fosters 1,2-methyl and hydride shifts. Transition states leading to rings D and E were bridged cyclopropane/carbonium ions, which allow ring expansion/annulation to bypass formation of undesirable anti-Markovnikov cations. Similar bridged species are also involved in many cation rearrangements. Our calculations revealed systematic errors in DFT cyclization energies. A spectacular example was the B3LYP/6-311+G(2d,p)//B3LYP/6-31G* prediction of endothermicity for the strongly exothermic cyclization of squalene to hopene. DFT cyclization energies for the 6-311+G(2d,p) basis set ranged from reasonable accuracy (mPW1PW91, TPSSh with 25% HF exchange) to underestimation (B3LYP, HCTH, TPSS, O3LYP) or overestimation (MP2, MPW1K, PBE1PBE). Despite minor inaccuracies, B3LYP/6-31G* geometries usually gave credible mPW1PW91 single-point energies. Nevertheless, DFT energies should be used cautiously until broadly reliable methods are established.     

类锗烯H2C=GeLiCl与RH(R=F,OH, NH2)的插入反应

采用DFT B3LYP和QCISD方法研究了不饱和类锗烯H2C=GeLiCl与RH(R=F, OH, NH2)的插入反应. 在B3LYP/6-311+G(d,p)水平上优化了反应势能面上的驻点构型. 结果表明, H2C=GeLiCl与HF、H2O 或NH3发生插入反应的机理相同. QCISD/6-311++G(d,p)//B3LYP/6-311+G(d,p)计算的三个反应的势垒分别为173.53、194.48和209.05 kJ·mol-1, 反应热分别为60.18、72.93和75.34 kJ·mol-1. 相同条件下发生插入反应时, 反应活性顺序都是H—F>H—OH>H—NH2.     

复合物HNO···H2O2内N-H···O蓝移氢键的理论研究

A theoretical study on the blue-shifted H-bond N-H…O and red-shifted H-bond O-H…O in the complexHNO…H_2O_2 was conducted by employment of both standard and counterpoise-corrected methods to calculate thegeometric structures and vibrational frequencies at the MP2/6-31G(d),MP2/6-31 G(d,p),MP2/6-311 q G(d,p),B3LYP/6-31G(d),B3LYP/6-31 G(d,p) and B3LYP/6-311 G(d,p) levels.In the H-bond N-H…O,the calcu-lated blue shift of N-H stretching frequency is in the vicinity of 120 cm~(-1) and this is indeed the largest theoreticalestimate of a blue shift in the X-H…Y H-bond ever reported in the literature.From the natural bond orbital analy-sis,the red-shifted H-bond O-H…O can be explained on the basis of the dominant role of the hyperconjugation.For the blue-shifted H-bond N-H…O,the hyperconjugation was inhibited due to the existence of significant elec-tron density redistribution effect,and the large blue shift of the N-H stretching frequency was prominently due tothe rehybridization of sp~n N-H hybrid orbital.     

Investigation of 1H NMR chemical shifts of organic dye with hydrogen bonds and ring currents

The (1)H NMR chemical shifts were theoretically computed for the organic dyes 2-(2,6-dimethyl-4H-pyran-4-ylidene)-malononitrile (1), cyano-(2,6-dimethyl-4H-pyran-4-ylidene)-acetic acid methyl ester (2), 2-(2,6-bis(4-(dimethylamino)styryl)-4H-pyran-4-ylidene)-malononitrile (3), and methyl 2-(2,6-bis(4-(dimethylamino)styryl)-4H-pyran-4-ylidene)-2-cyanoacetate (4) at the GIAO/B3LYP/6-311++G(d,p)//B3LYP/6-311++G(d,p) level of theory. Moreover, the intramolecular rotational barriers of the molecules were calculated to evaluate the internal flexibility with respect to the torsional degrees of freedom, and the nuclear-independent chemical shifts (NICS) were employed to analyze the ring currents. The difference was explained in terms of intramolecular hydrogen bonds and ring currents of the molecules. The (1)H NMR spectra were reproduced by experiments for the comparison with computationally constructed data. Our results suggest a good guideline in interpreting (1)H NMR chemical shifts using computational methods and furthermore a reliable perspective for designing molecular structures.     

FT-IR and FT-Raman spectroscopic investigation, computed vibrational frequency analysis and IR intensity and Raman activity peak resemblance analysis on 2-nitroanisole using HF and DFT (B3LYP and B3PW91) calculations

Fourier-transform Raman and infrared spectra of 2-nitroanisole are recorded (4000-100 cm(-1)) and interpreted by comparison with respective theoretical spectra calculated using HF and DFT method. The geometrical parameters with C(S) symmetry, harmonic vibrational frequencies, infrared and Raman scattering intensities are determined using HF/6-311++G (d, p), B3LYP/6-311+G (d, p), B3LYP/6-311++G (d, p) and B3PW91/6-311++G (d, p) level of theories. A detailed vibrational spectral analysis has been carried out and assignments of the observed fundamental bands have been proposed on the basis of peak positions and relative intensities. The results of the calculations have been used to simulate IR and Raman spectra for the molecule that showed good agreement with the observed spectra. The SQM method, which implies multiple scaling of the DFT force fields has been shown superior to the uniform scaling approach. The vibrational frequencies and the infrared intensities of the C-H modes involved in back-donation and conjugation are also investigated.     

The molecular structure of fluoromalononitrile as determined by gas-phase electron diffraction and ab initio calculations

The molecular structure of fluoromalononitrile was studied by means of gas-phase electron diffraction and quantum mechanical methods using HF/6-31G(d), MP2/6-311++G(2df,2pd) and DFT/B3LYP/6-31G(d), B3PW91/6-31G(d), B3LYP/6-311++G(2df,2pd) and B3PW91/6-311++G(2df,2pd). The r(g) and angle(alpha) structural parameters we obtained from the present analysis are: CC=1.487(5) A, CN=1.157(3) A, CF=1.386(5) A, CH=1.096 A (ass.), angleCCC=106.7(1.0) degrees , angleCCF=108.0(0.7) degrees , angleCCN=177.6(2.0) degrees . Uncertainties in parenthesis are 3sigma.