Au12M (M=Na,Mg, Al,Si, P,S, Cl)团簇的结构、稳定性和电子性质

采用基于密度泛函理论的第一性原理方法系统地研究了Au12M(M=Na,Mg,Al,Si,P,S,Cl)团簇的结构、稳定性和电子性质.对团簇的平均结合能、镶嵌能、垂直离化势、最高占据分子轨道(HOMO)和最低未占据分子轨道(LUMO)的能级差、电荷布居分析、自然键轨道(NBO)进行了计算和讨论.对于Au12M(M=Na,Mg,Al)团簇,它们形成了内含M原子的最稳定的笼状结构.然而对于Au12M(M=Si,P,S,Cl)团簇,它们却形成了以M元素为顶点的稳定锥形结构.在这些团簇中发现Au12S团簇相对是最稳定的,这是由于Au12S团簇形成了稳定的满壳层的电子结构.自然电荷布居分析表明:对于所有的Au12M(M=Na,Mg,Al,Si,P,S,Cl)团簇电荷总是从Au原子转向M原子.自然键轨道和HOMO分析表明Au12M团簇中发生了Au原子的s-d轨道和M原子的p轨道间的杂化现象.     

Dynamic relaying properties of a β-turn peptide in long-range electron transfer

The relay stations play a significant role in long-range charge hopping transfer in proteins. Although studies have clarified that many more protein structural motifs can function as relays in charge hopping transfers by acting as intermediate charge carriers, the relaying properties are still poorly understood. In this work, taking a β-turn oligopeptide as an example, we report a dynamic character of a relay with tunable relaying properties using the density functional theory calculations. Our main finding is that a β-turn peptide can serve as an effective electron relay in facilitating long-range electron migration and its relay properties is vibration-tunable. The vibration-induced structural transient distortions remarkably affect the lowest occupied molecular orbital (LUMO) energy, vertical electron affinity and electron-binding mode of the β-turn oligopeptide and the singly occupied molecular orbital (SOMO) energy of the corresponding electron adduct and thus the relaying properties. Different vibration modes lead to different structural distortions and thus have different effects on the relaying properties and ability of the β-turn peptide. For the relaying properties, there approximately is a linear negative correlation of electron affinity with the LUMO energy of the β-turn or the SOMO energy of its electron adduct. Besides, such relaying properties also vary in the vibration evolution process, and the electron-binding modes may be tunable. As an important addition to the known static charge relaying properties occurring in various protein structural motifs, this work reports the dynamic electron-relaying characteristics of a β-turn oligopeptide with variable relaying properties governed by molecular vibrations which can be applied to different proteins in mediating long-range charge transfers. Clearly, this work reveals molecular vibration effects on the electron relaying properties of protein structural motifs and provides new insights into the dynamics of long-range charge transfers in proteins. © 2018 Wiley Periodicals, Inc.     

Molecular electronegativity in density functional theory (VI)

Based on the density functional theory and partitioning the molecular electron density ρ (r) into atomic electronic densities and bond electronic densities, the expressions of the total molecular energy and the “effective electronegativity” of an atom or a bond in a molecule are obtained. The atom-bond electronegativity equalization model is then proposed for the direct calculation of the total molecular energy and the charge distribution of large molecules. Practical calculations show that the atom-bond electronegativity equalization model can reproduce the correspondingab initio values of the total molecular energies and charge distributions for a series of large molecules with a very satisfactory accuracy.     

Computational study of iron(II) systems containing ligands with nitrogen heterocyclic groups

Density functional theory (DFT) calculations show the higher energy HOMO (highest occupied molecular orbital) orbitals of four iron(II) diimine complexes are metal centered and the lower energy LUMO (lowest unoccupied molecular orbitals) are ligand centered. The energy of the orbitals correlates with electrochemical redox potentials of the complexes. Time-dependent density functional theory (TDDFT) calculations reveal ligand centered (LC) and metal-to-ligand charge transfer (MLCT) at higher energy than experimentally observed. TDDFT calculations also reveal the presence of d-d transitions which are buried under the MLCT and LC transitions. The difference in chemical and photophysical behavior of the iron complexes compared to that of their ruthenium analogues is also addressed.     

Stability of atomic oxygen chemisorption on Pt‐alloy surfaces

