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.

Using a prototype of the Tomo-e Gozen wide-field CMOS mosaic camera, we acquire wide-field optical images at a cadence of 2 Hz and search them for transient sources of duration 1.5 to 11.5 seconds. Over the course of eight nights, our survey encompasses the equivalent of roughly two days on one square degree, to a fluence equivalent to a limiting magnitude about $V = 15.6$ in a 1-second exposure. After examining by eye the candidates identified by a software pipeline, we find no sources which meet all our criteria. We compute upper limits to the rate of optical transients consistent with our survey, and compare those to the rates expected and observed for representative sources of ephemeral optical light.

We report observations of a stellar occultation by the classical Kuiper belt object (50000) Quaoar occurred on 28 June 2019. A single-chord high-cadence (2 Hz) photometry dataset was obtained with the Tomo-e Gozen CMOS camera mounted on the 1.05 m Schmidt telescope at Kiso Observatory. The obtained ingress and egress data do not show any indication of atmospheric refraction and allow to set new $1\sigma$ and $3\sigma$ upper limits of 6 and 16 nbar, respectively, for the surface pressure of a pure methane atmosphere. These upper limits are lower than the saturation vapor pressure of methane at Quaoar's expected mean surface temperature ($T \sim 44$ K) and imply the absence of a $\sim$10 nbar-level global atmosphere formed by methane ice on Quaoar's surface.

We present a fireball detected in the night sky over Kyoto, Japan on UT 2017 April 28 at ${\rm 15^{h}\,58^{m}\,19^{s}}$ by the SonotaCo Network. The absolute visual magnitude is $M_{\rm v}$=$-$4.10$\pm$0.42mag. Luminous light curves obtain a meteoroid mass $m$=29$\pm$1g, corresponding to the size $a_{\rm s}$=2.7$\pm$0.1cm. Orbital similarity assessed by D-criterions (cf. $D_{\rm SH}$=0.0079) has identified a likely parent, the binary near-Earth asteroid (164121) 2003 YT$_1$. The suggested binary formation process is a YORP-driven rotational disintegration (Pravec & Harris 2007). The asynchronous state indicates the age of $<$10$^4$yr, near or shorter than the upper limit to meteoroid stream lifetime. We examine potential dust production mechanisms for the asteroid, including rotational instability, resurfacing, impact, photoionization, radiation pressure sweeping, thermal fracture and sublimation of ice. We find some of them capable of producing the meteoroid-scale particles. Rotational instability is presumed to cause mass shedding, in consideration of the recent precedents (e.g. asteroid (6478) Gault), possibly releasing mm-cm scale dust particles. Impacts by micrometeorites with size $\simeq$1mm could be a trigger for ejecting the cm-sized particles. Radiation pressure can sweep out the mm-sized dust particles, while not sufficient for the cm-sized. For the other mechanisms, unprovable or unidentified. The feasibility in the parental aspect of 2003 YT$_1$ is somewhat reconciled with the fireball observation, yielding an insight into how we approach potentially hazardous objects.

Imaging observations of faint meteors were carried out on April 11 and 14, 2016 with a wide-field CMOS mosaic camera, Tomo-e PM, mounted on the 105-cm Schmidt telescope at Kiso Observatory, the University of Tokyo. Tomo-e PM, which is a prototype model of Tomo-e Gozen, can monitor a sky of ${\sim}1.98\,\mathrm{deg^2}$ at 2\,Hz. The numbers of detected meteors are 1514 and 706 on April 11 and 14, respectively. The detected meteors are attributed to sporadic meteors. Their absolute magnitudes range from $+4$ to $+10\,\mathrm{mag}$ in the $V$-band, corresponding to about $8.3{\times}10^{-2}$ to $3.3{\times}10^{-4}\,\mathrm{g}$ in mass. The present magnitude distributions we obtained are well explained by a single power-law luminosity function with a slope parameter $r = 3.1{\pm}0.4$ and a meteor rate $\log_{10}N_0 = -5.5{\pm}0.5$. The results demonstrate a high performance of telescopic observations with a wide-field video camera to constrain the luminosity function of faint meteors. The performance of Tomo-e Gozen is about two times higher than that of Tomo-e PM. A survey with Tomo-e Gozen will provide a more robust measurement of the luminosity function.

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.

We present a unique and significant polarimetric result regarding the near-Earth asteroid (152679) 1998 KU$_\mathrm{2}$ , which has a very low geometric albedo. From our observations, we find that the linear polarization degrees of 1998 KU$_\mathrm{2}$ are 44.6 $\pm$ 0.5\% in the R$_\mathrm{C}$ band and 44.0 $\pm$ 0.6\% in the V band at a solar phase angle of 81.0\degr. These values are the highest of any known airless body in the solar system (i.e., high-polarization comets, asteroids, and planetary satellites) at similar phase angles. This polarimetric observation is not only the first for primitive asteroids at large phase angles, but also for low-albedo (< 0.1) airless bodies. Based on spectroscopic similarities and polarimetric measurements of materials that have been sorted by size in previous studies, we conjecture that 1998 KU$_\mathrm{2}$ has a highly microporous regolith structure comprising nano-sized carbon grains on the surface.

Lightning and thundercloud are the most dramatic natural particle accelerators on the Earth. Relativistic electrons accelerated by electric fields therein emit bremsstrahlung gamma rays, which have been detected at ground observations, by airborne detectors, and as terrestrial gamma-ray flashes (TGFs) from space. The energy of the gamma rays is sufficiently high to potentially invoke atmospheric photonuclear reactions 14N(gamma, n)13N, which would produce neutrons and eventually positrons via beta-plus decay of generated unstable radioactive isotopes, especially 13N. However, no clear observational evidence for the reaction has been reported to date. Here we report the first detection of neutron and positron signals from lightning with a ground observation. During a thunderstorm on 6 February 2017 in Japan, a TGF-like intense flash (within 1 ms) was detected at our monitoring sites 0.5-1.7 km away from the lightning. The subsequent initial burst quickly subsided with an exponential decay constant of 40-60 ms, followed by a prolonged line emission at about 0.511 megaelectronvolt (MeV), lasting for a minute. The observed decay timescale and spectral cutoff at about 10 MeV of the initial emission are well explained with de-excitation gamma rays from the nuclei excited by neutron capture. The centre energy of the prolonged line emission corresponds to the electron-positron annihilation, and hence is the conclusive indication of positrons produced after the lightning. Our detection of neutrons and positrons is unequivocal evidence that natural lightning triggers photonuclear reactions. No other natural event on the Earth is known to trigger photonuclear reactions. This discovery places lightning as only the second known natural channel on the Earth after the atmospheric cosmic-ray interaction, in which isotopes, such as 13C, 14C, and 15N, are produced.

