Density functional theory of transition metal phthalocyanines,II: electronic structure of MnPc and FePc—symmetry and symmetry breaking

We present a two-part systematic density functional theory (DFT) study of the electronic structure of selected transition metal phthalocyanines. We use a semi-local generalized gradient approximation (GGA) functional, as well as several hybrid exchange-correlation functionals, and compare the results to experimental photoemission data. Here, we study the intermediate spin systems MnPc and FePc. We show that DFT calculations of these systems are extremely sensitive to the choice of functional and basis set with respect to the obtained electronic configuration and to symmetry breaking. Interestingly, all simulated spectra are in good agreement with experiment despite the differences in the underlying electronic configurations.     

Density functional study of α-graphyne derivatives: Energetic stability,atomic and electronic structure

The energetic stability, atomic and electronic structures of α-graphyne and its derivatives (α-GYs) with extended carbon chains were investigated by density functional (DF) calculations in this work. The studied α-GYs consist of hexagon carbon rings sharing their edges with carbon atoms N=1–10. The structure and energy analyses show that α-GYs with even-numbered carbon chains have alternating single and triple C–C bonds (polyyne), energetically more stable than those with odd-numbered carbon chains possessing continuous double C–C bonds (polycumulene). The calculated electronic structures indicate that α-GYs can be either metallic (odd N) or semiconductive (even N) depending on the parity of number of atoms on hexagon edges despite the edge length. The semiconducting α-graphyne derivatives are found to possess Dirac cones (DC) with small direct band gaps 2–40 meV and large electron velocities 0.554×10⁶–0.671×10⁶ m/s, 70–80% of that of graphene. Our DF studies suggest that introducing sp carbon atoms into the hexagon edges of graphene opens up an avenue to switch between metallic and DC electronic structures via tuning the parity of the number of hexagon edge atoms.     

Density functional theory study on the structure and properties of Si3N4 clusters

The hybrid density functional theory B3LYP with basis sets 6-31G＊ has been used to study on the equilibrium geometries and electronic structures of possible isomers of Si3N4 clusters. 24 possible isomers are obtained. The most stable isomer of Si3N4 is a 3D structure with 7 Si-N bonds and 2 N-N bonds that could beformed by 3 quadrangles. The bond properties of the most stable isomer was analyzed by using natural bond orbital method （NBO）, the results suggest that the charges on Si and N atoms in Si-N bonds are quite large, so theinteraction of N-Si atoms in Si3N4 cluster is of strongly electric interaction. The primary IR and Raman vibrational frequency located at 1033.40 cm^-1, 473.63 cm^-1 respectively. The polarizabilities and hyperpolarizabilities of the most stable isomer are also analyzed.     

Quantum-size effects in the electronic structure of low-dimensional metallic systems

By angle-resolved photoemission the electronic structure of quantum films of Mg, Ag, and Au has been compared on W(110) and Mo(110) substrates which are structurally and electronically very similar but differ in atomic number. In all cases, substrate-induced states with characteristic dispersions are observed in the region of a bulk band gap of the substrate. Based on the comparison between Mo and W, we can exclude that previously observed Mg states are spin-orbit split but observe a spin-orbit splitting in Ag and Au monolayers. This splitting is mainly caused by the substrate because it does not differ much between Ag and Au overlayers despite the large difference in atomic number.     

Structures,stabilities, and electronic properties of F-doped Sin （n=1～12） clusters：Density functional theory investigation

The geometries, stabilities, and electronic properties of FSin （n=1～12） clusters are systematically investigated by using first-principles calculations based on the hybrid density-functional theory at the B3LYP/6-311G level. The geometries are found to undergo a structural change from two-dimensional to three-dimensional structure when the cluster size n equals 3. On the basis of the obtained lowest-energy geometries, the size dependencies of cluster properties, such as averaged binding energy, fragmentation energy, second-order energy difference, HOMO–LUMO （highest occupied molecular orbital–lowest unoccupied molecular orbital） gap and chemical hardness, are discussed. In addition, natural population analysis indicates that the F atom in the most stable FSin cluster is recorded as being negative and the charges always transfer from Si atoms to the F atom in the FSin clusters.     

Density functional theory of a Gay—Berne film between aligning walls

