Radiationless decay mechanism of cytosine: an ab initio study with comparisons to the fluorescent analogue 5-methyl-2-pyrimidinone

The ultrafast radiationless decay mechanism of photoexcited cytosine has been theoretically supported by exploring the important potential energy surfaces using multireference configuration-interaction ab initio methods for the gas-phase keto-tautomer free base. At vertical excitation, the bright state is S1 (pipi*) at 5.14 eV, with S2 (nNpi*) and S3 (nOpi*) being dark states at 5.29 and 5.93 eV, respectively. Minimum energy paths connect the Franck-Condon region to a shallow minimum on the pipi* surface at 4.31 eV. Two different energetically accessible conical intersections with the ground state surface are shown to be connected to this minimum. One pathway involves N3 distorting out of plane in a sofa conformation, and the other pathway involves a dihedral twist about the C5-C6 bond. Each of these pathways from the minimum contains a low barrier of 0.14 eV, easily accessed by low vibronic levels. The path involving the N3 sofa distortion leads to a conical intersection with the ground state at 4.27 eV. The other pathway leads to an intersection with the ground state at 3.98 eV, lower than the minimum by about 0.3 eV. Comparisons with our previously reported study of the fluorescent cytosine analogue 5-methyl-2-pyrimidinone (5M2P) reveal remarkably similar conformational distortions throughout the decay pathways of both bases. The different photophysical behavior between the two molecules is attributed to energetic differences. Vertical excitation in cytosine occurs at a much higher energy initially, creating more vibrational energy than 5M2P in the Franck-Condon region, and the minimum S1 energy for 5M2P is too low to access an intersection with the ground state, causing population trapping and fluorescence. Calculations of vertical excitation energies of 5-amino-2-pyrimidinone and 2-pyrimidinone reveal that the higher excitation energy of cytosine is likely due to the presence of the amino group at the 4-position.     

Catalytic reaction mechanism of L-lactate dehydrogenase: an ab initio study

Studies on the catalytic reaction mechanism of L-lactate dehydrogenase have been carried out by using quantum chemical ab initio calculation at HF/6-31G* level. It is found that the interconversion reaction of pyruvate to L-lactate is dominated by the hydride ion HR- transfer, and the transfers of the hydride ion HR and proton HR are a quasi-coupled process, in which the energy barrier of the transition state is about 168.37 kJ/mol. It is shown that the reactant complex is 87.61 kJ/mol lower, in energy, than the product complex. The most striking features in our calculated results are that pyridine ring of the model cofactor is a quasi-boat-like configuration in the transited state, which differs from a planar conformation in some previous semiempirical quantum chemical studies. On the other hand, the similarity in the structure and charge between the HR transfer process and the hydrogen bonding with lower barrier indicates that the HR transfer process occurs by means of an unusual manner. In addition,     

Size-expanded yDNA bases: an ab initio study

xDNA and yDNA are new classes of synthetic nucleic acids characterized by having base-pairs with one of the bases larger than the natural congeners. Here these larger bases are called x- and y-bases. We recently investigated and reported the structural and electronic properties of the x-bases (Fuentes-Cabrera et al. J. Phys. Chem. B 2005, 109, 21135-21139). Here we extend this study by investigating the structure and electronic properties of the y-bases. These studies are framed within our interest that xDNA and yDNA could function as nanowires, for they could have smaller HOMO-LUMO gaps than natural DNA. The limited amount of experimental structural data in these synthetic duplexes makes it necessary to first understand smaller models and, subsequently, to use that information to build larger models. In this paper, we report the results on the chemical and electronic structure of the y-bases. In particular, we predict that the y-bases have smaller HOMO-LUMO gaps than their natural congeners, which is an encouraging result for it indicates that yDNA could have a smaller HOMO-LUMO gap than natural DNA. Also, we predict that the y-bases are less planar than the natural ones. Particularly interesting are our results corresponding to yG. Our studies show that yG is unstable because it is less aromatic and has a Coulombic repulsion that involves the amino group, as compared with a more stable tautomer. However, yG has a very small HOMO-LUMO gap, the smallest of all the size-expanded bases we have considered. The results of this study provide useful information that may allow the synthesis of an yG-mimic that is stable and has a small HOMO-LUMO gap.     

