• A precise description of neutrino-nucleus 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 next-generation neutrino oscillation experiments. The neutrino energy relevant to the neutrino-nucleus reactions spans a broad range and, accordingly, the dominant reaction mechanism varies across the energy region from quasi-elastic 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 meson-baryon to quarks. The main purpose of this review is to report our recent efforts towards a unified description of the neutrino-nucleus reactions over the wide energy range; recent overall progress in the field is also sketched. Starting with an overview of the current status of neutrino-nucleus scattering experiments, we formulate the cross section to be commonly used for the reactions over all the energy regions. A description of the neutrino-nucleon reactions follows and, in particular, a dynamical coupled-channels model for meson productions in and beyond the $\Delta$(1232) region is discussed in detail. We then discuss the neutrino-nucleus reactions, putting emphasis on our theoretical approaches. We start the discussion with electroweak processes in few-nucleon systems studied with the correlated Gaussian method. Then we describe quasi-elastic 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 charged-hadron production data in electron-positron annihilation. Accurate measurements were reported by the Belle and BaBar collaborations for the fragmentation functions at the center-of-mass 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 B-factory 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 high-statistical 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 light-quark 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 single-inclusive $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 charged-lepton 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 short-range nuclear correlation which is caused by high densities due to cluster formation.
  • Next generation neutrino oscillation experiments will need a quantitative understanding of neutrino-nucleus interaction far better than ever. Kinematics covered by the relevant neutrino-nucleus interaction spans wide region, from the quasi-elastic, through the resonance region, to the deeply inelastic scattering region. The neutrino-nucleus 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 neutrino-nucleus 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.
  • Semi-inclusive hadron-production processes are becoming important in high-energy hadron reactions. They are used for investigating properties of quark-hadron matters in heavy-ion collisions, for finding the origin of nucleon spin in polarized lepton-nucleon and nucleon-nucleon reactions, and possibly for finding exotic hadrons. In describing the hadron-production cross sections in high-energy 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 (Dokshitzer-Gribov-Lipatov-Altarelli-Parisi) evolution equations, which are often used in theoretical and experimental analyses of the fragmentation functions and in calculating semi-inclusive cross sections. The DGLAP equations are complicated integro-differential equations, which cannot be solved in an analytical method. In this work, a simple method is employed for solving the evolution equations by using Gauss-Legendre 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 next-to-leading 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 Be-9 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 Be-9 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, tensor-polarized 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 tensor-structure study in 2010's. Furthermore, the antiquark tensor polarization could be measured by polarized deuteron Drell-Yan processes at hadron facilities such as J-PARC and GSI-FAIR. Third, the recent CDF dijet anomaly is investigated within the standard model by considering possible modifications of the strange-quark distribution. We find that the shape of a dijet-mass spectrum changes depending on the strange-quark 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 high-energy 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 long-standing NuTeV sin^2 theta_w anomaly. There is also an issue of nuclear modification differences between the structure functions of charged-lepton 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 Be-9 nucleus at the Thomas Jefferson National Accelerator Facility (JLab). Last, tensor-polarized 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 spin-one hadrons.
  • An anomalous nuclear modification was reported by JLab measurements on the beryllium-9 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 Be-9 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 Be-9. The JLab measurement is possibly the first signature of clustering effects in high-energy 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 beryllium-9 nucleus at the Thomas Jefferson National Accelerator Facility (JLab), clustering aspects are studied in structure functions of deep inelastic lepton-nucleus 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 Be-9 nucleus consists of two alpha-like clusters with a surrounding neutron. The clustering produces high-momentum 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 Be-9 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 Be-9. It indicates that the anomalous Be-9 result should be explained by a different mechanism from the nuclear binding and Fermi motion. If nuclear-modification slopes d(F_2^A/F_2^D)/dx are shown by the maximum local densities, the Be-9 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 Be-9 should be explained by large densities in the clusters. For example, internal nucleon structure could be modified in the high-density 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 K-doped 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 K-doped picene has three structures similar to those of pristine picene with new states near EF, consistent with the semiconductor-metal transition. The K-induced change cannot be explained with a simple rigid-band model of pristine picene, but can be interpreted by molecular orbital calculations considering electron-intramolecular-vibration interaction. Excellent agreement of the K-doped spectrum with the calculations points to importance of electron-intramolecular-vibration interaction in K-doped picene.
  • We discuss studies on determination of fragmentation functions and an application to exotic-hadron 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 next-to-leading order (NLO) of the running coupling constant alpha_s by global analyses of hadron-production data in electron-positron annihilation. Various parametrization results are much different in disfavored-quark 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 disfavored-quark 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 neutrino-nucleus 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 charged-lepton and neutrino scattering measurements. Nuclear medium effects could be also an origin for the NuTeV anomaly in the weak-mixing angle. In addition, the modifications could affect neutrino-oscillation experiments because some DIS events of neutrino-oxygen 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 E07-011 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 lepton-nucleon 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 node-type 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 JLab-E07-011 measurements for $g_1^d$ should be also able to reduce the gluon uncertainty, and the reduction is comparable to the one by RUN-5 $\pi^0$-production data at RHIC. The reduction is caused by the error correlation between polarized antiquark and gluon distributions and by a next-to-leading-order (NLO) gluonic effect in the structure function $g_1^d$. We find that the JLab-E07-011 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 s-quark fragmentation functions can distinguish the tetraquark structure from $q\bar q$. By the global analysis of f_0 (980) production data in electron-positron 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 structure-function ratios F_2^A/F_2^{A'} and Drell-Yan cross-section ratios \sigma_{DY}^A/\sigma_{DY}^{A'}. The analyses are done in the leading order (LO) and next-to-leading order (NLO) of running coupling constant \alpha_s. Uncertainties of the NPDFs are estimated in both LO and NLO for finding possible NLO improvement. Valence-quark 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 scaling-violation 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 proton-deuteron 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 structure-function ratios F_2^A/F_2^{A'} and Drell-Yan 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 next-to-leading order (NLO) of \alpha_s. As a result, valence-quark 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 hadron-production data in $e^+e^-$ annihilation. It is particularly important that uncertainties of the fragmentation functions are estimated for the first time. We found that light-quark and gluon fragmentation functions have large uncertainties, so that one should be careful in using these functions for hadron-production processes in heavy-ion 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 charged-hadron production data in electron-positron annihilation. It is important that uncertainties of the determined fragmentation functions are estimated in this analysis. Analysis results indicate that gluon and light-quark functions have large uncertainties especially at small Q^2. We find that next-to-leading-order (NLO) uncertainties are significantly reduced in comparison with leading-order (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 charged-hadron production data in electron-positron 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 light-quark fragmentation functions, have large uncertainties at small $Q^2$. There are large differences between widely-used 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 next-to-leading-order (NLO) analyses for the pion and kaon in comparison with leading-order 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 hadron-production data at small $Q^2$ or $p_T$ ($Q^2, p_T^2 << M_Z^2$) in lepton scattering and hadron-hadron 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 RHIC-Spin 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 charged-hadron production data in electron-positron annihilation. The optimum functions are obtained in both leading order (LO) and next-to-leading 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 RHIC-Spin 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 next-to-leading 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 double-spin 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.