Steric effects of polar substituents evaluated in terms of energy by means of isodesmic reactions

Steric effects of various polar and some charged groups were estimated on sterically crowded cyclopropane cis-1,2-bis derivatives 2 or 3, in which the variable substituent is in the proximity of a t-butyl group or of a methyl group. The steric energy was evaluated with reference to the pertinent mono derivatives, that is as reaction energy of an isodesmic reaction, in which the crowded compound is formally synthesized from simple derivatives. Energies were calculated within the framework of the density functional theory at level B3LYP/6-311+G(d,p)//B3LYP/6-311+G(d,p) for 11 dipolar and 5 charged substituents. Interaction of charged substituents is not only steric (destabilizing) but also inductive (stabilizing). The steric effects evaluated in this way differ distinctly from the standard steric constants derived purely from the van der Waals radii of the substituents.     

Steric effects of the alkyl groups: evaluation on isolated molecules by means of isodesmic reactions

Several possible scales of steric effects of the alkyl groups were suggested on the basis of isodesmic model reactions, in which a sterically crowded compound is formally synthesized from simpler derivatives. The reaction energies were calculated within the framework of the density functional theory at the level B3LYP/6-311+G(d.p)//B3LYP/6-311+G(d.p) for 6 model systems and 7 various alkyl groups. The most important systems were cis-1,2-dialkylcyclopropanes 1 synthesized from two mono derivatives and sterically crowded derivatives of bicyclo[2.2.2]octane 2 with C(3) symmetry. The scales of steric effects evaluated from the two models were rather different: the first scale depended in effect only on the C atoms in the alpha and beta positions and the effects were almost equal for all primary alkyls. The second scale depended also on the gamma position and the effect of the CH(2)-t-Bu group was much greater than that of the ethyl group. Any relationship between various systems was found rarely, only in the case of very similar reaction series; even in such cases the relationship was sometimes linear, sometimes distinctly curvilinear. It is concluded that any universal scale of steric effects is in principle not possible since these effects depend specifically on the surroundings of the substituent in a particular reaction. Nevertheless, there is a similarity between various scales; a bulky group appears as bulky in any scale. Therefore, very rough correlations of steric effects are possible.     

Strong intramolecular hydrogen bonds and steric effects involving C═S groups: An NMR and computational study

A number of 5-acyl rhodanines and thiorhodanines with bulky acyl groups (pivaloyl and adamantoyl), not previously available, have been synthesized. The compounds are shown to exist in the enol form. Structures have been calculated using both the MP2 approach and the B3LYP-GD3BJ functional and the 6-311++G(d,p) basis set. Hydrogen bond energies are estimated by subtracting energies of a structure with the OH group turned 180° from those of the intramolecularly hydrogen-bonded one. Properties such as OH chemical shifts, two-bond isotope effects on ¹³C chemical shifts, electron densities at the bond critical point from atoms in molecules analysis, and the hydrogen bond energies show that the sterically hindered compounds have stronger hydrogen bonds than methyl or isopropyl derivatives. The combination of oxygen and sulfur derivatives enables a detailed analysis of hydrogen bond energies.     

Photooxidation of methylnaphthalenes

Studies on the photooxidation of methyl-substituted aromatic hydrocarbons have revealed that whereas electron density is a determinant of endoperoxide formation, steric factors are most important in influencing the stability of the endoperoxide. Additional information on the energetics of the reactions and on the magnitude of the steric interactions was obtained using calculations at the B3LYP/6-311+G* level of theory.     

GIAO/DFT evaluation of 13C NMR chemical shifts of selected acetals based on DFT optimized geometries