A density functional theory calculation is used to investigate the atomic oxygen (O) stability over platinum (Pt) and Pt‐based alloy surfaces. Here, the stability is connected with the preferential adsorption sites for O chemisorptions and the adsorption energy. Thus, the interaction mechanism between atomic O and metal surfaces is studied by using charge transfer analysis. In this present paper, atomic structure and binding energy of oxygen adsorption on the Pt(111) are in a very good agreement with experiment and previous density functional theory calculations. Furthermore, we obtained that the addition of ruthenium (Ru) and molybdenum (Mo) on the pure Pt surface enhances the adsorption energy. Our charge transfer analysis shows that the largest charge transfer contributing to the metal‐O bonding formation is observed in the case of O/PtRuMo surface followed by O/PtRu surface. This is in consistency with metal d‐orbital characteristic, where Mo has much more empty d‐orbital than Ru in correspondence to accept electrons from atomic oxygen. Copyright © 2016 John Wiley & Sons, Ltd.     

Towards a force field based on density fitting

Total intermolecular interaction energies are determined with a first version of the Gaussian electrostatic model (GEM-0), a force field based on a density fitting approach using s-type Gaussian functions. The total interaction energy is computed in the spirit of the sum of interacting fragment ab initio (SIBFA) force field by separately evaluating each one of its components: electrostatic (Coulomb), exchange repulsion, polarization, and charge transfer intermolecular interaction energies, in order to reproduce reference constrained space orbital variation (CSOV) energy decomposition calculations at the B3LYP/aug-cc-pVTZ level. The use of an auxiliary basis set restricted to spherical Gaussian functions facilitates the rotation of the fitted densities of rigid fragments and enables a fast and accurate density fitting evaluation of Coulomb and exchange-repulsion energy, the latter using the overlap model introduced by Wheatley and Price [Mol. Phys. 69, 50718 (1990)]. The SIBFA energy scheme for polarization and charge transfer has been implemented using the electric fields and electrostatic potentials generated by the fitted densities. GEM-0 has been tested on ten stationary points of the water dimer potential energy surface and on three water clusters (n = 16,20,64). The results show very good agreement with density functional theory calculations, reproducing the individual CSOV energy contributions for a given interaction as well as the B3LYP total interaction energies with errors below kBT at room temperature. Preliminary results for Coulomb and exchange-repulsion energies of metal cation complexes and coupled cluster singles doubles electron densities are discussed.     

Density functional theory studies of spin, charge, and orbital ordering in YBaT2O5 (T = Mn, Fe, Co)

Spin, charge, and orbital orderings are influenced by electron/hole doping, cation radii, oxygen stoichiometry, temperature, magnetic field, and so on. In order to understand the role of electron/hole doping, we have studied variations in spin, charge, and orbital ordering in terms of d-band filling for YBaT 2O 5 (T = Mn, Fe, Co). The calculations were performed using density functional theory as implemented in the full-potential linearized augmented-plane-wave method. We have carried out calculations for nonmagnetic, ferromagnetic, and antiferromagnetic configurations. A ferrimagnetic ground state was established for YBaMn 2O 5, whereas YBaFe 2O 5 and YBaCo 2O 5 have antiferromagnetic ground states; all of these results are in agreement with experimental findings. The effects of spin-orbit coupling, the Hubbard U parameter, and orbital polarization on the magnetic properties were also analyzed. The electronic band characteristics were analyzed using total as well as site- and orbital-projected densities of states. Inclusion of spin-orbit coupling and Coulomb correlation effects in the calculations was found to be important in order to reproduce the experimentally established semiconducting behaviors of YBaFe 2O 5 and YBaCo 2O 5. In order to quantify the charges at each atomic site, we made use of the Bader "atom-in-molecule" concept and Born effective-charge (BEC) analyses. The structural optimizations and BEC tensor calculations were performed using the VASP-PAW method. The different types of charge and orbital orderings in these compounds were visualized using the energy-projected density matrices of the d electrons. Substantial differences in ordering patterns with respect to d-band filling emerged. Ordering of the d z (2) orbital of Mn in YBaMn 2O 5 gave rise to G-type ferrimagnetic spin ordering along the c direction and checkerboard-type charge ordering, whereas ordering of the d x (2) - y (2) orbital of Fe in YBaFe 2O 5 caused Wollan-Koehler G-type antiferromagnetic spin ordering along the b direction and stripe-type charge ordering. Similarly, a complex pattern of orbital ordering in YBaCo 2O 5 activated spin and charge orderings similar to those in YBaFe 2O 5.     

Orbital contributions to the molecular charge and energy density distributions in Co2(CO)8

