Searches for the Most Metal-Poor Candidates from SDSS and SEGUE

Abstract: We report on efforts to identify large samples of very and extremely metal-poor stars based on medium-resolution spectroscopy and ugriz photometry obtained during the course of the Sloan Digital Sky Survey (SDSS), and its extension, SDSS-II, which includes the program SEGUE: Sloan Extension for Galactic Understanding and Exploration. To date, over 8000 stars with [Fe/H] <= -2.0 and effective temperatures in the range 4500K < T_eff < 7000K have been found, with the expected numbers in this temperature range to be well over 10,000 once SEGUE is completed. The numbers roughly double when one includes warmer blue stragglers and Blue Horizontal-Branch (BHB) stars in these counts. We show the observed low-metallicity tails of the Metallicity Distribution Functions for the cooler SDSS/SEGUE stars obtained thus far. We also comment on the confirmation of an inner/outer halo dichotomy in the Milky Way, and on how this realization may be used to direct searches for even more metal-poor stars in the near future.

\section{Introduction}

The primary recent searches for metal-poor stars in the Galaxy (the HK survey of
Beers and colleagues [1], and the Hamburg/ESO Survey of Christlieb and
colleagues [2]) were based on objective-prism techniques, followed by many years
of follow-up medium-resolution spectroscopy with intermediate-aperture
telescopes. Although this approach has been spectacularly successful, and has
resulted in the identification of several thousand stars with [Fe/H] $< -2.0$
(VMP stars according to [3]), and several hundred stars with [Fe/H] $< -3.0$
(EMP stars according to [3]), only three stars have been found to date with
[Fe/H] $< -4.0$ (UMP according to [3]), several of which are discussed in detail
by other contributions in this volume.

Clearly, much has already been learned from high-resolution spectroscopic
studies of many of the HK and HES stars, but many questions concerning the
nature of the first generations of objects born in the Galaxy
remain -- questions that can only be answered by obtaining a significant
increase in the numbers of VMP, EMP, and especially UMP and HMP ([Fe/H] $< -5.0$
according to [3]) stars. New methods and new strategies are required. Here we
report on one such effort, based on observations obtained during the SDSS and
SEGUE [4]. Other survey methods and approaches are discussed by Christlieb et
al. in this volume.

Kinematic studies of the SDSS stellar calibration objects, by Carollo et al. [5],
has confirmed the (long suspected) dichotomy of the halo of the Milky Way. Aside
from the obvious importance of this result for understanding the assembly and
evolution history of the halo itself, the fact that the outer-halo population
appears to exhibit a Metallicity Distribution Function (MDF) that peaks around
[Fe/H] $= -2.2$, roughly 0.5 dex lower than that of the inner-halo population,
opens up the possibility of specifically targeting the lowest metallicity stars
in the Galaxy {\it if one could efficiently pre-select} likely members
of the outer halo, as is discussed in more detail below.

\section{Metal-Poor Stars Discovered by the SDSS and SEGUE Surveys}

Metal-poor stars are identified from several sources within the SDSS and SEGUE.
In the original SDSS, numerous metal-poor stars are found among the 16 (color-
and apparent-magnitude selected) spectrophotometric and telluric calibration
stars that are targeted each time an exposure of a spectroscopic plug-plate
(containing a total of 640 fibers) is obtained. In addition, QSO candidates that
turn out not to be quasars, and BHB candidates that turn out to be somewhat
cooler than expected, based on spectroscopic follow-up, are often
low-metallicity stars. Additional metal-poor stars are found among the
calibration objects being observed during the ongoing LEGACY spectroscopy
program. The continuing SEGUE program specifically targets metal-poor stars, among
the 13 or so categories it considers. There are also many metal-poor stars that
are found among the F turnoff, G dwarf, and K giant target categories. As in
SDSS, each SEGUE plug-plate obtains 16 calibration stars as well. The total
number of stars observed spectroscopically by either SDSS or SEGUE is currently
over 350,000. It should be noted that the color selection, either for the
calibration objects or the low-metallicity candidates, is not capable of
discriminating between stars of metallicity lower than [Fe/H] = $-2.0$, since
the effect of declining metallicity on broadband stellar colors is minimal in
this regime. Hence, the MDF of such stars should be relatively unbiased below
this metallicity.

Stellar atmospheric parameters ($T_{\rm eff}$, log g, and [Fe/H]) are determined
for SDSS/SEGUE stars by application of the SEGUE Stellar Parameter Pipeline
(SSPP), which is described in detail by Lee at al. [6,7]. Tests of the SSPP
indicate that, over the temperature range 4500 $\le T_{\rm eff} \le 7500$~K,
external accuracies in the derived parameters are on the order 120~K, 0.25 dex,
and 0.2 dex, for $T_{\rm eff}$, log g, and [Fe/H], respectively.

