AM1 and PM3 Calculations of the Potential Energy Surfaces for CH2OH Reactions with NO and NO2

The AM1 and PM3 molecular orbital methods have been utilized to investigate the reactions of CH20H with NO and NO2 PM3 and AM1 calculated heats of formation differ from experimental values by 8.6 and 18.8 kcal mol-', respectively. The dominant reaction of CH20H with NO is predicted to produce the adduct HOCH2N0, supporting the hypothesis of Pagsberg, Munk, Anastasi, and Simpson. Calculated activation energies for the NO2 system predict the formation of the adducts HOCH2N02 and HOCH20N0. In addition, the PM3 calculations predict that the abstraction reaction producing CH20 and HN02 is more likely than one producing CH20 and HONO from reactions of CH20H with NO2.

Disjoint nonclassical hydrocarbons have very unstable lowest-lying singlet states: a PM3 study

Earlier workers have suggested that disjoint hydrocarbons have nearly-degenerate lowest-lying singlet and triplet states while non-disjoint (or joint) hydrocarbons should be ground-state triplets. PM3 results for an appropriate selection of alternant hydrocarbons are inconsistent with that generalization: disjoint, nonclassical, alternant hydrocarbons show the strongest predilection for triplet ground states.

Hydrogen bonding of nucleotide base pairs: Application of the PM3 method

The ability of the pm3 semiempirical quantum mechanical method to reproduce hydrogen bonding in nucleotide base pairs was assessed. Results of pm3 calculations on the nucleotides 2′-deoxyadenosine 5′-monophosphate (pdA), 2′-deoxyguanosine 5′-monophosphate (pdG), 2′-deoxycytidine 5′-monophosphate (pdC), and 2′-deoxythymidine 5′-monophosphate (pdT) and the base pairs pdA–pdT, pdG–pdC, and pdG(syn)–pdC are presented and discussed. The pm3 method is the first of the parameterized nddo quantum mechanical models with any ability to reproduce hydrogen bonding between nucleotide base pairs. Intermolecular hydrogen bond lengths between nucleotides displaying Watson–Crick base pairing are 0.1–0.2 Å less than experimental results. Nucleotide bond distances, bond angles, and torsion angles about the glycosyl bond (χ), the C4′C5′ bond (γ), and the C5′O5′ bond (β) agree with experimental results. There are many possible conformations of nucleotides. pm3 calculations reveal that many of the most stable conformations are stabilized by intramolecular CHO hydrogen bonds. These interactions disrupt the usual sugar puckering. The stacking interactions of a dT–pdA duplex are examined at different levels of gradient optimization. The intramolecular hydrogen bonds found in the nucleotide base pairs disappear in the duplex, as a result of the additional constraints on the phosphate group when part of a DNA backbone. Sugar puckering is reproduced by the pm3 method for the four bases in the dT–pdA duplex. pm3 underestimates the attractive stacking interactions of base pairs in a B-DNA helical conformation. The performance of the pm3 method implemented in SPARTAN is contrasted with that implemented in MOPAC. At present, accurate ab initio calculations are too timeconsuming to be of practical use, and molecular mechanics methods cannot be used to determine quantum mechanical properties such as reaction-path calculations, transition-state structures, and activation energies. The pm3 method should be used with extreme caution for examination of small DNA systems. Future parameterizations of semiempirical methods should incorporate base stacking interactions into the parameterization data set to enhance the ability of these methods.

Sulfoxides as an intramolecular sulfenylating agent for indoles and diverse applications of the sulfide-sulfoxide redox cycle in organic chemistry

This dissertation involves study of various aspects of sulfoxide chemistry. Specifically designed t-butyl and propanenitrile sulfoxides tethered to indole-2-carboxamide were used as a source of intramolecular sulfenylating agents to synthesize novel indolo[3,2-b]-1-5-benzothiazepinones which are structurally analogous to the other biologically active benzothiazepinones. This study reveals that the intramolecular cyclization of sulfoxide follows an electrophilic sulfenylation (Sulfoxide Electrophilic Sulfenylation, SES) reaction pathway. Evidence of the absence of sulfenic acid as a transient reactive intermediate in such intramolecular cyclization is also provided. In another study, sulfoxide was used as a “protecting group” of thioether to synthesize 8-membered, indole substituted, thiazocine-2-acetic acid derivative via Ring Closing Metathesis (RCM). Protection (oxidation) of inert (to RCM) sulfide to sulfoxide followed by RCM produced cyclized product in good yields. Deprotection (reduction) of sulfoxide was achieved using Lawessons Reagent (L.R.). Application of the sulfide-sulfoxide redox cycle to solve the existing difficulties in using RCM methodology to thioethers is illustrated. A new design of a “molecular brake”, based on the sulfide-sulfoxide redox cycle is described. N-Ar rotation in simple isoindolines is controlled by the oxidation state of the proximate sulfur atom. Sulfide [S(II)] shows “free” [brake OFF] N-Ar rotation whereas sulfoxide displayed hindered [brake ON] N-Ar rotation. The semi-empirical molecular orbital (PM3) calculations revealed concerted pyramidalization of amidic nitrogen with N-Ar rotation.

