Conformational study of jet-cooled L-phenylglycine

We investigated the conformational structures of L-phenylglycine in the gas phase by photoionization and double resonance spectroscopy techniques as well as high-level ab initio calculations. The UV-UV and IR-UV double resonance spectroscopy suggested that there exists only one conformer that has a free OH band for the carboxyl group. Rotational contour analysis combined with ab initio calculation indicated that the conformer we detected by resonant two-photon ionization was not one of those found by Sanz et al. in their microwave spectroscopic study [Chem. Eur. J. 12, 2564 (2006)]. Different methods of vaporization along with different expansion and cooling conditions and different detection methods are believed to be the culprit for such intriguing discrepancy. The identical hydrogen bonding structure of our phenylglycine conformer with the most abundant conformer of glycine found in helium droplets and their nearly identical OH frequencies suggest that the skeletal structure of glycine is not significantly altered by phenyl substitution.     

Theoretical study of the reaction mechanism of proton transfer in glycinamide

To investigate the tautomerism of glycinamide that is induced by proton transfer, we present detailed theoretical studies on the reaction mechanism of both the isolated gas phase and H₂O‐assisted proton transfer process of glycinamide, using density functional theory calculations by means of the B3LYP hybrid functional. Twenty‐six geometries, including 10 significant transition states, were optimized, and these geometrical parameters are discussed in detail. The relative order of the activation energy for hydrogen atom transfer of all the conformers has been systematically explored in this essay. For the amido hydrogen atom transfer process, the relative order of the activation energy is: IV < II < III < I, while in the carbonic hydrogen atom transfer process, the relative order is IV > II > III > I. Meanwhile, the most favorable structure for both the amido hydrogen atom transfer and the carbonic hydrogen atom transfer has been found. The involvement of the water molecule not only can stabilize the transition states and the ground states, but can also reduce the activation energy greatly. The superior catalytic effect of H₂O has been discussed in detail. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Polymorphs of the antiviral drug ganciclovir

Ganciclovir (GCV; systematic name: 2‐amino‐9‐{[(1,3‐dihydroxypropan‐2‐yl)oxy]methyl}‐6,9‐dihydro‐1H‐purin‐6‐one), C₉H₁₃N₅O₄, an antiviral drug for treating cytomegalovirus infections, has two known polymorphs (Forms I and II), but only the structure of the metastable Form II has been reported [Kawamura & Hirayama (2009). X‐ray Struct. Anal. Online , 25 , 51–52]. We describe a successful preparation of GCV Form I and its crystal structure. GCV is an achiral molecule in the sense that its individual conformers, which are generally chiral objects, undergo fast interconversion in the liquid state and cannot be isolated. In the crystalline state, GCV exists as two inversion‐related conformers in Form I and as a single chiral conformer in Form II. This situation is similar to that observed for glycine, also an achiral molecule, whose α‐polymorph contains two inversion‐related conformers, while the γ‐polymorph contains a single conformer that is chiral. The hydrogen bonds are exclusively intermolecular in Form I, but both inter‐ and intramolecular in Form II, which accounts for the different molecular conformations in the two polymorphs.     

The Shape of the Simplest Non-proteinogenic Amino Acid α-Aminoisobutyric Acid (Aib)

The simplest non-proteinogenic amino acid α-aminoisobutyric acid (Aib), an analogue of glycine and alanine, has been vaporized by laser ablation and probed by high-resolution Fourier transform microwave spectroscopic techniques. Comparison of the experimental rotational and ¹⁴N nuclear quadrupole constants with that predicted ab initio has allowed the identification of three conformers of Aib exhibiting three types of hydrogen-bond interactions I (NH⋅⋅⋅O=C, cis-COOH), II (OH⋅⋅⋅N, trans-COOH), and III (N−H⋅⋅⋅O−H, cis-COOH) within the amino acid backbone. The observation of conformer III, not detected previously for related proteinogenic amino acids with a nonpolar side chain in a supersonic expansion, indicates that the presence of the methyl groups should restrict the conformational relaxation from conformer Aib-III to Aib-I. For conformer Aib-II, the rotational spectra of the ¹³C isotopomers reveal a tunneling motion arising from the two equivalent methyl groups in the molecule. The observation of a single spectrum at the midpoint between those predicted for the two ¹³C of the methyl groups has been explained by considering a double-minimum potential function with a low-energy interconversion barrier for a large amplitude internal motion. This singular fact has been corroborated by the anomalous centrifugal distortion effects determined in conformer Aib-II.     