Figure 1 shows the MDF for the over 4200 stars from SDSS with adopted [Fe/H]
$\le -2.0$ (and with $S/N \ge 10/1$). Figure 2 shows the same distribution for the
over 2400 similar stars observed during the first year of SEGUE. The addition of
observations from LEGACY and from the second year of SEGUE increases the total
numbers of stars with [Fe/H] $\le -2.0$ to well over 8000 objects in the
temperature range considered herein. This already
represents roughly a tripling of the total numbers of VMP stars obtained
by the sum of all previous surveys conducted in the past half century since the
low-metallicity star phenomenon was recognized.

\begin{figure}[!t]
\includegraphics[height=2.75in]{fig1.eps}
\caption{The MDF for VMP stars with 4500~K $ < T_{\rm eff} < 7000$~K
observed during the course of SDSS-I.  The great majority of these stars
are objects taken for spectrophotometric and telluric calibrations.}
\end{figure}

\begin{figure}[!t]
\includegraphics[height=2.75in]{fig2.eps}
\caption{The MDF for VMP stars with 4500~K $ < T_{\rm eff} < 7000$~K observed
during the course of the first year of SEGUE.  These
stars include objects specifically targeted to be likely VMP stars, as well as
calibration stars.}
\end{figure}

It is clear from the shapes of the MDFs in Figures 1 and 2 that
the total number of EMP stars (with [Fe/H] $\le -3.0$) is declining very
rapidly, with only a handful of stars being identified to date with [Fe/H] $\sim
-4.0$, and none below this value.  We address the reasons for this below.

\begin{figure}[!t]
\includegraphics[scale=0.7]{fig3.eps}
\caption{The components of space motion for the local sample of stars in
Carollo et al., shown as a function of [Fe/H].  The dashed lines indicate
the location of the Local Standard of Rest (LSR).  As can be seen in the middle panel,
the thick-disk and metal-weak thick disk population exhibits a small lag with
respect to the LSR.  Note also, in this panel, the large numbers of VMP stars
with large retrograde motions ($ V < -220$ km/s) with respect to the center of
the Galaxy.}
\end{figure}

\section{The Dichotomy of the Galactic Halo}

The region of the Milky Way halo, beyond the thin- and thick-disk systems, has long
been thought of as a single entity, comprising old stars and globular clusters
that represent the earliest populations of objects to have formed in our Galaxy.
Although many previous studies have speculated that the halo may in fact be
partitioned into inner and outer components, possibly with different kinematics
and spatial distributions, these studies have been limited by the generally
small numbers of stars or globular clusters upon which this split could be made.
Hence, it has always remained a suspicion, not a certainty.  Until now.
Based on medium-resolution spectroscopy for an initial sample of 20,366 SDSS
calibration stars located within roughly 20 kpc of the Sun, Carollo et al. [5] have
shown that the halo is clearly divisible into two broadly overlapping structural
components. These are: (1) the inner halo, which is dominated by stars on highly
eccentric orbits and exhibits a peak metallicity [Fe/H] $= -1.6$, as well as a
flattened density distribution with a modest net prograde rotation, and (2) the
outer halo, which includes stars that possess a wide range of orbital
eccentricities (including many on low eccentricity orbits), exhibits a peak
metallicity [Fe/H] $= -2.2$, and a spherical density distribution with a highly
statistically significant net retrograde rotation (between $V_{\phi}$ = $-50$ and
$-80$ km/s, depending on which subset of the data is considered).

The kinematic analysis of the Carollo et al. data set is carried out for the
subset of over 11,000 stars located within 4 kpc of the Sun, where one can take
advantage of the existence of reasonably accurate proper motions obtained from a
re-calibration of the USNO-B catalog, as described by Munn et al. [8]. Typical
errors in the derived proper motions are on the order of 2.5 to 3.5 mas/year,
which is more than adequate for this exercise.  Figure 3 shows the distribution
of derived $UVW$ components of the space motions for these stars as a function of
[Fe/H], as determined by the SSPP.  The impact of the size of this sample can be
appreciated by comparing with similar diagrams in the article by Frebel et al.
in this volume.

As is clear from inspection of Figure 3, one can clearly make out the presence of
the concentrations of stars with relatively low space motions and moderately high
metallicities, which we associate with membership in the thick-disk and metal-weak
thick-disk populations.  The very large velocity dispersions at lower
metallicity are of course associated with the halo population(s) (which do not
clearly separate from one another in such diagrams).

