
A precise description of neutrinonucleus reactions will play a key role in
addressing fundamental questions such as the leptonic CP violation and the
neutrino mass hierarchy through analyzing data from nextgeneration neutrino
oscillation experiments. The neutrino energy relevant to the neutrinonucleus
reactions spans a broad range and, accordingly, the dominant reaction mechanism
varies across the energy region from quasielastic scattering through nucleon
resonance excitations to deep inelastic scattering. This corresponds to
transitions of the effective degree of freedom for theoretical description from
nucleons through mesonbaryon to quarks. The main purpose of this review is to
report our recent efforts towards a unified description of the neutrinonucleus
reactions over the wide energy range; recent overall progress in the field is
also sketched. Starting with an overview of the current status of
neutrinonucleus scattering experiments, we formulate the cross section to be
commonly used for the reactions over all the energy regions. A description of
the neutrinonucleon reactions follows and, in particular, a dynamical
coupledchannels model for meson productions in and beyond the $\Delta$(1232)
region is discussed in detail. We then discuss the neutrinonucleus reactions,
putting emphasis on our theoretical approaches. We start the discussion with
electroweak processes in fewnucleon systems studied with the correlated
Gaussian method. Then we describe quasielastic scattering with nuclear
spectral functions, and meson productions with a $\Delta$hole model. Nuclear
modifications of the parton distribution functions determined through a global
analysis are also discussed. Finally, we discuss issues to be addressed for
future developments.

Fragmentation functions are determined for the pion and kaon by global
analyses of chargedhadron production data in electronpositron annihilation.
Accurate measurements were reported by the Belle and BaBar collaborations for
the fragmentation functions at the centerofmass energies of 10.52 GeV and
10.54 GeV, respectively, at the KEK and SLAC B factories, whereas other
available $e^+e^$ measurements were mostly done at higher energies, mainly at
the $Z$ mass of 91.2 GeV. There is a possibility that gluon fragmentation
functions, as well as quark fragmentation functions, are accurately determined
by scaling violation. We report our global analysis of the fragmentation
functions especially to show impacts of the Bfactory measurements on the
fragmentation function determination. Our results indicate that the
fragmentation functions are determined more accurately not only by the scaling
violation but also by highstatistical nature of the Belle and BaBar data.
However, there are some tensions between the Belle and BaBar data in comparison
with previous measurements. We also explain how the flavor dependence of quark
fragmentation functions and the gluon function are separated by using
measurements at different $Q^2$ values. In particular, the electric and weak
changes are different depending on the quark type, so that a lightquark flavor
separation also became possible in principle due to the precise data at both
$\sqrt{s}\simeq 10.5$ GeV and 91.2 GeV.

We perform the first iterative Monte Carlo (IMC) analysis of fragmentation
functions constrained by all available data from singleinclusive $e^+ e^$
annihilation into pions and kaons. The IMC method eliminates potential bias in
traditional analyses based on single fits introduced by fixing parameters not
well contrained by the data and provides a statistically rigorous determination
of uncertainties. Our analysis reveals specific features of fragmentation
functions using the new IMC methodology and those obtained from previous
analyses, especially for light quarks and for strange quark fragmentation to
kaons.

A clustering aspect is explained for the $^9$Be nucleus in chargedlepton
deep inelastic scattering. Nuclear modifications of the structure function
$F_2$ are studied by the ratio $R_{\rm EMC} = F_2^A /F_2^D$, where $A$ and $D$
are a nucleus and the deuteron, respectively. In a JLab experiment, an
unexpectedly large nuclear modification slope $dR_{\rm EMC}/dx$ was found for
$^9$Be, which could be related to its clustering structure. We investigated a
mean conventional part of a nuclear structure function $F_2^A$ by a convolution
description with nucleon momentum distributions calculated by antisymmetrized
(or fermionic) molecular dynamics (AMD) and also by a simple shell model. We
found that clustering effects are small in the conventional part, so that the
JLab result could be associated with an internal nucleon modification or a
shortrange nuclear correlation which is caused by high densities due to
cluster formation.

