Recent exoplanet observations reported a large number of multiple-planet systems, in which some of the planets are in a chain of resonances. The fraction of resonant systems to non-resonant systems provides clues about their formation history. We investigated the orbital stability of planets in resonant chains by considering the long-term evolution of planetary mass and stellar mass and using orbital calculations. We found that while resonant chains were stable, they can be destabilized by a change of $\sim$10% in planetary mass. Such a mass evolution can occur by atmospheric escape due to photoevaporation. We also found that resonant chains can be broken by a stellar mass loss of $\lesssim1$%, which would be explained by stellar winds or coronal mass ejections. The long-term mass change of planets and stars plays an important role in the orbital evolutions of planetary systems including super-Earths.

Fully kinetic two-dimensional particle-in-cell simulations are used to study electron acceleration at high-Mach-number nonrelativistic perpendicular shocks. SNR shocks are mediated by the Weibel instability which is excited because of an interaction between shock-reflected and upstream ions. Nonlinear evolution of the Weibel instability leads to the formation of current sheets. At the turbulent shock ramp the current sheets decay through magnetic reconnection. The number of reconnection sites strongly depends on the ion-to-electron mass ratio and the Alfv\'enic Mach number of the simulated shock. Electron acceleration is observed at locations where magnetic reconnection operates. For the highest mass ratios almost all electrons are involved in magnetic reconnection, which makes the magnetic reconnection the dominant acceleration process for electrons at these shocks. We discuss the relevance of our results for 3D systems with realistic ion-to-electron mass ratio.

The first-order Fermi acceleration of electrons requires an injection of electrons into a mildly relativistic energy range. However, the mechanism of injection has remained a puzzle both in theory and observation. We present direct evidence for a novel stochastic shock drift acceleration theory for the injection obtained with Magnetospheric Multiscale (MMS) observations at Earth's bow shock. The theoretical model can explain electron acceleration to mildly relativistic energies at high-speed astrophysical shocks, which may provide a solution to the long-standing issue of electron injection.

Major components of chondrites are chondrules and matrix. Measurements of the volatile abundance in Semarkona chondrules suggest that chondrules formed in a dense clump that had a higher solid density than the gas density in the solar nebula. We investigate collisions between chondrules and matrix in the surface region of dense clumps using fluffy aggregate growth models. Our simulations show that the collisional growth of aggregates composed of chondrules and matrix takes place in the clumps well before they experience gravitational collapse. The internal structure of chondrite parent bodies (CPBs) can be thereby determined by aggregate growth. We find that the aggregate growth generates two scales within CPBs. The first scale is involved with the small scale distribution of chondrules and determined by the early growth stage, where chondrules accrete aggregates composed of matrix grains. This accretion can reproduce the thickness of the matrix layer around chondrules found in chondrites. The other scale is related to the large scale distribution of chondrules. Its properties (e.g., the abundance of chondrules and the overall size) depend on the gas motion within the clump, which is parameterized in this work. Our work thus suggests that the internal structure of CPBs may provide important clues about their formation conditions and mechanisms.

We explore electron pre-acceleration at high Mach-number nonrelativistic perpendicular shocks at, e.g., young supernova remnants, which are a prerequisite of further acceleration to very high energies via diffusive shock acceleration. Using fully kinetic particle-in-cell simulations of shocks and electron dynamics in them, we investigate the influence of shock-surfing acceleration at the shock foot on the nonthermal population of electrons downstream of the shock. The shock-surfing acceleration is followed by further energization at the shock ramp where the Weibel instability spawns a type of second-order Fermi acceleration. The combination of these two processes leads to the formation of a nonthermal electron population, but the importance of shock-surfing acceleration becomes smaller for larger ion-to-electron mass ratio in the simulation. We discuss the resulting electron spectra and the relevance of our results to the physics of systems with real ion-to-electron mass ratio and fully three-dimensional behavior.

