The Road to Quasars

Abstract: Although the extragalactic nature of 3C 48 and other quasi stellar radio sources was discussed as early as 1960 by John Bolton and others, it was rejected largely because of preconceived ideas about what appeared to be unrealistically high radio and optical luminosities. Not until the 1962 occultations of the strong radio source 3C 273 at Parkes, which led Maarten Schmidt to identify 3C 273 with an apparent stellar object at a redshift of 0.16, was the true nature understood. Successive radio and optical measurements quickly led to the identification of other quasars with increasingly large redshifts and the general, although for some decades not universal, acceptance of quasars as the very luminous nuclei of galaxies. Curiously, 3C 273, which is one of the strongest extragalactic sources in the sky, was first cataloged in 1959 and the magnitude 13 optical counterpart was observed at least as early as 1887. Since 1960, much fainter optical counterparts were being routinely identified using accurate radio interferometer positions which were measured primarily at the Caltech Owens Valley Radio Observatory. However, 3C 273 eluded identification until the series of lunar occultation observations led by Cyril Hazard. Although an accurate radio position had been obtained earlier with the OVRO interferometer, inexplicably 3C 273 was initially misidentified with a faint galaxy located about an arc minute away from the true quasar position.

\section{Introduction}

The discovery of radio galaxies and then quasars had a profound impact on development of astrophysics and cosmology during the later part of the 20th century.  Suddenly astronomers were faced with previously unimaginable luminosities coming from unimaginably small volumes.  Following the 1963 identification of 3C 273 at z=0.16, redshifts were quickly extended to z $\sim$ 2 and more.  But quasars could have, and indeed should have, been recognized earlier.  Indeed 3C 48 had been identified two years before 3C 273, but the apparent  redshift of 0.37 was rejected in favor of a galactic interpretation because of preconceived ideas about the limits to radio and optical luminosity.

Early suggestions about unusual activity in the nuclei of galaxies go back to the work of Sir James Jeans \cite{J61}.  But, the modern understanding of the important role of galactic nuclei probably began with the famous paper by Karl
 \cite{Seyfert43}
who reported the presence of broad strong emission lines in the nuclei of seven spiral nebulae.  Interestingly, although Seyfert's name ultimately became attached to the general category of galaxies with broad nuclear emission  lines associated with highly ionized elements, his 1943 paper apparently went unnoticed until Baade and Minkowski \cite[(1954)]{BM54} drew attention to the similarity of the optical emission line spectrum of the galaxies studied by Seyfert with the galaxy that they had identified with the Cygnus A radio source.

Earlier John Bolton, Gordon Stanley, and Bruce Slee
\cite[(1949)]{BSS49}
had identified three of the strongest discrete radio sources with the Crab Nebula, M87, and NGC 5128,  Until that time the discrete radio sources were widely thought to be associated with stars.  This was understandable, as Karl Jansky and Grote Reber had observed radio emission from the Milky Way.  The Milky Way is composed of stars, so it was natural to assume that the discrete radio sources had a stellar origin.  Bolton et al. understood the importance of the identification of the Taurus A radio source with the Crab Nebula which was widely recognized as the remnants of the 1054 supernova reported by Chinese observers. They correctly identified two other strong sources with M87 and NGC 5128, but realizing that their absolute radio luminosity would need to be $10^6$ times more luminous than that of the Crab Nebula, they  concluded that

 \begin{quote}
\it{NGC 5128 and NGC 4486 (M87) have not been resolved into stars, so there is little direct evidence that they are true galaxies. If the identification of the radio sources are accepted, it would indicate that they are [within our own Galaxy]}.
\end{quote}

They entitled their paper simply ``Positions of Three Discrete Radio Sources of Galactic Radio Frequency Radiation.''  John Bolton later explained that he really did understand that M87 and NGC 5128 were extraordinarily luminous extragalactic radio sources, but that he was concerned that that  a conservative Nature referee might hold up publication (personal communication from Bolton to the author, August, 1989).

The following years saw the identification of more radio galaxies, and the paradigm which had previously considered all discrete radio sources to be galactic stars quickly changed to one where to most high latitude sources were  assumed to be extragalactic.  The energy requirements were exacerbated with the identification of Cygnus A, the second strongest radio source in the sky with a magnitude 18 galaxy  which was at what was then considered a high redshift of 0.056 and a corresponding radio luminosity about $10^3$ times more luminous than M87 or NGC 5128 \cite[(Baade and Minkowski 1954)]{BM54}.

By the end of the decade, many radio sources had been identified with optical counterparts, but the redshift of Cygnus A remained the largest of any known radio galaxy.  Typically the optical counterparts of strong radio sources were elliptical galaxies \cite{B60} that were the brightest member of a cluster.  In 1960, Rudolph Minkowski \cite[(1960)]{M60} identified 3C 295 with a magnitude 20 galaxy at z=0.46, by far the largest known redshift of any galaxy.  3C~295 is about ten times smaller than Cygnus A and ten times more distant, consistent with the idea that the smallest radio sources might be path finders to finding very distant galaxies. But a few months later Caltech radio astronomers identified the first of several very small radio sources with what appeared to be galactic stars, thus again raising questions about the extragalactic nature of other small diameter radio sources.

\section{The First Quasars}

In late 1960, while searching for ever more distant radio galaxies, Caltech radio astronomers, John Bolton and Tom Matthews identified 3C 48 with a faint magnitude 16 apparent stellar object.  Apparently Bolton suggested that 3C 48 might not be a galactic star but rather extragalactic with z=0.37.  But, shortly later Bolton left Caltech to take up his new responsibilities as the person in charge of completing the construction of the new 210-ft radio telescope at Parkes, which completely occupied him for the next few years.  At the 107th meting of the American Astronomical Society held in December, 1960, Allan Sandage \cite{S60} gave a late paper on "The First True Radio Star."

In an exhaustive unpublished  analysis of the complex 3C 48 emission line spectrum, Greenstein \cite{G62}  interpreted the observed emission lines in terms of highly ionized states of various rare earth elements.  He briefly speculated on a possible redshift of 0.37, but then dismissed the possibility that 3C 48 was extragalactic.   Nearly two years would pass, and two other compact radio sources, 3C 196 and 3C 286 were also identified as galactic stars.  Based on OVRO interferometer measurements which had a nominal position accuracy of  five or ten arcsec, in his 1962 PhD thesis, Caltech graduate student Dick  \cite{R62}  suggested an identification for 3C 273 with a faint galaxy,  which later turned to be located about an arc minute east of the currently accepted quasar position.  In May 1962, Maarten Schmidt tried to measure the position of Read's faint galaxy thought to be identified with 3C 273.  However, in his subsequent publication in the Astrophysical Journal submitted in December, 1962,  \cite{R63} dropped the suggested 3C 273 identification.    Although Read and other Caltech radio astronomers had measured the position of 3C 273 good to better than 10 arcsec, and the  existence of star-like optical counterparts was already well established, not until after the 1962 occultations at Parkes \cite{H63} showed the double source structure, would 3C 273 be correctly identified by \cite{S63} with a star-like object at a redshift of 0.16. Like 3C~48, the optical counterpart of 3C~273 contained a faint jet-like feature extending some 20 arcseconds away from the bright stellar object.

Interestingly, the first occultation positions  that were sent by John Bolton to Schmidt on August 20, just two weeks after the occultation, were in error by about 15 arcseconds. The correct radio positions were not known by Schmidt until after he received Bolton's letter of January 26, 1963, more than a month after Schmidt had observed the 3C 273 spectrum at Palomar.   Further details of the events surrounding the Parkes occultation observations and the subsequent interaction with Caltech and Palomar astronomers that led to the identification of the 3C 273 and the determination of its redshift are given in a companion paper in this volume by \cite{HJG15}, by Schmidt \cite{S83,S11} and by \cite{K13, K14}.

Applying the derived 3C~273 Balmer series redshift of 0.16, Schmidt also recognized the Mg II line at a rest wavelength of 2798 Angstroms, leading he and Greenstein  to determine that 3C 48 indeed was also extragalactic and had a  a redshift of 0.37 \cite{GM63}, just the value that was suggested two years earlier by Bolton, but rejected by Greenstein and Bowen.  The announcement of the 3C 48 redshift appeared as a companion paper to the Hazard et al. occultation position and the \cite{S63}  identification  with the quasar and its jet-like extension.  Although Schmidt's measurement and interpretation of the 3C 273 spectrum appears to have been critical in recognizing the extragalactic nature of 3C 48,  the 3C 48 paper appeared without Schmidt's or Bolton's names, but included Tom Matthews as an author, as Greenstein remembered that it was Matthews who first drew his attention to the optical counterpart \cite{G96}.

In a 1989 lecture, John  \cite{B90} recalled that before leaving for Australia in November 1960, he had argued that 3C 48 had a redshift of 0.37 but that  Ira Bowen and Jesse Greenstein had convinced him that a 3 or 4 Angstrom discrepancy in the calculated rest wavelength of different lines was unacceptable.   Bolton's 1989 claim was immediately attacked  \cite{G96} as an attempt to rewrite history.  However, Bolton's assertion was substantiated by published memories of Fred \cite{H89}  and the much later discovery of a handwritten letter from Bolton to Joe Pawsey dated November 16, 1960 stating that 3C 48 had a redshift of 0.37.  However, a few weeks later when his ship stopped in Hawaii, Bolton again wrote to Pawsey saying "It's most likely a star."

\section{Radio Quiet Quasars}

