The Kohn-Sham density of states and band gap of water: from small clusters to liquid water

Electronic properties of water clusters (H2O)(n), with n=2, 4, 8, 10, 15, 20, and 30 molecules were investigated by sequential Monte Carlo/density-functional theory (DFT) calculations. DFT calculations were carried out over uncorrelated configurations generated by Monte Carlo simulations of liquid water with a reparametrized exchange-correlation functional that reproduces the experimental information on the electronic properties (first ionization energy and highest occupied molecular orbital-lowest unoccupied molecular orbital gap) of the water dimer. The dependence of electronic properties on the cluster size (n) shows that the density of states (DOS) of small water clusters (n>10) exhibits the same basic features that are typical of larger aggregates, such as the mixing of the 3a1 and 1b1 valence bands. When long-ranged polarization effects are taken into account by the introduction of embedding charges, the DOS associated with 3a1 orbitals is significantly enhanced. In agreement with valence-band photoelectron spectra of liquid water, the 1b1, 3a1, and 1b2 electron binding energies in water aggregates are redshifted by approximately 1 eV relative to the isolated molecule. By extrapolating the results for larger clusters the threshold energy for photoelectron emission is 9.6+/-0.15 eV (free clusters) and 10.58+/-0.10 eV (embedded clusters). Our results for the electron affinity (V0=-0.17+/-0.05 eV) and adiabatic band gap (E(G,Ad)=6.83+/-0.05 eV) of liquid water are in excellent agreement with recent information from theoretical and experimental works.     

Performance of numerical basis set DFT for aluminum clusters

We have investigated and compared the ability of numerical and Gaussian-type basis sets combined with density functional theory (DFT) to accurately describe the geometries, binding energies, and electronic properties of aluminum clusters, Al12XHn (X = Al, Si; n = 0, 1, 2). DFT results are compared against high-level benchmark calculations and experimental data where available. Properties compared include geometries, binding energies, ionization potentials, electron affinities, and HOMO-LUMO gaps. Generally, the PBE functional with the double numerical basis set with polarization (DNP) performs very well against experiment and the analytical basis sets for considerably less computational expense.     

Geometrical Structures and Electronic Properties of Copper-Doped Aluminum Clusters

Using density function theory (DFT), the Cu-doped Aln (n=1?15) clusters have been stud-ied. The electron a±nity, ionization potential, Mulliken population analysis of Cu, mean polarizability, polarizability anisotropy, dipole moments and HOMO-LUMO gaps have also been calculated on the basis of optimized geometries. The results indicate that there is magic numbers in copper-doped aluminum clusters and electronic characteristic depended on the size of clusters. As n=13, the electron affinity and ionization potential of cluster changed more than 0.3 and 0.6 eV respectively, compared with neighborhood clusters.     

Electron affinities of 1,1-diaryl-2,3,4,5-tetraphenylsiloles: direct measurements and comparison with experimental and theoretical estimates

We present a comprehensive experimental and theoretical characterization of the electronic structure of four 1,1-diaryl-2,3,4,5-tetraphenylsiloles (aryl = phenyl, 2-(9,9-dimethylfluorenyl), 2-thienyl, pentafluorophenyl). Solid-state electron affinities and ionization potentials of these siloles were measured using inverse-photoelectron spectroscopy (IPES) and photoelectron spectroscopy (PES), respectively; the density of electronic states obtained from calculations performed at the density functional theory (DFT) level corresponds very well to the PES and IPES data. The direct IPES measurements of electron affinity were then used to assess alternative estimates based on electrochemical and/or optical data. We also used DFT to calculate the reorganization energies for the electron-transfer reactions between these siloles and their radical anions. Additionally, optical data and ionization potential and electron affinity data were utilized to estimate the binding energies of excitons in these siloles.     

Joint photoelectron and theoretical study of (Rh(m)Co(n))⁻ (m=1-5, n=1-2) cluster anions and their neutral counterparts

Anion photoelectron spectroscopic experiments and calculations based on density functional theory have been used to investigate and uniquely identify the structural, electronic, and magnetic properties of both neutral and anionic (Rh(m)Co(n)) and (Rh(m)Co(n))(-) (m=1-5, n=1-2) clusters, respectively. Negative ion photoelectron spectra are presented for electron binding energies up to 3.493 eV. The calculated electron affinities and vertical detachment energies are in good agreement with the measured values. Computational results for geometric structures and magnetic moments of both cluster anions and their neutrals are presented.     

Small ScCn cyclic clusters: a density functional study of their structure and stability

