The oxygen migration in the apatite-type lanthanum silicate with the cation substitution

A theoretical model is proposed to determine the effects of Si substitution with Al on the oxygen diffusion in apatite-type lanthanum silicates based on density-functional theory (DFT) calculations for La10(SiO 4)4(AlO4)2O2. Substitution changes the stable configuration for excess oxygen from the split interstitial to a new cluster form with the original cluster. Al doping completely changes the migration mechanism from the interstitialcy one, which was proposed for the La9.33(SiO4)6O2 starting material, to a mechanism which contains an interstitial process. Nevertheless, the migration barrier is calculated to be 0.81 eV, which indicates small changes in oxygen conduction and is consistent with the observations. The present study indicates that the cation substitution on silicon site alone does not promise the improvement of the oxide ion conduction in the lanthanum silicate.

Position preference and diffusion path of an oxygen ion in apatite-type lanthanum silicate La9.33Si6O26: A density functional study

Using density functional theory, we investigated the position preference and diffusion mechanisms of interstitial oxygen ions in lanthanum silicate La9.33Si6O26, which is an apatite-structured oxide and a promising candidate electrolyte material for solid oxide fuel cells. The reported lanthanum vacancies were explicitly taken into account by theoretically determining their arrangement with a supercell model. The most stable structures and the formation energies of oxygen interstitials were determined for each charged state. It was found that the double-negatively charged state is stable over a wide range of the Fermi level, and that the excess oxygen ions form split interstitials with the original oxygen ions, while the neutral and the single-negatively charged states preferably form molecular oxygen. These species were found near the lanthanum vacancy site. The theoretically determined migration pathway along the c-axis essentially follows an interstitialcy mechanism. The obtained migration barrier is sensitive to the charge state, and is also affected by the lanthanum vacancy. The barrier height of the double-negatively charged state was calculated to be 0.58 eV for the model structure, which is consistent with the measured activation energy.

γ-Y2Si2O7, a machinable silicate ceramic: Mechanical properties and machinability

In this paper, the mechanical properties of bulk single-phase γ-Y2Si2O7 ceramic are reported. γ-Y2Si2O7 exhibits low shear modulus, excellent damage tolerance, and thus has a good machinability ready for metal working tools. To understand the underlying mechanism of machinability, drilling test, Hertzian contact test, and density functional theory (DFT) calculation are employed. Hertzian contact test demonstrates that γ-Y2Si2O7 is a "quasi-plastic" ceramic and the intrinsically weak interfaces contribute to its machinability. Crystal structure characteristics and DFT calculations of γ-Y2Si2O7 suggest that some weakly bonded planes, which involve Y-O bonds that can be easily broken, are the sources of the low shear deformation resistance and good machinability.

Predicting single-layer technetium dichalcogenides (TcX2, X = S, Se) with promising applications in photovoltaics and photocatalysis

One of the least known compounds among transition metal dichalcogenides (TMDCs) is the layered triclinic technetium dichalcogenides (TcX2, X = S, Se). In this work, we systematically study the structural, mechanical, electronic, and optical properties of TcS2 and TcSe2 monolayers based on density functional theory (DFT). We find that TcS2 and TcSe2 can be easily exfoliated in a monolayer form because their formation and cleavage energy are analogous to those of other experimentally realized TMDCs monolayer. By using a hybrid DFT functional, the TcS2 and TcSe2 monolayers are calculated to be indirect semiconductors with band gaps of 1.91 and 1.69 eV, respectively. However, bilayer TcS2 exhibits direct-bandgap character, and both TcS2 and TcSe2 monolayers can be tuned from semiconductor to metal under effective tensile/compressive strains. Calculations of visible light absorption indicate that 2D TcS2 and TcSe2 generally possess better capability of harvesting sunlight compared to single-layer MoS2 and ReSe2, implying their potential as excellent light-absorbers. Most interestingly, we have discovered that the TcSe2 monolayer is an excellent photocatalyst for splitting water into hydrogen due to the perfect fit of band edge positions with respect to the water reduction and oxidation potentials. Our predictions expand the two-dimensional (2D) family of TMDCs, and the remarkable electronic/optical properties of monolayer TcS2 and TcSe2 will place them among the most promising 2D TMDCs for renewable energy application in the future.

