Molecular dynamics simulation of hydrogen enhancing dislocation emission

The interactive pair potential between Al and H is obtained based on the ab initio calculation and the Chen-Mobius 3D lattice inversion formula. By utilizing the pair potentials calculated, the effects of hydrogen on the dislocation emission from crack tip have been studied. The simulated result shows that hydrogen can reduce the cohesive strength for Al single crystal, and then the critical stress intensity factor for partial dislocation emission decreases from 0.11 MPa root m (C-H = 0) to 0.075 MPa root m (C-H=0.72%) and 0.06 MPa root m (C-H = 1.44%). This indicates thar hydrogen can enhance the dislocation emission. The simulation also shows that atoms of hydrogen can gather and turn into small bubbles, resulting in enhancement of the equilibrium vacancy concentration.

Surface atomic properties of tetrahedral amorphous carbon

The surface energy and surface atomic structure of tetrahedral amorphous carbon has been calculated by an ab-initio method. The surface atoms are found to reconstruct into sp2 sites often bonded in graphitic rings. Placing the dangling bonds on adjacent surface atoms lower their energy by π-bonding and this is the source of the low surface energy. The even lower surface energy of hydrogenated amorphous carbon (a-C:H) is due to the hydrogenation of all broken surface bonds. © 2005 Elsevier B.V. All rights reserved.

Theoretical investigation of epitaxial deformation and the hcp-bcc transition of alkali metals

The deformation of alkali metals K, Rb, and Cs under epitaxial deformation is studied via the ab initio pseudopotential plane wave method using the local-density approximation. Under loading from the stable fee phase, metastable stares along directions [001], [111], and [201] are identified. One metastable state, presented at direction [201], has a very low symmetry in contrast to the planes [001] and [201]. Our results show that the softening direction and sequences of growth is significantly affected by the existence of the metastable states and magnitude of the energy barrier. The resulting softening sequences from soft to hard are [201], [110], [001], and [111] under biaxial compression and [001], [111], [201], and [110] under biaxial tension. An orthorhombic deformation path is used to investigate the fact, that the structure of the alkali films K and Cs evolve from the quasihexagonal structure into the (110)-oriented bcc structure, observed by experiments.

Estudio de la factibilidad de aplicación, a escala planta piloto, del catalizador Cr2O3/Al2O3 en la reducción de emisión de SO2 a la atmósfera, proveniente de fuentes fijas : avances 2013

Resumen: Se planificaron las experiencias con el objeto de analizar el comportamiento del catalizador en la columna metálica de mayor diámetro. Se modificaron las masas usadas para verificar la eficiencia de retención respecto de la masa. Se realizaron ciclos de adsorción, desorción y readsorción sobre una misma muestra para determinar variaciones en la eficiencia del catalizador. En otra fase, en colaboración con el Dr. V. A. Ranea y el Dr. E. E. Mola (INIFTA, UNLP), se desarrolló el estudio teórico de la adsorción de moléculas de SO2, CH4, CO2, O2 y CO sobre Cr2O3(0001) mediante Teoría del Funcional Densidad (programa VASP, Vienna Ab-initio Simulation Package), y el estudio de la cinética de la reacción entre CH4, SO2 y el O2 junto con la presencia de especies sulfito y sulfato. Este estudio permitió hallar los sitios preferenciales de adsorción de S0 y la posible competencia con SO2 experimentalmente y por cálculos teóricos. Experimentalmente, se observa que la eficiencia de adsorción del catalizador respecto al SO2 es cercana al 100%. Se observa un pico de termodesorción a 1120 K. Luego, se estudió la oxidación de CH4 con SO2. Se observa que hay producción de CO2 desde temperatura inicial, seguida de un aumento significativo en la formación de CO2 hasta 330-340 K. Luego, la producción de CO2 se mantiene aproximadamente constante. Mediante el empleo de la ecuación de Arrhenius y resultados experimentales, se obtuvo la energía de activación de la reacción global, de 7 Kcal/mol. Mediante estudios teóricos, se determinó que la energía de quimisorción del SO2 sobre el Cr2O3 es de -3.09 eV para la configuración más estable, una energía de adsorción de O2 en estado disociativo de -1.567 eV, una energía para CH4 sobre O2 adsorbido previamente de -0.335 eV, y -0.812 eV para la configuración más estable de CO2 sobre el sustrato.