It is shown that the simple Onsager second-virial approximation of density functional theory can successfully describe the orientational structure of a Gay-Berne film confined between aligning plates. The theory takes as input the density profile as determined by computer simulation and semi-quantitatively reproduces the variation in the nematic order parameter throughout the film, at different temperatures and for different surface potential strengths, without any adjustable parameters. In this context the validity of an earlier analytical approach is discussed where the density, order parameter and tilt angle profiles were assumed to be step functions.     

Density functional investigation of structural and electronic properties of small bimetallic silver–gold clusters

Structural and electronic properties of bimetallic silver–gold clusters up to eight atoms are investigated by the density functional theory using Wu and Cohen generalized gradient approximation functional. By substitution of Ag and Au atoms, in the optimized lowest energy structures of pure gold and silver clusters, we determine the ground state conformations of the bimetallic silver–gold ones. We reveal that Ag atoms prefer internal positions whereas Au atoms prefer exposed ones favoring charge transfer from Ag to Au atoms. For each size and composition, binding energy, HOMO–LUMO gap, magnetic moment, vertical ionization potential, electron affinity and chemical hardness were calculated. On increasing the size of the cluster by varying number of Ag atoms with fixed number of Au ones, vertical ionization potential and electron affinity show obvious odd–even oscillations consistent with the pure Ag and Au clusters. Au atoms inclusion in the cluster increases the binding energy and vertical ionization potential, indicating higher stability as the number of Au atoms grows. The variation of chemical hardness with the composition in a cluster with the same size shows peaks when the number of Ag atoms is greater than or equal to Au ones, corresponding to transition from planar to tri-dimensional structures. For clusters with even number of atoms, the peaks indicate that the clusters with the same number of Ag and Au atoms are the most stable ones. Analyzing the density of states, we found that increasing the concentration of Ag atoms affects the energy separation between the HOMO and the low lying occupied states.     

Consistent LDA’ + DMFT approach to the electronic structure of transition metal oxides: Charge transfer insulators and correlated metals

We discuss the recently proposed LDA’ + DMFT approach providing a consistent parameter-free treatment of the so-called double counting problem arising within the LDA + DMFT hybrid computational method for realistic strongly correlated materials. In this approach, the local exchange-correlation portion of the electron-electron interaction is excluded from self-consistent LDA calculations for strongly correlated electronic shells, e.g., d-states of transition metal compounds. Then, the corresponding double-counting term in the LDA’ + DMFT Hamiltonian is consistently set in the local Hartree (fully localized limit, FLL) form of the Hubbard model interaction term. We present the results of extensive LDA’ + DMFT calculations of densities of states, spectral densities, and optical conductivity for most typical representatives of two wide classes of strongly correlated systems in the paramagnetic phase: charge transfer insulators (MnO, CoO, and NiO) and strongly correlated metals (SrVO₃ and Sr₂RuO₄). It is shown that for NiO and CoO systems, the LDA’ + DMFT approach qualitatively improves the conventional LDA + DMFT results with the FLL type of double counting, where CoO and NiO were obtained to be metals. Our calculations also include transition-metal 4s-states located near the Fermi level, missed in previous LDA + DMFT studies of these monoxides. General agreement with optical and the X-ray experiments is obtained. For strongly correlated metals, the LDA’ + DMFT results agree well with the earlier LDA + DMFT calculations and existing experiments. However, in general, LDA’ + DMFT results give better quantitative agreement with experimental data for band gap sizes and oxygen-state positions compared to the conventional LDA + DMFT method.     

Angle-resolved photoemission of chain-like molecules: the electronic band structure of sexithiophene and sexiphenyl

Here we report the electronic π-band structure of sexithiophene obtained from 6T(010) oriented films. The angle-resolved valence band photoemission results taken parallel and perpendicular to the molecular axis are compared to those of sexiphenyl and interpreted in terms of intra- and inter-molecular band dispersion. We show that the strong photoemission intensity variations with emission angle parallel to the molecular axis are well reproduced by the Fourier transforms of the molecular orbitals of the isolated molecules. These results imply that ARUPS can yield quite detailed information about molecular geometry, both in terms of molecular orientation and internal structure.     

Density functional theory study on the structure and properties of Si3N4clusters