Thermal decomposition of decalin: an ab initio study

Density functional theory calculations (B3LYP and BH&HLYP functionals) of the potential energy surface have been performed to investigate the mechanisms of decalin breakdown, and the Rice-Ramsperger-Kassel-Marcus and transition state theory methods have been used to compute the high-pressure limit thermal rate constants for the new reaction pathways. The new pathways connect decalin to five primary monoaromatic species: benzene, toluene, styrene, ethylbenzene, and xylene. The reactions used for the new routes are carbon-carbon bond cleavage reaction, dissociation reaction, and hydrogen abstraction and addition reactions. A kinetic analysis was performed for pyrolytic conditions, and benzene, toluene, and xylene were identified as major products.     

Vibrational spectroscopic studies and ab initio calculations of 5-methyl-2-(p-fluorophenyl)benzoxazole

FT-Raman and FT-IR spectra of 5-methyl-2-(p-fluorophenyl)benzoxazole were recorded and analysed. The vibrational frequencies of the compound have been computed using the Hartree-Fock/6-31G* basis and compared with the experimental values.     

Vibrational spectroscopic studies and ab initio calculations of 5-methyl-2-(p-methylaminophenyl)benzoxazole

FT-IR spectra of 5-methyl-2-(p-methylaminophenyl)benzoxazole was recorded and analysed. The vibrational frequencies of the compound have been computed using the Hartree-Fock/6-31G* basis and compared with the experimental values.     

Wheland intermediates: an ab initio valence bond study

The traditional resonance model for electrophilic attacks on substituted aromatic rings is revisited using high level valence bond (VB) calculations. A large π-donation is found in the X = NH(2) case and a weaker one for the X = Cl case, not only for ortho and para isomers but also for the meta case, which can be explained by considering five (not three) fundamental VB structures. No substantial π-effect is found in the X = NO(2) case, generally viewed as π-attractive.     

Conformational stability of a model macrocycle tetraamide: an ab initio study

Ab initio calculations are carried out to investigate the conformational stability of a model macrocyle tetraamide. The four amide groups in the selected model are present in the sequence: -(O=CNH)-Ph-(NHC=O)-CH=CH-(O=CNH)-Ph-(NHC=O)-CH=CH-. In this sequence, two phenyl rings and two ethene groups act as bridges between the amide units. Each amide motif bonds to a phenyl ring through its amide nitrogen and to an ethene group through its amide carbon. Four clearly distinct minimum-energy conformations are found upon full geometry optimization using the B3LYP/6-31+G(d) method. Frequency calculations using the same method confirm that the four conformations are indeed minima in the macrocycle potential energy surface. Relative to the most stable conformer, the other conformations are higher in energy by 0.86, 2.09, and 9.17 kcal/mol, respectively, at the MP2/6-31+G(d,p) level. The stability of the macrocycle conformations is correlated primarily to the existence and strength of intramolecular N-H...O=C hydrogen bonds. Additional stability to the conformations is found to come from weak Ph-H...O=C hydrogen bonding between a carbonyl oxygen and a hydrogen atom of a phenyl group. Solvent effects play an important role in the relative energies of the various conformations, as indicated by the simple SCRF = dipole model calculations for the case of aqueous solution.     

Pseudorotation motion in tetrahydrofuran: an ab initio study

The use of different models based on experimental information about the observed level splitings, rotational constants, and far-infrared transition frequencies leads to different predictions on the equilibrium geometry for tetrahydrofuran. High-level ab initio calculations [coupled cluster singles, doubles (triples)/complete basis set (second order Moller-Plesset triple, quadrupole, quintuple)+zero-point energy(anharmonic)] suggest that the equilibrium conformation of tetrahydrofuran is an envelope C(s) structure. The theoretical geometrical parameters might be helpful to plan further microwave spectroscopic studies in order to get a physical interpretation of the measurements.     