The years following the recognition of the 3C 48 and 3C 273 redshifts led to the identification of more quasars at ever larger redshifts and the recognition that quasars are the extremely bright nuclei of active galaxies.  Generally, the identified quasars had a significant UV excess compared with stars, so due to the redshift of their spectrum, they appeared blue on photographic plates facilitating their identification with radio sources with even modest position accuracy \cite{RS64}.

   In 1965, Sandage noted that the areal density of blue stellar objects was some thousand times greater than that of 3C radio sources.  Sandage argued that what he called quasi stellar galaxies were related to quasars, except that they were ``radio quiet.''  Sandage's paper was widely attacked, perhaps in part because it was received on May 15, 1965 at the Astrophysical Journal, but S. Chandrasekar, the ApJ  editor was apparently so impressed by Sandage's claim for a ``New Constituent of Universe'' that he held up publication of the Journal, so Sandage's paper would appear in the May 15 issue.  Tom Kinman \cite{K65} along with Lynds and Villere \cite{LV65} argued that most of Sandage's Blue Stellar Objects are just that, blue stellar objects, while Fritz Zwicky \cite{Z65} pointed out that he had previously called attention to this phenomena which he had called ``compact galaxies'' and he later accused Sandage of ``one of the most outstanding feats of plagerism'' \cite{Z71}.

As it turned out, most of Sandage's quasi stellar galaxies were indeed galactic and there are only some ten times more radio quiet quasars than radio loud quasars.   But, it has now been more than half a century since we have divided quasars into the two classes of radio loud and radio quiet quasars, and it is still has not been clear if there are two distinct populations or  whether the radio loud population is merely the extreme end of a continuous distribution of radio luminosity. Proponents of each interpretation claim that the other interpretation is due to selection effects.  Many of the previous investigations designed to distinguish between radio loud and radio quiet quasars were limited by contamination from low luminosity AGN with absolute optical magnitudes greater than -23, by the use of samples based on radio rather than optical selection criteria, and by inadequate sensitivity to detect radio emission from most of the radio quiet population.

In an attempt to overcome these limitations, Kimball et al. \cite{K11} used the Jansky Very Large Array to observe 179 quasars selected from the SDSS. All of the  quasars were within the redshift range 0.2 to 0.3 and were brighter than $M_{i}$ = -23 so were legitimate quasars.  The radio observations were made at 6 GHs reaching an rms noise of $\sim$6 $\mu$Jy. All but about 6 quasars were detected as radio sources with an observed radio luminosity sharply peaked between $10^{22}$ and $10^{23}$ Watts/Hz characteristic of the radio luminosity typically observed from star forming galaxies.  About ten percent of the SDSS sample were found to be strong radio sources with radio luminosities ranging up to $10^{27}$ W/Hz.  Kimball et al.  concluded that the radio emission from radio quiet quasars is due to star formation in the host galaxy and is not directly related to the SMBH assumed to drive the large optical luminosity of quasars.

\section{Summary and Conclusions}

Quasars and AGN are now part of the normal lexicon of astronomy. Although the 1962 occultations of 3C 273 were crucial in bringing 3C 273 to the attention of Maarten Schmidt, the occultation position used by Schmidt for his Palomar observations was in error by about 15 seconds of arc.  A more accurate position had been established from OVRO interferometer observations but apparently was not accurately conveyed to Schmidt.

Prior to the Parkes occulation, 3C 273 was not on anyone's list of potentially interesting sources to be further studied, as it was already known not to be very small and thus presumably (and correctly) not very distant.  Schmidt's December 1962  spectroscopic observations  of 3C 273 were made  only because he had received the occultation position from Bolton.  However, since the correct radio position was not available until after Schmidt had recognized that the optical counterpart was not galactic, but had a redshift, z=0.16, the occultation positions did not directly result in the identification of the optical counterparts to the radio positions. Rather, not having an accurate radio position, Schmidt \cite{S83,S11}  believing that the correct optical counterpart was the jet, apparently took a spectrum of the quasar itself only to ``get it out of the way.''  Ironically, it was possible for Schmidt to recognize the ``large'' redshift of 3C 273, because, unlike for 3C 48, the redshift of 3C 273 was ``small,'' and the familiar Balmer spectrum was not shifted out of the optical window.

The earlier identification of 3C 48, and later the two other apparently stellar counterparts of 3C 196 and 3C 286 played no role in the recognition of quasars.  The 1960 rejection by Greenstein, Bowen, and ultimately also by Bolton of the suggested 3C 48 redshift  was explained by a 3 to 4 angstrom discrepancy in the derived rest wavelengths, an amount that might easily have been understood as the result of Doppler motions, but more likely was their unwillingness to accept the implications of the extraordinary luminosity  and small dimensions implied by the optical variability.

Following the recognition of the 3C 48 and 3C 273 redshifts, new identifications quickly led to redshifts  up to z $\sim$2.  But, contrary to hopes, because they are not standard candles, quasars have had little impact to classical cosmology and the quest to define $H_0$ or $q_0$.  However, discovery of quasars has led to the recognition that they are driven by a central engine, widely believed to be due to accretion onto a SMBH of $\sim 10^9$ or more solar masses.

Over the half century since their discovery, quasars have had a profound impact to the sociology of astronomy and astronomers as well as to astrophysics and cosmology.  For the first decade or more, there was vocal minority that continued to advocate a non cosmological interpretation of quasar redshifts.  The controversy was intense and sometimes very personal.  The arguments never died, until their proponents died.

\section{Unanswered Questions}

$\bullet$  Why did  Greenstein, Matthews, Sandage and others wait two years to submit papers interpreting  3C 48 as a galactic star. Why did John Bolton wait nearly thirty years before going public with his claim to the discovery of the 3C 48 redshift?

$\bullet$  Why did it take so long to identify 3C 273 with such a bright magnitude 13 optical counterpart, even though by the time of the 1962 Parkes occultation, at least three other sources, 3C 48, 3C 147, and 3C 196 had already been identified with apparent stellar counterparts, and a Caltech interferometer position was apparently known to an accuracy of 5 or 10 arcseconds  by the time of Schmidt's May 1962 Palomar observations.

$\bullet$ Why was the extragalactic nature of 3C 48 not accepted earlier as suggested by John Bolton in November, 1960 and considered by Greenstein in his unpublished manuscript submitted to the Astrophysical Journal?   Why was a 3 or 4 Angstrom discrepancy among the derived rest wavelengths considered fatal to the z=0.37 interpretation when just a few months earlier Minkowski had gone out on a limb with the z=0.46 redshift based on only one spectral line and a comparable radio luminosity?  Why were Bolton and Schmidt ignored in the 1963 publications by Greenstein and Matthews (1963) and by \cite{MS63}?

$\bullet$  Why did it take Maarten Schmidt six weeks to decipher the simple Balmer spectrum of 3C 273?  Was it possibly due to the intervention of  holiday activities?  Or, perhaps, it was the receipt of Bolton's January 26 letter, with improved radio positions, that inspired Schmidt to re-inspect and finally understand the 3C 273 spectrum six weeks after he had observed it.

$\bullet$  Although all quasars presumably contain a SMBH which drives the excessive luminosity, why are only about 10\% of optically selected quasars also strong radio sources?

\section{Acknowledgment}  The National Radio Astronomy Observatory is operated by Associated Universities Inc. under Cooperative Agreement with the National Science Foundation.   I am grateful to many colleagues, especially Maarten Schmidt, Tom Matthews, Jesse Greenstein, Allan Sandage, Cyril Hazard, and John Bolton who have shared with me their recollections of the events surrounding the discovery of quasars.   Miller Goss and Ron Ekers brought the 1960 and 1961 letters from John Bolton to Joe Pawsey to my attention.  Dave Jauncey provided valuable comments, and John Faulkner alerted me to the 1981 recollection by Hoyle of Bolton's 1960 claim that 3C 48 was extragalactic with z=0.37.

Instrumentation and Methods for Astrophysics (astro-ph.IM); Astrophysics of Galaxies (astro-ph.GA); High Energy Astrophysical Phenomena (astro-ph.HE)

The Road to Quasars


Graphane Nanoribbons: A Theoretical Study

Abstract: In this study, we investigate the electronic and magnetic properties of graphane nanoribbons. We find that zigzag and armchair graphane nanoribbons with H-passivated edges are nonmagnetic semiconductors. While bare armchair ribbons are also nonmagnetic, adjacent dangling bonds of bare zigzag ribbons have antiferromagnetic ordering at the same edge. Band gaps of the H-passivated zigzag and armchair nanoribbons exponentially depend on their width. Detailed analysis of adsorption of C, O, Si, Pt, Ti, V and Fe atoms on the graphane ribbon surface reveal that functionalization of graphane ribbons is possible via these adatoms. It is found that C, O, V and Pt atoms have tendency to replace H atoms of graphane. We showed that significant spin polarizations in graphane can be achieved through creation of domains of H-vacancies and CH-divacancies.

\section{Introduction}\label{Introduction}

After the synthesis of two dimensional (2D)
graphene\cite{novo1,novo2,zhang,berger,kats,geim}, its
nanoribbons (NRs) have been a subject of interest, both
experimentally\cite{han,li} and
theoretically.\cite{lee,barone,abanin,cresti,ezawa0,gunlycke2,hod1} Advances in experimental techniques are paving the way for
integrating the exceptional electrical, optical and magnetic
functionalities of these nanometer-sized materials into future
electronic technology. In addition to numerous theoretical and experimental works on
graphene and its NRs, research efforts have been also
devoted to synthesize its various derivatives. Notably,
the synthesis of a 2D hydrocarbon in honeycomb
structure,\cite{elias} namely graphane, followed its prediction through theoretical works.\cite{sofo,boukhvalov}
Detailed analysis of hydrogenation processes of graphene leading
to graphane and the existence of hydrogen (H) frustrations were discussed very
recently.\cite{Graphene to graphane} Furthermore, it was reported
that graphane NRs have more favorable formation
energies than experimentally available graphene ribbons.
\cite{yafei} Recently, we reported the possibility of obtaining
magnetization through dehydrogenation of domains on 2D
graphane.\cite{hasan-can} Modification of electronic structure of graphane by introducing either a hydroxyl group or a H-vacancy
was also investigated with GW self-energy calculations.
\cite{lebeque} Although there are several recent studies on
graphane, electrical and magnetic properties of its NRs
have remained unexplored.

