
This paper presents a detailed account of evaluation of the electron
anomalous magnetic moment a_e which arises from the gaugeinvariant set, called
Set V, consisting of 6354 tenthorder Feynman diagrams without closed lepton
loops. The latest value of the sum of Set V diagrams evaluated by the
MonteCarlo integration routine VEGAS is 8.726(336)(\alpha/\pi)^5, which
replaces the very preliminary value reported in 2012. Combining it with other
6318 tenthorder diagrams published previously we obtain
7.795(336)(\alpha/\pi)^5 as the complete massindependent tenthorder term.
Together with the improved value of the eighthorder term this leads to
a_e(theory)=1 159 652 181.643(25)(23)(16)(763) \times 10^{12}, where first
three uncertainties are from the eighthorder term, tenthorder term, and
hadronic and elecroweak terms. The fourth and largest uncertainty is from
\alpha^{1}=137.035 999 049(90), the finestructure constant derived from the
rubidium recoil measurement. Thus, a_e(experiment)  a_e(theory)= 0.91(0.82)
\times 10^{12}. Assuming the validity of the standard model, we obtain the
finestructure constant \alpha^{1}(a_e)=137.035 999 1570(29)(27)(18)(331),
where uncertainties are from the eighthorder term, tenthorder term, hadronic
and electroweak terms, and the measurement of a_e. This is the most precise
value of \alpha available at present and provides a stringent constraint on
possible theories beyond the standard model.

We show experimentally that a stable wave propagating into a region
characterized by an opposite current may become modulationaly unstable.
Experiments have been performed in two independent wave tank facilities; both
of them are equipped with a wavemaker and a pump for generating a current
propagating in the opposite direction with respect to the waves. The
experimental results support a recent conjecture based on a currentmodified
Nonlinear Schr\"odinger equation which establishes that rogue waves can be
triggered by nonhomogeneous current characterized by a negative horizontal
velocity gradient.

This paper reports the tenthorder QED contribution to the lepton g2 from
the gaugeinvariant set, called Set III(c), which consists of 390 Feynman
vertex diagrams containing an internal fourthorder lightbylightscattering
subdiagram. The massindependent contribution of Set III(c) to the electron g2
(a_e) is 4.9210(103) in units of (alpha/pi)^5. The massdependent contributions
to a_e from diagrams containing a muon loop is 0.00370(37) (alpha/pi)^5. The
taulepton loop contribution is negligible at present. Altogether the
contribution of Set III(c) to a_e is 4.9247 (104) (alpha/pi)^5. We have also
evaluated the contribution of the closed electron loop to the muon g2 (a_mu).
The result is 7.435(134) (alpha/pi)^5. The contribution of the taulepton loop
to a_mu is 0.1999(28)(alpha/pi)^5. The total contribution of variousleptonic
loops (electron, muon, and taulepton) of Set III(c) to a_mu is 12.556 (135)
(alpha/pi)^5.

We present a study of the evolution of magnetism from the quantum critical
system YbRh2Si2 to the stable trivalent Yb system YbCo2Si2. Single crystals of
Yb(Rh_(1x)Co_x)2Si2 were grown for 0 < x < 1 and studied by means of magnetic
susceptibility, electrical resistivity, and specific heat measurements, as well
as photoemission spectroscopy. The results evidence a complex magnetic phase
diagram, with a nonmonotonic evolution of T_N and two successive transitions
for some compositions resulting in two tricritical points. The strong
similarity with the phase diagram of YbRh2Si2 under pressure indicates that Co
substitution basically corresponds to the application of positive chemical
pressure. Analysis of the data proves a strong reduction of the Kondo
temperature T_K with increasing Co content, T_K becoming smaller than T_N for x
~ 0.5, implying a strong localization of the 4f electrons. Furthermore,
lowtemperature susceptibility data confirm a competition between ferromagnetic
and antiferromagnetic exchange. The series Yb(Rh_(1x)Co_x)2Si2 provides an
excellent experimental opportunity to gain a deeper understanding of the
magnetism at the quantum critical point in the vicinity of YbRh2Si2 where the
antiferromagnetic phase disappears (T_N=>0).

