Cellulose pyrolysis and quantum chemistry

Cellulose is the major constituent of most plants of interest as renewable sources of energy and is the most extensively studied form of biomass or biomass constituent. Predicting the mass loss and product yields when cellulose is subjected to increased temperature represents a fundamental problem in the thermal release of biomass energy. Unfortunately, at this time, there is no internally consistent model of cellulose pyrolysis that can organize the varied experimental data now available or provide a guide for additional experiments. Here, we present a model of direct cellulose pyrolysis using a multistage decay scheme that we first presented in the IJQC in 1984. This decay scheme can, with the help of an inverse method of assigning reaction rates, provide a reasonable account of the direct fast pyrolysis yield measurements. The model is suggestive of dissociation states of d-glucose (C6H10O5,), the fundamental cellulose monomer. The model raises the question as to whether quantum chemistry could now provide the dissociation energies for the principal breakup modes of glucose into C-1, C-2, C-3, C-4, and C-5 compounds. These calculations would help in achieving a more fundamental description of volatile generation from cellulose pyrolysis and could serve as a guide for treating hemicellulose and lignin, the other major biomass constituents. Such advances could lead to the development of a predictive science of biomass pyrolysis that would facilitate the design of liquifiers and gasifiers based upon renewable feedstocks. (C) 1998 John Wiley & Sons, Inc.

Quantum gauge boson propagators in the light front

Gauge fields in the light front are traditionally addressed via, the employment of an algebraic condition n·A = 0 in the Lagrangian density, where Aμ is the gauge field (Abelian or non-Abelian) and nμ is the external, light-like, constant vector which defines the gauge proper. However, this condition though necessary is not sufficient to fix the gauge completely; there still remains a residual gauge freedom that must be addressed appropriately. To do this, we need to define the condition (n·A) (∂·A) = 0 with n·A = 0 = ∂·A. The implementation of this condition in the theory gives rise to a gauge boson propagator (in momentum space) leading to conspicuous nonlocal singularities of the type (k·n)-α where α = 1, 2. These singularities must be conveniently treated, and by convenient we mean not only mathemathically well-defined but physically sound and meaningful as well. In calculating such a propagator for one and two noncovariant gauge bosons those singularities demand from the outset the use of a prescription such as the Mandelstam-Leibbrandt (ML) one. We show that the implementation of the ML prescription does not remove certain pathologies associated with zero modes. However we present a causal, singularity-softening prescription and show how to keep causality from being broken without the zero mode nuisance and letting only the propagation of physical degrees of freedom.

Quantum confinement in hydrogen bond

In this work, the quantum confinement effect is proposed as the cause of the displacement of the vibrational spectrum of molecular groups that involve hydrogen bonds. In this approach, the hydrogen bond imposes a space barrier to hydrogen and constrains its oscillatory motion. We studied the vibrational transitions through the Morse potential, for the NH and OH molecular groups inside macromolecules in situation of confinement (when hydrogen bonding is formed) and nonconfinement (when there is no hydrogen bonding). The energies were obtained through the variational method with the trial wave functions obtained from supersymmetric quantum mechanics formalism. The results indicate that it is possible to distinguish the emission peaks related to the existence of the hydrogen bonds. These analytical results were satisfactorily compared with experimental results obtained from infrared spectroscopy. (c) 2015 Wiley Periodicals, Inc.

Distribution of neutral prilocaine in a phospholipid bilayer: Insights from molecular dynamics simulations

In this work, we report a 20-ns constant pressure molecular dynamics simulation of prilocaine (PLC), in amine-amide local anesthetic, in a hydrated liquid crystal bilayer of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine. The partition of PLC induces the lateral expansion of the bilayer and a concomitant contraction in its thickness. PLC molecules are preferentially found in the hydrophobic acyl chains region, with a maximum probability at similar to 12 angstrom from the center of the bilayer (between the C(4) and C(5) methylene groups). A decrease in the acyl chain segmental order parameter, vertical bar S-CD vertical bar, compared to neat bilayers, is found, in good agreement with experimental H-2-NMR studies. The decrease in vertical bar S-CD vertical bar induced by PLC is attributed to a larger accessible volume per lipid in the acyl chain region. (C) 2008 Wiley Periodicals, Inc.