Estudio de la factibilidad de aplicación, a escala planta piloto, del catalizador Cr2O3 soportado en alúmina en la reducción de emisión de SO2 a la atmósfera, proveniente de fuentes fijas : avances 2012

Resumen: Se realizaron las experiencias planificadas con el objeto de analizar el comportamiento del catalizador en la columna metálica, para simular condiciones de planta piloto. Se modificaron las masas de catalizador y se realizaron ciclos de adsorción, desorción y readsorción sucesivos sobre una misma muestra, con lo que se determinaron variaciones en la eficiencia del mismo. En otra fase se desarrolló el estudio teórico de la adsorción de O2 y CO2 sobre el mismo sustrato, y el estudio de la cinética de la reacción entre CH4, SO2 y el O2 por medio del programa VASP (Vienna Ab-initio Simulation Package). Se verificó, a través de los datos experimentales y teóricos y en colaboración con el Dr. V. A. Ranea y el Prof. E. E. Mola (INIFTA, UNLP), la presencia de especies sulfito y sulfato sobre la superficie del soporte. Experimentalmente, se observa que la eficiencia de adsorción del catalizador respecto al SO2 es cercana al 100%. Se observa un pico de termodesorción a 1120 K. Luego, se estudió la oxidación de CH4 con SO2. Se observa que hay producción de CO2 desde temperatura inicial, seguida de un aumento en la formación de CO2 hasta 330-340 K. Luego, la producción de CO2 se mantiene aproximadamente constante. Mediante el empleo de la ecuación de Arrhenius y resultados experimentales, se obtuvo la energía de activación de la reacción global, de 7 Kcal/mol. También se observó que el incremento del flujo de SO2 a valores superiores a 200 ml/min no incrementa la cantidad de SO2 retenida en el rango de 923-1023K. Para un incremento de masa de sustrato catalítico de 0,025 a 0,050 g, la masa retenida de SO2 se incrementa un 70,61%. Mediante estudios teóricos, se determinó que la energía de quimisorción del SO2 sobre el Cr2O3 es de -3.09 eV para la configuración más estable.

Estudio de la factibilidad de aplicación, a escala planta piloto, del catalizador Cr2O3, soportado en alúmina, en la reducción de emisión de SO2 a la atmósfera proveniente de fuentes fijas

Resumen: Se propone utilizar un óxido como el Cr2O3 como catalizador ya que se ha determinado anteriormente, en la primera etapa de esta investigación, (“Estudio comparativo de la retención de SO2 sobre óxidos de metales de transición soportados en alúmina”), que la retención de SO2 sobre su superficie es un proceso de quimisorción con formación de especies sulfito superficiales sobre sitios básicos y un proceso de óxido reducción del ión metálico. Apoya este mecanismo el hecho de que la cantidad de SO2 adsorbido es función de la temperatura. La mayor eficiencia del Cr2O3 puede explicarse en base a sus propiedades superficiales, lo cual ha sido utilizado en la segunda etapa de reacción de reducción, ya que se ha completado la etapa inicial de quimisorción. En la segunda etapa de esta investigación (“Estudio de la reacción de reducción de SO2 con CH4 a altas temperaturas sobre catalizador de Cr2O3 soportado en alúmina”), se apuntó al estudio de un nuevo tipo de sinergia entre propiedades ácido-base y propiedades redox en una misma superficie. La tercera etapa apuntó a determinar la influencia que tiene el O2 en este proceso, ya que el O2 se encuentra presente en las chimeneas industriales en las condiciones de reacción entre el SO2 y el CH4, y produce modificaciones en los parámetros de reacción. Se experimentó con diferentes masas de catalizador y flujos de los distintos gases, y se estudió la influencia de la presencia de oxígeno en la reacción y particularmente con diferentes flujos del mismo, y la posibilidad de regeneración del catalizador.En esta cuarta y última etapa se están estudiando los cambios que se producen en la reacción al pasar de escala laboratorio a planta piloto utilizando una columna de mayor diámetro construída en metal. A través de los datos experimentales se está estudiando, en conjunto con el INIFTA, la presencia de especies sulfito y sulfato sobre la superficie del soporte. Adicionalmente, por medio del programa VASP (Vienna Ab-initio Simulation Package), se analiza la interacción entre los reactivos gaseosos y el soporte.

Quantum Monte Carlo study of the CO interaction with a dimer model surface for Cr(110)

