
We present a new magnetohydrodynamic (MHD) simulation code with the aim of
providing accurate numerical solutions to astrophysical phenomena where
discontinuities, shock waves, and turbulence are inherently important. The code
implements the HLLD approximate Riemann solver, the
fifthordermonotonicitypreserving interpolation (MP5) scheme, and the
hyperbolic divergence cleaning method for a magnetic field. This choice of
schemes significantly improved numerical accuracy and stability, and saved
computational costs in multidimensional problems. Numerical tests of one and
twodimensional problems showed the advantages of using the highorder scheme
by comparing with results from a standard secondorder TVD MUSCL scheme. The
present code enabled us to explore longterm evolution of a threedimensional
accretion disk around a black hole, in which compressible MHD turbulence caused
continuous mass accretion via nonlinear growth of the magnetorotational
instability (MRI). Numerical tests with various computational cell sizes
exhibited a convergent picture of the early nonlinear growth of the MRI in a
global model, and indicated that the MP5 scheme has more than twice the
resolution of the MUSCL scheme in practical applications.

We report the discovery of giant and anisotropic magnetoresistance due to the
orbital rearrangement in a nonmagnetic correlated metal. In particular, we
measured the magnetoresistance under fields up to 31.4 T in the cubic Prbased
heavy fermion superconductor PrV$_2$Al$_{20}$ with a nonmagnetic $\Gamma _3$
doublet ground state, exhibiting antiferroquadrupole ordering below 0.7 K. For
the [100] direction, we find that the highfield phase appears between 12 T and
25 T, accompanied by a large jump at 12 T in the magnetoresistance ($\Delta MR
\sim $ 100 $\% $) and in the anisotropic magnetoresistivity (AMR) ratio by
$\sim $ 20 $\% $. These observations indicate that the strong hybridization
between the conduction electrons and anisotropic quadrupole moments leads to
the Fermi surface reconstruction upon crossing the fieldinduced
antiferroquadrupole (orbital) rearrangement.

We investigated the precursor wave emission efficiency in magnetized purely
perpendicular relativistic shocks in pair plasmas. We extended our previous
study to include the dependence of upstream magnetic field orientations. We
performed twodimensional particleincell simulations and focused on two
magnetic field orientations: the magnetic field to be in the simulation plane
(i.e., inplane configuration) and perpendicular to the simulation plane (i.e.,
outofplane configuration). Our simulations in the inplane configuration
demonstrated that not only extraordinary but also ordinary mode waves are
excited. We quantified the emission efficiency as a function of the
magnetization parameter $\sigma_e$ and found that the largeamplitude precursor
waves are emitted for a wide range of $\sigma_e$. We found that especially at
low $\sigma_e$, the magnetic field generated by Weibel instability amplifies
the ordinary mode wave power. The amplitude is large enough to perturb the
upstream plasma, and transverse density filaments are generated as in the case
of the outofplane configuration investigated in the previous study. We
confirmed that our previous conclusion holds regardless of upstream magnetic
field orientations with respect to the twodimensional simulation plane. We
discuss the precursor wave emission in three dimensions and the feasibility of
wakefield acceleration in relativistic shocks based on our results.

Conventionally, the parallelism between the physical properties of Ce and Yb
based magnets and heavy fermions has been considered to arise as a consequence
of the electronhole symmetry known for the band theory in condensed matter
physics. Namely, the trivalent state of Ce has $4f^{1}$ configuration, while
$4f^{13}$ configuration for Yb$^{3+}$ can be viewed as a state carrying one
4$f$ hole. However, in the light of this parallelism, the magnetism of Yb based
intermetallic compounds has been a mystery over decades; the transition
temperature of the Yb based compounds is normally very small, as low as $\sim$
1 K or even lower, while Ce counterparts may often have the transition
temperature well exceeding 10 K. Here, we report our experimental discovery of
the transition temperature reaching 20 K for the first time in an Yb based
compound. Our single crystal study using transport, thermomagnetic and
photoemission spectroscopy measurements reveals that the Mn substitution at the
Al site in an intermediate valence state of $\alpha$YbAlB$_{4}$ not only
induces antiferromagnetic transition at a record high temperature of 20 K but
also transform the heavy fermion liquid state in $\alpha$YbAlB$_{4}$ into a
highly resistive metallic state proximate to a Kondo insulator. The high
temperature antiferromagnetism in a system close to a Kondo insulator is
unusual and suggest a novel type of itinerant magnetism in the Yb based heavy
fermion systems.

