Técnicas de análise do perfil de execução e otimização de programas em química computacional

In this paper we review the basic techniques of performance analysis within the UNIX environment that are relevant in computational chemistry, with particular emphasis on the execution profile using the gprof tool. Two case studies (in ab initio and molecular dynamics calculations) are presented in order to illustrate how execution profiling can be used to effectively identify bottlenecks and to guide source code optimization. Using these profiling and optimization techniques it was possible to obtain significant speedups (of up to 30%) in both cases.

Análise teórica da interação de CO, CO2 e NH3 com ZnO

This paper presents a study of the interaction of small molecules with ZnO surfaces by means of theoretical methods. The AM1 semi-empirical method was used for optimizing the geometric parameters of adsorbed molecules. The optimized AM1 structures were used in the calculations of the ab initio RHF method with the 3-21G* basis set. The interaction of CO, CO2 and NH3 molecules were studied with (ZnO)22 and (ZnO)60 cluster models. We have analyzed the interaction energy, SCF orbital energies, Mulliken charges and the density of states (DOS).

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.

Polarizabilidades e primeira hiperpolarizabilidades elétricas dipolares de ésteres E-4-amino-trans-1,3-butadienil-(1',2'-di-hidroxibenzeno) de boro, alumínio e gálio

Static electric dipole polarizabilities and first hyperpolarizabilites have been calculated for the title molecules and their 3' and 4'-nitro derivatives at ab-initio Hartree- Fock/6-31G(d, p) level. The influence of the pivotal p vacant 3A elements (B, Al or Ga) substitution on the electrical properties of these molecules is detailed. The axial vector components of the first hyperpolarizabilities β(0) of the push-pull 4'-nitro derivatives, -18.2×10-32 esu (B), -21.1×10-32 esu (Al) and -20.8×10-32 esu (Ga) are calculated to be as much as fourfold larger then that calculated for the p-nitroaniline, a reference organic molecule for comparison for this type of molecular property.

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.

Estudo das ligações de hidrogênio para dímeros formados pelas moléculas de H2O, NH3, HF, HCl e HBr através de cálculos baseados em primeiros princípios

Hydrogen bond energies of fifteen dimers were calculated using the large basis set 6-311++G(3df,3pd), at Hartree-Fock (HF) level including Møller-Plesset (MP2) calculations. The procedure for obtaining such energies were based on the dimer's energy rise provoked by increasing in intermolecular distance of the system component units. Deviations from a strictly linear hydrogen bond were investigated and rotational barriers were also computed allowing the calculation of the second order attractive interactions. In order to provide a more objective definition of hydrogen bond, a lower energy limit was proposed in place of the merely empirical parameters employed in the classical definition

First principles modeling of metallic alloys and alloy surfaces

By alloying metals with other materials, one can modify the metal’s characteristics or compose an alloy which has certain desired characteristics that no pure metal has. The field is vast and complex, and phenomena that govern the behaviour of alloys are numerous. Theories cannot penetrate such complexity, and the scope of experiments is also limited. This is why the relatively new field of ab initio computational methods has much to give to this field. With these methods, one can extend the understanding given by theories, predict how some systems might behave, and be able to obtain information that is not there to see in physical experiments. This thesis pursues to contribute to the collective knowledge of this field in the light of two cases. The first part examines the oxidation of Ag/Cu, namely, the adsorption dynamics and oxygen induced segregation of the surface. Our results demonstrate that the presence of Ag on the Cu(100) surface layer strongly inhibits dissociative adsorption. Our results also confirmed that surface reconstruction does happen, as experiments predicted. Our studies indicate that 0.25 ML of oxygen is enough for Ag to diffuse towards the bulk, under the copper oxide layer. The other part elucidates the complex interplay of various energy and entropy contributions to the phase stability of paramagnetic duplex steel alloys. We were able to produce a phase stability map from first principles, and it agrees with experiments rather well. Our results also show that entropy contributions play a very important role on defining the phase stability. This is, to the author’s knowledge, the first ab initio study upon this subject.

Defects in Persistent Luminescence Materials