This paper reports the tenthorder contributions to the g2 of the electron
a_e and those of the muon a_mu from the gaugeinvariant Set II(c), which
consists of 36 Feynman diagrams, and Set II(d), which consists of 180 Feynman
diagrams. Both sets are obtained by insertion of sixthorder
vacuumpolarization diagrams in the fourthorder anomalous magnetic moment. The
massindependent contributions from Set II(c) and Set II(d) are 0.116 489
(32)(alpha/pi)^5 and 0.243 00 (29)(alpha/pi)^5, respectively. The leading
contributions to a_mu, which involve electron loops only, are 3.888 27
(90)(alpha/pi)^5 and 0.4972 (65)(alpha/pi)^5 for Set II(c) and Set II(d),
respectively. The total contributions of the electron, muon, and taulepton
loops to a_e are 0.116 874 (32) (alpha/pi)^5 for Set II(c) and 0.243 10 (29)
(alpha/pi)^5 for Set II(d). The contributions of electron, muon, and taulepton
loops to a_mu are 5.5594 (11) (alpha/pi)^5 for Set II(c) and 0.2465 (65)
(alpha/pi)^5 for Set II(d).

This paper reports the Feynmanparametric representation of the
vacuumpolarization function consisting of 105 Feynman diagrams of the eighth
order, and its contribution to the gaugeinvariant set called Set I(i) of the
tenthorder lepton anomalous magnetic moment. Numerical evaluation of this set
is carried out using FORTRAN codes generated by an automatic code generation
system gencodevpN developed specifically for this purpose. The contribution of
diagrams containing electron loop to the electron g2 is 0.017 47 (11)
(alpha/pi)^5. The contribution of diagrams containing muon loop is 0.000 001 67
(3) (alpha/pi)^5. The contribution of taulepton loop is negligible at present.
The sum of all these terms is 0.017 47 (11) (alpha/pi)^5. The contribution of
diagrams containing electron loop to the muon g2 is 0.087 1 (59) (alpha/pi)^5.
That of taulepton loop is 0.000 237 (1) (alpha/pi)^5. The total contribution
to a_mu, the sum of these terms and the massindependent term, is 0.104 8 (59)
(alpha/pi)^5.

This paper reports the tenthorder QED contribution to lepton g2 from
diagrams of three gaugeinvariant sets VI(d), VI(g), and VI(h), which are
obtained by including various fourthorder radiative corrections to the
sixthorder g2 containing lightbylightscattering subdiagrams. In the case
of electron g2, they consist of 492, 480, and 630 vertex Feynman diagrams,
respectively. The results of numerical integration, including massdependent
terms containing muon loops, are 1.8418(95) (alpha/pi)^5 for the Set VI(d),
1.5918(65) (alpha/pi)^5 for the Set VI(g), and 0.1797(40) (alpha/pi)^5 for the
Set VI(h), respectively. We also report the contributions to the muon g2,
which derive from diagrams containing an electron, muon or tau lepton loop:
Their sums are 5.876(802) (alpha/pi)^5 for the Set VI(d), 5.710(490)
(alpha/pi)^5 for the Set VI(g), and 8.361(232) (alpha/pi)^5 for the Set VI(h),
respectively.

This paper reports the tenthorder QED contribution to the g2 of electron
and muon from two gaugeinvariant sets, Set I(g) and Set I(h). In the case of
electron g2 Set I(g) consists of 9 Feynman diagrams which have a 4thorder
vacuumpolarization loop containing another 4thorder vacuumpolarization loop.
Set I(h) consists of 30 Feynman diagrams which have a proper 6thorder
vacuumpolarization loop containing a secondorder vacuumpolarization loop.
The results of numerical integration, including massdependent terms containing
one closed loop of muon, are 0.028 597 (4) (alpha/pi)^5 for Set I(g) and 0.001
685 (13) (alpha/pi)^5 for Set I(h), respectively. We also report the
contributions of Set I(g) and Set I(h) to the muon anomaly. Diagrams included
are those containing electron, muon, and taulepton loops. Their sums are 2.640
9 (4) (alpha/pi)^5 and 0.564 8 (11) (alpha/pi)^5, respectively. The sum of
contributions of Sets I(g) and I(h) containing only electron loops are in fair
agreement with the recently obtained asymptotic analytic results.