The chemisorption of CO on a Cr( 110) surface is investigated using the quantum Monte Carlo method in the diffusion Monte Carlo (DMC) variant and a model Cr2CO cluster. The present results are consistent with the earlier ab initio HF study with this model that showed the tilted/ near-parallel orientation as energetically favoured over the perpendicular arrangement. The DMC energy difference between the two orientations is larger (1.9 eV) than that computed in the previous study. The distribution and reorganization of electrons during CO adsorption on the model surface are analysed using the topological electron localization function method that yields electron populations, charge transfer and clear insight on the chemical bonding that occurs with CO adsorption and dissociation on the model surface.

First-principles study of the electronic and magnetic properties of defective carbon nanostructures

149 p.: graf.

有分子内态激发和化学反应的流动的数值模拟

本文旨在研究双原子分子化学反应的微观模型、表面复合反应机制和分子 相继碰撞中的能量传递过程等问题。这些问题的本质是如何考虑分子和原子之 间的相互作用，传统的研究途径是从能量观点出发的。本文尝试将力学的观点 和方法，与分子模型和统计平均手段结合起来，在计算机的帮助下，分析和研 究这些问题。 本文所建立的依赖空间取向化学反应模型认为：化学键断裂的充要条件是 作用于化学键的拉伸力超过它的抗拉极限。由此导出的双原子分子的反应微观 判据及反应速率常数Arrhenius形式的表达式，反映了碰撞空间取向以及反应势 能面的影响。对这个模型的检验和应用表明：它能够解释了置换反应速率常数 实验数据拟合公式k_f=aTbexp(-Ts/T)中Ts<ab initio)方法。它的突出特点是能够模拟分子之间的相继碰撞，直接给出了分子转动能或振动能随碰撞次数的变化过程，这对于了解和认识热力学非平衡现象特别有帮助口对同核双原子分子和惰性气体分子碰撞 过程的模拟结果，与实验或理论结果符合，说明单个靶分子方法的基本想法和 实施步骤，以及我们对惰性气体分子和同核双原子分子及其原子相互作用的考 虑，是现实可行的。 此外，本文附录还简要介绍了作者在读博士期间完成的其它工作，包括一 种自然的宏观流动量分解方法及封闭的平均流动方程、过渡领域三维复杂外形 绕流场的直接统计模拟以及过渡领域高超声速圆柱绕流的直接统计模拟等。

Hybrid Qm/Mm Simulation Of The Hydration Phenomena Of Dipalmitoylphosphatidylcholine Headgroup

The polar headgroup of dipalmitoylphosphatidylcholine (DPPC) molecule both in gas phase and aqueous Solution is investigated by the hybrid quantum mechanical/molecular mechanical (QM/MM) method, in which the polar head of DPPC molecule and the bound water molecules are treated with density functional theory (DFT), while the apolar hydrocarbon chain of DPPC molecule is treated with MM method. It is demonstrated that the hybrid QM/MM method is both accurate and efficient to describe the conformations of DPPC headgroup. Folded structures of headgroup are found in gas phase calculations. In this work, both monohydration and polyhydration phenomena are investigated. In monohydration, different water association sites are studied. Both the hydration energy and the quantum properties of DPPC and water molecules are calculated at the DFT level of theory after geometry optimization. The binding force of monohydration is estimated by using the scan method. In polyhydration, more extended conformations are found and hydration energies in different polyhydration styles are estimated. (C) 2008 Elsevier Inc. All rights reserved.

Vibrational predissociation spectroscopy of mass-selected ionic clusters

This thesis presents structural investigations of molecular ions and ionic clusters using vibrational predissociation spectroscopy. Experimentally, a pulsed beam of the mass-selected ion is crossed by a tunable infrared laser beam generated by a Nd:YAG pumped LiNbO_3 optical parametric oscillator. The resulting fragment ion is mass-analyzed and detected, with its intensity as a function of the laser wavelength being the "action" spectrum of the parent ion. In the case of SiH_7^+, we observed a vibrational band centered at 3866 cm^(-1) with clear P, Q, R branches, which is assigned as a perturbed H_2 stretch. The absence of a second H_2 band suggests that the ion forms a symmetric complex with a structure H_2•SiH_3^+•H_2 , in contrast to the species CH_7^+, which has the structure CH_5^+•H_2. The infrared spectra of NO_2^+(H_2O)_n clusters exhibit a marked change with cluster size, indicating that an intracluster reaction occurs with sufficient solvation. Specifically, in NO_2^+(H_2O)_n clusters where n≤3, H_2O binds to a nitronium ion core; but at n=4 the NO_2^+ reacts, transforming the cluster to a new structure of H_3O^+•(H_2O)_(n_2)•HNO_3. For protonated chlorine nitrate, we have observed two distinct isomers previously predicted by ab initio calculations: NO_2^+•(HOC1), the lowest energy isomer, and (ClO)(HO)NO^+, a covalently bonded isomer about 20 kcal/mol higher in energy. Both isomers decompose to NO_2^+ and HOCl upon photo-excitation. These results for HClONO_2^+ lend strong support to the involvement of an ionic mechanism in the reaction of ClONO_2 on polar stratospheric cloud surfaces, a critical step in the dramatic springtime depletion of ozone over Antarctica. Current research activities on halide-solvent clusters and metal-ligand complexes as well as technological improvements of the apparatus are also discussed.