Next generation neutrino oscillation experiments will need a quantitative
understanding of neutrinonucleus interaction far better than ever. Kinematics
covered by the relevant neutrinonucleus interaction spans wide region, from
the quasielastic, through the resonance region, to the deeply inelastic
scattering region. The neutrinonucleus interaction in each region has quite
different characteristics. Obviously, it is essential to combine different
expertise to construct a unified model that covers all the kinematical region
of the neutrinonucleus interaction. Recently, several experimentalists and
theorists got together to form a collaboration to tackle this problem. In this
contribution, we report the collaboration's recent activity and a goal in near
future.

Semiinclusive hadronproduction processes are becoming important in
highenergy hadron reactions. They are used for investigating properties of
quarkhadron matters in heavyion collisions, for finding the origin of nucleon
spin in polarized leptonnucleon and nucleonnucleon reactions, and possibly
for finding exotic hadrons. In describing the hadronproduction cross sections
in highenergy reactions, fragmentation functions are essential quantities. A
fragmentation function indicates the probability of producing a hadron from a
parton in the leading order of the running coupling constant $\alpha_s$. Its
$Q^2$ dependence is described by the standard DGLAP
(DokshitzerGribovLipatovAltarelliParisi) evolution equations, which are
often used in theoretical and experimental analyses of the fragmentation
functions and in calculating semiinclusive cross sections. The DGLAP equations
are complicated integrodifferential equations, which cannot be solved in an
analytical method. In this work, a simple method is employed for solving the
evolution equations by using GaussLegendre quadrature for evaluating
integrals, and a useful code is provided for calculating the $Q^2$ evolution of
the fragmentation functions in the leading order (LO) and nexttoleading order
(NLO) of $\alpha_s$. The renormalization scheme is $\overline{MS}$ in the NLO
evolution. Our evolution code is explained for using it in one's studies on the
fragmentation functions.

We report recent studies on structure functions of the nucleon and nuclei.
First, clustering effects are investigated in the structure function F_2 of
Be9 for explaining an unusual nuclear correction found in a JLab experiment.
We propose that high densities created by formation of clustering structure
like 2*alpha+neutron in Be9 is the origin of the unexpected JLab result by
using the antisymmetrized molecular dynamics (AMD). There is an approved
proposal at JLab to investigate the structure functions of light nuclei
including the cluster structure, so that much details will become clear in a
few years. Second, tensorpolarized quark and antiquark distributions are
obtained by analyzing HERMES measurements on the structure function b_1 for the
deuteron. The result suggests a finite tensor polarization for antiquark
distributions, which is an interesting topic for further theoretical and
experimental investigations. An experimental proposal exists at JLab for
measuring b_1 of the deuteron as a new tensorstructure study in 2010's.
Furthermore, the antiquark tensor polarization could be measured by polarized
deuteron DrellYan processes at hadron facilities such as JPARC and GSIFAIR.
Third, the recent CDF dijet anomaly is investigated within the standard model
by considering possible modifications of the strangequark distribution. We
find that the shape of a dijetmass spectrum changes depending on the
strangequark distribution. It indicates that the CDF excess could be partially
explained as a PDF effect, particularly by the strangeness in the nucleon,
within the standard model if the excess at m_{jj}~140 GeV is not a sharp peak.

We report current status of global analyses on nuclear parton distribution
functions (NPDFs). The optimum NPDFs are determined by analyzing highenergy
nuclear reaction data. Due to limited experimental measurements, antiquark
modifications have large uncertainties at x>0.2 and gluon modifications cannot
be determined. A nuclear modification difference between u and d quark
distributions could be an origin of the longstanding NuTeV sin^2 theta_w
anomaly. There is also an issue of nuclear modification differences between the
structure functions of chargedlepton and neutrino reactions. Next, nuclear
clustering effects are discussed in structure functions F_2^A as a possible
explanation for an anomalous result in the Be9 nucleus at the Thomas Jefferson
National Accelerator Facility (JLab). Last, tensorpolarized quark and
antiquark distribution functions are extracted from HERMES data on the
polarized structure function b_1 of the deuteron, and they could be used for
testing theoretical models and for proposing future experiments, for example,
the one at JLab. Such measurements could open a new field of spin physics in
spinone hadrons.