Dynamics of meso and thermo citrate synthases with implicit solvation

The dynamics of hydration of meso and thermo citrate synthases has been investigated using the EEF1 methodology implemented with the CHARNM program. The native enzymes are composed of two identical subunits, each divided into a small and large domain. The dynamics behavior of both enzymes at 30 degrees C and 60 degrees C has been compared. The results of simulations show that during the hydration process, each subunit follows a different pathway of hydration, in spite of the identical sequence. The hydrated structures were compared with the crystalline structure, and the root mean square deviation (RMSD) of each residue along the trajectory was calculated. The regions with larger and smaller mobility were identified. In particular, helices belonging to the small domain are more mobile than those of the large domain. In contrast, the residues that constitute the active site show a much lower displacement compared with the crystalline structure. Hydration free energy calculations point out that Thermoplasma acidophilum citrate synthase (TCS) is more stable than chicken citrate synthase (CCS), at high temperatures. Such result has been ascribed to the higher number of superficial charges in the thermophilic homologue, which stabilizes the enzyme, while the mesophilic homologue denatures. These results are in accord with the experimental found that TCS keeps activity at temperatures farther apart from the catalysis regular temperature than the CCS. (c) 2005 Wiley Periodicals, Inc.

Computer-aided design of a novel ligand for retinoic acid receptor in cancer chemotherapy

The isotypes of RAR and RXR are retinoic acid and retinoid X acid receptors, respectively, whose ligand-binding domain contains the ligand-dependent activation function, with distinct pharmacological targets for retinoids, involved in the treatment of various cancers and skin diseases. Due to the major challenge which cancer treatment and cure still imposes after many decades to the international scientific community, there is actually considerable interest in new ligands with increased bioactivity. We have focused on the retinoid acid receptor, which is considered an interesting target for drug design. In this work, we carried out density functional geometry optimizations, and different docking procedures. We performed screening in a large database (hundreds of thousands of molecules which we optimized at the AM1 level) yielding a set of potential bioactive ligands. A new ligand was selected and optimized at the B3LYP/6-31G* level. A flexible docking program was used to investigate the interactions between the receptor and the new ligand. The result of this work is compared with several crystallographic ligands of RAR. Our theoretically more bioactive new-ligand indicates stronger and more hydrogen bonds as well as hydrophobic interactions with the receptor. (c) 2005 Wiley Periodicals, Inc.

Theoretical study of omeprazole behavior: Racemization barrier and decomposition reaction

omeprazole is a substituted benzimidazole which suppresses gastric-acid secretion by means of H+, K+-ATPase inhibition. It is an optically active drug with the sulfur of the sulfoxide being the chiral center. This pro-drug can be easily converted into its respective sulfenamide at low pH. In this work, omeprazole has been studied in relation to racemization barrier and decomposition reaction. Quantum chemistry coupled to PCA chemometric method were used to find all minimum energy structures. Conformational analysis and calculation of racemization barriers were carried out by PM3 semiempirical method (Gaussian 98). The average racemization energy barrier for all minimum energy structures (43.56 kcal mol(-1)) can be related to the velocity constant in Eyring's equation. The enormous half-life time at 100 degrees C (9.04 x 10(4) years) indicates that the process cannot be observed in human time scale. on the other hand, the difference of free energy change (Delta(Delta G) = -266.78 kcal mol(-1)) for the decomposition reaction shows that the process is favorable to the sulfenamide formation. The highly negative Delta(Delta G) obtained for the decomposition reaction shows that this process is extremely exothermic. This result explains why omeprazole decomposes and does not racemize. (C) 2008 Wiley Periodicals, Inc.

Dimensional Compactification and Two-Particle Binding

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Basicities of primary arylamines and calculated amine nitrogen electronic charges

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

An ab initio study of oxygen vacancies and doping process of Nb and Cr atoms on TiO2 (110) surface models