Towards SiC surface functionalization: an ab initio study

We present a microscopic model of the interaction and adsorption mechanism of simple organic molecules on SiC surfaces as obtained from ab initio molecular-dynamics simulations. Our results open the way to functionalization of silicon carbide, a leading candidate material for biocompatible devices.     

Molecular structure of 2-methylamino-5-[(5-methyl-2-benzoxazolinone-3-yl)methyl]-1,3,4-thiadiazole dihydrophosphate: a combined X-ray crystallographic and ab initio study

Crystal structure of the compound entitled 2-methylamino-5-[(5-methyl-2-benzoxazolinone-3-yl)methyl]-1,3,4-thiadiazole dihydrophosphate is determined using X-ray analysis and compared with the structure obtained from semiempirical and RHF methods at various levels of theory. RHF/6-31G(d) calculations offer the best conformity with X-ray results for bond lengths and bond angles. Moreover, at the result of the comparison of various combinations of basis sets and methods, it appears that there is not much gain in accuracy by using sophisticated methods.     

The structural effects of fluorine substitution in pyridine,pyrimidine, and s-triazine: An ab initio study

The structures of 2-fluoropyridine, 3-fluoropyridine, 4-fluoropyridine, pentafluoropyridine, 2-fluoropyrimidine, 5-fluoropyrimidine, and fluoro-s-triazine have been evaluated by the ab initio gradient method with a 4-21 basis set augmented as needed with polarization functions on nitrogen. The structural effects of fluorination on the parent heterocycle are very similar to the effects of fluorination on benzene studied earlier. The ring angle is enlarged by about 2° at the point of fluorination and the adjacent ring bonds are shortened, much more for an adjacent C-N than for C-C. Fluorination of pyridine at the C₂ position causes an in-crease of bond localization in the ring. Rotational constants calculated from the structural parameters corrected with standard offset values are in exceptionally close agreement with experimental constants where these are known.     

Electron attachment in ice-HCl clusters: an ab initio study

Experimental work has shown that small amounts of HCl strongly enhance electron capture in ice films. The purpose of the present study was to investigate the effect of adsorbed HCl on the interaction of electrons with small clusters of water. Studies were made with clusters of 6 and 12 water molecules with various geometries both with and without one HCl attached. A number of distinct HCl coordination motifs were examined. All of the neutral structures with HCl exhibited zero thresholds for electron attachment and formed dipole bound anionic states (DBS). The relaxation processes for these "initial DBS" depended on the number of H(2)O (n) and on the number and type of H-bonds to the HCl (x). The initial DBS of systems with only O-H...Cl H-binding underwent dissociative electron attachment (DEA), forming H atoms. Relaxation for systems with ClH...OH(2) bonds was more complex. For the two layer n = 12 systems with x = 2 or 3 the HCl proton moved to the nearest oxygen to form H(3)O(+). Then rearrangement of the proton network occurred, and the Cl(-) became solvated by three HO-H...Cl(-) bonds. The presence of Cl(-) and H(3)O(+) increases the dipole moment and the electron binding energy (EBE) of the network. Further stabilization is achieved by decay into deeper DBS electron traps and/or by reaction of the excess electron with H(3)O(+) to form H(*) atoms. The HCl(H(2)O)(6) clusters with a single Cl-H...OH(2) bond behaved differently. They increased their stability by becoming more linear. This raised the dipole moment and the EBE therefore increased, reducing the total energy. None of these species showed any signs of increasing the number of H-bonds to Cl. The implication of these observations for the interpretation of the results of the experiments with 0.2 monolayer of HCl on 5 monolayer of H(2)O at 20 K, and on the possible role of cosmic ray-induced ionization in polar stratospheric clouds in ozone depletion is discussed.     

The low-lying electronic states of H2CN+: a preliminary ab initio study

