Nature of the electronic band gap in lanthanide oxides

Accurate electronic structures of the technologically important lanthanide/rare-earth sesquioxides (Ln2O3, with Ln=La, ⋯,Lu) and CeO2 have been calculated using hybrid density functionals HSE03, HSE06, and screened exchange (sX-LDA). We find that these density functional methods describe the strongly correlated Ln f electrons as well as the recent G0W0@LDA+U results, generally yielding the correct band gaps and trends across the Ln period. For HSE, the band gap between O 2p states and lanthanide 5d states is nearly independent of the lanthanide, while the minimum gap varies as filled or empty Ln 4f states come into this gap. sX-LDA predicts the unoccupied 4f levels at higher energies, which leads to a better agreement with experiments for Sm2O 3, Eu2O3, and Yb2O3. © 2013 American Physical Society.

Regioselective Bis-Selenation of Allenes Catalyzed by Palladium Complexes: A Theoretical Study

The reaction mechanism of Pd(O)-catalyzed allene bis-selenation reactions is investigated by using density functional methods. The overall reaction mechanism has been examined. It is found that with the bulkier PMe3 ligand, the rate-determining step is the reductive elimination process, while allene insertion and reductive elimination processes are competitive for the rate-determining step with the PH3 ligand, indicating the importance of the ligand effect. For both cis and trans palladium complexes, allene insertion into the Pd-Se bond of the trans palladium complex using the internal carbon atom attached to the selenyl group is prefer-red among the four pathways of allene insertion processes. The formation of sigma-allyl and pi-allyl palladium complexes is favored over that of the sigma-vinyl palladium species. By using methylallene, the regioselectivity of monosubstituted allene insertion into the Pd-Se bond is analyzed.

Structural stability and electronic state of transition metal trimers

Ground state geometries were searched for transition metal trimers Sc-3, Y-3, La-3, Lu-3, Ti-3, Zr-3, and Hf-3 by density functional methods. For all the studied trimers, our calculation indicates that the ground state geometries are either equilateral triangle (Zr-3 and Hf-3) or near equilateral triangle (Ti-3, Sc-3, Y-3, La-3, and Lu-3). For rare earth trimers Sc-3, Y-3, La-3, and Lu-3, isosceles triangle (near equilateral triangle) at quartet state is the ground state. Isosceles triangle at doublet state is the competitive candidate for the ground state. For Zr-3 and Hf-3, equilateral triangle at singlet state is the most stable. For Ti-3, isosceles triangle (near equilateral triangle) at quintet state gives the ground state. For Sc-3, Zr-3, and Hf-3, where experimental results are available, the predicted geometries are in agreement with experiment in which the ground state is equilateral triangle (Zr-3) or fluxional (Sc-3 and Hf-3). For Y-3, the calculated geometry is in agreement with experimental observation and previous theoretical study that Y-3 is a bent molecule for the ground state.

Electronic structures of 3d-metal mononitrides

Bond distances, vibrational frequencies, electron affinities, ionization potentials, and dissociation energies of the title molecules in neutral, positively, and negatively charged ions were studied by use of density functional methods B3LYP, BLYP, BHLYP, BPW91, and B3PW91. The calculated results are compared with experiments and previous theoretical studies. It was found that the calculated properties are highly dependent on the functionals employed, in particular for the dissociation energy and vibrational frequency. For neutral species, pure density functional methods BLYP and BPW91 have relatively good performance in reproducing the experimental bond distance and vibrational frequency. For cations, hybrid exchange functional methods B3LYP and B3PW91 are good in predicting the dissociation energy. For both neutral and charged species, BHLYP tends to give smaller dissociation energy.

New insights on surface-enhanced Raman scattering based on controlled aggregation and spectroscopic studies, DFT calculations and symmetry analysis for 3,6-bi-2-pyridyl-1,2,4,5-tetrazine adsorbed onto citrate-stabilized gold nanoparticles

Asystematic study on the surface-enhanced Raman scattering (SERS) for 3,6-bi-2-pyridyl-1,2,4,5-tetrazine (bptz) adsorbed onto citrate-modified gold nanoparticles (cit-AuNps) was carried out based on electronic and vibrational spectroscopy and density functional methods. The citrate/bptz exchange was carefully controlled by the stepwise addition of bptz to the cit-AuNps, inducing flocculation and leading to the rise of a characteristic plasmon coupling band in the visible region. Such stepwise procedure led to a uniform decrease of the citrate SERS signals and to the rise of characteristic peaks of bptz, consistent with surface binding via the N heterocyclic atoms. In contrast, single addition of a large amount of bptz promoted complete aggregation of the nanoparticles, leading to a strong enhancement of the SERS signals. In this case, from the distinct Raman profiles involved, the formation of a new SERS environment became apparent, conjugating the influence of the local hot spots and charge-transfer (CT) effects. The most strongly enhanced vibrations belong to a(1) and b(2) representations, and were interpreted in terms of the electromagnetic and the CT mechanisms: the latter involving significant contribution of vibronic coupling in the system. Copyright (C) 2010 John Wiley & Sons, Ltd.

