• 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 fifth-order-monotonicity-preserving 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 two-dimensional problems showed the advantages of using the high-order scheme by comparing with results from a standard second-order TVD MUSCL scheme. The present code enabled us to explore long-term evolution of a three-dimensional accretion disk around a black hole, in which compressible MHD turbulence caused continuous mass accretion via nonlinear growth of the magneto-rotational 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 non-magnetic correlated metal. In particular, we measured the magnetoresistance under fields up to 31.4 T in the cubic Pr-based heavy fermion superconductor PrV$_2$Al$_{20}$ with a non-magnetic $\Gamma _3$ doublet ground state, exhibiting antiferro-quadrupole ordering below 0.7 K. For the [100] direction, we find that the high-field 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 field-induced antiferro-quadrupole (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 two-dimensional particle-in-cell simulations and focused on two magnetic field orientations: the magnetic field to be in the simulation plane (i.e., in-plane configuration) and perpendicular to the simulation plane (i.e., out-of-plane configuration). Our simulations in the in-plane 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 large-amplitude 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 out-of-plane configuration investigated in the previous study. We confirmed that our previous conclusion holds regardless of upstream magnetic field orientations with respect to the two-dimensional 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 electron-hole 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 high-Mach-number quasi-perpendicular shock is presented by means of ab initio particle-in-cell simulations in three dimensions. We found that coherent electrostatic Buneman waves and ion-Weibel magnetic turbulence coexist in a strong-shock structure whereby particles gain energy during shock-surfing and subsequent stochastic drift accelerations. Energetic electrons that initially experienced the surfing acceleration undergo pitch-angle diffusion by interacting with magnetic turbulence and continuous acceleration during confinement in the shock transition region. The ion-Weibel turbulence is the key to the efficient nonthermal electron acceleration.
  • We investigated the efficiency of coherent upstream large-amplitude electromagnetic wave emission via synchrotron maser instability at relativistic magnetized shocks by using two-dimensional particle-in-cell simulations. We considered the purely perpendicular shock in an electron-positron 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 one-dimensional 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 non-thermal electrons in a relativistic magnetized ion-electron 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 power-law 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 time-delayed radio brightenings of Sgr A* triggered by the pericenter passage of the G2 cloud is studied by carrying out global three-dimensional 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 re-organized after the passage of G2. We have found that the magnetic energy in the accretion flow increases by a factor 3-4 in 5-10 years after the pericenter passage of G2 by a dynamo mechanism driven by the magneto-rotational instability. Since this B-field 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 time-delay of the radio brightening enables us to determine the rotation axis of the preexisting disk.
  • We report the results of 27Al-NMR 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 low-temperature phase, indicating the breaking of the three fold rotation symmetry due to the field-induced magnetic dipole perpendicular to the field. This provides microscopic evidence for a ferro order of the O20 quadrupole. Although a first-order 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 well-characterized sample of "stuffed" (Pr-rich) Pr$_{2+x}$Ir$_{2-x}$O$_{7-\delta}$. Elastic neutron scattering shows the onset of long-range spin-ice "2-in/2-out" 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 long-range spin-ice ordering of 1.7$\mu_B$ moments (120--270 mT, depending on muon site). This shortfall is due in part to splitting of the non-Kramers crystal-field ground-state doublets of near-neighbor 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 Pr-moment 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 particle-in-cell simulations. Using the implicit finite-difference time-domain method to solve Maxwell equations, we found that nonphysical instability was greatly inhibited with a Courant-Friedrichs-Lewy (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.5-0.7 were obtained with explicit field solvers. In addition, we found employing higher-order shape functions and an optimal implicitness factor further suppressed long-wavelength modes of the instability. The findings allowed the examination of the long-term 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, gamma-ray 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 high-quality 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 high-precision 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 zero-field quantum criticality (QC) without tuning. We examined other possible scaling relationship such as $T/(B-B_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 Fermi-liquid component of the magnetization after subtracting the QC component estimated based on the $T/B$ scaling. The temperature dependence of the heavy Fermi-liquid component is found very similar to the magnetization of the polymorph $\alpha$-YbAlB$_4$. In addition, the heavy Fermi-liquid 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 zero-field 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 high-quality single crystal of PrV$_2$Al$_{20}$ exhibits superconductivity at $T_{\rm c}=$ 50 mK in the antiferroquadrupole-ordered 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 high-quality 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 Fermi-liquid breakdown in an intermediate valence system. We employ high precision magnetization measurements of the valence fluctuating superconductor beta-YbAlB4 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 beta-YbAlB4 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 multi-moment advection scheme (Minoshima et al., 2011, J. Comput. Phys.) to the three-dimensional 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 two-dimensional scheme, the three-dimensional 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.1-10 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 Kelvin-Helmholtz 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 growth-rates 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 solar-system.
  • Electron accelerations at high Mach number collision-less shocks are investigated by means of two-dimensional electromagnetic Particle-in-Cell simulations with various Alfven Mach numbers, ion-to-electron 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 super-high 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 super-high Mach number shock results in a coherent electrostatic potential structure. While multi-dimensionality 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 extremely-high 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 two-dimensional schemes and show that they provide quite accurate solutions compared to other existing schemes with the same memory usage. The two-dimensional 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 two-dimensional schemes and show that they provide quite accurate solutions within reasonable usage of computational resources compared to other existing schemes. The two-dimensional 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 two-dimensional 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 Yb-based heavy fermion superconductor which exhibits pronounced non-Fermi-liquid 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 c-axis. 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 fine-tuning 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 f-electron based superconductor {\beta}-YbAlB4 demonstrate a scaling of its free energy indicative of zero-field quantum criticality without tuning in a metal. The breakdown of Fermi-liquid behavior takes place in a mixed-valence state, in sharp contrast with other known examples of quantum critical f-electron systems that are magnetic Kondo lattice systems with integral valence.