Understanding chondrule formation provides invaluable clues about the origin of the solar system. Recent studies suggest that planetesimal collisions and the resulting impact melts are promising for forming chondrules. Given that the dynamics of planetesimals is a key in impact-based chondrule formation scenarios, we here perform direct $N$-body simulations to examine how the presence of Jupiter affects the properties of chondrule-forming collisions. Our results show that the absence/presence of Jupiter considerably changes the properties of high velocity collisions whose impact velocities are higher than 2.5 km s$^{-1}$; high velocity collisions occur due to impacts between protoplanets and planetesimals for the case without Jupiter; for the case with Jupiter, eccentricities of planetesimals are pumped up by the secular and resonant perturbations from Jupiter. We also categorize the resulting planetesimal collisions and find that most of high velocity collisions are classified as grazing ones for both cases. To examine the effect of Jupiter on chondrule formation directly, we adopt the impact jetting scenario and compute the resulting abundance of chondrules. Our results show that for the case without Jupiter, chondrule formation proceeds in the inside-out manner, following the growth of protoplanets. If Jupiter is present, the location and timing of chondrule formation are determined by Jupiter's eccentricity, which is treated as a free parameter in our simulations. Thus, the existence of Jupiter is the key parameter for specifying when and where chondrule formation occurs for impact-based scenarios.

We investigated electromagnetic precursor wave emission in relativistic shocks by using two-dimensional particle-in-cell simulations. We found that the wave amplitude is significantly enhanced by a positive feedback process associated with ion-electron coupling through the wakefields for high magnetization. The wakefields collapse during the nonlinear process of the parametric decay instability in the near-upstream region, where nonthermal electrons and ions are generated. The intense coherent emission and the particle acceleration may opperate in high-energy astrophysical objects.

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.

Electron injection at high Mach-number nonrelativistic perpendicular shocks is studied here for parameters that are applicable to young SNR shocks. Using high-resolution large-scale two-dimensional fully kinetic particle-in-cell (PIC) simulations and tracing individual particles we in detail analyze the shock surfing acceleration (SSA) of electrons at the leading edge of the shock foot. The central question is to what degree the process can be captured in 2D3V simulations. We find that the energy gain in SSA always arises from the electrostatic field of a Buneman wave. Electron energization is more efficient in the out-of-plane orientation of the large-scale magnetic field because both the phase speed and the amplitude of the waves are higher than for the in-plane scenario. Also, a larger number of electrons is trapped by the waves compared to the in-plane configuration. We conclude that significant modifications of the simulation parameters are needed to reach the same level of SSA efficiency as in simulations with out-of-plane magnetic field or 3D simulations.

We design a conservative finite difference scheme for ideal magnetohydrodynamic simulations that attains high-order accuracy, shock-capturing, and divergence-free condition of the magnetic field. The scheme interpolates pointwise physical variables from computational nodes to midpoints through a high-order nonlinear weighted average. The numerical flux is evaluated at the midpoint by a multi-state approximate Riemann solver for correct upwinding, and its spatial derivative is approximated by a high-order linear central difference to update the variables with designed order of accuracy and conservation. The magnetic and electric fields are defined at staggered grid points employed in the Constrained Transport (CT) method by Evans & Hawley (1988). We propose a new CT variant, in which the staggered electric field is evaluated so as to be consistent with the base one-dimensional Riemann solver and the staggered magnetic field is updated to be divergence-free as designed high-order finite difference representation. We demonstrate various benchmark tests to measure the performance of the present scheme. We discuss the effect of the choice of interpolation methods, Riemann solvers, and the treatment for the divergence-free condition on the quality of numerical solutions in detail.

In this paper, we report our recent findings on the relativistic Weibel instability and its nonlinear saturation by performing numerical simulations of collisionless plasmas. Analysis of the obtained numerical results revealed that the nonlinear phase of the Weibel instability can be described by characteristic phases based on the Weibel filaments' current density in terms of particle and Alfv\'en limit currents. We also analyzed the relativistic kink instability based on the energy principle in the magnetohydrodynamic (MHD) regime, and found that the Weibel filaments do not suffer from the kink-type instability in the MHD regime up to 1000 $\omega_{p,i}^{-1}$. This finding allowed a magnetic field to be sustained by relativistic Weibel instability that was stable enough to be a seed for MHD dynamos.