Effects of nanoscale spatial inhomogeneity in strongly correlated systems

We calculate ground-state energies and density distributions of Hubbard superlattices characterized by periodic modulations of the on-site interaction and the on-site potential. Both density-matrix renormalization group and density-functional methods are employed and compared. We find that small variations in the on-site potential v(i) can simulate, cancel, or even overcompensate effects due to much larger variations in the on-site interaction U-i. Our findings highlight the importance of nanoscale spatial inhomogeneity in strongly correlated systems, and call for a reexamination of model calculations assuming spatial homogeneity.

Electronic Properties of Water in Liquid Environment. A Sequential QM/MM Study Using the Free Energy Gradient Method

There is a continuous search for theoretical methods that are able to describe the effects of the liquid environment on molecular systems. Different methods emphasize different aspects, and the treatment of both the local and bulk properties is still a great challenge. In this work, the electronic properties of a water molecule in liquid environment is studied by performing a relaxation of the geometry and electronic distribution using the free energy gradient method. This is made using a series of steps in each of which we run a purely molecular mechanical (MM) Monte Carlo Metropolis simulation of liquid water and subsequently perform a quantum mechanical/molecular mechanical (QM/MM) calculation of the ensemble averages of the charge distribution, atomic forces, and second derivatives. The MP2/aug-cc-pV5Z level is used to describe the electronic properties of the QM water. B3LYP with specially designed basis functions are used for the magnetic properties. Very good agreement is found for the local properties of water, such as geometry, vibrational frequencies, dipole moment, dipole polarizability, chemical shift, and spin-spin coupling constants. The very good performance of the free energy method combined with a QM/MM approach along with the possible limitations are briefly discussed.

Structural Analysis of alpha-Fetoprotein (AFP)-like Peptides with Anti-Breast-Cancer Properties

The abundance of alpha-fetoprotein (AFP), a natural protein produced by the fetal yolk sac during pregnancy, correlates with lower incidence of estrogen receptor positive (ER+) breast cancer. The pharmacophore region of AFP has been narrowed down to a four amino acid (AA) region in the third domain of the 591 AA peptide. Our computational study focuses on a 4-mer segment consisting of the amino acids threonine-proline-valine-asparagine (TPVN). We have run replica exchange molecular dynamics (REMD) simulations and used 120 configurational snapshots from the total trajectory as starting configurations for quantum chemical calculations. We optimized structures using semiempirical (PM3, PM6, PM6-D2, PM6-H2, PM6-DH+, PM6-DH2) and density functional methods (TPSS, PBE0, M06-2X). By comparing the accuracy of these methods against RI-MP2 benchmarks, we devised a protocol for calculating the lowest energy conformers of these peptides accurately and efficiently. This protocol screens out high-energy conformers using lower levels of theory and outlines a general method for predicting small peptide structures.

Adatom CCV Auger rates via the local density of states

In this letter the core-core-valence Auger transitions of an atomic impurity, both in bulk or adsorbed on a jellium-like surface, are computed within a DFT framework. The Auger rates calculated by the Fermi golden rule are compared with those determined by an approximate and simpler expression. This is based on the local density of states (LDOS) with a core hole present, in a region around the impurity nucleus. Different atoms, Na and Mg, solids, Al and Ag, and several impurity locations are considered. We obtain an excellent agreement between KL1V and KL23V rates worked out with the two approaches. The radius of the sphere in which we calculate the LDOS is the relevant parameter of the simpler approach. Its value only depends on the atomic species regardless of the location of the impurity and the type of substrate. (C) 2003 Elsevier B.V. All rights reserved.

Investigation of inter-ion interactions in N,N,N',N'-tetramethylethylenediammonium dithiocyanate via experimental and theoretical charge density studies

The crystal structure of the N,N,N',N'-tetramethylethylenediammonium dithiocyanate salt has been examined by experimental charge density studies from high-resolution X-ray diffraction data. The corresponding results are compared with multipole refinements, using theoretical structure factors obtained from a periodic density functional theory calculation at the B3LYP level with a 6-31G** basis set. The salt crystallizes in space group P (1) over bar and contains only a single ion pair with an inversion center in the cation. The salt has thus one unique classical N+-H center dot center dot center dot(NCS)(-) hydrogen bond but also has six other weaker interactions: four C-H center dot center dot center dot S, one C-H center dot center dot center dot N, and one C-H center dot center dot center dot C-pi. The nature of all these interactions has been examined topologically using Bader's quantum theory of "atoms in molecules" and all eight of the Koch-Popelier criteria. The experimental and theoretical approaches agree well and both show that the inter-ion interactions, even in this simplest of systems, play an integrated and complex role in the packing of the ions in the crystal. Electrostatic potential maps are derived from experimental charge densities. This is the first time such a system has been examined in detail by these methods.

