Gamma-ray observations have revealed strong variability in blazar luminosities in the gamma-ray band over time scales as short as minutes. We show, for the first time, that the correlation of the spectrum with intensity is consistent with the behavior of the luminosity variation of blazar SEDs along a blazar sequence for low synchrotron peak blazars. We show that the observational signatures of variability with flux are consistent with wakefield acceleration of electrons initiated by instabilities in the blazar accretion disk. This mechanism reproduces the observed time variations as short as 100 seconds. The wakefield mechanism also predicts a reduction of the electron spectral index with increased gamma-ray luminosity, which could be detected in higher energy observations well above the inverse Compton peak.

Compilation of papers presented by the JEM-EUSO Collaboration at the 36th International Cosmic Ray Conference (ICRC), held July 24 through August 1, 2019 in Madison, Wisconsin.

Mini-EUSO will observe the Earth in the UV range (300 - 400 nm) offering the opportunity to study a variety of atmospheric events such as Transient Luminous Events (TLEs), meteors and marine bioluminescence. Furthermore it aims to search for Ultra High Energy Cosmic Rays (UHECR) above $10^{21}$ eV and Strange Quark Matter (SQM). The detector is expected to be launched to the International Space Station in August 2019 and look at the Earth in nadir mode from the UV-transparent window of the Zvezda module of the International Space Station. The instrument comprises a compact telescope with a large field of view ($44^{\circ}$), based on an optical system employing two Fresnel lenses for light collection. The light is focused onto an array of 36 multi-anode photomultiplier tubes (MAPMT), for a total of 2304 pixels and the resulting signal is converted into digital, processed and stored via the electronics subsystems on-board. In addition to the main detector, Mini-EUSO contains two ancillary cameras for complementary measurements in the near infrared (1500 - 1600 nm) and visible (400 - 780 nm) range and also a 8x8 SiPM imaging array.

Space Debris (SD) consist of non-operational artificial objects orbiting around the Earth, which could possibly damage space vehicles, such as the International Space Station (ISS) or other manned spacecrafts. The vast majority of such objects are cm-sized, not catalogued and usually the tracking data are not precise enough. Here we present the feasibility study of SD detection and tracking with techniques usually employed in cosmic-ray physics. For this purpose, we have evaluated the possibility of using Mini-EUSO, a space-borne fluorescence telescope to be deployed on the ISS, to track SD illuminated by the Sun. By means of ESAF (EUSO Simulation and analysis Framework) simulation and by developing the trigger algorithms, we estimated the minimum size and maximum distances of detectable SD. We then studied the number of possible SD detections using an ESA software called MASTER (Meteoroid and SD Terrestrial Environment Reference). With the Mini-EUSO Engineering Model (Mini-EUSO EM), we performed some measurements to estimate the reflectance of the most common SD materials and to demonstrate the ability of Mini-EUSO to detect SD events. We also performed some tests in open-sky conditions, identifying and tracking fast-moving objects. In particular, the detection of a rocket body allowed us to confirm the simulation outcomes predictions and the expected performance of the detector.

Mini-EUSO is a UV telescope that will look downwards to the Earth's atmosphere onboard the International Space Station. With the design of the ultra-high energy cosmic ray fluorescence detectors belonging to the JEM-EUSO program, it will make the first UV map of the Earth by observing atmospheric phenomena such as transient luminous events, sprites and lightning, as well as meteors and bioluminescence from earth. Diffused light from laser shots from the ground, which mimic the fluorescence light emitted by Nitrogen molecules when extensive air showers pass through the atmosphere, can be used to verify the capability of this kind of detector to observe ultra-high energy cosmic rays. To validate the electronics and the trigger algorithms developed for Mini-EUSO, a scaled down version of the telescope with 1:9 of the original focal surface and a lens of 2.5 cm diameter has been built. Tests of the Mini-EUSO engineering model have been made in laboratory and in open sky condition. In this paper, we report results of observations of the night sky, which include the detection of stars, meteors, a planet and a rocket body reflecting the sunlight. Interesting results of the observation of city lights are also reported.

EUSO-TA is a ground-based florescence detector built to validate the design of an ultra-high energy cosmic ray fluorescence detector to be operated in space. EUSO-TA detected the first air shower events with the technology developed within the JEM-EUSO program. It operates at the Telescope Array (TA) site in Utah, USA. With the external trigger provided by the Black Rock Mesa fluorescence detectors of Telescope Array (TA-FDs), EUSO-TA observed nine ultra-high energy cosmic ray events and several laser events from the Central Laser Facility of Telescope Array and portable lasers like the JEM-EUSO Global Light System prototype. The reconstruction parameters of the cosmic ray events which crossed the EUSO-TA field of view (both detected and not detected by EUSO-TA), were provided by the Telescope Array Collaboration. As the TA-FDs have a wider field of view than EUSO-TA ($\sim$30 times larger), they allow the cosmic ray energy reconstruction based on the observation of most of the extensive air-shower profiles, including the shower maximum, while EUSO-TA only observes a portion of the showers, usually far from the maximum. For this reason, the energy of the cosmic rays corresponding to the EUSO-TA signals appear lower than the actual ones. In this contribution, the analysis of the cosmic-ray events detected with EUSO-TA is discussed.

