\section{Introduction}
Recent advances in materials science and technology have witnessed a rapid expansion of research towards
discoveries of novel materials for sustainable energy. In this cause, hydrogen has been identified as an
alternative energy carrier because of its high abundance, lightweight, and environmental
friendliness~.\cite{sch-01} Yet, there are still some important challenges to be addressed prior to possible
actual spread of hydrogen-based infrastructure. One of the critical problems is lack of suitable hydrogen
storage materials which can store H$_{2}$ with high volumetric and gravimetric density, and also desorbs
H$_{2}$ through easy means. Furthermore, hydrogen adsorption and desorption should be operated under near
ambient conditions in cost effective ways.~\cite{meng-07, xiong-08} Since the first experiment of hydrogen
adsorption on carbon nanotubes by Dillon \textit{et al.}~\cite{nature-Dil97} and announcement of hydrogen
storage target of 6$\%$ by weight for on-board automotive applications by the US Department of
Energy~(DOE)\cite{heden-2000}, carbon nanotubes are considered as suitable candidates because of their light
weight, high strength, high adsorption characteristics, and also large surface to volume ratio. However, experimental
findings showed that carbon nanotubes are not suitable for the H$_{2}$ adsorption due to weak van der Waal
interaction between H$_{2}$ and carbon surface.\cite{nature-Dil97} These works have triggered extensive research
on the functionalization of the carbon nanotubes by incorporating transition metals as well as alkali and
alkaline earth metals.

The study of hydrogen storage using transition (Sc, Ti, and Ni), alkali (Li, K, Na) and alkaline earth
(Mg, Ca) metals doped carbon nanotubes have been reported by many research
groups~\cite{pavel-07,prl-Cir05,jpc-Liu09,science-Che99, jcp-Cab05,jacs-sun06,prl-yoon08}. A single
transition-metal (Sc, Ti, and Ni) atom coated on C$_{60}$ or carbon nanotube surface can bind up to
three to five H$_{2}$ molecules~\cite{prl-zhao05, prl-Cir05,apl-shin06}.  Thus, the dispersed C$_{60}$ and
nanotubes can adsorb up to 6.8-10~wt~$\%$ hydrogen with high metal coverage. Carbon nanomaterials decorated
with alkali and alkaline-earth metal atoms can also be used for hydrogen storage.  Chen~\textit{et al.} have
reported a hydrogen adsorption capacity of 20~wt\%\ and 14~wt\%\ in lithium and potassium doped carbon
nanotubes\cite{science-Che99} respectively.  Later on, it was found that observation of such high adsorption
capacity was due to the presence of water impurities\cite{carbon-Yan00}. Cabria~\textit{et al.}~\cite{jcp-Cab05}
have shown that the binding energy of hydrogen molecule on Li-doped carbon nanotubes is two times higher than
that on pristine carbon nanotubes.  Recently, Liu and coworkers\cite{jpc-Liu09} found that up to eight Li metal
atoms dispersed on single walled carbon nanotubes (SWCNTs) can adsorb 64~H$_{2}$ molecules with an average
binding energy -0.17~eV/H$_{2}$ leading to 13.5~wt\%\ of storage capacity\cite{jpc-Liu09}.
Ataca~\textit{et al.}~\cite{apl-Ataca08} have studied the dispersion of Li atoms on graphene surface and
predicted the gravimetric density of 12.8~wt\%. The isolated clusters (Li$_{12}C_{60}$, Na$_{8}C_{60}$ and
Ca$_{32}C_{60}$) where metal atoms are capped onto the pentagonal and hexagonal faces of C$_{60}$, can
store H$_{2}$ up to 8.4-9.5~wt~$\%$\cite{jacs-sun06, nanolett-Gho08, prl-yoon08}. The charge transfer from
metal atoms to C$_{60}$ gives rise to the electric field surrounding the coated fullerene and enhances the
H$_{2}$ adsorption. Recently, it is found that the sodium ions have been strongly absorbed onto a hydrophobic
graphite surface via cation-$\pi$ interactions\cite{chinphys-shi11}.

We take the sodium as a coating material on SWCNTs since the sodium ion is one of the most popular ions
in medical science, chemistry, physics and particularly in biological systems. It is sixth most abundant
element in nature and most studied in the periodic table due to its simplicity towards the experimentation.
On comparing the earlier experimental~\cite{langchak-95,assglee-04} and theoretical reports with the present
study, it can be mentioned that the sodium-doped carbon nanotubes have many advantages and satisfy some of
the important requirements for the hydrogen storage materials; namely,(a) strong interaction between the
sodium ions and H$_{2}$ molecules, (b) remarkably enhanced molecular H$_{2}$ adsorption, and (c) light-weight
material with higher capacity for the adsorption of H$_{2}$ molecules. The sodium-doped SWCNTs can act as a
better hydrogen storage material at lower temperature with better hydrogen adsorption efficiency. Further, the
adsorption can even be tuned by functionalizing with other suitably chosen organic functional groups~\cite{rscsri-12}.
There are many reports on hydrogen storage using sodium for e.g., the metal complex of sodium (sodium borohydrates,
sodium alanates\cite{chemmater-Zheng08, jamchemsoc-cornelis08}), sodium-doped fullerenes~\cite{nanolett-Gho08},
adsorption of sodium ions on hydrophobic graphite surfaces\cite{chinphys-shi11}, curvature-induced hydrogen
adsorption in sodium-doped carbon nanomaterials\cite{chandra-08}. All above studies have been carried out at
molecular levels. To the best of our knowledge, there is no report on sodium doped carbon nanotube at solid
state. Despite of many reports for hydrogen storage on Na-doped carbon materials at molecular level,
there is no attempt to explain the basic understanding of the stability of sodium atoms on carbon surfaces and
its high storage capacity towards H$_{2}$ molecules. So, the interesting issues we address are; (i). will the
high storage capacity of Na-doped carbon nanotubes sustain at solid state? (ii). binding mechanism and stability
of sodium atoms on carbon nanotube surfaces, (iii). evolution of the electronic properties of Na-SWCNT-H$_{2}$ systems,
(iv). adsorption mechanism of H$_{2}$ molecules on Na-SWCNTs, (v). effects of sodium clustering on the stability
and adsorption of H$_{2}$ molecules, (vi). electrostatic interaction of sodium and H$_{2}$ molecules on carbon surface,
(vii). effects of chirality of carbon nanotubes on the stability and adsorption of sodium and H$_{2}$ molecules,
(viii). transition states search of composite systems (Na/SWCNT-H$_{2}$). To address the above issues, we carry
out extensive search for high-capacity hydrogen storage material consisting of individual sodium atom doped carbon
nanotubes at zero temperature using density functional theory (DFT). We employ both local density approximation
(LDA)~\cite{prb-Per81,prl-Cep80} and generalized gradient approximation (GGA)~\cite{prl-Per96} to estimate the
binding strength. We have also verified the LDA results with hybrid B3LYP functional~\cite{jcpbeck-88}.

Earlier studies\cite{jcp-cab08,jpcb-oka01} have found that the LDA based predictions of the physisorption
energies of H$_{2}$ on the surface of graphite and carbon nanotubes were in good agreement with experiments.
Furthermore, quantum Monte Carlo calculations show that the adsorption energy of hydrogen in carbon materials
typically lies between the calculated adsorption energy based on GGA and LDA\cite{prl-zhao05}.  A previous
\textit{ab initio} study using the M{\o}ller-Plesset second-order perturbation~(MP2) method\cite{jpcb-oka01} has
shown that GGA tends to underestimate the adsorption energy of H$_{2}$ on SWCNT whereas LDA gives reasonable
adsorption energy as compared to the MP2 calculation\cite{jcp-cab08,prl-zhao05,jpcb-zho06}. The reliability
of LDA can be ascribed to the following facts\cite{jcp-cab08}. Firstly, when the electron densities of H$_{2}$
and carbon surface overlap weakly, the nonlinearity of the exchange-correlation energy density functional produces
an attractive interaction even in the absence of electron density redistribution. Also the overestimated binding
energy by LDA\cite{prb-lug07, prb-lee08} may compensate for the insufficient account of van der Waals
interactions.\cite{jcp-cab08} In contrast, DFT calculation using a uniform GGA produced a purely repulsive
interaction.  Calculation based on the GGA-PW91 functional has predicted a repulsive interaction between H$_{2}$
and a graphene layer and also between H$_{2}$ and a (6,6) carbon nanotube\cite{prb-tad01}, which contradicts the
experimental findings\cite{jmat-sah08}. It has also been noted that LDA calculations well reproduce the empirical
interaction potentials between graphitic layers and also in the other graphitic systems for distances near
to the equilibrium separation although the LDA is not able to reproduce the long-range dispersion
interaction\cite{prb-gir02}.