An anomalous nuclear modification was reported by JLab measurements on the
beryllium9 structure function F_2. It is unexpected in the sense that a
nuclear modification slope is too large to be expected from its average nuclear
density. We investigated whether it is explained by a nuclear clustering
configuration in Be9 with two \alpha nuclei and surrounding neutron clouds.
Such clustering aspects are studied by using antisymmetrized molecular dynamics
(AMD) and also by a simple shell model for comparison. We consider that nuclear
structure functions F_2^A consist of a mean conventional part and a remaining
one depending on the maximum local density. The first mean part does not show a
significant cluster effect on F_2. However, we propose that the remaining one
could explain the anonymous JLab slope, and it is associated with high
densities created by the cluster formation in Be9. The JLab measurement is
possibly the first signature of clustering effects in highenergy nuclear
reactions. A responsible physics could be an internal nucleon modification,
which is caused by the high densities due to the cluster configuration.

For understanding an anomalous nuclear effect experimentally observed for the
beryllium9 nucleus at the Thomas Jefferson National Accelerator Facility
(JLab), clustering aspects are studied in structure functions of deep inelastic
leptonnucleus scattering by using momentum distributions calculated in
antisymmetrized (or fermionic) molecular dynamics (AMD) and also in a simple
shell model for comparison. According to the AMD, the Be9 nucleus consists of
two alphalike clusters with a surrounding neutron. The clustering produces
highmomentum components in nuclear wave functions, which affects nuclear
modifications of the structure functions. We investigated whether clustering
features could appear in the structure function F_2 of Be9 along with studies
for other light nuclei. We found that nuclear modifications of F_2 are similar
in both AMD and shell models within our simple convolution description although
there are slight differences in Be9. It indicates that the anomalous Be9
result should be explained by a different mechanism from the nuclear binding
and Fermi motion. If nuclearmodification slopes d(F_2^A/F_2^D)/dx are shown by
the maximum local densities, the Be9 anomaly can be explained by the AMD
picture, namely by the clustering structure, whereas it certainly cannot be
described in the simple shell model. This fact suggests that the large nuclear
modification in Be9 should be explained by large densities in the clusters.
For example, internal nucleon structure could be modified in the highdensity
clusters. The clustering aspect of nuclear structure functions is an unexplored
topic which is interesting for future investigations.

We use photoemission spectroscopy to study electronic structures of pristine
and Kdoped solid picene. The valence band spectrum of pristine picene consists
of three main features with no state at the Fermi level (EF), while that of
Kdoped picene has three structures similar to those of pristine picene with
new states near EF, consistent with the semiconductormetal transition. The
Kinduced change cannot be explained with a simple rigidband model of pristine
picene, but can be interpreted by molecular orbital calculations considering
electronintramolecularvibration interaction. Excellent agreement of the
Kdoped spectrum with the calculations points to importance of
electronintramolecularvibration interaction in Kdoped picene.

We discuss studies on determination of fragmentation functions and an
application to exotichadron search by using characteristic differences between
favored and disfavored functions. The optimum fragmentation functions are
determined for pion, kaon, and proton in the leading order (LO) and
nexttoleading order (NLO) of the running coupling constant alpha_s by global
analyses of hadronproduction data in electronpositron annihilation. Various
parametrization results are much different in disfavoredquark and gluon
fragmentation functions; however, we show that they are within uncertainties of
the determined functions by using the Hessian method for uncertainty
estimation. We find that the uncertainties are especially large in the
disfavoredquark and gluon fragmentation functions. NLO improvements are
explicitly shown in the determination by comparing uncertainties of the LO and
NLO functions. Next, we propose to use differences between favored and
disfavored fragmentation functions for determining internal quark
configurations of exotic hadrons. We make a global analysis for f_0 (980) for
finding its internal configuration; however, uncertainties are too large to
specify the structure at this stage.

We explain the current status of nuclear parton distribution functions in
connection with neutrinonucleus interactions. Neutrino deep inelastic
scattering (DIS) measurements have been done for heavy nuclear targets such as
iron and lead. In order to extract structure functions of the nucleon, one
needs to remove nuclear effects from the data. However, recent studies indicate
that there are inconsistencies in nuclear modifications between chargedlepton
and neutrino scattering measurements. Nuclear medium effects could be also an
origin for the NuTeV anomaly in the weakmixing angle. In addition, the
modifications could affect neutrinooscillation experiments because some DIS
events of neutrinooxygen nucleus interactions are contained. On the other
hand, the nuclear medium effects themselves are interesting and important for
describing nuclei in terms of quark and gluon degrees of freedom.

