Vibrational spectroscopy and DFT calculations of di-amino acid cyclic peptides. Part I: cyclo(Gly-Gly), cyclo(L-Ala-L-Ala) and cyclo(L-Ala-Gly) in the solid state and in aqueous solution

Investigations of the vibrational spectra of cyclo(Gly-Gly), cyclo(L-Ala-L-Ala) and cyclo(t-Ala-Gly) are reported. Raman scattering and Fourier transform infrared (FTIR) spectra of solid-state and aqueous protonated samples, as well as their corresponding N-deuterated isotopomers, have been examined. In addition, density functional theory (DFT) (B3-LYP/cc-pVDZ) calculations of molecular structures and their associated vibrational modes were carried out. In each case, the calculated structures of lowest energy for the isolated gas-phase molecules have boat conformations. Assignments have been made for the observed Raman and FTIR vibrational bands of the cyclic di-amino acid peptides (CDAPs) examined. Raman polarization studies of aqueous phase samples are consistent with C-2 and C-1 symmetries for the six-membered rings of cyclo(L-Ala-L-Ala) and cydo(L-Ala-Gly), respectively. There is a good correlation between experimental and calculated vibrational bands for the three CDAPs. These data are in keeping with boat conformations for cydo(L-Ala-L-Ala) and cyclo(L-Ala-Gly) molecules, predicted by the ab initio calculations, in both the solid and aqueous solution states. However, Raman spectroscopic results might infer that cyclo(L-AlaGly) deviates only slightly from planarity in the solid state. The potential energy distributions of the amide I and II modes of a cis-peptide linkage are shown to be significantly different from those of the trans-peptides. For example, deuterium shifts have shown that the cis-amide I vibrations found in cyclo(Gly-Gly), cyclo(L-Ala-L-Ala), and cyclo(L-Ala-Gly) have larger N-H contributions compared to their trans-amide counterparts. Compared to trans-amide II vibrations, cis-amide II vibrations show a considerable decrease in N-H character.

IR/Raman spectroscopy and DFT calculations of cyclic di-amino acid peptides. Part III: comparison of solid state and solution structures of cyclo(L-Ser-L-Ser)

B3-LYP/cc-pVDZ calculations of the gas-phase structure and vibrational spectra of the isolated molecule cyclo(L-Ser-L-Ser), a cyclic di-amino acid peptide (CDAP), were carried out by assuming C-2 symmetry. It is predicted that the minimum-energy structure is a boat conformation for the diketopiperazine (DKP) ring with both L-Beryl side chains being folded slightly above the ring. An additional structure of higher energy (15.16 kJ mol(-1)) has been calculated for a DKP ring with a planar geometry, although in this case two fundamental vibrations have been calculated with imaginary wavenumbers. The reported X-ray crystallographic structure of cyclo(L-Ser-L-Ser), shows that the DKP ring displays a near-planar conformation, with both the two L-Beryl side chains being folded above the ring. It is hypothesized that the crystal packing forces constrain the DKP ring in a planar conformation and it is probable that the lower energy boat conformation may prevail in the aqueous environment. Raman scattering and Fourier-transform infrared (FT-IR) spectra of solid state and aqueous solution samples of cyclo(L-Ser-L-Ser) are reported and discussed. Vibrational band assignments have been made on the basis of comparisons with the calculated vibrational spectra and band wavenumber shifts upon deuteration of labile protons. The experimental Raman and IR results for solid-state samples show characteristic amide I vibrations which are split (Raman:1661 and 1687 cm(-1), IR:1666 and 1680 cm(-1)), possibly due to interactions between molecules in a crystallographic unit cell. The cis amide I band is differentiated by its deuterium shift of ~ 30 cm(-1), which is larger than that previously reported for trans amide I deuterium shifts. A cis amide II mode has been assigned to a Raman band located at 1520 cm(-1). The occurrence of this cis amide II mode at a wavenumber above 1500 cm(-1) concurs with results of previously examined CDAP molecules with low molecular weight substituents on the C-alpha atoms, and is also indicative of a relatively unstrained DKP ring.

Vibrational spectra and crystal structure of the di-amino acid peptide cyclo(L-Met-L-Met): Comparison of experimental data and DFT calculations

