New insights into NO-carbon and N2O-carbon reactions from quantum mechanical calculations

Density functional theory calculations were used to investigate the mechanisms of NO-carbon and N2O-carbon reactions. It was the first time that the importance of surface nitrogen groups was addressed in the kinetic behaviors of the NO-carbon reaction. It was found that the off-plane nitrogen groups that are adjacent to the zigzag edge sites and in-plane nitrogen groups that are located on the armchair sites make the bond energy of oxygen desorption even ca. 20% lower than that of the off-plane epoxy group adjacent to zigzag edge sites and in-plane o-quinone oxygen atoms on armchair sites; this may explain the reason why the experimentally obtained activation energy of the NO-carbon reaction is ca. 20% lower than that of the O-2-carbon reaction over 923 K. A higher ratio of oxygen atoms can be formed in the N2O-carbon reaction, because of the lower dissociation energy of N2O, which results in a higher ratio of off-plane epoxy oxygen atoms. The desorption energy of semiquinone with double adjacent off-plane oxygen groups is ca. 20% less than that of semiquinone with only one adjacent off-plane oxygen group. This may be the reason why the activation energy of N2O is also ca. 20% less than that of the O-2-carbon reaction. The new mechanism can also provide a good qualitative comparison for the relative reaction rates of NO-, N2O-, and O-2-carbon reactions. The anisotropic characters of these gas-carbon reactions can also be well explained.

Transition-state-theory calculations for reactions of O(3P) with halogenated olefins

The potential energy surfaces for the reactions of atomic oxygen in its ground electronic state, O(P-3), with the olefins: CF2=CCl2 and CF2=CF - CF3, have been characterized using ab initio molecular orbital calculations. Geometry optimization and vibrational frequency calculations were performed for reactants, transition states and products at the MP2 and QCISD levels of theory using the 6-31G(d) basis set. This database was then used to calculate the rate constants by means of Transition-State-Theory. To obtain a better reference and to test the reliability of the activation barriers we have also carried out computations using the CCSD(T)(fc)/6-311Gdagger, MP4(SDQ)(fc)/CBSB4 and MP2(fc)/CBSB3 single point energy calculations at both of the above levels of theory, as well as with the composite CBS-RAD procedure ( P. M. Mayer, C. J. Parkinson, D. M. Smith and L. Radom, J. Chem. Phys., 1998, 108, 604) and a modi. cation of this approach, called: CBS-RAD( MP2, MP2). It was found that the kinetic parameters obtained in this work particularly with the CBS-RAD ( MP2, MP2) procedure are in reasonable agreement with the experimental values. For both reactions it is found that the channels leading to the olefin double-bond addition predominates with respect to any other reaction pathway. However, on account of the different substituents in the alkenes we have located, at all levels of theory, two transition states for each reaction. Moreover, we have found that, for the reactions studied, a correlation exists between the activation energies and the electronic structure of the transition states which can explain the influence of the substituent effect on the reactivity of the halo-olefins.

Gas-phase solvolysis type reactions of SiCl3+ cations

Gas-phase SiCl3+ ions undergo sequential solvolysis type reactions with water, methanol, ammonia, methylamine and propylene. Studies carried out in a Fourier Transform mass spectrometer reveal that these reactions are facile at 10-8 Torr and give rise to substituted chlorosilyl cations. Ab initio and DFT calculations reveal that these reactions proceed by addition of the silyl cation to the oxygen or nitrogen lone pair followed by a 1,3-H migration in the transition state. These transition states are calculated to lie below the energy of the reactants. By comparison, hydrolysis of gaseous CCl3+ is calculated to involve a substantial positive energy barrier.

A theoretical study of the tautomerism and vibrational spectra of 4,5-diamine-2,6-dimercaptopyrimidine

The 4,5-diamine-2,6-dimercaptopyrimidine (DADMcP) compound is an interesting multifunctional species exhibiting a rather complex tautomerism, encompassing nine tautomeric forms. Investigation of tautomerism in this compound has been carried out by means of FTIR spectroscopy, in association with ab-initio HF/SCF and DFT calculations. According to this study three tautomers are energetically favored; the thione form being the most stable one. The theoretical vibrational spectra of such tautomeric forms have been successfully simulated by means of DFT calculations, allowing the elucidation and assignment of the complex composition of the vibrational bands observed for the mixture of isomers.

Study of the optical and magnetic properties of pyrimidine in water combining PCM and QM/MM methodologies