How electrons get accelerated to relativistic energies in a highMachnumber
quasiperpendicular shock is presented by means of ab initio particleincell
simulations in three dimensions. We found that coherent electrostatic Buneman
waves and ionWeibel magnetic turbulence coexist in a strongshock structure
whereby particles gain energy during shocksurfing and subsequent stochastic
drift accelerations. Energetic electrons that initially experienced the surfing
acceleration undergo pitchangle diffusion by interacting with magnetic
turbulence and continuous acceleration during confinement in the shock
transition region. The ionWeibel turbulence is the key to the efficient
nonthermal electron acceleration.

We investigated the efficiency of coherent upstream largeamplitude
electromagnetic wave emission via synchrotron maser instability at relativistic
magnetized shocks by using twodimensional particleincell simulations. We
considered the purely perpendicular shock in an electronpositron plasma. The
coherent wave emission efficiency was measured as a function of the
magnetization parameter {\sigma}, which is defined by the ratio of the Poynting
flux to the kinetic energy flux. The wave amplitude was systematically smaller
than that observed in onedimensional simulations. However, it continued to
persist, even at a considerably low magnetization rate, where the Weibel
instability dominated the shock transition. The emitted electromagnetic waves
were sufficiently strong to disturb the upstream medium, and transverse
filamentary density structures of substantial amplitude were produced. Based on
this result, we discuss the possibility of the wakefield acceleration model for
the production of nonthermal electrons in a relativistic magnetized
ionelectron shock.

A distinct thermal Hall signal is observed in a quantum spin liquid candidate
Ba$_3$CuSb$_2$O$_9$. The transverse thermal conduction shows a powerlaw
temperature dependence below 40 K where a spin gap opens. We suggest that,
through the very low longitudinal thermal conductivity and the thermal Hall
signals, a phonon Hall effect is induced by strong phonon scatterings by orphan
Cu$^{2+}$ spins formed in random domains of Cu$^{2+}$Sb$^{5+}$ dumbbells.

A possibility of timedelayed radio brightenings of Sgr A* triggered by the
pericenter passage of the G2 cloud is studied by carrying out global
threedimensional magnetohydrodynamic simulations taking into account the
radiative cooling of the tidal debris of the G2 cloud. Magnetic fields in the
accretion flow are strongly perturbed and reorganized after the passage of G2.
We have found that the magnetic energy in the accretion flow increases by a
factor 34 in 510 years after the pericenter passage of G2 by a dynamo
mechanism driven by the magnetorotational instability. Since this Bfield
amplification enhances the synchrotron emission from the disk and the outflow,
the radio and the infrared luminosity of Sgr A* is expected to increase around
A.D. 2020. The timedelay of the radio brightening enables us to determine the
rotation axis of the preexisting disk.

We report the results of 27AlNMR measurements on a single crystal of
PrTi2Al20. This compound shows a phase transition near 2 K associated with the
quadrupole degree of freedom of Pr ions, for which the ground state of the
crystalline electric field is a nonmagnetic doublet. When a magnetic field is
applied along the <111> direction, the NMR lines split upon entering the
lowtemperature phase, indicating the breaking of the three fold rotation
symmetry due to the fieldinduced magnetic dipole perpendicular to the field.
This provides microscopic evidence for a ferro order of the O20 quadrupole.
Although a firstorder transition is expected theoretically, the line splitting
evolves continuously within the experimental resolution.

We report on the evidence for the multiband electronic transport in
$\alpha$YbAlB$_{4}$ and $\alpha$Yb$_{0.81(2)}$Sr$_{0.19(3)}$AlB$_{4}$.
Multiband transport reveals itself below 10 K in both compounds via Hall effect
measurements, whereas anisotropic magnetic ground state sets in below 3 K in
$\alpha$Yb$_{0.81(2)}$Sr$_{0.19(3)}$AlB$_{4}$. Our results show that Sr$^{2+}$
substitution enhances conductivity, but does not change the quasiparticle mass
of bands induced by heavy fermion hybridization.

