Electronic structure and conformational properties of (carboxy-alkenyl)-phosphonic acids

The general conformational properties and electronic structure of (carboxy-alkenyl)-phosphonic derivatives were determined at RHF/STO-3G^* level. In all the series, low rotation barriers were found for the two C=C/P=O conformers. In the compounds in which the interactions between the carboxylic and phosphonic moieties are smaller, the most stable conformers are the C=C/P=O s-cis ones. In most of the conformers, the C=C/C=O system presents the disposition s-cis. The Z-(2-carboxy-vinyl) and Z-(2-carboxy-propenyl) phosphonic acids present intramolecular hydrogen bonds, existing in at least four conformer with internal hydrogen bonds. These last compounds were more rigorously studied at RHF/3-21G^* and RHF/6-31G^** levels. The most stable conformer shows a trans structure for the C=C/P=O angle, with an intramolecular hydrogen bond located between the hydroxylic hydrogen of phosphonic group and the carbonyl oxygen of carboxylic moiety. A secondary conformer is found with a double intramolecular hydrogen bond between two hydroxylic hydrogens of the phosphonic moiety and the oxygen of carboxylic bond. Another secondary conformer appears with an intramolecular hydrogen bond between the oxygen of the phosphoryl bond and the hydroxylic hydrogen of the carboxylic group. A study of the topology of charge densities is carried out. This analysis reveals bonds with an ionic participation. A very weak π conjugation, variable with the conformers, is found in the C=C/P=O system, as well as a strongly polarized P=O partial triple bond. The intramolecular hydrogen bonds give rise to cyclic structures.     

Gas phase and water solution conformational analysis of the herbicide diuron (DCMU): an ab initio study

