Estudo da adsorção de hidrogênio e sulfeto na superfície de paládio: aspectos experimentais (eletroquímica) e teóricos (ab initio e Teoria do Funcional da Densidade)

The adsorption of H and S2- species on Pd (100) has been studied with ab initio, density-functional calculations and electrochemical methods. A cluster of five Pd atoms with a frozen geometry described the surface. The computational calculations were performed through the GAUSSIAN94 program, and the basis functions adapted to a pseudo-potential obtained by using the Generator Coordinate Method adapted to the this program. Using the cyclic voltammetry technique through a Model 283 Potentiostat/Galvanostat E.G.&G-PAR obtained the electrochemical results. The calculated chemisorption geometry has a Pd-H distance of 1.55Å, and the potential energy surface was calculated using the Becke3P86//(GCM/DFT/SBK) methodology. The adsorption of S2- ions on Pd surface obtained both through comparison between the experimental and theoretical results, at MP2 level, suggest a S2- absorption into the metallic cluster. The produced Pd-(S2-) system was show to be very stable under the employed experimental conditions. The paper has shows the powerful aid of computational methods to interpret adsorption experimental data.

Aspectos relacionados à utilização da equação logística quadrática em processos eletroquímicos

The concepts of dissipation and feedback are contained in the behavior of many natural dynamical systems. They have been used to predict the evolution of populations leading to the formulation of the quadratic logistic equation (QLE). More recently, the QLE has been used to provide a better understanding of physicochemical systems with promising results. Many physical, chemical and biological dynamic phenomena can be understood on the basis of the QLE and this work describes the main aspects of this equation and some recent applications, with emphasis on electrochemical systems. Also, it is illustrated the concept of potential energy as a convenient way of describing the stability of the fixed points of the QLE.

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.

Coordenadas cartesianas moleculares a partir da geometria dos modos normais de vibração

A simple method to obtain molecular Cartesian coordinates as a function of vibrational normal modes is presented in this work. The method does not require the definition of special matrices, like the F and G of Wilson, neither of group theory. The Eckart's conditions together with the diagonalization of kinetic and potential energy are the only required expressions. This makes the present approach appropriate to be used as a preliminary study for more advanced concepts concerning vibrational analysis. Examples are given for diatomic and triatomic molecules.

Operador dipolo-dipolo na base de momento angular: um complemento ao estudo de ressonância magnética nuclear

The relationship between the magnetic dipole-dipole potential energy function and its quantum analogue is presented in this work. It is assumed the reader is familiar with the classical expression of the dipolar interaction and has basic knowledge of the quantum mechanics of angular momentum. Except for these two points only elementary steps are involved.

Estudio teórico de la isomerización del ácido maleico en ácido fumárico: un enfoque basado en el concepto de superficie de energía potencial

A comparative study based on potential energy surfaces (PES) of 2-butanedioic and hypothetic 2-butanedioic/HCl acids is useful for understanding the maleic acid isomerization. The PES enables locating conformers of minimum energy, intermediates of reactions and transition states. From contour diagrams, a set of possible reaction paths are depicted interconnecting the proposed structures. The study was carried out in absentia and in the presence of the catalyst (HCl), using an solvatation model provided by the Gaussian software package. Clearly, the effect of HCl is given by new reaction paths with lower energetic barriers in relation to the reaction without catalyzing.

Heteropentámeros (etanol)4-agua: estudio estructural y termodinámico

Stochastic exploration of the potential energy surface of (ethanol)4-water heteropentamers through simulated annealing calculations was used to find probable structures of these clusters. Subsequent geometry optimization with the B3LYP/6-31+G(d) approach of these initial structures led to 13 stable heteropentamers. The strength of the hydrogen bonds of the type O"H-O (primary) and their spatial arrangements seem to be responsible for the geometric preferences and the high stability of these heteropentamers. This result is a consequence of the presence of the cooperative effects among such interactions. There is no significant influence of the secondary hydrogen bonds (C"H-O) on the stability of the heteropentamers.

