Electronic structure of spheroidal metal containing carbon shells: study of the LaC60 and C60 clusters and their ions within the local density approximation

Ionization thresholds and electron affinities have been calculated within the Local Density Approximation (LDA) for the neutral, positively and negatively charged clusters of LaC₆₀ and C₆₀. The evaluated energies are found to be in good agreement with available experimental data. More accurate measurements are however necessary to verify the suggested spheroidal cage structure for these molecules.     

Accurate energies calculated by empirical corrections to the local spin density approximation

We tested the efficacy of three empirical correction schemes on atomization energies calculated by the following theories: Kohn–Sham density functional theory (KS-DFT) with local spin density approximation (LSDA), two KS-DFT gradient-corrected methods, one hybrid Hartree–Fock/KS-DFT method similar to B3LYP, and the ab initio extrapolation procedures G1 and G2. Empirical corrections improved the LSDA results greatly, while the other theories were improved slightly or not at all. The best procedure for correcting LSDA atomization energies brings the mean absolute deviation from experiment from 38.3 to 4.0 kcal/mol on a subset of 44 molecules in the G1 dataset that were not used in deriving the empirical parameters. This corrected LSDA is interesting for three reasons: it could be a useful computational tool in some cases, it implies that the LSDA itself gives accurate energies for reactions where atomic coordinations stay unchanged, and it gives insight into the search of better functionals.     

The self-interaction-corrected local spin density approximation. a case of too much correlation

We have carried out total-energy calculations for atoms from H to Ar using the exchange-only (XO) version of Perdew? orbital self-interaction-correction (SIC) functional. We find that for atoms beyond oxygen, this functional yields better total energies than the full exchange-correlation (XC) version.     

Theory for spin relaxation in small magnetic metal clusters

We discuss the magnetic properties of small neutral transition-metal clusters like Fe _n and Co _n deduced from Stern-Gerlach deflection experiments. We claim that the asymmetric Stern-Gerlach deflection profiles are due to a transfer from electronical angular momentum to the cluster rotation, allowing for a depopulation of the high energy magnetic levels. For finite temperatures we consider two limiting cases. First, the cluster magnetization is assumed to be tied to the random orientation of the cluster easy axes due to the lattice anisotropy. This causes a surprisingly small magnetization for small external magnetic fields. For larger fields and also for increasing temperatures the magnetization is released from the cluster geometry and allowed to align itself parallel to the field. In the second case the clusters are treated as an ensemble of superparamagnetic particles. Here, the effect of the anisotropy is less visible. The cluster lattice anisotropy per atom is expected to decrease for increasing cluster size. Preliminary results support this.     

Some efforts beyond the local density approximation

The exchange-correlation density functional can be expressed as a many-body perturbation series in terms of the Coulomb interaction using the exact Kohn-Sham orbitals as the basis. A self-consistent equation is derived for the exact exchangecorrelation potential. This perturbation approach forms a basis for going beyond the local density approximation (LDA ). The discontinuity in the exchange-correlation potential for semiconductors calculated by the perturbative approach gives a good account of the discrepancy of the band gap calculated in LDA . The discontinuity also plays an important role in the interface band diagrams. A theory to account for the interaction effects of localized d or f orbitals is reviewed and the physics of the applications to a model test, to some 3d transition metals, and to heavy fermions is discussed. The perturbative approach to improvement beyond LDA tends to be computation-intensive and to be system-specific. © 1995 John Wiley & Sons, Inc.     

BFW: a density functional for transition metal clusters

Ionization potentials (IPs) or electron affinities (EAs) for transition metal clusters are an important property that can be used to identify and differentiate between clusters. Accurate calculation of these values is therefore vital. Previous attempts using a variety of DFT models have correctly predicted trends, but have relied on the use of scaling factors to compare to experimental IPs. In this paper, we introduce a new density functional (BFW) that is explicitly designed to yield accurate, absolute IPs for transition metal clusters. This paper presents the numerical results for a selection of transition metal clusters and their carbides, nitrides, and oxides for which experimental IPs are known. When tested on transition metal clusters, the BFW functional is found to be significantly more accurate than B3LYP and B3PW91.     

Generalized gradient approximation adjusted to transition metals properties: Key roles of exchange and local spin density

Perdew-Burke-Ernzerhof (PBE) and PBE adapted for solids (PBEsol) are exchange-correlation (xc) functionals widely used in density functional theory simulations. Their differences are the exchange, μ, and correlation, β, coefficients, causing PBEsol to lose the Local Spin Density (LSD) response. Here, the μ/β two-dimensional (2D) accuracy landscape is analyzed between PBE and PBEsol xc functional limits for 27 transition metal (TM) bulks, as well as for 81 TM surfaces. Several properties are analyzed, including the shortest interatomic distances, cohesive energies, and bulk moduli for TM bulks, and surface relaxation degree, surface energies, and work functions for TM surfaces. The exploration, comparing the accuracy degree with respect experimental values, reveals that the found xc minimum, called VV, being a PBE variant, represents an improvement of 5% in mean absolute percentage error terms, whereas this improvement reaches ~11% for VVsol, a xc resulting from the restoration of LSD response in PBEsol, and so regarded as its variant.     

