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.

A microscopic view of substitution reactions solvated by ionic liquids

The solvation effect of the ionic liquid 1-N-butyl-3-methylimidazolium hexafluorophosphate on nucleophilic substitution reactions of halides toward the aliphatic carbon of methyl p-nitrobenzenesulfonate (pNBS) was investigated by computer simulations. The calculations were performed by using a hybrid quantum-mechanical/molecular-mechanical (QM/MM) methodology. A semiempirical Hamiltonian was first parametrized on the basis of comparison with ab initio calculations for Cl(-) and Br(-) reaction with pNBS at gas phase. In condensed phase, free energy profiles were obtained for both reactions. The calculated reaction barriers are in agreement with experiment. The structure of species solvated by the ionic liquid was followed along the reaction progress from the reagents, through the transition state, to the final products. The simulations indicate that this substitution reaction in the ionic liquid is slower than in nonpolar molecular solvents proper to significant stabilization of the halide anion by the ionic liquid in comparison with the transition state with delocalized charge. Solute-solvent interactions in the first solvation shell contain several hydrogen bonds that are formed or broken in response to charge density variation along the reaction coordinate. The detailed structural analysis can be used to rationalize the design of new ionic liquids with tailored solvation properties. (c) 2008 American Institute of Physics.

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]

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.

Solvent Effects in Optical Spectra of ortho-Aminobenzoic Acid Derivatives