A theoretical study of the ScCn, ScCn+, and ScCn- (n = 1-10) cyclic clusters has been carried out employing the B3LYP density functional method. Predictions for several molecular properties that could help in their possible experimental characterization, such as equilibrium geometries, electronic structures, dipole moments, and vibrational frequencies, are reported. All ScCn cyclic clusters are predicted to have doublet ground states. For cationic clusters the ground state is alternate between singlets (n-even species) and triplets (n-odd members). In the case of anionic clusters the singlet-triplet separation is relatively small, with the singlets being favored in most cases. In general, even-odd parity effects are also observed for different properties, such as incremental binding energies, ionization energies, and electron affinities. For all neutral, cationic, and anionic clusters it is found that cyclic species are more stable than their open-chain counterparts. Therefore, cyclic structures are the most interesting possible targets for an experimental search of scandium-doped carbon clusters.     

A study of the reactions of molecular hydrogen with small gold clusters

This work presents a study of reactions between neutral and negatively charged Au(n) clusters (n=2,3) and molecular hydrogen. The binding energies of the first and second hydrogen molecule to the gold clusters were determined using density functional theory (DFT), second order perturbation theory (MP2) and coupled cluster (CCSD(T)) methods. It is found that molecular hydrogen easily binds to neutral Au(2) and Au(3) clusters with binding energies of 0.55 eV and 0.71 eV, respectively. The barriers to H(2) dissociation on these clusters with respect to Au(n)H(2) complexes are 1.10 eV and 0.59 eV for n=2 and 3. Although negatively charged Au(n) (-) clusters do not bind molecular hydrogen, H(2) dissociation can occur with energy barriers of 0.93 eV for Au(2) (-) and 1.39 eV for Au(3) (-). The energies of the Au(2)H(2) (-) and Au(3)H(2) (-) complexes with dissociated hydrogen molecules are lower than the energies of Au(2) (-)+H(2) and Au(3) (-)+H(2) by 0.49 eV and 0.96 eV, respectively. There is satisfactory agreement between the DFT and CCSD(T) results for binding energies, but the agreement is not as good for barrier heights.     

Electronic and structural evolution and chemical bonding in ditungsten oxide clusters: W2O(n)- and W2O(n) (n = 1-6)

We report a systematic and comprehensive investigation of the electronic structures and chemical bonding in a series of ditungsten oxide clusters, W2O(n)- and W2O(n) (n = 1-6), using anion photoelectron spectroscopy and density functional theory (DFT) calculations. Well-resolved photoelectron spectra were obtained at several photon energies (2.331, 3.496, 4.661, 6.424, and 7.866 eV), and W 5d-based spectral features were clearly observed and distinguished from O 2p-based features. More complicated spectral features were observed for the oxygen-deficient clusters because of the W 5d electrons. With increasing oxygen content in W2O(n)-, the photoelectron spectra were observed to shift gradually to higher binding energies, accompanied by a decreasing number of W 5d-derived features. A behavior of sequential oxidation as a result of charge transfers from W to O was clearly observed. A large energy gap (2.8 eV) was observed in the spectrum of W2O6-, indicating the high electronic stability of the stoichiometric W2O6 molecule. Extensive DFT calculations were carried out to search for the most stable structures of both the anion and neutral clusters. Time-dependent DFT method was used to compute the vertical detachment energies and compare to the experimental data. Molecular orbitals were used to analyze the chemical bonding in the ditungsten oxide clusters and to elucidate their electronic and structural evolution.     

Electron solvation in water-ammonia mixed clusters: Structure, energetics, and the nature of localization states of the excess electron

The structure and energetics of water-ammonia mixed clusters with an excess electron, [(H2O)n(NH3)m]- with m=1, n=2-6 and m=2, n=2, and also the corresponding neutral clusters are investigated in detail by means of ab initio quantum chemical calculations. The authors focus on the localization structure of the excess electron with respect to its surface versus interiorlike states, its binding to ammonia versus water molecules, the spatial and orientational arrangement of solvent molecules around the excess electron, the changes of the overall hydrogen-bonded structure of the clusters as compared to those of the neutral ones and associated dipole moment changes, vertical detachment energies of the anionic clusters, and also the vertical attachment energies of the neutral clusters. It is found that the hydrogen-bonded structure of the anionic clusters are very different from those of the neutral clusters unlike the case of water-ammonia dimer anion, and these changes in structural arrangements lead to drastically different dipole moments of the anionic and the neutral clusters. The spatial distribution of the singly occupied molecular orbital holding the excess electron shows only surface states for the smaller clusters. However, for n=5 and 6, both surface and interiorlike binding states are found to exist for the excess electron. For the surface states, the excess electron can be bound to the dangling hydrogens of either an ammonia or a water molecule with different degrees of stability and vertical detachment energies. The interiorlike states, wherever they exist, are found to have a higher vertical detachment energy than any of the surface states of the same cluster. Also, for interiorlike states, the ammonia molecule with its dangling hydrogens is always found to stay on top or on a far side of the charge density of the excess electron without participating in the hydrogen bond network of the cluster; the intermolecular hydrogen bonds are formed by the water molecules only which add to the overall stability of these anionic clusters.     