Specific heat, elastic neutron scattering, and muon spin rotation ($\mu$SR)
experiments have been carried out on a wellcharacterized sample of "stuffed"
(Prrich) Pr$_{2+x}$Ir$_{2x}$O$_{7\delta}$. Elastic neutron scattering shows
the onset of longrange spinice "2in/2out" magnetic order at $T_M = 0.93$ K,
with an ordered moment of 1.7(1)$\mu_\mathrm{B}$/Pr ion at low temperatures.
Approximate lower bounds on the correlation length and correlation time in the
ordered state are 170 \AA\ and 0.7 ns, respectively. $\mu$SR experiments yield
an upper bound 2.6(7) mT on the local field $B_\mathrm{loc}^{4f}$ at the muon
site, which is nearly two orders of magnitude smaller than the expected dipolar
field for longrange spinice ordering of 1.7$\mu_B$ moments (120270 mT,
depending on muon site). This shortfall is due in part to splitting of the
nonKramers crystalfield groundstate doublets of nearneighbor Pr$^{3+}$ ions
by the $\mu^+$induced lattice distortion. For this to be the only effect,
however, $\sim$160 Pr moments out to a distance of $\sim$14 \AA\ must be
suppressed. An alternative scenario, which is consistent with the observed
reduced nuclear hyperfine Schottky anomaly in the specific heat, invokes slow
correlated Prmoment fluctuations in the ordered state that average
$B_\mathrm{loc}^{4f}$ on the $\mu$SR time scale (${\sim}10^{7}$ s), but are
static on the time scale of the elastic neutron scattering experiments
(${\sim}10^{9}$ s). In this picture the dynamic muon relaxation suggests a
Pr$^{3+}$ $4f$ correlation time of a few nanoseconds, which should be
observable in a neutron spin echo experiment.

We studied the stability property of numerical Cherenkov radiation in
relativistic plasma flows employing particleincell simulations. Using the
implicit finitedifference timedomain method to solve Maxwell equations, we
found that nonphysical instability was greatly inhibited with a
CourantFriedrichsLewy (CFL) number of 1.0. The present result contrasts with
recently reported results (Vay, J. L., et al. 2011, J. Comp. Phys, 230, 5908;
Godfrey. B., & Vay, J. L. 2013, J. Comp. Phys, 243, 260; Xu, X., et al. 2013,
Comput. Phys. Commun, 184, 2503) in which magical CFL numbers in the range
0.50.7 were obtained with explicit field solvers. In addition, we found
employing higherorder shape functions and an optimal implicitness factor
further suppressed longwavelength modes of the instability. The findings
allowed the examination of the longterm evolution of a relativistic
collisionless shock without the generation of nonphysical wave excitations in
the upstream. This achievement will allow us to investigate particle
accelerations in relativistic shocks associated with, for example, gammaray
bursts.

We have measured the specific heat of PrV$_2$Al$_{20}$ at very low
temperatures, using high quality single crystals with the residual resistivity
ratio ~ 20. The highquality single crystals exhibit clear double transitions
at $T_{\rm Q} =$ 0.75 K and $T^* =$ 0.65 K. These transitions are clearer and
shift to higher temperature in higher quality single crystals. Besides, there
was no hysteresis in those transitions in warming and cooling process of the
heat capacity measurements. In the ordered state below $T^*$, the specific heat
shows a $T^4$ dependence, indicating the gapless mode associated with the
quadrupole and/or octupole ordering.

We measured the transverse magnetoresistivity of the mixed valence compound
$\alpha$YbAlB$_4$. Two configurations were used where current was applied
along [110] direction for both and magnetic field was applied along [110] and
$c$axis. We found the transverse magnetoresistivity is highly anisotropic. In
the weak field below 1 T, it is consistent with stronger $c$$f$ hybridization
in the $ab$ plane which was suggested from the previous zero field resistivity
measurements. At the higher field above 3 T, we observed a negative transverse
magnetoresistivity for the field applied along the $c$axis. The temperature
dependences of the resistivity measured at several different fields suggest the
suppression of the heavy fermion behavior at the characteristic field of $\sim
5.5$ T.

We present the highprecision magnetization data of the valence fluctuating
heavy fermion superconductor $\beta$YbAlB$_4$ in a wide temperature range from
0.02 K to 320 K spanning four orders of magnitude. We made detailed analyses of
the $T/B$ scaling of the magnetization, and firmly confirmed the unconventional
zerofield quantum criticality (QC) without tuning. We examined other possible
scaling relationship such as $T/(BB_c)^{\delta}$ scaling, and confirmed that
$\delta = 1$ provides the best quality of the fit with an upper bound on the
critical magnetic field $\vert B_c \vert <0.2$~mT. We further discuss the heavy
Fermiliquid component of the magnetization after subtracting the QC component
estimated based on the $T/B$ scaling. The temperature dependence of the heavy
Fermiliquid component is found very similar to the magnetization of the
polymorph $\alpha$YbAlB$_4$. In addition, the heavy Fermiliquid component is
suppressed in the magnetic field above $\sim$ 5 T as in $\alpha$YbAlB$_4$.
This was also confirmed by the magnetization measurements up to $\sim 50$ T for
both $\alpha$ and $\beta$YbAlB$_4$. Interestingly, the detailed analyses
revealed that the only a part of $f$ electrons participates in the zerofield
QC and the heavy fermion behavior. We also present a temperature  magnetic
field phase diagram of \ybal to illustrate how the characteristic temperature
and field scales evolves near the QC.