We investigated three amino derivatives of ortho-aminobenzoic or anthranilic acid (o-Abz): a) 2-Amino-benzamide (AbzNH(2)); b) 2-Amino-N-methyl-benzamide (AbzNHCH(3)) and c) 2-Amino-N-N`-dimethyl-bezamide (AbzNH(CH(3))(2)), see Scheme 1. We describe the results of ab-initio calculations on the structural characteristics of the compounds and experimental studies about solvent effects in their absorption and steady-state and time-resolved emission properties. Ab-initio calculations showed higher stability for the rotameric conformation in which the oxygen of carbonyl is near to the nitrogen of ortho-amino group. The derivatives present decrease in the delocalization of pi electron, and absorption bands are blue shifted compared to the parent compound absorption, the extent of the effect increasing from to Abz-NH(2) to Abz-NHCH(3) Abz-NH(CH(3))(2). Measurements performed in several solvents have shown that the the dependence of Stokes shift of the derivatives with the orientational polarizability follows the Onsager-Lippert model for general effects of solvent. However deviation occurred in solvents with properties of Bronsted acids, or electron acceptor characteristics, so that hydrogen bonds formed with protic solvents predominates over intramolecular hydrogen bond. In most solvents the fluorescence decay of AbzNH(2) and AbzNHCH(3) was fitted to a single exponential with lifetimes around 7.0 ns and no correlation with polarity of the solvent was observed. The fluorescence decay of AbzN(CH(3))(2) showed lifetimes around 2.0 ns, consistent with low quantum yield of the compound. The spectroscopic properties of the monoamino derivative AbzNHCH(3) are representative of the properties presented by Abz labelled peptides and fatty acids previously studied.

Local character of magnetic coupling in ionic solids

Magnetic interactions in ionic solids are studied using parameter-free methods designed to provide accurate energy differences associated with quantum states defining the Heisenberg constant J. For a series of ionic solids including KNiF3, K2NiF4, KCuF3, K2CuF4, and high- Tc parent compound La2CuO4, the J experimental value is quantitatively reproduced. This result has fundamental implications because J values have been calculated from a finite cluster model whereas experiments refer to infinite solids. The present study permits us to firmly establish that in these wide-gap insulators, J is determined from strongly local electronic interactions involving two magnetic centers only thus providing an ab initio support to commonly used model Hamiltonians.

Electronic structure of a neutral oxygen vacancy in SrTiO3

The electronic structure of an isolated oxygen vacancy in SrTiO3 has been investigated with a variety of ab initio quantum mechanical approaches. In particular we compared pure density functional theory (DFT) approaches with the Hartree-Fock method, and with hybrid methods where the exchange term is treated in a mixed way. Both local cluster models and periodic calculations with large supercells containing up to 80 atoms have been performed. Both diamagnetic (singlet state) and paramagnetic (triplet state) solutions have been considered. We found that the formation of an O vacancy is accompanied by the transfer of two electrons to the 3d(z2) orbitals of the two Ti atoms along the Ti-Vac-Ti axis. The two electrons are spin coupled and the ground state is diamagnetic. New states associated with the defect center appear in the gap just below the conduction band edge. The formation energy computed with respect to an isolated oxygen atom in the triplet state is 9.4 eV.

Role of structural saturation and geometry in the luminescence of silicon-based nanostructured materials

The structural saturation and stability, the energy gap, and the density of states of a series of small, silicon-based clusters have been studied by means of the PM3 and some ab initio (HF/6-31G* and 6-311++G**, CIS/6-31G* and MP2/6-31G*) calculations. It is shown that in order to maintain a stable nanometric and tetrahedral silicon crystallite and remove the gap states, the saturation atom or species such as H, F, Cl, OH, O, or N is necessary, and that both the cluster size and the surface species affect the energetic distribution of the density of states. This research suggests that the visible luminescence in the silicon-based nanostructured material essentially arises from the nanometric and crystalline silicon domains but is affected and protected by the surface species, and we have thus linked most of the proposed mechanisms of luminescence for the porous silicon, e.g., the quantum confinement effect due to the cluster size and the effect of Si-based surface complexes.

Termodinâmica estatística de líquidos com o método de Monte Carlo. II. O líquido clorofórmio

Thermodynamic properties and radial distribution functions for liquid chloroform were calculated using the Monte Carlo method implemented with Metropolis algorithm in the NpT ensemble at 298 K and 1 atm. A five site model was developed to represent the chloroform molecules. A force field composed by Lennard-Jones and Coulomb potential functions was used to calculate the intermolecular energy. The partial charges needed to represent the Coulombic interactions were obtained from quantum chemical ab initio calculations. The Lennard-Jones parameters were adjusted to reproduce experimental values for density and enthalpy of vaporization for pure liquid. All thermodynamic results are in excelent agreement with experimental data. The correlation functions calculated are in good accordance with theoretical results avaliable in the literature. The free energy for solvating one chloroform molecule into its own liquid at 298 K and 1 atm was computed as an additional test of the potential model. The result obtained compares well with the experimental value. The medium effects on cis/trans convertion of a hypotetical solute in water TIP4P and chloroform solvents were also accomplished. The results obtained from this investigation are in agreement with estimates of the continuous theory of solvation.

Técnicas para análise de desempenho de computadores

A series of open source benchmarks for computer performance analysis of personal computers with a focus on computational chemistry calculations is presented. The results returned by these tests are discussed and used to correlate with the actual performance of a set of computers available for research on two computing intensive fields of chemistry, quantum chemical and molecular simulation calculations.

Estudo QSPR sobre os coeficientes de partição: descritores mecânico-quânticos e análise multivariada

Quantum chemistry and multivariate analysis were used to estimate the partition coefficients between n-octanol and water for a serie of 188 compounds, with the values of the q 2 until 0.86 for crossvalidation test. The quantum-mechanical descriptors are obtained with ab initio calculation, using the solvation effects of the Polarizable Continuum Method. Two different Hartree-Fock bases were used, and two different ways for simulating solvent cavity formation. The results for each of the cases were analised, and each methodology proposed is indicated for particular case.

Análise estrutural de ciclodextrinas: um estudo comparativo entre métodos teóricos clássicos e quânticos

In the present work, we analyzed the accuracy of distinct theoretical methods to reproduce the solid state structures of cyclodextrins. The a, b and g-cyclodextrins (CD) were considered and also their hydrates with included water molecules: a-CD.2H2O, b-CD.10H2O and g-CD.12H2O. The geometries were fully optimized using Molecular Mechanics (MM2), semiempirical (AM1 and PM3) and ab initio (HF/3-21G) methods and quantitatively compared with experimental data from X ray diffraction. The results obtained from the classical MM2 method were in best agreement with the experiment. The semiempirical and ab initio structures were also in satisfactory accordance with the experimental data. In general, the PM3 method was found to be more suitable than the AM1 to describe the CD geometries, mainly when the intramolecular hydrogen bonds are considered.

Estrutura e propriedades de elipticinas

The ellipticines constitute a broad class of molecules with antitumor activity. In the present work we analyzed the structure and properties of a series of ellipticine derivatives in the gas phase and in solution using quantum mechanical and Monte Carlo methods. The results showed a good correlation between the solvation energies in water obtained with the continuum model and the Monte Carlo simulation. Molecular descriptors were considered in the development of QSAR models using the DNA association constant (log Kapp) as biological data. The results showed that the DNA binding is dominated by electronic parameters, with small contributions from the molecular volume and area.

Implementação computacional do modelo carga-fluxo de carga-fluxo de dipolo para cálculo e interpretação das intensidades do espectro infravermelho

The first computational implementation that automates the procedures involved in the calculation of infrared intensities using the charge-charge flux-dipole flux model is presented. The atomic charges and dipoles from the Quantum Theory of Atoms in Molecules (QTAIM) model was programmed for Morphy98, Gaussian98 and Gaussian03 programs outputs, but for the ChelpG parameters only the Gaussian programs are supported. Results of illustrative but new calculations for the water, ammonia and methane molecules at the MP2/6-311++G(3d,3p) theoretical level, using the ChelpG and QTAIM/Morphy charges and dipoles are presented. These results showed excellent agreement with analytical results obtained directly at the MP2/6-311++G(3d,3p) level of theory.