Electronic structure of Nd1-xYxMnO3 from Mn K edge absorption spectroscopy and DFT methods

The electronic structure of Nd1-xYxMnO3 (x-0-0.5) is studied using x-ray absorption near-edge structure (XANES) spectroscopy at the Mn K-edge along with the DFT-based LSDA+U and real space cluster calculations. The main edge of the spectra does not show any variation with doping. The pre-edge shows two distinct features which appear well-separated with doping. The intensity of the pre-edge decreases with doping. The theoretical XANES were calculated using real space multiple scattering methods which reproduces the entire experimental spectra at the main edge as well as the pre-edge. Density functional theory calculations are used to obtain the Mn 4p, Mn 3d and O 2p density of states. For x=0, the site-projected density of states at 1.7 eV above Fermi energy shows a singular peak of unoccupied e(g) (spin-up) states which is hybridized Mn 4p and O 2p states. For x=0.5, this feature develops at a higher energy and is highly delocalized and overlaps with the 3d spin-down states which changes the pre-edge intensity. The Mn 4p DOS for both compositions, show considerable difference between the individual p(x), p(y) and p(z)), states. For x=0.5, there is a considerable change in the 4p orbital polarization suggesting changes in the Jahn-Teller effect with doping. (C) 2013 Elsevier Ltd. All rights reserved.

Intermolecular interactions, charge-density distribution and the electrostatic properties of pyrazinamide anti-TB drug molecule: an experimental and theoretical charge-density study

An experimental charge-density analysis of pyrazinamide (a first line antitubercular drug) was performed using high-resolution X-ray diffraction data (sin theta/lambda)(max) = 1.1 angstrom(-1)] measured at 100 (2) K. The structure was solved by direct methods using SHELXS97 and refined by SHELXL97. The total electron density of the pyrazinamide molecule was modeled using the Hansen-Coppens multipole formalism implemented in the XD software. The topological properties of electron density determined from the experiment were compared with the theoretical results obtained from CRYSTAL09 at the B3LYP/6-31G** level of theory. The crystal structure was stabilized by N-H center dot center dot center dot N and N-H center dot center dot center dot O hydrogen bonds, in which the N3-H3B center dot center dot center dot N1 and N3-H3A center dot center dot center dot O1 interactions form two types of dimers in the crystal. Hirshfeld surface analysis was carried out to analyze the intermolecular interactions. The fingerprint plot reveals that the N center dot center dot center dot H and O center dot center dot center dot H hydrogen-bonding interactions contribute 26.1 and 18.4%, respectively, of the total Hirshfeld surface. The lattice energy of the molecule was calculated using density functional theory (B3LYP) methods with the 6-31G** basis set. The molecular electrostatic potential of the pyrazinamide molecule exhibits extended electronegative regions around O1, N1 and N2. The existence of a negative electrostatic potential (ESP) region just above the upper and lower surfaces of the pyrazine ring confirm the pi-electron cloud.

Methods for melting temperature calculation

Melting temperature calculation has important applications in the theoretical study of phase diagrams and computational materials screenings. In this thesis, we present two new methods, i.e., the improved Widom's particle insertion method and the small-cell coexistence method, which we developed in order to capture melting temperatures both accurately and quickly.

We propose a scheme that drastically improves the efficiency of Widom's particle insertion method by efficiently sampling cavities while calculating the integrals providing the chemical potentials of a physical system. This idea enables us to calculate chemical potentials of liquids directly from first-principles without the help of any reference system, which is necessary in the commonly used thermodynamic integration method. As an example, we apply our scheme, combined with the density functional formalism, to the calculation of the chemical potential of liquid copper. The calculated chemical potential is further used to locate the melting temperature. The calculated results closely agree with experiments.

We propose the small-cell coexistence method based on the statistical analysis of small-size coexistence MD simulations. It eliminates the risk of a metastable superheated solid in the fast-heating method, while also significantly reducing the computer cost relative to the traditional large-scale coexistence method. Using empirical potentials, we validate the method and systematically study the finite-size effect on the calculated melting points. The method converges to the exact result in the limit of a large system size. An accuracy within 100 K in melting temperature is usually achieved when the simulation contains more than 100 atoms. DFT examples of Tantalum, high-pressure Sodium, and ionic material NaCl are shown to demonstrate the accuracy and flexibility of the method in its practical applications. The method serves as a promising approach for large-scale automated material screening in which the melting temperature is a design criterion.