The hybrid density functional theory B3LYP with basis sets 6-31G· has been used to study on the equilibrium geometries and electronic structures of possible isomers of Si3N4 clusters. 24 possible isomers are obtained. The most stable isomer of Si3N4 is a 3D structure with 7 Si-N bonds and 2 N-N bonds that could be formed by 3 quadrangles. The bond properties of the most stable isomer was analyzed by using natural bond orbital method (NBO), the results suggest that the charges on Si and N atoms in Si-N bonds are quite large, so the interaction of N-Si atoms in Si3N4 cluster is of strongly electric interaction. The primary IR and Raman vibrational frequency located at 1033.40 cm-1, 473.63 cm-1 respectively. The polarizabilities and hyperpolarizabilities of the most stable isomer are also analyzed.     

Surface electronic structure of Ge(001)2 × 1: experiment and theory

Detailed studies of the surface electronic structure of Ge(001)2 × 1 have been performed with angle-resolved photoemission. Five surface structures are identified in the spectra and their initial-energy dispersions E(k_∥) are presented along the [010] direction. Employing the local density approximation and Green's functions, the surface band structure along this direction has been calculated self-consistently for an asymmetric dimer model of the 2 × 1-reconstructed surface. Four of the five experimental surface structures are identified with calculated surface states or resonances, i.e. the dangling-bond state and two different back-bond resonances. Excellent agreement is obtained between the experimental and theoretical dispersions for these surface electronic bands.     

Algebraic connectivity analysis in molecular electronic structure theory I: coulomb potential,tensor connectivity, ε-approximation

In this (first) paper we attempt to generalize the notion of tensor connectivity, subsequently studying how this property is affected in different tensorial operations. We show that the often implied corollary of the linked diagram theorem, namely individual size-extensivity of arbitrary connected closed diagrams, can be violated in Coulomb systems. In particular, the assumption of the existence of localized Hartree–Fock orbitals is generally incompatible with the individual size-extensivity of connected closed diagrams when the interaction tensor is generated by the true two-body part of the electronic Hamiltonian. Thus, in general, size-extensivity of a many-body method may originate in specific cancellations of super-extensive quantities, breaking the convenient one-to-one correspondence between the connectivity of arbitrary many-body equations and the size-extensivity of the expectation values evaluated by those equations (for example, when certain diagrams are discarded from the method). Nevertheless, assuming that many-body equations are evaluated for a stable many-particle system, it is possible to introduce a workaround, called the ε-approximation, which restores the individual size-extensivity of an arbitrary connected closed diagram, without qualitatively affecting the asymptotic behavior of the computed expectation values. No assumptions concerning the periodicity of the system and its strict electrical neutrality are made.     

Density functional theory studies on the electronic, structural, phonon dynamical and thermo-stability properties of bicarbonates MHCO(3), M?=??Li, Na, K

The structural, electronic, phonon dispersion and thermodynamic properties of MHCO(3) (M = Li, Na, K) solids were investigated using density functional theory. The calculated bulk properties for both their ambient and the high-pressure phases are in good agreement with available experimental measurements. Solid phase LiHCO(3) has not yet been observed experimentally. We have predicted several possible crystal structures for LiHCO(3) using crystallographic database searching and prototype electrostatic ground state modeling. Our total energy and phonon free energy (F(PH)) calculations predict that LiHCO(3) will be stable under suitable conditions of temperature and partial pressures of CO(2) and H(2)O. Our calculations indicate that the [Formula: see text] groups in LiHCO(3) and NaHCO(3) form an infinite chain structure through O?H?O hydrogen bonds. In contrast, the [Formula: see text] anions form dimers, [Formula: see text], connected through double hydrogen bonds in all phases of KHCO(3). Based on density functional perturbation theory, the Born effective charge tensor of each atom type was obtained for all phases of the bicarbonates. Their phonon dispersions with the longitudinal optical-transverse optical splitting were also investigated. Based on lattice phonon dynamics study, the infrared spectra and the thermodynamic properties of these bicarbonates were obtained. Over the temperature range 0-900 K, the F(PH) and the entropies (S) of MHCO(3) (M =Li, Na, K) systems vary as F(PH)(LiHCO(3)) > F(PH)(NaHCO(3)) > F(PH)(KHCO(3)) and S(KHCO(3)) > S(NaHCO(3)) > S(LiHCO(3)), respectively, in agreement with the available experimental data. Analysis of the predicted thermodynamics of the CO(2) capture reactions indicates that the carbonate/bicarbonate transition reactions for Na and K could be used for CO(2) capture technology, in agreement with experiments.     

Density functional theory study on LaNi4.5Al0.5 hydride phase： electronic properties and sites occupation＊