Accordingly, in the present work, the conventional DFT calculations are directed toward obtaining structural
stability and electronic structure properties of  Na-SWCNT-H$_{2}$ systems at absolute temperature and pressure.
The transition states of the reaction are also studied in order to determine the energy barrier of the reaction
and the reaction rate. Semiconducting zigzag nanotubes show stronger binding energies for the Na atom, as
compared to metallic SWCNTs with armchair and chiral geometries. We show that the charge-induced dipole interaction
is responsible for the higher hydrogen uptake of Na-coated SWCNTs. We find that the storage capacity of hydrogen
depends strongly on the clustering of Na atoms on SWCNTs.
\section{Computational Details}
All the energy calculations and geometry optimizations were performed at absolute temperature by using density
functional theory as implemented in the CASTEP code~\cite{zfk-Cla05}. We use the local density
approximation~(LDA) with Perdew-Zunger parameterization~\cite{prb-Per81} of Ceperley-Alder data~\cite{prl-Cep80}
as well as generalized gradient approximation~(GGA) by Perdew-Burke-Ernzerhof~(PBE)~\cite{prl-Per96}. Ultrasoft
pseudopotential with a plane wave basis cutoff of 380~eV was used. To ensure the calculated results being
mutually comparable, the identical conditions are employed for the isolated H$_{2}$ molecules, Na atom,
the carbon nanotube, and also the adsorbed carbon nanotube system. The Monkhorst-Pack scheme with 1x1x7 special
K-points is used for Brillouin zone sampling for all the systems. The geometries were considered converged
when the change in energy was of the order of 10$^{-5}$~Hartrees and the absolute maximum force was
0.002~Hartree/\AA.
\begin{figure}
\center
\includegraphics[scale=0.2]{fig1.eps}
\caption{Four possible Na adsorption sites on~(8,0) SWCNT. Symbols A, C, T, and H
denote the bridge site over an axial bond, the bridge over a zigzag C-C bond,
top of the carbon atom, and the hollow site of the hexagonal C-ring.
The light green and pink spheres represent carbon and sodium atoms.}
\label{fig-cnt-site}
\end{figure}

We have selected several carbon nanotubes with different chiral indices as the model of carbon substrate.
All simulations were carried out in periodicity along Z-axis. The dimensions of simulation cells ranged
from 16$\times$16$\times$8.5 for smallest system to 27.4$\times$28$\times$7.38~\AA$^{3}$
for the largest systems respectively. The simulation cell were taken so that the empty space between
the nanotubes and the X and Y directions were $\sim$~7.5~\AA. Further increase in this vacuum space
leads to negligible change in the total energy of the systems and forces on the atoms. For most part of
simulations, we have selected (8,0) SWCNT because of its feasibility for synthesis.
\begin{figure}
\center
\includegraphics[scale=0.15]{fig2.eps}
\caption{Optimized geometries of Na adsorbed at hollow sites on (3,3), (4,2),
(5,0), (5,5) (6,3), and (8,0) SWCNTs. E$_{b}$ represents the binding energy of Na.
The light green and pink spheres represent carbon and Na atoms, respectively.}
\label{fig-cnt-na}
\end{figure}

A storage model composed of Na-doped (8,0) SWCNT is built by considering four possible Na adsorption
sites. These four possible sites are the top of the carbon atom (T), the hollow site of the hexagon
C-ring (H), the bridge site over an axial bond (A), and the bridge over a zigzag C-C bond (C) as shown
in figure 1.  After full structural optimization, the Na atom always locate on hollow site (H) shown
in figure 2, regardless of the initial location. The binding energy of the Na atom is given by,
\begin{equation}
BE_{Na} = E_{Na-SWCNT}-(E_{SWCNT} + E_{Na})
\end{equation}
where E$_{Na-SWCNT}$, E$_{SWCNT}$ and E$_{Na}$ are the total energies of Na-doped SWCNT, SWCNT and
isolated Na atom respectively. Next, we used lowest energy structure of Na-SWCNT for H$_{2}$ adsorption.
The binding energy per H$_{2}$ molecule is given by,
\begin{equation}
BE_{H_2} = [E_{(H_2)_n Na-SWCNT} - E_{Na-SWCNT}-nE_{H_2}]/n
\end{equation}
where n indicates the number of H$_{2}$ molecules, E$_{(H_2)_n -Na-SWCNT}$, E$_{Na-SWCNT}$ and E$_{H_2}$ are the
total energy of (H$_{2})_n $-Na-SWCNT complex, Na-doped SWCNT and H$_{2}$ molecule respectively. The adsorption
mechanisms are analyzed in terms of the Mulliken population analysis, electrostatic potential and partial
density of states (PDOS) which can provide a clear and definitive description for charge distribution. We
have verified the accuracy of our calculations by calculating the density of states (DOS) of (8,0)SWCNT and
the results agree well with reported results of Li \textit{et al}.~\cite{jcp-li03}
\section{Results and Discussions}
\label{sec-results}
First, we consider the adsorption energetics of a single Na atom on the (8,0) SWCNT as shown in figure 2.
We have performed structural optimizations starting with four possible metal adsorption sites (A, C, H, and
T) on carbon nanotubes. A single Na atom is found to always occupy the hollow site of a hexagonal carbon
ring in all the six different CNTs considered, with Na-C bond length of 2.4-2.9 \AA.
\begin{figure}
\center
\includegraphics[scale=0.8]{fig3.eps}
\caption{The density of states of (8,0) SWCNT and SWCNT-Na are shown. The partial density of states of
carbon~(2$p$) and sodium~(3$s$) atoms of (8,0)~SWCNT-Na system are also shown. The Fermi level is set to
zero and indicated by the dotted line.}
\label{fig-na-pdos}
\end{figure}
When a Na atom is adsorbed on SWCNT, there are two nearest and four next-nearest neighbors. It was believed
that the stability of adsorption was enhanced by the curvature of the carbon nanotube surface\cite{chandra-08}.
Here, we have considered six SWCNTs [(3,3), (4,2), (5,0), (5,5), (6,3) and (8,0)]. It is observed that
E$_{b}$ for Na atoms on (5,0) and (8,0) tubes are -2.59~eV and -2.05~eV, while E$_{b}$ for Na atoms on the
(3,3), (4,2), (5,5) and (6,3) are -1.54~eV, -2.51~eV, -2.11~eV and -2.65~eV respectively. The binding energy
difference comes from the configuration of the surface. It is observed that Na atoms on semiconducting SWCNTs
[(4,2), (5,0),(6,3)and (8,0)] are more strongly adsorbed by $\sim$0.4~eV as compared with metallic
SWCNTs [(3,3), (5,5)] due to strong electrostatic interaction between Na and SWCNTs.
\begin{figure}
\center
\includegraphics[scale=0.2]{fig4.eps}
\caption{The coverage of the H$_{2}$ molecules on (8,0) SWCNT-Na.
The light green and pink spheres have the same meaning Figure 2,
and the white one shows H atom.}
\label{fig-na-h2}
\end{figure}

Now, we address the underlying binding mechanism between Na atom and (8,0) SWCNT. A single Na atom prefers to adsorb
strongly at the H site, as shown in figure 2. The Na-C distances are 2.46~\AA~ and 2.70~\AA~for nearest and next-nearest
neighbor carbon atoms with the binding strength of E$_{b}$=-2.05~eV. The Mulliken charge analysis shows that Na
carries a 0.82$e$ positive charge, indicating that the Na atom is ionized. The two nearest and four next-nearest
carbon atoms are negatively charged with 0.18$e$ and 0.035$e$ respectively, whereas remaining charges are
distributed over the whole SWCNT. The partial density of states (PDOS) in figure 3 shows a lowering of carbon
eigenvalues by 0.77~eV for all the carbon 2$p$ orbitals. The charge around Na and the PDOS implies a simple
electrostatic attraction as a binding mechanism between charged Na atom and SWCNT.
\begin{figure}
\center
\includegraphics[scale=0.8]{fig5.eps}
\caption{The density of states and the partial density of states for C~(2$p$), Na(3s) and H$_{2}$(1s)
of (8,0) SWCNT-Na-H$_{2}$ are shown. The Fermi level is set to zero and indicated by the dotted line.}
\label{fig-na-h2-pdos}
\end{figure}