Persistent luminescence materials can store energy from solar radiation or artificial lighting and release it over a period of several hours without a continuous excitation source. These materials are widely used to improve human safety in emergency and traffic signalization. They can also be utilized in novel applications including solar cells, medical diagnostics, radiation detectors and structural damage sensors. The development of these materials is currently based on methods based on trial and error. The tailoring of new materials is also hindered by the lack of knowledge on the role of their intrinsic and extrinsic lattice defects in the appropriate mechanisms. The goal of this work was to clarify the persistent luminescence mechanisms by combining ab initio density functional theory (DFT) calculations with selected experimental methods. The DFT approach enables a full control of both the nature of the defects and their locations in the host lattice. The materials studied in the present work, the distrontium magnesium disilicate (Sr₂MgSi₂O₇) and strontium aluminate (SrAl₂O₄) are among the most efficient persistent luminescence hosts when doped with divalent europium Eu²⁺ and co-doped with trivalent rare earth ions R³⁺ (R: Y, La-Nd, Sm, Gd-Lu). The polycrystalline materials were prepared with the solid state method and their structural and phase purity was confirmed by X-ray powder diffraction. Their local crystal structure was studied by high-resolution transmission electron microscopy. The crystal and electronic structure of the nondoped as well as Eu²⁺, R^2+/3+ and other defect containing materials were studied using DFT calculations. The experimental trap depths were obtained using thermoluminescence (TL) spectroscopy. The emission and excitation of Sr₂MgSi₂O₇:Eu²+,Dy³⁺ were also studied. Significant modifications in the local crystal structure due to the Eu²⁺ ion and lattice defects were found by the experimental and DFT methods. The charge compensation effects induced by the R³⁺ co-doping further increased the number of defects and distortions in the host lattice. As for the electronic structure of Sr₂MgSi₂O₇ and SrAl₂O₄, the experimental band gap energy of the host materials was well reproduced by the calculations. The DFT calculated Eu²⁺ and R^2+/3+ 4fn as well as 4f^n-15d¹ ground states in the Sr₂MgSi₂O₇ band structure provide an independent verification for an empirical model which is constructed using rather sparse experimental data for the R³⁺ and especially the R²⁺ ions. The intrinsic and defect induced electron traps were found to act together as energy storage sites contributing to the materials’ efficient persistent luminescence. The calculated trap energy range agreed with the trap structure of Sr₂MgSi₂O₇ obtained using TL measurements. More experimental studies should be carried out for SrAl₂O₄ to compare with the DFT calculations. The calculated and experimental results show that the electron traps created by both the rare earth ions and vacancies are modified due to the defect aggregation and charge compensation effects. The relationships between this modification and the energy storage properties of the solid state materials are discussed.

Optical spectroscopy of organic material based on C60

The structure and optical properties of thin films based on C60 materials are studied. Reproducible vacuum method of thin fullerene films production with Cd impurity on Si, glass and mica surfaces is developed. Surface morphology of the films are investigated by AFM and SEM methods. The ab initio quantum - chemical calculations of the geometry, total energy and excited energy states of complex fullerene- cadmium telluride supramolecules are performed. Photoluminescence spectra of composite thin films based on C60 before and after X-ray irradiation were measured. The intensity of additional peaks is defined as the charge composition due to the type of substrate. These results are interpreted as an appearance of the dipole-allowed transitions in the fullerene excited singlet states spectrum cause of an interference with cadmium telluride. X-ray irradiated films were investigated, and additional peaks in photoluminescence spectra were detected. These peaks appear as a result of molecular complexes formation from C60CdTe mixture and dimerization of the films. Density functional B3LYP quantum-chemical calculations for C60CdTe, molecular complexes, (C60)2 and C120O dimers were performed to elucidate some experimental results.

Mass spectrometric studies on aryltin compounds and alkali halides /

The fragmentation behavior of aryltin compounds [(p-ThAnis)nSnPh4.n (n=l-4); (p-ThAnis)3SnX (X=C1, Br, I); (o-CH30C6H4)3SnCl; Ph3Sn(o-pyr)] have been studied comparatively under EI and FAB ionization modes. Alkali halides were run under FAB mode. For the aryltin compounds, the effect of ligand type on the spectra have been explored in both EI and FAB modes. The fragmentation mechanisms have been examined with linked scans, such as fragment ion scans (B/E) and parent ion scans (B^/E). Ab Initio molecular orbital calculations were used to determine the structures of the fragments by comparing their relative stabilities. In the EI MS studies, negative ion EI mode has also been used for some of the aryltin compounds, to examine the possible ion molecule reactions under low pressures at 70eV. In the positive ion FAB MS studies, matrix optimization experiments have been carried out. Negative ion FAB experiments of all the compounds have been done in two different ways. Finally, the comparison of the two methods, EI MS and FAB MS, have been made.For alkali halides, the studies focused on the FAB MS behavior under different conditions. The intensities of cluster ions were reported, and the anomalies in the intensity distribution was also discussed.