Potential surface crossing in classical trajectory calculations: Application to collision-induced predissociation in electronically excited I2

Classical trajectory calculations of inelastic collisions involving electronically excited iodine molecules with rare-gas atoms have been extended to include the effects of potential surface crossing. Adiabatic potential surfaces were generated assuming an r^?6 dependence of the energy splitting on the molecule—atom distance. Nonadiabatic transition probabilities were calculated using the Landau—Zener formula. The reduced-mass and final-state dependence of the calculated induced predissociation probabilities are in good agreement both with experiment and with earlier optical-model calculations.     

Semiclassical theory for atom—surface scattering: Calculations on Ar + W(110)

By using a semiclassical approach we derive an effective potential which governs the motion of the colliding atom/molecule. The effective potential depends on the surface temperature and the phonon excitation during the collision. Numerical results on Ar + W(110) are discussed and compared with experimental data.     

A microwave and quantum chemical study of the conformational properties and intramolecular hydrogen bonding of 1-fluorocyclopropanecarboxylic acid

The structural and conformational properties of 1-fluorocyclopropanecarboxylic acid have been explored by microwave spectroscopy and a series of ab initio (MP2/6-311++G(d,p) level), density functional theory (B3LYP/aug-cc-pVTZ level), and G3 quantum chemical calculations. Four "stable" conformers, denoted conformers I-IV, were found in the quantum chemical calculations, three of which (conformers I -III) were predicted to be low-energy forms. Conformer I was in all the quantum chemical calculations predicted to have the lowest energy, conformer III to have the second lowest energy, and conformer II to have the third lowest energy. Conformers II and III were calculated to have relatively large dipole moments, while conformer I was predicted to have a small dipole moment. The microwave spectrum was investigated in the 18-62 GHz spectral range. The microwave spectra of conformers II and III were assigned. Conformer I was not assigned presumably because its dipole moment is comparatively small. Conformer II is stabilized by an intramolecular hydrogen bond formed between the fluorine atom and the hydrogen atom of the carboxylic acid group. Conformer III has a synperiplanar orientation for the F-C-C=O and H-O-C=O chains of atoms. Its dipole moment is: mua = 3.4(10), mub = 10.1(13), and muc = 0.0 (assumed) and mu(tot) = 10.6(14) x 10(-30) C m [3.2(4) D]. Several vibrationally excited states of the lowest torsional mode of each of II and III were also assigned. The hydrogen-bonded conformer II was found to be 2.7(2) kJ/mol less stable than III by relative intensity measurements. Absolute intensity measurements were used to show that the unassigned conformer I is the most abundant form present at a concentration of roughly 65% at room temperature. Conformer I was estimated to be ca. 5.0 kJ/mol more stable than the hydrogen-bonded rotamer (conformer II) and ca. 2.3 kJ/mol more stable than conformer III. The best agreement with the theoretical calculations is found in the MP2 calculations, which predict conformer I to be 5.1 kJ/mol more stable than III and 1.7 kJ/mol more stable than II.     

Vibrational spectrum of glycine molecule

The infrared and Raman spectra of glycine molecule has been studied in spectral region 400-4000 cm(-1) in solid form as well as in water. The vibrational frequencies for the fundamental modes of the glycine in neutral and its zwitterionic form have also been calculated using AM1 semiempirical method as well as ab initio method with minimal basis set. The reliability of the minimal basis set and AM1 method with higher basis sets, for IR spectra of the neutral glycine conformers were examined. We find that the 6-21G basis set calculation yields structural parameters, rotational constant and dipole moment of glycine conformers, which are very similar to those obtained from extended basis set calculation as well as experimental values. IR frequencies for glycine conformer I are also calculated in water using SCRF=PCM model and compared with experimental values. A comparison between calculated frequencies for neutral glycine, and its zwitterionic form with observed IR and Raman bands have been made. The total energies for gas phase glycine and its zwitterionic form along with those of hydrated forms were also calculated. It is found from the calculations that in the gas phase neutral glycine is more stable as compared to its zwitterionic form.     