In this paper the crystal structure, electronic structure and hydrogen site occupation of LaNi4.5Al0.5Hy hydride phase （y = 5.0, 6.0） have been investigated by using full-potential linearized augmented plane wave method. The hydrogen atoms were found to prefer the 6m, 12o and 12n sites, while no 4h sites were occupied. A narrowed Ni-d band is found due to the lattice expansion, the total density of states at EF increases with y, which indicates that the compounds become less stable. The interaction between Al and Ni, H plays a dominant role in the stability of LaNi4.5Al0.5 hydride phase. The smaller the shift of EF towards the higher energy region, the more stable the compounds will be. The obtained results are compared with experimental data and discussed in the light of previous works.     

Computationally simple,analytic, closed form solution of the Coulomb self-interaction problem in Kohn–Sham density functional theory

We have developed and tested in terms of atomic calculations an exact, analytic and computationally simple procedure for determining the functional derivative of the exchange energy with respect to the density in the implementation of the Kohn–Sham formulation of density functional theory (KS-DFT), providing an analytic, closed-form solution of the self-interaction problem in KS-DFT. We demonstrate the efficacy of our method through ground-state calculations of the exchange potential and energy for atomic He and Be atoms, and comparisons with experiment and the results obtained within the optimized effective potential (OEP) method.     

Calculation of the structure and IR spectrum of methyl-β-D-glucopyranoside by density functional theory

Structural-dynamic models of methyl-β-D-glucopyranoside have been constructed by a density functional method using a B3LYP functional in bases 6-31G(d) and 6-31+G(d,p). Energies have been minimized. Structures, dipole moments, polarizabilities, frequencies of normal modes in the harmonic approximation, and the intensity distribution in the molecular IR spectrum have been calculated. The calculation results have been compared with both the experimental spectra of methyl-β-D-glucopyranoside in the region 400–3700 cm^–1 and data obtained within the framework of an approach that uses the classical valence-force method to calculate normal mode frequencies and the quantum-chemical CNDO/2 technique to calculate the electronic structure.     

Mechanical effects in the n → s transition of work-hardened metal and alloy superconductors

In the superconducting transitions, the dissolution of filamentary forms of intermediate phase, retained at edge-type dislocations to temperatures somewhat below (ns) and above (sn) the critical temperature, induces (+ or –) thermal spikes at the dislocation cores. In a recent model, these have been shown to trigger transients in creep, and to give rise to yield-point effects in tensile tests, such as are observed in the ns transition in metals and undeformed alloys. The model is now extended to cover corresponding observations in type I work-hardened metal and alloy superconductors.     

Parallelization of analytical Hartree—Fock and density functional theory Hessian calculations. Part I: parallelization of coupled-perturbed Hartree—Fock equations

Solving the coupled-perturbed Hartree-Fock (CPHF) equations is the most time consuming part in the analytical computation of second derivatives of the molecular energy with respect to the nuclei. This paper describes a unique parallelization approach for solving the CPHF equations. The computational load is divided by the nuclear perturbations and distributed evenly among the computing nodes. The parallel algorithm is scalable with respect to the size of the molecule, i.e. the larger the molecule, the greater the parallel speedup. The memory storage requirements are also distributed among the processors, with little communication among the processors. The method is implemented in the Q-Chem software package and its performance is discussed. This work represents the first step in a research project to parallelize analytical frequency calculations at Hartree-Fock and density functional theory levels.     

Symmetry dependence of x—ray absorption near—edge structure at rhe metal K edge of 3d transition metal compounds

The pre-edge features in system with “even” symmetry,apart from quadrupolar tranition contribution,are mainly dipolar in character,associated with the existence of unoccupied states made up of mixed cation np with higher-neighboring cation-(n-l)d orbitals,and reflect the density of states due to the medium-range order of the system,while in “ood” symmetry materials these pre-edge features are the result of a transition from the 1s to final denisty of states of p symmetry due to an unsymmetrical mixing of the ligand wave functions with the central cation d orbitals.In the latter case,they contain not only the p but also the d base of orbitals,similar to a tetrahedral configuration.These results are validated for Fe as a photoabsorber by comparing x-ray absorption near-edge spectra of Fe2SiO4(fayalite) and Fe2O3(hematite) to ab initio full multiple scattering calculations at the Fe K edge,but pertain to all systems containing sixfold-coordinated cations.