We investigate impact of $\pi^0$production data at Relativistic Heavy Ion
Collider (RHIC) and future E07011 experiment for the structure function $g_1$
of the deuteron at the Thomas Jefferson National Accelerator Facility (JLab) on
studies of nucleonic spin structure, especially on the polarized gluon
distribution function. By global analyses of polarized leptonnucleon
scattering and the $\pi^0$production data, polarized parton distribution
functions are determined and their uncertainties are estimated by the Hessian
method. Two types of the gluon distribution function are investigated. One is a
positive distribution and the other is a nodetype distribution which changes
sign at $x \sim 0.1$. Although the RHIC $\pi^0$ data seem to favor the node
type for $\Delta g(x)$, it is difficult to determine a precise functional form
from the current data. However, it is interesting to find that the gluon
distribution $\Delta g(x)$ is positive at large $x$ ($>0.2$) due to constraints
from the scaling violation in $g_1$ and RHIC $\pi^0$ data. The JLabE07011
measurements for $g_1^d$ should be also able to reduce the gluon uncertainty,
and the reduction is comparable to the one by RUN5 $\pi^0$production data at
RHIC. The reduction is caused by the error correlation between polarized
antiquark and gluon distributions and by a nexttoleadingorder (NLO) gluonic
effect in the structure function $g_1^d$. We find that the JLabE07011 data
are accurate enough to probe the NLO gluonic term in $g_1$. Both RHIC and JLab
data contribute to better determination of the polarized gluon distribution in
addition to improvement on polarized quark and antiquark distributions.

We discuss internal structure of an exotic hadron by using fragmentation
functions. The fragmentation functions for the f_0(980) meson are obtained by a
global analysis of e^++e^ \to f_0+X data. Quark configuration of the f_0(980)
could be determined by peak positions and second moments of the obtained
fragmentation functions.

It is proposed that fragmentation functions should be used to identify exotic
hadrons. As an example, fragmentation functions of the scalar meson f_0(980)
are investigated. It is pointed out that the second moments and functional
forms of the u and squark fragmentation functions can distinguish the
tetraquark structure from $q\bar q$. By the global analysis of f_0 (980)
production data in electronpositron annihilation, its fragmentation functions
and their uncertainties are determined. It is found that the current available
data are not sufficient to determine its internal structure, while precise data
in future should be able to identify exotic quark configurations.

