Ab initio many-body calculations on infinite carbon and boron-nitrogen chains

In this paper we report first-principles calculations on the ground-state electronic structure of two infinite one-dimensional systems: (a) a chain of carbon atoms and (b) a chain of alternating boron and nitrogen atoms. Meanfield results were obtained using the restricted Hartree-Fock approach, while the many-body effects were taken into account by second-order M{\o}ller-Plesset perturbation theory and the coupled-cluster approach. The calculations were performed using 6-31$G^{**}$ basis sets, including the d-type polarization functions. Both at the Hartree-Fock (HF) and the correlated levels we find that the infinite carbon chain exhibits bond alternation with alternating single and triple bonds, while the boron-nitrogen chain exhibits equidistant bonds. In addition, we also performed density-functional-theory-based local density approximation (LDA) calculations on the infinite carbon chain using the same basis set. Our LDA results, in contradiction to our HF and correlated results, predict a very small bond alternation. Based upon our LDA results for the carbon chain, which are in agreement with an earlier LDA calculation calculation [ E.J. Bylaska, J.H. Weare, and R. Kawai, Phys. Rev. B 58, R7488 (1998).], we conclude that the LDA significantly underestimates Peierls distortion. This emphasizes that the inclusion of many-particle effects is very important for the correct description of Peierls distortion in one-dimensional systems.Comment: 3 figures (included). To appear in Phys. Rev.

Work functions, ionization potentials, and in-between: Scaling relations based on the image charge model

We revisit a model in which the ionization energy of a metal particle is associated with the work done by the image charge force in moving the electron from infinity to a small cut-off distance just outside the surface. We show that this model can be compactly, and productively, employed to study the size dependence of electron removal energies over the range encompassing bulk surfaces, finite clusters, and individual atoms. It accounts in a straightforward manner for the empirically known correlation between the atomic ionization potential (IP) and the metal work function (WF), IP/WF$\sim$2. We formulate simple expressions for the model parameters, requiring only a single property (the atomic polarizability or the nearest neighbor distance) as input. Without any additional adjustable parameters, the model yields both the IP and the WF within $\sim$10% for all metallic elements, as well as matches the size evolution of the ionization potentials of finite metal clusters for a large fraction of the experimental data. The parametrization takes advantage of a remarkably constant numerical correlation between the nearest-neighbor distance in a crystal, the cube root of the atomic polarizability, and the image force cutoff length. The paper also includes an analytical derivation of the relation of the outer radius of a cluster of close-packed spheres to its geometric structure.Comment: Original submission: 8 pages with 7 figures incorporated in the text. Revised submission (added one more paragraph about alloy work functions): 18 double spaced pages + 8 separate figures. Accepted for publication in PR

Nucleation of a sodium droplet on C60

We investigate theoretically the progressive coating of C60 by several sodium atoms. Density functional calculations using a nonlocal functional are performed for NaC60 and Na2C60 in various configurations. These data are used to construct an empirical atomistic model in order to treat larger sizes in a statistical and dynamical context. Fluctuating charges are incorporated to account for charge transfer between sodium and carbon atoms. By performing systematic global optimization in the size range 1<=n<=30, we find that Na_nC60 is homogeneously coated at small sizes, and that a growing droplet is formed above n=>8. The separate effects of single ionization and thermalization are also considered, as well as the changes due to a strong external electric field. The present results are discussed in the light of various experimental data.Comment: 17 pages, 10 figure

Theoretical Evidence for Temperature-induced Proton Mobility in Isolated Lysine-rich Polyalanines

A multistate molecular mechanics method is introduced to model the possible competition between various protonation sites in gas-phase biomolecules with excess protons. The method relies on the Amber force field for each site and is calibrated against density-functional theory benchmark calculations at the 6−31+G(d,p) level. In its adiabatic version, where it has similarities with constant-pH algorithms, the model predicts that the small protonated Ala4−Lys peptide, unprotected at the N-terminus, changes protonation site above 400 K. In the larger [Ala9−Lys+H ]+ peptide, the proton remains at the lysine amine group in a favored charge/electric dipole conformation. In the three-state Ala4−Lys−Ala4−Lys peptide, the excess proton is found to be partially delocalized on the amine group of the first lysine and on the N-terminus. The statistical properties of the protonated peptides are found to significantly depend on the localized character of the proton. Finally, the model is extended by considering possible couplings between the protonation sites, in an empirical valence-bond version. Strong couplings can stabilize the peptides into unexpected proton-bound conformations over broad ranges of temperature

Photoionization of LinHm clusters

Ionization potentials of LinHm clusters have been measured from bare Lin clusters to hydrogen saturated clusters. The evolution of electronic properties with the number of H is discussed. We found that LinHm clusters behave like Lin-m clusters. This similarity may be due to a segregation between a metallic part and an insulator part inside the cluster