A theoretical investigation of the low energy conformers of the isomers glycine and methylcarbamic acid and their role in the interstellar medium

We have theoretically investigated the low energy conformers of neutral glycine (NH(2)CH(2)COOH) and its isomer methylcarbamic acid (CH(3)NHCOOH) in the gas phase. A total of 16 different levels of the theory, including CCSD(T), MP2 and B3LYP methods with various Pople and Dunning type basis sets with and without polarization and diffuse functions were used. We found eight low energy glycine conformers, where the heavy atoms in three have a planar backbone, and four low energy methylcarbamic acid conformers all with non-planar backbones. Interestingly at all levels of theory, we found that the most stable methylcarbamic acid conformer is significantly lower in energy than the lowest energy glycine conformer. The MP2 level and single point CCSD(T) calculations show the lowest energy methylcarbamic acid conformer to be between 31 to 37 kJ mol(-1) lower in energy than the lowest energy glycine conformer. These calculations suggest that methylcarbamic acid might serve as a precursor to glycine formation in the Interstellar Medium (ISM). We also report the theoretical harmonic vibrational frequencies, infrared intensities, moment of inertia, rotational constants and dipole moments for all of the conformers. In order to understand how glycine or methylcarbamic acid might be formed in the ISM, larger calculations which model glycine or its isomer interacting with several surrounding molecules, such as water, are needed. We demonstrate that B3LYP method should provide a reliable and computationally practical approach to modeling these larger systems.     

Modeling the fragmentation dynamics of ionic clusters inside helium nanodroplets: the case of He100Ne4+

We present simulation results on the effect of a helium nanodroplet environment on the fragmentation dynamics of embedded molecular systems. The helium atoms are treated explicitly, with zero-point effects taken into account through an effective helium-helium interaction potential. The ionized neon tetramer is used as a model molecular system because, like all the small rare-gas clusters, it fragments extensively upon ionization. All the nonadiabatic effects between electronic states of the ionized neon cluster are taken into account. The results reveal a predominance of Ne2+ and HepNe2+ fragments and the absence of bare Ne+ fragments, in agreement with available experimental data. The neutral monomer fragments exhibit a rather wide kinetic energy distribution that can be fitted to the sum of two Boltzmann distributions, one with a low kinetic energy and the other with a higher kinetic energy. This indicates that cooling by helium atom evaporation is more efficient than was believed so far, as suggested by recent experimental results. Purely classical calculations are shown to strongly overestimate the amount of cage effect (cooling), clearly indicating the need to take into account zero-point effects.     

A conformational analysis of the six isomers of oxybispyridine

Conformational energy contour maps of the six isomers of oxybispyridine have been constructed using the ab initio STO-3G molecular orbital method. The calculations (employing a partial rigid rotor) for all six isomers indicate that the minimum energy conformers are not planar and that energy barriers between 70–1000 kJ mol^?1 restrict interconversion to planar structures, thereby preventing conjugation between the p-electrons of the oxygen atom with the π system of the pyridine rings. It is postulated that of the three mechanisms used to explain conformer interconversion about the C? O bond, the disrotatory one-ring flip mechanism is the most appropriate, since the “Morino's” structures are all within 2.5 kJ mol^?1 of the minimum. Furthermore, room temperature accessibility of the “Morino's” structures suggests that the Smiles rearrangement would be possible for suitably substituted derivatives of these isomers.     

Collision-energy-resolved penning ionization electron spectroscopy of HCOOH, CH3COOH, and HCOOCH3 by collision with He*(2(3)S) metastable atoms

Penning ionization of formic acid (HCOOH), acetic acid (CH3COOH), and methyl formate (HCOOCH3) upon collision with metastable He*(2(3)S) atoms was studied by collision-energy/electron-energy-resolved two-dimensional Penning ionization electron spectroscopy (2D-PIES). Anisotropy of interaction between the target molecule and He*(2(3)S) was investigated based on the collision energy dependence of partial ionization cross sections (CEDPICS) obtained from 2D-PIES as well as ab initio molecular orbital calculations for the access of a metastable atom to the target molecule. For the interaction potential calculations, a Li atom was used in place of He*(2(3)S) metastable atom because of its well-known similarity in interaction with targets. The results indicate that in the studied collision energy range the attractive potential localizes around the oxygen atoms and that the potential well at the carbonyl oxygen atom is at least twice as much as that at the hydroxyl oxygen. Moreover we can notice that attractive potential is highly anisotropic. Repulsive interactions can be found around carbon atoms and the methyl group.     

