Search

In the ground state, the pair density n can be determined by solving a single auxiliary equation of a two-particle problem. Electron-electron cusp condition and asymptotic behavior for the Pauli potential of the effective potential of the two-particle equation are presented. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

We propose a direct method for reducing the dimension of the space of orbital products that occur, for example, in the calculation of time dependent density functional theory linear response and in Hedin’s GW approximation to the electron propagator. We do this by defining, within the linear space of orbital products, a subspace of dominant directions that are associated with a certain eigenvalue problem. These directions span the entire linear space of products with an error that decreases approximately exponentially with their number. Our procedure works best for atomic orbitals of finite range and it avoids the use of extra sets of auxiliary fit functions. [ABSTRACT FROM AUTHOR]

Stable and metastable fluid density distributions (FDDs) in a closed nanoslit between two identical parallel solid walls have been identified on the basis of a nonlocal canonical ensemble density functional theory. Similar to Monte Carlo simulations, periodicity of the FDD in one of the lateral (parallel to the walls surfaces) directions, denoted as the x direction, was assumed. In the other lateral direction, y direction, the FDD was considered uniform. It was found that depending on the average fluid density in the slit, both uniform as well as nonuniform FDDs in the x direction can occur. The uniform FDDs are either symmetric or asymmetric about the middle plane between walls; the latter FDD being the consequence of a symmetry breaking across the slit. The nonuniform FDDs in the x direction occur either in the form of a bump on a thin liquid film covering the walls or as a liquid bridge between those walls and provide symmetry breaking in the x direction. For small and large average densities, the stable state is uniform in the x direction and is symmetric about the middle plane between walls. In the intermediate range of the average density and depending on the length Lx of the FDD period, the stable state can be represented either by a FDD, which is uniform in the x direction and asymmetric about the middle of the slit (small values of Lx), or by a bump- and bridgelike FDD for intermediate and large values of Lx, respectively. These results are in agreement with the Monte Carlo simulations performed earlier by other authors. Because the free energy of the stable state decreases monotonically with increasing Lx, one can conclude that the real period is very large (infinite) and that for the values of the parameters employed, a single bridge of finite length over the entire slit is generated. [ABSTRACT FROM AUTHOR]

Soft x-ray natural circular dichroism of amino acids is studied by means of ab initio methods. Several approaches to evaluate the oscillator and rotary strengths of core-to-valence excitations are compared from the viewpoint of basis set dependence: ground-state Hartree-Fock (HF) orbital set employed in (i) random phase approximation (RPA), (ii) static exchange approach (STEX) (unrelaxed), (iii) core-ionized state HF orbital set applied in STEX (relaxed), and (iv) HF excited state orbital set for each core-to-valence excited state. Furthermore in (i) the PRA in the framework of the density functional method (DFT) is compared with the RPA where the ab initio HF orbital set is used. In (iv), the oscillator and rotary strengths evaluated by different orbital sets for the initial and final states, namely, nonorthogonal ground-state and core-excited HF orbitals, are compared with those evaluated by using the core-excited HF orbital set to describe the initial (ground) state. It was shown that, among considered methods, the RPA provides most consistent and less time-consuming results for circular dichroism core excitation spectra. Discussion of the low energy part of K edge circular dichroism spectra of five common amino acids obtained with the help of RPA is presented. [ABSTRACT FROM AUTHOR]

The density profiles in a fluid interacting with the two identical solid walls of a closed long slit were calculated for wide ranges of the number of fluid molecules in the slit and temperature by employing a nonlocal density functional theory. Using argon as the sample fluid and considering the walls composed of solid carbon dioxide, it is shown that the density profile corresponding to the stable state of the fluid considerably changes its shape with increasing average density ρav of the fluid inside the slit. Temperature dependent critical values ρsb1 and ρsb2 of ρav were identified, such that for ρsb1≤=ρav≤=ρsb2 the stable state of the system is described by an asymmetric density profile whereas outside this range it is described by a symmetric one. Hence a spontaneous symmetry breaking of the fluid density distribution in a closed slit with identical walls can take place. On the basis of the results obtained for closed slits, the symmetry breaking in open slits was also examined. [ABSTRACT FROM AUTHOR]

The common density functionals for the exchange-correlation energy make serious self-interaction errors in the molecular dissociation limit when real or spurious noninteger electron numbers N are found on the dissociation products. An “M-electron self-interaction-free” functional for positive integer M is one that produces a realistic linear variation of total energy with N in the range of M-1≤Na=M, and so can avoid these errors. This desideratum is a natural generalization to all M of the more familiar one of one-electron self-interaction freedom. The intent of this paper is not to advocate for any functional, but to understand what is required for a functional to be M-electron self-interaction-free and thus correct even for highly stretched bonds. The original Perdew-Zunger self-interaction correction (SIC) and our scaled-down variant of it are exactly one- and nearly two-electron self-interaction-free, but only the former is nearly so for atoms with M>2. Thus all these SIC’s produce an exact binding energy curve for H2+, and an accurate one for He2+, but only the unscaled Perdew-Zunger SIC produces an accurate one for Ne2+, where there are more than two electrons on each fragment Ne+0.5. We also discuss LiH+, which is relatively free from self-interaction errors. We suggest that the ability of the original and unscaled Perdew-Zunger SIC to be nearly M-electron self-interaction-free for atoms of all M stems in part from its formal resemblance to the Hartree-Fock theory, with which it shares a sum rule on the exchange-correlation hole of an open system. [ABSTRACT FROM AUTHOR]

In previous work the authors have presented a highly accurate two-spinor fully relativistic solution of the two-center Coulomb problem utilizing the finite-element method (FEM) and furthermore developed a relativistic minimax two-spinor linear combination of atomic orbitals (LCAO). In the present paper the authors present Dirac-Fock-Slater (DFS-) density functional calculations for two-atomic molecules up to super heavy systems using the fully nonlinear minimax FEM and the minimax LCAO in its linearized approximation (linear approximation to relativistic minimax). The FEM gives highly accurate benchmark results for the DFS functional. Especially considering molecules with up to super heavy atoms such as UubO and Rg2, the authors found that LCAO fails to give the correct systematic trends. The accurate FEM results shed a new light on the quality of the DFS-density functional. [ABSTRACT FROM AUTHOR]