The TurLab facility is a laboratory, equipped with a 5 m diameter and 1 m depth rotating tank, located in the Physics Department of the University of Turin. Originally, it was mainly built to study systems of different scales where rotation plays a key role in the fluid behavior such as in atmospheric and oceanic flows. In the past few years the TurLab facility has been used to perform experiments related to the observation of Extreme Energy Cosmic Rays (EECRs) from space using the fluorescence technique. For example, in the case of the JEM-EUSO mission, where the diffuse night brightness and artificial light sources can vary significantly in time and space inside the Field of View of the telescope. The Focal Surface of Mini-EUSO Engineering Model (Mini-EUSO EM) with the level 1 (L1) and 2 (L2) trigger logics implemented in the Photo-Detector Module (PDM) has been tested at TurLab. Tests related to the possibility of using an EUSO-like detector for other type of applications such as Space Debris (SD) monitoring and imaging detector have also been pursued. The tests and results obtained within the EUSO@TurLab Project on these different topics are presented.

Mini-EUSO is a space experiment selected to be installed inside the International Space Station. It has a compact telescope with a large field of view ($44 $\times$ 44$ sq. deg.) focusing light on an array of photo-multipliers tubes in order to observe UV emission coming from Earth's atmosphere. Observations will be complemented with data recorded by some ancillary detectors. In particular, the Mini-EUSO Additional Data Acquisition System (ADS) is composed by two cameras, which will allow us to obtain data in the near infrared, and in the visible range. These will be used to monitor the observation conditions, and to acquire useful information on several scientific topics to be studied with the main instrument, such as the physics of atmosphere, meteors, and strange quark matter. Here we present the ADS control software developed to stream cameras together with the UV main instrument, in order to grab images in an automated and independent way, and we also describe the calibration activities performed on these two ancillary cameras before flight.

We consider that electromagnetic pulses produced in the jets of this innermost part of the accretion disk accelerate charged particles (protons, ions, electrons) to very high energies including energies above $10^{20}$ eV for the case of protons and nucleus and $10^{12-15}$ eV for electrons by electromagnetic wave-particle interaction. The episodic eruptive accretion in the disk by the magneto-rotational instability gives rise to the strong Alfvenic pulses, which acts as the driver of the collective accelerating pondermotive force. This pondermotive force drives the wakes. The accelerated hadrons (protons and nuclei) are released to the intergalactic space to be ultra-high energy cosmic rays. The high-energy electrons, on the other hand, emit photons in the collisions of electromagnetic perturbances to produce various non-thermal emissions (radio, IR, visible, UV, and gamma-rays) of active galactic nuclei. Applying the theory to M82 X-1, we find that it can explain the northern hot spot of ultra high energy cosmic rays above $6\times 10^{19}$ eV. We also discuss astrophysical implications for other nearby active galactic nuclei, neutron star mergers, and high energy neutrinos.

The new technique of measuring frequency by optical lattice clocks now approaches to the relative precision of $(\Delta f/f)=O(10^{-18})$. We propose to place such precise clocks in space and to use Doppler tracking method for detecting low-frequency gravitational wave below 1 Hz. Our idea is to locate three spacecrafts at one A.U. distance (say at L1, L4 & L5 of the Sun-Earth orbit), and apply the Doppler tracking method by communicating "the time" each other. Applying the current available technologies, we obtain the sensitivity for gravitational wave with three or four-order improvement ($h_{\rm n}\sim 10^{-17}$ or $10^{-18}$ level in $10^{-5}$Hz -- $1$ Hz) than that of Cassini spacecraft in 2001. This sensitivity enables us to observe black-hole mergers of their mass greater than $10^5 M_\odot$ in the cosmological scale. Based on the hierarchical growth model of black-holes in galaxies, we estimate the event rate of detection will be 20-50 a year. We nickname "INO" (Interplanetary Network of Optical Lattice Clocks) for this system, named after Tadataka Ino (1745--1818), a Japanese astronomer, cartographer, and geodesist.

EUSO-Balloon is a pathfinder mission for the Extreme Universe Space Observatory onboard the Japanese Experiment Module (JEM-EUSO). It was launched on the moonless night of the 25$^{th}$ of August 2014 from Timmins, Canada. The flight ended successfully after maintaining the target altitude of 38 km for five hours. One part of the mission was a 2.5 hour underflight using a helicopter equipped with three UV light sources (LED, xenon flasher and laser) to perform an inflight calibration and examine the detectors capability to measure tracks moving at the speed of light. We describe the helicopter laser system and details of the underflight as well as how the laser tracks were recorded and found in the data. These are the first recorded laser tracks measured from a fluorescence detector looking down on the atmosphere. Finally, we present a first reconstruction of the direction of the laser tracks relative to the detector.

The episodic dynamics of the magnetic eruption of a spinning black hole (BH) accretion disks and its associated intense shapeup of their jets is studied via three-dimensional general-relativistic magnetohydrodynamics (GRMHD). The embedded magnetic fields in the disk get amplified by the magnetorotational instability (MRI) so large as to cause an eruption of magnetic field (recconection) and large chunks of matter episodically accrete toward the roots of the jets upon such an event. We also find that the eruption events produce intensive Alfv\'en pulses, which propagate through the jets. After the eruption, the disk backs to the weakly magnetic states. Such disk activities cause short time variabilities in mass accretion rate at the event horizon as well as electromagnetic luminosity inside the jet. Since the dimensionless strength parameter $a_0=eE/m_e \omega c$ of these Alfv\'en wave pulses is extremely high for a substantial fraction of Eddington accretion rate accretion flow onto a supermassive black hole, the Alfv\'en shocks turn into ultrarelativistic $(a_0\gg 1)$ bow wake acceleration, manifesting into the ultra-high energy cosmic rays and electrons which finally emit gamma-rays. Since our GRMHD model has universality in its spatial and temporal scales, it is applicable to a wide range of astrophysical objects ranging from those of AGN (which is the primary target of this research), to micro-quasars. Properties such as time variabilities of blazar gamma-ray flares and spectrum observed by {\it Fermi} Gamma-ray Observatory are well explained by linear acceleration of electrons by the bow wake.

