Modeling NMR parameters by DFT methods as an aid to the conformational analysis of cis-fused 7a(8a)-methyl octa(hexa)hydrocyclopenta[d][1,3]oxazines and [3,1]benzoxazines

The energies of the preferred conformations of four 7a-methyl octa(or hexa)hydrocyclopenta[d][1,3]oxazines, five 8a-methyl octa(or hexa)hydro[3,1]benzoxazines, and 8a-methyl hexahydro[1,3]benzoxazinone, all cis-fused, were investigated by DFT methods. Following geometry optimization at the B3LYP/6-31G(d,p) level, both the proton chemical shifts and the vicinal coupling constants between H-4a and the H-4 and H-5 protons were calculated at the B3LYP/cc-pVTZ level and compared to the previously experimentally measured values. The agreement between the calculated and the experimental chemical shifts was found to be good. Similarly, the agreement between the calculated and the experimental vicinal coupling constants was also found to be good, thus providing a methodology for determining the conformational equilibria of such systems that is comparable in many respects to experimental approaches such as variable-temperature NMR or to the use of model coupling constant values, when available, from analogous compounds.     

已烯自由基阳离子的密度泛函理论研究

使用密度泛函理论B3LYP方法和6-31G(d,p),6-31+G(d,p),6-311G(d,p)及6- 311+G(d,p)基组,分别对1-C_6H_(12)~+,2-C_6H_(12)~+和3-C_6H_(12)~+的各种构 象进行了几何构型优化,并在B3LYP/6-311G(d,p)水平上进行了频率分析计算,在 各优化构型上,使用B3LYP和MP2(full)方法进行了超精细结构的计算。计算的3- C_6H_(12)~+的超精细偶合常数比以往的计算结果更好;1-C_6H_(12)~+和2-C_6H_ (12)~+的超精细偶合常数目前尚无实验数据报道,本计算预言了它们的超精细偶合 常数和最稳定构型。     

丁烯和丁烷自由基阳离子的分子结构和超精细结构的 理论研究

用密度函数理论B3LYP方法和6-31G(d,p),6-311G(d,p)及6-311+G(d,p)基组，分别对1-C4H^+~8,2-C4H^+~8和C4H^+~10进行了构型优化和频率分析计算，预言1-C4H^+~8具有非平面构型，与以往报道的从头算和密度函数理论计算结果不同。在各自由基阳离子的B3LYP构型上，进行了B3LYP、MP2及MRSDCI方法的超精细偶合常数计算，得到了比以往更好的结果，特别是MP2/B3LYP计算值是至今与实验值符合得最好的理论计算结果。     

Experimental and DFT 1H NMR study of conformational equilibria in trans-4',7-dihydroxyisoflavan-4-ol and trans-isoflavan-4-ol

The solution-state conformational equilibria of trans-4',7-dihydroxyisoflavan-4-ol (1) and trans-isoflavan-4-ol (2) were assessed based on the temperature dependence of their vicinal coupling constants J(H)(-)(2)(alpha)(,H)(-)(3) and J(H-3,H)(-)(4) in comparison to values calculated with density functional theory (DFT) methods at the B3LYP/cc-pVTZ//B3LYP/6-31G(d,p) level of theory. For each half-chair conformer, several rotamers with respect to the C-4 hydroxyl and C-3 phenyl were calculated and the overall diequatorial-to-diaxial ratio at 298 K was assessed as 66:34 for 1 and 73:27 for 2. The syntheses of 1 and 2 are described.     

Conformational analysis of sulfur-containing 6-deoxy-l-hexose derivatives by molecular modeling and NMR spectroscopy. A theoretical study and experimental evidence of intramolecular nonbonded interactions between sulfur and oxygen