Experimental Raman and FT-IR spectra of solid-state non-deuterated and N-deuterated samples of cyclo(L-Met-L-Met) are reported and discussed. The Raman and FT-IR results show characteristic amide I vibrations (Raman: 1649 cm-1, infrared: 1675 cm-1) for molecules exhibiting a cis amide conformation. A Raman band, assigned to the cis amide II vibrational mode, is observed at sim1493 cm-1 but no IR band is observed in this region. Cyclo(L-Met-L-Met) crystallises in the triclinic space group P1 with one molecule per unit cell. The overall shape of the diketopiperazine (DKP) ring displays a (slightly distorted) boat conformation. The crystal packing employs two strong hydrogen bonds, which traverse the entire crystal via translational repeats. B3-LYP/cc-pVDZ calculations of the structure of the molecule predict a boat conformation for the DKP ring, in agreement with the experimentally determined X-ray structure. Copyright © 2009 John Wiley & Sons, Ltd.

DFT calculations of the structures and vibrational spectra of the [Fe(bpy)(3)](2+) and [Ru(bpy)(3)](2+) complexes

Structures of the [M(bpy)(3)](2+) complexes (M = Fe and Ru) have been calculated at the B3-LYP/DZVP level. IR and Raman spectra were calculated using the optimised geometries, employing a scaled quantum chemical force field, and compared with an earlier normal coordinate analysis of [Ru(bpy)(3) ](2+) which was based upon experimental data alone, and the use of a simplified model. The results of the calculations provide a highly satisfactory fit to the experimental data and the normal coordinate analyses, in terms of potential energy distributions, allow a detailed understanding of the vibrational spectra of both complexes. Evidence is presented for Jahn-Teller distortion in the E-1 MLCT excited state. (C) 2008 Elsevier B.V. All rights reserved.

Vibrational spectra of α-amino acids in the zwitterionic state in aqueous solution and the solid state: DFT calculations and the influence of hydrogen bonding

The zwitterionic forms of the two simplest alpha-amino acids, glycine and l-alanine, in aqueous solution and the solid state have been modeled by DFT calculations. Calculations of the structures in the solid state, using PW91 or PBE functionals, are in good agreement with the reported crystal structures, and the vibrational spectra computed at the optimized geometries provide a good fit to the observed IR and Raman spectra in the solid state. DFT calculations of the structures and vibrational spectra of the zwitterions in aqueous solution at the B3-LYP/cc-pVDZ level were found to require both explicit and implicit solvation models. Explicit solvation was modeled by inclusion of five hydrogen-bonded water molecules attached to each of the five possible hydrogen-bonding sites in the zwitterion and the integration equation formalism polarizable continuum model (IEF-PCM) was employed, providing a satisfactory fit to observed IR and Raman spectra. Band assignments are reported in terms of potential-energy distributions, which differ in some respects to those previously reported for glycine and l-alanine.

DFT and in situ EXAFS investigation of gold/ceria-zirconia low-temperature water gas shift catalysts: Identification of the nature of the active form of gold

A combined experimental and theoretical investigation of the nature of the active form of gold in oxide-supported gold catalysts for the water gas shift reaction has been performed. In situ extended X-ray absorption fine structure (EXAFS) and X-ray absorption near-edge structure (XANES) experiments have shown that in the fresh catalysts the gold is in the form of highly dispersed gold ions. However, under water gas shift reaction conditions, even at temperatures as low as 100 degrees C, the evidence from EXAFS and XANES is only 14 consistent with rapid, and essentially complete, reduction of the gold to form metallic clusters containing about 50 atoms. The presence of Au-Ce distances in the EXAFS spectra, and the fact that about 15% of the gold atoms can be reoxidized after exposure to air at 150 degrees C, is indicative of a close interaction between a fraction (ca. 15%) of the gold atoms and the oxide support. Density functional theory (DFT) calculations are entirely consistent with this model and suggest that an important aspect of the active and stable form of gold under water gas shift reaction conditions is the location of a partially oxidized gold (Audelta+) species at a cerium cation vacancy in the surface of the oxide support. It is found that even with a low loading gold catalysts (0.2%) the fraction of ionic gold under water gas shift conditions is below the limit of detection by XANES (<5%). It is concluded that under water gas shift reaction conditions the active form of gold comprises small metallic gold clusters in intimate contact with the oxide support.

DFT studies of long-chain FAMEs: theoretical justification for determining chain length and unsaturation from experimental Raman spectra