Semi-classical study of the molecular dissociation: H2 + He

Accurate vibrational wave functions and a state-dependent model interaction potential were used in the study, within the framework of a semi-classical theory, of the vibrational excitation and dissociation of the hydrogen molecule in collinear collisions with the helium atom. A molecule initially in the excited state is shown to be very efficient in energy transfer and twice more likely to be further excited than to be de-excited. The change in the population distribution among the vibrational states at the first few collisions was analyzed. It is shown that the population of the first vibrational excited state ψ₁ reaches its maximum after the very first collision and that of ψ₂ after the second. It is also found that at a sufficiently high collision energy, ψ₅ is the most efficient state in dissociation at the second collision while ψ₆ contributes most at the third collision.     

Dissociation of carbonic acid: gas phase energetics and mechanism from ab initio metadynamics simulations

A comprehensive metadynamics study of the energetics, stability, conformational changes, and mechanism of dissociation of gas phase carbonic acid, H2CO3, yields significant new insight into these reactions. The equilibrium geometries, vibrational frequencies, and conformer energies calculated using the density functional theory are in good agreement with the previous theoretical predictions. At 315 K, the cis-cis conformer has a very short life time and transforms easily to the cis-trans conformer through a change in the O=C-O-H dihedral angle. The energy difference between the trans-trans and cis-trans conformers is very small (approximately 1 kcal/mol), but the trans-trans conformer is resistant to dissociation to carbon dioxide and water. The cis-trans conformer has a relatively short path for one of its hydroxyl groups to accept the proton from the other end of the molecule, resulting in a lower activation barrier for dissociation. Comparison of the free and potential energies of dissociation shows that the entropic contribution to the dissociation energy is less than 10%. The potential energy barrier for dissociation of H2CO3 to CO2 and H2O from the metadynamics calculations is 5-6 kcal/mol lower than in previous 0 K studies, possibly due to a combination of a finite temperature and more efficient sampling of the energy landscape in the metadynamics calculations. Gas phase carbonic acid dissociation is triggered by the dehydroxylation of one of the hydroxyl groups, which reorients as it approaches the proton on the other end of the molecule, thus facilitating a favorable H-O-H angle for the formation of a product H2O molecule. The major atomic reorganization of the other part of the molecule is a gradual straightening of the O=C=O bond. The metadynamics results provide a basis for future simulation of the more challenging carbonic acid-water system.     

Theoretical investigation of rotationally inelastic collisions of CH2(ã) with helium

Rotationally inelastic collisions of the CH(2) molecule in its a?(1)A(1) electronic state have been investigated. We have determined a potential energy surface (PES) for the interaction of rigid CH(2)(a?), frozen at its equilibrium geometry, with a helium atom, using a coupled-cluster method that includes all single and double excitations, as well as perturbative contributions of connected triple excitations [RSSCD(T)]. The PES is quite anisotropic, due to lack of electron density in the unoccupied CH(2) non-bonding orbital perpendicular to the molecular plane. Quantum scattering calculations have been carried out to compute state-to-state rotational energy transfer and elastic depolarization cross sections at collision energies up to 2400 cm(-1). These cross sections were thermally averaged to derive room-temperature rate constants. The total removal and elastic depolarization rate constants for the ortho k(a) = 1 levels agree well with recent experimental measurements by Hall, Sears, and their co-workers. We observe a strong even-odd alternation in the magnitude of the total rate constants which we attribute to the asymmetry splitting of the k(a) = 1 levels.     

The structure and electronic property of the smallest C20‐glycine and Gd‐encapsulated C20‐glycine derivatives with potentially biological activity