Irradiation-mediated tailoring of carbon nanotubes

The ever-increasing demand for faster computers in various areas, ranging from entertaining electronics to computational science, is pushing the semiconductor industry towards its limits on decreasing the sizes of electronic devices based on conventional materials. According to the famous law by Gordon E. Moore, a co-founder of the world s largest semiconductor company Intel, the transistor sizes should decrease to the atomic level during the next few decades to maintain the present rate of increase in the computational power. As leakage currents become a problem for traditional silicon-based devices already at sizes in the nanometer scale, an approach other than further miniaturization is needed to accomplish the needs of the future electronics. A relatively recently proposed possibility for further progress in electronics is to replace silicon with carbon, another element from the same group in the periodic table. Carbon is an especially interesting material for nanometer-sized devices because it forms naturally different nanostructures. Furthermore, some of these structures have unique properties. The most widely suggested allotrope of carbon to be used for electronics is a tubular molecule having an atomic structure resembling that of graphite. These carbon nanotubes are popular both among scientists and in industry because of a wide list of exciting properties. For example, carbon nanotubes are electronically unique and have uncommonly high strength versus mass ratio, which have resulted in a multitude of proposed applications in several fields. In fact, due to some remaining difficulties regarding large-scale production of nanotube-based electronic devices, fields other than electronics have been faster to develop profitable nanotube applications. In this thesis, the possibility of using low-energy ion irradiation to ease the route towards nanotube applications is studied through atomistic simulations on different levels of theory. Specifically, molecular dynamic simulations with analytical interaction models are used to follow the irradiation process of nanotubes to introduce different impurity atoms into these structures, in order to gain control on their electronic character. Ion irradiation is shown to be a very efficient method to replace carbon atoms with boron or nitrogen impurities in single-walled nanotubes. Furthermore, potassium irradiation of multi-walled and fullerene-filled nanotubes is demonstrated to result in small potassium clusters in the hollow parts of these structures. Molecular dynamic simulations are further used to give an example on using irradiation to improve contacts between a nanotube and a silicon substrate. Methods based on the density-functional theory are used to gain insight on the defect structures inevitably created during the irradiation. Finally, a new simulation code utilizing the kinetic Monte Carlo method is introduced to follow the time evolution of irradiation-induced defects on carbon nanotubes on macroscopic time scales. Overall, the molecular dynamic simulations presented in this thesis show that ion irradiation is a promisingmethod for tailoring the nanotube properties in a controlled manner. The calculations made with density-functional-theory based methods indicate that it is energetically favorable for even relatively large defects to transform to keep the atomic configuration as close to the pristine nanotube as possible. The kinetic Monte Carlo studies reveal that elevated temperatures during the processing enhance the self-healing of nanotubes significantly, ensuring low defect concentrations after the treatment with energetic ions. Thereby, nanotubes can retain their desired properties also after the irradiation. Throughout the thesis, atomistic simulations combining different levels of theory are demonstrated to be an important tool for determining the optimal conditions for irradiation experiments, because the atomic-scale processes at short time scales are extremely difficult to study by any other means.

Synergistic crystal facet engineering and structural control of WO3 films exhibiting unprecedented photoelectrochemical performance

WO3 nanoplate arrays with (002) oriented facets grown on fluorine doped SnO2 (FTO) glass substrates are tailored by tuning the precursor solution via a facile hydrothermal method. A 2-step hydrothermal method leads to the preferential growth of WO3 film with enriched (002) facets, which exhibits extraordinary photoelectrochemical (PEC) performance with a remarkable photocurrent density of 3.7 mA cm–2 at 1.23 V vs. revisable hydrogen electrode (RHE) under AM 1.5 G illumination without the use of any cocatalyst, corresponding to ~93% of the theoretical photocurrent of WO3. Density functional theory (DFT) calculations together with experimental studies reveal that the enhanced photocatalytic activity and better photo-stability of the WO3 films are attributed to the synergistic effect of highly reactive (002) facet and nanoplate structure which facilitates the photo–induced charge carrier separation and suppresses the formation of peroxo-species. Without the use of oxygen evolution cocatalysts, the excellent PEC performance, demonstrated in this work, by simply tuning crystal facets and nanostructure of pristine WO3 films may open up new opportunities in designing high performance photoanodes for PEC water splitting.