Ab initio molecular-orbital calculations have been carried out on the low-lying triplet and singlet electronic states of the H₂CN⁺ cation, at the SCF and Møller-Plesset levels of theory. Both triplet ³A₂ and ³B₂ electronic states have similar energies. The barriers to isomerization to the ³A″ and ³A' electronic states of HCNH⁺ are estimated. It appears that ³A₂ and ³B₂ states are stable towards both isomerization and dissociation. The results of mass spectroscopic experiments involving H₂CN⁺ are discussed.     

Current rectification by asymmetric molecules: an ab initio study

Current rectification effect in an asymmetric molecule HCOO-C6H4-(CH2)n sandwiched between two aluminum electrodes has been studied using an ab initio nonequilibrium Green's function method. The conductance of the system decreases exponentially with the increasing number n of CH2. The phenomenon of current rectification is observed such that a very small current appears at negative bias and a sharp negative differential resistance at a critical positive bias when n>or=2. The rectification effect arises from the asymmetric structure of the molecule and the molecule-electrode couplings. A significant rectification ratio of approximately 38 can be achieved when n=5.     

Glycine interaction with carbon nanotubes: an ab initio study

The interaction of the glycine radical on the side walls of both armchair and zigzag single walled carbon nanotubes is investigated by density functional theory. It is found that the interaction potential of the N-centered glycine radical with the tubes has a minimum of 16.9 (armchair) and 20.2 (zigzag) kcal/mol with respect to the dissociation products. In contrast, the C-centered radical, which is 22.7 kcal/mol lower in energy than the N-centered radical, does not form stable complexes with both types of carbon nanotubes.     

Solvent-dependent stereoselectivity in a Still-Wittig rearrangement: an experimental and ab initio study

[see reaction]. The Still-Wittig rearrangement gave opposite selectivities for (Z:E)-alkenes in THF (3:1) vs toluene (1:3) in the synthesis of serine-proline dipeptide amide isosteres. Four transition states leading to (Z)-and (E)-alkenes with THF and without (representing toluene) were identified by ab initio calculations at the 3-21G* level. The calculated (Z:E)-ratios with THF (4.7:1) and without THF (1:3.2) suggested that the transition state geometries and energies were well-represented by the calculations.     

Addition of anions to carbonyl compounds: an ab initio study

Energies for the addition of anionic nucleophiles, Z(-), to carbonyl compounds, XYCO, are calculated at the G2(MP2) level of theory. The substituents X, Y, and Z are taken from the set {H, CH3, NH2, OH, F, CF3, CHCH2, CHO, CCH, and CN}. The basicity and, to a lesser extent, ionization potential of Z(-) were found to correlate with the enthalpy of addition of Z(-) to H2CO. The enthalpy of addition of Z(-) to XYCO relative to H2CO is largely independent of Z. The ordering of the enthalpies of addition for the series of XYCO's is rationalized. By using a thermodynamic cycle, the independence of this ordering from Z is attributed to the additivity of the inductive stabilization of XYZCO(-) by X and Y versus H2ZCO(-). A method for estimating the enthalpy of addition for nucleophile/carbonyl combinations not studied above is described and shown to give good results on a model system.     

First ultraviolet absorption band of methane: an ab initio study

Quantum mechanical calculations of the cross sections for photodissociation of CH4 and CD4 in the 1t2-->3s band are presented. The potential energy surfaces for the three states correlating with the 1 1T2 state at tetrahedral geometries are calculated. The elements of the (3x3) matrix representing the electronic Hamiltonian in the diabatic basis are expanded in powers of nuclear coordinates, up to the second order. The expansion coefficients are based on accurate multireference configuration interaction calculations. The electronically nonadiabatic dynamics is treated with the multiconfiguration time-dependent Hartree approach. All nine internal degrees of methane are included in the quantum dynamics simulations. The calculated cross section agrees well with experiment. Semiclassical calculations using the reflection principle suggest that the peaks in the spectrum correspond to the three adiabatic electronic states correlating with the 1 1T2 state at Td geometries. However, the non-Born-Oppenheimer terms in the Hamiltonian have a strong effect on the positions of the peaks in the absorption spectrum. The results of semiclassical calculations, which neglect these terms, are therefore quite different from the accurate quantum results and experiment.