6-Deoxy-l-mannose diphenyldithioacetal (1) unexpectedly gave the rearranged products phenyl 3,4-di-O-acetyl-2-S-phenyl-1,2-dithio-6-deoxy-beta-l-glucopyranoside (9) and 3,4-di-O-acetyl-2,5-anhydro-6-deoxy-l-glucose diphenyldithioacetal (10) upon treatment with acetyl chloride, while 6-deoxy-l-mannose ethylenedithioacetal (3) yielded (4aR,6S,7S,8R,8aS)-7,8-diacetyloxy-6-methylhexahydro-4aH-[1,4]dithiino[2,3b]pyran (11), whose structure was further confirmed by X-ray diffraction, and 3,4-di-O-acetyl-2,5-anhydro-l-rhamnose ethylenedithioacetal (12). The geometry of the four rearranged products as well as that of 1-thio-6-deoxy-l-mannopyranosides 5 and 7 and their acetyl derivatives 6 and 8 was studied by density functional theory (B3LYP/6-31G) molecular models, in combination with a Karplus-type analysis of the NMR vicinal coupling constants, revealing that the six-membered ring of pyranosides 5-9 and 11 exists in a slightly distorted chair conformation (6-13% distortion) and that the conformational behavior of the 2,5-anhydro-6-deoxy-l-glucose dithioacetals 10 and 12 is strongly influenced by the presence of stabilizing intramolecular nonbonded sulfur-oxygen 1,4- and 1,5-interactions. Compounds 9-12 were formed by a molecular rearrangement via sulfonium ion intermediates followed by stereoselective intramolecular cyclizations as formulated by the quantum chemical calculations performed in the present study.     

Synthesis,Crystal Structure and Spectroscopic Properties of 1,2-Benzothiazine Derivatives:An Experimental and DFT Study

1,2-Benzothiazine derivatives methyl 3-methoxy-4-oxo-3,4-dihydro-2H-benzo[e] [1,2]thiazine-3-carboxylate 1,1-dioxide(1) and methyl 2-ethyl-3-hydroxy-4-oxo-3,4-dihydro-2Hbenzo[e][1,2]thiazine-3-carboxylate 1,1-dioxide(2) were synthesized, and characterized by spectroscopic techniques; 1H-NMR and infrared(IR) spectroscopy. Crystals of 1 and 2 were grown by slow evaporation of methanol and ethyl acetate, respectively and their crystal structures were investigated by single-crystal X-ray diffraction analysis. Geometric properties were calculated by the B3 LYP method of density functional theory(DFT) at the 6-31G+(d) basis set to compare with the experimental data. Simulated properties were found in strong agreement with the experimental ones. Intermolecular forces have also been modeled in order to investigate the strength of packing and strong hydrogen bonding was observed in both compounds 1 and 2. Electronic properties such as Ionization Potential(IP), Electron Affinities(EA) and coefficients of the highest occupied molecular orbital(HOMO) and the lowest unoccupied molecular orbital(LUMO) of com- pounds 1 and 2 were simulated for the first time.     

