deletions | additions       diff --git a/Introduction.tex b/Introduction.tex  index fe81112..9d205bd 100644  --- a/Introduction.tex  +++ b/Introduction.tex 
...
\section{Introduction}  The development of theoretical methods for the simulation of electronic system is an active area of study. This interest has been fueled, on one hand, by the development success  of theoretical tools, tools  like density functional theory (DFT)~\cite{Hohenberg_1964,Kohn_1965}, that can predict many properties with good accuracy at a relatively modest computational cost. On the other hand,because  these same tools are not good enough for many applications, and better more accurate and more efficient  methodsin terms of new properties, higher accuracy or smaller numerical cost  are required. Current research is targeted on a broad range of aspects of electronic structure simulations: the development of new theoretical frameworks, new or improved methods to calculate properties within existing theories, or even new, more efficienent, efficient,  algorithms. In most cases, this theoretical work requires the development of test implementations to assess the properties and predictive power of the new methods.While the implementation of established methods is focused towards robustness, efficiency and generality. Given the experimentative nature of this work, the translation to code of new theory needs to be easy to implement and to modify. This is a factor that is not usually considered when analyzing and comparing the broad range of methods and codes used by chemists, phyisicst and material scientist to simulate electrons.  When simulating electrons different fields needs to be represented numerically, Given the experimentative nature of the development of methods  forexample  the ionic potential, simulations of electrons,  the single-particle orbitals or states, or translation to code of new theory needs to be easy to implement and to modify. This is a factor that is not usually considered when analyzing and comparing  the electronic density. broad range of methods and codes used by chemists, phyisicst and material scientist.  The most popular representations rely on basis sets, that usually have a certain physical connection to the system being simulated. In chemistry the method of choice is to use atomic orbitals as a basis to describe the orbitals of a molecule. When these atomic orbitals are expanded in Gaussian functions, it leads to an efficient method as many integrals can be calculated from analytical formulae~\cite{szabo1996modern}. In condensed matter physics, on the other hand, the traditional basis is a set of plane waves, that correspond to the eigenstates of a homogeneous electron gas. These physics-inspired basis sets have, however, some limitations. For example, it is not trivial to simulate crystalline systems using atomic orbitals~\cite{Dovesi_2014}, and, on the other hand, in plane wave approaches finite systems must be approximated as periodic system using a super cell approach. Several alternatives to atomic-orbital and plane-wave basis sets exist~\cite{Harrison_2004,Pask_2005,16008435,Genovese_2011}. One particular approach that does not depend on a basis set is to represent fields directly on a real-space mesh. The method was pioneered by Becke~\cite{Becke_1989}, who used a combination of radial-grids centered around each atom. In 1994 Chelikowsky, Troullier and Saad~\cite{Chelikowsky_1994} presented a practical approach for the solution of the Kohn-Sham equations using uniform grids combined with pseudo-potentials. What made the approach competitive was the use of high-order finite differences, that keep control the error of the Laplacian without requiring very dense meshes. From that moment, several real-space implementation have been presented~\cite{Briggs_1996,Fattebert_1996,Beck_1997,Nardelli_2001,Marques_2003,Kronik_2006,Enkovaara_2010}.