Layered palladates and their relation to nickelates and cuprates

Abstract: We explore the layered palladium oxides La$_2$PdO$_4$, LaPdO$_2$ and La$_4$Pd$_3$O$_8$ via ab initio calculations. La$_2$PdO$_4$, being low spin $d^8$, is quite different from its high spin nickel analog. Hypothetical LaPdO$_2$, despite its $d^9$ configuration, has a paramagnetic electronic structure very different from cuprates. On the other hand, the hypothetical trilayer compound La$_4$Pd$_3$O$_8$ ($d^{8.67}$) is more promising in that its paramagnetic electronic structure is very similar to that of overdoped cuprates. But even in the $d^9$ limit (achieved by partial substitution of La with a 4+ ion), we find that an antiferromagnetic insulating state cannot be stabilized due to the less correlated nature of Pd ions. Therefore, this material, if it could be synthesized, would provide an ideal platform for testing the validity of magnetic theories for high-temperature superconductivity.

\section{Introduction}

Cuprates are unique transition metal oxides in that the active transition metal $3d$ orbital ($d_{x^2-y^2}$) has a comparable energy to the $2p$ orbitals of the oxygen
ligands.  This leads to strong bonding-antibonding splitting, with a half-filled antibonding state characterizing its electronic structure \cite{pickett_cuprates}.
Chemical doping away from the stoichiometric $d^9$ configuration leads to high temperature superconductivity with a novel $d$-wave symmetry \cite{cuprates_r}.
From a materials design perspective, it is interesting to explore whether other transition metal oxides could have similar characteristics,
thus possibly leading to a new family of high temperature superconductors \cite{norman_RPP}.
In this context, layered nickelates have been an obvious target \cite{nickelate_analogs}.
The best known of the layered nickelates are the  Ruddlesden-Popper (RP) series, Ln$_{n+1}$Ni$_n$O$_{3n+1}$ (Ln being a lanthanide), composed of $n$-NiO$_2$ layers along the $c$-axis
separated by LnO spacer layers.  Their $d$ filling ranges from $d^8$ for the n=1 member to $d^7$ for the n=$\infty$ one. Single layer (n=1) Ln$_{2}$NiO$_4$
has been intensively studied \cite{214}. Upon Sr substitution for La, charge and spin
stripe order, a pseudogap phase, and
an insulator-to-metal transition have been reported,
but no superconductivity has been found.
Lower valences (higher $d$ occupations) can be achieved in RP phases by topotactic reduction.
This was used to create infinite-layer $d^9$ LaNiO$_2$ from the cubic perovskite LaNiO$_3$ \cite{hayward}.
Although isostructural to CaCuO$_2$, its behavior is very different from that of its cuprate counterpart given the reduced $d$-$p$ hybridization
and $d_{z^2}$ intermixing with La $d$ states \cite{lanio2_pickett}.
More promising has been the topotactic reduction of the n=3 member, Ln$_4$Ni$_3$O$_{10}$ (Ln= La, Pr).  Although La$_4$Ni$_3$O$_8$ ($d^{8.67}$)
is a charge ordered antiferromagnetic (AFM) insulator \cite{zhang},
Pr$_4$Ni$_3$O$_8$ is a metal that shares many of the characteristics of overdoped cuprates \cite{nickelates}.
However, electron doping this material has proven to be a challenge \cite{nickelate_doping}.

\begin{figure}
\includegraphics[width=0.9\columnwidth]{fig1.pdf}
\caption {Structure of (a) infinite layer (112), (b) single layer (214), and (c) trilayer (438) phases, with oxygen atoms in red, nickel/palladium atoms in blue, and lanthanum atoms in green. The square planar environment of the Ni/Pd ions is highlighted.}
\label{fig1}
\end{figure}

A different approach is to move down the periodic table and consider instead $4d$ analogs.  But none of the Ruddlesen-Popper phases of Ag exist.  So, we turn to Pd,
where the n=1 member (La$_2$PdO$_4$) \cite{la2pdo4} and the n=$\infty$ member (LaPdO$_3$) \cite{lapdo3} have been synthesized. Ln$_2$PdO$_4$ has been reported for a
variety of lanthanide ions \cite{la2pdo4} and has the T$^\prime$ structure found in Nd$_2$CuO$_4$ \cite{la2pdo4_2}.
They are insulators and show diamagnetic behavior consistent with a low spin (LS) $d^8$ configuration. Electron doping by substituting La$^{3+}$ by Ce$^{4+}$ has been reported in hopes
that metallic behavior would arise \cite{suzuki}. However, the highest doping concentrations reached (20\%) are not enough to induce metallicity,
though the resistivity is reduced.
No other Ruddlesen-Popper phases are known to exist.  This is likely because La$_4$PdO$_7$ and La$_2$Pd$_2$O$_5$ dominate the thermodynamic phase diagram \cite{jacob}.
But LaPdO$_3$ has been synthesized under pressure \cite{lapdo3}. Topotactic reduction would then give the possibility of $d^9$ Pd.  This would be significant, since
the only formally $d^9$ Pd compounds are the delafossites such as PdCoO$_2$ \cite{mackenzie}.
But there, the Pd actually forms metallic triangular sheets and exhibits free electron-like behavior.

Here, we use density functional theory (DFT) to study the electronic and magnetic properties of La$_2$PdO$_4$ (214 phase) as well as hypothetical LaPdO$_2$ (112 phase) and La$_4$Pd$_3$O$_8$ (438 phase),
and contrast them with their nickel analogs and with cuprates (Fig.~\ref{fig1}).  We find that the 214 and 112 phases are not promising as cuprate analogs.
However, the 438 phase has a remarkably similar paramagnetic (PM) band structure to that of overdoped cuprates.  On the other hand, we could not stabilize an antiferromagnetic insulating state in the $d^9$ limit (achieved by substituting one La$^{3+}$ per formula unit with a ${4+}$ ion). Hence, this material would be an ideal platform to test the importance of having a parent insulating phase with strong antiferromagnetic correlations for achieving high-T$_c$ superconductivity.

\section{Computational Methods}

Our electronic structure calculations were performed using the all-electron, full potential
code WIEN2k \cite{wien2k} based on the augmented plane wave
plus local orbitals (APW + lo) basis set.
Both LDA and LDA+$U$ \cite{sic} calculations were performed. For the latter, a Hund's rule $J$ is included with a typical value of 0.7 eV. The $U$ values for each calculation are specified below. For the structural optimizations, the Perdew-Burke-Ernzerhof version of the generalized
gradient approximation was used \cite{gga}.

For the calculations, we converged using R$_{mt}$K$_{max}$ = 7.0, a $k$
mesh of 14$\times$14$\times$14 for the 214 materials, 17$\times$17$\times$17 for the 438 materials, and 14$\times$14$\times$17 for the 112 ones. Muffin-tin radii of 2.38 a.u.~for La, 2.5 for Th,
1.97 a.u.~for Ni, 2.07 a.u.~for Pd, and 1.75 a.u.~for O were chosen.  Unless otherwise stated, experimental structural data were used.

\begin{figure}
\includegraphics[width=0.90\columnwidth]{fig2.pdf}
\caption {Atom-resolved density of states (DOS) and paramagnetic GGA band structures for La$_2$NiO$_4$ (upper panels) and La$_2$PdO$_4$ (lower panels).}
\label{fig2}
\end{figure}

\section{${\rm \bf La_2PdO_4}$: Single layer}

We start with single-layer La$_2$PdO$_4$ that displays a tetragonal $I_4/mmm$ structure with
PdO$_2$ layers in which the $d^8$ ions are in a square planar environment, similar to its nickel analog
(recognizing that La$_2$NiO$_4$ has the related K$_2$NiF$_4$ structure in which the Ni ions sit in an elongated octahedral environment of oxygens \cite{la2nio4_struct}).
Let us recall some basics about the paramagnetic GGA electronic structure
of La$_2$NiO$_4$ (Fig.~\ref{fig2}).  The La-$4f$ bands do not play a role in the vicinity of the Fermi level, appearing 2.5 eV above it. The O-$2p$ bands are between -7 to -3 eV.
The Ni-$3d$ bands extend from -2.5 to 2 eV. For the Ni ions, the crystal
field splitting in a distorted octahedral environment breaks the degeneracy of the
e$_g$ states, with the d$_{x^2-y^2}$ being higher in energy with a bandwidth of $\sim$2.5 eV crossing the Fermi energy,
and a narrower $d_{z^2}$ band below it.
The on-site energy of the O-$2p$ levels is shifted down in energy by 2.5 eV with respect to the Ni-$3d$ states,
in contrast to cuprates where they are nearly degenerate.

The paramagnetic GGA band structure of La$_2$PdO$_4$ is also shown in Fig.~\ref{fig2}.  The La-$4f$ bands are at higher energies. The O-$2p$ levels are shifted up in energy as well with respect to the nickel analog as can be seen in the DOS, with the hybridized O-$2p$ bands and Pd-$4d$ bands between -7 to 3 eV.
The much wider ($\sim$5 eV) $d_{x^2-y^2}$ band crossing the Fermi level reflects the stronger covalent bonding between the Pd-4$d_{x^2-y^2}$ and O-$2p$ orbitals
due to the larger radial extent of the $4d$ orbitals. For both Pd and Ni 214 phases, the involvement of $d_{z^2}$ bands gives rise to a Fermi surface different from that observed in cuprates (Fig.~\ref{fig3}).