Density functional calculations, using B3LPY/6-31G(d) methods, have been used to investigate the conformations and vibrational (Raman) spectra of a series of long-chain, saturated fatty acid methyl esters (FAMEs) with the formula CH2nO2 (n = 5-21) and two series of unsaturated FAMEs. The calculations showed that the lowest energy conformer within the saturated FAMEs is the simple (all-trans) structure and, in general, it was possible to reproduce experimental data using calculations on only the all-trans conformer. The only exception was C6H12O2, where a second low-lying conformer had to be included in order to correctly simulate the experimental Raman spectrum. The objective of the work was to provide theoretical justification for the methods that are commonly used to determine the properties of the fats and oils, such as chain length and degree of unsaturation, from experimental Raman data. Here it is shown that the calculations reproduce the trends and calibration curves that are found experimentally and also allow the reasons for the failure of what would appear to be rational measurements to be understood. This work shows that although the assumption that each FAME can simply be treated as a collection of functional groups can be justified in some cases, many of the vibrational modes are complex motions of large sections of the molecules and thus would not be expected to show simple linear trends with changes in structure, such as increasing chain length and/or unsaturation. Simple linear trends obtained from experimental data may thus arise from cancellation of opposing effects, rather than reflecting an underlying simplicity.

Conformations, vibrational frequencies and Raman intensities of short-chain fatty acid methyl esters using DFT with 6-31G(d) and Sadlej pVTZ basis sets

Density functional calculations, using B3LPY/6-31G(d) methods, have been used to investigate the conformations and vibrational (Raman) spectra of three short-chain fatty acid methyl esters (FAMEs) with the formula CnH2nO2 (n = 3-5). In all three FAMEs, the lowest energy conformer has a simple 'all-trans' structure but there are other conformers, with different torsions about the backbone, which lie reasonably close in energy to the global minimum. One result of this is that the solid samples we studied do not appear to consist entirely of the lowest energy conformer. Indeed, to account for the 'extra' bands that were observed in the Raman data but were not predicted for the all-trans conformer, it was necessary to add-in contributions from other conformers before a complete set of vibrational assignments could be made. Provided this was done, the agreement between experimental Raman frequencies and 6-31G(d) values (after scaling) was excellent, RSD = 12.6 cm(-1). However, the agreement between predicted and observed intensities was much less satisfactory. To confirm the validity of the approach followed by the 6-3 1 G(d) basis set, we used a larger basis set, Sadlej pVTZ, and found that these calculations gave accurate Raman intensities and simulated spectra (summed from two different conformers) that were in quantitative agreement with experiment. In addition, the unscaled Sadlej pVTZ, and the scaled 6-3 1 G(d) calculations gave the same vibrational mode assignments for all bands in the experimental data. This work provides the foundation for calculations on longer-chain FAMEs (which are closer to those found as triglycerides in edible fats and oils) because it shows that scaled 6-3 1 G(d) calculations give equally accurate frequency predictions, and the same vibrational mode assignments, as the much more CPU-expensive Sadlej pVTZ basis set calculations.

A DFT study of the chain growth probability in Fischer-Tropsch synthesis

The chain growth probability (alpha value) is one of the most significant parameters in Fischer-Tropsch (FT) synthesis. To gain insight into the chain growth probability, we systematically studied the hydrogenation and C-C coupling reactions with different chain lengths on the stepped Co(0001) surface using density functional theory calculations. Our findings elucidate the relationship between the barriers of these elementary reactions and the chain length. Moreover, we derived a general expression of the chain growth probability and investigated the behavior of the alpha value observed experimentally. The high methane yield results from the lower chain growth rate for C-1 + C-1 coupling compared with the other coupling reactions. After C-1, the deviation of product distribution in FT synthesis from the Anderson-Schulz-Flory distribution is due to the chain length-dependent paraffin/olefin ratio. (C) 2008 Elsevier Inc. All rights reserved.

Chain growth mechanism in Fischer-Tropsch synthesis: A DFT study of C-C coupling over Ru, Fe, Rh, and Re surfaces

A quantitative approach is used to understand the chain growth mechanism in FT synthesis on the Ru, Fe, Rh, and Re surfaces. The C-C coupling reactions are extensively calculated on the stepped metal surfaces. Combining the coupling barriers and reactant stabilities, we investigate the reaction rates of all possible C, + C-1 coupling pathways on the metal surfaces. It is found that (i) all the transition-state structures are similar on these surfaces, while some coupling barriers are very different; (ii) the dominant chain growth pathways on these surfaces are different: C + CH and CH + CH on Rh and Ru surfaces, C + CH3 on Fe surface, and C + CH on Re surface. The common features of the major coupling reactions together with those on the Co surface are also discussed.

Reactivity of the 4d transition metals toward N hydrogenation and NH dissociation: A DFT-based HSAB analysis