Stars in the disk population(s) can be effectively eliminated by considering
only those stars on retrograde orbits. Furthermore, as one takes subsets of the
stars that achieve maximum distances above the plane in the course of their
orbits about the center of the Galaxy (Z$_{\rm max}$) of 5 kpc or more, one can
address whether or not they exhibit properties consistent with that expected
from a single, or a more complex, halo population. Figure 4 shows the distribution
of [Fe/H] for stars with varying levels of increasingly retrograde orbits, and
with different cuts on Z$_{\rm max}$. As is clear from inspection of this
Figure, as one sweeps to more retrograde orbits, and to greater Z$_{\rm max}$,
the nature of the MDF changes over to favor stars at lower [Fe/H], and exhibits
fewer stars with higher [Fe/H]. This is due to the change in the relative
dominance of the outer-halo population over the inner-halo population at low
metallicity and at distances greater than roughly 15-20 kpc from the Galactic
center. This dichotomy may have a profound influence on searches for the most
metal-poor stars in the Galaxy in the future, as described below.

\begin{figure}
\includegraphics[scale=1.0] {fig4.eps}
\caption{Plots of the MDFs for stars on highly retrograde orbits with respect to the
Galactic center.  The left hand column corresponds to stars with $V_{\phi} < -100$
km/s, while the middle and right-hand columns correspond to stars with more extreme
retrograde motions, $V_{\phi} < -150$ km/s and $V_{\phi} < -200$ km/s,
respectively.  From top to bottom in each column, the panels show the MDFs for
stars with increasing cuts on Z$_{\rm max}$.  The vertical lines at [Fe/H]$ =
-2.0$ are shown for reference. Note that the nature of the MDFs changes within
each column, gradually favoring stars with [Fe/H] $< -2.0$ at higher Z$_{\rm
max}$. This change is even more evident (even though the total numbers of stars
are fewer) for the most extremely retrograde columns. All of the cuts on
individual panels with Z$_{\rm max} > 5$ kpc are highly inconsistent, according to a K-S
test, with being drawn from the same parent MDF as that of the low Z$_{\rm
max}$ cuts shown in the upper panels.}
\end{figure}

\section{Future Searches for the Most Metal-Poor Stars}

There are two primary techniques that have been exploited to date for
the identification of low-metallicity stars in the Galaxy.  The first
used selection of likely halo members based on observations of stars in the
local volume of the Galaxy with relatively high proper motions, e.g., Ryan and
Norris [9], or Carney et al. [10]. The second was based on the identification of
stars from $in-situ$ surveys (such as the HK survey and the HES) that sought to
find halo members based on low-resolution objective-prism spectroscopy to a
limiting magnitude between $B \sim 15.5$ (HK survey) and $B \sim 17.5$ (HES).
The SDSS/SEGUE efforts are color-selected $in-situ$ surveys, but many of
the VMP stars identified to date are still relatively local.

The interpretations of (in particular) the kinematics of stars selected by these
two techniques has been the subject of much debate in the literature over the
years. Now we understand why. Stars selected on the basis of high proper motions
favor the identification of members of the outer-halo population, while the
apparent-magnitude-limited prism (and color-selected) techniques favor stars
from the inner-halo population. We presently estimate that no more than 20\% of
halo stars within 4 kpc from the sun are members of the outer-halo population,
but this number has to be confirmed from more detailed modeling. The dominance
of the inner halo continues to at least 10-15 kpc from the Sun; beyond that the
outer-halo population begins to take over. Although there certainly exist
EMP stars that are members of this inner-halo population, they are
completely swamped by the huge numbers of stars near the peak of the
inner-halo MDF around [Fe/H] $= -1.6$. It is surely no accident that the three
UMP (which includes two HMP) stars known at present all are either located at
distances greater than 10 kpc from the Sun, or, in the case of HE~1327-2326 (see
Frebel et al. [8], and also Frebel et al., this volume), exhibit kinematics that
place them firmly in the outer halo.

The clear recognition that outer-halo stars are drawn from a population with
a lower peak in its MDF ([Fe/H] $= -2.2$) than that of the inner-halo stars
([Fe/H]$ = -1.6$), and the smaller number of stars with [Fe/H] $> -2.0$ that
appear among the outer-halo object,s suggests that one might seek to exploit
these facts in future searches for the most metal-deficient stars. One could, for
example, look to survey fainter stars, in particular distant giants, in order to
reach stars well outside the region where the inner-halo dominates.
Alternatively, one could make use of proper-motion surveys (now and in the
future) to identify stars in the local volume with large components of their
motions in the retrograde direction (or very high $U$ or $W$ motions), since
this appears to strongly favor outer-halo membership. In this manner one might
hope to efficiently identify sufficiently bright UMP and HMP stars that can be
studied at high resolution with current generation 8m-10m class telescopes, and
perhaps one day identify the long sought Mega Metal-Poor (MMP, with [Fe/H] $<
-6.0$, according to [3]) stars.