Recent advances in graphane have motivated us to explore the
properties of zigzag and armchair graphane NRs. In this
paper, using first-principles plane wave method within the
density functional theory (DFT) we investigated the electronic
and magnetic properties of bare and H- passivated graphane
NRs. We also explored the effects of specific
imperfections on these properties. These are various foreign
atoms adsorbed on the surfaces of graphane NRs, vacancies and
edge roughness. We found that these imperfections can attribute
interesting functionalities by modifying the electronic and
magnetic properties of graphane NRs.

\begin{figure}
\includegraphics[width=8.5cm]{fig-1.png}
\caption{(Color online) (a) Top and perspective view of the atomic
structure of infinite 2D graphane sheet having honeycomb structure.
Two sublattices of graphane are indicated by A and B. Black (dark)
and blue (light) balls are for carbon and hydrogen atoms,
respectively. (b) Charge density contour plots of diamond and
graphane are shown on a plane passing through C-C-C-C and H-C-C-H
bonds, respectively. The tetrahedral angle of the diamond $\theta_C
= 109.47^o$. Arrows indicate the direction of increasing charge
density. The calculated values of $\theta_C$ and $\theta_H$, namely
C-C-C and H-C-C bond angles in graphane respectively, are given in
Table~\ref{tab:graphene-graphane}. Contour spacings are 0.0286
e/\AA$^3$. (c) The LDA energy band structure where the orbital
character of specific bands is also given. The band gap is shaded
yellow/gray.} \label{2d-graphane}
\end{figure}

\section{Methods}\label{Methods}

We perform first-principles calculations\cite{kresse} within DFT using
projector augmented wave (PAW) potentials\cite{paw} and
approximate exchange-correlation functional by spin polarized
local density approximation\cite{lda} (LDA). Kinetic energy
cutoff $ \hbar^2 |\mathbf{k}+\mathbf{G}|^2 / 2m $ for plane-wave
basis set is taken as 500 eV. In the
self-consistent potential and total energy calculations of 2D graphane a set
of (35x35x1) \textbf{k}-point sampling is used for Brillouin
zone (BZ) integration. The convergence
criterion of self consistent calculations for ionic relaxations
is $10^{-5}$ eV between two consecutive steps. By using the
conjugate gradient method,
all atomic positions and unit cells are optimized until the
atomic forces are less than 0.03 eV/\AA. Pressures on the
lattice unit cell are decreased to values less than 1.0 kBar.

In order to correct the energy bands and band gap values obtained
by LDA, frequency-dependent GW$_0$  calculations are carried
out.\cite{gw} Screened Coulomb potential W, is kept fixed
at initial DFT value W$_0$ and Green's function G, is
iterated five times. GW$_0$ self-energy calculations
are carried out with $15 $ \AA ~ vacuum  spacing,
default cut-off potential, 160 bands and 64 grid points.

Graphane NRs are treated by the supercell geometry
within periodic boundary conditions. Vacuum spacing of at least $15 $ \AA~ is placed between adjacent graphane NRs to hinder the  interactions. In specific cases double unit cell is
used in our calculations to allow the possible
antiferromagnetic (AFM) orderings along the ribbon edges.
Reciprocal space integrations are carried out with 1x1x15
Monkhorst-Pack \textbf{k}-point grids.

\section{Two Dimensional Graphane}\label{2D}

For a better understanding of graphane NRs we start with 2D graphane
which is derived from the hydrogenation of graphene, where each
carbon atom is saturated by a single hydrogen atom. Accordingly, the
primitive cell of graphane consists of two carbon and two hydrogen
atoms. Chair like configuration of infinite sheet of 2D graphane is
formed, whereby each carbon atoms of A- and B-type sublattices are
saturated by hydrogen atoms from above and below, respectively, as
described in Fig.~\ref{2d-graphane}(a). This is known as the most
favorable and stable hydrocarbon conformation.\cite{sofo, hasan-can}
The planar graphene honeycomb structure, which is stabilized by
planar $sp^2$-hybrid orbitals and $\pi$-bonding between adjacent
perpendicular $p_z$-orbitals is puckered (buckled) as a result of
the adsorption of H atoms, whereby a single C atom is bound to three
adjacent C atoms and a single H atom through 'tetrahedrally
coordinated' $sp^{3}$-like hybrid orbitals. The buckling, $\delta$,
i.e. the perpendicular distance between A-type and B-type carbon
sublattices, is calculated to be 0.45 \AA. As a result, the lattice
constant of graphene increases from 2.46~\AA~ to 2.51~\AA~ and hence
the C-C and C-H bond lengths become 1.52 and 1.12~\AA,
respectively, as indicated in Fig.~\ref{2d-graphane}(a). The charge
density contour plots given in Fig.~\ref{2d-graphane}(b) indicate
that the charge distribution in graphane especially around the C-C
bond is reminiscent of that in tetrahedrally coordinated diamond. In
graphane, the angle $\theta_H$ between H-C and C-C bonds and
$\theta_C$ between adjacent C-C bonds are $107.35^o$ and $111.51^o$,
respectively. The mean value of these angles is equal to the
tetrahedral angle $\theta_C$ of diamond. The maximum of the C-C bond
charge in diamond and in graphane are found to be $\sim 0.282$ and
$\sim 0.286$ e/\AA$^3$, respectively.

\begin{table*}
\caption{Comparison of the calculated quantities of graphene and
graphane. Lattice constant, $a$; C-C bond distance, $d_{C-C}$; C-H
bond distance, $d_{C-H}$; the buckling, $\delta$ [see
Fig.~\ref{2d-graphane}(a)]; angle between adjacent C-C bonds,
$\theta_{C}$ [see Fig.~\ref{2d-graphane}(b)]; angle between adjacent
C-H and C-C bonds, $\theta_{H}$;  energy band gap calculated by LDA,
$E_{g}^{LDA}$; energy band gap corrected by GW$_0$, $E_{g}^{GW_0}$;
cohesive energy $E_{c}^{nm}$, ($E_{c}^{m}$) obtained with respect to
nonmagnetic (magnetic) free atom energies; the C-H bond energy,
$E_{C-H}^{nm}$ ($E_{C-H}^{m}$) obtained with respect to nonmagnetic
(magnetic) free atom energies; photoelectric threshold (work
function), $\Phi$; in-plane stiffness, $C$ and Poison ratio, $\nu$.
} \label{tab:graphene-graphane}
\begin{center}
\begin{tabular}{cccccccccccccccc}
\hline  \hline
Material &       $a$    & $d_{C-C}$  & $d_{C-H}$ & $\delta$ & $\theta_{C}$ & $\theta_{H}$ &$E_{g}^{LDA}$ &  $E_{g}^{GW_0}$  &    $E_{c}^{nm}$    &   $E_{c}^{m}$    &     $E_{C-H}^{nm}$     & $E_{C-H}^{m}$     & $\Phi$     & $C$     & $\nu$\\
(1x1) unit cell &(\AA)     & (\AA)   & (\AA)   &    (\AA) &  (deg) & (deg) &($eV$)  &  ($eV$)   &     ($eV$)           &  ($eV$)             &      ( $eV$ )          &    ($eV$)             &   ($eV$)    &  ($J/m^{2}$)   &\\
\hline \hline
Graphene & 2.46 & 1.42 & - & - & 120 & - &  0.00 & 0.00 &20.16 & 17.87& -& - &4.77 & 335\cite{lee-science} &0.16\\
\hline
Graphane & 2.51 & 1.52& 1.12 & 0.45 & 111.51 &107.35 &3.42& 5.97 & 27.65 & 23.57 &3.74 &2.84 & 4.97& 243\cite{topsakal-graphane} &0.07\\
\hline  \hline
\end{tabular}
\end{center}
\end{table*}

The cohesive energy of graphane (per unit cell) relative to free C
and H atoms is obtained from $E_c = 2\times E_T^H+2\times
E_T^C-E_T^{Gra}$ , where $E_T^H$ and $E_T^C$ are the total energies
of single free C, and free H, whereas $E_T^{Gra}$ is the total
energy of graphane. The cohesive energy depends whether the energies
of magnetic or nonmagnetic (NM) states of free C and free H atoms
are considered. Here the cohesive energy per unit cell of graphane
is calculated to be 23.57 (27.65) eV by considering the magnetic
(nonmagnetic) states of free atoms. As for the average C-H bond
energy, one can use the formula $E_{C-H}=
(-E_T^{Gra}+E_T^{Gre}+2\times E_T^H)/2$ where $E_T^{Gre}$ is the
total energy of graphene. The calculated value is 2.84 (3.74 for
nonmagnetic case). We reported\cite{hasan-can} that the desorption
of a single H atom from the graphane is an endothermic reaction with
an energy of 4.79 eV. When compared with the energy of C-H bond, the
desorption energy of single H is larger, since the single H removal
creates reconstruction of nearby atoms which further reduces the
total energy of the system.

Due to $sp^{3}$-saturation of C atoms, 2D graphane is a NM
semiconductor with a direct band gap of 3.42 eV as shown in
Fig.~\ref{2d-graphane}(c). However, this band gap, which is
underestimated within DFT, is corrected by GW$_{0}$ approximation to
become 5.97 eV.\cite{hasan-can, lebeque} While the top of the
valence band originates from $p_{x}$- and $p_{y}$-orbitals of C
atoms, the bottom of the conduction band has C-$p_{z}$ orbital
character. Earlier, we also calculate\cite{hasan-can} the phonon
modes of infinite 2D graphane, which yield all real frequencies in
BZ. Having all frequencies positive indicate the stability of
graphane structure.