Next, we investigate the interaction between Na-SWCNT with H$_{2}$ molecules. A single H$_{2}$ molecule
sits atop the Na atom in (8,0) SWCNT-Na-H$_{2}$. The optimized structure is shown in figure 4. The molecular
hydrogen is found to interact with Na atoms as a dihydrogen molecule with a binding energy -0.18~eV by LDA
and Na-H bond length 2.44~\AA. Mulliken population analysis shows  charges of +0.76$e$ on the sodium and
+0.03$e$ on the hydrogen. From the total density of states (TDOS) and PDOS given in figure 5, the hydrogen
bonding orbitals at $\sim$-9.5~eV far below the entire carbon 2$p$ band, hence there is no further lowering
of carbon 2$p$ eigenvalues as compared to the SWCNT-Na system.
\begin{table}
\center
\small
\begin{center}
\caption{The bond length, charges on Na and binding Energy (eV)
of Na and H$_{2}$ on (8,0) SWCNT-Na}
\begin{tabular}{|c|c|c|c|c|c|c|c|c|}
\hline
No. of H$_2$ &  \multicolumn{2}{c|}{C-Na bond length (\AA)} &  \multicolumn{2}{c|}{Na-H bond length (\AA)}
& Charge on Na (e)  &       \multicolumn{3}{c|}{B.E. (eV)}       \\      \hline
molecules  & Minimum       & Maximum      & Minimum      &       Maximum &       &       LDA     &
GGA     &       B3LYP   \\      \hline
0     &       2.50   &       2.70   &       -       &       -       &       0.82
&       -2.05   &       -1.32   &       -1.38   \\
1    &       2.44   &       2.69   &       2.44   &       2.44   &
0.76   &       -0.18   &       -0.16   &       -0.23   \\
2    &       2.46   &       2.70   &       2.31   &       2.34   &
0.68  &       -0.27   &       -0.16   &       -0.21   \\
3    &       2.48   &       2.72   &       2.33   &       2.35   &
0.60   &       -0.26   &       -0.14   &       -0.21   \\
4    &       2.48   &       2.76   &       2.34   &       2.41   &
0.55   &       -0.26   &       -0.13   &       -0.21   \\
5    &       2.49   &       2.78   &       2.34   &       2.46   &
0.46   &       -0.26   &       -0.12   &       -0.21   \\
6    &       2.59   &       2.87   &       2.34   &       2.52   &
0.36   &       -0.26   &       -0.12   &       -0.20   \\       \hline
\end{tabular}
\end{center}
\end{table}

We also compute the electrostatic potential of SWCNT-Na-H$_{2}$ system to understand the binding mechanism
of H$_{2}$. The substantial charge redistribution upon Na coating leads to high a electric field near the Na
atoms. We measure the electric fields, induced by the charge redistribution along the outward radial direction of
the complex. Figure 6 shows the radial component of the electric field from the center of the SWCNT to the
hollow site of SWCNT, Na atom and H$_{2}$ molecule. For the SWCNT, the electric field strength is
8.3x10$^{10}$ V/m, at the hollow site; with Na adsorbed, E(Na site) is 29.2x10$^{10}$ V/m; with single
H$_{2}$ coverage atop Na, the E~(Na site) drops to 24x10$^{10}$~V/m and at the center of H$_{2}$ molecule
E(H$_{2}$) is 8.6x10$^{10}$ V/m. Therefore, we assume that the electric field produced by SWCNT-Na polarizes
the H$_{2}$ molecules and the induced dipole-dipole as well as the charge-dipole interactions, consequently,
leads to binding of H$_{2}$ around SWCNT-Na. The PDOS of SWCNT-Na-H$_{2}$ shows that there is no bonding
orbital between H$_{2}$ and Na, which also supports this argument.
\begin{figure}
\center
\includegraphics[scale=0.9]{fig6.eps}
\caption{The radial component of the electric field of (8,0) SWCNT, SWCNT-Na and SWCNT-Na-H$_{2}$ systems.
The dotted vertical lines indicate the position of Carbon ring, Na atom and H$_{2}$ molecules. The insets
shows the isosurfaces of the electrostatic potentials of SWCNT, SWCNT-Na and SWCNT-Na-H$_{2}$.}
\label{fig-na-h2-esp}
\end{figure}

The optimized geometries of higher coverage of H$_{2}$ molecules on the (8,0) SWCNT-Na are shown in figure 4.
The hydrogen binding energies (E$_{b}$) are summarized in Table 1, which show the GGA, LDA and B3LYP results
for both Na atom and H$_{2}$ molecules. A single Na atom is found to adsorb up to six~H$_{2}$~molecules with
a relatively constant binding energy of~-0.26~eV/H$_{2}$. It may be interesting to note that the maximum number
of H$_{2}$ adsorbed at the Na site is very much similar to the case of Na-doped fullerene. Electrons in a
carbon nanotube undergo major changes in their charge distribution upon Na doping. The transferred charges
are spread over the hexagon formed out of carbon atoms nearest to the Na atom. The significant charge
redistribution results in an electric field strong enough to attract up to six H$_{2}$ molecules with a
binding energy of $\sim$-0.26~eV per H$_{2}$ with LDA and -0.14 with GGA functionals. It is known that GGA
generally underestimates the binding strengths whereas LDA overestimates them. To verify our results, we have
carried out calculations using the hybrid functional, B3LYP, since it has been successfully used to study the
hydrogen storage in alkali metal doped C$_{60}$. We find that the LDA energies are comparable to those obtained
by B3LYP, which are almost twice of the same computed by GGA.~\cite{nanolett-Gho08} The optimized structures
show that none of the relevant physical quantities, such as average H$_{2}$ binding energy, the Na-H bondlength
and the H-H distance undergoes significant changes as the number of adsorbed H$_{2}$ molecules increases.
However, there is monotonous decrease in the net charge on the sodium atom up to about 0.36$e$ for six H$_{2}$
molecules. We have observed that all the six H$_{2}$ molecules are absorbed at the optimum distance as shown
in Table 1. The optimum distance is an important parameter for estimating the hydrogen storage capacity of any
material. It determines a finite-scale region for storing more H$_{2}$ molecules around the doped atoms. So a
large distance is advantageous to reduce the repulsion between the adsorbed molecules. It may be noted that the
hydrogen molecules are found to interact with ionized Na atom as a dihydrogen molecule, forming a T-shaped complex,
and thus retaining their molecular identity irrespective of the type of SWCNTs.
\begin{table}
\small
\center
\begin{center}
\caption{Binding Energies (in~eV) of Na and H$_{2}$ molecules on (3,3), (4,2), (5,0), (5,5), (6,3) and (8,0)
SWCNT-Na systems}
\begin {tabular}{|c|c|c|c|c|c|c|c|c|c|c|}
\hline
No. of H$_2$ &
\multicolumn{2}{c|}{(3,3)SWCNT-Na}    &       \multicolumn{2}{c|}{(4,2)SWCNT-Na}      &
\multicolumn{2} {c|}{(5,0)SWCNT-Na}   &       \multicolumn{2} {c|}{(5,5)SWCNT-Na}     &
\multicolumn{2}{c|}{(6,3)SWCNT-Na}    \\      \hline
molecules&       LDA     &       GGA     &       LDA     &       GGA     &       LDA  &       GGA
        &  LDA  &       GGA     &       LDA     &       GGA     \\      \hline
0       &       -1.54   &       -1.04   &       -2.51   &       -1.71   &       -2.59 &
-1.69   &       -2.11   &       -1.66   &       -2.65   &       -1.83   \\
1       &       -0.22   &       -0.10   &       -0.24   &       -0.12   &       -0.20 &
-0.11   &       -0.30   &       -0.13   &       -0.25   &       -0.13   \\
2       &       -0.20   &       -0.10   &       -0.24   &       -0.13   &       -0.25 &
-0.12   &       -0.28   &       -0.14   &       -0.26   &       -0.13   \\
3       &       -0.20   &       -0.15   &       -0.25   &       -0.10   &       -0.24 &
-0.11   &       -0.28   &       -0.11   &       -0.26   & -0.11 \\
4       &       -0.21   &       -0.14   &       -0.25   &       -0.11   &       -0.25 &
-0.13   &       -0.28   &       -0.10   &       -0.27   &       -0.10   \\
5       &       -0.23   &       -0.11   &       -0.25   &       -0.10   &       -0.25 &
-0.10   &       -0.26   &       -0.10   &       -0.27   &       -0.11   \\
6       &       -0.22   &       -0.10   &       -0.25   &       -0.09   &       -0.23 &
-0.08   &       -0.25   &       -0.09   &       -0.26   & -0.10 \\      \hline
\end{tabular}
\end{center}
\end{table}