For Co2(CO)8, the representative of a whole class of bridged cobalt complexes, the 18-electron rule predicts a direct metal-metal bond in addition to the metal-bridge bonds. By intuition, this bond should have bent-bond character. However, it is well-known from charge density analyses that no bond critical point exists in the corresponding spatial region. Otherwise, the energy density distribution points to a certain stabilizing contribution of this local area to the total molecular energy. It is shown that a partitioning of the total charge and energy densities into orbital contributions can lead to a deeper insight into complex bonding properties.     

Structural and spectroscopic characterization of 2,3-difluorobenzoic acid and 2,4-difluorobenzoic acid with experimental techniques and quantum chemical calculations

In this study, the molecular conformation, vibrational and electronic transition analysis of 2,3-difluorobenzoic acid and 2,4-difluorobenzoic acid (C7H4F2O2) were presented using experimental techniques (FT-IR, FT-Raman and UV) and quantum chemical calculations. FT-IR and FT-Raman spectra in solid state were recorded in the region 4000-400 cm(-1) and 4000-5 cm(-1), respectively. The UV absorption spectra of the compounds that dissolved in ethanol were recorded in the range of 200-800 nm. The structural properties of the molecules in the ground state were calculated using density functional theory (DFT) and second order M?ller-Plesset perturbation theory (MP2) employing 6-311++G(d,p) basis set. Optimized structure of compounds was interpreted and compared with the earlier reported experimental values. The scaled vibrational wavenumbers were compared with experimental results. The complete assignments were performed on the basis of the experimental data and total energy distribution (TED) of the vibrational modes, calculated with scaled quantum mechanics (SQM) method. A study on the electronic properties, such as absorption wavelength, excitation energy, dipole moment and frontier molecular orbital energy, were performed by time dependent DFT (TD-DFT) approach. Based on the UV spectra and TD-DFT calculations, the electronic structure and the assignments of the absorption bands of steady compounds were discussed. The calculated HOMO and LUMO energies show that charge transfer occurs within the molecules.     

二苯二硫和二苄二硫润滑性能的量子化学研究

The molecular geometries optimization and electronic structures of diphenyl disulfide (DPDS) and dibenzyl disulfide (DBDS) compounds were investigated by density functional theory (DFT) and ab initio method at the 6-31G basis set level. The active atoms and bonds of reaction were provided by frontier molecular orbital theory. The molecular orbital parameters of DPDS and DBDS compounds and iron atom cluster were calculated by using density functional theory. The interaction pattern between the organic disulfide compounds and iron atom cluster was discussed based on the approximate rule of orbital energy. Some parameters characterizing the action strength between the organic disulfide compounds and iron atom cluster, including the bonding strength, reactive strength and static action strength, were analyzed by using frontier electron density, super de-localizability, net atomic charge and the interaction energy of chemical adsorption as criteria. The results indicate that S-S chemical bond and C-S chemical bond of the compounds are inclined to be broken when DPDS and DBDS interact with the metal. The anti wear ability order of DPDS and DBDS compounds is DPDS>DBDS, and the extreme pressure ability order of DPDS and DBDS compounds is DBDS>DPDS, and the prediction results based on quantum chemistry calculations are in good accordance with the friction and wear test results.     

Dicyano and pyridine derivatives of β-carotene: synthesis and vibronic, electronic, and photophysical properties

Density functional theory and time-dependent density functional theory calculations provide pictures of the molecular orbitals involved in the ground and excited states of two cyano derivatives of 8'-apo-β-caroten-8'-al synthesized via an acid-base-catalyzed Knoevenagel condensation reaction. Population analysis shows that the symmetry-allowed transition, S(0) ((1)A(g)) → S(2) ((1)B(u)) based on the C(2h) symmetry is a HOMO (highest occupied molecular orbital) to LUMO (lowest unoccupied molecular orbital) π → π* transition with electron densities located mostly on the polyene chain. Calculated and actual steady-state absorption spectra show similar features with low-energy peak maxima between 550 and 600 nm.     

Intermolecular electrostatic energies using density fitting

A method is presented to calculate the electron-electron and nuclear-electron intermolecular Coulomb interaction energy between two molecules by separately fitting the unperturbed molecular electron density of each monomer. This method is based on the variational Coulomb fitting method which relies on the expansion of the ab initio molecular electron density in site-centered auxiliary basis sets. By expanding the electron density of each monomer in this way the integral expressions for the intermolecular electrostatic calculations are simplified, lowering the operation count as well as the memory usage. Furthermore, this method allows the calculation of intermolecular Coulomb interactions with any level of theory from which a one-electron density matrix can be obtained. Our implementation is initially tested by calculating molecular properties with the density fitting method using three different auxiliary basis sets and comparing them to results obtained from ab initio calculations. These properties include dipoles for a series of molecules, as well as the molecular electrostatic potential and electric field for water. Subsequently, the intermolecular electrostatic energy is tested by calculating ten stationary points on the water dimer potential-energy surface. Results are presented for electron densities obtained at four different levels of theory using two different basis sets, fitted with three auxiliary basis sets. Additionally, a one-dimensional electrostatic energy surface scan is performed for four different systems (H2O dimer, Mg2+-H2O, Cu+-H2O, and n-methyl-formamide dimer). Our results show a very good agreement with ab initio calculations for all properties as well as interaction energies.     