We theoretically investigated how the formation of oxygen vacancies and the addition of niobium and chromium atoms as dopants modify the varistor properties of TiO2. The calculations were carried out at the HF level using a contracted basis set, developed by Huzinaga et al.. to represent the atomic centers on the (110) surface for the large (TiO2)(15) cluster model. The change of the values for the net atomic charges and band gap after oxygen vacancy formation and the presence of dopants in the lattice are analyzed and discussed. It is shown that the formation of oxygen vacancies decreases the band gap while an opposite effect is found when dopants are located in the reduced surface. The theoretical results are compared with available experimental data. A plausible explanation of the varistor behavior of this system is proposed. (C) 1997 John Wiley & Sons, Inc.

Relation between antioxidant activity and electronic structure of Phenols

We describe a new physicochemical descriptor of the antioxidant activity of phenols, the energy difference between the two highest occupied molecular orbitals, which we believe will improve quantitative structure-activity relationship studies about these compounds. (C) 2003 Wiley Periodicals, Inc.

Theoretical study of MgO(001) surfaces: Pure, doped with Fe, Ca, and Al, and with and without adsorbed water

Ab initio calculations of large cluster models have been performed in order to study water adsorption at the five-fold coordinated adsorption site on pure Mg(001) and MgO(001) surfaces doped with Fe, Ca, and Al. The geometric parameters of the adsorbed water molecule have been optimized preparatory to analysis of binding energies, charge transfer, preferential sites of interaction, and bonding distances. We have used Mulliken population analysis methods in order to analyze charge distributions and the direction of charge transfer. We have also investigated energy gaps, HOMO energies, and SCF orbital energies as well as the acid-base properties of our cluster model. Numerical results are compared, where possible, with experiment and interpreted in the framework of various analytical models. (C) 2001 John Wiley & Sons, Inc.

Time-dependent non-Hermitian Hamiltonians with real energies

In this work we intend to study a class of time-dependent quantum systems with non-Hermitian Hamiltonians, particularly those whose Hermitian counterparts are important for the comprehension of posed problems in quantum optics and quantum chemistry. They consist of an oscillator with time-dependent mass and frequency under the action of a time-dependent imaginary potential. The wave functions are used to obtain the expectation value of the Hamiltonian. Although it is neither Hermitian nor PT symmetric, the Hamiltonian under study exhibits real values of energy.

C-H center dot center dot center dot O and N-H center dot center dot center dot O hydrogen bonds in liquid amides investigated by Monte Carlo simulation

Monte Carlo simulations of liquid formamide, N-methylformamide (MF), and N,N-dimethylformamide (DMF) have been performed in the isothermal and isobaric ensemble at 298 K and 1 atm, aiming to investigate the C-H ... O and N-H ... O hydrogen bonds. The interaction energy was calculated using the classical 6-12 Lennard-Jones pairwise potential plus a Coulomb term on a rigid six-site molecular model with the potential parameters being optimized in this work. Theoretical values obtained for heat of vaporization and liquid densities are in good agreement with the experimental data. The radial distribution function [RDF, g(r)] obtained compare well with R-X diffraction data available. The RDF and molecular mechanics (MM2) minimization show that the C-H ... O interaction has a significant role in the structure of the three liquids. These results are supported by ab initio calculations. This Interaction is particularly important in the structure of MF. The intensity of the N-H ... O hydrogen bond is greater in the MF than formamide. This could explain some anomalous properties verified in MF. (C) 1997 John Wiley & Sons, Inc.

Theoretical analysis on TiO(2)(110)/V surface

Theoretical analysis based on the Hartree-Fock method were performed in order to study the stoichiometric TiO(2) (110) surface and the vanadium substituted system. The Pople with polarization 3-21G* basis set level was used. The TiO(2) (110) surface was modeled using a (TiO(2))(15) cluster model. In order to take into account the finite size of the cluster, we have studied two different models: the point charge and the hydrogen saturated methodologies. The charge values used in the point charge calculations were optimized. The density of states, orbital self-consistend field (SCF) energies, and Mulliken charge values were analyzed. The method and model's dependence on the analyzed results are discussed. The theoretical results are compared with available experimental data. (C) 2001 John Wiley & Sons, Inc.