First-principles local density approximation (generalized gradient approximation) +U study of catalytic CenOm clusters: U value differs from bulk

Ceria possesses strong catalytic properties for CONO(x) removal and H(2) production. Clusters often show more intriguing functionalities than their bulk counterparts. Here, the geometric and electronic structures of Ce(n)O(m) (n=1-4,m=2n-1,2n) clusters are studied for the first time using the projected augmented wave method in density functional theory with detailed assessment of the exchange-correlation functional and the Hubbard parameter U. We note that the U value strongly affects the electronic structures of the oxygen-deficient Ce(n)O(2n-1) clusters, though less so on the stoichiometric Ce(n)O(2n). Furthermore, the local density approximation (LDA)+U method is more accurate than the generalized gradient approximation+U in describing the localization of the 4f electrons of the Ce(n)O(m) clusters. The calculated vibration frequency of the CeO molecule with the LDA+U (U=4 eV) is 818.4 cm(-1), in close agreement with experimental values of 820-825 cm(-1) for the low lying states. Different optimal U values were noted for the ceria cluster (4 eV) and its bulk (6 eV), due to quantum-size and geometric effects. The largely reduced formation energy of an oxygen vacancy indicates that the catalytic effect of the Ce(n)O(m) clusters are far greater than bulk CeO(2).     

Time‐dependent local‐density approximation in real time: Application to conjugated molecules

The time‐dependent local‐density approximation (TDLDA) is applied to the optical response of conjugated carbon molecules in the energy range of 0–30 eV, with calculations given for carbon chains, polyenes, retinal, benzene, and C₆₀. The major feature of the spectra, the collective π–π* transition, is seen at energies ranging from below 2 to 7 eV and is reproduced by the theory to a few tenths of an electron volt with a good account of systematic trends. However, there is some indication that TDLDA predicts too much fragmentation of the strength function in large molecules. Transition strengths are reproduced with a typical accuracy of 20%. The theory also predicts a broad absorption peak in the range of 15–25 eV, and this feature agrees with experiment in the one case where quantitative data is available (benzene). ©1999 John Wiley & Sons, Inc. Int J Quant Chem 75: 55–66, 1999     

Shell structure and response properties of metal clusters

Electronic shell structure, which was first recognized in sodium clusters, has been observed in alkali and noble metals, as well as in divalent and trivalent metals. Shell structure with modifications is expected to be broadly applicable to most metals. Features in the cluster abundance spectra and in the experimental dipole polarizabilities and ionization potentials correlate well with predictions of electronic level filling in spherical and spheroidal potential wells. The lack of precise quantitative agreement between experiment and theory for the response properties indicates necessary refinements in the self-consistent uniform background jellium model for clusters.     

Nonlocal density functional calculations of the jellium metal surface

The Becke exchange functional is used for calculation of properties of the jellium model using the slab geometry inside a box with the infinite potential barriers at the boundaries. We simplify semianalytical representation of matrix elements. We calculate the surface energies and work functions with self-consistent electron densities. For all densities (here, we give results in erg/cm² for r_s = 2.07 bohr, in comparison with the LSD approximation (?602)) and the uncorrected Pw GGA -II (?730), the Becke-II exchange only (?1212), and the Becke-II exchange with Perdew86 correlation (?830) [always close to Pw GGA -I (?814)] give smaller surface energies. The most important factor determining values of surface energies from different GGAS seems to be a form of a correlation potential. We also calculate the effect of finite slab thickness and the vacuum region thickness on the surface energy at the LSD level and indicate its importance in various jellium model calculations. © 1995 John Wiley & Sons, Inc.     

Synthetic strategies for the surface functionalisation of gold nanoparticles with metals and metal clusters

The reaction of the new ditopic thiol-phosphine compound HS(CH(2))(11)OOCC(6)H(4)PPh(2) (L) with an excess of dodecanethiol-protected gold nanoparticles gave the asymmetric gold complex [CH(3)(CH(2))(11)SAuPPh(2)C(6)H(4)COO(CH(2))(11)SH] (4), but no phosphine-protected gold nanoparticles were formed. However, by blocking the phosphine function in L with metal fragments, we have been able to produce gold nanoparticles functionalised with AuCl- and cluster [Fe(2)(CO)(7)Au] units on the surface by the method of ligand-place exchange reaction.     

The response of metal clusters toq- andL-dependent external fields