Theoretical calculations on interaction of the C₂₀ fullerene (consists solely by pentagons) with the smallest amino acid (glycine) were carried out using density‐functional theory method. The glycine molecule energetically prefers to interact with the Top‐site on the C₂₀ cage via its amino nitrogen (N) active site. The stable ordering of three active sites on glycine molecule is NH₂‐site > O‐site > OH‐site. Moreover, when the Gd atom is encapsulated to the center of C₂₀‐glycine, the cage volume obviously increase ～24.8%; and the endohedral atom induces the generation of two strong bands in the partial density of states spectra, which could cause the effect on optical properties. Additionally, it is also found that the modified C₂₀‐glycine derivative by Gd atom can reduce the thermodynamic and kinetic stabilities. It could be expected that the study may provide a theoretical reference in exploring their intrinsic feature structurally to antitumor activity. © 2012 Wiley Periodicals, Inc.     

Time-dependent density functional theory study of charge transfer in collisions

We study the charge transfer between colliding ions, atoms, or molecules, within time-dependent density functional theory. Two particular cases are presented, the collision between a proton and a Helium atom, and between a gold atom and a butane molecule. In the first case, proton kinetic energies between 16?keV and 1.2?MeV are considered, with impact parameters between 0.31 and 1.9 ?. The partial transfer of charge is monitored with time. The total cross-section is obtained as a function of the proton kinetic energy. In the second case, we analyze one trajectory and discuss spin-dependent charge transfer between the different fragments.     

Matrix isolation and low temperature solid state FTIR spectroscopic study of alpha-furil

Alpha-furil [C(4)H(3)O-C(=O)-C(=O)-C(4)H(3)O] has been isolated in argon and xenon matrices and studied by FTIR spectroscopy, supported by DFT(B3LYP)/6-311++G(d,p) calculations. The obtained spectra were fully assigned and revealed the presence in the matrices of three different conformers, all of them exhibiting skewed conformations around the intercarbonyl bond with the two C(4)H(3)O-C(=O) fragments nearly planar. The three conformers differ in the orientation of the furan rings relative to the carbonyl groups: the most stable conformer, I (C(2) symmetry; O=C-C=O intercarbonyl dihedral equal to 153.1 degrees), has both furan rings orientated in such a way that one of their beta-hydrogen atoms approaches the oxygen atom of the most distant carbonyl group, forming two H-C=C-C-C=O six-membered rings; the second most stable conformer, II (C(1) symmetry; O=C-C=O intercarbonyl dihedral equal to 126.9 degrees ), has one furan ring orientated as in I, while the second furan group is rotated by ca. 180 degrees (resulting in an energetically less favourable H-C=C-C=O five-membered ring); in the third conformer, III (C(2) symmetry; O=C-C=O dihedral equal to 106.2 degrees ), both furan rings assume the latter orientation relative to the dicarbonyl group. The theoretical calculations predicted the two higher energy forms being 5.85 and 6.22 kJ mol(-1) higher in energy than the most stable form, respectively, and energy barriers for conformational interconversion higher than 40 kJ mol(-1). These barriers are high enough to prevent observation of conformational isomerization for the matrix isolated compound. The three possible conformers of alpha-furil were also found to be present in CCl(4) solution, as well as in a low temperature neat amorphous phase of the compound prepared from fast condensation of its vapour onto a suitable 10 K cooled substrate. On the other hand, in agreement with the available X-ray data [S. C. Biswas, S. Ray and A. Podder, Chem. Phys. Lett., 1987, 134, 541], the IR spectra obtained for the neat low temperature crystalline state reveals that, in this phase, alpha-furil exists uniquely in its most stable conformational state, I.     

Controlling Metalation Reaction of Phthalocyanine with Cobalt at Single-Molecule Level on Au(111) Surface

Metalation reaction of metal-free phthalocyanine molecule with Co atom adsorbed on Au(111) surface has been studied in situ at single atom/molecule scale by low-temperature scanning tunneling microscopy (STM) experiment combined with simulations based on density function theory calculations. Through manipulations using STM tip, we showed a controlled manner to have a single metal-free phthalocyanine molecule react with a Co atom to form Co phthalocyanine molecule. In this reaction process, an intermediate state originating from $\pi$-d interaction between the metal-free phthalocyanine molecule and Co atom has been identified. Moreover, we also revealed that the redox reaction represented as bond breaking and bond forming relative to the Co and pyrrolic N atoms, not pyrrolic H atoms, is a key process for dehydrogenation and metalation reaction. Our DFT calculations provided theoretical supporting for the above conclusions, and further understanding of the related mechanisms.