Based on a dynamical formation model of a supermassive black hole (SMBH), we estimate the expected observational profile of gravitational wave at ground-based detectors, such as KAGRA or advanced LIGO/VIRGO. Noting that the second generation of detectors have enough sensitivity from 10 Hz and up (especially with KAGRA owing to its location at less seismic noise), we are able to detect the ring-down gravitational wave of a BH with the mass $M < 2\times 10^3 M_\odot $. This enables us to check the sequence of BH mergers to SMBHs via intermediate-mass BHs. We estimate the number density of galaxies from the halo formation model and estimate the number of BH mergers from the giant molecular cloud model assuming hierarchical growth of merged cores. At the designed KAGRA (and/or advanced LIGO/VIRGO), we find that the BH merger of its total mass $M\sim 60M_\odot$ is at the peak of the expected mass distribution. With its signal-to-noise ratio $\rho=10 (30)$, we estimate the event rate $R \sim 200 (20)$ per year in the most optimistic case, and we also find that BH mergers in the range $M < 150 M_\odot$ are $R>1$ per year for $\rho=10$. Thus, if we observe a BH with more than $100 M_\odot$ in future gravitational-wave observations, our model naturally explains its source.

Population synthesis studies into planet formation have suggested that distributions consistent with observations can only be reproduced if the actual Type I migration timescale is at least an order of magnitude longer than that deduced from linear theories. Although past studies considered the effect of the Type I migration of protoplanetary embryos, in most cases they used a conventional formula based on static torques in isothermal disks, and employed a reduction factor to account for uncertainty in the mechanism details. However, in addition to static torques, a migrating planet experiences dynamic torques that are proportional to the migration rate. These dynamic torques can impact on planet migration and predicted planetary populations. In this study, we derived a new torque formula for Type I migration by taking into account dynamic corrections. This formula was used to perform population synthesis simulations with and without the effect of dynamic torques. In many cases, inward migration was slowed significantly by the dynamic effects. For the static torque case, gas giant formation was effectively suppressed by Type I migration; however, when dynamic effects were considered, a substantial fraction of cores survived and grew into gas giants.

We have identified iridium in an ~5 m-thick section of pelagic sediment cored in the deep sea floor at Site 886C, in addition to a distinct spike in iridium at the K-Pg boundary related to the Chicxulub asteroid impact. We distinguish the contribution of the extraterrestrial matter in the sediments from those of the terrestrial matter through a Co-Ir diagram, calling it the "extraterrestrial index" fEX. This new index reveals a broad iridium anomaly around the Chicxulub spike. Any mixtures of materials on the surface of the Earth cannot explain the broad iridium component. On the other hand, we find that an encounter of the solar system with a giant molecular cloud can aptly explain the component, especially if the molecular cloud has a size of ~100 pc and the central density of over 2000 protons/cm^3. Kataoka et al. (2013, 2014) pointed that an encounter with a dark cloud would drive an environmental catastrophe leading to mass extinction. Solid particles from the hypothesized dark cloud would combine with the global environment of Earth, remaining in the stratosphere for at least several months or years. With a sunshield effect estimated to be as large as -9.3 W m^-2, the dark cloud would have caused global climate cooling in the last 8 Myr of the Cretaceous period, consistent with the variations of stable isotope ratios in oxygen (Barrera and Huber, 1990; Li and Keller, 1998; Barrera and Savin, 1999; Li and Keller, 1999) and strontium (Barrera and Huber, 1990; Ingram, 1995; Sugarman et al., 1995). The resulting growth of the continental ice sheet also resulted in a regression of the sea level. The global cooling, which appears to be associated with a decrease in the diversity of fossils, eventually led to the mass extinction at the K-Pg boundary.

We have obtained a steady-state, 1-D model of the accretion disk of a protostar taking into account the magneto-rotational instability (MRI). We find that the disk is divided into an outer turbulent region (OTR), a MRI suppressed region (MSR), and an inner turbulent region (ITR). The outer turbulent region is fully turbulent because of MRI. However, in the range, r_{out} (= 8 - 60 AU) from the central star, MRI is suppressed around the midplane of the gas disk and a quiet area without turbulence appears, because the degree of ionization of gas becomes low enough. The disk becomes fully turbulent again in the range r in (= 0.2 - 1 AU), which is called the inner turbulent region, because the midplane temperature become high enough (\gt 1000 K) due to gravitational energy release. Planetesimals are formed through gravitational instability at the two distinct sites, outer and inner MRI fronts (the boundaries between the MRI suppressed region (MSR) and the outer and inner turbulent regions), because of the radial concentration of the solid particles. At the outer MRI front, icy particles grow through low-velocity collisions into porous aggregates with low densities. They eventually undergo gravitational instability to form icy planetesimals. On the other hand, rocky particles accumulate at the inner MRI front, since their drift velocities turn outward due to the local maximum in gas pressure. They undergo gravitational instability in a sub-disk of pebbles to form rocky planetesimals at the inner MRI front. The tandem regime is consistent with the ABEL model, in which the Earth was initially formed as a completely volatile-free planet. The water and other volatile elements came later through the accretion of icy particles by the occasional scatterings in the outer regions.