An analysis of the methyl rotation dynamics in the So (x' A) and T (a A) states of thioacetaldehyde from Pyrolysis jet spectra /

Jet-cooled, laser-induced phosphorescence excitation spectra (LIP) of thioacetaldehyde CH3CHS, CH3CDS, CD3CHS and CD3CDS have been observed over the region 15800 - 17300 cm"^ in a continuous pyrolysis jet. The vibronic band structure of the singlet-triplet n -* n* transition were attributed to the strong coupling of the methyl torsion and aldehydic hydrogen wagging modes . The vibronic peaks have been assigned in terms of two upper electronic state (T^) vibrations; the methyl torsion mode v^g, and the aldehydic hydrogen wagging mode v^^. The electronic origin O^a^ is unequivocally assigned as follows: CH3CHS (16294.9 cm"'' ), CH3CDS (16360.9 cm"'' ), CD3CHS (16299.7 cm"^ ), and CD3CDS (16367.2 cm"'' ). To obtain structural and dynamical information about the two electronic states, potential surfaces V(e,a) for the 6 (methyl torsion) and a (hydrogen wagging) motions were generated by ab initio quantum mechanical calculations with a 6-3 IG* basis in which the structural parameters were fully relaxed. The kinetic energy coefficients BQ(a,e) , B^(a,G) , and the cross coupling term B^(a,e) , were accurately represented as functions of the two active coordinates, a and 9. The calculations reveal that the molecule adopts an eclipsed conformation for the lower Sq electronic state (a=0°,e=0"') with a barrier height to internal rotation of 541.5 cm"^ which is to be compared to 549.8 cm"^ obtained from the microwave experiment. The conformation of the upper T^ electronic state was found to be staggered (a=24 . 68° ,e=-45. 66° ) . The saddle point in the path traced out by the aldehyde wagging motion was calculated to be 175 cm"^ above the equilibrium configuration. The corresponding maxima in the path taken by methyl torsion was found to be 322 cm'\ The small amplitude normal vibrational modes were also calculated to aid in the assignment of the spectra. Torsional-wagging energy manifolds for the two states were derived from the Hamiltonian H(a,e) which was solved variationally using an extended two dimensional Fourier expansion as a basis set. A torsionalinversion band spectrum was derived from the calculated energy levels and Franck-Condon factors, and was compared with the experimental supersonic-jet spectra. Most of the anomalies which were associated with the interpretation of the observed spectrum could be accounted for by the band profiles derived from ab initio SCF calculations. A model describing the jet spectra was derived by scaling the ab initio potential functions. The global least squares fitting generates a triplet state potential which has a minimum at (a=22.38° ,e=-41.08°) . The flatter potential in the scaled model yielded excellent agreement between the observed and calculated frequency intervals.

Etude théorique des propriétés électroniques et optiques des super-réseaux de Si/SiO2

FRANCAIS: L'observation d'une intense luminescence dans les super-réseaux de Si/SiO2 a ouvert de nouvelles avenues en recherche théorique des matériaux à base de silicium, pour des applications éventuelles en optoélectronique. Le silicium dans sa phase cristalline possède un gap indirect, le rendant ainsi moins intéressant vis-à-vis d'autres matériaux luminescents. Concevoir des matériaux luminescents à base de silicium ouvrira donc la voie sur de multiples applications. Ce travail fait état de trois contributions au domaine. Premièrement, différents modèles de super-réseaux de Si/SiO2 ont été conçus et étudiés à l'aide de calculs ab initio afin d'en évaluer les propriétés structurales, électroniques et optiques. Les deux premiers modèles dérivés des structures cristallines du silicium et du dioxyde de silicium ont permis de démontrer l'importance du rôle de l'interface Si/SiO2 sur les propriétés optiques. De nouveaux modèles structurellement relaxés ont alors été construits afin de mieux caractériser les interfaces et ainsi mieux évaluer la portée du confinement sur les propriétés optiques. Deuxièmement, un gap direct dans les modèles structurellement relaxés a été obtenu. Le calcul de l'absorption (par l'application de la règle d'or de Fermi) a permis de confirmer que les propriétés d'absorption (et d'émission) du silicium cristallin sont améliorées lorsque celui-ci est confiné par le SiO2. Un décalage vers le bleu avec accroissement du confinement a aussi été observé. Une étude détaillée du rôle des atomes sous-oxydés aux interfaces a de plus été menée. Ces atomes ont le double effet d'accroître légèrement le gap d'énergie et d'aplanir la structure électronique près du niveau de Fermi. Troisièmement, une application directe de la théorique des transitions de Slater, une approche issue de la théorie de la fonctionnelle de la densité pour des ensembles, a été déterminée pour le silicium cristallin puis comparée aux mesures d'absorption par rayons X. Une très bonne correspondance entre cette théorie et l'expérience est observée. Ces calculs ont été appliqués aux super-réseaux afin d'estimer et caractériser leurs propriétés électroniques dans la zone de confinement, dans les bandes de conduction.

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.

Extent and limitations of density functional theory in describing magnetic systems

The performance of density-functional theory to solve the exact, nonrelativistic, many-electron problem for magnetic systems has been explored in a new implementation imposing space and spin symmetry constraints, as in ab initio wave function theory. Calculations on selected systems representative of organic diradicals, molecular magnets and antiferromagnetic solids carried out with and without these constraints lead to contradictory results, which provide numerical illustration on this usually obviated problem. It is concluded that the present exchange-correlation functionals provide reasonable numerical results although for the wrong physical reasons, thus evidencing the need for continued search for more accurate expressions.

Magnetic coupling in the weac ferromagnet CuF2

CuF2 is known to be an antiferromagnetic compound with a weak ferromagnetism due to the anisotropy of its monoclinic unit cell (Dzialoshinsky-Moriya mechanism). We investigate the magnetic ordering of this compound by means of ab initio periodic unrestricted Hartree-Fock calculations and by cluster calculations which employ state-of-the-art configuration interaction expansions and modern density functional theory techniques. The combined use of periodic and cluster models permits us to firmly establish that the antiferromagnetic order arises from the coupling of one-dimensional subunits which themselves exhibit a very small ferromagnetic coupling between Cu neighbor cations. This magnetic order could be anticipated from the close correspondence between CuF2 and rutile crystal structures.