We have evaluated the contribution to the anomalous magnetic moment of the
electron from six tenthorder Feynman diagrams which contain eighthorder
vacuumpolarization function formed by two lightbylight scattering diagrams
connected by three photons. The integrals are constructed by two different
methods. In the first method the subtractive counter terms are used to deal
with ultraviolet (UV) singularities together with the requirement of
gaugeinvariance. In the second method, the WardTakahashi identity is applied
to the lightbylight scattering amplitudes to eliminate UV singularities.
Numerical evaluation confirms that the two methods are consistent with each
other within their numerical uncertainties. Combining the two results
statistically and adding small contribution from the muons and/or tau leptons,
we obtain $ 0.000 399 9 (18) (\alpha/\pi)^5$. We also evaluated the
contribution to the muon $g2$ from the same set of diagrams and found $ 1.263
44 (14) (\alpha/\pi)^5$.

We have carried out a new evaluation of the eighthorder contribution to the
electron g2 using FORTRAN codes generated by an automatic code generator
gencodeN. Comparison of the "new" result with the "old" one has revealed an
inconsistency in the treatment of the infrared divergences in the latter. With
this error corrected we now have two independent determinations of the
eighthorder term. This leads to the revised value 1 159 652 182.79 (7.71) x
10^{12} of the electron g2, where the uncertainty comes mostly from that of
the best nonQED value of the fine structure constant alpha. The new value of
alpha derived from the revised theory and the latest experiment is alpha^{1} =
137.035 999 084 (51) [0.37 ppb], which is about 4.7 ppb smaller than the
previous alpha^{1}.

We present angleresolved photoelectron spectroscopy data probing the
electronic structure of the Ndsubstituted high$T_c$ cuprate
La$_{1.48}$Nd$_{0.4}$Sr$_{0.12}$CuO$_4$ (NdLSCO). Data have been acquired at
low and high photon energies, $h\nu$ = 55 and 500 eV, respectively. Earlier
comparable lowenergy studies of La$_{1.4x}$Nd$_{0.6}$Sr$_{x}$CuO$_4$ ($x =
0.10, 0.12, 0.15$) have shown strongly suppressed photoemission intensity, or
absence thereof, in large parts of the Brillouin zone. Contrary to these
findings we observe spectral weight at all points along the entire Fermi
surface contour at low and high photon energies. No signs of strong charge
modulations are found. At high photon energy, the Fermi surface shows obvious
differences in shape as compared to the lowenergy results presented here and
in similar studies. The observed difference in shape and the high
bulksensitivity at this photon energy suggest intrinsic electronic structure
differences between the surface and bulk regions.

We have developed an efficient algorithm for the subtraction of infrared
divergences that arise in the evaluation of QED corrections to the anomalous
magnetic moment of lepton (g2). By incorporating this new algorithm, we have
extended the automated codegenerating system developed previously to deal with
diagrams without internal lepton loops (called qtype), which produced
convergent integrals when applied to diagrams that have only
ultravioletdivergent subdiagrams of vertex type. The new system produces
finite integrals for all qtype diagrams, including those that contain
selfenergy subdiagrams and thus exhibit infrareddivergent behavior. We have
thus far verified the system for the sixth and eighthorder cases. We are now
evaluating 6354 vertex diagrams of qtype that contribute to the tenthorder
lepton g2.

The contribution to the eighthorder anomalous magnetic moment (g2) of the
electron from a set of diagrams without closed lepton loops is recalculated
using a new FORTRAN code generated by an automatic code generator. Comparing
the contributions of individual diagrams of old and new calculations, we found
an inconsistency in the old treatment of infrared subtraction terms in two
diagrams. Correcting this error leads to the revised value 1.9144 (35)
(alpha/pi)^4 for the eighthorder term. This theoretical change induces the
shift of the inverse of the fine structure constant by 6.41180(73)x10^{7}.

Kaon photoproduction process on $^{12}$C has been studied by measuring
neutral kaons in a photon energy range of 0.8$$1.1 GeV. Neutral kaons were
identified by the invariant mass constructed from two charged pions emitted in
the $K^{0}_{S}\to\pi^{+}\pi^{}$ decay channel. The differential cross sections
as well as the integrated ones in the threshold photon energy region were
obtained. The obtained momentum spectra were compared with a Spectator model
calculation using elementary amplitudes of kaon photoproduction given by
recent isobar models. Present result provides, for the first time, the
information on $n(\gamma,K^{0})\Lambda$ reaction which is expected to play an
important role to construct models for strangeness production by the
electromagnetic interaction. Experimental results show that cross section of
$^{12}{\rm C}(\gamma,K^0)$ is of the same order to that of $^{12}{\rm
C}(\gamma,K^+)$ and suggest that slightly backward $K^0$ angular distribution
is favored in the $\gamma n\to K^0\Lambda$ process.