Electronic Structure and Optical Properties of the Surface F-Centers in MgO: A Theoretical Analysis by DFT Approach

Electronic structure of Mg9O9 and Mg9O8 clusters modeling nano-crystalline powders of magnesium oxide has been analyzed within the frames of the density functional theory (DFT). In the framework of time-dependent DFT method (TD-DFT), the relationship between the surface and bulk properties of nano-crystals is analyzed based on variations in the density of electronic states (DOS) and changes of electronic spectra. The spectroscopy of spatial defects like low-coordinated oxygen ions and of surface point defects like F⁺- and F-centers is investigated. Optical properties of the nano-sized crystalline magnesium oxide are characterized by a spectrum of absorption bands in the range of 1-5 eV. Point defects such as F-centers absorb light in the range of 1.2-1.5 eV. Spatial defects O_LC in nano-crystals generate absorption bands in the range of 2.5-5.0 eV. According to calculations, there is no direct relation between coordination numbers of surface ions and excitation energies. Theoretical excitation energies are compared with experimental optical properties of the F⁺- and F-centers.     

Vacuum ultraviolet mass-analyzed threshold ionization spectroscopy of hexafluorobenzene: the Jahn-Teller effect and vibrational analysis

One-photon mass-analyzed threshold ionization (MATI) spectrum of hexafluorobenzene was obtained by using vacuum ultraviolet radiation generated by four-wave difference frequency mixing in Kr. The ionization energy of hexafluorobenzene determined from the position of the 0-0 band was 9.9108+/-0.0006 eV. To aid the spectral analysis, the Jahn-Teller coupling parameters for four e(2g) modes of C(6)F(6) (+) in the ground electronic state were calculated from the topographical data of the potential energy surface obtained at the density functional theory (DFT) level. These were used in the initial calculation of the energies of the Jahn-Teller states and upgraded through the multimode fit to the experimental data. Excellent agreement between the experimental and calculated frequencies was achieved. The vibrations which are not linear Jahn-Teller active were observed and could be assigned by referring to the frequencies obtained at the DFT level.     

Electronic and optical properties of edge modified peritetracene: a DFT study

Peritetracene (PTA) molecules have promising applications in organic electronics and organic light-emitting diodes, but the major constraints come from their poor stability with higher energy gap. We have investigated the stability, electronic, and optical properties of different electron-donating- and electron-withdrawing-substituted PTA molecule groups using density functional theory (DFT) and time-dependent density functional theory (TDDFT) methods. On substituting suitable electron donor and acceptor groups, the energy gap of PTA molecules reduces up to 0.37 eV with an increase in their stability. The stabilities of considered PTA molecules have been investigated using the distribution pattern of frontier molecular orbital energies. The charge transfer properties with smaller ionization potential and larger electron affinity for PTA molecule have been evaluated using Koopmans’ theorem. Enhancement in dipole moment, polarizability, and hyperpolarizability of PTA substituted with electron-donating groups shows the presence of NLO properties. Also, we have investigated the spectroscopic properties of different electron-donating- and electron-withdrawing-substituted PTA molecule groups. Our spectroscopic properties show the bathochromic shift and small hypsochromic shifts in wavelengths of PTA with substituents groups. It is therefore concluded that the –NCH₃ and –NO₂ substituent groups on PTA are observed to have the strongest and highest stability than the other substituent groups considered.

     

Theoretical study of the adsorption and dissociation of oxygen on Pt(111) in the presence of homogeneous electric fields