We have investigated the effects of nanotube chiral index (metallic vs semiconducting) of SWCNTs on H$_{2}$
adsorption. Accordingly, we consider (8,0) and (5,5) SWCNT which are of nearly same diameter. It is observed
that the calculated binding energies per H$_{2}$ molecule on Na-doped (8,0) SWCNT are very close to that for
(5,5) SWCNT for all $n$=1,...,6 coverage of H$_{2}$. This indicates that the binding of H$_{2}$ to the
SWCNT-Na complex has little to do with the chirality of nanotubes. It is, however, important to know if the
results reported above for (8,0) and (5,5)~SWCNT hold for other nanotubes as well and if they depend on the
chirality. Here, besides (8,0) and (5,5), we have studied four other different SWCNT-Na-(H$_{2}$)$_n$, $i.e.$,
corresponding to (3,3), (4,2), (5,0), and (6,3)~SWCNT. The optimized geometries of SWCNT-Na-(H$_{2}$)$_n$
show some regularity. We observe that the H$_{2}$ molecules get adsorbed on and around the Na atom site.
Thus, one H$_{2}$ molecule is always found to sit atop the Na atom, while the remaining five H$_{2}$ molecules
sit around the Na atom and the H-H bonds are oriented tangential to SWCNT surface which form T-shape complex
(see figure 4). This holds true for all the SWCNT-Na-(H$_2$)$_n$ systems. Table 2 shows the binding energy
data for Na and H$_{2}$ molecules. These results show that although the binding energy of Na is different on
metallic and semiconducting nanotubes, it doesn't affect the binding and the uptake capacity of H$_{2}$ molecules.
In conclusion, the phenomenon that a single Na atom adsorbed on hollow site of SWCNT can bind up to six
molecular hydrogen is  very general and holds for all the geometries of nanotubes considered here and for
other carbon-based nanomaterials such as C$_{60}$ as reported earlier.~\cite{nanolett-Gho08}  In this
context, it may however be noted that different types of fullerenes doped with alkali metal atoms do not
necessarily adsorb the same number of H$_{2}$ molecules, where the number varies from 3 to 6. In theses
cases, curvature seems to play a more significant role.
\begin{figure}
\center
\includegraphics[scale=0.2]{fig7.eps}
\caption{The coverage of the H$_{2}$ molecules on (3,3), (4,2), (5,0), (5,5), (6,3) and (8,0) SWCNT-Na systems.}
\label{fig-cnt-all-na-h2}
\end{figure}

Now, we discuss the stable Na adsorbates with maximum coverage on SWCNT. If the sodium atoms are doped at
all available hexagonal sites, it is expected that 144H$_{2}$ molecules can be adsorbed by the (8,0)
SWCNT-24Na system. However, it is important to note that the adsorption of the maximum number of H$_{2}$
molecules at one cationic site is only possible when the cationic sites are separated by $\sim$5.0~\AA,
at which the system can experience minimum electrostatic repulsion and steric hindrance. Accordingly,
we have identified such sites. In figure 8, we show six cases where Na covers 1/3 area of the hollow
sites of (3,3), (4,2), (5,0), (5,5), (8,0) SWCNTs and 1/4 hollow sites of (6,3) SWCNT. The figure also
shows the binding energy per Na and the Na-Na distances. The calculated binding energy values show that
there is strong binding between the metal atoms and the SWCNT systems. We first discuss the (8,0)~SWCNT for
maximum coverage of Na. When Na atoms are adsorbed on the surfaces of SWCNT, there is significant charge
transfer from the sodium atoms to the carbon atoms. As the coverage become higher (maximum eight Na atoms),
the electrostatic repulsion increases due to shortening of the Na-Na distances which makes the system unstable.
(Recall that there is a charge transfer from sodium atoms to SWCNT making the sodium positively charged.)
For stable Na adsorbates on carbon nanotubes, the Na-Na distance should be larger than that in bulk Na
which is 3.12~\AA.
\begin{figure}
\center
\includegraphics[scale=0.9]{fig8.eps}
\caption{The coverage of Na atoms on (3,3), (4,2), (5,0), (5,5), (6,3) and (8,0) SWCNT-Na systems.}
\end{figure}
For the stable Na adsorbates on carbon nanotubes, the Na-Na distance is 5.82~\AA~ with a binding energy
of -1.22~eV per Na (In case of single Na, the binding energy is -2.05~eV). The decrease in binding energy
can be understood from the net charge on sodium atoms. It is found that as the number of Na atoms increases,
the charge on each Na atoms decreases to $\sim$0.48$e$ leading to decrease in binding energy per Na atom.
The calculated binding energy per Na atom is however larger than the cohesive energy of -1.13~eV of bulk Na.
This shows that clustering of Na atoms is unlikely to occur on SWCNT surface. It may be noted that when the
Na-Na distance is less than 5.1~\AA, the net hydrogen molecule adsorption on Na-SWCNT decreases which is
consistent with previous reported result\cite{nanolett-Gho08}.
\begin{figure}
\center
\includegraphics[scale=0.22]{fig9.eps}\\
\includegraphics[scale=0.20]{fig10.eps}
\caption{Optimized geometries of the two, three, four and eight Na atoms isolated and aggregated on a
pristine (8,0)SWCNT. ${\Delta}$E represent the energy difference with reference to ground state.}
\end{figure}

We now address the important of issue of Na clustering which directly determines the reversibility of hydrogen
adsorption and desorption. To address this, we have considered the clusters of 2-8 Na atoms on the surface
of (8,0) SWCNT as shown in figure 9. To check whether Na atoms form a cluster on the (8,0) SWCNT surface,
we have carried out two separate calculations. First, we placed two Na atoms at the hollow sites of next
nearest neighbor carbon hexagon such that their separation distance is 7.3~\AA ~(column two of figure 9).
In the second case, the Na atoms are on two adjacent carbon hexagon with a Na-Na separation distance of
4.3~\AA~(column 3 of figure 9) for $n$=2 and decreasing to 3.3 \AA~as the cluster size grows to $n$=8.
The total energy difference ($\Delta$E) relative to isolated Na coverage (column 1 of figure 9) shows that
for n=2, the coating structures are more stable than those involving clustering, and are 0.14 and 0.11~eV
lower in energy. The Na clustering geometries are however uniformly more stable than the isolated Na systems,
for n~$\ge$~2 and the stability increases with $n$, from 0.34~eV for $n$=3 to 2.21~eV for n~=~8; whereas the
Na coating/clustering structures of column 2 of figure 9 are less preferred to the isolated Na geometries. The
binding energies, 0.27~eV/H$_{2}$ for Na$_{2}$ and 0.22~eV/H$_{2}$ for Na$_{4}$ complexes hardly changes from
the saturation value $\sim$0.26~eV/H$_{2}$ for 6H$_{2}$ complexes. The above analysis shows that for the
clustering of Na atoms to occur, the Na-Na distance must be $\ge$~ 5.9~\AA. It is further to be noted
that our calculated binding energy for Na (for Na$_{8}$) on (8,0) SWCNT (-1.22~eV) is larger than the
cohesive energy of bulk Na, namely, -1.13~eV. Since the energy gained from the cohesion is lower than that
from binding to SWCNT, clustering of Na atoms is unlikely to occur on SWCNT.
\begin{figure}
\center
\includegraphics[scale=0.24]{fig11.eps}
\caption{Optimized geometries of the twelve and twenty four H$_{2}$ molecules on the two and four aggregated
Na atoms on (8,0)SWCNT.}
\end{figure}

Now we discuss in brief the effect of Na clustering on the adsorption of the H$_{2}$ molecules. In Figure
10, we show the optimized geometries for twelve and twenty four H$_{2}$ molecules on the aggregated Na$_{2}$
and Na$_{4}$ complexes. The optimized C-Na separation ranges from 2.5~\AA~to 3.0~\AA~. The clustering of Na has
thus little effect on the binding energy per H$_{2}$ molecules, and the uptake capacity of H$_{2}$ adsorption.
\begin{figure}
\center
\includegraphics[scale=0.8]{fig12.eps}
\caption{Energy Vs reaction path for successive transformation of Na atom from one C-C ring to next nearest
neighbor over a (8,0) SWCNT. Only one transition state is found during the reaction.}
\end{figure}