We discuss a mechanism for intense electromagnetic wave emission at an astrophysical relativistic shock in a magnetized collisionless plasma. At the magnetized shock, the particle reflection by a compressed magnetic field of the shock produces a ring-like distribution in momentum, which gives rise to plasma instabilities. Intense and coherent high-frequency electromagnetic waves will be emitted if the synchrotron maser instability (SMI) is excited, whereas non-propagating magnetic fluctuations will be generated when the Weibel instability (WI) is the dominant mode. The problem is of great astrophysical interest because if intense radiation is emitted, the interaction with the upstream medium induces a large-amplitude electrostatic field (or Wakefield), which may play a role for the acceleration of ultra-high-energy cosmic rays. We review our recent effort to measure the efficiency of the electromagnetic wave emission using fully self-consistent, two-dimensional (2D) particle-in-cell (PIC) simulations for pair plasmas. We found that the emission efficiency in 2D was systematically lower than one dimensional (1D) PIC simulation results. However, the power remains finite even when the WI is active to generate large-amplitude magnetic fluctuations. Astrophysical implications of the present results are briefly discussed.

The time evolution and saturation of the Weibel instability at the ion Alfv\'en current are presented by ab initio particle-in-cell simulations. We found that the ion Weibel current in 3D could evolve into the Alfv\'en current where the magnetic field energy is sustained at 1.5\% of the initial beam kinetic energy. The current filaments are no longer isolated at saturation, but rather connected to each other to form a network structure. Electrons are continuously heated during the coalescence of the filaments, which is crucial for obtaining sustained magnetic fields with much stronger levels than with 2D simulations. The results highlight again the importance of the Weibel instability in generating magnetic fields in laboratory, astrophysical, and cosmological situations.

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.

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.

Recent observations have revealed the eccentricity and inclination distributions of close-in super-Earths. These distributions have the potential to constrain their formation processes. In the in-situ formation scenario, the eccentricities and inclinations of planets are determined by gravitational scattering and collisions between protoplanets on the giant impact stage. We investigate the effect of the initial eccentricities and inclinations of protoplanets on the formation of close-in super-Earths. We perform $N$-body simulations of protoplanets in gas-free disks, changing the initial eccentricities and inclinations systematically. We find that while the eccentricities of protoplanets are well relaxed through their evolution, the inclinations are not. When the initial inclinations are small, they are not generally pumped up since scattering is less effective and collisions occur immediately after orbital crossing. On the other hand, when the initial inclinations are large, they tend to be kept large since collisional damping is less effective. Not only the resultant inclinations of planets, but also their number, eccentricities, angular momentum deficit, and orbital separations are affected by the initial inclinations of protoplanets.

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 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.

Chondrules are primitive materials in the Solar System. They are formed in the first about 3 Myr of the Solar System's history. This timescale is longer than that of Mars formation, and it is conceivable that protoplanets, planetesimals and chondrules might have existed simultaneously in the solar nebula. Due to protoplanets perturbation on the planetesimal dynamics and chondrule accretion on them, all the formed chondrules are unlikely to be accreted by planetesimals. We investigate the amount of chondrules accreted by planetesimals in such a condition. We assume that a protoplanet is in oligarchic growth, and we perform analytical calculations of chondrule accretion both by a protoplanet and by planetesimals. Through the oligarchic growth stage, planetesimals accrete about half of the formed chondrules. The smallest planetesimals get the largest amount of the chondrules, compared with the amount accreted by more massive planetesimals. We perform a parameter study and find that this fraction is not largely changed for a wide range of parameter sets.

Chondritic meteorites contain unique spherical materials named chondrules: sub-mm sized silicate grains once melted in a high temperature condition in the solar nebula. We numerically explore one of chondrule forming processes, planetesimal collisions. Previous studies found that impact jetting via protoplanet-planetesimal collisions make chondrules with an amount of 1 % of impactors' mass, when impact velocity exceeds 2.5 km s$^{-1}$. Based on the mineralogical data of chondrules, undifferentiated planetesimals would be more suitable for chondrule-forming collisions than potentially differentiated protoplanets. We examine planetesimal-planetesimal collisions using a shock physics code and find two things: one is that planetesimal-planetesimal collisions produce the nearly same amount of chondrules as protoplanet-planetesimal collisions ($\sim$ 1 %). The other is that the amount of produced chondrules becomes larger as the impact velocity increases when two planetesimals collide with each other. We also find that progenitors of chondrules can be ejected from deeper regions of large targets (planetesimals or protoplanets) than small impactors (planetesimals). The composition of targets is therefore important to fully account for the mineralogical data of currently sampled chondrules.