\begin{figure}
\includegraphics[width=0.80\columnwidth]{fig3.pdf}
\caption {Fermi surfaces for Pd and Ni 214, 438 and 112 phases from paramagnetic GGA calculations.}
\label{fig3}
\end{figure}

Checkerboard antiferromagnetic (AFM) ordering even at the GGA level is sufficient to open up a gap ($\sim$ 0.8 eV) in La$_2$NiO$_4$ (Fig.~\ref{fig4}).
This AFM state is lower in energy by 150 meV/Ni than the paramagnetic state. The significant Hund's rule
coupling stabilizes a high spin (HS) state for Ni$^{2+}$ (S = 1), with a moment of 1.31 $\mu_B$ per Ni.
The application of a Coulomb $U$ simply increases the value of the gap.

\begin{figure}
\begin{center}
\includegraphics[width=\columnwidth]{fig4.pdf}
\caption {Band structure of La$_2$NiO$_4$ (left) in an AFM state within GGA and La$_2$PdO$_4$ (middle) within LDA+$U$. Here, a $U$ of 8 eV is used to obtain an insulating state for La$_2$PdO$_4$ (though it is not magnetic).
Energy level schematic for $d^8$ high spin (HS) Ni (lower right) and low spin (LS) Pd (upper right).}
\label{fig4}
\end{center}
\end{figure}

In La$_2$PdO$_4$ the situation is different. The much larger crystal field	splitting stabilizes a low spin state for Pd$^{2+}$ (S=0). This should lead to a crystal field gap between occupied d$_{z^2}$ and unoccupied d$_{x^2-y^2}$ states. But as described above, the large bandwidth closes this gap unless a Coulomb $U$ that exceeds 6 eV is applied (Fig.~\ref{fig4}).  Even then, no magnetism is found.
This value of $U$ is larger than the 3.76 eV value determined by constrained RPA calculations for Pd metal \cite{rpa}, and also larger than the 4.5 eV value estimated
from spectroscopic data on PdO \cite{uozumi}. A similar observation about $U$ has been reported for PdO \cite{bennett}, and may be connected with issues associated
with using LDA+$U$ for 4$d$ ions (it has been claimed that hybrid functionals work better for PdO \cite{derzsi} but a very large fraction of exact exchange has to be introduced).
It should be noted that the susceptibility does not give any evidence for AFM order,
simply showing T independent diamagnetic behavior, consistent with a low spin state for Pd, with a low T upturn likely due to impurities.
Regardless, it is clear that La$_2$PdO$_4$ is very different from La$_2$NiO$_4$ due to its much larger $d_{x^2-y^2}$ bandwidth and its low spin nature.

\section {${\rm \bf LaPdO_2}$: Infinite layer}

We have also studied the infinite layer hypothetical compound LaPdO$_2$. It was modeled in analogy with LaNiO$_2$ (isostructural with CaCuO$_2$) with
an assumed space group $P4/mmm$, in which both volume and $c/a$ were optimized giving rise to lattice parameters $a$=4.15 \AA~and $c$=3.47 \AA. Here, the Ni/Pd$^{1+}$ $d^9$ ions
have a square planar coordination.

\begin{figure}
\includegraphics[width=0.9\columnwidth]{fig5.pdf}
\caption {Paramagnetic GGA band structure and atom-resolved DOS of LaPdO$_2$ (lower panels) and LaNiO$_2$ (upper panels).}
\label{fig5}
\end{figure}

We first present the GGA results (Fig.~\ref{fig5}).
For LaNiO$_2$, the La-$4f$ bands are located at 2.5 eV. The O-$2p$ bands
extend from about -8 eV to -3 eV. The Ni-$3d$ bands are distributed
from -2.5 to 2 eV. Unlike in CaCuO$_2$, there are two bands
crossing the Fermi level. One is $d_{x^2-y^2}$ in character, the other one is
a mixture of La-5$d_{z^2}$ and Ni-3$d_{z^2}$. In addition, as in La$_2$NiO$_4$, there is a smaller degree of $d-p$ hybridization as compared to cuprates.
As has been noted before, this $d^9$ metallic compound is very different from its insulating cuprate analogs \cite{lanio2_pickett}.

Is the Pd case different? In Fig.~\ref{fig5}, the band structure and DOS for LaPdO$_2$ are also shown. The La-$4f$ bands are shifted to lower energies relative to the nickelate, and hybridize with the much wider Pd-$4d_{x^2-y^2}$ band (with a bandwidth of 5 eV).  However, the electronic structure is similar in that two bands
still cross the Fermi level: the wide Pd-4$d_{x^2-y^2}$ and the La 5$d_{z^2}$-Pd 4d$_{z^2}$ one. As in the 214 phases, the extra involvement of a $d_z^2$ band makes the Fermi surface considerably different to that of the corresponding cuprate counterpart (Fig.~\ref{fig3}).  Moreover, strong three dimensional behavior is evident
in the Fermi surface, which we find to be the case as well for the nickel analog (Fig.~\ref{fig3}), again unlike the cuprates.

A stable checkerboard AFM metallic state can be obtained for LaNiO$_2$ within GGA with a spin moment of 0.73 $\mu_B$ per Ni. This state has a lower energy by 15 meV/Ni
than that of the PM state even though there is no experimental evidence for magnetic order in this material.
In LaPdO$_2$, an AFM metallic state can also be stabilized. The derived moment on Pd$^{1+}$ within GGA is lower than in the nickelate case (0.22 $\mu_B$). We find that the PM state is more stable than the AFM one by 4 meV/Pd.

Given that substitution of Ni by Pd was not promising in terms of eliminating the $d_{z^2}$ contribution around E$_F$, we turn back to LaNiO$_2$ where one could attempt to substitute La by another 3+ ion to eliminate the
problematic La-5d$_{z^2}$ states from the vicinity of the Fermi energy. One approach might be to substitute La by the smaller Y.
The $P4/mmm$ structure was assumed once again with the lattice constants being optimized.
The band structure of hypothetical YNiO$_2$ is quite similar to that of its La counterpart with Y-4$d_{z^2}$ states still crossing E$_F$, giving rise to even larger Fermi surface pockets than in LaNiO$_2$.  We also tried using Tl (with its filled d shell) in place of Y, but two bands still cross E$_F$.

\section{${\rm \bf La_4Pd_3O_8}$: Trilayer}

\begin{figure}
\includegraphics[width=0.92\columnwidth]{fig6.pdf}
\caption {Paramagnetic GGA band structure and atom-resolved DOS of La$_4$Pd$_3$O$_8$ (lower panels) and La$_4$Ni$_3$O$_8$ (upper panels).}
\label{fig6}
\end{figure}

After ruling out the 214 and 112 layered palladates as cuprate analogs, we turn our attention to the 438 phases.
Ln$_4$Ni$_3$O$_8$ (Ln= La, Pr) has a tetragonal $I_4/mmm$ structure in which the Ni is in a square planar environment.
Although the La438 phase is known to be a charge-ordered insulator \cite{zhang}, Pr438 is metallic \cite{nickelates}.
Therefore, we focus on the paramagnetic electronic structure of La438 (Fig.~\ref{fig6}).
Again, the La-$4f$ bands are located more than 3 eV above the Fermi level. The O-$2p$ bands
extend from about -5 eV to 2 eV. The Ni-$3d$ bands are distributed
from -2 to 2 eV.  A single band per Ni with d$_{x^2-y^2}$ character and a 3 eV bandwidth due to large $d$-$p$ hybridization crosses the Fermi level.

Although higher-order RP phases are not known for the palladates, it is possible that they could be stabilized under high oxygen partial pressures.
If the n=3 phase could be stabilized, then the resulting 438 phase could conceivably be made by reduction.  To explore this possibility,
the Pd analog was modeled from the Ni 438 structure by performing an optimization of both volume and $c/a$ giving rise to lattice parameters $a$=4.17 \AA~and $c$=26.31 \AA.
We find that the the electronic structure of the palladate is very similar to that of the nickelate (Fig.~\ref{fig6}), though with a larger bandwidth (4.4 eV) due to
increased $d$-$p$ hybridization. We also find that the Fermi surfaces of Pd and Ni-438 phases are very similar and cuprate-like (Fig.~\ref{fig3}).
Note the presence of two large hole-like cylinders and one electron-like cylinder, indicating that the doping level of the inner plane differs from that
of the two outer planes (Fig.~\ref{fig1}), as has been inferred in trilayer (and higher-layer) cuprates by NMR \cite{nmr}.
Given the similarity of the Fermi surfaces, one might then expect that an electron-doped version of this layered palladate might be a promising candidate
for high temperature superconductivity \cite{eriksson}.
This is particularly relevant, given that substituting La by a smaller 4+ ion in the nickel analog has proven to be difficult; the larger size of Pd versus Ni may make this easier for the Pd case.

To explore this further, we investigated whether a magnetic state could be obtained.  Previously, we had been able to stabilize magnetic solutions for
La$_4$Ni$_3$O$_8$ and Pr$_4$Ni$_3$O$_8$ if we include an on-site Coulomb $U$ \cite{nickelates,antia1,antia2}. In the 438 phases, the average Ni/Pd valence is +1.33 which leaves the e$_g$ orbitals with 2.67 electrons per Ni/Pd on average. If the Hund's rule coupling is larger than the splitting between the two $e_g$ orbitals, the transition metal ion would be in a HS state, in the opposite situation, in a LS state.
If the latter occurs, all the $z^2$ bands (majority and minority spin) and 2/3 of a majority spin $x^2-y^2$ band will be occupied. Due to this partially filled band crossing the Fermi level, the 438 materials in a LS state are metallic, with all unoccupied states having $x^2-y^2$ parentage. The resulting magnetic order in this case is ferromagnetic (FM). This is the magnetic ground state we find in Pr$_4$Ni$_3$O$_8$ and also in hypothetical La$_4$Pd$_3$O$_8$. The electronic structure of the majority spin channel in this FM-LS state is analogous to that of the PM state shown in Fig.~\ref{fig6}. As in the 112 case, the derived magnetic moments are smaller for Pd (0.50 $\mu_B$) than for Ni (0.70 $\mu_B$) at the same $U$ value (3.5 eV).
For the palladate, given the large crystal field splitting, a HS state could not be stabilized (unlike for the nickelate).