Field emission properties of N-doped capped single-walled carbon nanotubes: a first-principles density-functional study

The geometrical structures and field emission properties of pristine and N-doped capped (5,5) single-walled carbon nanotubes have been investigated using first-principles density-functional theory. The structures of N-doped carbon nanotubes are stable under field emission conditions. The calculated work function of N-doped carbon nanotube decreases drastically when compared with pristine carbon nanotube, which means the enhancement of field emission properties. The ionization potentials of N-doped carbon nanotubes are also reduced significantly. The authors analyze the field emission mechanism in terms of energy gap between the lowest unoccupied molecular orbital and the highest occupied molecular orbital, Mulliken charge population, and local density of states. Due to the doping of nitrogen atom, the local density of states at the Fermi level increases dramatically and donor states can be observed above the Fermi level. The authors' results suggest that the field emission properties of carbon nanotubes can be enhanced by the doping of nitrogen atom, which are consistent with the experimental results.     

A DFT/TDDFT study of porphyrazines and phthalocyanine oxo‐titanium derivatives as potential dyes in solar cells

Density functional theory and time dependent density functional theory calculations at the level of LDA/BP86/TZ2P were performed systematically on several Ti(IV) complexes of porphyrazines and one phthalocyanine. We performed an analysis of the frontier molecular orbitals of the ground state electronic structures and also discuss in particular the good concordance of our results with the experimental data, which affords to predict the geometrical and optical properties of new complexes ( 3 , 4 , and 7 ). We also emphasize the characterization of the UV–vis absorption spectra and propose transitions that contribute to the Q and B bands. Some useful calculated properties in complexes 2 , 3 , and 7 , like: high light absorption in the visible region of the spectra, transitions involved in these bands with a determined direction, charge separation, bigger highest occupied molecular orbital–lowest unoccupied molecular orbital (HOMO‐LUMO) gaps than complexes 4 and 5 , and the energy of their LUMO orbitals (that are higher than the lowest energy level of the conduction band of the TiO₂) indicate that system complexes 2 , 3 , and 7 could act as light‐harvesting sensitizers for dye‐sensitized solar cells (DSCs). These proposals were made using a model of the previously experimentally known phthalocyanine, which was used as sensitizer in DSCs devices, comparing its electronic properties with the herein proposed sensitizers. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Electronic Structure Insights into the Solvation of Magnesium Ions with Cyclic and Acyclic Carbonates

A computational framework to rank the solvation behavior of Mg²⁺ in carbonates by using molecular dynamics simulations and density functional theory is reported. Based on the binding energies and enthalpies of solvation calculated at the M06‐2X/6‐311++G(d,p) level of theory and the free energies of solvation from ABF‐MD simulations, we find that ethylene carbonate (EC) and the ethylene carbonate:propylene carbonate (EC:PC) binary mixture are the best carbonate solvents for interacting with Mg²⁺. Natural bond orbital and quantum theory of atoms in molecules analyses support the thermochemistry calculations with the highest values of charge transfer, perturbative stabilization energies, electron densities, and Wiberg bond indices being observed in the Mg²⁺(EC) and Mg²⁺(EC:PC) complexes. The plots of the noncovalent interactions indicate that those responsible for the formation of Mg²⁺ carbonate complexes are strong‐to‐weak attractive interactions, depending on the regions that are interacting. Finally, density of state calculations indicate that the interactions between Mg²⁺ and the carbonate solvents affects the HOMO and LUMO states of all carbonate solvents and moves them to more negative energy values.     

Charge transfer in the TCNQ-sexithiophene complex