Chondritic meteorites provide valuable opportunities to investigate the origins of the solar system. We explore impact jetting as a mechanism of chondrule formation and subsequent pebble accretion as a mechanism of accreting chondrules onto parent bodies of chondrites, and investigate how these two processes can account for the currently available meteoritic data. We find that when the solar nebula is $\le 5$ times more massive than the minimum-mass solar nebula at $a \simeq 2-3$ AU and parent bodies of chondrites are $\le 10^{24}$ g ($\le$ 500 km in radius) in the solar nebula, impact jetting and subsequent pebble accretion can reproduce a number of properties of the meteoritic data. The properties include the present asteroid belt mass, the formation timescale of chondrules, and the magnetic field strength of the nebula derived from chondrules in Semarkona. Since this scenario requires a first generation of planetesimals that trigger impact jetting and serve as parent bodies to accrete chondrules, the upper limit of parent bodies' masses leads to the following implications: primordial asteroids that were originally $\ge 10^{24}$ g in mass were unlikely to contain chondrules, while less massive primordial asteroids likely had a chondrule-rich surface layer. The scenario developed from impact jetting and pebble accretion can therefore provide new insights into the origins of the solar system.

Chondrules are one of the most primitive elements that can serve as a fundamental clue as to the origin of our Solar system. We investigate a formation scenario of chondrules that involves planetesimal collisions and the resultant impact jetting. Planetesimal collisions are the main agent to regulate planetary accretion that corresponds to the formation of terrestrial planets and cores of gas giants. The key component of this scenario is that ejected materials can melt when the impact velocity between colliding planetesimals exceeds about 2.5 km s$^{-1}$. The previous simulations show that the process is efficient enough to reproduce the primordial abundance of chondrules. We examine this scenario carefully by performing semi-analytical calculations that are developed based on the results of direct $N$-body simulations. As found by the previous work, we confirm that planetesimal collisions that occur during planetary accretion can play an important role in forming chondrules. This arises because protoplanet-planetesimal collisions can achieve the impact velocity of about 2.5 km s$^{-1}$ or higher, as protoplanets approach the isolation mass ($M_{p,iso}$). Assuming that the ejected mass is a fraction ($F_{ch}$) of colliding planetesimals' mass, we show that the resultant abundance of chondrules is formulated well by $F_{ch}M_{p,iso}$, as long as the formation of protoplanets is completed within a given disk lifetime. We perform a parameter study and examine how the abundance of chondrules and their formation timing change. We find that the impact jetting scenario generally works reasonably well for a certain range of parameters, while more dedicated work would be needed to include other physical processes that are neglected in this work and to examine their effects on chondrule formation.

Most of super-Earths detected by the radial velocity (RV) method have significantly smaller eccentricities than the eccentricities corresponding to velocity dispersion equal to their surface escape velocity ("escape eccentricities"). If orbital instability followed by giant impacts among protoplanets that have migrated from outer region is considered, it is usually considered that eccentricities of the merged bodies become comparable to those of orbital crossing bodies, which are excited up to their escape eccentricities by close scattering. However, the eccentricity evolution in the {\it in situ} accretion model has not been studied in detail. Here, we investigate the eccentricity evolution through {\it N}-body simulations. We have found that the merged planets tend to have much smaller eccentricities than the escape eccentricities due to very efficient collision damping. If the protoplanet orbits are initially well separated and their eccentricities are securely increased, an inner protoplanet collides at its apocenter with an outer protoplanet at its pericenter. The eccentricity of the merged body is the smallest for such configuration. Orbital inclinations are also damped by this mechanism and planets tend to share a same orbital plane, which is consistent with {\it Kepler} data. Such efficient collision damping is not found when we start calculations from densely packed orbits of the protoplanets. If the protoplanets are initially in the mean-motion resonances, which corresponds to well separated orbits, the {\it in situ} accretion model well reproduces the features of eccentricities and inclinations of multiple super-Earths/Earth systems discovered by RV and {\it Kepler} surveys.

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.