An exceptional feature of graphane is related with its
interesting charging configuration. While graphene has a
covalent bonding, upon the saturation of each C atom by a
single H atom, a charge of $\delta Q\sim$ 0.1 electrons is
transferred from H to C. At the end, the negatively charged
bilayer of carbon atoms becomes sandwiched between positively
charged H layers. Graphane having this charging has the photoelectric threshold (work
function) calculated to be $\Phi$=4.97 eV, which is 0.2 eV
larger than that of graphene. In Table~\ref{tab:graphene-graphane}, we presented the calculated values related with structure, energy bands, photoelectric threshold etc. of graphane and graphene for the sake of comparison. To be complete Table~\ref{tab:graphene-graphane} also includes calculated elastic constants such as in-plane stiffness\cite{brandbyge} $C ~[=\frac{1}{A}( \frac{\partial^2 E_T}{\partial \epsilon^2})$, where $E_T$ is the total energy of the system under strain, $\epsilon$ is the uniaxial strain and $A$ is the area of the unit cell] and Poisson's ratio\cite{hasan-ansiklopedi} $\nu$.

\section{Graphane Nanoribbons}\label{Graphane NRs}

The nanoribbons cut from a 2D graphane are structures providing
important features for various technological applications. Two major
families of graphane NRs are distinguished depending on their
orientations; namely armchair, and zigzag graphane NRs. Apart from
the orientation, graphane NRs are characterized by their widths. For
armchair graphane NRs, the width, $N$, is defined by the number of
C-C dimers in the unit cell which are parallel to the axis of the
nanoribbon; for zigzag graphane NRs, $N$ denotes the number of
zigzag carbon chains along the nanoribbon axis. The electronic
structure changes depending on whether the dangling bonds of carbon
atoms at the edge are free (bare) or passivated by single H atoms.
In Fig.~\ref{ribbon1}(a) and ~\ref{ribbon1}(b) we show the atomic
and electronic structure of armchair and zigzag graphane NRs with
$N$=11 and $N$=6, respectively.

Bare armchair graphane NR is a NM indirect band gap semiconductor,
since the electrons at the dangling orbitals are paired. The band
gap of 2D graphane is reduced due to the bands of edge states of
dangling bonds occurring in the band gap. The charge density
analysis of these edge state bands presented in
Fig.~\ref{ribbon1}(a) clearly shows that they are localized at the
edges of the NR. Upon hydrogen passivation of the threefold
coordinated edge atoms, the edge state bands in the gap disappear
and the direct band gap opens up.

Similar to armchair graphane NRs, the dangling bonds of the
threefold coordinated carbon atoms at both edges of zigzag graphane
NRs give rise to edge states in the band gap. Owing to relatively
large distance between adjacent dangling bonds and hence their
relatively smaller coupling, they have a minute dispersion. Charge
density isosurfaces of these edge states confirm their localization
at the edges. Also, owing to their weak coupling the spins of
adjacent dangling bonds cannot be paired. Consequently, each carbon
atom at the zigzag edge attains a magnetic moment of 1  Bohr
magneton ($\mu_B$). The ordering of these magnetic moments is
however dramatically different from that of the zigzag graphene NRs.
The significant coupling between two spin states at different edges
gives rise to AFM ordering along the graphene edge.

\begin{figure}
\includegraphics[width=8.5cm]{fig-2.png}
\caption{(Color online) (a) Atomic structure of bare armchair
graphane NR having $N$=11. The double unit cell of the nanoribbon is
delineated by red/dashed lines with the lattice constant $2a$.
Large/black and small/light blue balls indicate carbon and hydrogen
atoms. Energy band structure corresponding to the armchair graphane
NR and charge density of selected bands are shown in the panels at
the righthand side. (b) Atomic structure of bare zigzag graphane NR
having $N$=6 with double unit cell delineated by red/dashed lines
and with lattice constant $2a$. Energy band structure and isosurface
charge density of selected states corresponding to zigzag graphane
NR are indicated. Bands shown by red/dotted lines are derived from
edge states. Zero of energy is set to the Fermi level, shown by
dash-dotted lines, of the nanoribbons with H-passivated edges. In
spin polarized calculations double unit cell is used to allow
antiferromagnetic order along the edges.} \label{ribbon1}
\end{figure}

The physical origin of the magnetic interactions between the edges
of graphitic fragments was treated earlier.\cite{edge-states} For
zigzag graphene nanoribbon, collective alignment of the magnetic
moments along the same ribbon edge through the ferromagnetic (FM)
interaction, but AFM ordering between opposite edges is attributed
to magnetic tail interactions. The magnetic ordering of a wide
zigzag graphane NRs ($\sim$ 2 nm) are examined in a supercell
comprising two unit cells, whereby antiferromagnetic ordering
between adjacent dangling bonds at the same edge is allowed. We
consider four different orderings, such as AFM-AFM (where,
respectively, the ordering of magnetic moments at the atoms located
at different and same edge are antiferromagnetic), FM-AFM, AFM-FM
and FM-FM as described in Fig.~\ref{ribbon2}. Here the difference
charge density is defined as the $\Delta
\rho=\rho^{(\uparrow)}-\rho^{(\downarrow)}$ where $\rho^{\uparrow
(\downarrow )}$  is spin up (down) charge. Since the magnetic
interaction between two edges vanishes, AFM-AFM and FM-AFM
orderings have the same total energy corresponding to the ground
state. These orderings are 125 meV energetically more favorable than
the AFM-FM and FM-FM orderings. From this analysis it is revealed
that AFM spin alignment between adjacent atoms at the same edge is
preferred by the zigzag graphane NRs. FM ordering between two spins
has 62.5 meV higher energy. For zigzag graphane NR with $N$=6, the
magnetic coupling between two edges is negligible and hence either
AFM-AFM or FM-AFM orderings have the same energy and they are ground
states. However, when the width of the nanoribbon is less than
12~\AA~the degeneracy of AFM-AFM and FM-AFM states is broken. The
bare zigzag graphane NR is an indirect and antiferromagnetic
semiconductor with an indirect band gap relatively smaller than that
of 2D graphane. However, it becomes NM, direct band gap
semiconductor upon passivation of the dangling bonds. Also the
magnetic edge states disappear and band gap becomes larger.

The value of the band gap of graphane NRs can differ from that of parent 2D structure due to a combined effects of quantum confinement, edge state bands as well as folding in the direction perpendicular to the NR axis. We note that the band gap of graphane NR depends on its width given by N. However, the quantum confinement effect is not as emphasized as that found in graphene NR. In Fig.~\ref{ribbon3}, H-passivated armchair and zigzag graphane NRs display a quantum confinement effect, namely the band gap reduction with increasing width or $N$. This behavior is fitted to an expression as:

\begin{equation}
E_{gap}(N)=3.42~eV +\alpha\exp(-N\beta)
\end{equation}
Here $\alpha$ and $\beta$ are fitting parameters. The values of
$\alpha$ and $\beta$ are found to be 1.18 eV and 0.19 for zigzag
graphane NRs (2.15 eV and 0.14 for armchair graphane NRs). The band
gaps of both types of graphane NRs go to 3.42 eV as $N \rightarrow$
$\infty$. Considering the GW$_0$ corrected value of the band gap,
5.97 eV, the above fitting can be expressed as
$E_{gap}(N)=5.97+\alpha\exp(-N\beta)$ assuming that the same scissor
operation is applicable for graphane NRs. Apparently, the band gap
values of graphane NRs calculated using DFT are underestimated and
hence their values are expected to occur $\sim$ 2.5 eV larger than
presented in Fig.~\ref{ribbon3}. Unfortunately, GW$_{0}$ self-energy
corrections cannot be performed for graphane NRs due to large
computational time. We note that the variation of band gap is
relevant for $N<22$ for zigzag graphane NRs and $N<30$ for armchair
graphane NRs. In Ref.[\onlinecite{yafei}] the band gaps of hydrogen
saturated graphane NRs were calculated for $6 \leq N \leq 15$ using
Generalized Gradient Approximation (GGA). In the present study band
gaps are calculated for $10 \leq N \leq 40$ using LDA. The
hydrogenated graphane NRs included in both studies, namely those $10
\leq N \leq 15$, are in good agreement. Largest band gap difference
of 0.1 eV found between two studies is attributed to different
approximations (i.e LDA versus GGA) in representing
exchange-correlation potential.

\begin{figure}
\includegraphics[width=8.5cm]{fig-3.png}
\caption{(Color online) Total energies of possible magnetic
orderings at the edges of bare zigzag graphane nanoribbons.
Calculations are performed in double unit cell delineated by
red/dashed lines. Spin up and spin down states are shown by
green/dark and grey/light isosurfaces of the difference charge
density, $\Delta \rho$.}

\label{ribbon2}
\end{figure}

\section{Functionalization of Graphane NRs by Adatoms}\label{Adatoms}