\begin{theacknowledgments}

Funding for the SDSS and SDSS-II has been provided by the Alfred P. Sloan
Foundation, the Participating Institutions, the National Science Foundation, the
U.S. Department of Energy, the National Aeronautics and Space Administration,
the Japanese Monbukagakusho, the Max Planck Society, and the Higher Education
Funding Council for England. The SDSS Web Site is http://www.sdss.org/.

The SDSS is managed by the Astrophysical Research Consortium for the
Participating Institutions. The Participating Institutions are the American
Museum of Natural History, Astrophysical Institute Potsdam, University of Basel,
University of Cambridge, Case Western Reserve University, University of Chicago,
Drexel University, Fermilab, the Institute for Advanced Study, the Japan
Participation Group, Johns Hopkins University, the Joint Institute for Nuclear
Astrophysics, the Kavli Institute for Particle Astrophysics and Cosmology, the
Korean Scientist Group, the Chinese Academy of Sciences (LAMOST), Los Alamos
National Laboratory, the Max-Planck-Institute for Astronomy (MPIA), the
Max-Planck-Institute for Astrophysics (MPA), New Mexico State University, Ohio
State University, University of Pittsburgh, University of Portsmouth, Princeton
University, the United States Naval Observatory, and the University of
Washington.

T.~C.~B acknowledges support from grants AST 04-06784,  AST 07-07776, and PHY
02-16783, Physics Frontier Centers/JINA: Joint Institute for Nuclear
Astrophysics, awarded by the U.S. National Science Foundation.  D.~C. is
grateful to JINA for support of her long-term visitor status at Michigan State
University, where this analysis took place.

\end{theacknowledgments}

\bibliographystyle{aipprocl}


Emergence of Antiferromagnetic Correlation in LiTi2-xVxO4 via 7Li NMR

Abstract: We report 7Li NMR studies of V-substitution effects on spinel oxide superconductor LiTi2O4 (Tc = 13.4 K). In LiTi2-xVxO4 (x = 0-0.4), the V substitution for the Ti site suppressed the relative volume fraction of superconductivity faster than Tc. From the observation of fairly homogeneous enhancement in a 7Li nuclear spin-lattice relaxation rate, we conclude that the V substitution changes electron correlation effects through electron carrier doping from quarter electron filling 3d0.5 to 3d1.5 and then the antiferromagnetic correlation emerges.

\section{Introduction}
Transition-metal spinels have attracted great interests,
because intriguing electronic phases and their competition are associated with
frustration effects on spin, charge, orbital networks and the coupling with crystal lattice~\cite{JP2007}.
A normal spinel LiT$_2$O$_4$ is a type-II  oxide superconductor with
a high $T_\mathrm{c}$ $\sim$13 K~\cite{Johnston,Mos}.
The crystal structure is stable at low temperatures.
LiTi$_2$O$_4$ is an $s$-wave superconductor,
because of a Hebel-Slichter peak in $^{7}$Li nuclear spin-lattice relaxation rate~\cite{MItoh} and activation-type specific heat~\cite{Sun}.
The mystery is the unconventional break down of the superconductivity.
Although the coherence length $\xi$ is long enough~\cite{Sun},
the superconductivity is suppressed as the decrease in the superconducting volume fraction
for Li-deficient Li$_{1-\delta}$Ti$_2$O$_4$~\cite{Rygula} and Li-rich Li$_{1+y}$Ti$_{2-y}$O$_4$~\cite{Johnston,Ueda}.
The Li deficiency and Li substitution for Ti site make the conduction band from the quarter filling to
the band insulating states.

The full solid solution of LiTi$_{2-x}$V$_x$O$_4$ (0$\leq x\leq$ 2) is known~\cite{Tsuda}.
An itinerant-electron spinel LiV$_2$O$_4$ shows a heavy electronic specific heat and a Curie-Weiss spin susceptibility~\cite{Kondo}.
Although the ferromagnetic spin fluctuations were estimated from NMR data~\cite{Fujiwara},
the antiferromagnetic fluctuations were observed at low temperatures by neutron scattering experiments~\cite{ND1,ND2}.
The spin fluctuation spectrum may not be simple~\cite{ND2,PES}.
Ti substitution for LiV$_2$O$_4$ induces spin-glass-like behaviors and the microscopic effects were studied~\cite{Trinkl}.
The microscopic effects of V substitution for LiTi$_2$O$_4$, however, have been poorly understood.