This paper describes an empirical localized orbital correction model which improves the accuracy of density functional theory (DFT) methods for the prediction of thermochemical properties for molecules of first and second row elements. The B3LYP localized orbital correction version of the model improves B3LYP DFT atomization energy calculations on the G3 data set of 222 molecules from a mean absolute deviation (MAD) from experiment of 4.8 to 0.8 kcal/mol. The almost complete elimination of large outliers and the substantial reduction in MAD yield overall results comparable to the G3 wave-function-based method; furthermore, the new model has zero additional computational cost beyond standard DFT calculations. The following four classes of correction parameters are applied to a molecule based on standard valence bond assignments: corrections to atoms, corrections to individual bonds, corrections for neighboring bonds of a given bond, and radical environmental corrections. Although the model is heuristic and is based on a 22 parameter multiple linear regression to experimental errors, each of the parameters is justified on physical grounds, and each provides insight into the fundamental limitations of DFT, most importantly the failure of current DFT methods to accurately account for nondynamical electron correlation. [ABSTRACT FROM AUTHOR]

A theoretical description of the vibrational excitons in DNA is presented by using the vibrational basis mode theory developed in Papers I and II. The parameters obtained from the density functional theory calculations, such as vibrational coupling constants and basis mode frequencies, are used to numerically simulate two-dimensional (2D) IR spectra of dGn:dCn and dAn:dTn double helices with n varying from 1 to 10. From the molecular dynamics simulations of dG5C5 and dA5T5 double helices in D2O solution, it is found that the thermally driven internal motions of these systems in an aqueous solution do not induce strong fluctuations of basis mode frequencies nor vibrational couplings. In order to construct the two-exciton Hamiltonian, the vibrational anharmonicities of eight basis modes are obtained by carrying out B3LYP/6-31G* calculations for the nine basis modes. The simulated 2D IR spectra of dGn:dCn double helix in D2O solution are directly compared with closely related experimental results. The 2D IR spectra of dGn:dCn and dAn:dTn are found to be weakly dependent on the number of base pairs. The present work demonstrates that the computational procedure combining quantum chemistry calculation and molecular dynamics simulation methods can be of use to predict 2D IR spectra of nucleic acids in solutions. [ABSTRACT FROM AUTHOR]

In this paper we report recontracted LANL2DZ basis sets for first-row transition metals. The valence-electron shell basis functions were recontracted using the PWP86 generalized gradient approximation functional and the hybrid B3LYP one. Starting from the original LANL2DZ basis sets a cyclic method was used in order to optimize variationally the contraction coefficients, while the contraction scheme was held fixed at the original one of the LANL2DZ basis functions. The performance of the recontracted basis sets was analyzed by direct comparison between calculated and experimental excitation and ionization energies. Results reported here compared with those obtained using the original basis sets show clearly an improvement in the reproduction of the corresponding experimental gaps. [ABSTRACT FROM AUTHOR]

In this paper we consider the general form of the correlated-determinantal wave function functional of Colle and Salvetti (CS) for the He atom. The specific form employed by CS is the basis for the widely used CS correlation energy formula and the Lee-Yang-Parr correlation energy density functional of Kohn-Sham density functional theory. We show the following: (i) The key assumption of CS for the determination of this wave function functional, viz., that the resulting single-particle density matrix and the Hartree-Fock theory Dirac density matrix are the same, is equivalent to the satisfaction of the Coulomb hole sum rule for each electron position. The specific wave function functional derived by CS does not satisfy this sum rule for any electron position. (ii) Application of the theorem on the one-to-one correspondence between the Coulomb hole sum rule for each electron position and the constraint of normalization for approximate wave functions then proves that the wave function derived by CS violates charge conservation. (iii) Finally, employing the general form of the CS wave function functional, the exact satisfaction of the Coulomb hole sum rule at each electron position then leads to a wave function that is normalized. The structure of the resulting approximate Coulomb holes is reasonably accurate, reproducing both the short- and the long-range behavior of the hole for this atom. Thus, the satisfaction of the Coulomb hole sum rule by an approximate wave function is a necessary condition for constructing wave functions in which electron-electron repulsion is represented reasonably accurately. [ABSTRACT FROM AUTHOR]

We provide a rigorous proof that the Hohenberg-Kohn theorem holds for spin densities by extending Lieb’s Legendre-transform formulation to spin densities. The resulting spin-density-functional theory resolves several troublesome issues. Most importantly, the present paper provides an explicit construction for the spin potentials at any point along the adiabatic connection curve, thus providing a formal basis for the use of exchange-correlation functionals of the spin density in the Kohn-Sham density-functional theory (DFT). The practical implications of this result for unrestricted Kohn-Sham DFT calculations is considered, and the existence of holes below the Fermi level is discussed. We argue that an orbital’s energy tends to increase as its occupation number increases, which provides the basis for a computational algorithm for determining the occupation numbers in Kohn-Sham DFT and helps explain the origin of Hund’s rules and holes below the Fermi level. [ABSTRACT FROM AUTHOR]

In this paper, we examine the limits of accuracy of the single determinant approximations (Hartree-Fock, optimized effective potential, and density functional theory) to the exact energy eigenstates of many electron systems. We show that an approximate Slater determinant of Sz=M gives maximum accuracy for states with S=M, provided that perturbation theory for the spin up minus spin down potential is applicable. The overlap with the exact energy eigenstates with S≠M is much smaller. Therefore, for the case that the emphasis is on wave functions, one must use symmetry preserving theories, although this is at the expense of accuracy in energy. [ABSTRACT FROM AUTHOR]

Within the scope of studying the molecular implications of the Pb2+ cation in environmental and polluting processes, this paper reports Hartree-Fock and density functional theory (B3LYP) four-component relativistic calculations using an all-electron basis set applied to [Pb(H2O)]2+ and [Pb(OH)]+, two complexes expected to be found in the terrestrial atmosphere. It is shown that full-relativistic calculations validate the use of scalar relativistic approaches within the framework of density functional theory. [Pb(H2O)]2+ is found C2v at any level of calculations whereas [Pb(OH)]+ can be found bent or linear depending of the computational methodology used. When Cs is found the barrier to inversion through the C∞v structure is very low, and can be overcome at high enough temperature, making the molecule floppy. In order to get a better understanding of the bonding occurring between the Pb2+ cation and the H2O and OH- ligands, natural bond orbital and atoms-in-molecule calculations have been performed. These approaches are supplemented by a topological analysis of the electron localization function. Finally, the description of these complexes is refined using constrained-space orbital variation complexation energy decompositions. [ABSTRACT FROM AUTHOR]