Adsorption of adatoms is widely used and an efficient way to provide
new functionalities to structures in nanoscale
applications.\cite{hasan-prb,can1,cohen-adsorp} Adsorption of C, O,
Si, Ge, Pt, V, Fe and Ti adatoms on graphane nanoribbons is
investigated by using (1x1x6) supercell of H-passivated armchair
graphane NR having $N$=8. Adsorption geometry is determined by
calculating the lowest binding energy corresponding to the optimized
structures. To this end, the binding energies of adatoms at four
initial sites described in Fig.~\ref{adatom}: These are edge site E;
hollow site HO, top site of the carbon TC and top site of the
hydrogen atom TH in the middle of NR. The stable binding sites are
determined upon structure relaxation. In the same figure the
positions of adatoms after relaxation are also shown. Binding
energies ($E_b$) of these atoms are calculated by using the
expression
\begin{equation}\label{equ:binding}
E_b={-E_{T}[NR+A]} + E_{T}[NR] + E_{T}[A]
\end{equation}
in terms of the total energies of the optimized structure of
graphane NR with adatom, $E_{ T}[NR+A]$; bare NR, $E_{T}[NR]$ and
free adatom, $E_{T}[A]$. Owing to the large graphane nanoribbon unit
cell, \textbf{k}-point sampling is done only along the axis of the
ribbon. All the energies are calculated in the same unit cell and
obtained from the lowest ground-state total energies (either
magnetic or nonmagnetic). Initially, at least 2~\AA~ distance
between adatom and outermost graphane atom is provided to find the
relaxation site of adatom on the ribbon.

An atom adsorbed on graphane NR may give rise to resonance states in
the valence and conduction bands, and also localized states in the
band gap. Owing to the periodic boundary conditions, the energy
levels associated with these adatoms may form energy bands. Since we
consider a supercell consisting of six unit cell of H-passivated
armchair graphane NR, the adatom-adatom distance is large and
prevents coupling between adjacent adatoms. In this respect, the
flat bands associated with adatoms can mimic either resonance or
localized states relevant for dilute doping.

\begin{figure}
\includegraphics[width=8.5cm]{fig-4.png}
\caption{(Color online) Variation of the energy band gap of
H-passivated zigzag and armchair NRs of graphane as a function of
width given by $N$. The variation of the band gap with $N$ is fitted
to the curve given by continuous line. (See text)}

\label{ribbon3}
\end{figure}

\begin{figure}
\includegraphics[width=8.5cm]{fig-5.png}
\caption{(Color online) Schematic representations of possible
positions of adatoms on a large H-passivated armchair graphane
nanoribbon. Positions of adatoms obtained after the structure
optimization through minimization of total energy and forces
exerting on the atoms are also shown.}

\label{adatom}
\end{figure}

In Table~\ref{tab:adsorp}, we present nearest neighbor distances
between adatom and H ($d_{H}$) and C ($d_{C}$), binding energy
($E_{b}$), total magnetic moment of the system ($\mu_{T}$), charge
transfer between adatom and adsorbate ($\rho*$). Since the inversion
symmetry of graphane can be broken by the adsorption of an atom, a
net electric-dipole moment can be induced thereof. Thus electric
dipole moment values calculated in the direction normal to graphane
NR surface p, induced by adatom are also listed in
Table~\ref{tab:adsorp}. Photoelectric threshold of the
adatom+graphane system perpendicular to the plane of NR, $\Phi$, is
calculated by the difference of the electrostatic potential at
distances where the gradient of it is negligible and the Fermi level
of the adsorbed system. It should be noted that both p and $\Phi$ is
relevant only for uniform coverage of graphane. For a dilute
impurity, these calculated values converge to that of armchair
graphane. When the width of armchair graphane NR is large enough,
one expects the photoelectric threshold, $\Phi$, and dipole moment,
p, to converge to 4.97 eV and 0 e\AA, respectively.\cite{hasan-can}

\begin{table*}
\caption{Summary of the calculated quantities for adatoms adsorbed
on a H-passivated armchair graphane NR. The first and second lines
in each row associated with a given adatom adsorbed to edge site and
the sites near the center of the graphane NRs, respectively. $d_H$,
the adatom-nearest hydrogen distance; $d_C$, the nearest adatom
carbon distance; $E_b$, adatom binding energy; $\mu_{T}$, magnetic
moment per supercell; $\rho^{*}$, excess charge on the adatom (where
negative sign indicates excess electrons); $\Phi$, photoelectric
threshold (work function); p, dipole moment calculated in the
direction normal to graphane NR surface; $E_i$ energies of localized
states associated with adatoms. Localized states are measured from
the top of the valence bands in eV. The occupied ones are indicated
by bold numerals and their spin alignments are denoted by either
$\uparrow$ or $\downarrow$. Up to first seven adatom states of
$E_{i}$ are shown. } \label{tab:adsorp}
\begin{center}
\begin{tabular}{ccccccccccc}
\hline  \hline
$Atom$&$Site$ &$d_H$ & $d_C$ & $E_b$& $\mu_{T}$ & $\rho^{*}$ & $\Phi$& p & $E_{i}$  &\\
&&(\AA) & (\AA) & ($eV$)& $(\mu_{B})$ & (e)& $(eV)$& (e$\times$\AA)& $\uparrow$: Spin-up,  $\downarrow$: Spin-down States&\\
\hline \hline
\textbf{C}  & $E \rightarrow E1$   & 1.27 & 2.05 & 1.10 & 2.00& -0.28& 3.85&0.02& \textbf{0.45$\uparrow$}, \textbf{0.88$\uparrow$}, 2.17$\downarrow$, 2.88$\downarrow$, 3.53$\uparrow$,   3.64$\downarrow$ \\
\textbf{}   & $TH \rightarrow M1$  & 1.12 & 1.45 & 3.99 & 0.00& -0.18& 4.02& 0.34 &\textbf{0.72, 1.52}  \\

\hline
\textbf{O}  & $E \rightarrow E4$   & 0.98 & 1.42 & 5.67 & 0.00& -1.23& 4.63&0.19& \textbf{-0.87, -0.50, -0.22, 0.00}, 4.80  \\
\textbf{}   & $TH \rightarrow M1$  & 0.98 & 1.43 & 5.89 & 0.00& -1.27& 4.59& 0.11& \textbf{-1.16, -1.02, -0.66}, 4.08 \\

\hline
\textbf{Si}  & $E \rightarrow E3$   & 1.88 & 2.65 & 0.54 & 2.00& 0.03& 3.42&0.03& \textbf{1.20$\uparrow$, 1.34$\uparrow$}, 2.36$\uparrow$, 2.43$\downarrow$, 2.56$\downarrow$, 2.95$\downarrow$\\
\textbf{}    & $TH \rightarrow TH$  & 1.84 & 3.00 & 0.68 & 2.00& -0.02& 3.59& 0.42& \textbf{1.23$\uparrow$, 1.24$\uparrow$}, 1.93$\uparrow$, 2.42$\downarrow$, 2.45$\downarrow$, 2.46$\downarrow$ \\

\hline
\textbf{Ti}  & $E \rightarrow E2$   & 1.97 & 2.16 & 0.81 & 2.00& 0.19& 2.65&0.27& \textbf{1.77$\uparrow$, \textbf{2.30$\downarrow$}, 2.34$\uparrow$, 2.43$\uparrow$}, 2.81$\uparrow$, 2.89$\uparrow$, 3.11$\uparrow$\\
\textbf{}    & $TC \rightarrow TC$  & 1.95 & 2.52 & 0.89 & 3.91& 0.31& 2.16& 0.92& \textbf{2.18$\uparrow$, 2.48$\uparrow$, 2.49$\uparrow$, 2.95$\uparrow$}, 3.18$\downarrow$, 3.27$\uparrow$, 3.28$\uparrow$ \\

\hline
\textbf{V}  & $E \rightarrow E1$   & 1.99 & 2.40 & 0.54 & 5.00& 0.14& 2.56&0.22& \textbf{1.49$\uparrow$, 2.10$\uparrow$, 2.38$\uparrow$, 2.43$\uparrow$, 2.50$\uparrow$}, 2.70$\uparrow$, 3.16$\downarrow$\\
\textbf{}    & $TC \rightarrow M3$  & 1.71 & 2.07 & 1.24 & 3.00& 0.93& 3.13& 0.43& \textbf{0.10$\uparrow$, 0.49$\downarrow$, 1.45$\uparrow$, 1.60$\uparrow$, 1.84$\uparrow$}, 2.20$\uparrow$,  3.03$\downarrow$\\

\hline
\textbf{Fe}  & $E \rightarrow E2$   & 1.90 & 2.15 & 0.57 & 4.00& 0.05& 3.21&0.36& \textbf{-0.42$\uparrow$, -0.27$\uparrow$, -0.08$\uparrow$, -0.03$\uparrow$, 0.03$\uparrow$}, 1.48$\downarrow$, 1.59$\uparrow$ \\
\textbf{}   & $HO \rightarrow HO$  & 1.95 & 2.83 & 0.65 & 4.00& 0.07& 2.98& 0.53& \textbf{\textbf{-1.15$\downarrow$}, -0.75$\uparrow$, -0.57$\uparrow$, -0.41$\uparrow$, -0.34$\uparrow$, 0.09$\uparrow$ }, 1.39$\uparrow$\\

\hline
\textbf{Ge}  & $E \rightarrow E3$   & 1.96 & 2.71 & 0.48 & 2.00& 0.00& 3.38&0.03& \textbf{1.27$\uparrow$, 1.38$\uparrow$}, 2.25$\uparrow$, 2.41$\downarrow$, 2.52$\downarrow$, 2.81$\downarrow$\\
\textbf{}    & $TH \rightarrow TH$  & 1.91 & 3.06 & 0.62 & 2.00& -0.02&3.55 & 0.36& \textbf{1.30$\uparrow$, 1.31$\uparrow$}, 1.89$\uparrow$, 2.38$\uparrow$, 2.42$\uparrow$, 2.43$\uparrow$ \\

\hline
\textbf{Pt}  & $E \rightarrow E5$   & 1.56 & 2.00 & 3.79 & 0.00& -0.13& 4.57&-0.09& \textbf{-0.78, -0.36, -0.13}, 0.13, 2.76 \\
\textbf{}   & $TH \rightarrow M2$  & 1.54 & 2.04 & 3.50 & 0.00& -0.18& 4.30& 0.34&  \textbf{-0.36, -0.21, 0.38}, 2.58, 3.43 \\