Estudo teórico da estrutura eletrônica e da dinâmica induzida por lasers da molécula de HCl

Potential energy and dipole moment curves for the HCl molecule were computed. Calculations were performed at different levels of theory (DFT, MRCI). Spectroscopic properties are reported and compared with experimental data, for validating the theoretical approaches. Interaction of infrared radiation with HCl is simulated using the wave packet formalism. The quantum control model for population dynamics of the vibrational levels, based on pi-pulse theory, is applied. The results demonstrate that wavepackets with specific composition can be built with short infrared laser pulses and provide the basis for studies of H + HCl collision dynamics with infrared laser excitation.

Exploración estocástica de las superficies de energía potencial de dímeros cis-trans y trans-trans del ácido fórmico

Potential energy surface (PES) of cis-trans and trans-trans formic acid dimers were sampled using a stochastic method, and the geometries, energies, and vibrational frequencies were computed at B3LYP/6-311++G(3df,2p) level of theory. The results show that molar free energy of dimerization deviated up to 108.4% when basis set superposition error (BSSE) and zero-point energy (ZPE) were not considered. For cis-trans dimers, C=O and O - H bond weakened, whereas C - O bonds strengthened due to dimerization. Also, trans-trans FA dimers did not show a trend regarding strengthening or weakening of the C=O, O - H and C - O bonds.

Molecular modeling of asymmetric induction in heterogeneously catalyzed hydrogenation of the C=O bond

Modifiering av metallytor med starkt adsorberade kirala organiska molekyler är eventuellt den mest relevanta teknik man vet i dag för att skapa kirala ytor. Den kan utnyttjas i katalytisk produktion av enantiomeriskt rena kirala föreningar som behövs t.ex. som läkemedel och aromkemikalier. Trots många fördelar av asymmetrisk heterogen katalys jämfört med andra sätt för att få kirala föreningar, har den ändå inte blivit ett allmänt verktyg för storskaliga tillämpningar. Detta beror t.ex. på brist på djupare kunskaper i katalytiska reaktionsmekanismer och ursprunget för asymmetrisk induktion. I denna studie användes molekylmodelleringstekniker för att studera asymmetriska, heterogena katalytiska system, speciellt hydrering av prokirala karbonylföreningar till motsvarande kirala alkoholer på cinchona-alkaloidmodifierade Pt-katalysatorer. 1-Fenyl-1,2-propandion (PPD) och några andra föreningar, som innehåller en prokiral C=O-grupp, användes som reaktanter. Konformationer av reaktanter och cinchona-alkaloider (som kallas modifierare) samt vätebundna 1:1-komplex mellan dem studerades i gas- och lösningsfas med metoder som baserar sig på vågfunktionsteori och täthetsfunktionalteori (DFT). För beräkningen av protonaffiniteter användes också högst noggranna kombinationsmetoder såsom G2(MP2). Den relativa populationen av modifierarnas konformationer varierade som funktion av modifieraren, dess protonering och lösningsmedlet. Flera reaktant–modifierareinteraktionsgeometrier beaktades. Slutsatserna på riktning av stereoselektivitet baserade sig på den relativa termodynamiska stabiliteten av de diastereomeriska reaktant–modifierare-komplexen samt energierna hos π- och π*-orbitalerna i den reaktiva karbonylgruppen. Adsorption och reaktioner på Pt(111)-ytan betraktades med DFT. Regioselektivitet i hydreringen av PPD och 2,3-hexandion kunde förklaras med molekyl–yta-interaktioner. Storleken och formen av klustret använt för att beskriva Pt-ytan inverkade inte bara på adsorptionsenergierna utan också på de relativa stabiliteterna av olika adsorptionsstrukturer av en molekyl. Populationerna av modifierarnas konformationer i gas- och lösningsfas korrelerade inte med populationerna på Pt-ytan eller med enantioselektiviteten i hydreringen av PPD på Pt–cinchona-katalysatorer. Vissa modifierares konformationer och reaktant–modifierare-interaktionsgeometrier var stabila bara på metallytan. Teoretiskt beräknade potentialenergiprofiler för hydrering av kirala α-hydroxiketoner på Pt implicerade preferens för parvis additionsmekanism för väte och selektiviteter i harmoni med experimenten. De uppnådda resultaten ökar uppfattningen om kirala heterogena katalytiska system och kunde därför utnyttjas i utvecklingen av nya, mera aktiva och selektiva kirala katalysatorer.