One may expect two electronic effects of V substitution for Ti sites in LiTi$_2$O$_4$.
One is a simple pair-breaking effect of additional impurity electron spin $\Delta S$ = 1 on the superconductivity.
The V electrons are assumed to have localized moments.
The randomly distributed local moments can induce non-exponential NMR relaxation~\cite{Mc,ItohNi,ItohZn,JohnstonNMR}.
The other is the carrier doping effect on electron correlation through band filling from quarter electron filling 3$d^{0.5}$ to 3$d^{1.5}$.
The V electrons are assumed to hybridize the Ti conduction electrons and then to be itinerant.
Even for the random crystalline potentials introduced through the substitution,
NMR linewidths may be broader but
the NMR relaxation can have a single spin-lattice relaxation time~\cite{ItohSr}.

In this Letter, we report the $^{7}$Li NMR studies of LiTi$_{2-x}$V$_x$O$_4$ ($x$ = 0 $-$ 0.4).
The V substitution was found to suppress the relative volume fraction of superconductivity
faster than $T_\mathrm{c}$ and
to enhance fairly homogeneously and largely the $^{7}$Li nuclear spin-lattice relaxation but not so much the Knight shift.
The antiferromagnetic correlation was concluded to emerge in the V substituted samples.

\begin{figure}[bp]
	\begin{center}
		\includegraphics[width=7cm, clip]{f1.eps}
	\end{center}
	\caption{(a) Uniform magnetic susceptibility $M/H$ at 100 G after cooling in a zero field (ZFC) and in a field of 100 G (FC) for LiTi$_{2-x}$V$_x$O$_4$.
	(b) V concentration dependence of $T_\mathrm{c}$ and the relative volume fraction of superconductivity to $x$ = 0 for LiTi$_{2-x}$V$_x$O$_4$. The solid curves are visual guides.
	The bifurcation of temperature hysteresis and the onset of diamagnetic response of $M/H$ were defined as the superconducting transition temperature $T_\mathrm{c}$.
	}
		\label{fig_Tc}
\end{figure}%
Powder samples of LiTi$_{2-x}$V$_x$O$_4$ were synthesized by a solid state reaction method with a precursor of Li$_{4/3}$Ti$_{5/3}$O$_4$, because Ti is easily oxidized to be Ti$^{4+}$~\cite{Ueda}.
First, Li$_{4/3}$Ti$_{5/3}$O$_4$ was synthesized from the mixture of preheated dry Li$_2$CO$_3$ (99.9 $\%$) and TiO$_2$ (99.99 $\%$) after ref. 7.
Next, the mixtures of Li$_{4/3}$Ti$_{5/3}$O$_4$, TiO$_2$, Ti and V$_2$O$_3$ were sealed in evacuated quartz tubes and then fired at 760 $^\circ$C for 0$\leq x \leq$0.04, at 800 $^\circ$C for 0.08$\leq x \leq$0.16, and  at 850 $^\circ$C for 0.2$\leq x \leq$0.4 in a week.
The powder X-ray diffraction patterns indicated all the samples of $x \leq$ 0.04 in a single phase with spinel structure and
the samples of $x \geq$ 0.08 with a small amount of unreacted V$_2$O$_3$.
The cryopreservation method at 77 K in liquid nitrogen was employed to keep the samples fresh.

We performed high resolution Fourier-transformed $^{7}$Li (nuclear spin $I$ = 3/2 and nuclear gyromagnetic ratio $\gamma_{\rm n}$/2$\pi$ = 16.546 MHz/T) NMR measurements
of free induction decay signals or the nuclear spin-echoes at $H$ = 7.48414 T.
The applied magnetic field was estimated from the reference material LiCl$aq$.
Nuclear spin-lattice relaxation times were measured by an inversion recovery technique.

\begin{figure}[tp]
	\begin{center}
		\includegraphics[width=6.8cm, clip]{f2.eps}
	\end{center}
	\caption{V substitution effects on temperature dependence of magnetic susceptibility $\chi$ defined by $\Delta M/\Delta H$ at 4 and 5 T for $x$ = 0, 0.04, and 0.16.
	The inset shows $\Delta M/\Delta H$ for pure LiTi$_2$O$_4$ on an expanded scale.
	}
		\label{fig_XT}