Electron acceleration associated with various plasma kinetic instabilities in a nonrelativistic, very-high-Alfv\'en Mach-number ($M_A \sim 45$) shock is revealed by means of a two-dimensional fully kinetic PIC simulation. Electromagnetic (ion Weibel) and electrostatic (ion-acoustic and Buneman) instabilities are strongly activated at the same time in different regions of the two-dimensional shock structure. Relativistic electrons are quickly produced predominantly by the shock surfing mechanism with the Buneman instability at the leading edge of the foot. The energy spectrum has a high-energy tail exceeding the upstream ion kinetic energy accompanying the main thermal population. This gives a favorable condition for the ion acoustic instability at the shock front, which in turn results in additional energization. The large-amplitude ion Weibel instability generates current sheets in the foot, implying another dissipation mechanism via magnetic reconnection in a three-dimensional shock structure in the very-high-$M_A$ regime.

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.

Many extrasolar planetary systems containing multiple super-Earths have been discovered. N-body simulations taking into account standard type-I planetary migration suggest that protoplanets are captured into mean-motion resonant orbits near the inner disk edge at which the migration is halted. Previous N-body simulations suggested that orbital stability of the resonant systems depends on number of the captured planets. In the unstable case, through close scattering and merging between planets, non-resonant multiple systems are finally formed. In this paper, we investigate the critical number of the resonantly trapped planets beyond which orbital instability occurs after disk gas depletion. We find that when the total number of planets ($N$) is larger than the critical number ($N_{\rm crit}$), crossing time that is a timescale of initiation of the orbital instability is similar to non-resonant cases, while the orbital instability never occurs within the orbital calculation time ($10^8$ Kepler time) for $N\leq N_{\rm crit}$. Thus, the transition of crossing time across the critical number is drastic. When all the planets are trapped in 7:6 resonance of adjacent pairs, $N_{\rm crit} = 4$. We examine the dependence of the critical number of 4:3, 6:5 and 8:7 resonance by changing the orbital separation in mutual Hill radii and planetary mass. The critical number increases with increasing the orbital separation in mutual Hill radii with fixed planetary mass and increases with increasing planetary mass with fixed the orbital separation in mutual Hill radii. We also calculate the case of a system which is not composed of the same resonance. The sharp transition of the stability can be responsible for the diversity of multiple super-Earths (non-resonant or resonant), that is being revealed by $Kepler$ mission.

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.

Using fully relativistic GEANT4 simulation tool kit, the transport of energetic electrons generated in solar flares was Monte-Carlo simulated, and resultant bremsstrahlung gamma-ray spectra were calculated. The solar atmosphere was approximated by 10 vertically-stacked zones. The simulation took into account two important physical processes,that the bremsstrahlung photons emitted by precipitating relativistic electrons are strongly forward beamed toward the photosphere, and that the majority of these gamma-rays must be Compton back-scattered by the solar atmosphere in order to reach the observer. Then, the Compton degradation was found to make the observable gamma-ray spectra much softer than is predicted by simple analytic calculations. The gamma-ray signals were found to be enhanced by several conditions, including a broad pitch-angle distribution of the electrons, a near-limb flare longitude, and a significant tilt in the magnetic field lines if the flare longitude is rather small. These results successfully explain several important flare properties observed in the hard X-ray to gamma-ray range, including in particular those obtained with Yohkoh. A comparison of the Yohkoh spectrum from a GOES X3.7 class limb flare on 1998 November 22, with a simulation assuming a broad electron pitch-angle distribution, suggests that gamma-rays from this particular solar flare were a mixture of direct bremsstrahlung photons and their Comptonization.