DFT/B3LYP calculations of the ground-state conformation of eight cyclic and acyclic acetals are presented and compared with experimental data. Results of single-point GIAO/DFT calculations at five different levels of theory show that isotropic shieldings need to be empirically scaled to achieve agreement with experimental chemical shifts. Statistical evaluation of data indicates that the most accurate prediction of 13C chemical shifts is achieved at the MPW1PW91/6-311G** level of theory. An empirical equation describing the relationship between delta values and shielding constants is postulated. This equation has been applied to the non-chair ground-state conformation of the six-membered acetonide and to the conformationally flexible benzodioxonine derivative. The agreement observed between the experimental and predicted chemical shifts shows that calculations at the MPW1PW91/6-311G** level of theory are adequate for addressing questions of conformation.     

Theoretical investigation on multinuclear NMR chemical shifts of some tris(trifluoromethyl)boron complexes

Tris(trifluoromethyl)boron complexes have unusual properties and may find applications in many fields of chemistry, biology, and physics. To gain insight into their NMR properties, the isotropic ¹¹B, ¹³C, and ¹⁹F NMR chemical shifts of a series of tris(trifluoromethyl)boron complexes were systematically studied using the gauge‐included atomic orbitals (GIAO) method at the levels of B3LYP/6‐31 + G(d,p)//B3LYP/6‐31G* and B3LYP/6‐311 + G(d,p)//B3LYP/6‐311 + G(d,p). Solvent effects were taken into account by polarizable continuum models (PCM). The calculated results were compared with the experimental values. The reason that the structurally inequivalent fluorine atoms in a specific species give a same chemical shift in experimental measurements is attributed to the fast rotation of CF₃ group around the B? C(F₃) bond because of the low energy barrier. The calculated ¹¹B, ¹³C(F₃), and ¹⁹F chemical shifts are in good agreement with the experimental measurements, while the deviations of calculated ¹³C(X, X = O, N) chemical shifts are slightly large. For the latter, the average absolute deviations of the results from B3LYP/6‐311 + G(d,p)//B3LYP/6‐311 + G(d,p) are smaller than those from B3LYP/6‐31 + G(d,p)//B3LYP/6‐31G*, and the inclusion of PCM reduces the deviation values. The calculated ¹⁹F and ¹¹B chemical shieldings of (CF₃)₃BCO are greatly dependent on the optimized structures, while the influence of structural parameters on the calculated ¹³C chemical shieldings is minor. Copyright © 2009 John Wiley & Sons, Ltd.     

Anisotropic and hydrogen bonding effects in phenylglyoxamides and mandelamides: theoretical and NMR conformational evaluation