To use Na-doped SWCNT as a hydrogen storage materials, we have studied the transition state search using
nudged elastic band method (NEB) to calculate the diffusion barrier which directly determines the adsorption
and desorption mechanism of metal and hydrogen molecules. As the interaction between H$_{2}$ molecules and
Na is physisorption in nature, we consider the transition state search for Na atom only. We calculate the
diffusion barrier of Na atom between the two nearest neighbor C-C rings of the SWCNT. Through optimized search of the
transition state, only one transition state is found which is shown in figure~11. This transition state has
the Na atom on the top of the bridge site of the two C-C rings. The diffusion barrier (0.34~eV) being larger than
0.26~eV/H$_{2}$, which is the binding energy, implies that during the desorption of H$_{2}$ molecules from the
metal-nanotube complex, Na atom should remain in its equilibrium state. It may further be noted that the Na
atom can transit from site to site easily on the carbon nanotube which facilitates the ultimate formation of
high coverage.
\begin{table}
\center
\begin{center}
\caption{The binding energy per H$_{2}$ molecules (in~eV) and hydrogen weight percentage of the different SWCNTs}
\begin{tabular}{|l|c|c|c|c|c|}
\hline
Systems &       \multicolumn{2}{c|}{H$_{2}$} &  \multicolumn{2}{c|}{Na} &       \%~wt~  \\
\hline
&       LDA     &       GGA     &       LDA     &       GGA     &       \\      \hline
(3,3)SWCNT-Na$_{6}$-(H$_{2}$)$_{36}$    &       -0.24   &       -0.11   &       -1.26   &       -0.93
&       11.28 \\
(4,2)SWCNT-Na$_{7}$-(H$_{2}$)$_{42}$    &       -0.23   &       -0.10   &       -1.48   &       -1.18
&       9.20    \\
(5,0)SWCNT-Na$_{5}$-(H$_{2}$)$_{30}$    &       -0.24   &       -0.13   &       -1.73   &       -1.34
&       9.22   \\
(5,5)SWCNT-Na$_{10}$-(H$_{2}$)$_{60}$       &       -0.21   &       -0.09   &       -1.05    &      -0.85
&       11.28 \\
(6,3)SWCNT-Na$_{11}$-(H$_{2}$)$_{66}$       &       -0.23   &       -0.11   &       -1.12    &      -0.94
&       9.50 \\
(8,0)SWCNT-Na$_{8}$-(H$_{2}$)$_{48}$        &       -0.24   &       -0.08   &       -1.22    &      -0.86
&       9.20      \\   \hline
\end{tabular}
\end{center}
\end{table}

We have so far discussed the interaction of H$_{2}$ with a single atom bonded to a nanotube, but clearly one
imagines attaching more Na to a nanotube, thereby increasing the hydrogen storage capacity. Figure 12 shows
the optimized structures of all geometries considered here. Table 3 shows the binding energies of Na, H$_{2}$
molecules and weight percentage. We have observed that the binding energy~-0.24~eV of (8,0)SWCNT-Na$_{8}$-(H$_{2}$)$_{48}$
is larger than -0.17~eV/H$_{2}$ in (H$_{2}$)$_{64}$-Li$_{8}$-C$_{64}$ system with gravimetric density of
13.45~wt\%\cite{jpc-Liu09}. In fact, these configurations store H$_2$ molecules approximately in the range
of 9.2--11.28~wt~$\%$. However, these numbers are based on the assumption that SWCNT-Na-H$_{2}$ group will also
release the hydrogen molecule without difficulty. This assumption is perhaps a good one, since the binding
energy of H$_{2}$ in SWCNT-Na-H$_{2}$ is about 1/8 of the binding energy of SWCNT-Na.
\begin{figure}
\center
\includegraphics[scale=0.2]{fig13.eps}
\caption[small]{Optimized atomic structure of maximal number of adsorbed H$_{2}$ molecules on (3,3), (4,2), (5,0),
(5,5), (6,3) and (8,0) SWCNT-Na$_{n}$ systems.}
\end{figure}