The effect of homogeneous electric fields on the adsorption energies of atomic and molecular oxygen and the dissociation activation energy of molecular oxygen on Pt(111) were studied by density functional theory (DFT). Positive electric fields, corresponding to positively charged surfaces, reduce the adsorption energies of the oxygen species on Pt(111), whereas negative fields increase the adsorption energies. The magnitude of the energy change for a given field is primarily determined by the static surface dipole moment induced by adsorption. On 10-atom Pt(111) clusters, the adsorption energy of atomic oxygen decreased by ca. 0.25 eV in the presence of a 0.51 V/A (0.01 au) electric field. This energy change, however, is heavily dependent on the number of atoms in the Pt(111) cluster, as the static dipole moment decreases with cluster size. Similar calculations with periodic slab models revealed a change in energy smaller by roughly an order of magnitude relative to the 10-atom cluster results. Calculations with adsorbed molecular oxygen and its transition state for dissociation showed similar behavior. Additionally, substrate relaxation in periodic slab models lowers the static dipole moment and, therefore, the effect of electric field on binding energy. The results presented in this paper indicate that the electrostatic effect of electric fields at fuel cell cathodes may be sufficiently large to influence the oxygen reduction reaction kinetics by increasing the activation energy for dissociation.     

Negative ions of nitroethane and its clusters

Valence and dipole-bound negative ions of the nitroethane (NE) molecule and its clusters are studied using photoelectron spectroscopy (PES), Rydberg electron transfer (RET) techniques, and ab initio methods. Valence adiabatic electron affinities (EA(a)s) of NE, C(2)H(5)NO(2), and its clusters, (C(2)H(5)NO(2))(n), n=2-5, are estimated using vibrationally unresolved PES to be 0.3+/-0.2 eV (n=1), 0.9+/-0.2 eV (n=2), 1.5+/-0.2 eV (n=3), 1.9+/-0.2 eV (n=4), and 2.1+/-0.2 eV (n=5). These energies were then used to determine stepwise anion-neutral solvation energies and compared with previous literature values. Vertical detachment energies for (C(2)H(5)NO(2))(n)(-) were also measured to be 0.92+/-0.10 eV (n=1), 1.63+/-0.10 eV (n=2), 2.04+/-0.10 eV (n=3), and 2.3+/-0.1 eV (n=4). RET experiments show that Rydberg electrons can be attached to NE both as dipole-bound and valence bound anion states. The results are similar to those found for nitromethane (NM), where it was argued that the diffuse dipole state act as a "doorway state" to the more tightly bound valence anion. Using previous models for relating the maximum in the RET dependence of the Rydberg effective principle number n(max)(*), the dipole-bound electron affinity is predicted to be approximately 25 meV. However, a close examination of the RET cross section data for NE and a re-examination of such data for NM finds a much broader dependence on n(*) than is seen for RET in conventional dipole bound states and, more importantly, a pronounced [l] dependence is found in n(max)(*) (n(max)(*) increases with [l]). Ab initio calculations agree well with the experimental results apart from the vertical electron affinity value associated with the dipole bound state which is predicted to be 8 meV. Moreover, the calculations help to visualize the dramatic difference in the distributions of the excess electron for dipole-bound and valence states, and suggest that NE clusters form only anions where the excess electron localizes on a single monomer.     

Structural and electronic properties of reduced transition metal oxide clusters, M4O10 and M4O10- (M = Cr, W), from photoelectron spectroscopy and quantum chemical calculations

Anion photoelectron spectroscopy and quantum chemical calculations at the density functional theory (DFT), coupled cluster theory (CCSD(T)), and complete active space self-consistent field (CASSCF) theory levels are employed to study the reduced transition metal oxide clusters M(4)O(10)(-) (M = Cr, W) and their neutrals. Photoelectron spectra are obtained at 193 and 157 nm photon energies, revealing very different electronic structures for the Cr versus W oxide clusters. The electron affinity and HOMO-LUMO gap are measured to be 3.68 ± 0.05 and 0.7 eV, respectively, for the Cr(4)O(10) neutral cluster, as compared to 4.41 ± 0.04 and 1.3 eV for W(4)O(10). A comprehensive search is performed to determine the ground-state structures for M(4)O(10) and M(4)O(10)(-), in terms of geometry and electronic states by carefully examining the calculated relative energies at the DFT, CCSD(T), and CASSCF levels. The ground states of Cr(4)O(10) and Cr(4)O(10)(-) have tetrahedral structures similar to that of P(4)O(10) with the anion having a lower symmetry due to a Jahn-Teller distortion. The ground states of W(4)O(10) and W(4)O(10)(-) have butterfly shape structures, featuring two fused five-member rings with a metal-metal multiple bond between the central metal atoms. The much stronger WW bonding than the CrCr bonding is found to be the primary cause for the different ground state structures of the reduced Cr(4)O(10)(0/-) versus W(4)O(10)(0/-) oxide clusters. The photoelectron spectra are assigned by comparing the experimental and theoretical adiabatic and vertical electron detachment energies, further confirming the determination of the ground electronic states of M(4)O(10) and M(4)O(10)(-). The time-dependent DFT method is used to calculate the excitation energies of M(4)O(10). The TD-DFT results in combination with the self-consistently calculated vertical detachment energies for some of the excited states at the DFT and CCSD(T) levels are used to assign the higher energy bands. Accurate clustering energies and heats of formation of M(4)O(10) are calculated and used to calculate accurate reaction energies for the reduction of M(4)O(12) to M(4)O(10) by CH(3)OH, as well as for the oxidation of M(4)O(10) to M(4)O(12) by O(2). The performance of the DFT method with the B3LYP and BP86 functionals in the calculations of the relative energies, electron detachment energies, and excitation energies are evaluated, and the BP86 functional is found to give superior results for most of these energetic properties.     