Two types of orbital detectors of extreme energy cosmic rays are being developed nowadays: (i) TUS and KLYPVE with reflecting optical systems (mirrors) and (ii) JEM-EUSO with high-transmittance Fresnel lenses. They will cover much larger areas than existing ground-based arrays and almost uniformly monitor the celestial sphere. The TUS detector is the pioneering mission developed in SINP MSU in cooperation with several Russian and foreign institutions. It has relatively small field of view (+/-4.5 deg), which corresponds to a ground area of 6.4x10^3 sq.km. The telescope consists of a Fresnel-type mirror-concentrator (~2 sq.m) and a photo receiver (a matrix of 16x16 photomultiplier tubes). It is to be deployed on the Lomonosov satellite, and is currently at the final stage of preflight tests. Recently, SINP MSU began the KLYPVE project to be installed on board of the Russian segment of the ISS. The optical system of this detector contains a larger primary mirror (10 sq.m), which allows decreasing the energy threshold. The total effective field of view will be at least +/-14 degrees to exceed the annual exposure of the existing ground-based experiments. Several configurations of the detector are being currently considered. Finally, JEM-EUSO is a wide field of view (+/-30 deg) detector. The optics is composed of two curved double-sided Fresnel lenses with 2.65 m external diameter, a precision diffractive middle lens and a pupil. The ultraviolet photons are focused onto the focal surface, which consists of nearly 5000 multi-anode photomultipliers. It is developed by a large international collaboration. All three orbital detectors have multi-purpose character due to continuous monitoring of various atmospheric phenomena. The present status of development of the TUS and KLYPVE missions is reported, and a brief comparison of the projects with JEM-EUSO is given.

Reliability assessment in concerned with the analysis of devices and systems whose individual components are prone to fail. This reliability analysis documents the process and results of reliability determination of the JEM-EUSO photomultiplier tube component using the methods 217 Plus. Quantum efficiency degradation and radiation hardness assurance. In conclussion, the levels of damage suffered by the PMTs which comprise the focal surface of JEM-EUSO Space Telescope, are acceptable. The results show as well the greatest contribution to the failure is due to radiation SET. The guaranteed performance of this equipment is a 99.45 per cent, an accepted value of reliability thus fulfilling the objectives and technological challenges of JEM-EUSO.

An accreting supermassive blackhole, the central engine of active galactic nucleus (AGN), is capable of exciting extreme amplitude Alfven waves whose wavelength (wave packet) size is characterized by its clumpiness. The pondermotive force and wakefield are driven by these Alfven waves propagating in the AGN (blazar) jet and accelerate protons/nuclei to extreme energies beyond Zetta-electron volt (ZeV$=10^{21}$ eV). Such acceleration is prompt, localized, and does not suffer from the multiple scattering/bending enveloped in the Fermi acceleration that causes excessive synchrotron radiation loss beyond $10^{19}$ eV. The production rate of ZeV cosmic rays is found to be consistent with the observed gamma-ray luminosity function of blazars and their time variabilities.

We evaluate the exposure during nadir observations with JEM-EUSO, the Extreme Universe Space Observatory, on-board the Japanese Experiment Module of the International Space Station. Designed as a mission to explore the extreme energy Universe from space, JEM-EUSO will monitor the Earth's nighttime atmosphere to record the ultraviolet light from tracks generated by extensive air showers initiated by ultra-high energy cosmic rays. In the present work, we discuss the particularities of space-based observation and we compute the annual exposure in nadir observation. The results are based on studies of the expected trigger aperture and observational duty cycle, as well as, on the investigations of the effects of clouds and different types of background light. We show that the annual exposure is about one order of magnitude higher than those of the presently operating ground-based observatories.