Molecular crystals from thiophene molecules can be doped with TCNQ-F4 molecules for use in all-organic optoelectronic and semiconductor devices. The charge transfer and the molecular orbital energy level formation in between these two organic molecules are investigated here by density functional theory calculations. The isolated molecules are calculated nonbonded and bonded together, forming a charge transfer complex (CTC). The relaxed structure of the complex shows essentially coplanar and centered molecules with the alpha-sexithiophene rings tilted alternatingly by 4.8 degrees. The bond formation of these molecules results in a charge transfer of approximately 0.4 e from the alpha-sexithiophene to the TCNQ-F4 molecule. The highest occupied molecular orbital-lowest unoccupied molecular orbital gap width is reduced as compared to the isolated molecules due to the newly formed orbitals in the CTC. Upon adsorption on a Au(111) surface, electrons are transferred onto the molecule complex, thereby causing the molecular levels to align asymmetric with respect to the charge neutrality level. The theoretical results for the single molecule and CTC layer are compared to experimental photoemission and scanning tunneling spectroscopy results.     

Sulfur K-edge X-ray absorption spectroscopy as a probe of ligand-metal bond covalency: metal vs ligand oxidation in copper and nickel dithiolene complexes

A combination of Cu L-edge and S K-edge X-ray absorption data and density functional theory (DFT) calculations has been correlated with 33S electron paramagnetic resonance superhyperfine results to obtain the dipole integral (Is) for the S 1s-->3p transition for the dithiolene ligand maleonitriledithiolate (MNT) in (TBA)2[Cu(MNT)2] (TBA= tetra-n-butylammonium). The results have been combined with the Is of sulfide derived from XPS studies to experimentally obtain a relation between the S 1s-->4p transition energy (which reflects the charge on the S atom, QSmol) and the dipole integral over a large range of QSmol. The results show that, for high charges on S, Is can vary from the previously reported Is values, calculated using data over a limited range of QSmol. A combination of S K-edge and Cu K- and L-edge X-ray absorption data and DFT calculations has been used to investigate the one-electron oxidation of [Cu(MNT)2]2- and [Ni(MNT)2]2-. The conversion of [Cu(MNT)2]2- to [Cu(MNT)2]- results in a large change in the charge on the Cu atom in the molecule (QCumol) and is consistent with a metal-based oxidation. This is accompanied by extensive charge donation from the ligands to compensate the high charge on the Cu in [Cu(MNT)2]- based on the increased S K-edge and decreased Cu L-edge intensity, respectively. In contrast, the oxidation of [Ni(MNT)2]2- to [Ni(MNT)2]- results in a small change in QNimol, indicating a ligand-based oxidation consistent with oxidation of a molecular orbital, psiSOMO (singly occupied molecular orbital), with predominant ligand character.     

半胱氨酸在碳钢与硫酸界面的缓蚀行为

Interfacial behavior of cysteine (Cys) between mild steel and sulfuric acid solution as a corrosion inhibitor has been studied with electrochemical AC (alternating current) and DC (direct current) techniques at (25.0±0.1) ℃. The AC impedance results were evaluated using equivalent circuits in which a constant phase element (CPE) has been replaced with double layer capacitance (Cdl) to represent the frequency distribution of experimental data. Changes in impedance parameters (charge transfer resistance and double layer capacitance) indicated that cysteine molecules acted by accumulating at the metal/solution interface. The fractional coverage of the metal surface (θ) was determined using AC impedance results and it was found that the adsorption of cysteine on the mild steel surface followed a Langmuir isothermmodel with a standard free energy of adsorption (⊿G0ads) of -35.1 kJ·mol-1.  To clarify the type of interaction between mild steel surface and cysteine molecules with a molecular orbital approach, electronic properties, such as, the highest occupied molecular orbital (HOMO) energy, the lowest unoccupied molecular orbital (LUMO) energy, and the frontier molecular orbital coefficients have been calculated. Energy gaps for the interaction of mild steel surface and cysteine molecules (ELUMOFe-EHOMOCys and ELUMOCys-EHOMOFe) were used to determine whether cysteine molecules acted as electron donors or electron acceptors when they interacted with the mild steel surface. The local reactivity was evaluated through the condensed Fukui indices. Theoretical calculations were carried out using the density functional theory (DFT) at B3LYP level with the 6-311++G(d,p) basis set for all atoms by Gaussian 03W program.     

Isoelectronic changes in energy of quark atoms and molecules via the Levy equation

Summary Z-transition state calculations based on the Levy equation suggest that the isoelectronic changes in energy of quark atoms,Q, (ordinary atoms with extra nuclear charge in units of ±1/3 and/or ±2/3) can be expressed quantitatively in terms of the electrostatic potential at the nucleus of an isoelectronic ordinary atom. Numerical tests within the local density functional theory are presented for the quark atoms of Li-F. Theab initio MO (molecular orbital) calculations using STO-5G basis on the C₂ molecule and its quark derivatives lead to similar conclusions.