This work presents the derivation, implementation, and first applications of the generalized restricted-unrestricted method based on the density functional Kohn-Sham formalism. By using a spin-restricted Kohn-Sham representation for the reference state the well-known spin contamination problem is avoided, while the unrestricted representation of the perturbation response retains a proper description of spin polarization. The formulation is a generalization of our previous implementation of the restricted-unrestricted method [Z. Rinkevicius et al., J. Chem. Phys. 121, 7614 (2004)], as it accounts for the full unrestricted response instead of describing the spin polarization in terms of triplet operators only. The purpose of this paper is to investigate the role of the generalization employed and demonstrate its numerical performance. For this purpose we focus on isotropic hyperfine coupling constants of a set of organic radicals and transition metal compounds. For both classes of molecules we observe that the effect of neglecting singlet excitation operators in the response part of restricted-unrestricted formalism changes the calculated hyperfine coupling parameters by a few MHz. The obtained results confirm the validity of the approximation used in the simplified restricted-unrestricted approach, i.e. spin polarization can, in most cases, be adequately described by an account of only triplet operators in the response term. [ABSTRACT FROM AUTHOR]

We propose a general method of thermodynamic integration to find the free energy of a surface, where our integration parameter is taken to be the strain on the unit cell of the system (which in the example presented in this paper is simply the extension of the unit cell along the normal to the surface), and the integration is performed over the thermal average stress from a molecular dynamics run. In order to open up a vacuum gap in a continuous and reversible manner, an additional control interaction has been introduced. We also use temperature integration to find a linear relation for the temperature dependence of the free surface energy. These methods have been applied to the titanium dioxide (110) surface, using first principles density functional theory. A proof of principle calculation for zero temperature shows excellent agreement between the integral calculation and the difference in energy calculated by the DFT program. Calculations that have been performed at 295 and 1000 K give excellent agreement between the two integration methods. [ABSTRACT FROM AUTHOR]

First, the Thomas-Fermi-Amaldi (TFA) equation was formulated with a newly derived condition to remove the singularities at the nuclei, which coincided with the molecular cusp condition. Next, the collocation method was applied to the TFA equation using the grid-based density functional theory. In this paper, the electron densities and the radial probabilities for specific atoms (He, Be, Ne, Mg, Ar, Ca) were found to agree with those from the Thomas-Fermi-Dirac (TFD) method. Total energies for specific atoms (He, Ne, Ar, Kr, Xe, Rn) and molecules (H2,CH4) were also found to be close to those from the Hartree-Fock method using the Pople basis set 6-311G relative to the TFD method. In addition, the computational expense to determine the electron density and its corresponding energy for a large scale structure, such as a carbon nanotube, is shown to be much more efficient compared to the conventional Hartree-Fock method using the 6-31G Pople basis set. [ABSTRACT FROM AUTHOR]

In this paper a novel approach to study the formation and relaxation of excited states in solution is presented within the integral equation formalism version of the polarizable continuum model. Such an approach uses the excited state relaxed density matrix to correct the time dependent density functional theory excitation energies and it introduces a state-specific solvent response, which can be further generalized within a time dependent formalism. This generalization is based on the use of a complex dielectric permittivity as a function of the frequency, [variant_greek_epsilon]⁁(ω). The approach is here presented in its theoretical formulation and applied to the various steps involved in the formation and relaxation of electronic excited states in solvated molecules. In particular, vertical excitations (and emissions), as well as time dependent Stokes shift and complete relaxation from vertical excited states back to ground state, can be obtained as different applications of the same theory. Numerical results on two molecular systems are reported to better illustrate the features of the model. [ABSTRACT FROM AUTHOR]