The QED contribution to the anomalous magnetic moments of electron and muon
are known very precisely up to the order $\alpha^4$. However, the knowledge of
$\alpha^5$ term will also be required when the precision of measurement
improves further. This paper reports the first systematic attempt to evaluate
the $\alpha^5$ term. Feynman diagrams contributing to this term can be
classified into six gaugeinvariant sets which can be subdivided further into
32 gaugeinvariant subsets. Thus far we have numerically evaluated all
integrals of 17 gaugeinvariant subsets which contain lightbylightscattering
subdiagrams and/or vacuumpolarization subdiagrams. They cover most of leading
terms of muon $g2$ and lead to a preliminary result 663(20) $(\alpha/\pi)^5$,
which is 8.5 times more precise than the old estimate.

This article reports an automated approach to the evaluation of higherorder
terms of QED perturbation to anomalous magnetic moments of charged leptons by
numerical means. We apply this approach to tenthorder correction due to a
particular subcollection of Feynman diagrams, which have no virtual lepton
loops. This set of diagrams is distinctive in that it grows factorially in
number as the order increases, and also each of the diagrams holds quite a
large number of subtraction terms to be treated along renormalization
procedure. Thus some automated scheme has long been required to evaluate
correctly this class of diagrams. We developed a fast algorithm and an
implementation which automates necessary steps to generate from the
representation of each Feynman diagram the FORTRAN codes for numerical
integration. Currently those diagrams of tenth order are being evaluated.

We report a new value of electron $g2$, or $a_e$, from 891 Feynman diagrams
of order $\alpha^4$. The FORTRAN codes of 373 diagrams containing closed
electron loops have been verified by at least two independent formulations. For
the remaining 518 diagrams, which have no closed lepton loop, verification by a
second formulation is not yet attempted because of the enormous amount of
additional work required. However, these integrals have structures that allow
extensive crosschecking as well as detailed comparison with lowerorder
diagrams through the renormalization procedure. No algebraic error has been
uncovered for them. The numerical evaluation of the entire $\alpha^4$ term by
the integration routine VEGAS gives $1.7283 (35) (\alpha/\pi)^4$, where the
uncertainty is obtained by careful examination of error estimates by VEGAS.
This leads to $a_e = 1 159 652 175.86 (0.10) (0.26) (8.48) \times 10^{12}$,
where the uncertainties come from the $\alpha^4$ term, the estimated
uncertainty of $\alpha^5$ term, and the inverse fine structure constant,
$\alpha^{1} = 137.036 000 3 (10)$, measured by atom interferometry combined
with a frequency comb technique, respectively. The inverse fine structure
constant $\alpha^{1} (a_e)$ derived from the theory and the Seattle
measurement of $a_e$ is $137.035 998 83 (51)$.

The $\eta$meson photoproduction cross sections have been measured on the C
and Cu targets for the photon energies between 600 and 1100 MeV to investigate
the behavior of the $S_{11}$(1535) resonance in a nucleus. The excitation
functions of the cross section as well as the angular and momentum
distributions of $\eta$mesons are in quantitative agreement with the Quantum
Molecular Dynamics (QMD) model calculations in which the $\eta$meson emission
processes other than the $S_{11}$(1535) resonance are also incorporated as
proposed in the $\eta$MAID model. It is shown that the excitation of the
$D_{15}$(1675) resonance might play an important role for $E_{\gamma}>900$ MeV.

We have completed the evaluation of all massdependent $\alpha^4$ QED
contributions to the muon $g2$, or $a_\mu$, in two or more different
formulations. Their numerical values have been greatly improved by an extensive
computer calculation. The new value of the dominant $\alpha^4$ term $A_2^{(8)}
(m_\mu / m_e) $ is 132.6823 (72), which supersedes the old value 127.50 (41).
The new value of the threemass term $A_3^{(8)} (m_\mu / m_e, m_\mu / m_\tau) $
is 0.0376 (1). The term $A_2^{(8)} (m_\mu / m_\tau) $ is crudely estimated to
be about 0.005 and may be ignored for now. The total QED contribution to
$a_\mu$ is $116 584 719.58 (0.02)(1.15)(0.85) \times 10^{11}$, where 0.02 and
1.15 are uncertainties in the $\alpha^4$ and $\alpha^5$ terms and 0.85 is from
the uncertainty in $\alpha$ measured by atom interferometry. This raises the
Standard Model prediction by $13.9 \times 10^{11}$, or about 1/5 of the
measurement uncertainty of $a_\mu$. It is within the noise of current
uncertainty ($\sim 100 \times 10^{11}$) in the estimated hadronic
contributions to $a_\mu$.