Experimental and Theoretical Charge Density Analysis of Polymorphic Structures: The Case of Coumarin 314 Dye

Experimental charge density distributions in two known conformational polymorphs (orange and yellow) of coumarin 314 dye are analyzed based on multipole modeling of X-ray diffraction data collected at 100 K. The experimental results are compared with the charge densities derived from multipole modeling of theoretical structure factors obtained from periodic quantum calculation with density functional theory (DFT) method and B3LYP/6-31G(d,p) level of theory. The presence of disorder at the carbonyl oxygen atom of ethoxycarbonyl group in the yellow form, which was not identified earlier, is addressed here. The investigationof intermolecular interactions, based on Hirshfeld surface analysis and topological properties via quantum theory of atoms in molecule and total electrostatic interaction energies, revealed significant differences between the polymorphs. The differences of electrostatic nature in these two polymorphic forms were unveiled via construction of three-dimensional deformation electrostatic potential maps plotted over the molecular surfaces. The lattice energies evaluated from ab initio calculations on the two polymorphic forms indicate that the yellow form is likely to be the most favorable thermodynamically. The dipole moments derived from experimental and theoretical charge densities and also from Lorentz tensor approach are compared with the single-molecule dipole moments. In each case, the differences of dipole moments between the polymorphs are identified.

Molecular theory of solvation and solvation dynamics of a classical ion in a dipolar liquid

A microscopic theory of equilibrium solvation and solvation dynamics of a classical, polar, solute molecule in dipolar solvent is presented. Density functional theory is used to explicitly calculate the polarization structure around a solvated ion. The calculated solvent polarization structure is different from the continuum model prediction in several respects. The value of the polarization at the surface of the ion is less than the continuum value. The solvent polarization also exhibits small oscillations in space near the ion. We show that, under certain approximations, our linear equilibrium theory reduces to the nonlocal electrostatic theory, with the dielectric function (c(k)) of the liquid now wave vector (k) dependent. It is further shown that the nonlocal electrostatic estimate of solvation energy, with a microscopic c(k), is close to the estimate of linearized equilibrium theories of polar liquids. The study of solvation dynamics is based on a generalized Smoluchowski equation with a mean-field force term to take into account the effects of intermolecular interactions. This study incorporates the local distortion of the solvent structure near the ion and also the effects of the translational modes of the solvent molecules.The latter contribution, if significant, can considerably accelerate the relaxation of solvent polarization and can even give rise to a long time decay that agrees with the continuum model prediction. The significance of these results is discussed.

First-principles analysis of electron correlation, spin ordering and phonons in the normal state of FeSe1-x

We present first-principles density-functional-theory-based calculations to determine the effects of the strength of on-site electron correlation, magnetic ordering, pressure and Se vacancies on phonon frequencies and electronic structure of FeSe1-x. The theoretical equilibrium structure (lattice parameters) of FeSe depends sensitively on the value of the Hubbard parameter U of on-site correlation and magnetic ordering. Our results suggest that there is a competition between different antiferromagnetic states due to comparable magnetic exchange couplings between first- and second-neighbor Fe sites. As a result, a short range order of stripe antiferromagnetic type is shown to be relevant to the normal state of FeSe at low temperature. We show that there is a strong spin-phonon coupling in FeSe (comparable to its superconducting transition temperature) as reflected in large changes in the frequencies of certain phonons with different magnetic ordering, which is used to explain the observed hardening of a Raman-active phonon at temperatures (similar to 100 K) where magnetic ordering sets in. The symmetry of the stripe antiferromagnetic phase permits an induced stress with orthorhombic symmetry, leading to orthorhombic strain as a secondary order parameter at the temperature of magnetic ordering. The presence of Se vacancies in FeSe gives rise to a large peak in the density of states near the Fermi energy, which could enhance the superconducting transition temperature within the BCS-like picture.

Magnetically induced currents and magnetic response properties of molecules