To investigate further, we explored the effects of electron doping the 438 palladate.
For ThLa$_3$Pd$_3$O$_8$, Pd should have a $d^9$ configuration, the same as in parent cuprates. We chose Th over Ce to avoid dealing with $f$ states and because Ce has two plausible oxidation states (3+ and 4+). Here, an AFM state could be stabilized within LDA+$U$, but a gap cannot be opened regardless of the $U$ value due to the much stronger hybridization of Pd-$d$ and Th-$d$ states as compared to the nickelate \cite{antia2}. Though this metallic AFM state is more stable than a FM one by 27 meV/Pd  for $U$=3.5 eV, this energy difference is one order of magnitude smaller than in the nickelate.  The derived magnetic moment per Pd is 0.50 $\mu_B$ for this $U$ value.
Regardless, this AFM metallic solution differs from the AFM insulating solution found in the $3d$ nickelates and cuprates.

La$_4$Pd$_3$O$_8$ is then a cuprate analog in that it has a quasi-2D crystal structure, with strong $p$-$d$ hybridization, a single band of $d_{x^2-y^2}$ character per Pd crossing the Fermi level, and a similar Fermi surface topology \cite{eriksson,fs_wep}. Yet, the material is different in that a parent insulating antiferromagnetic phase is predicted to be absent. Hence, this material would be a good platform to investigate the importance of having an insulating phase with strong antiferromagnetic correlations for high-T$_c$ superconductivity, which would test whether superexchange is playing a key role as advocated by a number of authors \cite{pwa,lee,scalapino}.

\section{summary}

To summarize, the substitution of Ni with Pd in 214-single and 112-infinite layer compounds (aside from their formal structural similarities) induces drastic changes in the
electronic structure and magnetic properties, mainly due to the much larger crystal field splitting and bandwidth of the Pd compounds.
Therefore, we argue that the route to search for new high temperature superconductors by doping these layered palladates is not a promising one. However, the
hypothetical 438 material shares many cuprate-like features in terms of its paramagnetic electronic structure and fermiology, yet it is different in that a parent insulating antiferromagnetic phase could not be stabilized.
Therefore, if this material could be synthesized, it would be an ideal platform to test the importance of antiferromagnetic correlations for high temperature superconductivity.

\section{acknowledgments}

This work was supported by the Materials Sciences and Engineering
Division, Basic Energy Sciences, Office of Science, US DOE.
We acknowledge the computing resources provided on Blues, a high-performance computing clusters operated by the Laboratory
Computing Resource Center at Argonne National Laboratory.


Piezoresistance characterization of silicon nanowires in uniaxial and isostatic pressure variation

Abstract: Silicon nanowires (SiNWs) are known to exhibit large piezoresistance (PZR) effect, making it suitable for various sensing applications. Here, we report the results of a PZR investigation on randomly distributed and interconnected vertical silicon nanowire arrays as a pressure sensor. The samples were produced from p-type (100) Si wafers using a silver catalysed top-down etching process. The piezoresistance response of these SiNW arrays was analysed by measuring their I-V characteristics under applied uniaxial as well as isostatic pressure. The interconnected SiNWs exhibit increased mechanical stability in comparison with separated or periodic nanowires. The repeatability of the fabrication process and statistical distribution of measurements were also tested on several samples from different batches. A sensing resolution down to roughly \SI{1}{\milli\bar} pressure was observed with uniaxial force application, and more than two orders of magnitude resistance variation was determined for isostatic pressure below atmospheric pressure.

\section{Introduction}

Low dimensional structures may possess unique mechanical, electrical, optical, and thermoelectric properties.
Particularly, silicon nanowires (SiNWs), have demonstrated properties suitable for various advanced applications,
\cite{Peng2013,Heris2020,Heris2022}, including also low-cost thermoelectric devices and chemo-biological
sensors with ultrahigh sensitivity \cite{Zhou2003,Peng2009}. The integration of SiNWs in electronic devices is favored
by their compatibility with the well established Si-SiO\textsubscript{2} electronic industrial technology. Bulk silicon
has been known for a while to exhibit high piezo resistance (PZR) effect \cite{Smith1954}. In bulk semiconductors the
PZR-effect takes place, in principle, due to a change in the electronic structure and modification of the
charge-carriers effective masses. This phenomenon has found practical applications in many Si-based devices, such as
pressure transducers \cite{Tufte1962}, cantilevers for atomic force microscopy \cite{Tortonese1993}, accelerometers
\cite{Ning1995}, biosensors \cite{Wee2005}, and multi-axis force sensing tools \cite{Tiwari2021}.

Recently, nanowires have been shown to possess the ability to significantly increase the PZR response. A giant PZR was observed in p-doped SiNWs with diameters of \SIrange{50}{350}{\nano\meter} and length of microns, initially under tensile uniaxial stress \cite{He2006}. However, the PZR effect in n-doped nanowires was found comparable to that in the bulk counterpart, both for tensile and compressive uniaxial stress  \cite{Gao2017}.

On the theoretical side, the origin of the PZR effect in SiNWs has long been under debate and most frequently it is
referred to as anomalous PZR \cite{Zhang2012}. It has been related to quantum confinement effects \cite{Cheng2018},
surface charge effects \cite{Nguyen2021,Kim2020,Ghosh2021}, strain-induced bandgap shift \cite{shiri2008}, or changes of
the charge carriers effective masses \cite{Zhang2016}. A complex model incorporating these mechanisms has been proposed
in order to analytically quantify the PZR effect in silicon \cite{Rowe2014}.

A survey of several PZR sensors based on SiNWs, e.g. cantilever \cite{Toriyama2002}, opto-mechanical sensor
\cite{Toriyama2003}, flexible pressure sensor \cite{Kim2020}, or breath detector \cite{Ghosh2021}, shows that different
methods have been used for fabricating the SiNWs, such as vapor-liquid-solid (VLS), laser ablation, and metal assisted
catalysed etching (MACE) \cite{Schmidt2009}. Among these methods, MACE is the simplest and most versatile one. It relies
on catalyzed etching with assistance of a perforated metal template film (typically gold or silver)
\cite{svavarsson2016} or randomly distributed metallic nanoparticles (typically gold or silver)
\cite{Fakhri2021,plugaru2022structure} spread  on the Si-wafer. To date, studies have been focused on uniaxial load and
have neglected the SiNWs response under isostatic pressure, which creates a load uniformly distributed on the
sample surface. Here, we report on the PZR effect in SiNWs obtained by MACE, under uniaxial compression load, as well
as isostatic pressure in a vacuum chamber. We find that the interconnected SiNWs are mechanically stronger and
functionally more stable comparing with the arrays of separated wires under applied uniaxial pressure. They show
higher PZR sensitivity under isostatic pressure variation. We also demonstrate a simple, low cost and reproducible
fabrication method, to fabricate a robust and sensitive pressure sensor.

\section{Materials and Methods}

\subsection{Fabrication of SiNWs}

Arrays of interconnected SiNWs were fabricated by silver (Ag) MACE in a three step process, from p-type,
single-side polished \SI{525}{\micro\meter} thick Si wafers, with resistivity, $\rho$, of
\SIrange{10}{20}{\ohm\centi\meter}. The nanowire patterns were made on areas of about \SI{1}{\centi\meter\squared}, on
the polished side of the wafers. The sequence of steps used to prepare the SiNW areas is as follows:

\begin{enumerate}
\item Deposition of metal catalyst: Ag nano-particles were deposited on the surface of the Si-wafers by immersing the
wafers in a solution of \SI{3}{M} HF and \SI{1.5}{M} AgNO\textsubscript{3} for \SI{60}{\second}.

\item Wire etching: The samples were etched by immersing them in HF:H\textsubscript{2}O\textsubscript{2} (5M:0.4M) solution to obtain vertically aligned SiNWs.

\item Removal of residual Ag nano-particles: Samples were immersed in \si{20}{\%} w/v HNO\textsubscript{3} to remove residual silver particles. More detailed description can be found in Refs.~\cite{Fakhri2021,plugaru2022structure}.
\end{enumerate}

Subsequently, \SI{150}{\nano\meter} thick aluminum electrodes were deposited on the samples by electron beam evaporation (Polyteknik Cryofox Explorer 600 LT). For the uniaxial measurements, the electrodes were deposited on top and backside of the samples, while for the isostatic pressure measurements, the electrodes were made co-planar.

Four sets of interconnected SiNWs samples, denoted as follows, were prepared by varying the etching time: A
(\SI{1}{\minute}), B (\SI{3}{\minute}), C (\SI{5}{\minute}), D (\SI{7}{\minute}) were made for uniaxial pressure
application and sample E (\SI{40}{\minute}) was made for isostatic pressure testing.