The role of van der Waals interactions in the case of H2 adsorption on ruthenium (0001) surface

Computational material science with the Density Functional Theory (DFT) has recently gained a method for describing, for the first time the non local bonding i.e., van der Waals (vdW) bonding. The newly proposed van der Waals-Density Functional (vdW-DF) is employed here to address the role of non local interactions in the case of H2 adsorption on Ru(0001) surface. The later vdW-DF2 implementation with the DFT code VASP (Vienna Ab-initio Simulation Package) is used in this study. The motivation for studying H2 adsorption on ruthenium surface arose from the interest to hydrogenation processes. Potential energy surface (PES) plots are created for adsorption sites top, bridge, fcc and hcp, employing the vdW-DF2 functional. The vdW-DF yields 0.1 eV - 0.2 eV higher barriers for the dissociation of the H2 molecule; the vdW-DF seems to bind the H2 molecule more tightly together. Furthermore, at the top site, which is found to be the most reactive, the vdW functional suggests no entrance barrier or in any case smaller than 0.05 eV, whereas the corresponding calculation without the vdW-DF does. Ruthenium and H2 are found to have the opposite behaviors with the vdW-DF; Ru lattice constants are overestimated while H2 bond length is shorter. Also evaluation of the CPU time demand of the vdW-DF2 is done from the PES data. From top to fcc sites the vdW-DF computational time demand is larger by 4.77 % to 20.09 %, while at the hcp site it is slightly smaller. Also the behavior of a few exchange correlation functionals is investigated along addressing the role of vdW-DF. Behavior of the different functionals is not consistent between the Ru lattice constants and H2 bond lengths. It is thus difficult to determine the quality of a particular exchange correlation functional by comparing equilibrium separations of the different elements. By comparing PESs it would be computationally highly consuming.

Production of Mg(OH)2 from Mg-silicate rock for CO2 mineral sequestration