Electronic structure and phonon thermodynamics of iron alloys

Inelastic neutron scattering (INS) and nuclear-resonant inelastic x-ray scattering (NRIXS) were used to measure phonon spectra of FeV as a B2- ordered compound and as a bcc solid solution. Contrary to the behavior of ordering alloys studied to date, the phonons in the B2-ordered phase are softer than in the solid solution. Ordering increases the vibrational entropy, which stabilizes the ordered phase to higher temperatures. Ab initio calculations show that the number of electronic states at the Fermi level increases upon ordering, enhancing the screening between ions, and reducing the interatomic force constants. The effect of screening is larger at the V atomic sites than at the Fe atomic sites.

The phonon spectra of Au-rich alloys of fcc Au-Fe were also measured. The main effect on the vibrational entropy of alloying comes from a stiffening of the Au partial phonon density of states (DOS) with Fe concentration that increases the miscibility gap temperature. The magnitude of the effect is non- linear and it is reduced at higher Fe concentrations. Force constants were calculated for several compositions and show a local stiffening of Au–Au bonds close to Fe atoms, but Au–Au bonds that are farther away do not show this effect. Phonon DOS curves calculated from the force constants reproduced the experimental trends. The Au–Fe bond is soft and favors ordering, but a charge transfer from the Fe to the Au atoms stiffens the Au–Au bonds enough to favor unmixing. The stiffening is attributed to two main effects comparable in magnitude: an increase in electron density in the free-electron-like states, and stronger sd-hybridization.

INS and NRIXS measurements were performed at elevated temperatures on B2-ordered FeTi and NRIXS measurements were performed at high pressures. The high-pressure behavior is quasi- harmonic. The softening of the phonon DOS curves with temperature is strongly nonharmonic. Calculations of the force constants and Born-von Karman fits to the experimental data show that the bonds between second nearest neighbors (2nn) are much stiffer than those between 1nn, but fits to the high temperature data show that the former softens at a faster rate with temperature. The Fe–Fe bond softens more than the Ti–Ti bond. The unusual stiffness of the 2nn bond is explained by the calculated charge distribution, which is highly aspherical and localized preferentially in the t2g orbitals. Ab initio molecular dynamics (AIMD) simulations show a charge transfer from the t2g orbitals to the eg orbitals at elevated temperatures. The asphericity decreases linearly with temperature and is more severe at the Fe sites.

Theoretical studies of novel organic systems : biradicals, polyradicals, and conducting polymers

Chapter 1

Cyclobutanediyl has been studied in both its singlet and triplet states by ab initio electronic structure theory. The triplet, which is the ground state of the molecule, exists in both C_(2h) and C_(2v) forms, which interconvert via a C_s transition state. For the singlet, only a C_(2h) form is found. It passes, via a C_s transition state, onto the C_(2v) surface on which bicyclobutane is the only minimum. The ring-flipping (inversion) process in bicyclobutane includes the singlet biradical as an intermediate, and involves a novel, nonleast motion pathway. Semiclassical periodic orbit theory indicates that the various minima on both the singlet and triplet surfaces can interconvert via quantum mechanical tunneling.

Chapter 2

The dimethylenepolycyclobutadienes (n) are the non-Kekulé analogues of the classical acenes. Application of a variety of theoretical methods reveals several novel features of such structures. Most interesting is the emergence of a parity rule. When n is even, n is predicted to be a singlet, with n disjoint NBMOs. When n is odd, theory predicts a triplet ground state with (n+1) NBMOs that are not fully disjoint.