We conducted a polarimetric observation of the fast-rotating near-Earth asteroid (1566) Icarus at large phase (Sun-asteroid-observer's) angles $\alpha$= 57 deg--141deg around the 2015 summer solstice. We found that the maximum values of the linear polarization degree are $P_\mathrm{max}$=7.32$\pm$0.25 % at phase angles of $\alpha_\mathrm{max}$=124$\pm$8 deg in the $V$-band and $P_\mathrm{max}$=7.04$\pm$0.21 % at $\alpha_\mathrm{max}$=124$\pm$6 deg in the $R_\mathrm{C}$-band. Applying the polarimetric slope-albedo empirical law, we derived a geometric albedo of $p_\mathrm{V}$=0.25$\pm$0.02, which is in agreement with that of Q-type taxonomic asteroids. $\alpha_\mathrm{max}$ is unambiguously larger than that of Mercury, the Moon, and another near-Earth S-type asteroid (4179) Toutatis but consistent with laboratory samples with hundreds of microns in size. The combination of the maximum polarization degree and the geometric albedo is in accordance with terrestrial rocks with a diameter of several hundreds of micrometers. The photometric function indicates a large macroscopic roughness. We hypothesize that the unique environment (i.e., the small perihelion distance $q$=0.187 au and a short rotational period of $T_\mathrm{rot}$=2.27 hours) may be attributed to the paucity of small grains on the surface, as indicated on (3200) Phaethon.

Recent submillimeter and far-infrared wavelength observations of absorption in the rotational ground-state lines of various simple molecules against distant Galactic continuum sources have opened the possibility of studying the chemistry of diffuse molecular clouds throughout the Milky Way. In order to calculate abundances, the column densities of molecular and atomic hydrogen, HI, must be known. We aim at determining the atomic hydrogen column densities for diffuse clouds located on the sight lines toward a sample of prominent high-mass star-forming regions that were intensely studied with the HIFI instrument onboard Herschel. Based on Jansky Very Large Array data, we employ the 21 cm HI absorption-line technique to construct profiles of the HI opacity versus radial velocity toward our target sources. These profiles are combined with lower resolution archival data of extended HI emission to calculate the HI column densities of the individual clouds along the sight lines. We employ Bayesian inference to estimate the uncertainties of the derived quantities. Our study delivers reliable estimates of the atomic hydrogen column density for a large number of diffuse molecular clouds at various Galactocentric distances. Together with column densities of molecular hydrogen derived from its surrogates observed with HIFI, the measurements can be used to characterize the clouds and investigate the dependence of their chemistry on the molecular fraction, for example.

We present the early-phase spectra and the light curves of the broad-lined supernova (SN) 2016coi from $t=7$ to $67$ days after the estimated explosion date. This SN was initially reported as a broad-lined Type SN Ic (SN Ic-BL). However, we found that spectra up to $t=12$ days exhibited the He~{\sc i} $\lambda$5876, $\lambda$6678, and $\lambda$7065 absorption lines. We show that the smoothed and blueshifted spectra of normal SNe Ib are remarkably similar to the observed spectrum of SN 2016coi. The line velocities of SN 2016coi were similar to those of SNe Ic-BL and significantly faster than those of SNe Ib. Analyses of the line velocity and light curve suggest that the kinetic energy and the total ejecta mass of SN 2016coi are similar to those of SNe Ic-BL. Together with broad-lined SNe 2009bb and 2012ap for which the detection of He~{\sc i} were also reported, these SNe could be transitional objects between SNe Ic-BL and SNe Ib, and be classified as broad-lined Type `Ib' SNe (SNe `Ib'-BL). Our work demonstrates the diversity of the outermost layer in broad-lined SNe, which should be related to the variety of the evolutionary paths.

Multiple outbursts of a Jupiter-family comet, 15P/Finlay, occurred from late 2014 to early 2015. We conducted an observation of the comet after the first outburst and subsequently witnessed another outburst on 2015 January 15.6-15.7. The gas, consisting mostly of C2 and CN, and dust particles expanded at speeds of 1,110 +/- 180 m/s and 570 +/- 40 m/s at a heliocentric distance of 1.0 AU. We estimated the maximum ratio of solar radiation pressure with respect to the solar gravity beta_max = 1.6 +/- 0.2, which is consistent with porous dust particles composed of silicates and organics. We found that 10^8-10^9 kg of dust particles (assumed to be 0.3 micron - 1 mm) were ejected through each outburst. Although the total mass is three orders of magnitude smaller than that of the 17P/Holmes event observed in 2007, the kinetic energy per unit mass (104 J/kg) is equivalent to the estimated values of 17P/Holmes and 332P/2010 V1 (Ikeya-Murakami), suggesting that the outbursts were caused by a similar physical mechanism. From a survey of cometary outbursts on the basis of voluntary reports, we conjecture that 15P/Finlay-class outbursts occur >1.5 times annually and inject dust particles from Jupiter-family comets and Encke-type comets into interplanetary space at a rate of ~10 kg/s or more.

Remote sensing data from orbiter missions have proposed that ground ice may currently exist on Mars, although the volume is still uncertain. Recent analyses of Martian meteorites have suggested that the water reservoirs have at least three distinct hydrogen isotope compositions (D/H ratios): primordial and high D/H ratios, which are approximately the same and six times that of ocean water on Earth, respectively, and a newly identified intermediate D/H ratio, which is approximately two to three times higher than that in ocean water on Earth. We calculate the evolution of the D/H ratios and the volumes of the water reservoirs on Mars by modeling the exchange of hydrogen isotopes between multiple water reservoirs and the atmospheric escape. The D/H ratio is slightly higher in the topmost thin surface-ice layer than that in the atmosphere because of isotopic fractionation by sublimation, whereas the water-ice reservoir just below the exchangeable topmost surface layer retains the intermediate D/H signature found in Martian meteorites. We propose two possible models for constraining the volume of the ground ice considering the observed D/H ratios and geomorphological estimates of Paleo-oceans. The first assumes that the atmospheric loss is dominated by the Jeans escape. In this case, the volume of ground ice should be larger than the total volume of the observable surface ice that mainly occurs as polar layered ice deposits. The other model assumes diffusion-limited atmospheric loss in which the interactive evolution of the multiple water reservoirs naturally accounts for the observed D/H ratios. In this scenario, a large volume of ground ice does not necessarily exist currently on Mars as opposed to the perspective view proposed on the basis of recent orbiter missions.

We report measurements of trigonometric parallaxes for six high-mass star-forming regions in the Scutum spiral arm of the Milky Way as part of the BeSSeL Survey. Combining our measurements with 10 previous measurements from the BeSSeL Survey yields a total sample of 16 sources in the Scutum arm with trigonometric parallaxes in the Galactic longitude range from 5 deg to 32 deg. Assuming a logarithmic spiral model, we estimate a pitch angle of 19.8 (+/-3.1) deg for the Scutum arm, which is larger than pitch angles reported for other spiral arms. The high pitch angle of the arm may be due to the arm's proximity to the Galactic bar. The Scutum arm sources show an average peculiar motion of 4 km/s slower than the Galactic rotation and 8 km/s toward the Galactic center. While the direction of this non-circular motion has the same sign as determined for sources in other spiral arms, the motion toward the Galactic center is greater for the Scutum arm sources.

We report measurements of parallaxes and proper motions of ten high-mass star-forming regions in the Sagittarius spiral arm of the Milky Way as part of the BeSSeL Survey with the VLBA. Combining these results with eight others from the literature, we investigated the structure and kinematics of the arm between Galactocentric azimuth around -2 and 65 deg. We found that the spiral pitch angle is 7.3 +- 1.5 deg; the arm's half-width, defined as the rms deviation from the fitted spiral, is around 0.2 kpc; and the nearest portion of the Sagittarius arm is 1.4 +- 0.2 kpc from the Sun. Unlike for adjacent spiral arms, we found no evidence for significant peculiar motions of sources in the Sagittarius arm opposite to Galactic rotation.

Martian surface morphology implies that Mars was once warm enough to maintain persistent liquid water on its surface. While the high D/H ratios (~6 times the Earth's ocean water) of the current martian atmosphere suggest that significant water has been lost from the surface during martian history, the timing, processes, and the amount of the water loss have been poorly constrained. Recent technical developments of ion-microprobe analysis of martian meteorites have provided accurate estimation of hydrogen isotope compositions (D/H) of martian water reservoirs at the time when the meteorites formed. Based on the D/H data from the meteorites, this study demonstrates that the water loss during the pre-Noachian (>41-99 m global equivalent layers, GEL) was more significant than in the rest of martian history (>10-53 m GEL). Combining our results with geological and geomorphological evidence for ancient oceans, we propose that undetected subsurface water/ice (~100-1000 m GEL) should have existed, and it exceeds the observable present water inventory (~20-30 m GEL) on Mars.

Over 100 trigonometric parallaxes and proper motions for masers associated with young, high-mass stars have been measured with the BeSSeL Survey, a VLBA key science project, the EVN, and the Japanese VERA project. These measurements provide strong evidence for the existence of spiral arms in the Milky Way, accurately locating many arm segments and yielding spiral pitch angles ranging from 7 to 20 degrees. The widths of spiral arms increase with distance from the Galactic center. Fitting axially symmetric models of the Milky Way with the 3-D position and velocity information and conservative priors for the solar and average source peculiar motions, we estimate the distance to the Galactic center, Ro, to be 8.34 +/- 0.16 kpc, a circular rotation speed at the Sun, To, to be 240 +/- 8 km/s, and a rotation curve that is nearly flat (a slope of -0.2 +/- 0.4 km/s/kpc) between Galactocentric radii of 5 and 16 kpc. Assuming a "universal" spiral galaxy form for the rotation curve, we estimate the thin disk scale length to be 2.44 +/- 0.16 kpc. The parameters Ro and To are not highly correlated and are relatively insensitive to different forms of the rotation curve. Adopting a theoretically motivated prior that high-mass star forming regions are in nearly circular Galactic orbits, we estimate a global solar motion component in the direction of Galactic rotation, Vsun = 14.6 +/- 5.0 km/s. While To and Vsun are significantly correlated, the sum of these parameters is well constrained, To + Vsun = 255.2 +/- 5.1 km/s, as is the angular speed of the Sun in its orbit about the Galactic center, (To + Vsun)/Ro = 30.57 +/- 0.43 km/s/kpc. These parameters improve the accuracy of estimates of the accelerations of the Sun and the Hulse-Taylor binary pulsar in their Galactic orbits, significantly reducing the uncertainty in tests of gravitational radiation predicted by general relativity.

We report trigonometric parallax measurements of masers in the massive star forming complex W43 from VLBA observations as part of the BeSSeL Survey. Based on measurements of three 12 GHz methanol maser sources (G029.86-00.04, G029.95-00.01 and G031.28+00.06) and one 22 GHz water maser source (G031.58+00.07) toward W43, we derived a distance of $5.49^{+0.39}_{-0.34}$ kpc to W43. By associating the masers with CO molecular clouds, and associating the clouds kinematically with CO longitude-velocity spiral features, we assign W43 to the Scutum spiral arm, close to the near end of the Galactic bar. The peculiar motion of W43 is about 20 km/s toward the Galactic Center and is very likely induced by the gravitational attraction of the bar.

We report four trigonometric parallaxes for high-mass star-forming regions within 4 kpc of the Galactic center. These measurements were made with the VLBA as part of the BeSSeL Survey. By associating these sources kinematically with large-scale features in CO and HI longitude-velocity diagrams, we begin to outline some major features of the inner Milky Way: the Connecting arm, the near and far 3 kpc arms, and the Norma arm. The Connecting arm in the first Galactic quadrant lies closer to the Galactic center than the far 3 kpc arm and is offset by the long-bar's major axis near its leading edge, supporting the presence of an inner Lindblad resonance. Assuming the 3 kpc arms are a continuous physical structure, the relative Galactocentric distance of its near and far sides suggests highly elliptical streamlines of gas around the bar(s) and a bar corotation radius, r_CR > 3.6 kpc. At a Galactic longitude near 10{\deg} and a heliocentric distance of about 5 kpc, the near 3 kpc arm and the Norma arm intersect on a face-on view of our Galaxy, while passing at different Galactic latitudes. We provide an accurate distance measurement to the W31 star-forming complex of 4.95 (+0.5;-0.43) kpc from the Sun, which associates it with a bright CO feature belonging to the near 3 kpc arm.

We present the first measurement of the absolute proper motions of IRAS 00259+5625 (CB3, LBN594) associated with the HI loop called the "NGC281 superbubble" that extends from the Galactic plane over ~300 pc toward decreasing galactic latitude. The proper motion components measured with VERA are (mu_alpha cos(delta), mu_delta) = (-2.48 +/- 0.32, -2.85 +/- 0.65) mas yr^{-1}, converted into (mu_l cos(b), mu_b) = (-2.72 +/- 0.32, -2.62 +/- 0.65) mas yr^{-1} in the Galactic coordinates. The measured proper motion perpendicular to the Galactic plane (mu_b) shows vertical motion away from the Galactic plane with a significance of about ~4-sigma. As for the source distance, the distance measured with VERA is marginal, 2.4^{+1.0}_{-0.6} kpc. Using the distance, an absolute vertical motion (v_{b}) of -17.9 +/- 12.2 km s^{-1} is determined with ~1.5-sigma significance. The tendency of the large vertical motion is consistent with previous VLBI results for NGC 281 associated with the same superbubble. Thus, our VLBI results indicate the superbubble expansion motion whose origin is believed to be sequential supernova explosions.

It is crucial to understand and model a behavior of galaxy biasing for future ambitious galaxy redshift surveys. Using 40 large cosmological N-body simulations for a standard LambdaCDM cosmology, we study the cross-correlation coefficient between matter and the halo density field, which is an indicator of the stochasticity of bias, over a wide redshift range 0\le z \le 3. The cross-correlation coefficient is important to extract information on the matter density field, e.g., by combining galaxy clustering and galaxy-galaxy lensing measurements. We compare the simulation results with integrated perturbation theory (iPT) proposed by one of the present authors and standard perturbation theory (SPT) combined with a phenomenological model of local bias. The cross-correlation coefficient derived from the iPT agrees with N-body simulation results down to r~15 (10) h^{-1}Mpc within 0.5 (1.0) % for all redshifts and halo masses we consider. The SPT with local bias does not explain complicated behaviors on quasilinear scales at low redshifts, while roughly reproduces the general behavior of the cross-correlation coefficient on fully nonlinear scales. The iPT is powerful to predict the cross-correlation coefficient down to quasilinear regimes with a high precision.

We study how well the Gaussian approximation is valid for computing the covariance matrices of the convergence power and bispectrum in weak gravitational lensing analyses. We focus on its impact on the cosmological parameter estimations by comparing the results with and without non-Gaussian error contribution in the covariance matrix. We numerically derive the covariance matrix as well as the cosmology dependence of the spectra from a large set of N-body simulations performed for various cosmologies and carry out Fisher matrix forecasts for tomographic weak lensing surveys with three source redshifts. After showing the consistency of the power and bispectra measured from our simulations with the state-of-the-art fitting formulas, we investigate the covariance matrix assuming a typical ongoing survey across 1500 deg^2 with the mean source number density of 30 arcmin^{-2} at the mean redshift z_s=1.0. Although the shape noise contributes a significant fraction to the total error budget and it mitigates the impact of the non-Gaussian error for this source number density, we find that the non-Gaussian error degrades the cumulative signal-to-noise ratio up to the maximum multipole of 2000 by a factor of about 2 (3) in the power (bi-) spectrum analysis. Its impact on the final cosmological parameter forecast with 6 parameters can be as large as 15% in the size of the one-dimensional statistical error. This can be a problem in future wide and deep weak lensing surveys for precision cosmology. We also show how much the dark energy figure of merit is affected by the non-Gaussian error contribution and demonstrate an optimal survey design with a fixed observational time.

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.

We report on a measurement of the temperature of the cosmic microwave background radiation field, T_CMB, at z = 0.88582 by imaging HC3N (3-2) and (5-4) absorption in the foreground galaxy of the gravitationally lens magnified radio source PKS 1830-211 using the Very Long Baseline Array and the phased Very Large Array. Low-resolution imaging of the data yields a value of Trot = 5.6+2.5-0.9 K, for the rotational temperature, Trot, which is consistent with the temperature of the cosmic microwave background at the absorber's redshift of 2.73(1+z) K. However, our high-resolution imaging reveals that the absorption peak position of the foreground gas is offset from the continuum peak position of the synchrotron radiation from PKS 1830-211 SW, which indicates that the absorbing cloud is covering only part of the emission from PKS 1830-211, rather than the entire core-jet region. This changes the line-to-continuum ratios, and we find Trot between 1.1 and 2.5 K, which is lower than the expected value. This shows that previous, Trot, measurements could be biased due to unresolved structure.

We present analyses to determine the fundamental parameters of the Galaxy based on VLBI astrometry of 52 Galactic maser sources obtained with VERA, VLBA and EVN. We model the Galaxy's structure with a set of parameters including the Galaxy center distance R_0, the angular rotation velocity at the LSR Omega_0, mean peculiar motion of the sources with respect to Galactic rotation (U_src, V_src, W_src), rotation-curve shape index, and the V component of the Solar peculiar motions V_sun. Based on a Markov chain Monte Carlo method, we find that the Galaxy center distance is constrained at a 5% level to be R_0 = 8.05 +/- 0.45 kpc, where the error bar includes both statistical and systematic errors. We also find that the two components of the source peculiar motion U_src and W_src are fairly small compared to the Galactic rotation velocity, being U_src = 1.0 +/- 1.5 km/s and W_src = -1.4 +/- 1.2 km/s. Also, the rotation curve shape is found to be basically flat between Galacto-centric radii of 4 and 13 kpc. On the other hand, we find a linear relation between V_src and V_sun as V_src = V_sun -19 (+/- 2) km/s, suggesting that the value of V_src is fully dependent on the adopted value of V_sun. Regarding the rotation speed in the vicinity of the Sun, we also find a strong correlation between Omega_0 and V_sun. We find that the angular velocity of the Sun, Omega_sun, which is defined as Omega_sun = Omega_0 + V_sun/R_0, can be well constrained with the best estimate of Omega_sun = 31.09 +/- 0.78 km/s/kpc. This corresponds to Theta_0 = 238 +/- 14 km/s if one adopts the above value of R_0 and recent determination of V_sun ~ 12 km/s.

Based on a suite of state-of-the-art high-resolution $N$-body simulations, we revisit the so-called halofit model (Smith et al. 2003) as an accurate fitting formula for the nonlinear matter power spectrum. While the halofit model has been frequently used as a standard cosmological tool to predict the nonlinear matter power spectrum in a universe dominated by cold dark matter, its precision has been limited by the low-resolution of $N$-body simulations used to determine the fitting parameters, suggesting the necessity of improved fitting formula at small scales for future cosmological studies. We run high-resolution $N$-body simulations for 16 cosmological models around the Wilkinson Microwave Anisotropy Probe (WMAP) best-fit cosmological parameters (1, 3, 5, and 7 year results), including dark energy models with a constant equation of state. The simulation results are used to re-calibrate the fitting parameters of the halofit model so as to reproduce small-scale power spectra of the $N$-body simulations, while keeping the precision at large scales. The revised fitting formula provides an accurate prediction of the nonlinear matter power spectrum in a wide range of wavenumber ($k \leq 30h$\,Mpc$^{-1}$) at redshifts $0 \leq z \leq 10$, with 5% precision for $k\leq1 h$ Mpc$^{-1}$ at $0 \leq z \leq 10$ and 10% for $1 \leq k\leq 10 h$ Mpc$^{-1} $ at $0 \leq z \leq 3$. We discuss the impact of the improved halofit model on weak lensing power spectra and correlation functions, and show that the improved model better reproduces ray-tracing simulation results.

We examine scatter and bias in weak lensing selected clusters, employing both an analytic model of dark matter haloes and numerical mock data of weak lensing cluster surveys. We pay special attention to effects of the diversity of dark matter distributions within clusters. We find that peak heights of the lensing convergence map correlates rather poorly with the virial mass of haloes. The correlation is tighter for the spherical overdensity mass with a higher mean interior density. We examine the dependence of the halo shape on the peak heights, and find that the rms scatter caused by the halo diversity scales linearly with the peak heights with the proportionality factor of 0.1-0.2. The noise originated from the halo shape is found to be comparable to the source galaxy shape noise and the cosmic shear noise. We find the significant halo orientation bias, i.e., weak lensing selected clusters on average have their major axes aligned with the line-of-sight direction. We compute the orientation bias using an analytic triaxial halo model and obtain results quite consistent with the ray-tracing results. We develop a prescription to analytically compute the number count of weak lensing peaks taking into account all the main sources of scatters in peak heights. We find that the improved analytic predictions agree well with the simulation results for high S/N peaks. We also compare the expected number count with our weak lensing analysis results for 4 sq deg of Subaru/Suprime-Cam observations and find a good agreement.

We investigate the performance of an analytic model of the 3D matter distribution, which combines perturbation theory with halo models, for weak-lensing statistics. We compare our predictions for the weak-lensing convergence power spectrum and bispectrum with numerical simulations and fitting formulas proposed in previous works. We find that this model provides better agreement with simulations than published fitting formulas. This shows that building on systematic and physically motivated models is a promising approach. Moreover, this makes explicit the link between the weak-lensing statistics and the underlying properties of the 3D matter distribution, as a function of scale $\ell$. Thus, we obtain the contributions to the lensing power spectrum and bispectrum that arise from perturbative terms (complete up to one-loop) and nonperturbative terms (e.g., "1-halo" term). Finally, we show that this approach recovers the dependence on cosmology (for realistic scenarios).

We investigate the performance of an analytic model of the 3D matter distribution, which combines perturbation theory with halo models, for weak-lensing configuration-space statistics. We compared our predictions for the weak-lensing convergence two-point and three-point correlation functions with numerical simulations and fitting formulas proposed in previous works. We also considered the second- and third-order moments of the smoothed convergence and of the aperture-mass. As in our previous study of Fourier-space weak-lensing statistics, we find that our model agrees better with simulations than previously published fitting formulas. Moreover, we recover the dependence on cosmology of these weak-lensing statistics and we can describe multi-scale moments. This approach allows us to obtain the quantitative relationship between these integrated weak-lensing statistics and the various contributions to the underlying 3D density fluctuations, decomposed over perturbative, two-halo, or one-halo terms.

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

We study the lensing convergence power spectrum and its covariance for a standard LCDM cosmology. We run 400 cosmological N-body simulations and use the outputs to perform a total of 1000 independent ray-tracing simulations. We compare the simulation results with analytic model predictions. The semi-analytic model based on Smith et al.(2003) fitting formula underestimates the convergence power by ~30% at arc-minute angular scales. For the convergence power spectrum covariance, the halo model reproduces the simulation results remarkably well over a wide range of angular scales and source redshifts. The dominant contribution at small angular scales comes from the sample variance due to the number fluctuations of halos in a finite survey volume. The signal-to-noise ratio for the convergence power spectrum is degraded by the non-Gaussian covariances by up to a factor 5 for a weak lensing survey to z_s ~1. The probability distribution of the convergence power spectrum estimators, among the realizations, is well approximated by a chi-square distribution with broadened variance given by the non-Gaussian covariance, but has a larger positive tail. The skewness and kurtosis have non-negligible values especially for a shallow survey. We argue that a prior knowledge on the full distribution may be needed to obtain an unbiased estimate on the ensemble averaged band power at each angular scale from a finite volume survey.

A logarithmic transform of the convergence field improves `the information content', ie., the overall precision associated with the measurement of the amplitude of the convergence power spectrum by improving the covariance matrix properties. The translation of this improvement in the information content to that in cosmological parameters, such as those associated with dark energy, requires knowing the sensitivity of the log-transformed field to those cosmological parameters. In this paper we use N-body simulations with ray tracing to generate convergence fields at multiple source redshifts as a function of cosmology. The gain in information associated with the log-transformed field does lead to tighter constraints on dark energy parameters, but only if shape noise is neglected. The presence of shape noise quickly diminishes the advantage of the log mapping, more quickly than we would expect based on the information content. With or without shape noise, using a larger pixel size allows for a more efficient log-transformation.

We perform high resolution ray-tracing simulations to investigate probability distribution functions (PDFs) of lensing convergence, shear, and magnification on distant sources up to the redshift of $z_s=20$. We pay particular attention to the shot noise effect in $N$-body simulations by explicitly showing how it affects the variance of the convergence. We show that the convergence and magnification PDFs are closely related with each other via the approximate relation $\mu=(1-\kappa)^{-2}$, which can reproduce the behavior of PDFs surprisingly well up to the high magnification tail. The mean convergence measured in the source plane is found to be systematically negative, rather than zero as often assumed, and is correlated with the convergence variance. We provide simple analytical formulae for the PDFs, which reproduce simulated PDFs reasonably well for a wide range of redshifts and smoothing sizes. As explicit applications of our ray-tracing simulations, we examine the strong lensing probability and the magnification effects on the luminosity functions of distant galaxies and quasars.

Understanding a behavior of galaxy biasing is crucial for future galaxy redshift surveys. One aim is to measure the baryon acoustic oscillations (BAOs) within the precision of a few percent level. Using 30 large cosmological N-body simulations for a standard LCDM cosmology, we study the halo biasing over a wide redshift range. We compare the simulation results with theoretical predictions proposed by Matsubara (2008) which naturally incorporate the halo bias and redshift-space distortions into their formalism of perturbation theory with a resummation technique via the Lagrangian picture. The power spectrum and correlation function of halos obtained from Lagrangian resummation theory (LRT) well agree with N-body simulation results on scales of BAOs. Especially nonlinear effects on the baryon acoustic peak of the halo correlation function are accurately explained both in real and redshift space. We find that nonlinearity and scale dependence of bias are fairly well reproduced by 1-loop LRT up to k=0.35hMpc^{-1} (z=2 and 3) within a few percent level in real space and up to k=0.1hMpc^{-1} (z=2) and 0.15hMpc^{-1} (z=3) in redshift space. Thus, the LRT is very powerful for accurately extracting cosmological information in upcoming high redshift BAO surveys.

Using 1000 ray-tracing simulations for a {\Lambda}-dominated cold dark model in Sato et al. (2009), we study the covariance matrix of cosmic shear correlation functions, which is the standard statistics used in the previous measurements. The shear correlation function of a particular separation angle is affected by Fourier modes over a wide range of multipoles, even beyond a survey area, which complicates the analysis of the covariance matrix. To overcome such obstacles we first construct Gaussian shear simulations from the 1000 realizations, and then use the Gaussian simulations to disentangle the Gaussian covariance contribution to the covariance matrix we measured from the original simulations. We found that an analytical formula of Gaussian covariance overestimates the covariance amplitudes due to an effect of finite survey area. Furthermore, the clean separation of the Gaussian covariance allows to examine the non-Gaussian covariance contributions as a function of separation angles and source redshifts. For upcoming surveys with typical source redshifts of z_s=0.6 and 1.0, the non-Gaussian contribution to the diagonal covariance components at 1 arcminute scales is greater than the Gaussian contribution by a factor of 20 and 10, respectively. Predictions based on the halo model qualitatively well reproduce the simulation results, however show a sizable disagreement in the covariance amplitudes. By combining these simulation results we develop a fitting formula to the covariance matrix for a survey with arbitrary area coverage, taking into account effects of the finiteness of survey area on the Gaussian covariance.

Astrometric Very Long Baseline Interferometry (VLBI) observations of maser sources in the Milky Way are used to map the spiral structure of our Galaxy and to determine fundamental parameters such as the rotation velocity ($\Theta_0$) and curve and the distance to the Galactic center (R$_0$). Here, we present an update on our first results, implementing a recent change in the knowledge about the Solar motion. It seems unavoidable that the IAU recommended values for R$_0$ and $\Theta_0$ need a substantial revision. In particular the combination of 8.5 kpc and 220 \kms\, can be ruled out with high confidence. Combining the maser data with the distance to the Galactic center from stellar orbits and the proper motion of Sgr\,A* gives best values of R$_0$ = 8.3 $\pm$ 0.23 kpc and $\Theta_0$ = 239 or 246 $\pm$ 7 \kms, for Solar motions of V$_ \odot$ = 12.23 and 5.25 \kms, respectively. Finally, we give an outlook to future observations in the Bar and Spiral Structure Legacy (BeSSeL) Survey.

To estimate cosmological parameters from a given dataset, we need to construct a likelihood function, which sometimes has a complicated functional form. We introduce the copula, a mathematical tool to construct an arbitrary multivariate distribution function from one-dimensional marginal distribution functions with any given dependence structure. It is shown that a likelihood function constructed by the so-called Gaussian copula can reproduce very well the n-dimensional probability distribution of the cosmic shear power spectrum obtained from a large number of ray-tracing simulations. This suggests that the Copula likelihood will be a powerful tool for future weak lensing analyses, instead of the conventional multivariate Gaussian likelihood.

The estimation of cosmological parameters from a given data set requires a construction of a likelihood function which, in general, has a complicated functional form. We adopt a Gaussian copula and constructed a copula likelihood function for the convergence power spectrum from a weak lensing survey. We show that the parameter estimation based on the Gaussian likelihood erroneously introduces a systematic shift in the confidence region, in particular for a parameter of the dark energy equation of state w. Thus, the copula likelihood should be used in future cosmological observations.

We report measurement of the trigonometric parallax of W51 Main/South using the Very Long Baseline Array (VLBA). We measure a value of 0.185 +/- 0.010 mas, corresponding to a distance of 5.41 (+0.31/-0.28) kpc. W51 Main/South is a well-known massive star-forming region near the tangent point of the Sagittarius spiral arm of the Milky Way. Our distance to W51 yields an estimate of the distance to the Galactic center of Ro = 8.3 +/- 0.46 (statistical) +/- 1.0 (systematic) kpc by simple geometry. Combining the parallax and proper motion measurements for W51, we obtained the full-space motion of this massive star forming region. We find W51 is in a nearly circular orbit about the Galactic center. The H2O masers used for our parallax measurements trace four powerful bipolar outflows within a 0.4 pc size region, some of which are associated with dusty molecular hot cores and/or hyper- or ultra-compact HII regions.

Gravitational lensing of distant galaxies can be exploited to infer the convergence field as a function of angular position on the sky. The statistics of this field, much like that of the cosmic microwave background (CMB), can be studied to extract information about fundamental parameters in cosmology, most notably the dark energy in the Universe. Unlike the CMB, the distribution of matter in the Universe which determines the convergence field is highly non-Gaussian, reflecting the nonlinear processes which accompanied structure formation. Much of the cosmic information contained in the initial field is therefore unavailable to the standard power spectrum measurements. Here we propose a method for re-capturing cosmic information by using the power spectrum of a simple function of the observed (nonlinear) convergence field. We adapt the approach of Neyrinck et al. (2009) to lensing by using a modified logarithmic transform of the convergence field. The Fourier transform of the log-transformed field has modes that are nearly uncorrelated, which allows for additional cosmological information to be extracted from small-scale modes.

Continuing work initiated in an earlier publication [Sato and Asada, PASJ, 61, L29 (2009)], we consider light curves influenced by the orbital inclination and eccentricity of a companion in orbit around a transiting extrasolar planet (in a planet-satellite system or a hypothetical true binary). We show that the semimajor axis, eccentricity and inclination angle of a `moon' orbit around the host planet can be determined by transit method alone. For this purpose, we present a formulation for the parameter determinations in a small-eccentricity approximation as well as in the exact form. As a result, the semimajor axis is expressed in terms of observables such as brightness changes, transit durations and intervals in light curves. We discuss also a narrow region of parameters that produce a mutual transit by an extrasolar satellite.

We report on trigonometric parallax measurements for the Galactic star forming region G14.33-0.64 toward the Sagittarius spiral arm. We conducted multi-epoch phase-referencing observations of an H2O maser source in G14.33-0.64 with the Japanese VLBI array VERA. We successfully detected a parallax of 0.893+/-0.101 mas, corresponding to a source distance of 1.12+/-0.13 kpc, which is less than half of the kinematic distance for G14.33-0.64. Our new distance measurement demonstrates that the Sagittarius arm lies at a closer distance of ~1 kpc, instead of previously assumed ~2-3 kpc from kinematic distances. The previously suggested deviation of the Sagittarius arm toward the Galactic center from the symmetrically fitted model (Taylor & Cordes 1993) is likely due to large errors of kinematic distances at low galactic longitudes. G14.33-0.64 most likely traces the near side of the Sagittarius arm. We attempted fitting the pitch angle of the arm with other parallax measurements along the arm, which yielded two possible pitch angles of i=34.7+/-2.7 degrees and i=11.2+/-10.5 degrees. Our proper motion measurements suggest G14.33-0.64 has no significant peculiar motion relative to the differential rotation of the Galaxy (assumed to be in a circular orbit), indicating that the source motion is in good agreement with the Galactic rotation.

We consider the effects of mutual transits by extrasolar planet-companion systems (in a true binary or a planet-satellite system) on light curves. We show that induced changes in light curves depend strongly on a ratio between a planet-companion's orbital velocity around their host star and a planet-companion's spin speed around their common center of mass. In both the slow and fast spin cases (corresponding to long and short distances between them, respectively), a certain asymmetry appears in light curves. We show that, especially in the case of short distances, occultation of one faint object by the other, while the transit of the planet-companion system occurs in front of its parent star, causes an apparent increase in light curves and characteristic fluctuations appear as important evidence of mutual transits. We show also that extrasolar mutual transits provide a complementary method of measuring the radii of two transiting objects, their separation and mass, and consequently identifying them as a true binary, planet-satellite system or others. Monitoring $10^5$ stars for three years with Kepler may lead to a discovery of a second Earth-Moon-like system if the fraction of such systems for an averaged star is larger than 0.05, or it may put upper limits on the fraction as f < 0.05.

We report Suzaku observations of the galaxy cluster Abell 1795 that extend to r_200 ~ 2 Mpc, the radius within which the mean cluster mass density is 200 times the cosmic critical density. These observations are the first to probe the state of the intracluster medium in this object at r > 1.3 Mpc. We sample two disjoint sectors in the cluster outskirts (1.3 < r < 1.9 Mpc) and detect X-ray emission in only one of them to a limiting (3-sigma) soft X-ray surface brightness of B(0.5-2 keV) = 1.8 x 10^-12 erg s^-1 cm^-2 deg^-2, a level less than 20% of the cosmic X-ray background brightness. We trace the run of temperature with radius at r > 0.4 Mpc and find that it falls relatively rapidly (T ~ r^-0.9), reaching a value about one third of its peak at the largest radius we can measure it. Assuming the intracluster medium is in hydrostatic equilibrium and is polytropic, we find a polytropic index of 1.3 +0.3-0.2 and we estimate a mass of 4.1 +0.5-0.3 x 10^14 M_solar within 1.3 Mpc, somewhat (2.7-sigma) lower than that reported by previous observers. However, our observations provide evidence for departure from hydrostatic equilibrium at radii as small as r ~ 1.3 Mpc ~ r_500 in this apparently regular and symmetrical cluster.

We are using the VLBA and the Japanese VERA project to measure trigonometric parallaxes and proper motions of masers found in high-mass star-forming regions across the Milky Way. Early results from 18 sources locate several spiral arms. The Perseus spiral arm has a pitch angle of 16 +/- 3 degrees, which favors four rather than two spiral arms for the Galaxy. Combining positions, distances, proper motions, and radial velocities yields complete 3-dimensional kinematic information. We find that star forming regions on average are orbiting the Galaxy ~15 km/s slower than expected for circular orbits. By fitting the measurements to a model of the Galaxy, we estimate the distance to the Galactic center R_o = 8.4 +/- 0.6 kpc and a circular rotation speed Theta_o = 254 +/- 16 km/s. The ratio Theta_o/R_o can be determined to higher accuracy than either parameter individually, and we find it to be 30.3 +/- 0.9 km/s/kpc, in good agreement with the angular rotation rate determined from the proper motion of Sgr A*. The data favor a rotation curve for the Galaxy that is nearly flat or slightly rising with Galactocentric distance. Kinematic distances are generally too large, sometimes by factors greater than two; they can be brought into better agreement with the trigonometric parallaxes by increasing Theta_o/R_o from the IAU recommended value of 25.9 km/s/kpc to a value near 30 km/s/kpc. We offer a "revised" prescription for calculating kinematic distances and their uncertainties, as well as a new approach for defining Galactic coordinates. Finally, our estimates of Theta_o and To/R_o, when coupled with direct estimates of R_o, provide evidence that the rotation curve of the Milky Way is similar to that of the Andromeda galaxy, suggesting that the dark matter halos of these two dominant Local Group galaxy are comparably massive.

We present a distance measurement for the semiregular variable S Crateris (S Crt) based on its annual parallax. With the unique dual beam system of the VLBI Exploration for Radio Astrometry (VERA) telescopes, we measured the absolute proper motion of a water maser spot associated with S Crt, referred to the quasar J1147-0724 located at an angular separation of 1.23$^{\circ}$. In observations spanning nearly two years, we have detected the maser spot at the LSR velocity of 34.7 km s$^{-1}$, for which we measured the annual parallax of 2.33$\pm$0.13 mas corresponding to a distance of 430$^{+25}_{-23}$ pc. This measurement has an accuracy one order of magnitude better than the parallax measurements of HIPPARCOS. The angular distribution and three-dimensional velocity field of maser spots indicate a bipolar outflow with the flow axis along northeast-southwest direction. Using the distance and photospheric temperature, we estimate the stellar radius of S Crt and compare it with those of Mira variables.

We report astrometric observations of H2O masers around the red supergiant VY Canis Majoris (VY CMa) carried out with VLBI Exploration of Radio Astrometry (VERA). Based on astrometric monitoring for 13 months, we successfully measured a trigonometric parallax of 0.88 +/- 0.08 mas, corresponding to a distance of 1.14 +0.11/-0.09 kpc. This is the most accurate distance to VY CMa and the first one based on an annual parallax measurement. The luminosity of VY CMa has been overestimated due to a previously accepted distance. With our result, we re-estimate the luminosity of VY CMa to be (3 +/- 0.5) x 10^5 L_sun using the bolometric flux integrated over optical and IR wavelengths. This improved luminosity value makes location of VY CMa on the Hertzsprung-Russel (HR) diagram much closer to the theoretically allowable zone (i.e. the left side of the Hayashi track) than previous ones, though uncertainty in the effective temperature of the stellar surface still does not permit us to make a final conclusion.

We present results of multi-epoch VLBI observations with VERA (VLBI Exploration of Radio Astrometry) of the 22 GHz H$_{2}$O masers associated with a young stellar object (YSO) IRAS 22198+6336 in a dark cloud L1204G. Based on the phase-referencing VLBI astrometry, we derive an annual parallax of IRAS 22198+6336 to be 1.309$\pm$0.047 mas, corresponding to the distance of 764$\pm$27 pc from the Sun. Although the most principal error source of our astrometry is attributed to the internal structure of the maser spots, we successfully reduce the errors in the derived annual parallax by employing the position measurements for all of the 26 detected maser spots. Based on this result, we reanalyze the spectral energy distribution (SED) of IRAS 22198+6336 and find that the bolometric luminosity and total mass of IRAS 22198+6336 are 450$L_{\odot}$ and 7$M_{\odot}$, respectively. These values are consistent with an intermediate-mass YSO deeply embedded in the dense dust core, which has been proposed to be an intermediate-mass counterpart of a low-mass Class 0 source. In addition, we obtain absolute proper motions of the H$_{2}$O masers for the most blue-shifted components. We propose that the collimated jets aligned along the east-west direction are the most plausible explanation for the origin of the detected maser features.

We present 1.4 and 2.5 GHz Australia Telescope Compact Array (ATCA) observations of the galaxy cluster A3158 (z=0.0597) which is located within the central part of the Horologium-Reticulum Supercluster (HRS). Spectroscopic data for the central part of the HRS suggest that A3158 is in a dynamically important position within the supercluster and that it is moving toward the double cluster system A3125/A3128 which marks the centre of the HRS. A total of 110 radio galaxies are detected in a 35 arcminute radius about the cluster at 1.4 GHz, of which 30 are also detected at 2.5 GHz.We examine the source counts and compute the Radio Luminosity Function (RLF) at 1.4 GHz from the subset of 88 sources found within the fullwidth half-power area of the ATCA beam. Comparison of the source counts in the area over the background, as computed by Prandoni et al. (2001), shows some evidence of an excess of galaxies with L(1.4) GHz < 2 x 10^22 W/Hz. This result seems to indicate a star forming population and is a result similar to that found recently by Owen et al. (2005) for the merging cluster A2125. In addition we find that the radio luminosity function for early-type galaxies (E and S0) below log P(1.4) = 22.5 is lower than that found for a composite cluster environment (Ledlow & Owen, 1996) but is similar to the early-type RLF for clusters in the centre of the Shapley Super cluster (Venturi et al. 2000) which are believed to be in the latter stages of merging. This result implies that the cores of superclusters are environments where radio emission, particularly resultant from AGN, is suppressed in the later stages of merging. Thus, radio observations of clusters might be sensitive indicators of the precise merger stage of the cluster but more observational evidence is still required to establish this trend.

We have used the Japanese VLBI array VERA to perform high-precision astrometry of an H2O maser source in the Galactic star-forming region NGC 281 West, which has been considered to be part of a 300-pc superbubble. We successfully detected a trigonometric parallax of 0.355+/-0.030 mas, corresponding to a source distance of 2.82+/-0.24 kpc. Our direct distance determination of NGC 281 has resolved the large distance discrepancy between previous photometric and kinematic studies; likely NGC 281 is in the far side of the Perseus spiral arm. The source distance as well as the absolute proper motions were used to demonstrate the 3D structure and expansion of the NGC 281 superbubble, ~650 pc in size parallel to the Galactic disk and with a shape slightly elongated along the disk or spherical, but not vertically elongated, indicating the superbubble expansion may be confined to the disk. We estimate the expansion velocity of the superbubble as ~20 km/s both perpendicular to and parallel to the Galactic disk with a consistent timescale of ~20 Myr.

We present the initial results of multi-epoch VLBI observations of the 22 GHz H2O masers in the Orion KL region with VERA (VLBI Exploration of Radio Astrometry). With the VERA dual-beam receiving system, we have carried out phase-referencing VLBI astrometry and successfully detected an annual parallax of Orion KL to be 2.29+/-0.10 mas, corresponding to the distance of 437+/-19 pc from the Sun. The distance to Orion KL is determined for the first time with the annual parallax method in these observations. Although this value is consistent with that of the previously reported, 480+/-80 pc, which is estimated from the statistical parallax method using proper motions and radial velocities of the H2O maser features, our new results provide the much more accurate value with an uncertainty of only 4%. In addition to the annual parallax, we have detected an absolute proper motion of the maser feature, suggesting an outflow motion powered by the radio source I along with the systematic motion of source I itself.

We report on absolute proper-motion measurements of an H2O maser source in the NGC 281 West molecular cloud, which is located ~320 pc above the Galactic plane and is associated with an HI loop extending from the Galactic plane. We have conducted multi-epoch phase-referencing observations of the maser source with VERA (VLBI Exploration of Radio Astrometry) over a monitoring period of 6 months since May 2006. We find that the H2O maser features in NGC 281 West are systematically moving toward the southwest and further away from the Galactic plane with a vertical velocity of ~20-30 km/s at its estimated distance of 2.2-3.5 kpc. Our new results provide the most direct evidence that the gas in the NGC 281 region on the HI loop was blown out from the Galactic plane, most likely in a superbubble driven by multiple or sequential supernova explosions in the Galactic plane.

The in-orbit performance and calibration of the Hard X-ray Detector (HXD) on board the X-ray astronomy satellite Suzaku are described. Its basic performances, including a wide energy bandpass of 10-600 keV, energy resolutions of ~4 keV (FWHM) at 40 keV and ~11% at 511 keV, and a high background rejection efficiency, have been confirmed by extensive in-orbit calibrations. The long-term gains of PIN-Si diodes have been stable within 1% for half a year, and those of scintillators have decreased by 5-20%. The residual non-X-ray background of the HXD is the lowest among past non-imaging hard X-ray instruments in energy ranges of 15-70 and 150-500 keV. We provide accurate calibrations of energy responses, angular responses, timing accuracy of the HXD, and relative normalizations to the X-ray CCD cameras using multiple observations of the Crab Nebula.

The relation between shapes of cosmic string and baryon number violation is investigated. If there exist fermionic zero modes on the string, using bosonization technique, it is possible to obtain the effective action which describes the fermion coupled to the arbitrarily shaped cosmic string. The relation between baryon number and the sum of the writhing number and linking number of the cosmic strings is rederived. Baryons are created on the strings as the shapes of the cosmic strings change. Furthermore we discuss implications of this baryon number violating process to baryogenesis.

We present a ``topological'' formulation of arbitrarily shaped vortex strings in four dimensional field theory. By using a large Higgs mass expansion, we then evaluate the effective action of the closed Abrikosov-Nielsen-Olesen vortex string. It is shown that the effective action contains the Nambu-Goto term and an extrinsic curvature squared term with negative sign. We next evaluate the topological $\FtF$ term and find that it becomes the sum of an ordinary self-intersection number and Polyakov's self-intersection number of the world sheet swept by the vortex string. These self-intersection numbers are related to the self-linking number and the total twist number, respectively. Furthermore, the $\FtF$ term turns out to be the difference between the sum of the writhing numbers and the linking numbers of the vortex strings at the initial time and the one at the final time. When the vortex string is coupled to fermions, the chiral fermion number of the vortex string becomes the writhing number (modulo $\bZ$) through the chiral anomaly. Our formulation is also applied to ``global'' vortex strings in a model with a broken global $U(1)$ symmetry.