The density functional theory (DFT) based hard-soft acid-base (HSAB) reactivity indices, including the electrophilicity index, have been successfully applied to many areas of molecular chemistry. In this work we test the applicability of such an approach to fundamental surface chemistry. We have considered, as prototypical surface reactions, both the hydrogenation of atomic nitrogen and the dissociative adsorption of the NH molecular radical. By use of a DFT methodology, the minimum energy reaction pathways, and corresponding reaction barriers, of the above reactions over Zr(001), Nb(110), Mo(110), Tc(001), Ru(001), Rh(111), and Pd(111) have been determined. By consideration of the chemical potential and chemical hardness of the surface metal atoms, and the principle of electronegativity equalization, it is found that the charge transferred to the NH radical during the process of dissociative adsorption correlates very well with that determined by Mulliken population analysis. Furthermore, it is found that the stability of the NH/surface transition state complex relates directly to this charge transfer and that the trend in transition state stability predicted by a HSAB; treatment correlates very strongly with that determined by DFT calculations. With regards to N hydrogenation, we find that during the course of the reaction, H loses cohesion to the surface, as it must migrate from a 3-fold hollow site to either a bridge or top site, to react with N. Partial density of states (PDOS) and Mulliken population analysis reveal that this loss of bonding is accompanied by charge transfer from H to the surface metal atoms. Moreover, by simple modeling, we show that the reaction barriers are directly proportional to this mandatory charge transfer. Indeed, it is found that the reaction barriers correlate very well with the electrophilicity index of the metal atoms.

The temperature dependence of the adsorption of NO on Pt{211}: A RAIRS and DFT investigation

RAIRS experiments have been performed to investigate the adsorption of NO on Pt{211}. Results show that adsorption is complex and strongly temperature dependent. At 307 K, three bands are seen at saturation with frequencies of 1801, 1609, and 1576 cm(-1). However, at 120 K only two bands, at 1688 and 1620 cm(-1), are observed. To help with the assignment of these vibrational bands, DFT calculations were also performed. The calculations show that a bridged NO species, bonded to the step edge, is the most stable species on the surface and gives rise to the band observed at 1610-1620 cm(-1). The calculations also suggest that the temperature dependence of NO adsorption on Ptf{211} can be assigned to NO dissociation which occurs at room temperature but not at 120 K. In particular, the RAIRS band observed at 1801 cm(-1), which is observed on adsorption at 307 K but not at 120 K, is tentatively assigned to the formation of an O-NO complex. This species forms when a NO molecule bonds on top of an O atom, which results from the dissociation of NO on the Pt{211} surface at room temperature.

Resonance Raman and DFT studies of tetra-tert-butyl porphine: Assignment of strongly enhanced distortion modes in a ruffled porphyrin

The free-base form of tetra-tert-butyl porphine (TtBP), which has extremely bulky meso substituents, is severely distorted from planarity, with a ruffling angle of 65.5degrees. The resonance Raman spectrum of TtBP (lambda(ex) = 457.9 nm) and its d(2), d(8), and d(10) isotopomers have been recorded, and while the spectra show high-frequency bands similar to those observed for planar meso-substituted porphyrins, there are several additional intense bands in the low-frequency region. Density functional calculations at the B3-LYP/6-31G(d) level were carried out for all four isotopomers, and calculated frequencies were scaled using a single factor of 0.98. The single factor scaling approach was validated on free base porphine where the RMS error was found to be 14.9 cm(-1). All the assigned bands in the high-frequency (> 1000 cm(-1)) region of TtBP were found to be due to vibrations similar in character to the in-plane skeletal modes of conventional planar porphyrins. In the low-frequency region, two of the bands, assigned as nu(8) (ca. 330 cm(-1)) and nu(16) (ca. 540 cm(-1)), are also found in planar porphyrins such as tetra-phenyl porphine (TPP) and tetra-iso-propyl porphine (IPP). Of the remaining three very strong bands, the lowest frequency band was assigned as gamma(12) (pyr swivel, obsd 415 cm(-1), calcd 407 cm(-1) in do). The next band, observed at 589 cm-1 in the do compound (calcd 583 cm(-1)), was assigned as a mode whose composition is a mixture of modes that were previously labeled gamma(13) (gamma(CmCaHmCa)) andy gamma(11) (pyr fold(asym)) in NiOEP. The final strong band, observed at 744 cm(-1) (calcd 746 cm(-1)), was assigned to a mode whose composition is again a mixture of gamma(11) and gamma(13), although here it is gamma(11) rather than gamma(13) which predominates. These bands have characters and positions similar to those of three of the four porphyrin ring-based, weak bands that have previously been observed for NiTPP. In addition there are several weaker bands in the TtBP spectra that are also