Earth-mass dark matter microhalos with size of ~100 AUs are the first structures formed in the universe, if the dark matter of the Universe are made of neutralino. Here, we report the results of ultra-high-resolution simulations of the formation and evolution of these microhalos. We found that microhalos have the central density cusps of the form $\rho \propto r^{-1.5}$, much steeper than the cusps of larger dark halos. The central regions of these microhalos survive the encounters with stars except in very inner region of the galaxy down to the radius of a few hundreds pcs from the galactic center. The annihilation signals from nearest microhalos are observed as gamma-ray point-sources (radius less than 1'), with unusually large proper motions of ~0.2 degree per year. Their surface brightnesses are ~10% of that of the galactic center. Their S/N ratios might be better if they are far from the galactic plane. Luminosities of subhalos are determined only by their mass, and they are more than one order of magnitude luminous than the estimation by Springel et al. (2008): A boost factor can be larger than 1000. Perturbations to the millisecond pulsars by gravitational attractions of nearby earth-mass microhalos can be detected by the observations of Parkes Pulsar Timing Array (PPTA).

This document contains a summary of the workshop which took place on 22 - 24 February 2012 at the Kavli Institute of Cosmological Physics in the University of Chicago. The goal of the workshop was to discuss the physics reach of the JEM-EUSO mission and how best to implement a global ground based calibration system for the instrument to realize the physics goal of unveiling the origin of the highest energy cosmic rays.

Contributions of the JEM-EUSO Collaboration to the 32nd International Cosmic Ray Conference, Beijing, August, 2011.

The JEM-EUSO mission explores the origin of the extreme energy cosmic rays (EECRs) above 100 EeV and explores the limits of the fundamental physics through the observations of their arrival directions and energies. It is designed to achieve an exposure larger than 1 million km^2 sr year to open a new particle astronomy channel. This super-wide-field (60 degrees) telescope with a diameter of about 2.5 m looks down from space onto the night sky to detect near UV photons (330-400nm, both fluorescent and Cherenkov photons) emitted from the giant air showers produced by EECRs. The arrival direction map with more than five hundred events will tell us the origin of the EECRs and allow us to identify the nearest EECR sources with known astronomical objects. It will allow them to be examined in other astronomical channels. This is likely to lead to an understanding of the acceleration mechanisms, perhaps producing discoveries in astrophysics and/or fundamental physics. The comparison of the energy spectra among the spatially resolved individual sources will help to clarify the acceleration/emission mechanism, and also finally confirm the Greisen-Zatsepin-Kuz'min process for the validation of Lorentz invariance up to {\gamma}~10^{11}. Neutral components (neutrinos and gamma rays) can also be detected as well, if their fluxes are high enough. The JEM-EUSO mission is planned to be launched by a H2B rocket about JFY 2015-2016 and transferred to ISS by H2 Transfer Vehicle (HTV). It will be attached to the Exposed Facility external experiment platform of "KIBO".

Structure of hot dense matter at subnuclear densities is investigated by quantum molecular dynamics (QMD) simulations. We analyze nucleon distributions and nuclear shapes using two-point correlation functions and Minkowski functionals to determine the phase-separation line and to classify the phase of nuclear matter in terms of the nuclear structure. Obtained phase diagrams show that the density of the phase boundaries between the different nuclear structures decreases with increasing temperature due to the thermal expansion of nuclear matter region. The critical temperature for the phase separation is $\agt 6$ MeV for the proton fraction $x=0.5$ and $\agt 5$ MeV for $x=0.3$. Our result suggests the existence of "spongelike" phases with negative Euler characteristic in addition to the simple "pasta" phases in supernova cores until $T \alt 3$ MeV.

We examine the model dependence of the phase diagram of inhomogeneous nulcear matter in supernova cores using the quantum molecular dynamics (QMD). Inhomogeneous matter includes crystallized matter with nonspherical nuclei - "pasta" phases - and the liquid-gas phase separating nuclear matter. Major differences between the phase diagrams of the QMD models can be explained by the energy of pure neutron matter at low densities and the saturation density of asymmetric nuclear matter. We show the density dependence of the symmetry energy is also useful to understand uncertainties of the phase diagram. We point out that, for typical nuclear models, the mass fraction of the pasta phases in the later stage of the collapsing cores is higher than 10-20 %.

JEM-EUSO, on board of the Japanese Exploration Module of the International Space Station, is being proposed as the first space observatory devoted to UHECR. Its privileged position at 430 km above the Earth surface, combined with a large field of view, innovative optics and a high efficiency focal surface, results in an unprecedented exposure which significantly surpasses that of the largest ground observatories. The large number of events expected above the GZK threshold for photo-pion production by protons will allow the directional identification of individual sources and the determination of their spectra, i.e., doing astronomy and astrophysics through the particle channel. Similar goals can be achieved in the case of light UHECR nuclei. Furthermore, the atmospheric target volume ($\sim 10^{12}$ ton) makes the possibility of neutrino observation a highlight of the mission. Other exploratory objectives include the detection of extreme energy gammas and the study of Galactic magnetic fields as well as global observations of the earth's atmosphere, including clouds, night-glows, plasma discharges, and meteors. In this contribution we will describe the scientific objectives of JEM-EUSO.

The identification of very high energy photons is of great importance for the understanding of the origin of extreme energy cosmic rays (EECR). Several can be the sources of high energy photons at Earth. A guaranteed component is the flux of high energy photons expected as a consequence of the interaction of cosmic rays with the cosmic photon background. Another contribution may be expected as by-product at the acceleration sites of protons and nuclei, although such flux should be strongly suppressed for distant sources. On the other hand, top-down scenarios involving the decay of super heavy relic particles or topological defects, even if not currently favored, have as a characteristic signature an increasingly dominant flux of photons at the highest energies. In this work we study the statistical separation between hadron and photon showers at energies where both, LPM effect and magnetospheric interactions are important for the development of the cascades. We consider a detector with the same orbital characteristics as JEM-EUSO, but disregard trigger and reconstruction efficiencies, in order to define the maximum ideal discrimination power attainable.

Neutrinos are a very promising messenger at tens of EeV and above. They can be produced by several channels, namely as by products of hadronic interactions at the sources, as the main products of the decay of super massive particles and, in a guaranteed way, as the result of the propagation of UHECR through the bath of microwave relic photons. A new era of very large exposure space observatories, of which the JEM-EUSO mission is a prime example, is on the horizon and, with it, it is even larger the possibility of astrophysical neutrino detection at the highest energies. In the present work we use a combination of the PYTHIA interaction code with the CONEX shower simulation package in order to produce fast one-dimensional simulations of neutrino initiated showers in air. We make a detail study of the structure of the corresponding longitudinal profiles, but focus our physical analysis mainly on the development of showers at mid and high altitudes, where they can be an interesting target for space fluorescence observatories.

In supernova cores, nuclear "pasta" phases such as triangular lattice of rod-like nuclei and layered structure of slab-like nuclei are considered to exist. However, it is still unclear whether or not they are actually formed in collapsing supernova cores. Using {\it ab-initio} numerical simulations called the Quantum Molecular Dynamics (QMD), we here solve this problem by demonstrating that a lattice of rod-like nuclei is formed from a bcc lattice by compression. We also find that, in the transition process, the system undergoes zigzag configuration of elongated nuclei, which are formed by a fusion of two original spherical nuclei.

Magnetic elements are thought to be described by flux tube models, and are well reproduced by MHD simulations. However, these simulations are only partially constrained by observations. We observationally investigate the relationship between G-band bright points and magnetic structures to clarify conditions, which make magnetic structures bright in G-band. The G-band filtergrams together with magnetograms and dopplergrams were taken for a plage region covered by abnormal granules as well as ubiquitous G-band bright points, using the Swedish 1-m Solar Telescope (SST) under very good seeing conditions. High magnetic flux density regions are not necessarily associated with G-band bright points. We refer to the observed extended areas with high magnetic flux density as magnetic islands to separate them from magnetic elements. We discover that G-band bright points tend to be located near the boundary of such magnetic islands. The concentration of G-band bright points decreases with inward distance from the boundary of the magnetic islands. Moreover, G-band bright points are preferentially located where magnetic flux density is higher, given the same distance from the boundary. There are some bright points located far inside the magnetic islands. Such bright points have higher minimum magnetic flux density at the larger inward distance from the boundary. Convective velocity is apparently reduced for such high magnetic flux density regions regardless of whether they are populated by G-band bright points or not. The magnetic islands are surrounded by downflows.These results suggest that high magnetic flux density, as well as efficient heat transport from the sides or beneath, are required to make magnetic elements bright in G-band.

We investigate the evolution of collisionally merged stars with mass of ~100 MSun which might be formed in dense star clusters. We assumed that massive stars with several tens Msun collide typically after ~1Myr of the formation of the cluster and performed hydrodynamical simulations of several collision events. Our simulations show that after the collisions, merged stars have extended envelopes and their radii are larger than those in the thermal equilibrium states and that their interiors are He-rich because of the stellar evolution of the progenitor stars. We also found that if the mass-ratio of merging stars is far from unity, the interior of the merger product is not well mixed and the elemental abundance is not homogeneous. We then followed the evolution of these collision products by a one dimensional stellar evolution code. After an initial contraction on the Kelvin-Helmholtz (thermal adjustment) timescale (~1000-10000 yr), the evolution of the merged stars traces that of single homogeneous stars with corresponding masses and abundances, while the initial contraction phase shows variations which depend on the mass ratio of the merged stars. We infer that, once runaway collisions have set in, subsequent collisions of the merged stars take place before mass loss by stellar winds becomes significant. Hence, stellar mass loss does not inhibit the formation of massive stars with mass of ~1000Msun.

Nuclear ``pasta'', nonspherical nuclei in dense matter, is predicted to occur in collapsing supernova cores. We show how pasta phases affect the neutrino transport cross section via weak neutral current using several nuclear models. This is the first calculation of the neutrino opacity of the phases with rod-like and slab-like nuclei taking account of finite temperature effects, which are well described by the quantum molecular dynamics. We also show that pasta phases can occupy 10--20% of the mass of supernova cores in the later stage of the collapse.

Using the luminosity function and the apparent GRB rate inferred from the spectral peak energy-luminosity relation, we investigate the absolute GRB formation rate, taking into account the effects of both the jet-luminosity evolution and the jet opening angle evolution. In the case that there is no jet opening angle evolution, the jet-corrected luminosity for high-redshift GRBs is much larger than typical values for low-z (z>~1) bursts. On the other hand, if the jet-luminosity does not evolve with time, the jet opening angle for high-z bursts should be much smaller than that for low-z GRBs. Therefore, it is preferable to take into account both evolution effects. We also estimate the local GRB event rate in a galaxy of ~10^{-7}-10^{-5}/yr.

We simulate the inner 100pc of the Milky-Way Galaxy to study the formation and evolution of the population of star clusters and intermediate mass black holes. For this study we perform extensive direct N-body simulations of the star clusters which reside in the bulge, and of the inner few tenth of parsecs of the super massive black hole in the Galactic center. In our N-body simulations the dynamical friction of the star cluster in the tidal field of the bulge are taken into account via (semi)analytic soluations. The N-body calculations are used to calibrate a (semi)analytic model of the formation and evolution of the bulge. We find that about 10% of the clusters born within 100pc of the Galactic center undergo core collapse during their inward migration and form intermediate-mass black holes (IMBHs) via runaway stellar merging. After the clusters dissolve, these IMBHs continue their inward drift, carrying a few of the most massive stars with them. We predict that region within about 10 parsec of the SMBH is populated by about 50IMBHs of some 1000Msun. Several of these are expected to be accompanied still by some of the most massive stars from the star cluster. We also find that within a few milliparsec of the SMBH there is a steady population of several IMBHs. This population drives the merger rate between IMBHs and the SMBH at a rate of about one per 10Myr, sufficient to build the accumulate majority of mass of the SMBH. Mergers of IMBHs with SMBHs throughout the universe are detectable by LISA, at a rate of about two per week.

Resent observations and theoretical interpretations suggest that IMBHs (intermediate-mass black hole) are formed in the centers of young and compact star clusters born close to the center of their parent galaxy. Such a star cluster would sink toward the center of the galaxy, and at the same time stars are stripped out of the cluster by the tidal field of the parent galaxy. We investigated the orbital evolution of the IMBH, after its parent cluster is completely disrupted by the tidal field of the parent galaxy, by means of large-scale N-body simulations. We constructed a model of the central region of our galaxy, with an SMBH (supermassive black hole) and Bahcall-Wolf stellar cusp, and placed an IMBH in a circular orbit of radius 0.086pc. The IMBH sinks toward the SMBH through dynamical friction, but dynamical friction becomes ineffective when the IMBH reached the radius inside which the initial stellar mass is comparable to the IMBH mass. This is because the IMBH kicks out the stars. This behavior is essentially the same as the loss-cone depletion observed in simulations of massive SMBH binaries. After the evolution through dynamical friction stalled, the eccentricity of the orbit of the IMBH goes up, resulting in the strong reduction in the merging timescale through gravitational wave radiation. Our result indicates that the IMBHs formed close to the galactic center can merge with the central SMBH in short time. The number of merging events detectable with DECIGO is estimated to be around 50 per year. Event rate for LISA would be similar or less, depending on the growth mode of IMBHs.

Calculations of equilibrium properties of dense matter predict that at subnuclear densities nuclei can be rodlike or slablike. To investigate whether transitions between phases with non-spherical nuclei can occur during the collapse of a star, we perform quantum molecular dynamic simulations of the compression of dense matter. We have succeeded in simulating the transitions between rodlike and slablike nuclei and between slablike nuclei and cylindrical bubbles. Our results strongly suggest that non-spherical nuclei can be formed in the inner cores of collapsing stars.

We have followed the evolution of multi-mass star clusters containing massive central black holes through collisional N-body simulations done on GRAPE6. Each cluster is composed of between 16,384 to 131,072 stars together with a black hole with an initial mass of M_BH=1000 Msun. We follow the evolution of the clusters under the combined influence of two-body relaxation, stellar mass-loss and tidal disruption of stars. The (3D) mass density profile follows a power-law distribution \rho \sim r^{-\alpha} with slope \alpha=1.55. This leads to a constant density profile of bright stars in projection, which makes it highly unlikely that core collapse clusters contain intermediate-mass black holes (IMBHs). Instead globular clusters containing IMBHs can be fitted with standard King profiles. The disruption rate of stars is too small to form an IMBH out of a M_BH \approx 50 Msun progenitor black hole, unless a cluster starts with a central density significantly higher than what is seen in globular clusters. Kinematical studies can reveal 1000 Msun IMBHs in the closest clusters. IMBHs in globular clusters are only weak X-ray sources since the tidal disruption rate of stars is low and the star closest to the IMBH is normally another black hole. For globular clusters, dynamical evolution can push compact stars near the IMBH to distances small enough that they become detectable through gravitational radiation. If 10% of all globular clusters contain IMBHs, extragalactic globular clusters could be one of the major sources for {\it LISA}. (abridged)

In this paper we report results of collisional N-body simulations of the dynamical evolution of equal-mass star clusters containing a massive central black hole. Each cluster is composed of between 5,000 to 180,000 stars together with a central black hole which contains between 0.2% to 10% of the total cluster mass. We find that for large enough black hole masses, the central density follows a power-law distribution with slope \rho \sim r^{-1.75} inside the radius of influence of the black hole, in agreement with predictions from earlier Fokker Planck and Monte Carlo models. The tidal disruption rate of stars is within a factor of two of that derived in previous studies. It seems impossible to grow an intermediate-mass black hole (IMBH) from a M \le 100 Msun progenitor in a globular cluster by the tidal disruption of stars, although M = 10^3 Msun IMBHs can double their mass within a Hubble time in dense globular clusters. The same is true for the supermassive black hole at the centre of the Milky Way. Black holes in star clusters will feed mainly on stars tightly bound to them and the re-population of these stars causes the clusters to expand, reversing core-collapse without the need for dynamically active binaries. Close encounters of stars in the central cusp also lead to an increased mass loss rate in the form of high-velocity stars escaping from the cluster. A companion paper will extend these results to the multi-mass case.

We have followed the evolution of multi-mass star clusters containing massive central black holes by N-body simulations on the GRAPE6 computers of Tokyo University. We find a strong cluster expansion and significant structural changes of the clusters. Star clusters with IMBHs have power-law density profiles $rho ~ r^{-alpha}$ with slopes $alpha=1.55$ inside the influence sphere of the central black hole. This leads to a constant density profile of bright stars in projection, which rules out the presence of intermediate mass black holes in core collapse clusters. If the star clusters are surrounded by a tidal field, a central IMBH speeds up the destruction of the cluster until a remnant of a few hundred stars remains, which stays bound to the IMBH for a long time. We also discuss the efficiency of different detection mechanisms for finding IMBHs in star clusters.

Structure of cold and hot dense matter at subnuclear densities is investigated by quantum molecular dynamics (QMD) simulations. Obtained phase diagrams show that the density of the phase boundaries between the different nuclear structures decreases with increasing temperature due to the thermal expansion of nuclear matter region.

Structure of cold nuclear matter at subnuclear densities for the proton fraction $x=0.5$, 0.3 and 0.1 is investigated by quantum molecular dynamics (QMD) simulations. We demonstrate that the phases with slablike and rodlike nuclei, etc. can be formed dynamically from hot uniform nuclear matter without any assumptions on nuclear shape, and also systematically analyze the structure of cold matter using two-point correlation functions and Minkowski functionals. In our simulations, we also observe intermediate phases, which has complicated nuclear shapes. It is found out that these phases can be characterized as those with negative Euler characteristic. Our result implies the existence of these kinds of phases in addition to the simple ``pasta'' phases in neutron star crusts and supernova inner cores. In addition, we investigate the properties of the effective QMD interaction used in the present work to examine the validity of our results. The resultant energy per nucleon $\epsilon_{n}$ of the pure neutron matter, the proton chemical $\mu_{p}^{(0)}$ in pure neutron matter and the nuclear surface tension $E_{\rm surf}$ are generally reasonable in comparison with other nuclear interactions.

Structure of cold dense matter at subnuclear densities is investigated by quantum molecular dynamics (QMD) simulations. We succeeded in showing that the phases with slab-like and rod-like nuclei etc. can be formed dynamically from hot uniform nuclear matter without any assumptions on nuclear shape. We also observe intermediate phases, which has complicated nuclear shapes. Geometrical structures of matter are analyzed with Minkowski functionals, and it is found out that intermediate phases can be characterized as ones with negative Euler characteristic. Our result suggests the existence of these kinds of phases in addition to the simple ``pasta'' phases in neutron star crusts.

Structure of cold dense stellar matter at subnuclear densities is investigated with quantum molecular dynamics (QMD). We succeeded in showing that the phases with slab-like and rod-like nuclei etc. can be formed dynamically without any assumptions on the nuclear shape. Our result suggests the existence of these kinds of phases in neutron star crusts.

The TAROT telescope has for primary goal the search for the prompt optical counterpart of Cosmic Gamma-Ray Bursts. It is a completely autonomous 25cm telescope installed near Nice (France), able to point any location of the sky within 1-2 seconds. The control, scheduling, and data processing activities are completely automated, so the instrument is completely autonomous. In addition to its un-manned modes, we added recently the possibility to remotely control the telescope, as a request of the "Hands-On Universe" (HOU) program for exchange of time within automatic telescopes for the education and public outreach. To this purpose we developed a simple control interface. A webcam was installed to visualize the telescope. Access to the data is possible through a web interface. The images can be processed by the HOU software, a program specially suited for use within the classroom. We experienced these feature during the open days of the University of California Berkeley and the Astronomy Festival of Fleurance (France). We plan a regular use for an astronomy course of the Museum of Tokyo, as well as for French schools. Not only does Hands-On TAROT gives the general public an access to professional astronomy, but it is also a more general tool to demonstrate the use of a complex automated system, the techniques of data processing and automation. Last but not least, through the use of telescopes located in many countries over the globe, a form of powerful and genuine cooperation between teachers and children from various countries is promoted, with a clear educational goal.

Hands-On Universe (HOU) is an educational program that enables students to investigate the Universe while applying tools and concepts from science, math, and technology. Using the Internet, HOU participants around the world request observations from an automated telescope, download images from a large image archive, and analyze them with the aid of user-friendly image processing software. This program is developing now in many countries, including the USA, France, Germany, Sweden, Japan, Australia, and others. A network of telescopes has been established among these countries, many of them remotely operated, as shown in the accompanying demo. Using this feature, students in the classroom are able to make night observations during the day, using a telescope placed in another country. An archive of images taken on large telescopes is also accessible, as well as resources for teachers. Students are also dealing with real research projects, e.g. the search for asteroids, which resulted in the discovery of a Kuiper Belt object by high-school students. Not only Hands-On Universe gives the general public an access to professional astronomy, but it is also a more general tool to demonstrate the use of a complex automated system, the techniques of data processing and automation. Last but not least, through the use of telescopes located in many countries over the globe, a form of powerful and genuine cooperation between teachers and children from various countries is promoted, with a clear educational goal.

Observations of stellar kinematics, gas dynamics and masers around galactic nuclei have now firmly established that many galaxies host central supermassive black holes (SMBHs) with masses in the range $10^6 \sim 10^9$M$_{\odot}$. However, how these SMBHs formed is not well understood. One reason for this situation is the lack of observations of intermediate-mass BHs (IMBHs), which could bridge the gap between stellar-mass BHs and SMBHs. Recently, this missing link (i.e., an IMBH) has been found in observations made by the ASCA and the Chandra of the central region of the starburst galaxy M82 \citep{MT99, PG99, MT01, Ka01}. Subsequent observations by SUBARU have revealed that this IMBH apparently coincides with a young compact star cluster. Based on these findings, we suggest a new formation scenario for SMBHs. In this scenario, IMBHs first form in young compact star clusters through runaway merging of massive stars. While these IMBHs are forming, the host star clusters sink toward the galactic nucleus through dynamical friction, and upon evaporation deposit their IMBHs near the galactic center. The IMBHs then form binaries and eventually merge via gravitational radiation, forming an SMBH.

The homogeneity of the cosmic microwave background radiation (CBR) is one of the most severe constraint for theories of the structure formation in the universe. We investigated the effect of the gravitational scattering (lensing) of galaxies, clusters of galaxies, and superclusters on the anisotropy of the CBR by numerical simulations. Although this effect was thought to be unimportant, we found that the gravitational scatterings by superclusters can significantly reduce the anisotropy of the CBR. We took into account the exponential growth of the distance between two rays due to multiple scatterings. The bending angle of each ray grows through the random walk process. On the other hand, difference between two rays grows exponentially while it is small. This exponential growth is caused by coherent scatterings that two rays suffer, and was neglected in the previous studies. The gravitational scattering by superclusters reduces the observed temperature anisotropy of the CBR at present time approximately by 40--60 $\%$ from that at the recombination time for angular scale up to a few degrees if the supercluster were formed at $z=$2--4.