Sequestration of carbon dioxide in mineral rocks, also known as CO2 Capture and Mineralization (CCM), is considered to have a huge potential in stabilizing anthropogenic CO2 emissions. One of the CCM routes is the ex situ indirect gas/sold carbonation of reactive materials, such as Mg(OH)2, produced from abundantly available Mg-silicate rocks. The gas/solid carbonation method is intensively researched at Åbo Akademi University (ÅAU ), Finland because it is energetically attractive and utilizes the exothermic chemistry of Mg(OH)2 carbonation. In this thesis, a method for producing Mg(OH)2 from Mg-silicate rocks for CCM was investigated, and the process efficiency, energy and environmental impact assessed. The Mg(OH)2 process studied here was first proposed in 2008 in a Master’s Thesis by the author. At that time the process was applied to only one Mg-silicate rock (Finnish serpentinite from the Hitura nickel mine site of Finn Nickel) and the optimum process conversions, energy and environmental performance were not known. Producing Mg(OH)2 from Mg-silicate rocks involves a two-staged process of Mg extraction and Mg(OH)2 precipitation. The first stage extracts Mg and other cations by reacting pulverized serpentinite or olivine rocks with ammonium sulfate (AS) salt at 400 - 550 oC (preferably < 450 oC). In the second stage, ammonia solution reacts with the cations (extracted from the first stage after they are leached in water) to form mainly FeOOH, high purity Mg(OH)2 and aqueous (dissolved) AS. The Mg(OH)2 process described here is closed loop in nature; gaseous ammonia and water vapour are produced from the extraction stage, recovered and used as reagent for the precipitation stage. The AS reagent is thereafter recovered after the precipitation stage. The Mg extraction stage, being the conversion-determining and the most energy-intensive step of the entire CCM process chain, received a prominent attention in this study. The extraction behavior and reactivity of different rocks types (serpentinite and olivine rocks) from different locations worldwide (Australia, Finland, Lithuania, Norway and Portugal) was tested. Also, parametric evaluation was carried out to determine the optimal reaction temperature, time and chemical reagent (AS). Effects of reactor types and configuration, mixing and scale-up possibilities were also studied. The Mg(OH)2 produced can be used to convert CO2 to thermodynamically stable and environmentally benign magnesium carbonate. Therefore, the process energy and life cycle environmental performance of the ÅAU CCM technique that first produces Mg(OH)2 and the carbonates in a pressurized fluidized bed (FB) were assessed. The life cycle energy and environmental assessment approach applied in this thesis is motivated by the fact that the CCM technology should in itself offer a solution to what is both an energy and environmental problem. Results obtained in this study show that different Mg-silicate rocks react differently; olivine rocks being far less reactive than serpentinite rocks. In summary, the reactivity of Mg-silicate rocks is a function of both the chemical and physical properties of rocks. Reaction temperature and time remain important parameters to consider in process design and operation. Heat transfer properties of the reactor determine the temperature at which maximum Mg extraction is obtained. Also, an increase in reaction temperature leads to an increase in the extent of extraction, reaching a maximum yield at different temperatures depending on the reaction time. Process energy requirement for producing Mg(OH)2 from a hypothetical case of an iron-free serpentine rock is 3.62 GJ/t-CO2. This value can increase by 16 - 68% depending on the type of iron compound (FeO, Fe2O3 or Fe3O4) in the mineral. This suggests that the benefit from the potential use of FeOOH as an iron ore feedstock in iron and steelmaking should be determined by considering the energy, cost and emissions associated with the FeOOH by-product. AS recovery through crystallization is the second most energy intensive unit operation after the extraction reaction. However, the choice of mechanical vapor recompression (MVR) over the “simple evaporation” crystallization method has a potential energy savings of 15.2 GJ/t-CO2 (84 % savings). Integrating the Mg(OH)2 production method and the gas/solid carbonation process could provide up to an 25% energy offset to the CCM process energy requirements. Life cycle inventory assessment (LCIA) results show that for every ton of CO2 mineralized, the ÅAU CCM process avoids 430 - 480 kg CO2. The Mg(OH)2 process studied in this thesis has many promising features. Even at the current high energy and environmental burden, producing Mg(OH)2 from Mg-silicates can play a significant role in advancing CCM processes. However, dedicated future research and development (R&D) have potential to significantly improve the Mg(OH)2 process performance.

Contact with friction using the augmented Lagrangian Method: a conditional constrained minimization problem

This work presents a formulation of the contact with friction between elastic bodies. This is a non linear problem due to unilateral constraints (inter-penetration of bodies) and friction. The solution of this problem can be found using optimization concepts, modelling the problem as a constrained minimization problem. The Finite Element Method is used to construct approximation spaces. The minimization problem has the total potential energy of the elastic bodies as the objective function, the non-inter-penetration conditions are represented by inequality constraints, and equality constraints are used to deal with the friction. Due to the presence of two friction conditions (stick and slip), specific equality constraints are present or not according to the current condition. Since the Coulomb friction condition depends on the normal and tangential contact stresses related to the constraints of the problem, it is devised a conditional dependent constrained minimization problem. An Augmented Lagrangian Method for constrained minimization is employed to solve this problem. This method, when applied to a contact problem, presents Lagrange Multipliers which have the physical meaning of contact forces. This fact allows to check the friction condition at each iteration. These concepts make possible to devise a computational scheme which lead to good numerical results.