Photoinduced charge-transfer dehydrogenation in a gas-phase metal-DNA base complex: Al-cytosine

An Al-cytosine association complex has been generated via laser ablation of a mixture of aluminum and cytosine powders that were pressed into a rod form. The ionization energy of the complex is found to be 5.16 +/- 0.01 eV. The photoionization efficiency spectrum of Al-cytosine has also been collected. DFT calculations indicate that binding of Al to cytosine manifests a significant weakening of the N-H bond, predicted to have a strength of 1.5 eV in the complex, and a significant stabilization of the oxo tautomeric form relative to the hydroxy forms. The predicted ionization energy of 5.2 eV agrees well with the experimental value. The threshold for dehydrogenation/ionization of Al-cytosine, forming (Al-cytosine-H)+, is found to occur at photoexcitation energies between 11.4 and 12.8 eV. This is a two-photon process that is proposed to occur via photoinduced electron transfer from Al to an antibonding (sigma) orbital localized on N-H. In the context of this mechanism, this work constitutes the first time charge transfer between a metal and DNA base has been photoinitiated in the gas phase.     

Norbornane: an investigation into its valence electronic structure using electron momentum spectroscopy, and density functional and Green's function theories

We report on the results of an exhaustive study of the valence electronic structure of norbornane (C(7)H(12)), up to binding energies of 29 eV. Experimental electron momentum spectroscopy and theoretical Green's function and density functional theory approaches were all utilized in this investigation. A stringent comparison between the electron momentum spectroscopy and theoretical orbital momentum distributions found that, among all the tested models, the combination of the Becke-Perdew functional and a polarized valence basis set of triple-zeta quality provides the best representation of the electron momentum distributions for all of the 20 valence orbitals of norbornane. This experimentally validated quantum chemistry model was then used to extract some chemically important properties of norbornane. When these calculated properties are compared to corresponding results from other independent measurements, generally good agreement is found. Green's function calculations with the aid of the third-order algebraic diagrammatic construction scheme indicate that the orbital picture of ionization breaks down at binding energies larger than 22.5 eV. Despite this complication, they enable insights within 0.2 eV accuracy into the available ultraviolet photoemission and newly presented (e,2e) ionization spectra, except for the band associated with the 1a(2) (-1) one-hole state, which is probably subject to rather significant vibronic coupling effects, and a band at approximately 25 eV characterized by a momentum distribution of "s-type" symmetry, which Green's function calculations fail to reproduce. We note the vicinity of the vertical double ionization threshold at approximately 26 eV.     

Electronic properties of liquid water by sequential molecular dynamics/density functional theory

Electronic properties of liquid water were analysed by a sequential molecular dynamics (MD)/density functional theory approach. MD simulations are based on a polarisable model for water. Emphasis was placed on the prediction of the water dipole moment, liquid state polarisability, ionisation potential (IP), and vertical electron affinity. The dipole moment of the water molecule in liquid water is not dependent on the number of molecules included in the quantum mechanical calculations. The polarisability of the water molecule in liquid water is 4% lower than its gas phase value. The IP of liquid water (9.7 ± 0.06 eV) is in good agreement with recent experimental data.     

用密度泛函和遗传算法研究Cun (n≤20) 团簇的尺寸效应

结合遗传算法和Gupta多体势系统地搜索金属团簇Cu_n (n≤20)的几何结构, 并利用密度泛函方法进一步确定最稳定构型. 分析了平均键长、平均配位数、结合能、二阶差分能、电离势和电子亲和势等性质随着尺寸的变化规律. 发现在Cu₇处团簇最稳定构型从二维结构转向三维结构, Cu_n (n≤20)团簇的幻数为8, 13, 20. 团簇的键长、配位数和结合能属性随着尺寸的增长而递增最终接近相应的体相值; 而二阶差分能、电离势和电子亲和势随着尺寸增加出现奇偶交替, 说明偶数电子形成闭壳层结构, 比相邻团簇更稳定.