\hline \hline

\end{tabular}
\end{center}
\end{table*}

We found that some atoms, such as C, O, V and Pt that are adsorbed
on the H-passivated armchair graphane NR surface show a tendency to
remove H atoms from graphane surface. Since the single atoms of C,
O, V and Pt have rather high binding energy, it seems possible to
create graphene domains on graphane by striping H atoms but covering
the domain by adsorbed atoms. Independent of their initial position,
C and O atoms are adsorbed to carbon atoms of graphane by replacing
H atoms. A typical binding configuration of C and O is indicated as
M1 in Fig.~\ref{adatom}. Adsorption of C and O atoms occur at the
middle of nanoribbon with strong binding energies of 5.89 and 3.99
eV, respectively. The resulting system is NM. Binding structures of
Pt and V, which also remove H atom from the host C atom of graphane,
are also shown in the same figure.

Although H has very low binding energy (and hence is not shown in Table~\ref{tab:adsorp}), it may be important for striping H atoms from graphane. Upon all
geometry relaxations from different initial conditions, single
adsorbed H atom prefers to bind on top of H atom of graphane and
forms H$_2$ molecule which is weakly bound to the nearest C atom.
The binding energy of H$_2$ molecule is calculated to be $\sim 36$
meV.

\section{Vacancies in Graphane NRs}\label{Vacancies}

In our recent study\cite{hasan-can}, we showed that a graphane,
which has a NM ground state, can be made magnetic simply by removing
hydrogen atoms from the uniform hydrogen layers adsorbed on its both
sides. It was found that the magnetization depends on whether the
defect region is one-sided or two-sided. It was also shown that in
certain circumstances remarkably large magnetic moments can be
attained in a small domain on the sheet of graphane. Experimentally,
removal of surface hydrogen atoms by using laser beam resonating
with surface-atom bond\cite{laser}, ionic vapor
application\cite{breaux} and subnanometer Pt clusters\cite{vajda}
are reported before. Since it is revealed that the double-sided
vacancies give higher magnetic moments we focus on these type of
vacancies on the nanoribbons.

In this section we consider six different types of vacancies in a
bare armchair graphane NR with $N$=15. Upon desorption of a single
hydrogen atom, as shown in Fig.~\ref{vacancy}, local bonding through
$sp^{3}$-hybrid orbital is retransformed into planar $sp^{2}$- and
perpendicular $p_{z}$-orbitals. At the vacancy site one electron
accommodated by the dangling $p_{z}$-orbital becomes unpaired and
hence contributes to the magnetization by 1 $\mu_{B}$. We also
calculated that a single H-vacancy located at the center of the NR
has 70 meV lower energy than that located near the edge of the
ribbon.

The situation is rather intriguing for the magnetization of large
domains of hydrogen vacancies. As in the case of single H-vacancy,
leading to transformation of $sp^{3}$-bonding into planar
$sp^{2}$-bonding, triangular $\Delta$, and rectangular $\square$
shaped large graphene islands can be created by removing H atoms
hence by creating H-vacancy domains. As shown in Fig.~\ref{vacancy},
hydrogen desorbed triangle-shaped island consisting of $n$ carbon
atom at its each edge is defined as $\Delta$n-type H-vacancy.
Therefore, a $\Delta$n-type vacancy domain can be obtained by
removing $n(n+1)/2$ H atom from top side  and $n(n-1)/2$ H atom from
down side  of the graphane nanoribbon. This means that $n(n+1)/2$
$p_{z}$-orbital electrons freed from A-type and $n(n-1)/2$
$p_{z}$-orbital electrons freed from B-type sublattice C atoms.
According to the Lieb's theorem \cite{Lieb} net magnetic moment of
such a system can be calculated from the difference between the
sublattice atoms. For triangle vacancies this difference is simply
equal to $n$ and thus the $\Delta$n-vacancy domain has net magnetic
moment with $n\mu_{B}$. However, the $\square$-type domains does not
yield large magnetic moments and have values between 0 and 2
$\mu_{B}$. It is seen that magnetic interactions in double-sided
triangular and square H-vacancy domains resulting in net magnetic
moments are relatively straightforward and are in good agreement
with Lieb's theorem.

We also consider CH- and C$_{2}$H$_{2}$-vacancies, the synthesis of
them are relatively easier than creating only H-vacancies. In our
study, these vacancies are first created in a NR by removing
involved atoms and subsequently their structures are optimized.
Calculated values of magnetic moments of the structure are presented
in Fig.~\ref{vacancy}. After the creation of a CH-vacancy on a
graphane nanoribbon, geometric structure rearranges itself via
Jahn-Teller-like distortion in the lattice. Charge density plot in
Fig.~\ref{vacancy} shows the unpaired up and down spin electron
states located on the region surrounding the vacancy. It is seen
that a bare armchair graphane NR, which is NM semiconductor when it
is defect free, attains permanent magnetic moment in the presence of
CH-vacancies. While CH-type vacancy makes small changes in the
geometric structure, Stone-Wales type deformation occurs after
formation of the C$_{2}$H$_{2}$-vacancy. Since electrons are paired
in C$_{2}$H$_{2}$-vacancy, the ground state is again NM and thus
there is no visible difference charge density in Fig.~\ref{vacancy}.
From the point view of Lieb's theorem, removal of a CH couple
results in an absence in the number of $p_{z}$-orbital electron
belonging to a A-type (or B-type) sublattice. Upon the creation of
C$_{2}$H$_{2}$-vacancy, both the A- and B-type sublattices lose one
electron and thus the difference is zero resulting in a NM ground
state. Our results are consistent with the recently reported values
for CH-vacancies in infinite sheets of graphane.\cite{ch-vacancy}

\begin{figure}
\includegraphics[width=8.5cm]{fig-6.png}
\caption{(Color online)  Atomic structures corresponding to
single-H, double-sided triangular shaped $\Delta$2 and $\Delta$3,
double-sided rectangular shaped $\square$2, CH and C$_{2}$H$_{2}$
vacancies and resulting difference charge density, $\Delta
\rho=\rho^{(\uparrow)}-\rho^{(\downarrow)}$, surfaces for a bare
armchair graphane NR. Large/black and small/light blue-gray balls
indicate C and H atoms, respectively. Only a small part which
includes vacancy region and its nearby atoms, of the armchair
graphane NR with $N$=15 is shown.} \label{vacancy}
\end{figure}

\section{Edge Roughness}\label{Roughness}

Earlier, it was shown that the edge roughness of graphene NRs can
affect their electronic\cite{Gunlycke, Yoon, Evaldsson} and
mechanical\cite{can-mechanical} properties. In particular, it was
revealed that periodically repeating edge profiles can be treated
within the superlattice structure, which result in confinement of
spin states in zigzag NRs.\cite{haldun-super} Band alignments of
these superlattice structures have been studied both experimentally
and theoretically.\cite{dai} Here we investigated the effect of
periodically repeating edge roughness of a bare zigzag graphane NR.
While one edge of NR is kept flat, periodically repeating
undulations are carved at the other edge. The same structure can be
viewed as periodically repeating heterostructures of wide ($N$=8)
and narrow ($N$=6) segments of zigzag graphane NRs. Atomic and
electronic band structure of this superlattice, which mimics the
edge roughness are shown in Fig.~\ref{roughness}. Here we note that
infinite bare zigzag graphane NR with $N$=8 has a band gap of 0.8
eV, while the band gap of infinite bare zigzag graphane NR with
$N$=6 is relatively wider and 0.9 eV. As pointed out in Sec. IV,
wide band gap of H-passivated graphane NRs is reduced to 0.8-0.9 eV
because of edge states of unpassivated dangling bonds of bare
graphane NRs, which appear in the band gap. These edge states in
Fig.~\ref{roughness} of the superlattice occur as several closely
lying flat bands below and above the Fermi level separated  by a
superlattice band gap of $\sim$0.8 eV. The character of various
bands are revealed by plotting the charge density isosurfaces of
various states. Charge density isosurfaces of the highest filled
spin up and spin down edge state indicate the localization at the flat
edge. As for the charge distribution of the lower lying valence
states is rather uniform along the narrow region. The states of the
flat band at the edge of the conduction band is confined in the wide
segment of the superlattice since they cannot find a matching state
in the narrow region. As a proof of the concept, it is shown that
electronic states can be confined at specific regions of periodic
edge roughness. Of course, confinement, superlattice band gap etc.
depend on the structural parameters of superlattices and need
detailed investigation.

\begin{figure}
\includegraphics[width=8.5cm]{fig-7.png}
\caption{(Color online) Energy band diagram and band projected
charge density isosurfaces of various states for bare zigzag
graphane NR including edge roughness. The band gap between edge
states are shaded yellow/gray. Zero of the band energy is set to the
Fermi level.} \label{roughness}
\end{figure}

\section{Discussions and Conclusions}\label{Discussion}

Research on recently synthesized graphane revealed interesting
electronic and magnetic properties of this two dimensional
honeycomb structure. Graphane attains a $sp^3$-like bonding
through covering of its both sides by hydrogen atom. This bonding
is rather different from the $sp^2$-bonding of graphene and
attributes a number of additional properties to graphane. For
example, in contrast to semimetallic graphene with linear band
crossing at the Fermi level, graphane is a nonmagnetic, wide band
semiconductor.