The solvent effects on the low-lying absorption spectrum and on the (15)N chemical shielding of pyrimidine in water are calculated using the combined and sequential Monte Carlo simulation and quantum mechanical calculations. Special attention is devoted to the solute polarization. This is included by an iterative procedure previously developed where the solute is electrostatically equilibrated with the solvent. In addition, we verify the simple yet unexplored alternative of combining the polarizable continuum model (PCM) and the hybrid QM/MM method. We use PCM to obtain the average solute polarization and include this in the MM part of the sequential QM/MM methodology, PCM-MM/QM. These procedures are compared and further used in the discrete and the explicit solvent models. The use of the PCM polarization implemented in the MM part seems to generate a very good description of the average solute polarization leading to very good results for the n-pi* excitation energy and the (15)N nuclear chemical shield of pyrimidine in aqueous environment. The best results obtained here using the solute pyrimidine surrounded by 28 explicit water molecules embedded in the electrostatic field of the remaining 472 molecules give the statistically converged values for the low lying n-pi* absorption transition in water of 36 900 +/- 100 (PCM polarization) and 36 950 +/- 100 cm(-1) (iterative polarization), in excellent agreement among one another and with the experimental value observed with a band maximum at 36 900 cm(-1). For the nuclear shielding (15)N the corresponding gas-water chemical shift obtained using the solute pyrimidine surrounded by 9 explicit water molecules embedded in the electrostatic field of the remaining 491 molecules give the statistically converged values of 24.4 +/- 0.8 and 28.5 +/- 0.8 ppm, compared with the inferred experimental value of 19 +/- 2 ppm. Considering the simplicity of the PCM over the iterative polarization this is an important aspect and the computational savings point to the possibility of dealing with larger solute molecules. This PCM-MM/QM approach reconciles the simplicity of the PCM model with the reliability of the combined QM/MM approaches.

Piezomagnetic behavior of Co-doped ZnO nanoribbons

Using ab initio total energy calculations, we show that bilayer systems of ZnO nanoribbons, (ZnO)(2)NR, doped with Co atoms exhibit a piezomagnetic behavior. We find the formation of energetically stable zigzag chains of Co atoms along the edge sites of (ZnO)(2)NR's, Co(Zn(chain))-(ZnO)(2)NR. At the ground state, the antiferromagnetic and the ferromagnetic states are very close in energy, whereas upon longitudinal stretch, parallel to the nanoribbon growth direction, it becomes ferromagnetic. Further electronic structure calculations indicate that not only the magnetic state but also the electronic structure of CoZn(chain)-(ZnO)(2)NR can be tuned by the mechanical stretch. In this case, we find that stretched NR's exhibit dispersive unpaired electronic states within the (ZnO)(2)NR band gap.

Splitting of the zero-energy edge states in bilayer graphene

Bilayer graphene nanoribbons with zigzag termination are studied within the tight-binding model. We also include single-site electron-electron interactions via the Hubbard model within the unrestricted Hartree-Fock approach. We show that either the interactions between the outermost edge atoms or the presence of a magnetic order can cause a splitting of the zero-energy edge states. Two kinds of edge alignments are considered. For one kind of edge alignment (?) the system is nonmagnetic unless the Hubbard parameter U becomes greater than a critical value Uc. For the other kind of edge alignment (?) the system is magnetic for any U>0. Our results agree very well with ab initio density functional theory calculations.

Atomistic molecular dynamics study of interface formation: Al on poly(p-phenylene vinylene)

We model interface formation by metal deposition on the conjugated polymer poly-para-phenylene vinylene, studying direct aluminum and layered aluminum-calcium structures Al/PPV and Al/Ca/PPV. To do that we use classical molecular dynamics simulations, checked by ab initio density-functional theory calculations, for selected relevant configurations. We find that Al not only migrates easily into the film, with a strong charge transfer to the neighboring chains, but also promotes rearrangement of the polymer in the interfacial region to the hexagonal structure. On the other hand, our results indicate that a thin Ca layer is sufficient to protect the film and maintain a well-defined metal/polymer interface, and that also a thin Al capping layer may protect the whole from environmental degradation.

Molecular dynamics simulation of liquid trimethylphosphine

Structural and dynamical properties of liquid trimethylphosphine (TMP), (CH(3))(3)P, as a function of temperature is investigated by molecular dynamics (MD) simulations. The force field used in the MD simulations, which has been proposed from molecular mechanics and quantum chemistry calculations, is able to reproduce the experimental density of liquid TMP at room temperature. Equilibrium structure is investigated by the usual radial distribution function, g(r), and also in the reciprocal space by the static structure factor, S(k). On the basis of center of mass distances, liquid TMP behaves like a simple liquid of almost spherical particles, but orientational correlation due to dipole-dipole interactions is revealed at short-range distances. Single particle and collective dynamics are investigated by several time correlation functions. At high temperatures, diffusion and reorientation occur at the same time range as relaxation of the liquid structure. Decoupling of these dynamic properties starts below ca. 220 K, when rattling dynamics of a given TMP molecules due to the cage effect of neighbouring molecules becomes important. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3624408]

Electron theoretical investigation of the stability of the B2-TiFe compound

The metastable phase diagram of the BCC-based ordering equilibria in the Ti-Fe system has been calculated using a truncated cluster expansion, through the combination of FP-LAPW and cluster variation method (CVM) in the irregular tetrahedron cluster approximation. The results are compared with phenomenological CVM assessments of the system and suggest that the value for the experimental formation enthalpy of the B2-TiFe compound should be significantly more negative than the currently assessed value. (C) 2008 Elsevier B.V. All rights reserved.

Synthesis and evaluation of sesquiterpene lactone inhibitors of phospholipase A(2) from Bothrops jararacussu