Additionally, periodic SiNWs, sample F, was made as described in the Appendix.
\begin{table}[H]
\caption{List of samples with different etching time and corresponding length. Samples A,B,C,D,E are random wires and sample F is with periodic wires (shown in the Appendix).}
\label{tb:tab1}
\begin{center}
\begin{tabular}{ c c c }
\hline
\textbf{Sample~~~} & \textbf{Etching time (min)~~}	& \textbf{SiNW length (\textmu m)}	\\
\hline
A & 1 & 0.7		   \\
B & 3 & 1.5		   \\
C & 5 & 2.2        \\
D & 7 & 3          \\
E & 40 & 10    \\
F & 5 & 0.65       \\
\hline
\end{tabular}
\end{center}
\end{table}

A scanning electron microscope (SEM, Zeiss Supra 35) was used to characterize the SiNWs geometry. Top and cross-sectional SEM images were used to estimate the average diameter of the wires.
The Gwyddion software for data visualization and analysis was applied to surface SEM images in order to estimate the total top-area of the wires.
\autoref{fig:SEM} shows (a) the cross-section of wire array, and (b) the surface area of SiNWs obtained after \SI{7}{\minute} etching. It is worth mentioning here that the SEM analysis was carried out on the sample prior to removal of Ag-nanoparticles which are visualized as bright spots at the base of the nanowires in the cross-section.
\begin{figure
    \centering
    \includegraphics[width=0.8\linewidth]{Fig/finalfinalsem.png}
    \caption{(a) Cross-sectional SEM image of SiNWs, (b) top-view SEM image of SiNWs etched for \SI{7}{\minute}}
    \label{fig:SEM}
\end{figure}
\begin{figure
  \centering
  \includegraphics[width=0.8\linewidth]{100mic1r.png}(a)
 \includegraphics[width=0.72\linewidth]{Fig/gwe-eps-converted-to.pdf}(b)\hspace{17mm}
  \caption{(a) Gwyddion analysis of top-view SEM image in (\SI{100}{\micro\meter\squared}). (b) Two Gwyddion line-scans of the figure in left  (red color - Horizontal line,  black color - Vertical line).}
   \label{fig:Gwyddion}
\end{figure}
As can be seen in the top-view image, the wires are partly interconnected, forming a continuous rigid structure. Also seen from  place to place are free standing nanowires forming bundles.
Such bundle formation may take place because of capillary forces acting during the drying process following the wet-etching step.
In the cross sectional image one may observe that the length of the wires is relatively homogeneous, around
\SI{3}{\micro\meter} and their typical diameter is approximately \SI{150}{\nano\meter}.
According to Peng et al. \cite{Peng2009}, porosity plays an essential role in the PZR response, that increases with increased porosity.
The porosity is most conveniently controlled by the concentration of Ag deposition solution and etching time.
In a previous study \cite{Fakhri2021}, it was demonstrated that the SiNWs porosity was highly affected by the concentration of the AgNO\textsubscript{3} during the Ag-deposition.
A \SI{1.5}{mM} AgNO\textsubscript{3} (as used for samples A-D) provided highly porous SiNWs with maximum photoluminescence spectra intensity.

From the SEM image analysis, the wires cross-sectional surface coverage was estimated as roughly \si{28}{\%} by using
the Gwyddion program, \autoref{fig:Gwyddion}.
By counting the average number of wires on several line-scans, the wire density was estimated to be \SI{1.6 e7}{\milli\meter\squared}, which corresponds to roughly \num{8 e8} wires under the force meter area (\SI{7 x 7}{\milli\meter}).\\

\subsection{Measurement setups}

For uniaxial PZR tests, the samples were clamped between a metal pin of a force meter (Mark-10, M5-012), touching the top-side of the samples, and a rigid copper (Cu) plate at the backside as shown in \autoref{fig:schem}. In order to improve the electrical contacts, the samples were glued to the Cu backside plate with silver paste. The force meter was mounted on a vertical rod allowing for a movement on the z-axis (vertical) by a manually adjusted screw. A uniaxial pressure was applied to the samples by pressing the pin of the force meter into the samples' surface with intensity determined by the aforementioned screw.
The applied force was in the range of \SIrange{100}{900}{\milli\newton} on an area of \SI{7 x 7}{\milli\meter} area (the cross sectional area of the pin). Taking the wire coverage (\SI{28}{\%}) into account, this corresponding to a gauge pressure of \SIrange{7}{66}{\kilo\pascal}, respectively.
\autoref{fig:schem} shows sketch of the experimental setup.
A Keithley 2400 SourceMeter was used to measure the resistance $R$ through the sample as a function of applied pressure, at a constant bias voltage of \SI{2}{\volt}.

For PZR tests under isostatic pressure variation a vacuum chamber was used, and the air was removed while the resistance was measured at fixed \SI{5}{\volt}.  In this case (for sample E) the contacts were made in co-planar configuration, with separate contacts on each side of the patterned area, such that the nanowires are maximally exposed to air, instead of top-and-bottom contacts as for the uniaxial measurements (performed with samples A-D).
Furthermore, the wire length was increased (to \SI{10}{\micro\meter}) to increase the surface area exposed to the air.
The samples were mounted on a fixed sample station inside the chamber with tungsten needle tip kelvin probes as a connection.

\begin{figure
    \centering
   \includegraphics[width=0.6\linewidth]{Fig/schem2.png}
    \caption{Schematic of the force meter setup for testing PZR characteristics.}
    \label{fig:schem}
\end{figure}

\section{Results}

\subsection{Electrical response under uniaxial force}

A maximum force of approximately \SI{900}{\milli\newton} was applied to the samples using the setup shown in \autoref{fig:schem}.
The maximum  vertical  force exerted on each wire thus corresponds to \SI{900e-3}{\newton}/\num{8e8} \SI{\approx 1.1e-9}{\newton/\wire} or \SI{6.6e 4}{\pascal} (\SI{0.65}{\bar}).

The force and the instantaneous resistance measured after applying the force on the sample with nanowires of length \SI{3}{\micro\meter} (etched for \SI{7}{\minute}),  are shown in \autoref{fig:Force}(a), as function of time. The maximum force was applied at the beginning (time zero) and then reduced step-wise while the resistance was being measured. Each force level was kept constant for roughly 50~s to confirm that the resistance value was stable with time. As clearly visible, the resistance changes significantly and inversely with pressure.
\begin{figure
  \centering
    \includegraphics[width=0.8\linewidth]{Fig/shape1-1-eps-converted-to.pdf}
   \includegraphics[width=0.85\linewidth]{allcurves3-3-eps-converted-to.pdf}
   \caption{(a) Instantaneous resistance (right y-axis, red line) in response to applied force (left y-axis, blue line) over time for sample D.  (b) Relative resistance change versus normalized force for SiNWs samples A (purple), B (green), C (red), D (blue).}
   \label{fig:Force}
\end{figure}

We define the relative resistance change vs force variation, for each step $k=1,2,...$ as
\begin{equation}
   \frac{\Delta R}{R}=\frac{R_k-R_1}{R_1},
   \label{eq:first}
\end{equation}
where $R_k$ is the average resistance obtained for the constant force step $k$.
The measurements were repeated four times on each sample (A, B, C, D), every time reproducing the force steps as well as possible with the adjustable screw.
A second average, this time of the relative resistance in step $k$, defined by \autoref{eq:first}, obtained for the four measurements, vs. the average normalized forces for each step, are presented in \autoref{fig:Force}(b).
Here, by normalized force we mean the ratio of the force in a particular step $k$ to the initial force, $F_k/F_1$.  Hence, the maximum force values for each step-wise measurement are normalized to unity.
Note that the resistance decreases when the force increases, and that in \autoref{fig:Force}(b) we show the negative of $\Delta R / R$.

Increasing the length of the nanowires (by increasing the etching time) leads to higher relative resistance, although the trend appears to saturate; only a relatively small difference between the \SI{5}{\minute} and \SI{7}{\minute} samples has been obtained.
The highest relative change (\SI{37.3}{\%}) in resistance is observed for the longest etching time (sample D) with regards to uniaxial pressure.
\begin{figure
    \includegraphics[width=0.8\linewidth]{Fig/errorbar2-eps-converted-to.pdf}
   \caption{Resistance vs. applied force for the D-like samples (blue line) and bulk-Si (red line). Data for SiNWs are average values for 4 distinct samples, error bars show the range in  which the measurements fell in.

   }
   \label{fig:RandF}
\end{figure}
In \autoref{fig:RandF} the resistance versus applied force is shown again for sample D (blue line), now compared with the case of the bulk Si (red line).
For this measurement we replicated the sample D four times, in the same conditions of 7 min etching time, and performed the measurements on each sample replica.
The data displayed are average values measured on these samples (D-like), shown with the error bars.  A striking difference between these two configurations (SiNWs vs. bulk Si) is observed. No measurable change in resistance is seen
for the bulk Si, in a very sharp contrast to the SiNWs sample.

Up to \SI{35}{\%} change in resistivity was observed over the pressure variation in the range of \SIrange{0.1}{0.9}{\newton}.
In order to test the stability of SiNWs PZR response, we perform several sets of measurements over different samples.
In this stability test the resistance shift was measured in periods of \SI{60}{\second} of loading and unloading
forcet, for all the samples.
All samples were loaded by \SI{860+-10}{\milli\newton} and then unloaded to \SI{220+-10}{\milli\newton}. We observe fast
response, good stability and high sensitivity to repeated pressure changes. The results of a
reproducibility test for sample D is presented in \autoref{fig:load}.

\begin{figure
    \centering
   \includegraphics[width=0.9\linewidth]{Fig/shape4-1-eps-converted-to.pdf}
    \caption{SiNWs PZR response due to repeated load-unload force in real time for sample D. }
    \label{fig:load}
\end{figure}
Hydrogenation (exposure to hydrogen plasma) is widely used in the electronic industry to increase the mobility of charge carriers in semiconductors.
It neutralizes deep and shallow defects and charged surface states \cite{danielsson2010effect,svavarsson2020hydrogenated}.
Because the PZR effect has been attributed to surface states, we applied a hydrogenation treatment in order to explore
the origin of the PZR effect in SiNWs.
For hydrogen plasma treatment we used a custom built inductively coupled discharge setup with cylindrical geometry (\SI{290}{\milli\meter} long quartz tube with diameter of \SI{34}{\milli\meter}).
The quartz tube was held inside a circular copper inductive coil with diameter of \SI{54}{\milli\meter}.
A radio frequency power generator CERSAR (c) (\SI{13.56}{\mega\hertz}) source coupled with impedance matching unit was utilized.
For hydrogenation a gas mixture of Ar/H\textsubscript{2} (\SI{30/70}{\%}) was used and the throttle valves were adjusted to stabilize the gas pressure of \SI{29}{\milli\bar}.
A more detailed description and schematic of the hydrogenation setup can be found elsewhere \cite{sultan2019enhanced,sultan2018effect}.
In \autoref{fig:hydrogen} we compare the behavior of sample D before and after the hydrogenation.
The PZR effect decreases dramatically, and can be attributed to passivation of the surface states.

\begin{figure
    \centering
    \includegraphics[width=0.9\linewidth]{hydrogen1-eps-converted-to.pdf}
    \caption{SiNWs PZR response for sample D before and after hydrogenation. }
    \label{fig:hydrogen}
\end{figure}

\subsection{Electrical response under isostatic pressure}
Sample E (\SI{10}{\micro\meter}) was used for the investigation of the PZR response to isostatic pressure in
a vacuum chamber.
The PZR response was tested in the chamber under re-pressurizing condition in the range of \SIrange{e-4}{e3}{\milli\bar}.
The results are shown in \autoref{fig:r_vs_p} where the resistance is plotted as a function of the pressure measured at a fixed bias voltage of \SI{5}{\volt}.
We observed a dramatic increase of the resistance, by more than two orders of magnitude, when pumping out the air from the
vacuum chamber. We believe this is a result of combined mechanical and chemical effects: the pressure drop removes
mechanical stress, and the lack of air and humidity suppresses the chemisorption.  Note that the resistance values for sample E are much higher than for the samples A-D because of the co-planar configuration of the contacts and also
because of the larger length of the nanowires.

In the Appendix we show, for comparison, the behavior of periodic arrays of SiNWs with increasing isostatic pressure.
\begin{figure
    \centering
 \includegraphics[width=1\linewidth]{Fig/Resistance_vs_pressure.png}
    \caption{The resistance of sample E (on a log 10 scale) as a function of pressure measured in the vacuum chamber, is shown with the blue curve. For comparison, in red color, the resistance of a bulk silicon sample (i.\ e. without nanowires)}
    \label{fig:r_vs_p}

\end{figure}

\section{Discussion}

For solid materials the inter-atomic spacing may be altered by strain.
Consequently, apart from geometric changes, in semiconductors, bandstructure related details, such as bandgap, or effective mass, may change and thereby the resistivity as well.
Within a certain range of strain this relationship is linear \cite{Matsuda1993,Rowe2014}.
When a uniaxial stress $X$ is applied, the piezoresistance coefficient of the resistivity $\rho$ in the direction of stress, is defined as
\begin{equation}
\pi_l=\frac{\Delta\rho}{\rho_0}\frac{1}{X}   \ ,
\label{eq:x1}
\end{equation}
where $\Delta \rho$ is the stress-induced change in the resistivity and $\rho_0$ is the reference resistivity of the unstressed material.

In our case a uniaxial compressive force $F$ was applied on SiNWs along their length by a force meter.
The stress is $X=F/A_{\mathrm t}$, $A_{\mathrm t}$ being the total cross sectional area of the nanowires, which is equal to $pA_{\mathrm m}$, where $A_{\mathrm m}$ is the cross-sectional area of the pin pressed into the wires, and $p$ is the relative cross-sectional area of the wires.
Assuming the electrical resistivity of the nanowires is proportional to their resistance \autoref{eq:x1} becomes

\begin{equation}
   \pi_{l}= \frac{p A_{\mathrm m}}{F}  \times \frac{\Delta R}{R_0}\ .
    \label{eq:our_final_formula}
\end{equation}

The structure of our SiNWs array is robust and stable as the wires are partially interconnected, which provides high structural strength.
Such stability and fast response of SiNWs is a desired  property for many devices such as solid state accelerometers and bipolar transistors \cite{Doll2013}.
Also our results are in agreement with the study of Ghosh et al.~\cite{Ghosh2021}, in which large diameter Si nanorod-based sensors was used for force detection.

He et al. \cite{He2006}, who measured PZR coefficients of single (p-doped) SiNWs, with diameters \SIrange{50}{300}{\nano\meter}, made from wafers with resistivities between \SIrange{0.003}{10}{\ohm\cm}, found that the PZR was roughly inversely proportional to the diameter and proportional to the wafer resistivity.
For a similar diameter and wafer resistivity as in our present work, \SI{\sim 150}{\nm} and \SI{10}{\ohm\cm}, respectively, the PZR coefficient $\pi_l$ for a single wire was of the order of \SI{e-7}{\per\pascal}.
(According to Figure 2d of Ref. \onlinecite{He2006}.)
Using \autoref{eq:our_final_formula}, with $p$\num{= 0.28}, $A_m$\SI{= 49}{\mm\squared}, $F$\SI{= 0.8}{\newton}, and $\Delta R /R$\num{= 0.35}, we obtain $\pi_l$\SI{\approx 6e-6}{\per\pascal}, i.\ e. almost two orders of magnitude higher.
Note that in principle \autoref{eq:x1} does not depend on sample details such as the number of nanowires or their configuration, which we believe plays a role in our case.
Therefore we attribute this higher value to a collective PZR effect brought about by the interaction between multiple and interconnected wires, rather than the response of a single nanowire.
\section{Conclusions}
In summary, large arrays of interconnected SiNWs were fabricated in a simple tree step wet chemical process and
used for testing piezoresistance effect in nanowires.
The interconnected structures of the SiNWs provide a great increase in mechanical stability.
A pressure change of \SI{100}{\pascal} could be measured with this robust device.
The calculated PZR coefficient based on SiNWs array with NWs length of \SI{3}{\micro\meter} as resulted after MACE
etched for \SI{7}{\minute}, sample D, was almost two orders of magnitude higher for our sensor than reported for a
single SiNW.
Repeated measurements for different samples fabricated with the same process demonstrated good reproducibility with
less than \si{5}{\%} deviation in pressure sensing.
The electrical resistance of SiNWs of \SI{10}{\micro\meter} length increased more than two orders of magnitude when
measured in vacuum.
This findings make the device based on random and interconnected SiNWs a strong candidate as a simple and inexpensive
alternative to various pressure sensing applications.


Na atomic order, Co charge disproportionation and magnetism in Na$_{x}$CoO$_{2}$ for large Na contents

Abstract: We have synthesized and characterized four different stable phases of Na ordered Na$_{x}$CoO$_{2}$, for $0.65<x<0.8$. Above 100 K they display similar Curie-Weiss susceptibilities as well as ferromagnetic $q=0$ spin fluctuations in the CoO$_{2}$ planes revealed by $^{23}$Na NMR data. In all phases from $^{59}$Co NMR data we display evidences that the Co disproportionate already above 300 K into non magnetic Co$^{3+}$ and magnetic $\approx $Co$^{3.5+}$ sites on which holes delocalize. This allows us to understand that metallic magnetism is favored for these large Na contents. Below 100 K the phases differentiate, and a magnetic order sets in only for $x\gtrsim 0.75$ at $T_{N}=$22 K. We suggest that the charge order also governs the low $T$ energy scales and transverse couplings.

\section{Introduction}

The interplay between magnetism and metallicity in the copper oxygen plane
structure of the cuprates is at the origin of the great interest which has
been devoted to their high temperature superconductivity. Such an interplay
has been also found recently in the lamellar cobaltates Na$_{x}$CoO$_{2}$,
which further display a high thermoelectric power together with good
electronic conductivity \cite{Terasaki}. In both systems metallicity is
induced by doping the square CuO$_{2}$ or triangular CoO$_{2}$ layers by
insertion of dopants in the ionic separation layers. Soon after the
discovery of the cuprate superconductivity, it has been suggested that the
electron gas could display intrinsic inhomogeneous charge structures, such
as linear ''stripes'' of Cu$^{2+}$ spins separated by doped metallic stripes
\cite{Emery}. However, such static structures have only been found in
specific cases \cite{Tranquada}, and the magnetic properties of the Cu sites
have been usually found rather homogeneous.

A distinct property of the Co ions in the large crystal field induced by
their oxygen octahedral environment in the CoO$_{2}$ structure, is that the $%
t_{2g}$ triplet of the Co site is much lower in energy than the $e_{g}$
doublet, so that the electronic structure of the Co ions is expected to
correspond to low spin configurations Co$^{3+}$ ($S=0$) or Co$^{4+}$ ($S=1/2$%
) obtained by filling only the $t_{2g}$ triplet states. Consequently ordered
charge structures are more frequently expected than for cuprates. Indeed,
using the quasi-unique site sensitivity of NMR techniques, a Na atomic order
\cite{NaPaper} and an associated Co ordered charge disproportionation (OCD)
\cite{CoPaper} has been revealed for a specific metallic $x_{0}\approx 0.7$
composition, which displays a Curie-Weiss susceptibility. However for the
metallic cobaltate antiferromagnetic (AF) phases found for $x\geqslant 0.75$
\cite{Sugiyama,Mendels}, neutron scattering experiments have been analyzed
in a uniform local moment picture \cite{Boothroyd,Bayrakci}. This is
unexpected for hole doping of Na$_{1}$CoO$_{2}$ which is a band insulator
built from the non magnetic Co$^{3+}$ state \cite{Lang}, for which the $%
t_{2g}$ multiplet is filled by the six $d$ electrons.

Here we present the first systematic effort undertaken so far to study with
local probes the variation with $x$ of the magnetic properties and
demonstrate that the OCD is generic for $x>0.65$ including the AF phase with
$T_{N}=$22~K. This OCD always results into non magnetic Co$^{3+}$ and a
second type of sites with formal valence of about 3.5 responsible for the
peculiar magnetic properties. We establish that the latter display above
100~K a generic nearly ferromagnetic behaviour for four phases displaying
distinct ordered Na structures, so that the OCD is an essential ingredient
to explain the magnetic properties of these high $x$ phases. Surprisingly
the low $T$ ground state magnetic properties differ markedly as, apart the
AF phase, the former $x_{0}\approx 0.7$ phase appears an experimental
realization of a nearly ferromagnetic 2-dimensional metal without static
magnetic order down to $T=0$. The hole doped AF phases have been shown
ferromagnetic in plane and AF between planes (A-type AF) \cite
{Boothroyd,Bayrakci}. So, the Na atomic order and the OCD of Co appear
essential as well in governing the low energy modifications of the band
structure of the correlated metallic state, which drive the coupling between
Co planes and the low $T$ magnetic and thermoelectric properties \cite{Ong}.

\begin{figure}[tbp]
\onefigure[width=1\linewidth]{FigXray.eps}%
\caption{Hexagonal (P63/mmc N$^\circ$ 196) indexation of part of the X ray
diffraction patterns. At these high angle values the Cu K$\alpha $1 and
K$\alpha $2 splits the reflections into two separate peaks with intensity
ratio 2/1. (*) In the 0.71 case an additional structural splitting separates
the hexagonal 110 (and 112) reflections into two distinct 310 and
020 orthorhombic (Ccmm N$^\circ$ 63) peaks.}%
\label{FigXray}
\end{figure}

\section{Samples}

We have synthesized our samples by standard solid state reaction of Na$_{2}$
CO$_{3}$ and Co$_{3}$O$_{4}$ powders in flowing oxygen, with nominal
concentrations $x$ increasing by increments of 2 to 3/1000 in the range
$0.65\leqslant x\lesssim 0.80$. X ray powder diffraction data always
exhibited the Bragg peaks corresponding to the two-layer Co structure with a
hexagonal unit cell. However we systematically detected weak additional
reflections indicative of three dimensional Na long range ordering. It has
been immediately clear that the corresponding Na order is complicated and
highly dependent on Na content. This ordering was even found for a specific
composition to drive an orthorhombic distortion of the average lattice as
shown for $x=0.71$ in fig.~\ref{FigXray}, on a detail of the Bragg peaks of
the X ray diffraction profile. From the Rietveld refinements of the X ray
data summarized in table~\ref{table1} we could isolate four distinct phases
for $0.65\leqslant x\lesssim 0.80$. Each phase exhibits a specific Na
ordering leading to characteristic additional diffractions: commensurate
reflections, or incommensurate superstructure satellites with $q(b^{\ast})$
as the component of the wave vector modulation.

The concentrations for which single phase samples could be stabilized are
sequenced in four distinct narrow $x$ non overlapping domains outside which
the synthesized samples were found as intermediate mixtures of these phases.
In addition, multiphasing has been found to occur quite commonly for any $x$
value if no particular care to homogeneize the samples was taken. We could
however synthesize these four phases reproducibly, forbidding any air
exposure of the samples. Confirmation that nearly pure phases could be
achieved on large mass samples ($\simeq $400~mg) required for NMR
measurements has been directly obtained from $^{23}$Na and $^{59}$Co NMR
data, as shown hereafter.

The orderings found for these four phases always correspond to a symmetry
lowering and is systematically based on an orthorhombic reference subcell
(Ccmm, N$^{\circ }$63): $a_{ort}=a_{hex}\sqrt{3}$; $b_{ort}=a_{hex}$; $%
c_{ort}=c_{hex}$. Owing to the orthorhombic distortion for $x=0.71$ the
labelling H67, O71, H72 and H75 has been used hereafter and in table~\ref
{table1}. The formerly studied $x_{0}\approx 0.7$ phase \cite
{NaPaper,CoPaper} is in fact that with the lowest $x$ value $x=0.67(1)$. We
indeed found that any sample with larger $x$ evolves towards this limiting
composition if it is kept in insufficiently dry atmosphere.

\begin{largetable}%
\caption{Parameters of the studied phases. The accuracy on their difference
of Na content is much better than that on $x$ itself ($\pm $ 0.01). The
lattice constants $a,b,c$ corresponding to fig.~\ref{FigXray} are given in
the hexagonal or orthorhombic reference cell. When determined, the
incommensurate $q(b^{\ast})$ or commensurate modulations $q(c^{\ast})$ are
given. The Co$^{3+}$ fraction $y$ is obtained from $^{59}$Co NMR intensity
data. The 0.67 sample displays a commensurate
orthorombic superstructure $(a_{hex}\sqrt{3},\;3a_{hex},\;3c_{hex})$.}%
\label{table1}.
\begin{center}
\begin{tabular}{cccccccccc}
Phase & $x$ & $a_{ort}/\surd{3}$ (\AA) & $a_{hex}$ or $b_{ort}$ (\AA) & $c$ (\AA) & $q(b^{\ast})$ & $q(c^{\ast})$ & $A_{eff}^{iso}$ (kG/$\mu _{B}$)& y (\%) & \\
\hline
H67 & $0.67$ & $a_{hex}$ & 2.82920(1) & 10.9387(4) & $1/3$ & $1/3$ & 9.1(3) & 26(4) &  \\
O71 & $0.71$ & 2.83931(2) & 2.83031(3) & 10.8929(2) & 0.2849(1) & 0 & 8.0(3) & 40(5) &  \\
H72 & $0.72$ & $a_{hex}$ & 2.83651(4) & 10.8770(2) & 0.2810(1) & 0 & 7.3(3) & 37(5) &  \\
H75 & $\geq 0.75$ & $a_{hex}$ & 2.84162(1) & 10.8058(3) & - & - & 7.8(3) & 33(4) & \\
\end{tabular}
\end{center}
\end{largetable}

\section{Magnetic susceptibility and $^{23}$Na NMR shifts.}

The single crystal grains of samples of these phases were oriented in the $%
H_{0}=7$ Tesla NMR field within Stycast or paraffin. SQUID measurements of
the macroscopic susceptibility $\chi _{m}$ taken in 5~T field allow us to
evidence that the different phases display different magnetic properties.
For instance, as evidenced in fig.~\ref{fig.2}a, the low T magnitude of
$\chi_{m}$ decreases progressively with increasing $\emph{x}$. The H75 phase
is furthermore found to be the only phase in which a magnetic order is
detected in low applied field. However, as minute amounts of impurities or
slight admixture of phases could spoil the bulk measurements, spectroscopic
measurements with local probes better determine the susceptibility of each
phase.

As reported in ref.~\cite{NaPaper} on H67, the $^{23}$Na NMR is an excellent
probe of both Na order and of the magnetic properties. Indeed the $^{23}$Na
NMR displays distinct quadrupole splittings for the different Na sites of
this structure. In view of the more complex structures of the new phases
there was no surprise in finding a larger number of resolved Na sites,
although with similar magnitudes of their quadrupole splittings. The
magnetic properties of the compounds are probed at the local scale through
the NMR shifts of the different Na sites resolved in the $(-\frac{1}{2}%
\leftrightarrow \frac{1}{2})$ transition of the $^{23}$Na spectra presented
in fig.~\ref{fig.1}. There one can see that the $T=$5~K spectra are quite
distinct for the four phases. For H75 a large broadening occurs in the AF
state below $T_{N}=22$~K, while the spectrum of H67 is much more shifted
than those of O71 and H72, which are distinct but display some overlap. The
H72 batch has been found to evolve in time at room $T$ towards O71, and a
slight mixture of the two pure phases might be unavoidable in the bulk
samples \footnote{In various published papers for x values in the range
studied, e.g. \cite{Sakurai}, we could find on the susceptibility data many
signs indicating that the samples were superposition of our phases.}. Quite
generally $^{23}$Na NMR allowed us to control the phase purity of the bulk
of the NMR sample, as multiphase samples display superimposed spectra, as
seen in fig.~\ref{fig.1}.

\begin{figure}[tbp]
\onefigure[width=1\linewidth]{Fig2epl.eps}%
\caption{(a) $T$-dependencies of the bulk susceptibility $\chi_{m}$ measured
with a DC SQUID in a 5~T field; (b) $T$ dependence of the mean $^{23}$Na NMR
shift. Identical behaviour above 100~K can be seen for the four phases with
remarkable differences at low $T$; (c) the linear variations of $K_{s}$
versus $\chi_{m}$ underlines the phase purity of the samples. The H67 data
for $T<$30~K \cite{NaPaper} has been omitted to better display
those of the new phases.}%
\label{fig.2}
\end{figure}

Let us recall, as detailed in ref.~\cite{NaPaper} that for a field $%
H_{0}\parallel \alpha $, the NMR shift $K_{\beta}^{\alpha}$ of a Na atomic
site $\beta$ probes the spin susceptibility $\chi _{s,i}^{\alpha }(T)$ of
the neighbouring Co sites $i$ through transferred hyperfine couplings $%
A_{\beta ,i}^{\alpha }$ with $K_{\beta }^{\alpha }=\sum_{i}A_{\beta
,i}^{\alpha }\;\chi _{s,i}^{\alpha }(T)$. The main result found for H67, and
verified as well here for the other phases, is that the $K_{\beta }^{\alpha
}(T)$ variations scale with eachother for all Na sites. This $T$ dependence
is associated with the average $\chi _{s}^{\alpha }(T)$ of the magnetic Co
sites of the structure. So, overlooking the diversity of Na sites, the first
moment (or center of gravity) of the $^{23}$Na NMR signal writes $%
K_{s}^{\alpha }=A_{eff}^{\alpha }\chi _{s}^{\alpha }(T)$ where $%
A_{eff}^{\alpha }$ is an effective hyperfine field per Co site. The $T$
variations of $K_{s}^{\alpha }$ are reported in fig.~\ref{fig.2}b, and are
shown to be quite identical for $T>100$~K with a unique Curie-Weiss
$(T-\Theta)^{-1}$ variation (with $\Theta\approx$-80~K). They differ
markedly below 100~K, as does the SQUID data for $\chi _{m}$, the low $T$
enhancement of $\chi _{s}^{\alpha }(T)$ observed for H67 being progressively
reduced for increasing $x$.

The usual comparison between the SQUID and Na NMR data, displayed in the
good linear $K_{s}$ versus $\chi _{m}$ plots of fig.~\ref{fig.2}c allows us
to confirm the purity of the isolated phases. The high-T slopes of these
plots yield similar values for the effective hyperfine coupling
$A_{eff}^{iso}$ (table~\ref{table1}), which could be expected as $^{23}$Na
sites are coupled with many cobalts \cite{NaPaper}. In all phases the
anisotropy of $\chi _{s}^{\alpha }(T)$, given by that of $K_{s}^{\alpha }$,
has been found $\lesssim $ $\pm $ 0.1.

\begin{figure}[tbp]
\onefigure[width=1\linewidth]{Fig1epl.eps}%
\caption{$^{23}$Na NMR central line spectra taken at 5~K. They are quite
distinct for the four nearly pure phases, with some overlap between O71 and
H72 spectra. That for a sample which is clearly a mixture of H67, O71 and
H72 is shown as well. $\nu_{ref}$ is the non magnetic $^{23}$Na NMR
reference.}%
\label{fig.1}
\end{figure}

In the H75 AF phase, the saturation of $K_{s}(T)$, that is $\chi
_{s}^{\alpha }(T)$, seen at low $T$ in fig.~\ref{fig.2}b should be
associated with the onset of AF correlations. In a uniform Heisenberg model,
one would then assign the progressive increase of $\chi _{s}^{\alpha }(T)$
at low $T$ with decreasing $x$ to a decrease of $T_{N}$ and of out of plane
AF coupling strength. However this primary interpretation fails as NMR data
taken down to 1.4~K (and $\mu $SR to 50~mK \cite{Mendels2}), did not
evidence any frozen magnetic state in the three other phases, which are then
\textit{paramagnets in their ground state}, most probably metallic, as no
electronic magnetic transition is detected.

\begin{figure}[tbp]
\onefigure[width=1\linewidth]{Fig3epl.eps}%
\caption{Spectra taken with a large pulse spacing $\tau =200~\mu s$ in a
spin echo sequence allowed us to isolate the narrow spectra of Co1 sites
with long $T_{2}$. The broad spectra of the magnetic Co2 sites with short
$T_{2}$ are obtained by subtracting the Co1 spectra from those taken with
$\tau =10~\mu s$. The average shift of the Co2 sites decreases with
increasing $x$ as does the $^{23}$Na shift in figs.~\ref
{fig.2}b and \ref{fig.1}.}%
\label{fig.3}
\end{figure}

\section{$^{59}$Co NMR}

The difference between phases appears as well in the $^{59}$Co NMR spectra,
which also display more Co sites than for H67 \cite{CoPaper}.

\subsection{Charge disproportionation}

As in H67 we identified two classes of Co sites, their $(-\frac{1}{2}%
\rightarrow \frac{1}{2})$ transitions being easily visualized when $H_{0}$
is applied at 54.7$^{\circ }$ with respect to the $c$ axis, the ''magic
angle'' for which quadrupole effects are reduced, as seen in
fig.~\ref{fig.3}. A first series which we label as the Co1 ''class'', is
associated with non magnetic Co$^{3+}$, with small spin lattice $T_{1}^{-1}$
and spin spin $T_{2}^{-1}$ relaxation rates, and occurs with similar NMR
shifts at low $T$ in all four phases. The more complex spectra of fast
relaxing magnetic Co sites, which we label as the Co2 ''class'' is seen in
fig.~\ref{fig.3} to include diverse sites with much larger shifts at low $T$
than the Co1 sites. The increase of average shift of these Co2 with
decreasing $x$ agrees perfectly with the SQUID and $^{23}$Na NMR data
(fig.~\ref{fig.2}).

The Co1 NMR shifts were found to vary with temperature, which is a sign that
these non magnetic sites are sensing the magnetism of the Co2 sites through
transferred hyperfine couplings. Indeed, in Fig.~\ref{fig.5}, we evidence
that the shifts $^{59}K_{1}^{\alpha }$ of the Co1 nuclei scale linearly at
low $T$ with that of $^{23}K$. The  $T$ independent orbital contribution to
$^{59}K_{1}^{\alpha }$ can be obtained by extrapolating this linear
dependence to $^{23}K=0$ (that is vanishing spin susceptibility $%
\chi _{s}(T)$).\ It is found isotropic and increases slightly from 1.95 to
2.05\% from H67 to H75. These values are quite comparable with the purely
orbital shift of 1.95\% found for Co$^{3+}$ in the band insulator Na$_{1}$CoO%
$_{2}$ \cite{Lang}. The fraction $y$ of Co$^{3+}$ sites, estimated from the
Co1 relative NMR intensity (corrected for $T_{2}$ decay), increases
slightly, but not regularly with $x$ (table~\ref{table1}), this overall
trend being expected as all Co sites become Co$^{3+}$ for $x=1$. However as
we find $y<x$, the average valence of the Co2 class of sites is always $<4+$%
, and the OCD detected in H67 is present in all phases, \textit{including
the AF ordered phase}.

\begin{figure}[tbp]
\onefigure[width=1\linewidth]{Fig5epl.eps}%
\caption{At low $T$, the Co1 sites NMR shifts display linear variations
versus the $^{23}$Na shifts $^{23}K$. For $H_{0}\parallel c$, the Co1
hyperfine coupling\ with the magnetic Co2 sites have similar magnitudes for
all samples. For $H_{0}\perp c$ it has a similar magnitude for H72 and H75
but nearly vanishes for H67 and O71. For H67 two Co1 sites are well
resolved, and their shifts depart slightly below 40~K (i.e. $^{23}K>0.2\%$),
as reported in ref.~\cite{Gavilano}. The upturns for small $^{23}K$ are due
to the onset of Co1-Co2 site exchange due to Na motion, as discussed in the
text.}%
\label{fig.5}
\end{figure}

\subsection{Na motion}

As seen in Fig.~\ref{fig.5} we could detect the Co1 sites up to room $T$ in
most phases, well above the onset of Na motion detected hereafter at $%
\approx $200~K from $^{23}$Na $T_{1}$ data \cite{MotionPaper}. This implies
that the OCD occurs already above room $T$, contrary to the proposal of
ref.~ \cite{Gavilano}. However in fig.~\ref{fig.5} significant increases of
$^{59}K_{1}^{\alpha }$ with respect to the low $T$ linear $^{23}K$
dependence are observed for $T\gtrsim 150$~K, that is for $^{23}K<0.08$. We
attribute this to the onset of Co1-Co2 site exchange which becomes
significant when Na motion begins to set in. Such site exchanges are usually
easily detected in NMR spectra as they yield a progressive reduction of the
line splitting between the two sites until they merge in a single exchange
narrowed line at very high $T$, when the exchange rate exceeds the line
splitting. Here the increase of the Co1 NMR shift merely corresponds to the
onset of Co1-Co2 site exchange. Indeed from ref.~\cite{CoPaper},
$K^{ab}$(Co2) is as large as $\approx $4\% above 200~K, so the tiny
0.1\%-0.2\% increase of $K^{ab}$(Co1) corresponds to a $\approx $10\%
decrease of the 1.5~MHz Co1-Co2 splitting, that is a very small exchange
rate $\tau _{ex}^{-1}\approx 150$~kHz. This analysis is validated by the
fact that weaker upturns of $K^{c}$(Co1) are seen for $H_{0}\parallel c$ for
which $K^{c}$(Co1) and $K^{c}$(Co2) happen to differ only slightly as seen
in ref.~\cite{CoPaper}.

Therefore this key observation that the rate of exchange $\tau _{ex}^{-1}$
of Co sites is very slow with respect to expectations for electronic
processes allows us to conclude that the Co1-Co2 site exchange is connected
with Na motion. This proves that \textit{the Co charge does correlate with
the Na environment} (e.g. Na1 sites being on top of Co$^{3+}$).

\section{$\chem{^{23}}$Na spin lattice relaxation and electron spin dynamics}

To search for differences in the dynamic electronic susceptibilities of
these phases we have taken extensive $^{23}$Na $T_{1}$ data. As $^{23}$Na
has a spin $I=3/2$, its nuclear magnetization recovery should be given by

\begin{equation}
M(t)/M_{0}=1-W\exp (-6t/T_{1})-(1-W)\exp (-t/T_{1}),
\end{equation}

with $W=0.9$ if only the central transition has been saturated. Since
condition being impossible to fulfill strictly experimentally, $W$ has been
left as an adjustable parameter, which was found to evolve between 0.9 and
0.7 depending on the sample and experimental conditions. The $T_{1}^{-1}$
data were found slightly anisotropic, i.e. $\approx $30\% larger for
$H_{0}\perp c$ than for $H_{0}\parallel c$, for $T<$200~K. So, in fig.~\ref
{fig.4} we only plotted the data for $H_{0}\parallel c$.

Although the low $T$ variations of $\left( T_{1}T\right)^{-1}$ differ for
the four phases, they do become identical above 100~K as seen in fig.~\ref
{fig.4}a. There, we assign the extra high $T$ contribution to the occurrence
of Na motion, their onset taking place at distinct $T\gtrsim 200$~K for the
different phases. As can be anticipated, we could check that the magnetic
behaviour is not that of a Fermi liquid, especially for the H67 sample, as
$\left(T_{1}T\right)^{-1}$ and $^{23}K(T)$ do not combine into a constant
$R=S_{0}/T_{1}T(^{23}K)^{2}$, where $S_{0}=(\hbar /4\pi k_{B})(\gamma
_{e}/\gamma _{n})^{2}$ is the universal Korringa ratio. In fig.~\ref{fig.4}b
$R$ is of course sample independent between 80 and 160~K and always smaller
than unity, as noticed in a $x\approx 0.7$ mixed phase sample by Ihara
\textit{et al.} \cite{Ishida}. This is expected if a ferromagnetic
quasi-elastic peak at \textbf{q}$\approx 0$ dominates the spin excitations
as revealed by Inelastic Neutron Scattering (INS) for the H75 phase above
$T_{N}$ \cite {Boothroyd}. Indeed, in such a case the \textbf{q}$\approx 0$
response enhances markedly $\chi _{s}($\textbf{q}$=0)$, that is
$^{23}K_{s}$, while $\left( T_{1}T\right) ^{-1}$\ is less enhanced as it
probes $\chi ^{\prime \prime }($\textbf{q}$,\omega)$ at all \textbf{q}
values. The identical behaviour found here above 100~K extends this result
to all our Curie-Weiss phases.

\begin{figure}[tbp]
\onefigure[width=1\linewidth]{Fig4.eps}%
\caption{$T$ variation of
$(T_{1}T)^{-1}$ (a) and the normalized Korringa product $R$ (b) for the four
phases. While the data are distinct below 100~K, in (c) a universal scaling
between $(T_{1}T)^{-1}$ and $^{23}K$ is shown to apply. The high $T$
deviations due to Na motion and the slight low $T$ increases for H72 and O71
are discussed in the text.}%
\label{fig.4}
\end{figure}

To characterize further the difference of magnetic excitations in the
diverse ground states we searched for relationships between $(T_{1}T)^{-1}$
and $^{23}K$, as done in fig.~\ref{fig.4}c. Spin fluctuation theories in
nearly ferromagnetic metals are expected to give $(T_{1}T)^{-1}=aK^{n}$
\cite{Moriya-Ueda}, with $n=1$ in 3D \cite{Alloul-Mihaly}, while $n=3/2$ is
expected for 2D \cite{Hatatani-Moriya,Kitagawa}. For H67, we do remarkably
find an accurate scaling, with $n=1.5\pm 0.1$, that is $%
(T_{1}T)^{-1}=aK^{3/2}$ over the entire range 1.5~K$<T<$300~K. This
interpretation of the data implies that  the Curie-Weiss temperature
$\Theta$=-80~K found in $K(T)$ and $\chi(T)$ is not associated with AF
in-plane couplings but results from the nearly ferromagnetic electronic band
behaviour \cite{Moriya-Ueda,Hatatani-Moriya} as seen for instance in
TiBe$_{2}$ \cite{Alloul-Mihaly}. Impressively this 2D \textit{ferromagnetic}
scaling extends down to low $T$ with a unique $a$ for all phases \footnote{
For H72 and O71, the small deviations for $T\lesssim 5$K are not intrinsic
as data for the resolved Na lines were not found to scale perfectly.
Purification of these phases is required for accurate low $T$ studies.},
including the AF phase down to $T_{N}$.

At this stage one might wonder whether the saturation of $^{23}K(T)$
observed below 100~K in H75 is related to AF plane to plane couplings. Such
AF fluctuations that enhance $\chi ^{\prime \prime }(\mathbf{q}=\mathbf{q}%
_{AF})$ should result in an increase of $(T_{1}T)^{-1}$ and $R$ with a
divergence at $T_{N}$. But they are not probed by the $^{23}$Na nuclei, as
the local fields induced by two adjacent Co layers cancel in the A type AF
structure, as confirmed by the weak $^{23}$Na NMR shift in the N\'{e}el
state (fig.~\ref{fig.1}). The $^{23}$Na $T_{1}$ only probes then the
strength of the ferro fluctuations, and the perfect scaling found above
proves that the main incidence of AF fluctuations is to reduce the ferro
ones below 100~K. Comparison with the numerical results of Hatatani \emph{et
al.} \cite{Hatatani-Moriya} allows us to point out that, for H67, the low
$T$ increase of $\chi (\mathbf{q}=0)$ with respect to the common Curie-Weiss
variation is that expected by these authors in the immediate vicinity of a
ferromagnetic instability. Therefore the H67 phase \textit{appears as an
ideal 2D nearly ferromagnetic metal without 3D ordering settling in at low}
$T$.

\section{Discussion}

The four phases exhibit specific Na orderings, but similar charge
disproportionation, and identical ferro in plane fluctuations above 100~K,
independent of the details of the OCD. It seems to us that this analogy of
properties relies heavily on the occurrence of non magnetic Co$^{3+}$ as it
reduces the number of hopping paths between magnetic sites with respect to a
homogeneous structure. The associated decrease of bandwidth $W$ and increase
of hole density on the magnetic sites magnify the role of correlations $U$,
which enhances the ferro tendency supported by LDA calculations \cite{Singh}
for the uniform case.

On the contrary the ground state magnetic properties are certainly driven by
the diverse Na atomic orders evidenced here, similar to those suggested, or
observed \cite{Zandbergen,Roger} as well for various other $x$ values. One
could expect then distinct ordered magnetic states as in the case of the $%
x=0.5$ phase for which an AF order inside the Co plane sets in at $T_{N}=$
86~K \cite{Mendels,Gasparovic} in the absence Co$^{3+}$ \cite{Bobroff}. But
for the analogous ferro in plane couplings evidenced here for large $x$, one
would always expect A type AF order to occur. In these AF phases, the 3D
dispersion of the spin wave excitations found by INS has been analyzed with
Heisenberg Co-Co AF coupling between planes either with nearest
\cite{Boothroyd,Bayrakci} or next nearest neighbour exchange through Na
orbitals \cite{Mazin}. For such Heisenberg transverse couplings, a hectic
evolution of $T_{N}$ versus $x$ should presumably result, depending on the
actual Na order, contrary to the smooth evolution of \textit{paramagnetism}
found for $x>0.65$ and the abrupt occurrence of AF above $x=0.75$.

This definitely allows us then to conclude that metallic magnetism is indeed
responsible as well for the low $T$ states for large $x$ in these Na
cobaltates. One might consider \cite{McKenzie} that a Fermi liquid state is
only reached below an energy scale given by the temperature $T^{\ast }$ at
which $\chi _{s}$ saturates, which increases from $\approx $1~K for H67
\cite {NaPaper} to $\approx$ 60~K for H75. The band parameters associated
with Na order would then be responsible for these $T^{\ast }$ values and for
the transverse couplings which drive the AF order. While some attempts have
been done to take into account both correlations and OCD \cite{Marianetti},
extensions to diverse Na atomic orders are required to explain the evolution
with $x$ of the ground state properties \footnote{After the disclosure of an
earlier version of these results \cite{CondMat}, a theoretical attempt to
explain them within a strong correlation approach has been done \cite{Gao}.
The existence of a critical doping value for the onset of in-plane
ferromagnetism has been reproduced, as well as the occurrence of Co$^{3+}$
in the presence of Na disorder, using heuristic values for the Na-Co
interaction potentials.}.

An important experimental aspect revealed by our thorough investigation is
that, contrary to the case of cuprates for which dopant induced disorder is
quite influential, and governs many aspects of the cuprates phase diagram
\cite{FRA-alloul}, the hole doping achieved in bulk cobaltate samples by
insertion of ordered Na planes corresponds to rather clean situations for
many Na concentrations.

\acknowledgments We thank J.~Bobroff, F.~Bert and P.~Mendels for performing
the $\mu$SR measurements on the paramagnetic phases and for helpful
discussions, as well as G.Kotliar, I.~Mazin, F.~Rullier-Albenque and
D.~Singh for their stimulating interest. We acknowledge financial support
from the ANR (NT05-4-41913), INTAS (04-83-3891) and RFBR (06-02-17197).

Strongly Correlated Electrons (cond-mat.str-el); Superconductivity (cond-mat.supr-con)