Chapter 3

Bi(cyclobutadienyl) (2), the cyclobutadiene analogue of biphenyl, and its homologues tri- (3) and tetra(cyclobutadienyl) (4) have been studied using electronic structure theory. Ab initio calculations on 2 reveal that the central bond is a true double bond, and that the structure is best thought of as two allyl radicals plus an ethylene. The singlet and triplet states are essentially degenerate. Trimer 3 is two allyls plus a dimethylenecyclobutanediyl, while 4 is two coplanar bi(cyclobutadienyl) units connected by a single bond. For both 3 and 4, the quintet, triplet, and singlet states are essentially degenerate, indicating that they are tetraradicals. The infinite polymer, polycyclobutadiene, has been studied by HMO, EHCO, and VEH methods. Several geometries based on the structures of 3 and 4 have been studied, and the band structures are quite intriguing. A novel crossing between the valence and conduction bands produces a small band gap and a high density of states at the Fermi level.

Chapter 4

At the level of Hückel theory, polyfulvene has a HOCO-LUCO degeneracy much like that seen in polyacetylene. Higher levels of theory remove the degeneracy, but the band gap (E_g) is predicted to be significantly smaller than analogous structures such as polythiophene and polypyrrole at the fulvenoid geometry. An alternative geometry, which we have termed quinoid, is also conceivable for polyfulvene, and it is predicted to have a much larger E_g. The effects of benzannelation to produce analogues of polyisothianaphthene have been evaluated. We propose a new model for such structures based on conventional orbital mixing arguments. Several of the proposed structures have quite interesting properties, which suggest that they are excellent candidates for conducting polymers.

Chapter 5

Theoretical studies of polydimethylenecyclobutene and polydiisopropylidene- cyclobutene reveal that, because of steric crowding, they cannot achieve a planar, fully conjugated structure in either their undoped or doped states. Rather, the structure consists of essentially orthogonal hexatriene units. Such a structure is incompatible with conventional conduction mechanisms involving polarons and bipolarons.

Spectroscopic studies of gas-phase molecular clusters

Spectroscopic investigations of hydrogen-bonding and van der Waals' interactions m molecular clusters were studied by the techniques of infrared predissociation and resonance-enhanced multiphoton ionization spectroscopies (REMPI). Ab initio calculations were applied in conjunction for data interpretation.

The infrared predissociation spectroscopy of CN^-•(H_2O)_n (n = 2 - 6) clusters was reported in the region of 2950 - 3850 cm^(-1). The hydrogen bondings for the C-site and N-site binding, and among the water molecules were identified for n = 2 to 4. A spectral transition was observed for n = 5 and 6, implying that the anion was surface-bound onto the water aggregates in larger clusters.

The infrared predissociation spectroscopy of Br^-•(NH_3) and I^-•(NH_3)_n (n =1-3) clusters was reported in the region of 3050-3450 cm^(-1). For the Br^-•(NH_3) complex, a dominating ionic NH stretch appeared at 3175 cm^(-1), and the weaker free NH stretch appeared at 3348 cm^(-1). The observed spectrum was consistent to the structure in which there was one nearly linear hydrogen bond between Br^- and the NH_3 moiety. For the I^- •(NH_3) complex, five distinct IR absorption bands were observed in the spectrum. The spectrum was not consistent with basic frequency patterns of three geometries considered in the ab initio calculations - complex with one, two and three hydrogen bondings between I^- and the NH_3 moiety. Substantial inhomogenous broadening were displayed in the spectra for I^-•(NH_3)_n (n =2-3), suggesting the presence of multiple isomers.

The REMPI spectroscopy of the bound 4p ^2П 1/2 and ^2П 3/2 states, and the dissociative 3d ^2Σ^+ 1/2 state in the Al•Ar complex was reported. The dissociative spectrum at Al^+ channel suggested the coupling of the 4p ^2П 1/2,3/2 states to the repulsive 3d ^2Σ^+ 1/2 state. The spin-electronic coupling was further manifested in the dissociative Al^+ spectrum of the 3d ^2Σ^+ 1/2 state. Using the potential energy curves obtained from ab initio calculations, a bound → continuum Franck-Condon-intensity simulation was performed and compared with the one-photon 3d ^2Σ^+ 1/2 profile. The agreement provided evidence for the petturbation above the Al(3d)Ar dissociation limit, and the repulsive character of the 3d ^2Σ^+ 1/2 state.