Spectroscopic analysis of C80 doping with group III and group V elements using semiempirical PM3 molecular modelling technique

Semiempirical calculations at the level of PM3 of theory were carried out to study the structural and electronic properties of C80 and some of its doped derivatives with the elements of group III and V at the level of PM3 of theory. We have selected these elements to be substituted in the fullerene-C80 cage in order to show the effect of such structural change on the electronic properties of the molecules studied. The theoretical IR spectra, some of physical and chemical properties of the molecules studied are obtained and discussed.

Ability of the PM3 Quantum Mechanical Method to Model Intermolecular Hydrogen Bonding between Neutral Molecules

The PM3 semiempirical quantum-mechanical method was found to systematically describe intermolecular hydrogen bonding in small polar molecules. PM3 shows charge transfer from the donor to acceptor molecules on the order of 0.02-0.06 units of charge when strong hydrogen bonds are formed. The PM3 method is predictive; calculated hydrogen bond energies with an absolute magnitude greater than 2 kcal mol-' suggest that the global minimum is a hydrogen bonded complex; absolute energies less than 2 kcal mol-' imply that other van der Waals complexes are more stable. The geometries of the PM3 hydrogen bonded complexes agree with high-resolution spectroscopic observations, gas electron diffraction data, and high-level ab initio calculations. The main limitations in the PM3 method are the underestimation of hydrogen bond lengths by 0.1-0.2 for some systems and the underestimation of reliable experimental hydrogen bond energies by approximately 1-2 kcal mol-l. The PM3 method predicts that ammonia is a good hydrogen bond acceptor and a poor hydrogen donor when interacting with neutral molecules. Electronegativity differences between F, N, and 0 predict that donor strength follows the order F > 0 > N and acceptor strength follows the order N > 0 > F. In the calculations presented in this article, the PM3 method mirrors these electronegativity differences, predicting the F-H- - -N bond to be the strongest and the N-H- - -F bond the weakest. It appears that the PM3 Hamiltonian is able to model hydrogen bonding because of the reduction of two-center repulsive forces brought about by the parameterization of the Gaussian core-core interactions. The ability of the PM3 method to model intermolecular hydrogen bonding means reasonably accurate quantum-mechanical calculations can be applied to small biologic systems.

Measurement of carotenoid levels in human serum and a catalog of the lutein conformation populations from semi-empirical calculations

Lutein is a principal constituent of the human macular pigment. This study is composed of two projects. The first studies the conformational geometries of lutein and its potential adaptability in biological systems. The second is a study of the response of human subjects to lutein supplements. Using semi-empirical parametric method 3 (PM3) and density functional theory with the B3LYP/6-31G* basis set, the relative energies of s- cis conformers of lutein were determined. All 512 s-cis conformers were calculated with PM3. A smaller, representative group was also studied using density functional theory. PM3 results were correlated systematically to B3LYP values and this enables the results to be calibrated. The relative energies of the conformers range from 1-30 kcal/mole, and many are dynamically accessible at normal temperatures. Four commercial formulations containing lutein were studied. The serum and macular pigment (MP) responses of human subjects to these lutein supplements with doses of 9 or 20 mg/day were measured, relative to a placebo, over a six month period. In each instance, lutein levels in serum increased and correlated with MP increases. The results demonstrate that responses are significantly dependent upon formulation and that components other than lutein have an important influence serum response.

Core-valence Correlation Effects On Ir Calculations: The Bf3 And Bcl3 Cases.

The first theoretical results of core-valence correlation effects are presented for the infrared wavenumbers and intensities of the BF3 and BCl3 molecules, using (double- and triple-zeta) Dunning core-valence basis sets at the CCSD(T) level. The results are compared with those calculated in the frozen core approximation with standard Dunning basis sets at the same correlation level and with the experimental values. The general conclusion is that the effect of core-valence correlation is, for infrared wavenumbers and intensities, smaller than the effect of adding augmented diffuse functions to the basis set, e.g., cc-pVTZ to aug-cc-pVTZ. Moreover, the trends observed in the data are mainly related to the augmented functions rather than the core-valence functions added to the basis set. The results obtained here confirm previous studies pointing out the large descrepancy between the theoretical and experimental intensities of the stretching mode for BCl3.

Estudo teórico da interação existente entre a artemisinina e o heme

The purpose of this work is to study theoretically stereoelectronic aspects of the interaction between heme and artemisinin in the transitional heme-artemisinin complex. Through semi-empirical calculations using the PM3 method, the potential energy barrier of artemisinin rotation relative to heme in the heme-artemisinin complex was studied in vacuum and in the partially solvated state. The minimum heat of formation obtained for the complex with and without water molecules is -702.39 and -100.86 kcal mol-1, respectively, which corresponds to the dihedral angle C-Fe-O1-O2 of 43.93º and 51.90º around the iron-oxygen O1 bond, respectively. The water molecules bind to heme via 13 hydrogen bonds and O-H O and 6 C-H O interactions, which accounts for -67.23 kcal mol-1. It is observed that the inclusion of water molecules does not affect significantly the stability of the heme-artemisinin complex.