Interesting anisotropic effects were observed for phenylglyoxamides and their respective mandelamides. Such effects were observed in experimental (1)H and (13)C NMR (in CDCl(3), CD(3)OD, and DMSO-d(6) solvents) and in some cases with good correlation to theoretical (1)H and (13)C NMR DFT-GIAO (B3LYP/6-311++G**//B3LYP/6-31G*) calculations. A systematic conformational analysis of these compounds was performed in a two-step methodology, using PM3 and DFT (B3LYP/6-31G*) calculations; with good accomplishment and computational time economy. It was observed that intramolecular hydrogen bonding plays a significant role in the conformation of such compounds. Finally, a geminal nonequivalence of an N-CH(2) moiety, in one of the alkyl side chain (R1 = R2), was found for the tertiary mandelamides studied.     

A GIAO/DFT study of 1H, 13C and 15N shieldings in amines and its relevance in conformational analysis

The 1H, 13C and 15N absolute shieldings of 13 amines were calculated at the GIAO/B3LYP/6-311++G** level. For some compounds (ethylamine, piperidine and 1-methylpiperidine) two conformations were calculated. The 13C and 15N data could be correctly correlated with experimental chemical shifts, allowing the conformation of 1-methylpiperidine to be established. The 1H NMR absolute shieldings, although less well correlated with delta values, were used to account for the anisotropy effects of the N lone pair.     

Conformational equilibria of trans -3-X-cyclohexanols (X = Cl, Br, CH3 and OCH3). A low temperature NMR study and theoretical calculations

The integration of 1H and 13C NMR spectra, at - 90 degrees C in CS2/CD2Cl2 (9:1), for the trans-3-chlorocyclohexanol (1), trans-3-bromocyclohexanol (2), and trans-3-methoxycyclohexanol (4) showed that the equatorial-axial (ea) conformer occurs as ca 63, 63, and 69% in the conformational equilibrium, respectively. This corresponds to the following DeltaG(ea-ae) values (from (1)H spectrum): - 0.32 +/- 0.01, - 0.32 +/- 0.04, - 0.48 +/- 0.05 kcal mol(-1); and to (from 13C spectrum): - 0.31 +/- 0.04, - 0.35 +/- 0.05, and - 0.44 +/- 0.01 kcal mol(-1), respectively, in very good agreement within both series. Thus, although bromine is bulkier than chlorine, the 1,3-diaxial steric effects are similar in these equilibria. However, the integration of (1)H NMR spectrum for the trans-3-methylcyclohexanol (3) gave 90% of the 3ae conformer in the equilibrium, at - 90 degrees C on CS2/CD2Cl2 (9:1), corresponding to a DeltaG(ea-ae) value of 1.31 +/- 0.02 kcal mol(-1). The values obtained through the additivity rule, with data from monosubstituted cyclohexanes (DeltaG(Ad) = DeltaG(X) + DeltaG(OH)), for compounds 1, 2, and 4 (-0.37 +/- 0.15, - 0.34 +/- 0.09, and - 0.46 +/- 0.04 kcal mol(-1), respectively) are in very good agreement with the experimental values, but it is significantly smaller for compound 3 (0.79 +/- 0.02 kcal mol(-1)). Theoretical calculations through different levels of theory (HF/6-311 + g**, B3LYP/6-311 + g**, MP2/6-31 + g**, and CBS-4M) showed that CBS-4M is the best method for the study of conformational equilibria for these systems, since it provides DeltaG(ea-ae) values similar to the experimental values.     

Conformations and rotational barriers of 2,2′‐bithiazole and 4,4′‐dimethyl‐2,2′‐bithiazole semiemperical,ab initio,and density functional theory calculations

Conformational properties of 2,2′‐bithiazole and 4,4′‐dimethyl‐2,2′‐ bithiazole have been studied by using AM1 and PM3 semiemperical methods and ab initio HF/6‐311+G* and B3LYP/6‐311+G* calculations. All methods agree that the planar s‐trans conformation is the global minimum and the perpendicular conformation is the transition state. Additional local minima were found using the Hartree–Fock (HF) and B3LYP levels for 2,2′‐bithiazole while for 4,4′‐dimethyl derivative the minima was located only at the MP2//B3LYP level. The barrier heights for rotation are 1.72, 7.69, and 7.88 kcal/mol at the PM3, HF, and B3LYP levels, respectively, and methyl substitution did not affect appreciably this value. Fourier expansion terms and bond orders were used to explain the origins of the rotational barrier in terms of π conjugation, electrostatic interaction, and steric effects, which represent the main factors in the shape of the rotational barrier. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 79: 367–377, 2000     

(13)C NMR spectral assignment of 1,4-diarylpiperazinones

The piperazinone derivatives have potential application in the pharmaceutical, polymer and textile fields. The present work describes the preparation of a series of new 1,4-diarylsubstituted-2-piperazinones by condensation of substituted N,N'-bis-(2-hydroxyphenyl)-ethylenediamines with glyoxal and the complete (13)C NMR spectral assignment accomplished using APT, HMQC and HMBC techniques. Substituent chemical-shift effects (SCS) were calculated, which showed different values for the lactam- and amine-substituted aromatic rings. The results show that predictions based on SCS effects are not simple for these molecules due to electronic and steric effects. Moreover, in the case of the ortho-substituted derivative 2 g, the NMR spectra reveal a dynamic behavior related to restricted rotation of the phenyl groups (atropisomerism).     

Conformation preference and related intramolecular noncovalent interaction of selected short chain chlorinated paraffins

Short chain chlorinated paraffins(SCCPs) are not only research focus of environmental issues but also interesting model molecules for organic chemistry which exhibit diverse conformation preference and intramolecular noncovalent interactions(NCIs). A systematic study was conducted to reveal the conformation preference and the related intramolecular NCIs in two C_(10)-isomers of SCCPs, 5,5,6,6-tetrachlorodecane and 4,4,6,6-tetrachlorodecane. The overall conformation profile was determined on the basis of relative energies calculated at the MP2/6-311++G(d,p) level with the geometries optimized by B3LYP/6-311++G(d,p) method. Then, quantum theory of atoms in molecules(QTAIM) has been adopted to identify the NCIs in the selected conformers of the model molecules at both B3LYP/6-311++G(d,p) and M06-2X/aug-cc-pvdz level. Different chlorine substitution modes result in varied conformation preference. No obvious gauche effect can be observed for the SCCPs with chlorination on adjacent carbon atoms. The most stable conformer of 5,5,6,6-tetrachlorodecane(t Tt) has its three dihedral angles in the T configuration, and there is no intramolecular NCIs found in this molecule. On the contrary, the chlorination on interval carbon atoms favors the adoption of gauche configuration for the H–C–C–Cl axis. Not only intramolecular H···Cl contacts but also H···H interactions have been identified as driving forces to compensate the instability from steric crowding of the gauche configuration. The gggg and g′g′g′g′ conformers are the most popular ones, while the populations of tggg and tg′g′g′ conformer are second to those of the gggg and g′g′g′g′ conformers. Meanwhile, the M06-2X method with large basis sets is preferred for identification of subtle intramolecular NCIs in large molecules like SCCPs.     

The Transition States for CO2 Capture by Substituted Ethanolamines

Quantum chemical studies are used to understand the electronic and steric effects on the mechanisms of the reaction of substituted ethanolamines with CO₂. SCS‐MP2/6‐311+G(2d,2p) calculations are used to obtain the activation energy barriers and reaction energies for both the carbamate and bicarbonate formation. Implicit solvent effects are included with the universal solvation model SMD. Carbamate formation is more favorable than bicarbonate formation for monoethanolamine (MEA) both kinetically and thermodynamically. Increase of the steric hindrance on the C atoms around the N atom in substituted ethanolamines favors bicarbonate formation over carbamate formation with lower activation barriers and thereby higher reaction rates. In contrast, substitution by an N‐methyl or N‐ethyl group on MEA leads to a lower activation barrier for both carbamate formation and bicarbonate formation. As a result, higher reaction rates are expected as compared to MEA, and therefore these compounds have significant potential as industrial CO₂ capturing solvents.     

Regio- and stereoselectivity in the Diels–Alder reaction of vinylallenes with acrolein: an ab initio study

To explore the regio- and stereoselectivity in Diels–Alder reactions of vinylallenes with acrolein, a parent vinylallene/acrolein system and a methyl-substituted vinylallene/acrolein system were studied. Ab initio calculations were used to identify eight transition state structures for each of the two Diels–Alder reactions at various computational levels (RHF/6-31G(d), RHF/6-311G(2d), B3LYP/6-311G(2d), and MP2/6-311G(2d)). The relative energies of the endo and exo transition states along with the regioselectivity have been determined from these calculations. In the parent vinylallene/acrolein system, the endo s-cis transition structure is the preferred stereoselectivity at all levels of theory, however, there is no regioselectivity. In the methyl substituted vinylallene/acrolein system, the endo s-trans transition state tends to compete with the endo s-cis transition state at the RHF levels of theory and is 1 kcal/mol more stable at the B3LYP/6-311G(2d) level of theory. Also, in the methyl-substituted system, there is now a definite preference for one regioisomer over the other. Both Diels–Alder reaction systems are asynchronous with the methyl-substituted system being more pronounced.     

Remarks on GIAO-DFT predictions of 13C chemical shifts

Predicting (13)C chemical shifts by GIAO-DFT calculations appears to be more accurate than frequently expected provided that: (a) the comparison between experimental and theoretical data is performed using the linear regression method, (b) a sufficiently high level of theory [e.g. B3LYP/6-311 + + G(2d,p)//B3LYP/6-311 + + G(2d,p) or PBE1PBE/6-311 + G(2df,p)//B3LYP/6-311 + + G(2d,p)] is used, (c) the experimental data originate from the measurements performed in one solvent whose influence is taken into account at the molecular geometry optimization step and, first of all, during the shielding calculation, (d) the experimental data are free of heavy atom effects or such effects are appropriately treated in calculations, and finally (e) the conformational compositions of the investigated objects are known.     

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.     

Substituent effects on 13C chemical shifts of ketenimines: a GIAO/HF and DFT study

GIAO/HF and DFT methods were utilized to predict the ¹³C chemical shifts of substituted ketenimines. GIAO HF/6–311+G(2d,p) and B3LYP/6–311+G(2d,p) methods were applied on the optimized B3LYP/6–31G(d) geometries and ¹³C chemical shifts of C_α and C_β of substituted ketenimines were correlated with group electronegativities. HF and DFT calculations indicated that increasing substituent group electronegativity leads to increasing chemical shift of C_β of substituted ketenimines, whereas the C_α values decrease. Copyright © 2003 John Wiley & Sons, Ltd.     

Quantum theory of atoms in molecules analysis on the conformational preferences of vinyl alcohol and related ethers

Vinyl alcohol, methyl vinyl ether, and tert-butyl vinyl ether were studied within the framework of the quantum theory of atoms in molecules at the B3LYP/6-311++G(2d,2p) level. Local and integrated properties of the charge density indicate that the anti conformational preference of the tert-butyl derivative is not due to a differing resonance contribution with regard to the less bulky vinyl ethers but to steric effects. There is a small delocalization of charge density, either total or pi, between oxygen and the terminal vinyl carbon, which does not support the resonance picture of vinyl compounds.     

Computation and analysis of 19F substituent chemical shifts of some bridgehead‐substituted polycyclic alkyl fluorides

The ¹⁹F NMR shieldings for several remotely substituted rigid polycyclic alkyl fluorides with common sets of substituents covering a wide range of electronic effects were calculated using the DFT‐GIAO theoretical model. The level of theory, B3LYP/6–311+G(2d,p), was chosen based on trial calculations which gave good agreement with experimental values where known. The optimized geometries were used to obtain various molecular parameters (fluorine natural charges, electron occupancies on fluorine of lone pairs and of the C? F bond, and hybridization states) by means of natural bond orbital (NBO) analysis which could help in understanding electronic transmission mechanisms underlying ¹⁹F substituent chemical shifts (SCS) in these systems. Linear regression analysis was employed to explore the relationship between the calculated ¹⁹F SCS and polar substituent constants and also the NBO derived molecular parameters. The ¹⁹F SCS are best described by an electronegativity parameter. The most pertinent molecular parameters appear to be the occupation number of the NBO p‐type fluorine lone pair and the occupation number of the C? F antibonding orbital. This trend suggests that in these types of rigid saturated systems hyperconjugative interactions play a key role in determining the ¹⁹F SCS. Electrostatic field effects appear to be relatively unimportant. Copyright © 2003 John Wiley & Sons, Ltd.     

2-丁基-四氢噻吩亚砜13C-NMR的理论研究

在RHF/6-31G和B3LYP/6-31G水平上对顺式(Cis-)与反式(Trans-)2-丁基-四氢噻吩亚砜(BTHTO)进行几何优化,应用规范不变原子轨道法(GIAO)在6-31G、6-31+G、6-31++G和6-31+G(2d,p)水平上计算了Cis-和Trans-BTHTO的13C-NMR,对13C-NMR谱进行了归属。结果表明,BTHTO噻吩五元环的稳定构象呈半椅式,Cis-和Trans-BTHTO中与硫原子直接碳原子13C-NMR的显著差异主要是由于空间构型不同引起分子的静电势场对相应碳原子的屏蔽作用不同所致。