The Hard X-ray Detector (HXD) on board Suzaku covers a wide energy range from 10 keV to 600 keV by combination of silicon PIN diodes and GSO scintillators. The HXD is designed to achieve an extremely low in-orbit back ground based on a combination of new techniques, including the concept of well-type active shield counter. With an effective area of 142 cm^2 at 20 keV and 273 cm2 at 150 keV, the background level at the sea level reached ~1x10^{-5} cts s^{-1} cm^{-2} keV^{-1} at 30 keV for the PI N diodes, and ~2x10^{-5} cts s^{-1} cm^{-2} keV^{-1} at 100 keV, and ~7x10^{-6} cts s^{-1} cm^{-2} keV^{-1} at 200 keV for the phoswich counter. Tight active shielding of the HXD results in a large array of guard counters surrounding the main detector parts. These anti-coincidence counters, made of ~4 cm thick BGO crystals, have a large effective area for sub-MeV to MeV gamma-rays. They work as an excellent gamma-ray burst monitor with limited angular resolution (~5 degree). The on-board signal-processing system and the data transmitted to the ground are also described.

Gamma-ray emitting solar flares observed with Yohkoh were analyzed from a statistical viewpoint. The four-band hard X-ray (15--95 keV) photometric data, taken with the Hard X-ray Telescope onboard Yohkoh, were utilized in combination with the spectro-photometric gamma-ray (0.2--30 MeV) data obtained with the Gamma-Ray Spectrometer. The GOES class was also incorporated. Out of 2788 X-ray flares observed from 1991 October to 2001 December, 178 events with strong hard X-ray emission were selected. Among them, 40 flares were further found to show significant gamma-ray emission. A fractal dimension analysis and multi-band color--color plots of the 40 flares suggest that their soft X-ray to MeV gamma-ray spectral energy distributions involve at least four independent parameters. These are: (1) the overall flare size; (2) the relative intensities of the thermal vs. non-thermal signals; (3) the gamma-ray to hard X-ray intensity ratio; and (4) the hard X-ray spectral slope. These results are examined for possible selection effects. Also, the meanings of the third parameter are briefly considered.

We observed the nearby galaxy M~51 (NGC 5194) with BeppoSAX. The X-ray properties of the nucleus below 10 keV are almost the same as the ASCA results regarding the hard component and the neutral Fe K$\alpha$ line, but the intensity is about half of the ASCA 1993 data. Beyond this, in the BeppoSAX PDS data, we detected a bright hard X-ray emission component which dominates above 10 keV. The 10 -- 100 keV flux and luminosity of this component are respectively $2\times10^{-11}$ erg s$^{-1}$ cm$^{-2}$ and $2\times10^{41}$ erg s$^{-1}$. These are about 10 times higher than the extrapolation from the soft X-ray band, and similar to the flux observed with Ginga, which found a bright power law component in 2 -- 20 keV band. Considering other wavelength properties and the X-ray luminosity, together with strong neutral Fe K line, the hard X-ray emission most likely arises from a low luminosity active nucleus, which is obscured with a column density of $\sim10^{24}$cm$^{-2}$. This suggests that hidden low luminosity AGNs may well be present in other nearby galaxies. We interpret the discrepancy between Ginga and other X-ray satellites to be due to a large variability of absorption column density toward the line of sight over several years, suggesting that the Compton thick absorption material may be present on a spatial scale of a parsec. Apart from the nucleus, several ultra-luminous off-nuclear X-ray sources detected in M~51 exhibit long-term time variability, suggesting the state transition similar to that observed in Galactic black hole candidates.

NGC 1316 (Fornax A) is a radio galaxy with prototypical double lobes, where the magnetic field intensity is accurately measured via the inverse-Compton technique. The radio-emitting electrons in the lobes are inferred to have a synchrotron life time of 0.1 Gyr. Considering the lobe energetics, we estimate the past nuclear X-ray luminosity of NGC 1316 to be at least 4 times 10^{34} W (4 times 10^{41} erg s^{-1}). Thus, the nucleus was rather active at least 0.1 Gyr ago. In contrast, we confirmed with ASCA and ROSAT that the nucleus of NGC 1316 is very faint in X-rays at present, with the 2--10 keV luminosity of any AGN-like hard component being < 2 times 10^{33} W (2 times 10^{40} erg s^{-1}) even assuming a nuclear obscuration up to 10^{28} m^{-2} (10^{24} cm^{-2}). This is at least an order of magnitude lower than the estimated past activity, indicating that the nucleus is presently very inactive. From these two results, we conclude that the nucleus of NGC 1316 has become dormant during the last 0.1 Gyr. This suggests the possible abundance of ``dormant'' quasars in nearby galaxies.