PrV$_2$Al$_{20}$ is a rare example of a heavy fermion system based on strong
hybridization between conduction electrons and nonmagnetic quadrupolar moments
of the cubic $\Gamma_3$ ground doublet. Here, we report that a highquality
single crystal of PrV$_2$Al$_{20}$ exhibits superconductivity at $T_{\rm c}=$
50 mK in the antiferroquadrupoleordered state under ambient pressure. The
heavy fermion character of the superconductivity is evident from the specific
heat jump of $\Delta C/T \sim 0.3$ J/mol K$^2$ and the effective mass $m^*/m_0
\sim 140$ estimated from the temperature dependence of the upper critical
field. Furthermore, the highquality single crystals exhibit double transitions
at $T_{\rm Q}$ = 0.75 K and $T^{*}$ = 0.65 K associated with quadrupole and
octapole degrees of freedom of the $\Gamma_3$ doublet. In the ordered state,
the specific heat $C/T$ shows a $T^3$ dependence, indicating the gapless mode
associated with the quadrupole and/or octapole order. The strong sensitivity to
impurity of the superconductivity suggests unconventional character due to
significant quadrupolar fluctuations.

Here we provide the first clear evidence of Fermiliquid breakdown in an
intermediate valence system. We employ high precision magnetization
measurements of the valence fluctuating superconductor betaYbAlB4 to probe the
quantum critical free energy down to temperatures far below the characteristic
energy scale of the valence fluctuations. The observed T/B scaling in the
magnetization over three decades not only indicates unconventional quantum
criticality, but places an upper bound on the critical magnetic field B_c <
0.2 mT, a value comparable with the Earth's magnetic field and six orders of
magnitude smaller than the valence fluctuation scale. This tiny value of the
upper bound on B_c, well inside the superconducting dome, raises the
fascinating possibility that valence fluctuating betaYbAlB4 is intrinsically
quantum critical, without tuning the magnetic field, pressure, or composition:
the first known example of such a phenomenon in a metal.

We present an extension of the multimoment advection scheme (Minoshima et
al., 2011, J. Comput. Phys.) to the threedimensional case, for full
electromagnetic Vlasov simulations of magnetized plasma. The scheme treats not
only point values of a profile but also its zeroth to second order piecewise
moments as dependent variables, and advances them on the basis of their
governing equations. Similar to the twodimensional scheme, the
threedimensional scheme can accurately solve the solid body rotation problem
of a gaussian profile with little numerical dispersion or diffusion. This is a
very important property for Vlasov simulations of magnetized plasma. We apply
the scheme to electromagnetic Vlasov simulations. Propagation of linear waves
and nonlinear evolution of the electron temperature anisotropy instability are
successfully simulated with a good accuracy of the energy conservation.

Hydrodynamical instabilities are believed to power some of the small scale
(0.110 pc) turbulence and chemical mixing in the interstellar medium.
Identifying such instabilities has always been difficult but recent
observations of a wavelike structure (the Ripples) in the Orion nebula have
been interpreted as a signature of the KelvinHelmholtz instability (KHI),
occurring at the interface between the HII region and the molecular cloud.
However, this has not been verified theoretically. In this letter, we
investigate theoretically the stability of this interface using observational
constraints for the local physical conditions. A linear analysis shows that the
HII/molecular cloud interface is indeed KH unstable for a certain range of
magnetic field orientation. We find that the maximal growthrates correspond to
typical timescales of a few 1e4 years and instability wavelengths of 0.06 to
0.6 pc. We predict that after 2e5 years the KHI saturates and forms a turbulent
layer of thickness ~0.5 pc. The KHI can remain in linear phase over a maximum
distance of 0.75 pc. These spatial and time scales are compatible with the
Ripples representing the linear phase of the KHI. These results suggest that
the KHI may be crucial to generate turbulence and to bring heavy elements
injected by the winds of massive stars in HII regions to colder regions where
planetary systems around low mass stars are being formed. This could apply to
the transport of 26Al injected by a massive star in an HII region to the
nascent solarsystem.