The magnetically induced currents in organic monoring and multiring molecules, in Möbius shaped molecules and in inorganic all-metal molecules have been investigated by means of the Gauge-including magnetically induced currents (GIMIC) method. With the GIMIC method, the ring-current strengths and the ring-current density distributions can be calculated. For open-shell molecules, also the spin current can be obtained. The ring-current pathways and ring-current strengths can be used to understand the magnetic resonance properties of the molecules, to indirectly identify the effect of non-bonded interactions on NMR chemical shifts, to design new molecules with tailored properties and to discuss molecular aromaticity. In the thesis, the magnetic criterion for aromaticity has been adopted. According to this, a molecule which has a net diatropic ring current might be aromatic. Similarly, a molecule which has a net paratropic current might be antiaromatic. If the net current is zero, the molecule is nonaromatic. The electronic structure of the investigated molecules has been resolved by quantum chemical methods. The magnetically induced currents have been calculated with the GIMIC method at the density-functional theory (DFT) level, as well as at the self-consistent field Hartree-Fock (SCF-HF), at the Møller-Plesset perturbation theory of the second order (MP2) and at the coupled-cluster singles and doubles (CCSD) levels of theory. For closed-shell molecules, accurate ring-current strengths can be obtained with a reasonable computational cost at the DFT level and with rather small basis sets. For open-shell molecules, it is shown that correlated methods such as MP2 and CCSD might be needed to obtain reliable charge and spin currents. The basis set convergence has to be checked for open-shell molecules by performing calculations with large enough basis sets. The results discussed in the thesis have been published in eight papers. In addition, some previously unpublished results on the ring currents in the endohedral fullerene Sc3C2@C80 and in coronene are presented. It is shown that dynamical effects should be taken into account when modelling magnetic resonance parameters of endohedral metallofullerenes such as Sc3C2@C80. The ring-current strengths in a series of nano-sized hydrocarbon rings are related to static polarizabilities and to H-1 nuclear magnetic resonance (NMR) shieldings. In a case study on the possible aromaticity of a Möbius-shaped [16]annulene we found that, according to the magnetic criterion, the molecule is nonaromatic. The applicability of the GIMIC method to assign the aromatic character of molecules was confirmed in a study on the ring currents in simple monocylic aromatic, homoaromatic, antiaromatic, and nonaromatic hydrocarbons. Case studies on nanorings, hexaphyrins and [n]cycloparaphenylenes show that explicit calculations are needed to unravel the ring-current delocalization pathways in complex multiring molecules. The open-shell implementation of GIMIC was applied in studies on the charge currents and the spin currents in single-ring and bi-ring molecules with open shells. The aromaticity predictions that are made based on the GIMIC results are compared to other aromaticity criteria such as H-1 NMR shieldings and shifts, electric polarizabilities, bond-length alternation, as well as to predictions provided by the traditional Hückel (4n+2) rule and its more recent extensions that account for Möbius twisted molecules and for molecules with open shells.

New computational approaches for inelastic x-ray scattering

Inelastic x-ray scattering spectroscopy is a versatile experimental technique for probing the electronic structure of materials. It provides a wealth of information on the sample's atomic-scale structure, but extracting this information from the experimental data can be challenging because there is no direct relation between the structure and the measured spectrum. Theoretical calculations can bridge this gap by explaining the structural origins of the spectral features. Reliable methods for modeling inelastic x-ray scattering require accurate electronic structure calculations. This work presents the development and implementation of new schemes for modeling the inelastic scattering of x-rays from non-periodic systems. The methods are based on density functional theory and are applicable for a wide variety of molecular materials. Applications are presented in this work for amorphous silicon monoxide and several gas phase systems. Valuable new information on their structure and properties could be extracted with the combination of experimental and computational methods.

Mechanistic Aspects of Nucleophilic Substitution at Half-Sandwich Metal Complexes