Edge effects in bilayer graphene nanoribbons: Ab initio total-energy density functional theory calculations

We show that the ground state of zigzag bilayer graphene nanoribbons is nonmagnetic. It also possesses a finite gap, which has a nonmonotonic dependence with the width as a consequence of the competition between bulk and strongly attractive edge interactions. All results were obtained using ab initio total-energy density functional theory calculations with the inclusion of parametrized van der Waals interactions.

Spin filtering and disorder-induced magnetoresistance in carbon nanotubes: Ab initio calculations

Nitrogen-doped carbon nanotubes can provide reactive sites on the porphyrin-like defects. It is well known that many porphyrins have transition-metal atoms, and we have explored transition-metal atoms bonded to those porphyrin-like defects inN-doped carbon nanotubes. The electronic structure and transport are analyzed by means of a combination of density functional theory and recursive Green's function methods. The results determined the heme B-like defect (an iron atom bonded to four nitrogens) is the most stable and has a higher polarization current for a single defect. With randomly positioned heme B defects in nanotubes a few hundred nanometers long, the polarization reaches near 100%, meaning they are effective spin filters. A disorder-induced magnetoresistance effect is also observed in those long nanotubes, and values as high as 20 000% are calculated with nonmagnectic eletrodes.

Barrier-free substitutional doping of graphene sheets with boron atoms: Ab initio calculations

Using ab initio methods, we propose a simple and effective way to substitutionally dope graphene sheets with boron. The method consists of selectively exposing each side of the graphene sheet to different elements. We first expose one side of the membrane to boron while the other side is exposed to nitrogen. Proceeding this way, the B atoms will be spontaneously incorporated into the graphene membrane without any activation barrier. In a second step, the system should be exposed to a H-rich environment, which will remove the CN radical from the layer and form HCN, leading to a perfect substitutional doping.

Optical absorption and electronic band structure first-principles calculations of alpha-glycine crystals

Light absorption of alpha-glycine crystals grown by slow evaporation at room temperature was measured, indicating a 5.11 +/- 0.02 eV energy band gap. Structural, electronic, and optical absorption properties of alpha-glycine crystals were obtained by first-principles quantum mechanical calculations using density functional theory within the generalized gradient approximation in order to understand this result. To take into account the contribution of core electrons, ultrasoft and norm-conserving pseudopotentials, as well as an all electron approach were considered to compute the electronic density of states and band structure of alpha-glycine crystals. They exhibit three indirect energy band gaps and one direct Gamma-Gamma energy gap around 4.95 eV. The optical absorption related to transitions between the top of the valence band and the bottom of the conduction band involves O 2p valence states and C, O 2p conduction states, with the carboxyl group contributing significantly to the origin of the energy band gap. The calculated optical absorption is highly dependent on the polarization of the incident radiation due to the spatial arrangement of the dipolar glycine molecules; in the case of a polycrystalline sample, the first-principles calculated optical absorption is in good agreement with the measurement when a rigid energy shift is applied.

Accurate band gaps of AlGaN, InGaN, and AlInN alloys calculations based on LDA-1/2 approach

We present parameter-free calculations of electronic properties of InGaN, InAlN, and AlGaN alloys. The calculations are based on a generalized quasichemical approach, to account for disorder and composition effects, and first-principles calculations within the density functional theory with the LDA-1/2 approach, to accurately determine the band gaps. We provide precise results for AlGaN, InGaN, and AlInN band gaps for the entire range of compositions, and their respective bowing parameters. (C) 2011 American Institute of Physics. [doi:10.1063/1.3576570]

Structure, stability, depolarized light scattering, and vibrational spectra of fullerenols from all-electron density-functional-theory calculations

We have investigated the stability, electronic properties, Rayleigh (elastic), and Raman (inelastic) depolarization ratios, infrared and Raman absorption vibrational spectra of fullerenols [C(60)(OH)(n)] with different degrees of hydroxylation by using all-electron density-functional-theory (DFT) methods. Stable arrangements of these molecules were found by means of full geometry optimizations using Becke's three-parameter exchange functional with the Lee, Yang, and Parr correlation functional. This DFT level has been combined with the 6-31G(d,p) Gaussian-type basis set, as a compromise between accuracy and capability to treat highly hydroxylated fullerenes, e.g., C(60)(OH)(36). Thus, the molecular properties of fullerenols were systematically analyzed for structures with n=1, 2, 3, 4, 8, 10, 16, 18, 24, 32, and 36. From the electronic structure analysis of these molecules, we have evidenced an important effect related to the weak chemical reactivity of a possible C(60)(OH)(24) isomer. To investigate Raman scattering and the vibrational spectra of the different fullerenols, frequency calculations are carried out within the harmonic approximation. In this case a systematic study is only performed for n=1-4, 8, 10, 16, 18, and 24. Our results give good agreements with the expected changes in the spectral absorptions due to the hydroxylation of fullerenes.