Electron accelerations at high Mach number collisionless shocks are
investigated by means of twodimensional electromagnetic ParticleinCell
simulations with various Alfven Mach numbers, iontoelectron mass ratios, and
the upstream electron beta_e (the ratio of the thermal pressure to the magnetic
pressure). We found electrons are effectively accelerated at a superhigh Mach
number shock (MA~30) with a mass ratio of M/m=100 and beta_e=0.5. The electron
shock surfing acceleration is an effective mechanism for accelerating the
particles toward the relativistic regime even in two dimensions with the large
mass ratio. Buneman instability excited at the leading edge of the foot in the
superhigh Mach number shock results in a coherent electrostatic potential
structure. While multidimensionality allows the electrons to escape from the
trapping region, they can interact with the strong electrostatic field several
times. Simulation runs in various parameter regimes indicate that the electron
shock surfing acceleration is an effective mechanism for producing relativistic
particles in extremelyhigh Mach number shocks in supernova remnants, provided
that the upstream electron temperature is reasonably low.

We present a new numerical scheme for solving the advection equation and its
application to the Vlasov simulation. The scheme treats not only point values
of a profile but also its zeroth to second order piecewise moments as dependent
variables, and advances them on the basis of their governing equations. We have
developed one and twodimensional schemes and show that they provide quite
accurate solutions compared to other existing schemes with the same memory
usage. The twodimensional scheme can solve the solid body rotation problem of
a gaussian profile with little numerical diffusion. This is a very important
property for Vlasov simulations of magnetized plasma. The application of the
scheme to the electromagnetic Vlasov simulation of collisionless shock waves is
presented as a benchmark test.

We present a new numerical scheme for solving the advection equation and its
application to Vlasov simulations. The scheme treats not only point values of a
profile but also its zeroth to second order piecewise moments as dependent
variables, for better conservation of the information entropy. We have
developed one and twodimensional schemes and show that they provide quite
accurate solutions within reasonable usage of computational resources compared
to other existing schemes. The twodimensional scheme can accurately solve the
solid body rotation problem of a gaussian profile for more than hundred
rotation periods with little numerical diffusion. This is crucially important
for Vlasov simulations of magnetized plasmas. Applications of the one and
twodimensional schemes to electrostatic and electromagnetic Vlasov simulations
are presented with some benchmark tests.

$\alpha$YbAlB$_4$ is the locally isostructural polymorph of
$\beta$YbAlB$_4$, the first example of an Ybbased heavy fermion
superconductor which exhibits pronounced nonFermiliquid behavior above
$T_{\rm c}$. Interestingly, both $\alpha$YbAlB$_4$ and $\beta$YbAlB$_4$ have
strongly intermediate valence. Our single crystal study of the specific heat,
magnetization and resistivity has confirmed the Fermi liquid ground state of
$\alpha$YbAlB$_4$ ~in contrast with the quantum criticality observed in
$\beta$YbAlB$_4$. Both systems exhibit Kondo lattice behavior with the
characteristic temperature scale $T^* \sim$ 8 K in addition to a valence
fluctuation scale $\sim 200$ K. Below $T^*$, $\alpha$YbAlB$_4$ a heavy Fermi
liquid state with an electronic specific heat coefficient $\gamma\sim$ 130
mJ/mol K$^2$ and a large Wilson ratio more than 7, which indicates
ferromagnetic correlation between Yb moments. A large anisotropy in the
resistivity suggests that the hybridization between 4$f$ and conduction
electrons is much stronger in the $ab$plane than along the caxis. The
strongly anisotropic hybridization as well as the large Wilson ratio is the key
to understand the unusual Kondo lattice behavior and heavy fermion formation in
mixed valent compounds.

Fermi liquid theory, the standard theory of metals, has been challenged by a
number of observations of anomalous metallic behavior found in the vicinity of
a quantum phase transition. The breakdown of the Fermi liquid is accomplished
by finetuning the material to a quantum critical point using a control
parameter such as the magnetic field, pressure, or chemical composition. Our
high precision magnetization measurements of the ultrapure felectron based
superconductor {\beta}YbAlB4 demonstrate a scaling of its free energy
indicative of zerofield quantum criticality without tuning in a metal. The
breakdown of Fermiliquid behavior takes place in a mixedvalence state, in
sharp contrast with other known examples of quantum critical felectron systems
that are magnetic Kondo lattice systems with integral valence.