The armchair and zigzag graphane NRs of graphane display important
number of properties and hence constitute basic structures to
fabricate various devices. Both zigzag and armchair graphane NRs are
wide band gap semiconductor when their edges are passivated with
hydrogen. The band gaps vary exponentially with their widths. For
narrow graphane NRs the band gap is large due to quantum confinement
effect, but approaches to the band gap of graphane as the width
increases. Unpaired, unpassivated dangling bonds at the edges have 1
Bohr magneton magnetic moment and have antiferromagnetic coupling
with adjacent dangling bonds. These graphane NRs can be
functionalized to attain additional properties through H-passivation
of their edges, adatom adsorption, vacancy creation, edge profiling
and superlattice formation. In particular, graphane NRs attain
magnetic moment through the creation of H vacancy at the surfaces of
graphane. This property can be utilized to achieve interesting
functionalities through decoration and patterning of H-vacancies on
the graphane NRs. The possibility of generation of large magnetic
moments at small domains of H-vacancies makes graphane based
structures promising for data storage and nanospintronic
applications. These functionalities can be further extended by
adsorbing adatoms to carbon atoms deprived from hydrogen.

In conclusion, the present study demonstrates that graphane
NRs can be an important basic nanomaterial and presents interesting properties for future technological applications.

\begin{acknowledgments}
Computing resources used in this work were partly provided by the
National Center for High Performance Computing of Turkey (UYBHM)
under grant number 2-024-2007. This work is partially supported by
the project of The State Planning Organization (DPT) of Turkey and
by Academy of Science of Turkey (T\"{U}BA).
\end{acknowledgments}

Materials Science (cond-mat.mtrl-sci); Mesoscale and Nanoscale Physics (cond-mat.mes-hall); Computational Physics (physics.comp-ph)


Magnetoresistance of metallic perovskite oxide LaNiO$_{3-δ}$

Abstract: We report a study of the magnetoresistance (MR) of the metallic perovskite oxide LaNiO$_{3-\delta}$ as a function of the oxygen stoichiometry $\delta$ ($\delta \leq$ 0.14), magnetic field (H $\leq 6T$) and temperature (1.5K $\leq $ T $\leq $ 25K). We find a strong dependence of the nature of MR on the oxygen stoichiometry. The MR at low temperatures change from positive to negative as the sample becomes more oxygen deficient (i.e, $\delta$ increases). Some of the samples which are more resistive, show a resistivity minima at $T_{min}$ $\approx$ 20K. We find that in these samples the MR is positive at T > $T_{min}$ and negative for T < $T_{min}$. We conclude that in the absence of strong magnetic interaction, the negative MR in these oxides can arise from weak localisation effects.

\section{Introduction}
In recent years the electrical conduction in transition metal
perovskite oxides at low temperature has attracted considerable
research interest\cite{AKR}. Most of these oxides studied belong to
the ABO$_3$ class or its closely linked derivatives. In recent years
the observation of colossal magnetoresistance (CMR) in rare-earth
maganates belonging to the ABO$_3$ class has added yet another new
phenomena \cite{CHAA}. These oxides have chemical formula of the type
La$_{1-x}$(Sr/Ca)$_{x}$MnO$_{3}$ with $x\approx 0.2-0.4$. The
colossal {\bf negative} MR seen in these oxides have become a topic
of intense research activity. In these mixed valent manganates (doped
to create Mn$^{4+}$) the CMR is closely linked to the occurence of
strong ferromagnetism which arises from the Double Exchange (DE)
interactions. The additional important contribution to the CMR
phenomenon comes from the strong Jahn-Teller effect of the
Mn$^{3+}$ ions. In a closely related compound
La$_{1-x}$Sr$_{x}$CoO$_{3}$ (with Co$^{3+}$ and Co$^{4+}$) which is
also a ferromagnetic metallic oxide for $x\geq 0.2$ the MR was found
to be much smaller except in the composition range $x$ $\approx$ 0.2
where the oxide undergoes a composition driven metal-insulator
transition \cite{MAHI2}. The negative MR in this material for $x$
$\approx$ 0.2 is believed to be related to the formation of
ferromagnetic clusters. In both these oxides the electron is on the
verge of getting localized. On the other hand there are metallic
oxides also belonging to ABO$_3$ class like Sr$_{x}$NbO$_{3}$ (x$\geq$
0.75) \cite{SRNB} and Na$_{0.9}$WO$_{3}$ \cite{NG} with relatively
high conductivity which have positive MR at all temperatures. In
these oxides the MR has a quadratic dependence on the applied field
$H$ much like a band metal. The principal differences between the
two classes of perovskite oxides seems to be the extent of electron
localization and the presence of ferromagnetic interactions. This
leads us to the following interesting question: {\it If we start with
a metallic perovskite oxide showing positive MR and gradually make
the electrons localized (by some means), will the positive MR change
to a negative MR as the electrons get localized?} We investigate this
question in the metallic perovskite LaNiO$_{3-\delta}$ in which we
find that the MR changes from positive to negative as the material is
made more resistive by increasing $\delta$. LaNiO$_3$ belonging to
the ABO$_3$ class is structurally similar to LaCoO$_3$ and LaMnO$_3$.
However unlike the cobaltate and the manganate, it is a fairly good
metal with low resistivity \cite{NYV,DESHI,RLG,REGI,SANCH}. However,
as the oxygen deficiency is increased the resistivity is enhanced.
Eventually the material becomes an insulator at $\delta$ > 0.25
\cite{SANCH}. The investigation of the MR in this oxide is thus of
importance because it allows us to investigate the MR in a class of
disordered oxides where a change in oxygen stoichiometry leads to a
qualitative change in the nature of electronic transport without
bringing in strong magnetic effects atleast for $\delta \leq$ 0.2.

LaNiO$_3$ is an interesting material because of its potential in
application as metallic interconnects (or electrodes) in thin film
oxide electronics particularly those needing epitaxial multilayer
perovskite oxide films \cite{HEGDE}. The MR data of this oxide will
thus be of importance technologically.

\section{Experimental details}
The samples of LaNiO$_{3-\delta}$ were prepared from co-precipitation
and subsequent decomposition of nitrates of La and Ni as described in
refs.~[8] and [12]. The samples were heat treated
under different oxygen atmosphere to get varying oxygen content and
the oxygen stoichiometry $\delta$ was determined by iodometric
titration. The grain sizes as measured by scanning electron
microscope were typically of the order of 2 $\mu$m and there was no
variation in the grain sizes of the samples with different $\delta$.
[See ref [13] for details of sample preparation and charactrisation].

The electrical resistivity $\rho$(T) and the MR were measured by a
high precision low frequency (20 Hz) ac bridge technique with a
precision of $\pm$10ppm.  The resistivity was measured in the range
0.4K -300K and the MR was measured in the range 1.5K- 30K in a field of
up to 6 Tesla.

\section{Results and Discussion}
In figure~1 we present the resistivity measurements for the three
samples with the general formula LaNiO$_{3-\delta}$.  We identify the
samples in decreasing order of $\delta$ and have called them A,B and
C for $\delta$ = 0.14, 0.08 and 0.02 respectively. For all the
samples the $\rho$(T) curve show a metallic behaviour although the
residual resistivity $\rho_{0}$ increases by a factor of 30 as
$\delta$ changes from 0.02 to 0.14. We restrict our investigations to
$\delta\leq$ 0.14 to avoid any complications arising from magnetic
clusters which appear for $\delta \geq$ 0.2 \cite{SANCH}.

For the sample with $\delta$=0.14 which has the highest $\rho_0$ the
mean free path ($\ell_{el}$), estimated from the simple free electron
approximation, is somewhat smaller than the Ni-O bond length
($\approx$ 2\AA). Interestingly, it can be seen that for this sample
there is a resistivity minima around 20K. This is very similar to the
behavior of highly disordered metallic alloys.

The exact temperature dependence of $\rho$(T) in LaNiO$_{3-\delta}$
with varying $\delta$ has been critically looked into \cite{GAYA}. We
do not discuss these in detail here. We just mention that in the samples
B and C which show no resistivity minima we have a power law dependence
of the resistivity with T given as \cite{GAYA}:

\begin{equation}
\rho(T) = \rho_0 + \alpha T^{1.5} + \beta f(T)  \label{eq:fullt}
\end{equation}

\noindent where $f$(T) is the function giving a linear temperature
dependence at high temperature (T~>~200K) and a T$^n$ (n$\geq$ 2)
dependence for low temperatures. $\rho_0$ is the residual resistivity
and $\alpha$ and $\beta$ are constants. For the sample A, which shows
the highest resistivity, we have a resistivity minima
at T$_{min}$~=~20K and $\rho$(T) for T~<~T$_{min}$ shows a shallow
rise (see figure~1). We find that this shallow rise of $\rho$ below
T$_{min}$ can be best accounted for by using expressions for weak
localization (WL). We also find below that the MR measurements
supports this view. We quantify the weak localization
\cite{TVR} contribution by an additional term -$\gamma$T$^{0.75}$
which is incorporated into eqn.\ref{eq:fullt}. Thus, we fit the
resistivity of sample A using the expression

\begin{equation}
\rho = \rho_0 - \gamma T^{0.75} + \alpha T^{1.5} + \beta f(T)
\label{eq:afullt}
\end{equation}

\noindent We find that this relation fits the resitivity curve with an error
better than $\pm$ 0.5\% over the entire temperature range. The fits
for the samples are shown as a solid line in figure~1.

In figure~2 we show the MR of the three samples as a function of the
applied field H at T = 4.2K. We define the MR as
$\Delta\rho$/$\rho_{0H}$ = ($\rho$(H) - $\rho$(0))/$\rho_{0H}$, where
$\rho_{0H}$ for the zero field resistivity at a given temperature T.
The magnetoresistance of the LaNiO$_{3-\delta}$ system is small but
has a distinct dependence on the oxygen stoichiometry. It is positive
for the least resistive sample C, it is negative for sample A which
has the highest resistivity and it is almost zero (< 0.05\%) for
sample B with intermediate resistivity. The MR thus becomes
progressively negative as the resistivity increases.

The field dependence of MR also changes as $\delta$ increases. For
sample C, MR has a quadratic dependence on H which is shown in figure~3.
The MR seems to follow K\"ohler rule \cite{KHOLER}. The MR for this
sample is positive at all T. According to Koehler rule the dependence
of MR on temperature, magnetic field and purity of the sample should
be a function of H$\tau$ only, where $\tau$ is the relaxation time.
Since $\rho_{0H} \propto
\tau^{-1}$, $\Delta\rho/\rho_{0H} \propto$ $\it{k}$(H/$\rho_{0H}$),
where $\it{k}$(H/$\rho_{0H}$) is a function of H/$\rho_{0H}$.  For
normal metals the MR shows a quadratic dependence with field due to
the band MR. Hence we can write:

\begin{eqnarray}
\frac{\Delta\rho}{\rho_{0H}} = \beta_\rho\left(\frac{H}{\rho_{0H}}\right)^2
\label{eq:kho}
\end{eqnarray}

\noindent where, $\beta_\rho$ is a constant for the material. Fitting
the MR of sample C with the above expression we get $\beta_\rho$ =
3.98$\times$10$^{-12}$ ($\Omega$cm/T)$^2$. For normal metals the value of
$\beta_\rho \approx$ 10$^{-15}$ ($\Omega$cm/T)$^2$. This difference
is very interesting because if the positive MR indeed arises from a
band mechanism then it is one of the largest MR arising from the band
mechanism.  It will be of interest to compare the MR of this oxide
with other oxides showing positive MR \cite{SRNB,NG}. The other known
metallic oxides showing positive MR, Sr$_{x}$NbO$_3$ (x $\leq$ 0.8)
\cite{SRNB} and Na$_{0.9}$WO$_3$ \cite{NG} shows a H$^2$ dependence
over a small field range (H $\leq$ 3 Tesla). The values of
$\beta_\rho$ for Sr$_{0.85}$NbO$_3$ ($\rho_{4.2K}$ = 1.75m$\Omega$cm)
and Na$_{0.9}$WO$_3$ ($\rho_{4.2K}$ = 3$\mu\Omega$cm) are
1.127$\times$10$^{-9}$ ($\Omega$cm/T)$^2$ and 4x10$^{-15}$
($\Omega$cm/T)$^2$ respectively.  It seems therefore in these oxides
the constant $\beta_\rho$ increases as $\rho$ increases. This
observation needs more investigation to establish whether this indeed
is a general feature of metallic oxides.

For sample A the MR is not only negative, the field dependence is
also qualitatively different. We will show below that it is possible
to explain the field dependence of the negative MR as arising from
destruction of weak localisation on application of magnetic field.
For sample B, the MR being almost zero we didn't attempt any fit.
Most likely for this sample we have almost equal contribution from
both the positive and negative terms.

The temperature dependence of the MR also carries a distinct mark of
the oxygen stoichiometry. As shown in figure~4, the MR of sample A is
negative for T < T$_{min}$ and positive for T > T$_{min}$, where
T$_{min}$ is the temperature of resistivity minima (see figure~1). It
has been discussed before that the resistivity minima can arise due
to the contribution by WL for T < T$_{min}$. For T > T$_{min}$ the MR
is positive as in sample C. For T < T$_{min}$ the WL is supressed by
the application of the magnetic field and hence there will be a
negative contribution to the MR. To look at this quantitatively we
have fitted the magnetoconductance curves to the WL expression
\cite{KAWA,BAX} with an additional H$^2$ term to account for the
K\"ohler type negative magnetoconductance (or positive MR). In the
following we used conductivity change ($\delta\sigma$) instead of
resistivity change ($\delta\rho$) to follow the customary procedure
for WL expressions.  The expression thus used was \cite{BAX}:

\begin{eqnarray}
\delta\sigma =
\frac{e^2}{2\pi^2\hbar}{\left(\frac{eH}{\hbar}\right)}^{1/2}
{\cal H}\left(\frac{4{\cal D}\tau_\phi eH}{\hbar}\right) -
\beta_\sigma H^2
\label{eq:mr}
\end{eqnarray}

\noindent where
\begin{eqnarray}
{\cal H}(x) = 2\left[\sqrt{2+\frac{1}{x}}-\sqrt{\frac{1}{x}}\right]
-\left(\frac{1}{2}+\frac{1}{x}\right)^{-1/2} -
\left(\frac{3}{2}+\frac{1}{x}\right)^{-1/2} +
\frac{1}{48}\left(2.03+\frac{1}{x}\right)^{-3/2} \nonumber
\end{eqnarray}

\noindent Here, ${\cal D}$ is the electron diffusivity, $\tau_\phi$
is the electron phase relaxation time and $\beta_\sigma$ is a
constant. In this expression $\beta_\sigma$ and $\cal D\tau_\phi$ can
be free parameters. However, to restrict the number of free
parameters we fixed the parameter $\beta_\sigma$ following K\"ohler's
rule using the following strategy.

The K\"ohler's rule (see eqn \ref{eq:kho}) expressed in terms of the
conductivity is given by $\delta\sigma_{Kohler}$ = -$\beta_\sigma$
H$^2$ where $\beta_\sigma$ = $\beta_\rho/\rho_{0H}^3$. Assuming the
validity of the K\"ohler's rule for the samples A and C would imply
that both samples will have the same $\beta_\rho$. Thus, using
$\beta_\rho$ as found from sample C in figure~3 we have estimated the
value of $\beta_\sigma$ for sample A. The values of $\beta_\sigma$
for T $\leq$ 4.2K are $\approx$ 0.01 (S/cm)/T$^2$. Using these
$\beta_\sigma$ we fit the magnetoconductance data to eqn
\ref{eq:mr} using $\cal D\tau_\phi$ as the only fit parameter. A good
fit can be obtained for the magnetoconductance data as shown in the
inset of figure~4. The $\cal D\tau_\phi$ values (fit parameters) for the
two temperatures were obtained as 1.26x10$^{-12}$ cm$^2$ for T=2.3K
and 1.12x10$^{-12}$ cm$^2$ for T=4.2K.  Using the above values of
$\cal D\tau_\phi$, we have estimated the value of the WL correction
to the zero field conductivity of the sample. We get
$\delta\sigma_{WL}$ ($\approx (e^2/2\pi^2\hbar) (1/\sqrt{\cal
D\tau_\phi})$ to be about 12.3 S/cm which is about a factor of 2.5
more than that estimated from the resistivity data of this sample at
low temperature ($\approx$ 4.8 S/cm). Within the limits of
experimental error we can say that this is a good agreement since,
the estimate of $\delta\sigma_{WL}$ from the resistivity involves
some uncertainities because of the presence of other terms with much
larger contributions (see eqn.~2). Also we have only an approximate
value for $\beta_{\sigma}$. We can thus conclude that the negative MR
seen in this oxide can arise from the weak localization contribution.
At T~=~25~K the positive contribution wins over. But given the
smallness of MR it is difficult to obtain any meaningful fit.

One important effect that we should address is the formation of
local magnetic moments and any ferromagnetic interaction arising from
it which can give rise to negative MR.
While stoichiometric LaNiO$_{3}$ is a Pauli paramagnet,
creation of Ni$^{2+}$ ions in oxygen deficient samples can lead to a
Curie-Weiss type contribution. We have measured the susceptibility for
sample A. The susceptibility and the analysis are reported in ref.[13].
We mention the important results, since they are relevant to this paper.
For T<100K we could fit the susceptibility to $\chi$ = $\chi_0$+
C/(T-$\theta$). A fit to the
expression gives $\chi_0$ = 5.95x10$^{-6}$ emu/gm, $\theta$ = -51~K
and C = 3.83x10$^{-4}$ emu~ K/gm. The large value of C can arise from
the Ni$^{2+}$ formation as $\delta$ increases \cite{GAYA}.  For
sample A, for which $\delta$ = 0.14, about 30\% of the Ni in the
system is in the Ni$^{2+}$ state. However, not all of the Ni$^{2+}$
carry moments because part of the Ni$^{2+}$ which are in square
planar cordination are diamagnetic. Though there is a net small
antiferromagnetic interaction ($\theta \approx$ -51~K), the
susceptibility shows no sign of any magnetic ordering or cluster
formation. It can be contrasted to the composition LaNiO$_{2.75}$
which shows the transition to an insulating state and a clear
signature of a spin-glass or cluster glass like transition near 150K
\cite{SANCH}. From these observations we conclude that the magnetic
interactions will not play an important role for the MR of the
samples ($\delta \leq$ 0.14) which we have investigated. It may be
that for $\delta \geq$ 0.2 the cluster formation can lead to MR
dominated by magnetic interaction effects.

In summary, we have measured the magnetoresistance of the metallic
perovskite oxide LaNiO$_{3-\delta}$ as a function of the oxygen
defect $\delta$. The most metallic sample has a clear positive MR
which progressively goes to the negative magnitude as the oxygen
defect increases. The positive MR seems to arise from a band-like
contribution as seen in normal metals but the magnitude is much
larger than that seen in conventional metals. The negative MR seen in
the oxygen deficient sample arises from weak-localisation
contribution and is not due to any magnetic interactions. Our
experiment shows that the effects like weak-localisation of electrons
can lead to small but finite negative MR in these oxides which may be
masked by stronger negative MR arising from magnetic origin in
other oxides.

\vspace{2 truecm}
{\bf Acknowledgements:} One of us (AKR) wants to thank the Department
of Science and Technology, Government of India for a sponsored
project.  NG wants to thank the Council for Scientific and Industrial
Research for a fellowship. RLG acknowledges the supported by NSF
under Grant No: DMR9510475.

\newpage

Strongly Correlated Electrons (cond-mat.str-el)