We present in detail two examples of refractory materials. First, we demonstrate how key material properties that provide guidance in the design of refractory materials can be accurately determined via ab initio thermodynamic calculations in conjunction with experimental techniques based on synchrotron X-ray diffraction and thermal analysis under laser-heated aerodynamic levitation. The properties considered include melting point, heat of fusion, heat capacity, thermal expansion coefficients, thermal stability, and sublattice disordering, as illustrated in a motivating example of lanthanum zirconate (La₂Zr₂O₇). The close agreement with experiment in the known but structurally complex compound La₂Zr₂O₇ provides good indication that the computation methods described can be used within a computational screening framework to identify novel refractory materials. Second, we report an extensive investigation into the melting temperatures of the Hf-C and Hf-Ta-C systems using ab initio calculations. With melting points above 4000 K, hafnium carbide (HfC) and tantalum carbide (TaC) are among the most refractory binary compounds known to date. Their mixture, with a general formula Ta_xHf_1-xC_y, is known to have a melting point of 4215 K at the composition Ta₄HfC₅, which has long been considered as the highest melting temperature for any solid. Very few measurements of melting point in tantalum and hafnium carbides have been documented, because of the obvious experimental difficulties at extreme temperatures. The investigation lets us identify three major chemical factors that contribute to the high melting temperatures. Based on these three factors, we propose and explore a new class of materials, which, according to our ab initio calculations, may possess even higher melting temperatures than Ta-Hf-C. This example also demonstrates the feasibility of materials screening and discovery via ab initio calculations for the optimization of "higher-level" properties whose determination requires extensive sampling of atomic configuration space.

Development and Application of Embedding Methods for the Simulation of Large Chemical Systems

The high computational cost of correlated wavefunction theory (WFT) calculations has motivated the development of numerous methods to partition the description of large chemical systems into smaller subsystem calculations. For example, WFT-in-DFT embedding methods facilitate the partitioning of a system into two subsystems: a subsystem A that is treated using an accurate WFT method, and a subsystem B that is treated using a more efficient Kohn-Sham density functional theory (KS-DFT) method. Representation of the interactions between subsystems is non-trivial, and often requires the use of approximate kinetic energy functionals or computationally challenging optimized effective potential calculations; however, it has recently been shown that these challenges can be eliminated through the use of a projection operator. This dissertation describes the development and application of embedding methods that enable accurate and efficient calculation of the properties of large chemical systems.

Chapter 1 introduces a method for efficiently performing projection-based WFT-in-DFT embedding calculations on large systems. This is accomplished by using a truncated basis set representation of the subsystem A wavefunction. We show that naive truncation of the basis set associated with subsystem A can lead to large numerical artifacts, and present an approach for systematically controlling these artifacts.

Chapter 2 describes the application of the projection-based embedding method to investigate the oxidative stability of lithium-ion batteries. We study the oxidation potentials of mixtures of ethylene carbonate (EC) and dimethyl carbonate (DMC) by using the projection-based embedding method to calculate the vertical ionization energy (IE) of individual molecules at the CCSD(T) level of theory, while explicitly accounting for the solvent using DFT. Interestingly, we reveal that large contributions to the solvation properties of DMC originate from quadrupolar interactions, resulting in a much larger solvent reorganization energy than that predicted using simple dielectric continuum models. Demonstration that the solvation properties of EC and DMC are governed by fundamentally different intermolecular interactions provides insight into key aspects of lithium-ion batteries, with relevance to electrolyte decomposition processes, solid-electrolyte interphase formation, and the local solvation environment of lithium cations.

The investigation of mechanisms in environmental catalysis using time-resolved methods

The formation of nitrogen oxides (NOx) during a combustion process is difficult to avoid because of the large exotherm and the consequent problem of avoiding local high-temperature spikes. Consequently, for many applications, such as for automotive power generation, there will be a continuing need to use catalytic after-treatment to reduce harmful emissions. The investigation of the mechanisms of the key catalytic reactions in environmental catalysis can provide an insight into the action of the catalyst, and time-resolved methods offer a powerful means to study these processes under realistic conditions. The use of Temporal Analysis of Products (TAP) and Steady State Isotopic Transient Kinetic Analysis (SSITKA) methods to investigate the reduction of NOx under various experimental conditions is described. From a detailed analysis of the SSITKA profiles, it is shown that at low temperatures the mechanism for the formation of N-2 and N2O from NO may differ from the conventional high-temperature mechanism. This is supported by density functional theory calculations, which show that the barrier to the formation of N2O from the reaction of N(ads) and NO(ads) may be too high to allow this process to occur at low temperatures. The alternative reaction of NO(ads) + NO(ads) = N2O(g) + O(ads) is shown to be much more favorable and is consistent with the SSITKA analysis. The remarkable effect of hydrogen as a reductant at low temperatures is described, and alternative interpretations of the role of hydrogen are discussed.