Thus far, for the various states of the chemisorptions process, we have obtained the energy and geometry
structure of atoms and molecules from first principle calculations at absolute temperature. A practical
question thus arises as to how pressure and temperature effects will influence these stable and transition
states, and how the associated molecular vibrations vary over these energy states at given pressure and
temperature. For an accurate understanding of the effects of pressure and temperature detailed thermodynamic
study is required. However, such a study is beyond the scope of present work. Hence, we give a qualitative
explanation of effect of pressure and temperature on the stability and storage capacity. It is noted that
under the influence of high pressure, the internuclear separation decreases. The decrease in internuclear
distances induces the confinement of electrons which result in localization and overlap of atomic orbitals,
so the binding strength and storage capacity increases. This is supported by recent report on pressure
dependence on hydrogen storage in Li-doped carbon nanotube\cite{catalysis-jung07}. The dependence of the
stability of composite system (Na$_{n}$-SWCNT-(H$_{2}$)$_{n}$ on temperature and pressure is also explained
with chemical potential which is given in terms of Gibbs free energy difference. The Gibbs free energy
difference for the formation of stable adsorbates is given as follows;
\begin{equation}
\Delta{G}=E_{ads-swcnt}-E_{swcnt}-x\mu_{ads}
\end{equation}
where {\it{x}}, $\mu_{ads}$ are the hydrogen coverage and chemical potential of adsorbates
respectively. The $\mu_{ads}$ depends on the temperature {\it{T}} and pressure {\it{P}} of the environment.
The stable states are obtained by determining the minima under the same chemical potential. For an ideal gas,
$\mu_{ads}$ is given as
\begin{equation}
\mu_{ads}=-k_{B}Tln\frac{(2\pi m)^\frac{3}{2}(kT)^\frac{5}{2}}{Ph^3}
\end{equation}
where $k_{B}$, m and {\it h} are Boltzmann constant, mass of adsorbate atoms, and Planck constant respectively.
From above equation, we found that as the temperature approaches zero, $\mu_{ads}$~$\rightarrow$~0 and adsorbates
coverage on the SWCNT's reaches a maximum. Here, we have $\mu_{ads}$$~<$~0 with increasing temperature, which
results in decrease in the stability and storage capacity of adsorbates\cite{prb-yang06}.
\section{Conclusions}
In summary, using the state of the art first-principles total energy calculations at absolute temperature and
pressure, we have shown that each Na atom adsorbed on SWCNT can bind up to six hydrogen molecules with an average
binding energy of -0.26~eV/H$_{2}$. It is important to note that hydrogen molecules are found to interact with Na
atom as a dihydrogen molecule, forming T-shaped complex and thus retaining their molecular identity irrespective
of the type of SWCNTs. Our calculations show that the electric field produced by SWCNT-Na polarizes the H$_{2}$
molecules and the resulting induced dipole-dipole as well as the charge-dipole interaction, consequently,
induces the binding of H$_{2}$ around SWCNT-Na. The PDOS of SWCNT-Na-H$_{2}$ shows that there is no bonding
orbital between H$_{2}$ and Na, which also supports this argument. We have also investigated the effects
of nanotube chiral index (metallic and semiconducting) on Na and H$_{2}$ adsorption and found that semiconducting
zigzag nanotubes, (8,0) and (5,0), show stronger binding energies of the Na atom, -2.1~eV and -2.6~eV
respectively, as compared to metallic SWCNTs. We have shown that the Na adsorbates on SWCNTs are highly stable
without the tendency for clustering at Na-Na distance of 5.9~\AA. Transition state calculations show that Na
atom can easily move on the carbon nanotube and fill the nanotube exterior with metal and retain its stability during
desorption of H$_{2}$ molecules. At full coverage this yields a storage capacity of $\sim$ 9.2-11.28~wt~$\%$,
which is higher than the value set for required gravimetric density of hydrogen storage. Our analysis shows that,
although indeed Na-coated SWCNTs presents a potential material for the hydrogen storage, utter care should be
taken to avoid Na atoms clustering on support material, to achieve higher hydrogen capacity.
\section{Acknowledgement}
This research is supported by Grant no. BRNS/BSC/37/38/3784 from the Department of Atomic Energy, Board of
Research in Nuclear Sciences (DAE-BRNS), India. It is also supported in part by the Department of Science and
Technology (DST) through computing resources provided by the High Performance Computing Facility
at the Inter University Accelerator Centre (IUAC), New Delhi. SC acknowledges DST for providing financial assistance.


\section*{Introduction}
In this paper we study the following  polynomial differential system in ${\mathbb{C}}^n$:
$$\dfrac{dX_1}{dt}=A_1(X_1,\ldots,X_n), \, \ldots, \, \dfrac{dX_n}{dt}=A_n(X_1,\ldots,X_n),$$
where $A_i \in {\mathbb{C}}[X_1,\ldots,X_n]$ and $\deg A_i \leq d$.  We associate  to this polynomial differential system the polynomial derivation $D=\sum_{i=1}^nA_i(X_1,\ldots,X_n)\partial_{X_i}$. \\

The computation of  first integrals of such  polynomial differential systems is an old and classical problem. The situation is the following: we want to compute a function $\mathcal{F}$ such that the hypersurfaces $\mathcal{F}(X_1,\ldots,X_n)=c$, where $c$ are constants, give  orbits of the differential system. Thus we want to find  a function $\mathcal{F}$ such that $D(\mathcal{F})=0$.\\

In 1878, G.~Darboux \cite{Darboux} has given a strategy to find first integrals. One of the tools developed by G.~Darboux is now called \emph{Darboux polynomials}. \\
A polynomial $f$ is said to be a \emph{Darboux polynomial},  if $D(f)=g.f$, where $g$ is a polynomial. The polynomial $g$ is called the cofactor. A lot of properties of a polynomial differential system are related to Darboux polynomials of the corresponding derivation $D$, see e.g.\cite{Goriely,Dumortier_Llibre_Artes}.\\
There exists a lot of different names in the literature for Darboux polynomials, for example we can find: special integrals, eigenpolynomials, algebraic invariant hypersurfaces, or special polynomials. \\

 G.~Darboux shows in \cite{Darboux2} that \emph{if the derivation  $D$ has at least $\binom{n+d-1}{n}+1$ Darboux polynomials then $D$ has a  first integral which can be expressed by means of these polynomials}. More precisely the first integral has the following form: $\prod_i f_i^{\lambda_i}$ where $f_i$ are Darboux polynomials and $\lambda_i$ are complex numbers. This kind of integral is called today a Darboux first integral.

In 1979,   J.-P.~Jouanolou shows in his book \cite{Joua_Pfaff}, that \emph{if a derivation has at least $\binom{n+d-1}{n}+n$ Darboux polynomials then the derivation has a rational first integral}. We recall that a rational first integral is a first integral which belongs to ${\mathbb{C}}(X_1,\ldots,X_n)$. \\
Several authors have given simplified proof for this result. M.~Singer proves this result in ${\mathbb{C}}^2$, see \cite{Singer}. This approach is based on a work of Rosenlicht \cite{Rosenlicht}. J.-A. Weil generalizes this strategy and gives a proof for derivations in ${\mathbb{C}}^n$, see \cite{Weil}. J. Llibre and X. Zhang gives a direct proof  of Jouanolou's result in \cite{LZ}. \\

Darboux and Jouanolou theorem are improved in \cite{LZ1,LZ2}. The authors show that we get the same kind of result if we take into account the multiplicity of Darboux polynomials. The multiplicity of a Darboux polynomial is defined and studied in \cite{Christopher_Llibre_Pereira}.\\

The Darboux theory of integrability has been successfully used in the study of some physical problems, see e.g. \cite{Valls,LlibreValls}, and in the study of  limit cycles and centers, see e.g. \cite{Christopher1994,Schlomiuk,Llibre-rodriguez}. Unfortunately, to the author knowledge, there do not exist example showing if these bounds are optimal. In this note we are going to study the situation in the sparse case. This means that we are going to consider  polynomials $A_i$ with some coefficients equals to zero. In this situation, the size of the polynomials $A_i$ is not measured by the degree but by the size of its Newton polytope. We recall that the Newton polytope of a Laurent polynomial $f(\underline{X})=\sum_{\alpha} c_{\alpha}X^{\alpha}$, where $\underline{X}=X_1,\dots,X_n$ and $\alpha$ is a multi-index $(\alpha_1,\ldots,\alpha_n) \in {\mathbb{Z}}^n$, is the convex hull in ${\mathbb{R}}^n$ of the exponent $\alpha$ of all  nonzero terms of $f$. We denote this polytope by $\mathcal{N}(f)$.\\

In this note we prove a result improving Darboux and Jouanolou theorem. Our result depends on the size of a Newton polytope associated to the derivation and not on the degree $d$. Furthermore, in this context we can give example showing that the bound is optimal. This is our result:
\begin{Thm}\label{Thm}
Let $D=\sum_{i=1}^nA_i(X_1,\ldots,X_n)\partial_{X_i}$ be a derivation.
Consider generic values $(x_1,\ldots,x_n)$ in ${\mathbb{C}}^n$ and the polytope $N_D=\mathcal{N}\Big(\sum_{i=1}^n x_i \dfrac{A_i}{X_i}\Big)$.\\
Let $B$ be the number of integer points in $N_D \cap {\mathbb{N}}^n$, then
\begin{enumerate}
\item \label{darbouxsparse} if $D$ has at least $B+1$ Darboux polynomials then $D$ has a Darboux first integral,
\item \label{Jouanolousparse} if $D$ has at least $B+n$ Darboux polynomials then $D$ has a rational first integral. Furthermore, this bound is optimal.
\end{enumerate}

\end{Thm}

We can remark that if we consider dense polynomials $A_i$ with degree $d$, that is to say each coefficient of $A_i$ is nonzero, then $B=\binom{n+d-1}{n}$. Thus Theorem \ref{Thm} gives the classical bounds in the dense case.\\

Now, we illustrate  why these bounds can be better than the classical ones.  We give an example with $n=2$ in order to give a picture. If $A_i$ has the following form: $A_i(X_1,X_2)=c_{e,e}X_1^eX_2^e+c_{e-1,e}X_1^{e-1}X_2^e+c_{e,e-1}X_1^{e}X_2^{e-1} +c_{0,0}$, then $B=3e+2$, and $d=\deg(A_i)=2e$. In this situation we have $\binom{n+d-1}{2}=2e(2e+1)/2$. Thus for such examples Theorem \ref{Thm} gives a linear bound instead of a quadratic bound.\\

Figure \ref{fig1} shows the Newton polygon of $A_i(X_1,X_2)$, when $e=3$. The triangle corresponds to the Newton polygon of dense polynomials with total degree $6$. In this situation, Jouanolou's theorem says that if we have $23$ Darboux polynomials then we have a rational first integral. Here, our bound improves this result and means that $13$ Darboux polynomials are sufficient to construct a rational first integral.

\begin{center}
\begin{figure}[!h]
\setlength{\unitlength}{.25cm}
\begin{picture}(10,10)
\put(0.5,1){$0$}
\put(1,2){\vector(1,0){7}}
\put(1,1.5){\vector(0,1){7}}
\put(-1,8){$X_2$}
\put(8.5,1.5){$X_1$}
\put(1,2){\circle*{.2}}
\put(2,3){\circle*{.2}}
\put(3,4){\circle*{.2}}
\put(4,4){\circle*{.2}}
\put(3,5){\circle*{.2}}
\put(3,2){\circle*{.1}}
\put(2,2){\circle*{.1}}
\put(4,2){\circle*{.1}}
\put(1,3){\circle*{.1}}
\put(1,4){\circle*{.1}}
\put(1,5){\circle*{.1}}
\put(1,6){\circle*{.1}}
\put(1,7){\circle*{.1}}
\put(2,4){\circle*{.1}}
\put(2,5){\circle*{.1}}
\put(2,6){\circle*{.1}}
\put(3,2){\circle*{.1}}
\put(3,3){\circle*{.1}}
\put(4,3){\circle*{.1}}
\put(5,3){\circle*{.1}}
\put(5,2){\circle*{.1}}
\put(6,2){\circle*{.1}}
\put(7,2){\circle*{.1}}
\put(6,3){\circle*{.1}}
\put(5,4){\circle*{.1}}
\put(4,5){\circle*{.2}}
\put(3,6){\circle*{.1}}
\put(2,7){\circle*{.1}}
\put(1,8){\circle*{.1}}
\thicklines\put(1,2){\line(3,2){3}}
\thicklines\put(1,2){\line(2,3){2}}
\thicklines\put(3,5){\line(1,0){1}}
\thicklines\put(4,5){\line(0,-1){1}}
\end{picture}
\caption{Newton polygon  $\mathcal{N}\big(A_i(X_1,X_2)\big)$.}
\label{fig1}
\end{figure}
\end{center}

The Newton polygon of $\dfrac{A_i(X_1,X_2)}{X_1}$ corresponds to a translation of the Newton polygon of  $A_i(X_1,X_2)$. In Figure \ref{fig2} we show $\mathcal{N}\Big(\dfrac{A_i(X_1,X_2)}{X_1}\Big)$.

\begin{center}
\begin{figure}[!h]
\setlength{\unitlength}{.25cm}
\begin{picture}(10,10)
\put(0.5,1){$0$}
\put(0,2){\vector(1,0){8}}
\put(1,1.5){\vector(0,1){7}}
\put(-1,8){$X_2$}
\put(8.5,1.5){$X_1$}
\put(1,2){\circle*{.1}}
\put(1,3){\circle*{.2}}
\put(2,4){\circle*{.2}}
\put(4,4){\circle*{.1}}
\put(3,5){\circle*{.2}}
\put(3,2){\circle*{.1}}
\put(2,2){\circle*{.1}}
\put(4,2){\circle*{.1}}
\put(1,3){\circle*{.1}}
\put(1,4){\circle*{.1}}
\put(1,5){\circle*{.1}}
\put(1,6){\circle*{.1}}
\put(1,7){\circle*{.1}}
\put(2,4){\circle*{.1}}
\put(2,5){\circle*{.2}}
\put(2,6){\circle*{.1}}
\put(3,2){\circle*{.1}}
\put(3,4){\circle*{.2}}
\put(4,3){\circle*{.1}}
\put(5,3){\circle*{.1}}
\put(5,2){\circle*{.1}}
\put(6,2){\circle*{.1}}
\put(7,2){\circle*{.1}}
\put(6,3){\circle*{.1}}
\put(5,4){\circle*{.1}}
\put(4,5){\circle*{.1}}
\put(3,6){\circle*{.1}}
\put(2,3){\circle*{.1}}
\put(3,3){\circle*{.1}}
\put(2,7){\circle*{.1}}
\put(1,8){\circle*{.1}}
\put(0,2){\circle*{.2}}
\put(2,5){\line(1,0){1}}
\put(3,5){\line(0,-1){1}}
\thicklines\put(0,2){\line(2,3){2}}
\thicklines\put(0,2){\line(3,2){3}}
\thicklines\put(2,5){\line(1,0){1}}
\thicklines\put(3,5){\line(0,-1){1}}
\end{picture}
\caption{Newton polygon  $\mathcal{N}\Big(\dfrac{A_i(X_1,X_2)}{X_1}\Big)$.}
\label{fig2}
\end{figure}
\end{center}
Figure \ref{fig3} shows the part of the Newton polygon of  $x_1\dfrac{A_1(X_1,X_2)}{X_1}+x_2\dfrac{A_2(X_1,X_2)}{X_2}$ in ${\mathbb{N}}^2$, when $x_1$, $x_2$ are generic.
\begin{center}
\begin{figure}[!h]
\setlength{\unitlength}{.4cm}
\begin{picture}(10,10)
\put(0.5,1){$0$}
\put(1,2){\vector(1,0){7}}
\put(1,1.5){\vector(0,1){7}}
\put(-0.5,8){$X_2$}
\put(8.5,1.5){$X_1$}
\put(1,2){\circle*{.2}}
\put(2,3){\circle*{.2}}
\put(3,4){\circle*{.2}}
\put(4,4){\circle*{.2}}
\put(3,5){\circle*{.2}}
\put(3,2){\circle*{.1}}
\put(2,2){\circle*{.2}}
\put(4,2){\circle*{.1}}
\put(1,3){\circle*{.2}}
\put(1,4){\circle*{.1}}
\put(1,5){\circle*{.1}}
\put(1,6){\circle*{.1}}
\put(1,7){\circle*{.1}}
\put(2,4){\circle*{.2}}
\put(2,5){\circle*{.2}}
\put(2,6){\circle*{.1}}
\put(3,2){\circle*{.1}}
\put(3,3){\circle*{.2}}
\put(4,3){\circle*{.2}}
\put(5,3){\circle*{.1}}
\put(5,2){\circle*{.1}}
\put(6,2){\circle*{.1}}
\put(7,2){\circle*{.1}}
\put(6,3){\circle*{.1}}
\put(5,4){\circle*{.1}}
\put(4,5){\circle*{.1}}
\put(3,6){\circle*{.1}}
\put(2,7){\circle*{.1}}
\put(1,8){\circle*{.1}}
\thicklines\put(3,5){\line(1,-1){1}}
\thicklines\put(2,2){\line(2,1){2}}
\thicklines\put(1,3){\line(1,2){1}}
\thicklines\put(2,5){\line(1,0){1}}
\thicklines\put(4,4){\line(0,-1){1}}
\thicklines\put(1,2){\line(1,0){1}}
\thicklines\put(1,2){\line(0,1){1}}

\end{picture}
\caption{Newton polygon $\mathcal{N}\Big(x_1\dfrac{A_1(X_1,X_2)}{X_1}+x_2\dfrac{A_2(X_1,X_2)}{X_2}\Big)\cap {\mathbb{N}}^2$.}
\label{fig3}
\end{figure}
\end{center}

\subsection*{Structure of the paper} In Section \ref{section:1} we give some results about Newton  polytopes and weighted degree. In Section \ref{section:proof} we prove Theorem \ref{Thm} and we show with an example that the bound is optimal.

 \section{ToolBox}\label{section:1}
In this section we introduce some notations and results that will be useful in Section \ref{section:proof}. These kinds of tools are already present in the work of Ostrowski, see \cite{Ostro}. For more results about  sparse polynomials, see e.g. \cite{GKZ,CLO}.

\begin{Def}
Let $P$ be a polytope, then $H$ is a supporting hyperplane of $P$ if
\begin{enumerate}
\item $H \cap P \neq \emptyset$,
\item $P$ is fully contained in one of the two halfspaces defined by $H$.
\end{enumerate}
\end{Def}

In our situation, as we consider Newton polytopes, the equation of $H$ has integer coefficients.  More precisely, the equation of $H$ is $\nu. m=a$, where $.$ denotes the usual scalar product,  $\nu$ is a vector with integer coefficients, and $a$ is an integer.\\

We can represent a Newton polytope with the equations of its supporting hyperplanes:\\
$$\mathcal{N}(f)=\{ m \in {\mathbb{Z}}^n \mid \nu_j.m\leq a_j, \textrm{ for } j=1,\ldots,k\},$$
where $k$ is the number of supporting hyperplanes.

Now we define  a degree related to a vector $\nu \in {\mathbb{Z}}^n$.
\begin{Def}
Let $\nu \in {\mathbb{Z}}^n$, we set $\deg_{\nu}(f)=\max_{m \in \mathcal{N}(f)} \nu.m$.
\end{Def}

Now, we explain why we can call $\deg_{\nu}(f)$ a degree.

\begin{Prop}\label{prop:degree}
Let $f$ and $g$ be two polynomials in ${\mathbb{C}}[X_1,\ldots,X_n]$, $\nu=(\nu_1,\ldots,\nu_n)$ in ${\mathbb{Z}}^n$ and $x_1,x_2$ two generic elements in ${\mathbb{C}}^2$.
\begin{enumerate}
\item $\deg_{\nu}(f+g)\leq \max\big(\deg_{\nu}(f),\deg_{\nu}(g)\big)$,
\item $\deg_{\nu}(x_1f+x_2g)= \max\big(\deg_{\nu}(f),\deg_{\nu}(g)\big)$,
\item $\deg_{\nu}(f.g)=\deg_{\nu}(f)+\deg_{\nu}(g)$,
\item $\deg_{\nu}(\partial_{X_i}f) \leq \deg_{\nu}(f)-\nu_i=\deg_{\nu}(f/X_i)$.
\end{enumerate}
\end{Prop}

\begin{proof}
\begin{enumerate}
\item As $\mathcal{N}(f+g)$ is included in the convex hull of  $\mathcal{N}(f) \cup \mathcal{N}(g)$, we have
$$\deg_{\nu}(f+g)=\max_{m \in \mathcal{N}(f+g)} \nu.m\leq \max_{m \in Conv( \mathcal{N}(f) \cup \mathcal{N}(g) )}\nu.m,$$
where $Conv(.)$ denotes the convex hull.\\
Furthemore, as we consider a convex set, we deduce that this maximum is reached on a point in $\mathcal{N}(f) \cup \mathcal{N}(g)$. Thus
\begin{eqnarray*}
\deg_{\nu}(f+g)&\leq &\max_{m \in \mathcal{N}(f) \cup \mathcal{N}(g) }\nu.m\\
&\leq&\max\big(   \max_{m \in \mathcal{N}(f)} \nu.m, \max_{m \in \mathcal{N}(g)} \nu.m \big)\\
&\leq&    \max\big(\deg_{\nu}(f),\deg_{\nu}(g)\big).
\end{eqnarray*}
\item As $x_1, x_2$ are generic then $\mathcal{N}(x_1f+x_2g)$ is equal to the convex hull of  $\mathcal{N}(f) \cup \mathcal{N}(g)$. Indeed, with generic $x_1$ and $x_2$ we avoid simplifications in the sum $x_1f+x_2g$. Then the proof in this case proceeds as before.
\item This result comes from the well-known result by Ostrowski, \cite{Ostro}, which gives: $\mathcal{N}(f.g)=\mathcal{N}(f)+\mathcal{N}(g)$, where $+$  in this situation is the Minkowski sum.
\item Let $e_i$ be the i-th vector of the canonical basis of ${\mathbb{R}}^n$, then we have
$$\max_{m \in \mathcal{N}(\partial_{X_i} f)} \nu.m \leq \max_{m \in \mathcal{N}(f)} \nu.(m-e_i)=\max_{m \in \mathcal{N}(f)} \nu.m -\nu_i =\deg_{\nu}(f)-\nu_i.$$
 We also have
$$\max_{m \in \mathcal{N}(f)} \nu.(m-e_i)=\max_{m \in \mathcal{N}(f/X_i)} \nu.m$$
this completes the proof.
\end{enumerate}
\end{proof}

Newton polytopes and degree $\deg_{\nu}$ are related by the following proposition.

\begin{Prop}\label{prop:inclusion}
Let $f \in  {\mathbb{C}}[X_1^{\pm 1},\ldots,X_n^{\pm 1}]$ be a Laurent polynomial with corresponding Newton polytope:
$$\mathcal{N}(f)=\{ m \in {\mathbb{Z}}^n \mid \nu_j.m\leq a_j, \textrm{ for } j=1,\ldots,k\},$$
where $\nu_j.m=a_j$, whith  $a_j \in{\mathbb{Z}}$, and $\nu_j \in {\mathbb{Z}}^n$ are the equations of the $k$  supporting hyperplanes of $\mathcal{N}(f)$.\\
Let $g \in {\mathbb{C}}[X_1^{\pm 1},\ldots,X_n^{\pm 1}]$ such that for $j=1,\ldots,k$, $\deg_{\nu_j}(g)\leq \deg_{\nu_j}(f)$ then $\mathcal{N}(g)$ is included in $\mathcal{N}(f)$.
\end{Prop}

\begin{proof}
We just have to remark that
$$\max_{n \in \mathcal{N}(g)}  \nu_j.n=\deg_{\nu_j}(g) \leq \deg_{\nu_j}(f) =\max_{m \in \mathcal{N}(f)} \nu_j.m = a_j.$$
Thus each element in $\mathcal{N}(g)$ satisfies the equations of $\mathcal{N}(f)$.
\end{proof}
\section{Proof of Theorem \ref{Thm}}\label{section:proof}
\subsection{Newton polytope and cofactors}
We are going to show that if $A_i$ are sparse then the cofactors are sparse. This property will be the main tool of the proof of Theorem \ref{Thm}.

\begin{Prop}\label{prop:cofsparse}
Let $D=\sum_{i=1}^nA_i(X_1,\ldots,X_n)\partial_{X_i}$ be a derivation. Let $f$ be a Darboux polynomial with corresponding cofactor $g$.\\
Consider generic values $(x_1,\ldots,x_n)$ in ${\mathbb{C}}^n$ and let $N_D$ be the convex set $\mathcal{N}\Big(\sum_{i=1}^n x_i \dfrac{A_i}{X_i}\Big)$
 then
$$\mathcal{N}(g) \subset  N_D \cap {\mathbb{N}}^n.$$
\end{Prop}

\begin{proof}
First, obviously $\mathcal{N}(g) \in {\mathbb{N}}^n$ since $g$ is a polynomial.\\
Second, we have $g.f=\sum_{i=1}^n A_i \partial_{X_i} f$ , thus for all $\nu \in {\mathbb{Z}}^n$ we have $\deg_{\nu}(g.f)=\deg_{\nu}(\sum_{i=1}^n A_i \partial_{X_i}f)$.  Thanks to Proposition \ref{prop:degree}, we deduce these inequalities
\begin{eqnarray*}
\deg_{\nu}(g)+\deg_{\nu}(f) & \leq & \max_i(\deg_{\nu}(A_i) \partial_{X_i}f)\\
& \leq & \max_i(\deg_{\nu}(A_i) +\deg_{\nu}\partial_{X_i}f)\\
& \leq & \max_i(\deg_{\nu}(A_i) +\deg_{\nu}(f) -\nu_i).
\end{eqnarray*}
Thus we have
\begin{eqnarray*}
\deg_{\nu}(g)&\leq& \max_i(\deg_{\nu} (A_i)  -\nu_i)\\
\deg_{\nu}(g)&\leq& \max_i\Big(\deg_{\nu} (\dfrac{A_i}{X_i})\Big)\\
\deg_{\nu}(g)&\leq& \deg_{\nu} \Big(\sum_{i=1}^n x_i.\dfrac{A_i}{X_i}\Big).\\
\end{eqnarray*}
Now, we apply Proposition \ref{prop:inclusion} with  $\nu\in {\mathbb{Z}}^n$ corresponding to supporting hyperplanes of $\mathcal{N}\Big(\sum_{i=1}^n x_i.\dfrac{A_i}{X_i}\Big)$ and we get the desired result.
\end{proof}

We can now prove easily Theorem \ref{Thm}. Indeed, we use the classical strategy to prove Darboux Theorem in our situation.\\
As for all cofactors $g_{f_i}$ associated to a Darboux polynomial $f_i$, we have, by Proposition \ref{prop:cofsparse},
$$\mathcal{N}(g_{f_i}) \subset  N_D \cap {\mathbb{N}}^n,$$
then all  cofactors belong to a ${\mathbb{C}}$-vector space of dimension  $B$, where  $B$ is  the number of integer points in $N_D \cap {\mathbb{N}}^n$. Thus if we have $B+1$ cofactors, then there exists a relation between them: $$(\star)\,\, \sum_{i \in I} \lambda_i g_{f_i}=0,$$
where $\lambda_i$ are complex numbers.
Now, we recall a fundamental and  straightforward result on Darboux polynomials: $g_{f_1.f_2}=g_{f_1}+g_{f_2}$.  Thus relation $(\star)$ gives the Darboux first integral $\prod_{i \in I}f_i^{\lambda_i}$. This proves the first part of Theorem \ref{Thm}.\\

Now, in order to prove the second part of our theorem, we can use the strategy proposed in \cite{LZ}. In \cite{LZ} the authors show that if we have $n$ relations of the type $(\star)$ then we can deduce a relation with integer coefficients, i.e. $\lambda_i \in {\mathbb{Z}}$. This gives a first integral of this kind: $\prod_{i \in I}f_i^{\lambda_i}$ with $\lambda_i \in {\mathbb{Z}}$, thus this first integral belongs to ${\mathbb{C}}(X_1,\ldots,X_n)$ and we have a rational first integral.\\
 In the sparse case as the cofactors belong to a ${\mathbb{C}}$-vector space of dimension $B$, if we have $B+n$ Darboux polynomials then we have $B+n$ cofactors and then we deduce $n$ relations between the cofactors. With the strategy used in \cite{LZ} we obtain that the derivation has a rational first integral.
\subsection{The bound is optimal}
Consider a polynomial $p(X_1) \in {\mathbb{C}}[X_1]$ with degree $d$. Let $\alpha$ be a root of $p$ and $\xi_2,\ldots,\xi_n$ be distinct complex numbers such that $p'(\alpha),\xi_2,\ldots,\xi_n$ are ${\mathbb{Z}}$-independent. We denote by $D$ the following derivation:
$$D=p(X_1)\partial_{X_1}+\xi_2X_2\partial_{X_2}+\cdots+\xi_nX_n\partial_{X_n}.$$

This derivation has no non-trivial rational first integrals, by Corollary 5.3 in \cite{Goriely}. Indeed, $(\alpha,0,\dots,0)$ is a fixed point of the polynomial vector field, and the corresponding eigenvalues are distinct and ${\mathbb{Z}}$ independent.\\

By Proposition \ref{prop:cofsparse}, if $g$ is a cofactor then $\mathcal{N}(g) \subset  N_D \cap {\mathbb{N}}^n$. Here, $N_D\cap {\mathbb{N}}^n$ is the set of univariate polynomials in $X_1$ with degree smaller than $d-1$. Thus $B=d$. \\

Furthermore, $(X_1-\alpha)$, $(X_1-\alpha_2),\ldots,(X_1-\alpha_d)$ where $\alpha_i$ are roots of $p$, and  $X_2,\ldots,X_n$ are Darboux polynomials. Thus we have $d+n-1$ Darboux polynomials.\\

In conclusion, we cannot improve the bound given in Theorem \ref{Thm}, since there exists a derivation without rational first integrals  which has  $B+n-1$ Darboux polynomials.