\end{figure
Magnetization was measured by a SQUID magnetometer (Quantum Design MPMS)
for $x$ = 0, 0.01, 0.02, 0.04, 0.08, 0.16, 0.2 and 0.4.
 Figure \ref{fig_Tc} (a) shows uniform magnetic susceptibility $M/H$ at 100 G after cooling in a zero field (ZFC) and in a field of 100 G (FC) for LiTi$_{2-x}$V$_x$O$_4$ with $x$ = 0, 0.08, 0.12, 0.16, 0.2. The temperature hysteresis and the onset of diamagnetic response in $M/H$ diminished for $x >$ 0.12.
Figure \ref{fig_Tc} (b) shows V concentration dependence of $T_\mathrm{c}$ and the relative volume fraction of superconductivity to $x$ = 0 at 5 K.
The V impurities suppress the relative volume fraction faster than $T_\mathrm{c}$~\cite{SunV}.

Low field magnetization curves were non-linear even at 300 K, but high field ones ($H >$ 0.5 T) were linear.
Thus, we defined the intrinsic magnetic susceptibility of  LiTi$_{2-x}$V$_x$O$_4$ in the normal states
by the difference in magnetization $M$ at $H$ = 4 and 5 T, $\Delta M/\Delta H$.
The non-linear magnetization at lower fields than 0.4 T might be due to unintentional magnetic impurities (minute impurity phase).
Figure \ref{fig_XT} shows paramagnetic susceptibility $\Delta M/\Delta H$ for $x$ = 0, 0.04 and 0.16.
The V impurities induce Curie-Weiss like behaviors~\cite{Tsuda}.

\begin{figure}[tbp]
	\begin{center}
		\includegraphics[width=6.8cm, clip]{f3.eps}
	\end{center}
	\caption{Fourier-transformed $^{7}$Li NMR spectra of LiTi$_{2}$O$_4$ (a) and  LiTi$_{1.92}$V$_{0.08}$O$_4$ (b). The peak frequency of a $^{7}$Li NMR spectrum of LiCl$aq$ is denoted by $\nu_\mathrm{ref}$ for reference.
	(c) Temperature dependences of $^{7}$Li Knight shifts for LiTi$_{2-x}$V$_{x}$O$_4$. Note the direction of the vertical axis to the negative values. The sold lines are visual guides.
	}
		\label{fig_FT}
\end{figure}%
Figure \ref{fig_FT} shows $^{7}$Li NMR frequency spectra as a function of temperature for $x$ = 0 (a) and $x$ = 0.08 (b).
All the $^{7}$Li NMR spectra except $x$ = 0 below $T_\mathrm{c}$ are symmetric.
No quadrupole splits are observed.
The peak frequencies show negative shifts and decrease as temperature is decreased.
The linewidths are broadened by the V substitution.
This evidences the actual substitution of the V ions for the Ti ions.

Figure \ref{fig_FT}(c) shows the temperature dependences of $^{7}$Li Knight shifts  in LiTi$_{2-x}$V$_{x}$O$_4$ ($x \leq$ 0.16). The negative shifts show Curies-Weiss behaviors. The  $^{7}$Li Knight shift of Curies-Weiss sort in LiV$_{2}$O$_{4}$ was positive~\cite{Fujiwara}.

\begin{figure}[tbp]
	\begin{center}
		\includegraphics[width=6.5cm, clip]{f4.eps}
	\end{center}
	\caption{$^{7}$Li Knight shifts $K$ are plotted against magnetic susceptibility $\chi$ defined by $\Delta M$/$\Delta H$ at $H$ = 4 and 5 T, where temperature is an implicit parameter, for LiTi$_{2-x}$V$_{x}$O$_4$ of $x$ = 0, 0.04 and 0.16. The solid lines are the results from least-squares fittings.
	The inset shows on an expanded scale the $K-\chi$ plot of pure LiTi$_2$O$_4$.
		}
		\label{fig_KX}
\end{figure
In Fig. \ref{fig_KX}, $^{7}$Li Knight shifts $K$ are plotted against magnetic susceptibility $\chi$ defined by $\Delta M$/$\Delta H$ at $H$ = 4 and 5 T, where temperature is an implicit parameter, for LiTi$_{2-x}$V$_{x}$O$_4$ with $x$ = 0, 0.04 and 0.16.
The inset shows on an expanded scale the $K-\chi$ plot of pure LiTi$_2$O$_4$.

The $^{7}$Li Knight shift consists of a spin shift $K_s$ and a chemical shift $\sigma$,
\begin{equation}
^{7}K = K_s + \sigma.
\label{eq_K1}
\end{equation}
The spin shift $K_s$ is expressed by a product of a hyperfine coupling constant $A_{\rm hf}$ and a temperature dependent spin susceptibility $\chi_s(T)$
\begin{equation}
K_s(T) = \frac{A_{\rm hf}}{N_{\rm A}\mu_{\rm B}}\chi_s(T),
\label{eq_K2}
\end{equation}
where $N_{\rm A}$ is the Avogadro number and $\mu_{\rm B}$ is the Bohr magneton.

The bulk magnetic susceptibility $\chi$ is given by
\begin{equation}
\chi = \chi_s(T) + \chi_{\rm orb} + \chi_{\rm dia},
\label{eq_K3}
\end{equation}
where $\chi_{\rm orb}$ is the Van Vleck orbital susceptibility and  $\chi_{\rm dia}$ is the diamagnetic susceptibility of inner core electrons.

In Fig. \ref{fig_KX}, the solid lines are the least-squares fitting results by $K(T)$ = $p\chi(T)$ + $q$ ($p$ and $q$ are the fit parameters).
The linear relation between $^{7}K$ and $\chi$ breaks down at lower temperatures for $x$ = 0.04 and 0.16. The bulk magnetic susceptibility must include Curie components being different from
the peak Knight shifts of the broadened $^{7}$Li NMR spectra.

From the fitting results, the hyperfine coupling constant $^{7}A_{\rm hf}$ was estimated to be $-$ 3.12, $-$ 1.05 and $-$ 0.74 kOe/$\mu_{\rm B}$Ti for $x$ = 0, 0.04 ($T >$ 100 K) and 0.16 ($T >$ 200 K), respectively.
In passing, the $^{7}A_{\rm hf}$ of LiV$_2$O$_4$ is positive~\cite{Fujiwara}.
From $\chi_{\rm dia}$ = $-$ 6.28$\times$10$^{-5}$ emu/f.u.mole and tentative $\chi_{\rm orb}$ = $+$ 3.33$\times$10$^{-5}$ emu/f.u.mole~\cite{TiO2},
the chemical shift $\sigma$ was estimated to be positive 0.002 $-$ 0.005 $\%$.
Although the magnitude of $\sigma$ depends on the choice of $\chi_{\rm orb}$,
the positive chemical shift $\sigma$ is unconventional.
\begin{figure}[htbp]
	\begin{center}
		\includegraphics[width=6.5cm, clip]{f5.eps}
	\end{center}
	\caption{V substitution effects on $^{7}$Li nuclear spin-lattice relaxation curves (recovery curves) in LiTi$_{2-x}$V$_x$O$_4$ ($x$ = 0$-$ 0.2).
	}
		\label{fig_rec}
\end{figure
\begin{figure}[htb]
	\begin{center}
		\includegraphics[width=6.5cm, clip]{f6.eps}
	\end{center}
	\caption{Temperature dependences of the exponent $\beta$ of $^{7}$Li nuclear spin-lattice relaxation curves (a) and 1/$T_{1}$ (b) in LiTi$_{2-x}$V$_x$O$_4$ ($x$ = 0$-$ 0.4), and LiV$_2$O$_4$ reproduced from ref. 10. The solid curves are visual guides.
		}
		\label{fig_T1}
\end{figure

Figure \ref{fig_rec} shows the $^{7}$Li nuclear spin-lattice relaxation curves (recovery curves) $p(t)\equiv 1-M(t)/M(\infty)$ ($t$ is the time after an inversion pulse to the observation pulse and $M(t)$ is the nuclear magnetization) for LiTi$_{2-x}$V$_x$O$_4$ ($x$ = 0 $-$ 0.4).
The solid curves are the results of the least-squares fitting using a stretched exponential function
\begin{equation}
p(t) = p(0){\rm exp}[-(t/T_1)^{\beta}],
\label{eq_rec}
\end{equation}
where $p(0)$, $\beta$ and $T_1$ are the fit parameters.
As seen in Fig. \ref{fig_rec}, the recovery curves for the V substitution are nearly single exponential functions except the low temperature for $x$ = 0.2.
The V substitution enhances 1/$T_1$ while keeping nearly the homogeneous spin-lattice relaxation.

Figure \ref{fig_T1}(a) shows the temperature dependence of $\beta$ for $x$ = 0 $-$ 0.4.
For $x <$ 0.2, the exponent $\beta >$ 0.8 is nearly independent of temperature,
while for $x$ = 0.2 and 0.4, the cooling down below 50 K leads to $\beta \rightarrow$ 0.5.
For the Li poor and rich samples, we observed $\beta \rightarrow$ 0.5 immediately when the Li deficiency and the Li substitution for Ti site~\cite{MItoh} are introduced into LiTi$_2$O$_4$.
Thus, nearly the single exponential functions exclude the deviation of Li composition.
The observed V impurity effect at the high magnetic field $H \sim$ 7.5 T for $x <$ 0.2 is in contrast to the conventional magnetic impurity effect on the NMR relaxation, where the non-exponential recovery curves are induced only at low fields and low temperatures ($\beta$ =1 at 300 K is reduced to 0.5 at low temperatures) and easily suppressed by the high magnetic field of $H\sim$ 7.5 T~\cite{Mc,ItohNi}.

Figure \ref{fig_T1}(b) shows the temperature dependence of 1/$T_1$ for $x$ = 0 $-$ 0.4.
For pure LiTi$_2$O$_4$, 1/$T_1$ shows $T^{0.7}$ dependence.
With the V substitution, 1/$T_1$ is highly enhanced
and the temperature dependence below 100 K is changed into $T^{n}$ with $n \leq$ 0.7.
Evidently, the enhancement of 1/$T_1$ due to the V substitution is larger than that of the magnitude of the Knight shift $K$. That is the emergence of the antiferromagnetic correlation.

Electron correlation changes the Korringa ratio~\cite{Moriya}.
The modified Korringa relation is characterized by
\begin{align}
K(\alpha)=\frac{\gamma_e}{\gamma_{\rm n}}\frac{\hbar}{4\pi k_{\rm B}}\frac{1}{T_1TK^2_s},
\label{eq:ka}
\end{align}
where $\gamma_e$ is the electron gyromagnetic ratio and
$\alpha$ is the exchange enhancement factor~\cite{Moriya}.
The value of $K(\alpha)$ is associated with the wave vector (${\bm q}$) dependence of a generalized spin susceptibility $\chi$($\bm q$, $\omega$), that is the ratio of the $\bm q$-averaged $\chi\prime$($\bm q$) to the uniform $\chi\prime$($\bm q$ = 0).
The ferromagnetic and antiferromagnetic $\chi\prime$($\bm q$) lead to $K(\alpha) <$ 1 and $K(\alpha) >$ 1, respectively. Incommensurately enhanced  $\chi\prime$($\bm q$) also leads to $K(\alpha) >$ 1.

One should note that the $\bm q$ dependence of the hyperfine coupling constant $A_{\rm hf}$ also plays a significant role.
A Li site (8a) has 12 nearest neighbor Ti sites (16d).
The 12 Ti ions located at the corners of a truncated tetrahedron surround the Li ion at the center of the tetrahedron in a cubic spinel.
Thus, the staggered magnetic fields from the Ti electrons on four sublattices may be cancelled out and masked at the Li site.
The $\bm q$ dependent $A_{\rm hf}({\bm q})$ can work as a filter to the staggered mode.

\begin{figure}[htp]
	\begin{center}
		\includegraphics[width=6.0cm, clip]{f7.eps}
	\end{center}
	\caption{V subsitution effect on $K(\alpha)$ of eq.(~\ref{eq:ka}).
	The broken line indicates $K(\alpha)$ = 1.  The solid curves are visual guides.
	}
		\label{fig_Ka}
\end{figure
Figure~\ref{fig_Ka} shows the V substitution effect on the modified Korriga ratio $K(\alpha)$ defined by eq.(~\ref{eq:ka}).
The V substitution changes $K(\alpha) <$ 1 for pure LiTi$_2$O$_4$ into $K(\alpha) >$ 1.
The emergence of the antiferromagnetic correlation due to the V substitution, not a conventional quantum phase transition, was seen through the Li NMR.
One should note that the enhanced 1/$T_1$ may be due to incommensurate magnetic correlation,
which is not filtered by $A_{\rm hf}({\bm q})$.

The V substitution effect is not a simple pair-breaking effect on the superconductivity
but also the change of the electron correlation effect of the 3$d$ $t_{2g}$ band
from 3$d^{0.5}$ filling to 3$d^{1.5}$.
The electron carrier doping via V substitution is consistent with the sign of Seebeck coefficient~\cite{Tsuda}.

In conclusion, we observed V-induced enhancement
in $^{7}$Li nuclear spin-lattice relaxation rates in LiTi$_{2-x}$V$_x$O$_4$ ($x$ = 0 $-$ 0.4),
which indicates the emergence of the antiferromagnetic correlation.
The V substitution for Ti ions changes the electron correlation effects by controlling the band filling from quarter electron filling 3$d^{0.5}$ to 3$d^{1.5}$.

We thank T. Waki, H. Chudo, H. Ohta, A. Tanizawa, C. Michioka for their experimental assistance and helpful discussions.
This work was supported in part by a Grant-in-Aid for Science Research on
Priority Area, ``Invention of Anomalous Quantum Materials," from the Ministry
of Education, Culture, Sports, Science and Technology of Japan (Grant No. 16076210)
and in part by a Grant-in-Aid for Scientific Research from the Japan Society
for the Promotion of Science (Grant No. 19350030).

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