The hydrolysis reactions of organometallic ruthenium(II) piano-stool complexes of the type Ru-II(eta(6)-cymene)(L)Cl](0/+) (1-5, where L = kappa(1)- or kappa(2)-1,1-bis(diphenylphosphino)methane,1,1bis-(diphenylphosphino)methane oxide, kappa(1)-mercaptobenzothiazole) have been studied using density functional theory at the B3LYP level. In addition to considering a syn attack in an associative fashion, where the nucleophile approaches from the same side as the leaving group, we have explored alternative paths such as an anti attack in an associative manner, where the nucleophile attacks from the opposite side of the leaving group. During the anti attack, an intermediate is formed and there is a coordination mode change of the arene ring from eta(6) to eta(2) along with its rotation. When the intermediate goes to the product, the arene ring slips back from eta(2) to eta(6) coordination. This coordinated movement of the arene ring makes the associative anti attack an accessible pathway for the substitution process. Our calculations predict very similar activation barriers for both syn and anti attacks. In the dissociative path, the rate-determining step is the generation of a coordinatively unsaturated 16-electron ruthenium species. This turns out to be viable once solvent effects are included. The large size of the ancillary ligands on Ru makes the dissociative process as favorable as the associative process. Activation energy calculations reveal that although the dissociative path is favorable for kappa(1) complexes, both dissociative and associative processes can have significant contribution to the hydrolysis reaction in kappa(2) complexes. Once activated by hydrolysis, these complexes react with guanine and adenine bases of DNA. The thermodynamic stabilities of complexes formed with the nucleobases are also presented.

Raman evidence for the superconducting gap and spin-phonon coupling in the superconductor Ca( Fe0.95Co0.05)(2)As-2

Inelastic light scattering studies on a single crystal of electron-doped Ca(Fe0.95Co0.05)(2)As-2 superconductor, covering the tetragonal-to-orthorhombic structural transition as well as the magnetic transition at T-SM similar to 140 K and the superconducting transition temperature T-c similar to 23 K, reveal evidence for superconductivity-induced phonon renormalization. In particular, the phonon mode near 260 cm(-1) shows hardening below T-c, signaling its coupling with the superconducting gap. All three Raman active phonon modes show anomalous temperature dependence between room temperature and T-c, i.e. the phonon frequency decreases with lowering temperature. Further, the frequency of one of the modes shows a sudden change in temperature dependence at TSM. Using first-principles density functional theory based calculations, we show that the low temperature phase (T-c < T < T-SM) exhibits short-ranged stripe antiferromagnetic ordering, and estimate the spin-phonon couplings that are responsible for these phonon anomalies.

Theoretical Studies on the Structure and Bonding of Metallacyclocumulenes, -cyclopentynes, and -cycloallenes

In this paper, we compare the electronic structure of the hafnacycloallene complex Cp(2)HfC(4)Rr'(2)R `' (5Hf), which was previously described by Erker et al., with those of the titanium, zirconium, and hafnium complexes Cp(2)M(eta(4)-RHC(4)HR) (3M; i.e. metallacyclopent-2,3,4-trienes, metallacyclocumulenes) and Cp(2)M(eta(2)-R(2)C(4)R(2)) (4M; i.e. 1-metallacyclopent-3-ynes) using density functional theory (BP86/LANL2DZ) calculations. Moreover, the eta(3)-phenylallenyl zirconocene complex 7Zr, which was synthesized by Wojcicki et al., is included for the comparison. These calculations and extended Huckel calculations show that the bonding in complex 5Hf is remarkably similar to that of complexes 4M and 7Zr. An analysis of the structural parameters and bonding reveals that the unique interaction of the internal carbon atoms along with the terminal carbon atoms with the bent-metallocene moiety is the reason for the unusual stability of these metallacycles. The molecular orbital analysis further suggests that complex 5Hf can react with another metal fragment to give the bimetallic complexes 9 and 10. The electronic structures of complexes 3M, 4M, 5Hf, and 7Zr have been comparatively studied to get a general understanding of the bonding in these metallacycles.

Surface effects on stacking fault and twin formation in fcc nanofilms: A first-principles study

The free surface effects on stacking fault and twin formation in fcc metals (Al, Cu, and Ni) were examined by first-principles calculations based on density functional theory (DFT). It is found that the generalized planar fault (GPF) energies of Ni are much larger than bulk Ni with respect to Al and Cu. The discrepancy is attributed to the localized relaxation of Ni nanofilm to accommodate the large expansion of the inter-planar separation induced at the fault plane. The localized relaxation can be coupled to the electronic structure of Ni nanofilms. (C) 2011 Elsevier B.V. All rights reserved.