Projection of plane-wave calculations into atomic orbitals

The projection of the eigenfunctions obtained in standard plane-wave first-principle electronic-structure calculations into atomic-orbital basis sets is proposed as a formal and practical link between the methods based on plane waves and the ones based on atomic orbitals. Given a candidate atomic basis, ({\it i}) its quality is evaluated by its projection into the plane-wave eigenfunctions, ({\it ii}) it is optimized by maximizing that projection, ({\it iii}) the associated tight-binding Hamiltonian and energy bands are obtained, and ({\it iv}) population analysis is performed in a natural way. The proposed method replaces the traditional trial-and-error procedures of finding appropriate atomic bases and the fitting of bands to obtain tight-binding Hamiltonians. Test calculations of some zincblende semiconductors are presented.Comment: RevTex. 4 pages. 3 uuencoded compressed (tared) postscript figs. To appear in Solid St. Commu

A new projection method for finding the closest point in the intersection of convex sets

In this paper we present a new iterative projection method for finding the closest point in the intersection of convex sets to any arbitrary point in a Hilbert space. This method, termed AAMR for averaged alternating modified reflections, can be viewed as an adequate modification of the Douglas--Rachford method that yields a solution to the best approximation problem. Under a constraint qualification at the point of interest, we show strong convergence of the method. In fact, the so-called strong CHIP fully characterizes the convergence of the AAMR method for every point in the space. We report some promising numerical experiments where we compare the performance of AAMR against other projection methods for finding the closest point in the intersection of pairs of finite dimensional subspaces

Simulations of water nano-confined between corrugated planes

Two-dimensionally nanoconfined water between ideal planar walls has been the subject of ample study, aiming at understanding the intrinsic response of water to confinement, avoiding the consideration of the chemistry of actual confining materials. In this work, we study the response of such nanoconfined water under a periodic confining potential by means of computer simulations, both using empirical potentials and from first-principles. We propose a periodic confining potential emulating the atomistic oscillation of the confining walls, which allows varying the lattice parameter and amplitude of the oscillation. We do it for a triangular lattice, with several values of the lattice parameter: one which is ideal for commensuration with layers of Ih ice, and other values that would correspond to more realistic substrates. For the former, an overall rise of the melting temperature is observed. The liquid maintains a bi-layer triangular structure, however, despite the fact that it is not favoured by the external periodicity. The first-principles liquid is significantly affected by the modulation in its layering and stacking even at relatively small amplitudes of the confinement modulation. Beyond some critical modulation amplitude the hexatic phase present in flat confinement is replaced by a trilayer crystalline phase unlike any of the phases encountered for flat confinement. For more realistic lattice parameters, the liquid does not display higher tendency to freeze, but it clearly shows inhomogeneous behaviour as the strength of the rugosity increases. In spite of this expected inhomogeneity, the structural and dynamical response of the liquid is surprisingly insensitive to the external modulation. Although the first-principles calculations give a more triangular liquid than the one observed with empirical potentials, both agree remarkably well for the main conclusions of the study.Comment: 13 pages, 16 figure

Short range repulsive interatomic interactions in energetic processes in solids

The repulsive interaction between two atoms at short distances is studied in order to explore the range of validity of standard first-principles simulation techniques and improve the available short-range potentials for the description of energetic collision cascades in solids. Pseudopotentials represent the weakest approximation, given their lack of explicit Pauli repulsion in the core-core interactions. The energy (distance) scale realistically accessible is studied by comparison with all-electron reference calculations in some binary systems. Reference calculations are performed with no approximations related to either core (frozen core, augmentation spheres) or basis set. This is important since the validity of such approximations, even in all-electron calculations, rely on the small core perturbation usual in low-energy studies. The expected importance of semicore states is quantified. We propose a scheme for improving the electronic screening given by pseudopotentials for very short distances. The results of this study are applied to the assessment and improvement of existing repulsive empirical potentials.Comment: 10 pages, 7 figure

Two exchange-correlation functionals compared for first-principles liquid water

The first-principles description of liquid water using ab initio molecular dynamics (AIMD) based on Density Functional theory (DFT) has recently been found to require long equilibration times, giving too low diffusivities and a clear over-structuring of the liquid. In the light of these findings we compare here the room-temperature description offered by two different exchange correlation functionals: BLYP, the most popular for liquid water so far, and RPBE, a revision of the widely used PBE. We find for RPBE a less structured liquid with radial distribution functions closer to the experimental ones than the ones of BLYP. The diffusivity obtained with RPBE for heavy water is still 20% lower than the corresponding experimental value, but it represents a substantial improvement on the BLYP value, one order of magnitude lower than experiment. These characteristics and the hydrogen-bond (HB) network imperfection point to an effective temperature ~3% lower than the actual simulation temperature for the RPBE liquid, as compared with BLYP's ~17% deviation. The too long O--O average nearest-neighbor distance observed points to an excessively weak HB, possibly compensating more fundamental errors in the DFT description.Comment: Jorunal reference adde

Resistive and rectifying effects of pulling gold atoms at thiol-gold nano-contacts

We investigate, by means of first-principles calculations, structural and transport properties of junctions made of symmetric dithiolated molecules placed between Au electrodes. As the electrodes are pulled apart, we find that it becomes energetically favorable that Au atoms migrate to positions between the electrode surface and thiol terminations, with junction structures alternating between symmetric and asymmetric. As a result, the calculated $\emph{IV}$ curves alternate between rectifying and non-rectifying behaviors as the electrodes are pulled apart, which is consistent with recent experimental results