We have calculated the static polarizability and mean excitation energy of metal clusters submitted toq-andL-dependent external fields ofj _L (qr)Y _L0(Ω) type. Use has been made of an Extended Random-Phase Approximation which includes exchange and correlation effects within a local model, and of the spherical jellium model to describe the neutralizing positive background.     

Adiabatic approximation of time-dependent density matrix functional response theory

Time-dependent density matrix functional theory can be formulated in terms of coupled-perturbed response equations, in which a coupling matrix K(omega) features, analogous to the well-known time-dependent density functional theory (TDDFT) case. An adiabatic approximation is needed to solve these equations, but the adiabatic approximation is much more critical since there is not a good "zero order" as in TDDFT, in which the virtual-occupied Kohn-Sham orbital energy differences serve this purpose. We discuss a simple approximation proposed earlier which uses only results from static calculations, called the static approximation (SA), and show that it is deficient, since it leads to zero response of the natural orbital occupation numbers. This leads to wrong behavior in the omega-->0 limit. An improved adiabatic approximation (AA) is formulated. The two-electron system affords a derivation of exact coupled-perturbed equations for the density matrix response, permitting analytical comparison of the adiabatic approximation with the exact equations. For the two-electron system also, the exact density matrix functional (2-matrix in terms of 1-matrix) is known, enabling testing of the static and adiabatic approximations unobscured by approximations in the functional. The two-electron HeH(+) molecule shows that at the equilibrium distance, SA consistently underestimates the frequency-dependent polarizability alpha(omega), the adiabatic TDDFT overestimates alpha(omega), while AA improves upon SA and, indeed, AA produces the correct alpha(0). For stretched HeH(+), adiabatic density matrix functional theory corrects the too low first excitation energy and overpolarization of adiabatic TDDFT methods and exhibits excellent agreement with high-quality CCSD ("exact") results over a large omega range.     

The density of surface ligands regulates the luminescence of thiolated gold nanoclusters and their metal ion response

The fascinating luminescence properties of gold nanoclusters(AuNCs) have drawn considerable research interests,and been widely harnessed for a wide range of applications.However,a fundamental understanding towards ligand density's role in the luminescence properties of these ultrasmall AuNCs remains unclear yet.In this communication,through systematic investigation of surface chemistries of glutathione-protected Au NCs(GSH-Au NCs) with diffe rent density of GSH as well as other thiolates,it is discovered that the density of surface ligands can significantly regulate the luminescence properties of AuNCs.Fluorescence lifetime spectroscopy and X-ray photoelectron spectroscopy showed that AuNCs with a higher density of electron-rich ligands facilitate their luminescence generation.Moreover,differences in the surface coverage of AuNCs can also affect their interactions with foreign species,as illustrated by significantly different fluorescence quenching capability of GSH-AuNCs with different ligand density towards Hg~(2+).This study provides new insight into the intriguing luminescence properties of metal NCs,which is hoped to stimulate further research on the design of metal NCs with strong luminescence and sensitive/specific responses for promising optoelectronic,sensing and imaging applications.     

Description of core excitations by time-dependent density functional theory with local density approximation, generalized gradient approximation, meta-generalized gradient approximation, and hybrid functionals

Time-dependent density functional theory (TDDFT) is employed to investigate exchange-correlation-functional dependence of the vertical core-excitation energies of several molecules including H, C, N, O, and F atoms. For the local density approximation (LDA), generalized gradient approximation (GGA), and meta-GGA, the calculated X1s-->pi* excitation energies (X = C, N, O, and F) are severely underestimated by more than 13 eV. On the other hand, time-dependent Hartree-Fock (TDHF) overestimates the excitation energies by more than 6 eV. The hybrid functionals perform better than pure TDDFT because HF exchange remedies the underestimation of pure TDDFT. Among these hybrid functionals, the Becke-Half-and-Half-Lee-Yang-Parr (BHHLYP) functional including 50% HF exchange provides the smallest error for core excitations. We have also discovered the systematic trend that the deviations of TDHF and TDDFT with the LDA, GGA, and meta-GGA functionals show a strong atom-dependence. Namely, their deviations become larger for heavier atoms, while the hybrid functionals are significantly less atom-dependent.     

First-principles local density approximation + U and generalized gradient approximation + U study of plutonium oxides

The electronic structure and properties of PuO2 and Pu2O3 have been studied from first principles by the all-electron projector-augmented-wave method. The local density approximation+U and the generalized gradient approximation+U formalisms have been used to account for the strong on-site Coulomb repulsion among the localized Pu 5f electrons. We discuss how the properties of PuO2 and Pu2O3 are affected by the choice of U as well as the choice of exchange-correlation potential. Also, oxidation reaction of Pu2O3, leading to formation of PuO2, and its dependence on U and exchange-correlation potential have been studied. Our results show that by choosing an appropriate U, it is promising to correctly and consistently describe structural, electronic, and thermodynamic properties of PuO2 and Pu2O3, which enable the modeling of redox process involving Pu-based materials possible.