Several sesquiterpene lactone were synthesized and their inhibitive activities on phospholipase A(2) (PLA(2)) from Bothrops jararacussu venom were evaluated. Compounds Lac01 and Lac02 were efficient against PLA(2) edema-inducing, enzymatic and myotoxic activities and it reduces around 85% of myotoxicity and around 70% of edema-inducing activity. Lac05-Lac08 presented lower efficiency in inhibiting the biological activities studied and reduce the myotoxic and edema-inducing activities around only 15%. The enzymatic activity was significantly reduced. The values of inhibition constants (K(1)) for Lac01 and Lac02 were approximately 740 mu M, and for compounds Lac05-Lac08 the inhibition constants were approximately 7.622-9.240 mu M. The enzymatic kinetic studies show that the sesquiterpene lactones inhibit PLA(2) in a non-competitive manner. Some aspects of the structure-activity relationships (topologic, molecular and electronic parameters) were obtained using ab initio quantum calculations and analyzed by chemometric methods (HCA and PCA). The quantum chemistry calculations show that compounds with a higher capacity of inhibiting PLA(2) (Lac01-Lac04) present lower values of highest occupied molecular orbital (HOMO) energy and molecular volume (VOL) and bigger values of hydrophobicity (LogP). These results indicate some topologic aspects of the binding site of sesquiterpene lactone derivatives and PLA(2). (C) 2010 Elsevier Ltd. All rights reserved.

Computational chemistry on Fujitsu vector-parallel processors: Development and performance of applications software

In this and a preceding paper, we provide an introduction to the Fujitsu VPP range of vector-parallel supercomputers and to some of the computational chemistry software available for the VPP. Here, we consider the implementation and performance of seven popular chemistry application packages. The codes discussed range from classical molecular dynamics to semiempirical and ab initio quantum chemistry. All have evolved from sequential codes, and have typically been parallelised using a replicated data approach. As such they are well suited to the large-memory/fast-processor architecture of the VPP. For one code, CASTEP, a distributed-memory data-driven parallelisation scheme is presented. (C) 2000 Published by Elsevier Science B.V. All rights reserved.

Solid-State Experimental and Theoretical Investigation of the Ammonium Salt of Croconate Violet, a Pseudo-Oxocarbon Ion

The present work describes the crystal structure, vibrational spectra, and theoretical calculations of ammonium salts of 3,5-bis-(dicyanomethylene)cyclopentane-1,2,4-trionate, (NH(4))(2)(C(11)N(4)O(3)) [(NH(4))(2)CV], also known as ammonium croconate violet. This compound crystallizes in triclinic P (1) over bar and contains two water molecules per unit formula. The crystal packing is stabilized by hydrogen bonds involving water molecules and ammonium cations, giving rise to a 3D polymeric arrangement. In this structure, a pi-stacking interaction is not observed, as the smaller centroid-centroid distance is 4.35 angstrom. Ab initio electronic structure calculations under periodic boundary conditions were performed to predict vibrational and electronic properties. The vibrational analysis was used to assist the assignments of the Raman and infrared bands. The solid structure was optimized and characterized as a minimum in the potential-energy surface. The stabilizing intermolecular hydrogen bonds in the crystal Structure were characterized by difference charge-density analysis. The analysis of the density of states of (NH(4))(2)CV gives an energy gap of 1.4 eV with a significant contribution of carbon and nitrogen 2p states for valence and conduction bands.

2-pyrazinylnitrene and 4-pyrimidylnitrene. Ring expansion to 1,3,5-triazacyclohepta-1,2,4,6-tetraene and ring opening to (2-isocyanovinyl)carbodiimide

Tetrazolo[1,5-a]pyrazine/2-azidopyrazine 9T/9A undergo photolysis in Ar matrix at cryogenic temperatures to yield 1,3,5-triazacyclohepta-1,2,4,6-tetraene 21 as the first observable intermediate, and 1-cyanoimidazole 11 and (2-isocyanovinyl)carbodiimide 22 as the final products. The latter tautomerizes to 2-(isocyanovinyl)cyanamide 23 on warming to 40 K. The same intermediate 21 and the same final products are obtained on matrix photolysis of the isomeric tetrazolo[1,5-c]pyrimidine/4-azidopyrimidine 24T/24A. These photolysis results as well as those of the previously reported thermal ring contraction of N-15-labeled 2-pyrazinyl- and 4-pyrimidylnitrenes to 1-cyanoimidazoles can all be rationalized in terms of selective ring opening of 21 or nitrine 10 to a nitrile ylide zwitterion 28 prior to formation of the final products, 11 and 22. The results are supported by high-level ab initio and DFT calculations (CASPT2-CASSCF(6,6), G3(MP2), and B3LYP/6-31+G*) of the energies and IR spectra of the intermediates and products.

2-cyano(2-pyridyl)ketene

Flash vacuum thermolysis of quinolizinones is a new way of generating ketenes. The title ketene is obtained from 1-cyano-2-hydroxyquinolizine-4-one and characterized by its Ar matrix infrared spectrum. (C) Wiley-VCH Verlag GmbH, 69451 Weinheim, Germany 2002.