Conformational study of (8alpha,8'beta)-bis(substituted phenyl)-lignano-9,9'-lactones by means of combined computational, database mining, NMR, and chemometric approaches

Beta-(3,4-Methylenedioxybenzyl)-gamma-butyrolactone (MDBL) and (-)-hinokinin (HK) were obtained by partial synthesis and characterized by 1H NMR and computational methods (conformational analysis, molecular modeling, structural data mining and chemometrics). Three conformers were detected for MDBL and nine were found for HK. The energy differences are around 1 and 2 kcal mol(-1) and rotation barriers are less than 3 and 5 kcal mol(-1) for MDBL and HK conformers, respectively. The geometries of these conformers, obtained from semiempirical PM3 and density functional theory (DFT) B3LYP 6-31G** calculations agree satisfactorily with 1H NMR data (vicinal proton-proton coupling constants) and structures retrieved from the Cambridge Structural Database (torsion angles). DFT combined with some variants of the Haasnoot-de Leeuuw-Altona equations gives the best predictions for the coupling constants. The molecular conformation of MDBL, of HK, and of related systems depends not only on intramolecular interactions but also on crystal packing forces and solvent-solute interactions, in particular hydrogen bonds and polar interactions. Hydration favors more stable HK conformers, which can be important for their behavior in chemical and biological systems.     

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.     

Quantum chemical computations of 1,3-phenylenediacetic acid

Molecular geometry, vibrational wavenumbers and gauge including atomic orbital (GIAO) 13C NMR and 1H NMR chemical shift values of 1,3-phenylenediacetic acid (C10H10O4), in the ground state have been calculated by using ab initio Hartree-Fock (HF) and density functional theory (DFT/B3LYP) methods with 6-311++G(d,p) basis set for the first time. Comparison of the observed fundamental vibrational modes of 1,3-phenylenediacetic acid and calculated results by density functional B3LYP and Hartree-Fock methods indicate that B3LYP is superior to the scaled Hartree-Fock approach for quantum chemical studies. Geometric parameters (bond lengths and bond angles) and vibrational wavenumbers obtained by the HF and DFT/B3LYP methods are in good agreement with the experimental data. Furthermore, this is the first time the results of the calculated JCH and JCC coupling constants of the C10H10O4 molecule are presented in this study.     

Solvent effects on isomerization and spectral properties of photochromic-switching diarythene derivatives in polar and apolar solutions

The photocyclization behavior and dynamic conformational transition of photochromic switches of diarythene derivatives in solutions are investigated by using the density functional theory (DFT) and molecular dynamics (MD) simulations. Three possible conformations, antiparallel (anti), parallel (para), and twist, for the open-ring isomers of 1,2-bis(2-methylbenzothiophene-3-yl)maleic anhydride are located. Both PCM-B3LYP/6-31G* calculations and MD simulations demonstrate that anti and twist open-ring isomers can interconvert freely in n-hexane and acetonitrile solutions at room temperature. The statistical ratio of twist to anti isomers from MD simulations is 2.09 in n-hexane and 1.07 in CH(3)CN, in qualitative agreement with those (1.18 in n-hexane and 1.05 in CH(3)CN) estimated from Arrhenius analysis of DFT activation energies. The solvent polarity has little influence on the isomerization of open-ring isomers in the ground state. Due to the evident charge transfer upon excitations, the solvent effects on the electronic structures and absorption spectra of low-lying excited states (S(1) and S(2)) are more significant. For such charge-transfer excited states, the long-range corrected functional CAM-B3LYP gives better agreement with the experimental spectra than B3LYP. The solvent polarity and polarization of the charge-transfer excited states are crucial for fabricating the novel functionalized photochromic molecular switches.     

Vicinal 1H-1H NMR coupling constants from density functional theory as reliable tools for stereochemical analysis of highly flexible multichiral center molecules

A protocol for stereochemical analysis, based on the systematic comparison between theoretical and experimental vicinal (1)H-(1)H NMR coupling constants, was developed and applied to a series of flexible compounds (1-8) derived from the 6-heptenyl-5,6-dihydro-2H-pyran-2-one framework. The method included a broad conformational search, followed by geometry optimization at the DFT B3LYP/DGDZVP level, calculation of the vibrational frequencies, thermochemical parameters, magnetic shielding tensors, and the total NMR spin-spin coupling constants. Three scaling factors, depending on the carbon atom hybridizations, were found for the (1)H-C-C-(1)H vicinal coupling constants: f((sp3)-(sp3)) = 0.910, f((sp3)-(sp2)) = 0.929, and f((sp2)-(sp2))= 0.977. A remarkable correlation between the theoretical (J(pre)) and experimental (1)H-(1)H NMR (J(exp)) coupling constants for spicigerolide (1), a cytotoxic natural product, and some of its synthetic stereoisomers (2-4) demonstrated the predictive value of this approach for the stereochemical assignment of highly flexible compounds containing multiple chiral centers. The stereochemistry of two natural 6-heptenyl-5,6-dihydro-2H-pyran-2-ones (14 and 15) containing diverse functional groups in the heptenyl side chain was also analyzed by application of this combined theoretical and experimental approach, confirming its reliability. Additionally, a geometrical analysis for the conformations of 1-8 revealed that weak hydrogen bonds substantially guide the conformational behavior of the tetraacyloxy-6-heptenyl-2H-pyran-2-ones.     

Indirect nuclear 57Fe-13C and 57Fe-1H spin-spin coupling in ferrocenes and cyclopentadienyliron complexes: measurements and DFT calculations

Several 57Fe-labeled ferrocene derivatives and other cyclopentadienyliron complexes were studied by 57Fe and 13C NMR with respect to isotope-induced chemical shifts 1Delta12/13C(57Fe) and the magnitude and sign of coupling constants 1J(57Fe,13C) and 2J(57Fe,1H). The geometries of the complexes were optimized by DFT methods [B3LYP/6-311+G(d,p)] and chemical shifts (GIAO) and coupling constants were calculated at the same level of theory. The trends in calculated 57Fe nuclear shielding agree well with the experimental data and, in the case of coupling constants 1J(57Fe,13C) and 2J(57Fe,1H), both sign and magnitude are correctly reproduced.     

戊烯自由基阳离子的密度泛函理论研究

使用密度泛函理论B3LYP方法和6-31G(d,p)、6-31＋G(d,p)、6-311G(d,p)及6-311＋G(d,p)基组，分别对2-C5H10＋和1-C5H10＋的各种构象进行了几何构型优化，并用B3LYP/6-311G(d,p)进行了频率分析计算.计算预言1-C5H10＋具有非平面构型，与以往报导的从头算计算结论相反.在两个自由基阳离子的各种构象的B3LYP几何构型上，进行了B3LYP和UMP2(full)方法的超精细偶合常数计算，得到了比以往更好的结果.     

Ab initio and density functional theory studies of the structure,gas-phase acidity and aromaticity of tetraselenosquaric acid

Results of ab initio self-consistent-field (SCF) and density functional theory (DFT) calculations of the gas-phase structure,acidity (free energy of deprotonation,G0) and aro-maticity of tetraselenosquaric acid (3,4-diselenyl-3-cy-dobutene-1,2-diselenone,H2C4Se4) are reported.The global minimum found on the potential energy surface of tetraselenosquaric acid presents a planar conformation.The ZZ iso-mer was found to have the lowest energy among the three planar conformers and the ZZ and ZE isomers are very dose in energy.The optimized geometric parameters exhibit a bond length equalization relative to reference compounds,cyclobu-tanediselenone,and cydobutenediselenol.The computed aromatic stabilization energy (ASE) by homodesmotic reaction is -77.4 (MP2(fu)/6 - 311 G //RHF/6 - 311 G) and - 54.8 kJ/mol (B3LYP/6 - 311 G //B3LYP/6 -311 G).The aromaticity of tetraselenosquaric add is indicated by the calculated diamagnetic susceptibility exaltation (A) - 19.13 (CSGT(IGAEM) - RHF/6 - 311 G// RHF/6-     

Anharmonic analysis of the vibrational spectra of some cyanides and related molecules of astrophysical importance

A detailed analysis of the vibrational spectra of carbonyl cyanide, diethynyl ketone and acetyl cyanide has been conducted in harmonic and anharmonic approximations. RHF, MP2 and density functional theory (DFT) methods with 6-311++G(2df,2p) basis sets and B3LYP functionals have been employed. Spectroscopic constants such as anharmonicity constants, rotational and centrifugal distortion constants, rotation-vibration coupling constants and Coriolis coupling coefficients have been calculated for each molecule and compared with the experimental data, where available. A close agreement between the calculated and experimental values of the spectroscopic constants has been obtained. Complete assignments have been provided to the fundamental bands, overtones and combination tones of the molecules. Density functional theory based anharmonic frequencies compare well with the experimental frequencies within +/-18 cm(-1) on an average. RHF and MP2 methods, however, give much higher values for the frequencies that need scaling even in the anharmonic approximation.     

Structure and ethanol complexation of cyclic tetrasaccharide in aqueous solution studied by NMR and molecular mechanics

The structure and ethanol complexation of a cyclic tetrasaccharide (CTS) in aqueous solution were investigated by proton NMR spectroscopy and molecular mechanics calculations. Two glucose units, A and B, of CTS are alternatively bonded by alpha-1,3 and alpha-1,6 linkages. The overlapped signals of protons A5, A6S, A6R, B3, B6S and B6R were resolved by spectral simulations to determine their chemical shifts and vicinal coupling constants. All vicinal coupling constants except for the A5-A6 spin system are consistent with the dihedral angles in the X-ray crystal structure. Each of protons A5, A6S, and A6R in the two units of A is equivalent with respect to the chemical shift. The vicinal coupling constants of (3)J(5-6S) and (3)J(5-6R) for unit A are close to the average of two rotamers that are present in crystals. The intensities of cross-peaks in the rotating frame nuclear Overhauser effect spectroscopy (ROESY) spectrum were rather well correlated with the effective distances calculated for the X-ray structure and molecular mechanics structures calculated in vacuo and water, although they are slightly better correlated with molecular mechanics structure in vacuo than with the other structures. From the changes of the chemical shifts of several CTS protons with increasing ethanol concentration, it was suggested that adsorption sites of ethanol on the plate structure of CTS are protons B2 and B4 (site B) in the concave face side and protons A1 and A2 (site A) in the convex back side. The binding constants for sites A and B are 0.0061 and 0.0176 M(-1), respectively. These binding constants are much smaller than a value of 4.1 M(-1) for the ethanol-alpha-cyclodextrin complex.     

Determination of the absolute configuration of [3(2)](1,4)barrelenophanedicarbonitrile using concerted time-dependent density functional theory calculations of optical rotation and electronic circular dichroism

The technique of time-dependent density functional theory (TDDFT) has very recently been applied to the calculation of both transparent spectral region optical rotations and electronic circular dichroism (CD). Here, we report the concerted application of the new methodologies to the determination of the absolute configuration (AC) of [3(2)](1,4)barrelenophanedicarbonitrile, 1, the first optically active barrelenophane. 1 is conformationally flexible: the two three-carbon bridges of 1 can each exhibit two conformations, leading to three inequivalent conformations of 1: a, b, and c. Conformational structures and energies are predicted using DFT at the B3LYP/6-31G level. Comparison of the calculated structures to structures obtained via X-ray crystallography of (+)-1 shows that (remarkably) all three conformations a-c are simultaneously present in crystalline (+)-1. The sodium D line specific rotations, [alpha](D), and CD spectra of a-c are calculated using TDDFT at the B3LYP/aug-cc-pVDZ level. Comparison of the conformationally averaged specific rotation and CD spectrum to the experimental data of Matsuda-Sentou and Shinmyozu leads to the AC 9S,12S(+)/9R,12R(-). The same AC is obtained both from [alpha](D) and from the CD, strongly supporting its reliability.     

Experimental and theoretical investigation of epoxide quebrachitol derivatives through spectroscopic analysis

Two synthetic epoxide derivatives, important intermediates in organic synthesis, were obtained from L-quebrachitol, and their conformations were proposed based on spectroscopic analysis. Density functional theory (DFT) calculations of infrared and NMR spectra were shown to be reliable enough for organic chemistry applications. The observed structures were determined with the aid of the DFT spectroscopic data, stressing the relevance and utility of combined experimental/theoretical studies and also the usefulness of the (13)C NMR B3LYP/6-31G(d,p) calculations.     

DFT and NMR parameterized conformation of valeranone

A Monte Carlo random search using molecular mechanics, followed by geometry optimization of each minimum energy structure employing density functional theory (DFT) calculations at the B3LYP/6-31G* level and a Boltzmann analysis of the total energies, generated accurate molecular models which describe the conformational behavior of the antispasmodic bicyclic sesquiterpene valeranone (1). The theoretical H-C-C-H dihedral angles gave the corresponding 1H, 1H vicinal coupling constants using a generalized Karplus-type equation. In turn, the 3J(H,H) values were used as initial input data for the spectral simulation of 1, which after iteration provided an excellent correlation with the experimental 1H NMR spectrum. The calculated 3J(H,H) values closely predicted the experimental values, excepting the coupling constant between the axial hydrogen alpha to the carbonyl group and the equatorial hydrogen beta to the carbonyl group (J(2beta, 3beta)). The difference is explained in terms of the electron density distribution found in the highest occupied molecular orbital (HOMO) of 1. The simulated spectrum, together with 2D NMR experiments, allowed the total assignment of the 1H and 13C NMR spectra of 1.     

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.