A new approach to evaluate imperfection sensitivity in asymmetric bifurcation buckling analysis

A direct procedure for the evaluation of imperfection sensitivity in bifurcation problems is presented. The problems arise in the context of the general theory of elastic stability for discrete structural systems, in which the energy criterion of stability of structures and the total potential energy formulation are employed. In cases of bifurcation buckling the sensitivity of the critical load with respect to an imperfection parameter e is singular at the state given by epsilon =0, so that, a regular perturbation expansion of the solution is not possible. In this work we describe a direct procedure to obtain the relations between the critical loads, the generalized coordinates at the critical state, the eigenvector, and the amplitude of the imperfection, using singular perturbation analysis. The expansions are assumed in terms of arbitrary powers of the imperfection parameter, so that both exponents and coefficients of the expansion are unknown. The solution of the series exponents is obtained by searching the least degenerate solution. The formulation is here applied to asymmetric bifurcations, for which explicit expressions of the coefficients are obtained. The use of the method is illustrated by a simple example, which allows consideration of the main features of the formulation.

Hydrologisen tilanteen ja tuulivoimatuotannon muutosten vaikutus vesiarvoihin pohjoismaisilla sähkön spot-markkinoilla

Vesivoimalla on merkittävä rooli pohjoismaisessa sähköntuotantojärjestelmässä. Spot-markkinoille tarjottavan vesivoiman hinta riippuu vaihtoehtoisen tuotannon hinnasta ja odotetusta vesivoimantuottajien käytettävissä olevasta veteen sitoutuneen potentiaalienergian määrästä. Hydrologisella tilanteella tarkoitetaan tässä tämän potentiaalienergian poikkeamaa normaalitasostaan. Viime vuosina tuulivoimatuotanto pohjoismaisella sähkömarkkina-alueella on kasvanut voimakkaasti, ja on tullut aiheelliseksi tarkastella, millaisia vaikutuksia tällä on vesivoimantuottajien toimintaan. Työssä määritellään vesivoimalle vaihtoehtoisen sähköntuotannon tuotantokustannus, joka pitkällä aikavälillä toimii vertailutasona, jonka perusteella vesivoimantuottajat määrittävät markkinoilla tarjontahinnan tuotannolleen. Tarkastellaan, kuinka hydrologisen tilanteen ja vaihtoehtoisen tuotannon tuotantokustannusten muutokset vaikuttavat vesiarvoon, joka on hinta, jolla hintariippuvaista eli säätyvää vesivoimaa tarjotaan spot-markkinoille. Todetaan, että hydrologisen tilanteen vahvistuminen ja vaihtoehtoisen tuotantokustannuksen aleneminen alentavat vesiarvoja. Todetaan lisäksi, että tuulivoima vaikuttaa sähkön hinnanmuodostukseen markkinoilla samankaltaisesti kuin hintariippumaton vesivoimatuotanto. Esitetään aikasarjamalli vesivoimatuotannon hintariippuvuuden mallintamiseksi. Vertaillaan vesivoimatuottajien toimintaa kahdella vesivoimatuotantoa sisältävällä hinta-alueella, joista toisella tuulivoimatuotanto on kasvanut voimakkaammin kuin toisella. Havaitaan, että molemmilla hinta-alueilla hydrologisen tilanteen vahvistuminen on alentanut ja heikkeneminen nostanut vesiarvoja. Lisäksi havaitaan, että alueella, jonka tuulivoimatuotanto on kasvanut enemmän, vesiarvot ovat laskeneet suhteessa alueen, jolla tuulivoimatuotanto on kasvanut vähemmän, vesiarvoihin. Tuulivoiman voidaan todeta syrjäyttäneen markkinoilta tuotantokustannuksiltaan kalliimpaa tuotantoa.