In this paper, the second and third order polarizabilities of small GanAsm (n+m=4–10) clusters are systematically investigated using the time dependent density functional theory (TDDFT)/6-311+G* combined with the sum-over-states method (SOS//TDDFT/6-311+G*). For the static second order polarizabilities, the two-level term (βvec.2) makes a significant contribution to the βvec for all considered GanAsm clusters except for the Ga3As4 cluster. And, for the static third order polarizabilities, the positive channel (<γ>II) makes a larger contribution to 〈γ〉tot than the negative channel (<γ>I). Similar to the cubic GaAs bulk materials, the small GanAsm cluster assembled materials exhibit large second order (1×10-6 esu) and third order susceptibilities (5×10-11 esu). The dynamic behavior of β(-2ω;ω,ω) and γ(-3ω;ω,ω,ω) show that the small GanAsm cluster will be a good candidate of nonlinear optical materials due to the avoidance of linear resonance photoabsorption. [ABSTRACT FROM AUTHOR]

In this paper we present the theory and implementation of analytic derivatives of time-dependent density functional theory (TDDFT) excited states energies, both in vacuo and including solvent effects by means of the polarizable continuum model. The method is applied to two case studies: p-nitroaniline and 4-(dimethyl)aminobenzonitrile. For both molecules PCM-TDDFT is shown to be successful in supporting the analysis of experimental data with useful insights for a better understanding of photophysical and photochemical pathways in solution. [ABSTRACT FROM AUTHOR]

Semilocal density functional approximations (DFAs) for the exchange-correlation energy suffer from self-interaction error, which is believed to be the cause of many of the failures of common DFAs, such as poor description of charge transfer and transition states of chemical reactions. The standard self-interaction correction (SIC) of Perdew and Zunger mends some of these failures but spoils such essential properties as thermochemistry and equilibrium bond lengths. The Perdew-Zunger SIC seems to overcorrect many-electron systems. In this paper, we propose a modified SIC, which is scaled down in many-electron regions. The new SIC has an improved performance for many molecular properties, including total energies, atomization energies, barrier heights of chemical reactions, ionization potentials, electron affinities, and bond lengths. The local spin-density approximation (LSDA) benefits from SIC more than higher-level functionals do. The scaled-down SIC has only one adjustable parameter. Rationalization of the optimal value of this parameter enables us to construct an almost-nonempirical version of the scaled-down SIC-LSDA, which is significantly better than uncorrected LSDA and even better than the uncorrected generalized gradient approximation. We present an analysis of the formal properties of the scaled-down SIC and define possible directions for further improvements. In particular, we find that exactness for all one-electron densities does not guarantee correct asymptotics for the exchange-correlation potential of a many-electron system. [ABSTRACT FROM AUTHOR]

Given a general, N-particle Hamiltonian operator, analogs of the Hohenberg-Kohn theorem are derived for functions that are more general than the particle density, including density matrices and the diagonal elements thereof. The generalization of Lieb’s Legendre transform ansatz to the generalized Hohenberg-Kohn functional not only solves the υ-representability problem for these entities, but, more importantly, also solves the N-representability problem. Restricting the range of operators explored by the Legendre transform leads to a lower bound on the true functional. If all the operators of interest are incorporated in the restricted maximization, however, the variational principle dictates that exact results are obtained for the systems of interest. This might have important implications for practical work not only for density matrices but also for density functionals. A follow-up paper will present a useful alternative approach to the v- and N-representability problems based on the constrained search formalism. [ABSTRACT FROM AUTHOR]

In this paper we apply the direct-mapping density-functional theory (DFT) to open-shell systems, in order to get many-electron wave functions having the same transformation properties as the eigenstates of the exact Hamiltonians. Such a case is that of spin, where in order to get the magnetic properties, the many-particle states must be eigenstates not only of Sz but also of S2. In this theory the Kohn and Sham [Phys. Rev. A 140, 1133 (1965)] potential is expressed directly as a mapping of the external potential. The total energies of the molecules calculated were satisfactory as their relative deviations (ΔE/E) from the exact Hartree-Fock ones were of the order of 10-4. This accuracy is much higher than that of the standard DFT in its local exchange potential approximation. This method does not need an approximate density as input, as the effective potential is derived directly from the external potential. [ABSTRACT FROM AUTHOR]

Recently, three of us have proposed a method [Phys. Rev. Lett. 91, 33201 (2003)] for an accurate calculation of the dispersion energy utilizing frequency-dependent density susceptibilities of monomers obtained from time-dependent density-functional theory (DFT). In the present paper, we report numerical calculations for the helium, neon, water, and carbon dioxide dimers and show that for a wide range of intermonomer separations, including the van der Waals and short-range repulsion regions, the method provides dispersion energies with accuracies comparable to those that can be achieved using the current most sophisticated wave-function methods. If the dispersion energy is combined with (i) the electrostatic and first-order exchange interaction energies as defined in symmetry-adapted perturbation theory (SAPT) but computed using monomer Kohn-Sham (KS) determinants, and (ii) the induction energy computed using the coupled KS static response theory, (iii) the exchange-induction and exchange-dispersion energies computed using KS orbitals and orbital energies, the resulting method, denoted by SAPT(DFT), produces very accurate total interaction potentials. For the helium dimer, the only system with nearly exact benchmark values, SAPT(DFT) reproduces the interaction energy to within about 2% at the minimum and to a similar accuracy for all other distances ranging from the strongly repulsive to the asymptotic region. For the remaining systems investigated by us, the quality of the SAPT(DFT) interaction energies is so high that these energies may actually be more accurate than the best available results obtained with wave-function techniques. At the same time, SAPT(DFT) is much more computationally efficient than any method previously used for calculating the dispersion and other interaction energy components at this level of accuracy. [ABSTRACT FROM AUTHOR]

It is still a challenge to accurately calculate the polarizabilities of large fullerene cages and nanotubes. In this paper, a simple bond order bond polarizability relationship for carbon was found, which allowed us to apply the bond polarizability model to any pentagon isolation rule (PIR) fullerene (cage or nanotube). Following this approach, the following simple equation, α=1.262n, was obtained relating the static dipole polarizability (α) of PIR fullerenes (cages or closed nanotubes) to their number (n) of carbon atoms. Furthermore, it was shown that the polarizabilities of C60 and C70, calculated on the basis of this model, are in excellent agreement with those obtained experimentally and by density-functional theory calculations. [ABSTRACT FROM AUTHOR]

An iterative approach for calculating the frequency domain linear response of molecular systems within time-dependent density-functional theory is presented. The method completely avoids computing the exchange-correlation kernel which is typically the most expensive step for large systems. In particular, virtual orbitals are not needed. This approach may be useful for treating the response of large systems. We give an outline of the theory and a demonstration on a jellium model of an elliptic gold cluster. A detailed theory is appended discussing the computation of conductance and ac impedance of molecular junctions under bias. [ABSTRACT FROM AUTHOR]

A derivation of the density-functional-theory- (DFT) based reactivity indices in the ensemble unrestricted Hartree-Fock (eUHF) method is presented. The comparison between the properties of the reactivity indices evaluated in one and two sets of spin-orbital approach of the eUHF and hyper-unrestricted Hartree-Fock (UHF) methods are shown. All approaches give similar Fukui function irrespective of methodology used, but significantly differ for the global indices, containing important chemical information, and so their interpretation in terms of DFT- based indices can be questionable. The calculation scheme for the indices using the first- and second-order coupled perturbed eHF equations is proposed. A method for the identification of the spinorbitals involved in the change of the total number of electrons is included. The illustrative examples (water and hydrogen cyanide) show that the ground-state (GS) properties of the (Z±1)-electron systems can be predicted from the GS properties of the Z-electron systems with an accuracy comparable with the UHF calculations. The relaxation effect, important for the HCN system in which a change in the symmetry of the highest-occupied spin-orbital occurs, is effectively predicted. [ABSTRACT FROM AUTHOR]

In this paper the structures of 4-aminophenol(H2O)1+ and 3-aminophenol(H2O)1+ clusters are investigated in molecular beam experiments by different IR/UV-double resonance techniques as well as the mass analyzed threshold ionization spectroscopy yielding both inter- and intramolecular vibrations of the ionic and neutral species. Possible structures are extensively calculated at the level of density functional theory (DFT) or at the ab initio level of theory. From the experimental and theoretical investigations it can be concluded that in the case of 4-aminophenol(H2O)1 one O[Single_Bond]H...O hydrogen-bonded structure exists in the neutral cluster but two structures containing either an O[Single_Bond]H...O or a N[Single_Bond]H...O hydrogen-bonded arrangement are observed in the spectra of the ionic species. This observation is a result of an intramolecular rearrangement reaction within the ion which can only take place if high excess energies are used. A reaction path via the CH bonds is calculated and explains the experimental observations. In the case of 3-aminophenol(H2O)1+ only one O[Single_Bond]H...O bound structure is observed both in the neutral and ionic species. Ab initio and DFT calculations show that due to geometrical and energetical reasons a rearrangement cannot be observed in the 3-aminophenol(H2O)1+ cluster ion. [ABSTRACT FROM AUTHOR]

We present the case for the nonempirical construction of density functional approximations for the exchange-correlation energy by the traditional method of “constraint satisfaction” without fitting to data sets, and present evidence that this approach has been successful on the first three rungs of “Jacob’s ladder” of density functional approximations [local spin-density approximation (LSD), generalized gradient approximation (GGA), and meta-GGA]. We expect that this approach will also prove successful on the fourth and fifth rungs (hyper-GGA or hybrid and generalized random-phase approximation). In particular, we argue for the theoretical and practical importance of recovering the correct uniform density limit, which many semiempirical functionals fail to do. Among the beyond-LSD functionals now available to users, we recommend the nonempirical Perdew–Burke–Ernzerhof (PBE) GGA and the nonempirical Tao–Perdew–Staroverov–Scuseria (TPSS) meta-GGA, and their one-parameter hybrids with exact exchange. TPSS improvement over PBE is dramatic for atomization energies of molecules and surface energies of solids, and small or moderate for other properties. TPSS is now or soon will be available in standard codes such as GAUSSIAN, TURBOMOLE, NWCHEM, ADF, WIEN, VASP, etc. We also discuss old and new ideas to eliminate the self-interaction error that plagues the functionals on the first three rungs of the ladder, bring up other related issues, and close with a list of “do’s and don’t’s” for software developers and users. [ABSTRACT FROM AUTHOR]

In this paper a new scheme was proposed to calculate the intramolecular hydrogen-bonding energies in peptides and was applied to calculate the intramolecular seven-membered ring N–H...O==C hydrogen-bonding energies of the glycine and alanine peptides. The density-functional theory B3LYP/6-31G(d) and B3LYP/6-311G(d,p) methods and the second-order Mo\ller-Plesset perturbation theory MP2/6-31G(d) method were used to calculate the optimal geometries and frequencies of glycine and alanine peptides and related structures. MP2/6-311++G(d,p), MP2/6-311++G(3df,2p), and MP2/aug-cc-pVTZ methods were then used to evaluate the single-point energies. It was found that the B3LYP/6-31G(d), MP2/6-31G(d), and B3LYP/6-311G(d,p) methods yield almost similar structural parameters for the conformers of the glycine and alanine dipeptides. MP2/aug-cc-pVTZ predicts that the intramolecular seven-membered ring N–H...O==C hydrogen-bonding strength has a value of 5.54 kcal/mol in glycine dipeptide and 5.73 and 5.19 kcal/mol in alanine dipeptides, while the steric repulsive interactions of the seven-membered ring conformers are 4.13 kcal/mol in glycine dipeptide and 6.62 and 3.71 kcal/mol in alanine dipeptides. It was also found that MP2/6-311++G(3df,2p) gives as accurate intramolecular N–H...O==C hydrogen-bonding energies and steric repulsive interactions as the much more costly MP2/aug-cc-pVTZ does. [ABSTRACT FROM AUTHOR]

Based upon the optimized effective potential with the self-interaction correction, we present in this paper an alternative gradient-corrected density-functional approximation with the proper long-range behavior of the effective potential. As applied to the study of the interconfigurational energies of the whole transition-metal atoms, the present combination of the gradient-corrected contribution and the modified optimized effective potential lead the s ionization to the excellent agreement with the experiment. The calculated d ionizations and s–d transition energies are also discussed. [ABSTRACT FROM AUTHOR]

We propose a simulation method for infinitely dilute polymer solutions. In this method, an effective Hamiltonian of the solvated polymer chain is introduced to eliminate the degree of freedom of the solvent particle. The effective Hamiltonian is coupled with the density-functional theory (DFT) that we have developed for a polymer-solvent pair correlation function. All the equations proposed in this paper are derived from the first principle. This simulation method was applied to polymer chains in supercritical solvents. We observed anomalous behaviors of polymer chains near the liquid-vapor critical point: both solvophilic and solvophobic polymers expand significantly near the critical point; this is in contrast to the behavior of polymer chains in vacuum. This expansion can be interpreted as a cooperative phenomenon, which enhances the large long-wavelength density fluctuation of the solvent. [ABSTRACT FROM AUTHOR]

Orbital-free density-functional theory (OF-DFT) with modern kinetic-energy density functionals (KEDFs) is a linear scaling technique that accurately describes nearly-free-electron-like (main group) metals. In an attempt towards extending OF-DFT to transition metals, here we consider whether OF-DFT can be used effectively to study Ag, a metal with a localized d shell. OF-DFT has two approximations: use of a KEDF and local pseudopotentials (LPSs). This paper reports construction of a reasonably accurate LPS for Ag by means of inversion of the Kohn–Sham (KS) DFT equations in a bulk crystal environment. The accuracy of this LPS is determined within KS-DFT (where the exact noninteracting kinetic energy is employed) by comparing its predictions of bulk properties to those obtained from a conventional (orbital-based) nonlocal pseudopotential (NLPS). We find that the static bulk properties of fcc and hcp Ag predicted within KS-DFT using this LPS compare fairly well to those predicted by an NLPS. With the transferability of the LPS established, we then use this LPS in OF-DFT, where several approximate KEDFs were tested. We find that a combination of the Thomas–Fermi (TTF) and von Weizsäcker (TvW) functionals (TvW+0.4TTF) produces better densities than those from the linear-response-based Wang–Teter KEDF. However, the equations of state obtained from both KEDFs in OF-DFT contain unacceptably large errors. The lack of accurate KEDFs remains the final barrier to extending OF-DFT to treat transition metals. [ABSTRACT FROM AUTHOR]

Using the fluctuation-dissipation theorem (FDT) in the context of density-functional theory (DFT), one can derive an exact expression for the ground-state correlation energy in terms of the frequency-dependent density response function. When combined with time-dependent density-functional theory, a new class of density functionals results that use approximations to the exchange-correlation kernel fxc as input. This FDT-DFT scheme holds promise to solve two of the most distressing problems of conventional Kohn–Sham DFT: (i) It leads to correlation energy functionals compatible with exact exchange, and (ii) it naturally includes dispersion. The price is a moderately expensive O(N6) scaling of computational cost and a slower basis set convergence. These general features of FDT-DFT have all been recognized previously. In this paper, we present the first benchmark results for a set of molecules using FDT-DFT beyond the random-phase approximation (RPA)—that is, the first such results with fxc≠0. We show that kernels derived from the adiabatic local-density approximation and other semilocal functionals suffer from an “ultraviolet catastrophe,” producing a pair density that diverges at small interparticle distance. Nevertheless, dispersion interactions can be treated accurately if hybrid functionals are employed, as is demonstrated for He2 and HeNe. We outline constraints that future approximations to fxc should satisfy and discuss the prospects of FDT-DFT. [ABSTRACT FROM AUTHOR]

A recently developed density functional theory (DFT) for tethered bead-spring chains is used to investigate colloidal forces for the good solvent case. A planar surface of tethered chains is opposed to a bare, hard wall and the force exerted on the bare wall is calculated by way of the contact density. Previously, the case of large wall separation was investigated. The density profiles of the unperturbed chains, in that case, were found to be neither stepfunctions nor parabolas and were shown to accurately predict computer simulation results. In the present paper, the surface forces that result from the distortion of these density profiles at finite wall separation is studied. The resulting force function is analyzed for varying surface coverages, wall separations, and chain lengths. The results are found to be in near quantitative agreement with the scaling predictions of Alexander [S. Alexander, J. Phys. (Paris) 38, 983 (1977)] when the layer thickness is “correctly” defined. Finally, a hybrid Alexander–DFT theory is suggested for the analysis of experimental results. [ABSTRACT FROM AUTHOR]

This paper describes a systematic study on the clathrate structure of H+(H2O)21 using tandem mass spectrometry, vibrational predissociation spectroscopy, Monte Carlo simulations, and density functional theory calculations. We produced H+(H2O)n from a continuous corona-discharged supersonic expansion and observed three anomalies simultaneously at the cluster temperature near 150 K, including (1) the peak at n=21 is more intense than its neighboring ions in the mass spectrum, (2) the size-dependent dissociation fractions show a distinct drop for the 21-mer, and (3) the infrared spectrum of H+(H2O)21 exhibits only a single feature at 3699 cm-1, corresponding to the free-OH stretching of three-coordinated water molecules. Interestingly, the anomalies appear or disappear together with cluster temperature, indicating close correlation of these three observations. The observations, together with Monte Carlo simulations and density functional theory calculations, corroborate the notion for the formation of a distorted pentagonal dodecahedral (512) cage with a H2O molecule in the cage and a H3O+ ion on the surface for this “magic number” water cluster ion. The dodecahedral cage melts at higher temperatures, as evidenced by the emergence of a free-OH stretching feature at 3717 cm-1 for the two-coordinated water in H+(H2O)21 produced in a warmer molecular beam. Extension of this study to larger clusters strongly suggests that the experimentally observed isomer of H+(H2O)28 is most likely to consist of a distorted protonated pentakaidecahedral (51263) cage enclosing two neutral water molecules. [ABSTRACT FROM AUTHOR]

This paper proposes methods for calculating the derivative couplings between adiabatic states in density-functional theory (DFT) and compares them with each other and with multiconfigurational self-consistent field calculations. They are shown to be accurate and, as expected, the costs of their calculation scale more favorably with system size than post-Hartree-Fock calculations. The proposed methods are based on single-particle excitations and the associated Slater transition-state densities to overcome the problem of the unavailability of multielectron states in DFT which precludes a straightforward calculation of the matrix elements of the nuclear gradient operator. An iterative scheme employing linear-response theory was found to offer the best trade-off between accuracy and efficiency. The algorithms presented here have been implemented for doublet-doublet excitations within a plane-wave-basis and pseudopotential framework but are easily generalizable to other excitations and basis sets. Owing to their fundamental importance in cases where the Born-Oppenheimer separation of motions is not valid, these derivative couplings can facilitate, for example, the treatment of nonadiabatic charge transfers, of electron-phonon couplings, and of radiationless electronic transitions in DFT.© 2005 American Institute of Physics. [ABSTRACT FROM AUTHOR]

The time-dependent density functional theory (TD-DFT) calculation of excitation spectra places certain demands on the DFT exchange-correlation potential, v[sub xc], that are not met by the functionals normally used in molecular calculations. In particular, for high-lying excitations, it is crucial that the asymptotic behavior of v[sub xc] be correct. In a previous paper, we introduced a novel asymptotic-correction approach which we used with the local density approximation (LDA) to yield an asymptotically corrected LDA (AC-LDA) potential [Casida, Casida, and Salahub, Int. J. Quantum Chem. 70, 933 (1998)]. The present paper details the theory underlying this asymptotic correction approach, which involves a constant shift to incorporate the effect of the derivative discontinuity (DD) in the bulk region of finite systems, and a spliced asymptotic correction in the large r region. This is done without introducing any adjustable parameters. We emphasize that correcting the asymptotic behavior of v[sub xc] is not by itself sufficient to improve the overall form of the potential unless the effect of the derivative discontinuity is taken into account. The approach could be used to correct v[sub xc] from any of the commonly used gradient-corrected functionals. It is here applied to the LDA, using the asymptotically correct potential of van Leeuwen and Baerends (LB94) in the large r region. The performance of our AC-LDA v[sub xc] is assessed for the calculation of TD-DFT excitation energies for a large number of excitations, including both valence and Rydberg states, for each of four small molecules: N[sub 2], CO, CH[sub 2]O, and C[sub 2]H[sub 4]. The results show a significant improvement over those from either the LB94 or the LDA functionals. This confirms that the DD is indeed an important element in the design of functionals. The quality of TDLDA/LB94 and TDLDA/AC-LDA oscillator strengths were also assessed in what we believe to be the first rigorous assessment ... [ABSTRACT FROM AUTHOR]

We reply to the comment on our paper "Four-component relativistic density functional calculations of heavy diatomic molecules" [J. Chem. Phys. 112, 3499 (2000)] in which the values of our calculated spectroscopic constants for heavy diatomic molecules were reinvestigated. We point out that the improved spectroscopic constants presented by Liu and van Wu¨llen confirm our basic conclusion that the generalized gradient approximations tend to overcorrect errors of the bond lengths calculated by the local density approximation for molecules with heavy constituents. © 2000 American Institute of Physics. [ABSTRACT FROM AUTHOR]

Complete and detailed experimental transition probability density functions P(E[sup ′],E) have been determined for the first time for collisions between a large, highly vibrationally excited molecule, toluene, and several bath gases. This was achieved by applying the method of kinetically controlled selective ionization (KCSI) (Paper I [J. Chem. Phys. 112, 4076 (2000), preceding article]). An optimum P(E[sup ′],E) representation is recommended (monoexponential with a parametric exponent in the argument) which uses only three parameters and features a smooth behavior of all parameters for the entire set of bath gases. In helium, argon, and CO[sub 2] the P(E[sup ′],E) show relatively increased amplitudes in the wings--large energy gaps |E[sup ′]-E|--which can also be represented by a biexponential form. The fractional contribution of the second exponent in these biexponentials, which is directly related to the fraction of the so-called "supercollisions," is found to be very small (<0.1%). For larger colliders the second term disappears completely and the wings of P(E[sup ′],E) have an even smaller amplitude than that provided by a monoexponential form. At such low levels, the second exponent is therefore of practically no relevance for the overall energy relaxation rate. All optimized P(E[sup ′],E) representations show a marked linear energetic dependence of the (weak) collision parameter α[sub 1](E), which also results in an (approximately) linear dependence of 〈ΔE〉 and of the square root of 〈ΔE[sup 2]〉. The energy transfer parameters presented in this study form a new benchmark class in certainty and accuracy, e.g., with only 2%-7% uncertainty for our 〈ΔE〉 data below 25 000 cm-1. They should also form a reliable testground for future trajectory calculations and theories describing collisional energy transfer of polyatomic molecules. © 2000 American Institute of Physics. [ABSTRACT FROM AUTHOR]

A coupled system of two molecules bearing spins of 1/2, which are allowed to diffuse relative to each other, is considered. By using a symmetry-adapted basis operator set, the overall density matrix equation is decoupled into two equations for the time-resolved isochromat components, the sum of which yields the observed signal. The appropriate stochastic Liouville equation is solved by a combination of eigenfunction expansion and finite-differences for the angular and radial relative motions, respectively. A full range of spectra from classic Pake patterns in the rigid limit to motionally narrowed Lorentzians is recovered. As an extension of the above approach, the solid-echo experiment is described in terms of the ensemble-averaged isochromats. Homogeneous transverse relaxation times (T[sub 2]) as a function of the translational diffusion coefficient (D[sub T]) are obtained from simulating SECSY (spin-echo correlation spectroscopy) signals, which show distinct T[sub 2] minima vs D[sub T]. The present method of separating the quantum and stochastic classical variables constitutes a useful approach for treating multiquantum statistical systems, and it can be generalized to an arbitrary number of spins as shown in a companion paper. In the present study we obtained the usual linear dependence of T[sub 2] on D[sub T] in the motional narrowing (or Redfield) limit, whereas in the slow motional regime a D[sub T][sup -1/2] dependence is observed, consistent with studies of rotational diffusion. Varying the distance of maximum separation between the two spins (r[sub max]) has virtually no effect on the location of the T[sub 2] minimum with respect to D[sub T], implying that the onset of slow motion is essentially independent of r[sub max]. © 2000 American Institute of Physics. [ABSTRACT FROM AUTHOR]

Hybrid molecular dynamics (MD) simulations, which combine density functional theory (DFT) descriptions of a molecule with a molecular mechanics (MM) modeling of its solvent environment, have opened the way towards accurate computations of solvation effects in the vibrational spectra of molecules. Recently, Wheeler et al. [ChemPhysChem 4, 382 (2002)] have suggested to compute these spectra from DFT/MM-MD trajectories by diagonalizing the covariance matrix of atomic fluctuations. This so-called principal mode analysis (PMA) allegedly can replace the well-established approaches, which are based on Fourier transform methods or on conventional normal mode analyses. By scrutinizing and revising the PMA approach we identify five conditions, which must be guaranteed if PMA is supposed to render exact vibrational frequencies. Besides specific choices of (a) coordinates and (b) coordinate systems, these conditions cover (c) a harmonic intramolecular potential, (d) a complete thermal equilibrium within the molecule, and (e) a molecular Hamiltonian independent of time. However, the PMA conditions [(c)–(d)] and [(c)–(e)] are generally violated in gas phase DFT-MD and liquid phase DFT/MM-MD trajectories, respectively. Based on a series of simple analytical model calculations and on the analysis of MD trajectories calculated for the formaldehyde molecule in the gas phase (DFT) and in liquid water (DFT/MM) we show that in both phases the violation of condition (d) can cause huge errors in PMA frequency computations, whereas the inevitable violations of conditions (c) and (e), the latter being generic to the liquid phase, imply systematic and sizable underestimates of the vibrational frequencies by PMA. We demonstrate that the huge errors, which are caused by an incomplete thermal equilibrium violating (d), can be avoided if one introduces mode-specific temperatures Tj and calculates the frequencies from a “generalized virial” (GV) expression instead from PMA. Concerning ways to additionally remove the remaining errors, which GV still shares with PMA, we refer to Paper II of this work [M. Schmitz and P. Tavan, J. Chem. Phys. 121, 12247 (2004)]. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

The midinfrared (MIR) spectra of molecules in polar solvents exhibit inhomogeneously broadened bands whose spectral positions are shifted as compared to the gas phase. The shifts are caused by interactions with structured solvation shells and the broadenings by fluctuations of these interactions. The MIR spectra can be calculated from hybrid molecular dynamics (MD) simulations, which treat the solute molecule by density functional theory and the solvent by molecular mechanics by the so-called instantaneous normal mode analysis (INMA) or by Fourier transforming the time correlation function (FTTCF) of the molecular dipole moment. In Paper I of this work [M. Schmitz and P. Tavan, J. Chem. Phys. 121, 12233 (2004)] we explored an alternative method based on generalized virial (GV) frequencies noting, however, that GV systematically underestimates frequencies. As shown by us these artifacts are caused by solvent-induced fluctuations of the (i) equilibrium geometry, (ii) force constants, and (iii) normal mode directions as well as by (iv) diagonal and (v) off-diagonal anharmonicities. Here we now show, by analyzing the time scales of fluctuations and sample MD trajectories of formaldehyde in the gas phase and in water, that all these sources of computational artifacts can be made visible by a Fourier analysis of the normal coordinates. Correspondingly, the error sources (i) and (iii)–(v) can be removed by bandpass filtering, as long as the spectral signatures of the respective effects are well separated from the fundamental band. Furthermore, the artifacts arising from effect (ii) can be strongly diminished by a time-resolved version of the GV approach (TF-GV). The TF-GV method then yields for each mode j a trajectory of the vibrational frequency ωj(t|τ) at a time resolution τ>τj, which is only limited by the corresponding oscillation time τj=2π/ωj and, thus, is in the femtosecond range. A correlation analysis of these trajectories clearly separates the librational motions from the conformational dynamics of the solvation shells and yields the inhomogeneously broadened MIR spectra, if the theory of motional narrowing is properly included. The MIR spectrum of formaldehyde in solution obtained by TF-GV agrees very well with the FTTCF result, if one applies the so-called “harmonic approximation” quantum correction factor and a temperature scaling to the FTTCF intensities. Also for INMA an excellent agreement is achieved if one disregards a slight INMA overestimate of linewidths.© 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

Recently an inorganic fullerine-like [As@Ni12@As20]3- onion with near-perfect icosahedral symmetry in the crystalline phase was reported [M. J. Moses, J. C. Fettinger, and B. W. Eichhorn, Science 300, 778 (2003)]. This paper presents a detailed computational study in the framework of density functional theory on various aspects of this molecule. The electronic structure of the As@Ni12@As20 is investigated in its neutral as well as -3 charged state together with its subunits As20 and As@Ni12 by the all electron linear combination of Gaussian-type orbitals method. The bonding is studied by examining the integrated charge within atomic sphere, the electron localization function, changes in the electron density distribution, and from vibrational modes. We find that strong covalent As-As bonds seen in isolated As20 become weaker in the As@Ni12@As20 and strong covalent As-Ni bonds are formed. The structural stability of all four clusters is examined by analyzing the energetics and by calculating the vibrational frequencies. Further, the infrared and Raman spectra is predicted for both the neutral and charged As@Ni12@As20 clusters. Finally, the energy barrier for removal of a single arsenic atom is calculated for the neutral As@Ni12@As20 cluster. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

Density functional theory has been used to investigate the binding of propene to small Ag clusters in the gas phase. The binding mechanism based on frontier orbital theory, which we used previously to describe the binding between propene and the Au clusters, works for the pure Ag clusters as well. Among other things, it explains the trends of the desorption energy of propene as a function of the Ag cluster size. We show that one can predict the binding site of propene by examining the shape of the lowest unoccupied molecular orbitals (LUMOs) of the bare clusters and correlate the strength of the bond to the orbital energies of the LUMOs of the bare cluster. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

Most present applications of time-dependent density functional theory use adiabatic functionals, i.e., the effective potential at time t is determined solely by the density at the same time. This paper discusses a method that aims to go beyond this approximation, by incorporating “memory” effects: the derived exchange-correlation potential will depend not only on present densities but also on the past. In order to ensure the potentials are causal, we formulate the action on the Keldysh contour for electrons in electromagnetic fields, from which we derive suitable Kohn–Sham equations. The exchange-correlation action is now a functional of the electron density and velocity field. A specific action functional is constructed which is Galilean invariant and yields a causal exchange-correlation vector potential for the Kohn–Sham equations incorporating memory effects. We show explicitly that the net exchange-correlation Lorentz force is zero. The potential is consistent with known dynamical properties of the homogeneous electron gas (in the linear response limit). © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

In the previous paper it was shown that in the ground state the diagonal of the spin independent second-order density matrix n can be determined by solving a single auxiliary equation of a two-particle problem. Thus the problem of an arbitrary system with even electrons can be reduced to a two-particle problem. The effective potential of the two-particle equation contains a term vp of completely kinetic origin. Virial theorem and hierarchy of equations are derived for vp and simple approximations are proposed. A relationship between the effective potential up of the shape function equation and the potential vp is established. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

Gas-phase hydrogen atoms are accelerated towards metallic surfaces in their vicinity. As it approaches the surface, the velocity of an atom increases and this motion excites the metallic electrons, causing energy loss to the atom. This dissipative dynamics is frequently described as atomic motion under friction, where the friction coefficient is obtained from ab initio calculations assuming a weak interaction and slow atom. This paper tests the aforementioned approach by comparing to a real-time Ehrenfest molecular dynamics simulation of such a process. The electrons are treated realistically using standard approximations to time-dependent density functional theory. We find indeed that the electronic excitations produce a frictionlike force on the atom. However, the friction coefficient strongly depends on the direction of the motion of the atom: it is large when the atom is moving towards the cluster and much smaller when the atom is moving away. It is concluded that a revision of the model for energy dissipation at metallic surfaces, at least for clusters, may be necessary. [ABSTRACT FROM AUTHOR]

Orbital-free density functional theory as an extension of traditional Thomas-Fermi theory has attracted a lot of interest in the past decade because of developments in both more accurate kinetic energy functionals and highly efficient numerical methodology. In this paper, we developed a conjugate-gradient method for the numerical solution of spin-dependent extended Thomas-Fermi equation by incorporating techniques previously used in Kohn-Sham calculations. The key ingredient of the method is an approximate line-search scheme and a collective treatment of two spin densities in the case of spin-dependent extended Thomas-Fermi problem. Test calculations for a quartic two-dimensional quantum dot system and a three-dimensional sodium cluster Na216 with a local pseudopotential demonstrate that the method is accurate and efficient. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]