Nuclear parton distribution functions (NPDFs) are determined by global
analyses of experimental data on structurefunction ratios F_2^A/F_2^{A'} and
DrellYan crosssection ratios \sigma_{DY}^A/\sigma_{DY}^{A'}. The analyses are
done in the leading order (LO) and nexttoleading order (NLO) of running
coupling constant \alpha_s. Uncertainties of the NPDFs are estimated in both LO
and NLO for finding possible NLO improvement. Valencequark distributions are
well determined, and antiquark distributions are also determined at x<0.1.
However, the antiquark distributions have large uncertainties at x>0.2. Gluon
modifications cannot be fixed at this stage. Although the advantage of the NLO
analysis, in comparison with the LO one, is generally the sensitivity to the
gluon distributions, gluon uncertainties are almost the same in the LO and NLO.
It is because current scalingviolation data are not accurate enough to
determine precise nuclear gluon distributions. Modifications of the PDFs in the
deuteron are also discussed by including data on the protondeuteron ratio
F_2^D/F_2^p in the analysis. A code is provided for calculating the NPDFs and
their uncertainties at given x and Q^2 in the LO and NLO.

Nuclear parton distribution functions (NPDFs) are determined by a global
analysis of experimental measurements on structurefunction ratios
F_2^A/F_2^{A'} and DrellYan cross section ratios
\sigma_{DY}^A/\sigma_{DY}^{A'}, and their uncertainties are estimated by the
Hessian method. The NPDFs are obtained in both leading order (LO) and
nexttoleading order (NLO) of \alpha_s. As a result, valencequark
distributions are relatively well determined, whereas antiquark distributions
at x>0.2 and gluon distributions in the whole x region have large
uncertainties. The NLO uncertainties are slightly smaller than the LO ones;
however, such a NLO improvement is not as significant as the nucleonic case.

Fragmentation functions of pion, kaon, and nucleon are determined by global
analyses of hadronproduction data in $e^+e^$ annihilation. It is particularly
important that uncertainties of the fragmentation functions are estimated for
the first time. We found that lightquark and gluon fragmentation functions
have large uncertainties, so that one should be careful in using these
functions for hadronproduction processes in heavyion collisions and lepton
scattering. The analysis is extended to possible exotic hadron search by
fragmentation functions. We found that internal structure of $f_0 (980)$, such
as $s\bar s$ or tetraquark configuration, can be determined by noting
differences between favored and disfavored fragmentation functions.

Fragmentation functions are determined for the pion, kaon, and proton by
analyzing chargedhadron production data in electronpositron annihilation. It
is important that uncertainties of the determined fragmentation functions are
estimated in this analysis. Analysis results indicate that gluon and
lightquark functions have large uncertainties especially at small Q^2. We find
that nexttoleadingorder (NLO) uncertainties are significantly reduced in
comparison with leadingorder (LO) ones in the pion and kaon. The fragmentation
functions are very different in various analysis groups. However, all the
recent functions are roughly within the estimated uncertainties, which
indicates that they are consistent with each other. We provide a code for
calculating the fragmentation functions and their uncertainties at a given
kinematical point of z and Q^2 by a user.

Fragmentation functions and their uncertainties are determined for pion,
kaon, and proton by a global $\chi^2$ analysis of chargedhadron production
data in electronpositron annihilation and by the Hessian method for error
estimation. It is especially important that the uncertainties of the
fragmentation functions are estimated in this analysis. The results indicate
that the fragmentation functions, especially gluon and lightquark
fragmentation functions, have large uncertainties at small $Q^2$. There are
large differences between widelyused functions by KKP (Kniehl, Kramer, and
P\"otter) and Kretzer; however, they are compatible with each other and also
with our functions if the uncertainties are taken into account. We find that
determination of the fragmentation functions is improved in
nexttoleadingorder (NLO) analyses for the pion and kaon in comparison with
leadingorder ones. Such a NLO improvement is not obvious in the proton. Since
the uncertainties are large at small $Q^2$, the uncertainty estimation is very
important for analyzing hadronproduction data at small $Q^2$ or $p_T$ ($Q^2,
p_T^2 << M_Z^2$) in lepton scattering and hadronhadron collisions. A code is
available for general users for calculating obtained fragmentation functions.

We determine the polarized gluon distribution \Delta g(x) by a global
analysis using current DIS and \pi^0 asymmetry data. The \pi^0 data from
RHICSpin experiments provide a strong constraint on \Delta g(x), so that its
uncertainty is reduced. However, a sign problem appears in the analysis using
\pi^0 data, which means that positive and negative distributions are allowed
for \Delta g(x). These two types of solutions are discussed.

Fragmentation functions are determined for pions, kaons, and nucleons by a
global analysis of chargedhadron production data in electronpositron
annihilation. The optimum functions are obtained in both leading order (LO) and
nexttoleading order (NLO) of alpha_s. It is important that uncertainties of
the fragmentation functions are estimated in this work by the Hessian method.
We found that the uncertainties are large at small Q^2 and that they are
generally reduced in the NLO in comparison with the LO ones. We supply a code
for calculating the fragmentation functions and their uncertainties for the
pions, kaons, and nucleons at given z and Q^2.

We investigate the polarized gluon distribution \Delta g(x) by a global
analysis of current DIS data and the \pi^0 data from RHICSpin experiments. The
\pi^0 data provide a strong constraint on \Delta g(x), so that its uncertainty
is reduced. Furthermore, new DIS data of COMPASS and HERMES play an important
role in determining \Delta g(x) at large x.

Global analysis has been performed within the nexttoleading order in
Quantum Chromodynamics (QCD) to determine polarized parton distributions with
new experimental data in spin asymmetries. The new data set includes JLab,
HERMES, and COMPASS measurements on spin asymmetry A_1 for the neutron and
deuteron in lepton scattering. Our new analysis also utilizes the doublespin
asymmetry for pi^0 production in polarized pp collisions, A_{LL}^{pi^0},
measured by the PHENIX collaboration. Because of these new data, uncertainties
of the polarized PDFs are reduced. In particular, the JLab, HERMES, and COMPASS
measurements are valuable for determining Delta d_v(x) at large x and Delta
qbar(x) at x~0.1. The PHENIX pi^0 data significantly reduce the uncertainty of
Delta g(x). Furthermore, we discuss a possible constraint on Delta g(x) at
large x by using the HERMES data on g_1^d in comparison with the COMPASS ones
at x~0.05.