The $\alpha^4$ contribution to the lepton $g2$ from a gaugeinvariant set of
18 Feynman diagrams containing a lightbylight scattering subdiagram
internally has been reevaluated by a method independent of the previous
approach. Comparison of two methods revealed a program error in the first
version. Correcting this error, the contributions of these 18 diagrams become
0.990 72 (10)$(\alpha/\pi)^4$ and 4.432 43 (58)$(\alpha/\pi)^4$ for the
electron and muon $g2$, respectively. The correction is not large enough to
affect the comparison between theory and experiment for the muon $g2$, but it
does alter the inferred value for the fine structure constant $\alpha^{1}$ by
6 ppb.

The quantum oscillatory effect and superconductivity in a
nontetrachalcogenafulvalene (TCF) structure based organic superconductor
beta(BDATTP)_2SbF_6 are studied. Here the Shubnikovde Haas effect (SdH) and
angular dependent magnetoresistance oscillations (AMRO) are observed. The
oscillation frequency associated with a cylindrical Fermi surface is found to
be about 4050 tesla, which is also verified by the tunnel diode oscillator
(TDO) measurement. The upper critical field Hc2 measurement in a tilted
magnetic field and the TDO measurement in the mixed state reveal a highly
anisotropic superconducting nature in this material. We compared physical
properties of beta(BDATTP)_2SbF_6 with typical TCF structure based quasi
twodimensional organic conductors. A notable feature of beta(BDATTP)_2SbF_6
superconductor is a large value of effective cyclotron mass m_c^*=12.4+/1.1
m_e, which is the largest yet found in an organic superconductor. A possible
origin of the enhanced effective mass and its relation to the superconductivity
are briefly discussed.

The contribution of the $\alpha^3$ single electronloop vacuumpolarization
diagrams to the Lamb shift of the muonic hydrogen has been evaluated recently
by two independent methods. One uses the exact parametric representation of the
vacuumpolarization function while the other relies on the Pad\'{e}
approximation method. High precision of these values offers an opportunity to
examine the reliability of the MonteCarlo integration as well as that of the
Pad\'{e} method. Our examination covers both muonic hydrogen atom and muon
$g2$. We tested them further for the cases involving twoloop vacuum
polarization, where an exact analytic result is known. Our analysis justifies
the result for the Lamb shift of the muonic hydrogen and also resolves the
longstanding discrepancy between two previous evaluations of the muon $g2$
value.

The sixthorder electronloop vacuumpolarization contribution to the
$2P_{1/2}  2S_{1/2}$ Lamb shift of the muonic hydrogen ($\mu^{} p^+$ bound
state) has been evaluated numerically. Our result is 0.007608(1) meV. This
eliminates the largest uncertainty in the theoretical calculation. Combined
with the proposed precision measurement of the Lamb shift it will lead to a
very precise determination of the proton charge radius.

The method of NRQED is an adaptation of QED to bound systems. While it is
fully equivalent to QED, it enables us to explore the QED bound states
systematically as an expansion in the fine structure constant $\alpha$ and
velocity ($\sim Z\alpha$) of the bound electron. I describe how to construct
the NRQED Hamiltonian choosing recent works on the $\alpha (Z\alpha)$,
$\alpha^2 (Z\alpha)$,$\alpha (Z\alpha)^2$, $\alpha (Z\alpha)^3$, and $\alpha^2
(Z\alpha)^2$ radiative corrections to the hyperfine structure of the muonium
ground state as examples.

This is the second of a series of papers on the radiative corrections of
order $\alpha^2 (Z\alpha)$, $\alpha (Z\alpha )^2$, and various logarithmic
terms of order $\alpha^4$, to the hyperfine structure of the muonium ground
state. This paper deals with the $\alpha (Z\alpha)^2$ correction. Based on the
NRQED bound state theory, we isolated the term of order $\alpha(Z\alpha)^2$
exactly. Our result $+16.904 2 (11) \alpha(Z\alpha)^2 E_F / \pi$ for the
nonlogarithmic part is consistent with the $\alpha (Z\alpha )^2$ part of
Sapirstein's calculation and the recent result of Pachucki, and reduces the
numerical uncertainty in the $\alpha (Z\alpha)^2$ term by two orders of
magnitude.