In the present work, the conformational equilibrium for the herbicide diuron (DCMU) has been investigated using high level ab initio calculations. The solvent effect was included through two different continuum models: (1) the real cavity IPCM method and (2) the standard dipole Onsager model SCRF. The effect due to solute-solvent hydrogen-bond interactions was analyzed considering a hybrid discreet-continuum model. At the Hartree-Fock level, the gas phase results showed that only the trans forms (A and B) are present in the equilibrium mixture, with the relative concentrations found to be 33% (A) and 67% (B) (HF/6-311+G**//6-31G**). When the electronic correlation effect is included (MP2/6-31G*//HF/6-31G*), a relative stabilization of the cis forms was observed, with the conformational distribution calculated as 38% (A), 50% (B), 6% (C) and 6% (D). The trans conformations were found to be completely planar, which has been considered to be a prerequisite for the herbicide binding. In water solution, the trans conformation A should be the most abundant conformer, the IPCM and SCRF values being ca. 100% and ca. 85% respectively. The IPCM calculations with the isodensity level set to 0.0005 present a conformational distribution close to that obtained from the hybrid model [92% (A) and 8% (B)], which has been considered our best solvent approach. Regarding the biological action of urea-type herbicides, the results presented here are important, because some QSAR studies have suggested that the partition coefficient is related to the herbicide activity, so the conformational equilibrium may play a role in the biological action. Received: 23 February 1998 / Accepted: 28 May 1998 / Published online: 19 August 1998     

Vibrational analysis of a rate-slowing conformational kinetic isotope effect

An enthalpy-entropy approach to analyzing a rate-slowing conformational kinetic isotope effect (CKIE) in a deuterated doubly-bridged biaryl system is described. The computed isotope effect (k_H/k_D?=?1.075, 368?K) agrees well with the measured value (k_H/k_D?=?1.06, 368?K). The rate-slowing (normal isotope effect) nature of the computed CKIE is shown to originate from a vibrational entropy contribution defined by the twenty lowest frequency normal modes in the ground state and transition state structures. This normal entropy contribution is offset by an inverse vibrational enthalpy contribution, which also arises from the twenty lowest frequency normal modes. Zero point vibrational energy contributions are found to be relatively small when all normal modes are considered. Analysis of the H_ZPE, H_vib, and S_vib energy terms arising from the low frequency vibrational modes reveals their signs and magnitudes are determined by larger vibrational energy differences in the labeled and unlabeled ground state structures.     

The influence of end-capping on the enantioselectivity of a chiral phase

Summary The effect of end-capping chiral stationary phases (CSP's) derived fromN-(2-naphthyl)alanine undecyl ester has been examined using either trimethylchlorosilane (TMCS), hexamethyldisilazane (HMDS), or bis(trimethylsilyl) trifluoroacetamide (BSTFA) as end-capping reagents. The separation factor () and capacity factor (k) of the enantiomers ofN-(3,5-dinitrobenzoyl)leucine octadecyl amide andN-(3,5-dinitrobenzoyl)alanine butyl ester were evaluated on three columns all packed with material from the same batch of stationary phase. These columns were essentially identical before, but not after end-capping with the above reagents. TMCS and HMDS were found to be superior to BSTFA, which appears to cause a significant loss of bonded phase from the silica surface. It seems that residual silanols affect the retention either by interacting with the analyte or by interacting with strands of stationary phase. End-capping usually increases enantioselectivity, sometimes by decreasing k for the first enantiomer and increasing k for the second enantiomer. The enhancement in enantioselectivity is greatest in relatively nonpolar mobile phases and occurs to a greater extent for phases having incomplete surface coverages.     

Ab initio conformational maps in the gas phase and aqueous solution for a prototype of the glycosidic linkage

Ab initio conformational maps for methoxyethoxymethane (MEM) in both the gas phase and aqueous solution have been constructed using two different approaches. The results obtained allow us to conclude that a rigid conformational map is able to predict the regions of the minima, in the potential energy surface of MEM, in full agreement with those found in the relaxed conformational map, in both phases studied. This is a good indication that ab initio rigid conformational maps may be reliably used to sort the stablest conformers of disaccharides in aqueous solution. Besides that, in the MEM case, the solvation effects do not give rise to any new local minimum in its potential energy surface, but just change the relative energies of the stablest conformers found in the gas phase. This may be an indication that even in aqueous solution the anomeric effect is still the determinant effect defining the conformation of the molecule.Proceedings of the 11th International Congress of Quantum Chemistry satellite meeting in honor of Jean-Louis Rivail     

Infrared and Raman spectra, conformational stability, ab initio calculations of structure, and vibrational assignment of 2-hexyne

The infrared spectra (3500–50 cm⁻¹) of the gas and solid and the Raman spectra (3500–50 cm⁻¹) of the liquid and solid have been recorded for 2-hexyne, CH₃–CC–CH₂CH₂CH₃. Variable temperature studies of the infrared spectrum (3500–400 cm⁻¹) of 2-hexyne dissolved in liquid krypton have also been recorded. Utilizing four anti/gauche conformer pairs, the anti(trans) conformer is found to be the lower energy form with an enthalpy difference of 74±8 cm⁻¹ (0.88±0.10 kJ/mol) determined from krypton solutions over the temperature range −105 to −150 °C. At room temperature it is estimated that there is 42% of the anti conformer present. Equilibrium geometries and energies of the two conformers have been determined by ab initio (HF and MP2) and hybrid DFT (B3LYP) methods using a number of basis sets. Only the HF and DFT methods predict the anti conformer as the more stable form as found experimentally. A vibrational assignment is proposed based on the force constants, relative intensities, depolarization ratios from the ab initio and DFT calculations and on rotational band contours obtained using the calculated equilibrium geometries. From calculated energies it is shown that the CH₃ group exhibits almost completely free rotation which is in agreement with the observation of sub-band structure for the degenerate methyl vibrations from which values of the Coriolis coupling constants, ζ, have been determined. The results are compared to similar properties of some corresponding molecules.     

Comparison of conformer distributions in the crystalline state with conformational energies calculated by ab initio techniques

Summary The conformational preferences of 12 molecular substructures in the crystalline state have been determined and compared with those predicted for relevant model compounds by ab initio molecular orbital calculations. Least-squares regression shows that there is a statistically significant correlation between the crystal-structure conformer distributions and the calculated potential-energy differences, even though the calculations relate to a gas-phase environment. Torsion angles associated with high strain energy (>1 kcal mol^-1) appear to be very unusual in crystal structures and, in general, high-energy conformers are underrepresented in crystal structures compared with a gas-phase, room-temperature Boltzmann distribution. It is concluded that crystal packing effects rarely have a strong systematic effect on molecular conformations. Therefore, the conformational distribution of a molecular substructure in a series of related crystal structures is likely to be a good guide to the corresponding gas-phase potential energy surface.     

Infrared, and Raman spectra, conformational stability, normal coordinate analysis, ab initio calculations, and vibrational assignment of 1-chlorosilacyclobutane

The infrared spectra (3500–40 cm⁻¹) of gaseous and solid and the Raman spectra (3500–30 cm⁻¹) of liquid and solid 1-chlorosilacyclobutane, c-C₃H₆SiClH, have been obtained. Both the axial and equatorial conformers with respect to the chlorine atom have been identified in the fluid phases. Variable temperature (−105 to −150°C) studies of the infrared spectra of the sample dissolved in liquid krypton have been carried out. From these data, the enthalpy difference has been determined to be 211±17 cm⁻¹ (2.53±0.21 kJ/mol), with the equatorial conformer being the more stable form and the only conformer remaining in the annealed solid. At ambient temperatures, approximately 26% of the axial conformers are present in the vapor phase. A complete vibrational assignment is proposed for the equatorial conformer, and many of the fundamentals of the axial conformers have also been identified. The vibrational assignments are supported by normal coordinate calculations utilizing ab initio force constants. Complete equilibrium geometries, conformational stabilities, harmonic force fields, infrared intensities, Raman activities, and depolarization ratios have been determined for both rotamers by ab initio calculations employing the 6-31G(d) basis set at the levels of restricted Hartree–Fock (RHF) and/or Moller–Plesset (MP) to second order. Structural parameters have also been obtained using MP2/6-311+G(d,p) ab initio calculations. The r₀ parameters for both conformers are obtained from a combination of the ab initio predicted values and the twelve previously reported microwave rotational constants. The results are discussed and compared to those obtained for some similar molecules.     

A theoretical investigation into chiral phosphoric acid-catalyzed asymmetric Friedel-Crafts reactions of nitroolefins and 4,7-dihydroindoles: reactivity and enantioselectivity

This article mainly focused on high level Density Functional Theory (DFT) studies on the chiral phosphoric acid-catalyzed Friedel-Crafts reactions between 4,7-dihydroindoles and nitroolefins. Firstly, the reactivities of 4,7-dihydroindole and indole in the chiral phosphoric acid-catalyzed Friedel-Crafts reactions with nitroolefin have been compared. The higher reactivity of 4,7-dihydroindole could be attributed to its higher HOMO energy as well as its more suitable trajectory to attack the nitroolefin in the transition state. Secondly, the origin of the enantioselectivity of the chiral phosphoric acid-catalyzed Friedel-Crafts reaction of 4,7-dihydroindole with nitroolefin has been studied using complete models on PBE1PBE/[6-311+G(d,p), 6-31G(d,p)] level. When (S)-1b was used as the catalyst, the enantioselectivity of the reaction is entirely controlled by the steric effect between the catalyst and the substrate. Whereas for catalyst (S)-1c the enantioselectivity is determined by the solvent effect.     

Infrared and Raman spectra, conformational stability, normal coordinate analysis, ab initio calculations, and vibrational assignment of 1-fluoro-1-methylsilacyclobutane

The infrared (3500–40 cm⁻¹) spectra of gaseous and solid 1-fluoro-1-methylsilacyclobutane, c-C₃H₆SiF(CH₃), have been recorded. Additionally, the Raman spectrum (3500–30 cm⁻¹) of the liquid has been recorded and quantitative depolarization values have been obtained. Both the axial and equatorial (with respect to the methyl group) conformers have been identified in the fluid phases. Variable temperature (−55–−100°C) studies of the infrared spectra of the sample dissolved in liquid xenon have been carried out. From these data, the enthalpy difference has been determined to be 267±10 cm⁻¹ (3.19±0.12 kJ mol⁻¹), with the axial conformer being the more stable form and the only conformer remaining in the polycrystalline solid. A complete vibrational assignment is proposed for the axial conformer and many of the fundamentals for the equatorial conformer have also been identified. The vibrational assignments are supported by normal coordinate calculations utilizing ab initio force constants. Complete equilibrium geometries have been determined for both rotamers by ab initio calculations employing the 6-31G* and 6-311++G** basis sets at the levels of restricted Hartree–Fock (RHF) and/or Moller–Plesset (MP) to second order. The results are discussed and compared to those obtained for some similar molecules.     

Theoretical investigation into the influence of conformational equilibria on the water-exchange process in magnetic resonance imaging contrast agents

The conformational behavior in aqueous solution of four complexes of the Eu(III) ion with bis (R–amide) derivatives (R=H, methyl, ethyl, butyl) of diethylentriamine pentacetate ligands has been characterized at the ab initio level to rationalize the experimentally observed influence of alkyl substituents on the rate of the exchange process of the water molecule coordinated to the ion with the bulk water. Calculations were performed in vacuo and for aqueous solution, the latter by using the polarizable continuum model. Geometry optimizations provide, for each system, four isomers as stable conformations, all presenting a distorted tricapped trigonal prism coordination geometry around the ion. No significant influence of the alkyl substitution on the coordination geometry, nor on the europium–water distance, was observed. Moreover, increasing the length of the alkyl chain had no significant effect on the relative isomer population in solution. Thus, these results lead us to suppose that other effects, like those deriving from lateral chain folding in solution, should be considered to explain the increased rate of the water-exchange process with alkyl chain lengthening.Contribution to the Jacopo Tomasi Honorary Issue     

Ab initio studies of conformational properties of dimethyl diphosphate dianion and its complex with magnesium

The conformational properties of the diphosphate linkage have been studied with ab initio methods using the dimethyl diphosphate dianion (1) and magnesium dimethyl diphosphate (2) as models. The ab initio energy and geometry of the conformers around the P-O bonds have been determined at the self-consistent-field (SCF) using the 6-31G* and the tzp basis sets; whereas, the 6-31G* basis set alone has been used for 2. In addition, the adiabatic connection method (ACM) of density functional theory (DFT) using the dzvp basis set has been employed for 1. The optimization of all possible staggered conformers assumed for the four P-O bonds, led to nine minima for 1. In agreement with the general anomeric effect, the sc conformation about the P-O bonds is clearly preferred over the ap one. Vibrational frequencies were calculated at the SCF level using the 6-31G* basis set and used to evaluate zero-point energies, thermal energies, and entropies for all minima of 1. The effect of zero-point energies and thermal energies is quite small. However, the effect of entropies, mainly resulting from a multiplicity contribution, changes the stability of the conformers. For each minimum of 1, up to six different arrangements of the Mg²⁺ were used to determine minima of 2. This procedure led to 21 distinct minima. The presence of the magnesium counter-ion appeared to completely change the structure and relative energy of the conformers. The preferred structures of the complex exhibit the (sc, ap) orientation around the two central P–O bonds and an arrangement in which the magnesium cation is coordinated by three phosphoryl oxygen atoms. The results of this work clearly demonstrate that interactions with the metal counter-ion can induce conformational changes in the overall 3D-shape adopted by molecules containing diphosphate linkages. The PM3 and MNDO quantum semi-empirical methods and molecular mechanics methods using the CVFF force field were tested and large differences in the minimum structures, as well as in the conformational energies between these and ab initio methods, are discussed.     

Ab initio studies of the conformers and conformational distribution of the gaseous hydroxyamino acid threonine

Systematic and extensive conformational search has been performed to characterize the gas-phase threonine structures. A total of 1296 unique trial structures were generated by allowing for all combinations of internal single-bond rotamers. All the trial structures were optimized at the B3LYP/6-311G* level of the theory and then subjected to further optimization at the B3LYP/6-311++G** level. A total of 71 conformers were found and their rotational constants, dipole moments, zero-point vibrational energies, harmonic frequencies and vertical ionization energies of all the conformers were determined. Single-point energies were also calculated at the MP2/6-311G(2df,p) and B3LYP/6-311G(2df,p) levels. Characteristic H-bonding types were classified and listed for all the conformers. The conformational distributions of gaseous threonine at various temperatures were calculated.     

Application of chiral dipyridylmethane ligands in the enantioselective palladium-catalyzed allylic alkylation

New chiral C₁-symmetric dipyridylmethane ligands have been prepared from naturally occurring monoterpenes according to a method based on two consecutive constructions of the pyridine rings. These ligands have been assessed in the enantioselective palladium catalyzed allylic alkylation of 1,3-diphenylprop-2-enyl acetate with dimethylmalonate. Enantioselectivity up to 68% ee has been obtained.     

On the vibrational spectra and conformational stability of 1,1-dimethylhydrazine from temperature dependent FT-IR spectra of krypton solutions and ab initio calculations

Variable temperature (−105 to −150 °C) studies of the infrared spectra (3500–400 cm⁻¹) of 1,1-dimethylhydrazine, (CH₃)₂NNH₂, in liquid krypton have been carried out. No convincing spectral evidence could be found for the trans conformer which is expected to be at least 600 cm⁻¹ less stable than the gauche form. The structural parameters, dipole moments, conformational stability, vibrational frequencies, and infrared and Raman intensities have been predicted from MP2/6-31G(d) ab initio calculations. The predicted infrared and Raman spectra are compared to the experimental ones. The adjusted r₀ parameters from MP2/6-311+G(d,p) calculations are compared to those reported from an electron diffraction study. The energy differences between the gauche and trans conformers have been obtained from MP2 ab initio calculations as well as from density functional theory by the B3LYP method calculations from a variety of basis sets. All of these calculations indicate an energy difference of 650–900 cm⁻¹ with the B3LYP calculations predicted the larger values. The potential function governing the conformational interchange has been predicting from both types of calculations and comparisons have been made. The barrier to internal rotation by the independent rotor model of the inner methyl group is predicted to have a value of 1812 cm⁻¹ and that of the outer one of 1662 cm⁻¹ from ab initio MP2/6-31G(d) calculations. These values agree well with the experimentally determined values of 1852±16 and 1558±12 cm⁻¹, respectively, from a fit of the torsional transitions with the coupled rotor model. For the coupled rotor model the predicted V′₃₃ (sin 3τ₀ sin 3τ₁ term) value which ranged from 190 to 232 cm⁻¹ is in reasonable agreement with the experimental value of 268±3 cm⁻¹ but the predicted V₃₃ (cos 3τ₀ cos 3τ₁ term) value of −73 to −139 cm⁻¹ is 25% smaller and of the opposite sign of the experimental value of 333±22 cm⁻¹. These theoretical and spectroscopy results are compared to similar quantities of some corresponding molecules.     

A Conformational Study on Silacyclohexane. Comparison of ab initio (HF,MP2), DFT,and Molecular Mechanics Calculations. Conformational Energy Surface of Silacyclohexane

The structures and relative energies for the basic conformations of silacyclohexane 1 have been calculated using HF, RI‐MP2, RI‐DFT and MM3 methods. All methods predict the chair form to be the dominant conformation and all of them predict structures which are in good agreement with experimental data. The conformational energy surface of 1 has been calculated using MM3. It is found that there are two symmetric lowest energy pathways for the chair‐to‐chair inversion. Each of them consists of two sofa‐like transition states, two twist forms with C₁ symmetry (twist‐C₁), two boat forms with Si in a gunnel position (C₁ symmetry), and one twist form with C₂ symmetry (twist‐C₂). All methods calculate the relative energy to increase in the order chair < twist‐C₂ < twist‐C₁ < boat. At the MP2 level of theory and using TZVP and TZVPP (Si atoms) basis sets the relative energies are calculated to be 3.76, 4.80, and 5.47 kcal mol^–1 for the twist‐C₂, twist‐C₁, and boat conformations, respectively. The energy barrier from the chair to the twisted conformations of 1 is found to be 6.6 and 5.7 kcal mol^–1 from MM3 and RI‐DFT calculations, respectively. The boat form with Si at the prow (C_s symmetry) does not correspond to a local minimum nor a saddle point on the MM3 energy surface, whereas a RI‐DFT optimization under C_s symmetry constraint resulted in a local minimum. In both cases its energy is above that of the chair‐to‐twist‐C₁ transition state, however, and it is clearly not a part of the chair‐to‐chair inversion.     

Molecular structure of free canonical 2′-deoxyribonucleosides: a density functional study

The molecular structure of isolated canonical 2′-deoxyrinobucleosides was calculated using the density functional theory. It was demonstrated that the geometry of the base unit (BU) is almost unchanged compared to free nucleobases. Only slight out-of-plane deformation of the pyrimidine ring in deoxy-cytidine is observed. The conformation of the furanose ring strongly depends on the nature and orientation of the nucleobase. All nucleosides possess different conformations of this ring. Significant influence of the steric repulsion between the nucleobase and the sugar unit (SU) on puckering of the furanose ring and variation of the C–O and glycosyl bond lengths was demonstrated. The C(3′)-endo conformer of the furanose ring is more stable at the anti-orientation BU with respect to SU. An opposite trend is observed for the syn-orientation which is additionally stabilized by an intramolecular hydrogen bond with participation of the C(5′)OH group.     

Raman and infrared spectra, conformational stability, barriers to internal rotation, vibrational assignment, and ab initio calculations of 3-methyl-3-butenenitrile

The Raman (3500-30 cm⁻¹) spectra of liquid and solid and the infrared (3500-40 cm⁻¹) spectra of gaseous and solid 3-methyl-3-butenenitrile, CH₂C(CH₃)CH₂CN, have been recorded. Both cis and gauche conformers have been identified in the fluid phases but only the cis form remains in the solid. Variable temperature (−55 to −100 °C) studies of the infrared spectra of the sample dissolved in liquid xenon have been carried out. From these data, the enthalpy difference has been determined to be 163±16 cm⁻¹ (1.20±0.19 kJ mol⁻¹), with the cis conformer the more stable rotamer. It is estimated that there is 48±2% of the gauche conformer present at  25°C. A complete vibrational assignment is proposed for the cis conformer based on infrared band contours, relative intensities, depolarization ratios and group frequencies. Several of the fundamentals for the gauche conformer have also been identified. The vibrational assignments are supported by normal coordinate calculations utilizing ab initio force constants. Complete equilibrium geometries have been obtained for both rotamers by ab initio calculations employing the 6-31G(d), 6-311G(d,p), 6-311+G(d,p) and 6-311+G(2d,2p) basis sets at the levels of restricted Hartree-Fock (HF) and/or Møller-Plesset perturbation theory to the second order (MP2). Only with the 6-311G(2d,2p) and 6-311G(2df,2pd) basis sets with or without diffuse functions is the cis conformer predicted to be more stable than the gauche form. The potential energy terms for the conformational interchange have been obtained at the MP2(full)/6-311+G(2d,2p) level, and compared to those obtained from the experimental data. The results are discussed and compared to the corresponding quantities obtained for some similar molecules.     

Conformational Analysis of 2,2′-bifuran: Correlated High-level Ab initio and DFT Results

The torsional potential for inter-ring rotation in 2,2′-bifuran has been systematically tackled using highly accurate ab initio calculations as well as cost-effective DFT methods. The successful convergence of the ab initio results allowed to confirm the presence of a shallow gauche minimum in the torsional potential curve. The standard DFT methods failed to capture such a tiny energy barrier but, interestingly, the results could be remarkably improved by a mixture of wavefunction and DFT energies in a multi-coefficient fashion; thus, accurate DFT-based and ab initio reference data also become available. Since the experimental evaluation of torsional potentials faces quantitative problems, the outcome of high-level theoretical calculations is expected to be reliably used in further investigation on structure and conformational distribution of this system.     

新型多齿手性氨基醇的合成和结构表征

以L-苯甘氨酸和N,N-二甲基苯胺为原料,合成了2种新型多齿手性氨基醇S-1,1,2-三苯基-2-[(2-二甲氨基-5-溴-1-苯基)甲氨基]-1-乙醇和S-2-苯基-2-[(2-二甲氨基-5-溴-1-苯基)甲氨基]-1-乙醇;利用红外光谱仪、核磁共振谱仪及质谱仪表征了合成产物的结构.