We report on high angular resolution polarimetric observations of the nearby radio galaxy M87 using the Very Long Baseline Array at 24 GHz ($\lambda=$1.3 cm) and 43 GHz ($\lambda=$7 mm) in 2017-2018. New images of the linear polarization substructure in the nuclear region are presented, characterized by a two-component pattern of polarized intensity and smooth rotation of the polarization plane around the 43 GHz core. From a comparison with an analogous dataset from 2007, we find that this global polarization pattern remains stable on a time interval of 11 yr, while showing smaller month-scale variability. We discuss the possible Faraday rotation toward the M87 nucleus at centimeter to millimeter wavelengths. These results can be interpreted in a scenario where the observed polarimetric pattern is associated with the magnetic structure in the confining magnetohydrodynamic wind, which also serves as the source of the observed Faraday rotation.

Planetary-scale waves are thought to play a role in powering the yet-unexplained atmospheric superrotation of Venus. Puzzlingly, while Kelvin, Rossby and stationary waves manifest at the upper clouds (65--70 km), no planetary-scale waves or stationary patterns have been reported in the intervening level of the lower clouds (48--55 km), although the latter are probably Lee waves. Using observations by the Akatsuki orbiter and ground-based telescopes, we show that the lower clouds follow a regular cycle punctuated between 30$^{\circ}$N--40$^{\circ}$S by a sharp discontinuity or disruption with potential implications to Venus's general circulation and thermal structure. This disruption exhibits a westward rotation period of $\sim$4.9 days faster than winds at this level ($\sim$6-day period), alters clouds' properties and aerosols, and remains coherent during weeks. Past observations reveal its recurrent nature since at least 1983, and numerical simulations show that a nonlinear Kelvin wave reproduces many of its properties.

Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

We study the kinematics of the M87 jet using the first year data of the KVN and VERA Array (KaVA) large program, which has densely monitored the jet at 22 and 43 GHz since 2016. We find that the apparent jet speeds generally increase from $\approx0.3c$ at $\approx0.5$ mas from the jet base to $\approx2.7c$ at $\approx20$ mas, indicating that the jet is accelerated from subluminal to superluminal speeds on these scales. We perform a complementary jet kinematic analysis by using archival Very Long Baseline Array monitoring data observed in $2005-2009$ at 1.7 GHz and find that the jet is moving at relativistic speeds up to $\approx5.8c$ at distances of $200-410$ mas. We combine the two kinematic results and find that the jet is gradually accelerated over a broad distance range that coincides with the jet collimation zone, implying that conversion of Poynting flux to kinetic energy flux takes place. If the jet emission consists of a single streamline, the observed trend of jet acceleration ($\Gamma\propto z^{0.16\pm0.01}$) is relatively slow compared to models of a highly magnetized jet. This indicates that Poynting flux conversion through the differential collimation of poloidal magnetic fields may be less efficient than expected. However, we find a non-negligible dispersion in the observed speeds for a given jet distance, making it difficult to describe the jet velocity field with a single power-law acceleration function. We discuss the possibility that the jet emission consists of multiple streamlines following different acceleration profiles, resulting in jet velocity stratification.

We investigated the jet width profile with distance along the jet in the nearby radio galaxy NGC 1052 at radial distances between $\sim300$ to $4 \times 10^7$ Schwarzschild Radii($R_{\rm S}$) from the central engine on both their approaching and receding jet sides. The width of jets was measured in images obtained with the VLBI Space Observatory Programme (VSOP), the Very Long Baseline Array (VLBA), and the Very Large Array (VLA). The jet-width profile of receding jets are apparently consistent with that of approaching jets throughout the measuring distance ranges, indicating symmetry at least up to the sphere of gravitational influence of the central black hole. The power-law index $a$ of the jet-width profile ($w_{\rm{jet}} \propto r^{a}$, where $w_{\rm jet}$ is the jet width, $r$ is the distance from the central engine in the unit of $R_{\rm S}$) apparently shows a transition from $a \sim 0$ to $a \sim 1$, i.e., the cylindrical-to-conical jet structures, at a distance of $\sim1\times10^{4} \ R_{\mathrm{S}}$. The cylindrical jet shape at the small distances is reminiscent of the innermost jets in 3C 84. Both the central engines of NGC 1052 and 3C 84 are surrounded by dense material, part of which is ionized and causes heavy free-free absorption.

Planetary-scale waves at the Venusian cloud-top cause periodic variations in both winds and ultraviolet (UV) brightness. While the wave candidates are the 4-day Kelvin wave and 5-day Rossby wave with zonal wavenumber 1, their temporal evolutions are poorly understood. Here we conducted a time series analysis of the 365-nm brightness and cloud-tracking wind variations, obtained by the UV Imager onboard the Japanese Venus Climate Orbiter Akatsuki from June to October 2017, revealing a dramatic evolution of planetary-scale waves and corresponding changes in planetary-scale UV features. We identified a prominent 5-day periodicity in both the winds and brightness variations, whose phase velocities were slower than the dayside mean zonal winds (or the super-rotation) by >35 m s$^{-1}$. The reconstructed planetary-scale vortices were nearly equatorially symmetric and centered at ~35{\deg} latitude in both hemispheres, which indicated that they were part of a Rossby wave. The amplitude of winds variation associated with the observed Rossby wave packet were amplified gradually over ~20 days and attenuated over ~50 days. Following the formation of the Rossby wave vortices, brightness variation emerges to form rippling white cloud belts in the 45{\deg}-60{\deg} latitudes of both hemispheres. ~3.8-day periodic signals were observed in the zonal wind and brightness variations in the equatorial region before the Rossby wave amplification. Although the amplitude and significance of the 3.8-day mode were relatively low in the observation season, this feature is consistent with a Kelvin wave, which may be the cause of the dark clusters in the equatorial region.

Seasonal cycle of Martian relative humidity shows that, in the southern and the northern RSL sites, humidity decreases when RSL increment while humidity is high on the fading season. We suggest that the recurrence of RSL is caused by the different angle of response between wet and dry grains. 2D DEM simulation indicates that, when the ratio of a liquid bridge volume to a particle volume is the order of 10^{-7}, the repose angle of Martian sand can increase to ~40 degrees due to cohesion, which is near the slope around RSL starting area. Thus, if this amount of liquid bridge is generated by deliquescence or capillary condensation of vapor and if the water evaporates following the humidity and temperature cycles, Martian sand can induce avalanche on RSL slope. When grains are wet again, sand and dust deposit on RSL, which may be the mechanism of RSL disappearance and recharge of grains for the next season. In our model, because liquid induces RSL indirectly (RSL itself is the granular flow), the amount of liquid can be drastically reduced compared with the wet mechanism. In this case, as well as temperature, the seasonal cycle of relative humidity plays an important role in the recurrence of RSL.

Longwave Infrared Camera (LIR) onboard Akatsuki first revealed the global structure of the thermal tides in the upper cloud layer of Venus. The data were acquired over three Venusian years, and the analysis was done over the areas from the equator to the mid-latitudes in both hemispheres and over the whole local time. Thermal tides at two vertical levels were analyzed by comparing data at two different emission angles. Dynamical wave modes consisting of tides were identified; the diurnal tide consisted mainly of Rossby-wave and gravity-wave modes, while the semidiurnal tide predominantly consisted of a gravity-wave mode. The revealed vertical structures were roughly consistent with the above wave modes, but some discrepancy remained if the waves were supposed to be monochromatic. In turn, the heating profile that excites the tidal waves can be constrained to match this discrepancy, which would greatly advance the understanding of the Venusian atmosphere.

We propose a new method of estimating a mass of a super massive black hole residing in the center of an active galaxy. The active galaxy M87 offers a convenient test case for the method due to the existence of a large amount of observational data on the jet and ambient environment properties in the central area of the object. We suggest that the observed transition of a jet boundary shape from a parabolic to a conical form is associated with the flow transiting from the magnetically dominated regime to the energy equipartition between plasma bulk motion and magnetic field. By coupling the unique set of observations available for the jet kinematics, environment and boundary profile with our MHD modelling under assumption on the presence of a dynamically important magnetic field in the M87 jet, we estimate the central black hole mass and spin. The method leads us to believe that the M87 super massive black hole has a mass somewhat larger than typically accepted so far.

We present ALMA observations at 233, 678, and 870 GHz of the Galactic Center black hole, Sagittarius A*. These observations reveal a flat spectrum over this frequency range with spectral index $\alpha \approx -0.3$, where the flux density $S \propto \nu^\alpha$. We model the submm and far infrared spectrum with a one zone synchrotron model of thermal electrons. We infer electron densities $n = (2-5) \times 10^6$ cm$^{-3}$, electron temperatures $T_e = (1-3) \times 10^{11}$ K, and magnetic field strength $B = 10-50$ G. The parameter range can be further constrained using the observed quiescent X-ray luminosity. The flat submm spectrum results in a high electron temperature and implies that the emitting electrons are efficiently heated. We also find that the emission is most likely optically thin at 233 GHz. These results indicate that millimeter and submillimeter wavelength very long baseline interferometry of Sgr A* including those of the Event Horizon Telescope should see a transparent emission region down to event horizon scales.

In the current paradigm, it is believed that the compact VLBI radio core of radio-loud active galactic nuclei (AGNi) represents the innermost upstream regions of relativistic outflows. These regions of AGN jets have generally been modeled by a conical outflow with a roughly constant opening angle and flow speed. Nonetheless, some works suggest that a parabolic geometry would be more appropriate to fit the high energy spectral distribution properties and it has been recently found that, at least in some nearby radio galaxies, the geometry of the innermost regions of the jet is parabolic. We compile here multi-frequency core sizes of archival data to investigate the typically unresolved upstream regions of the jet geometry of a sample of 56 radio-loud AGNi. Data combined from the sources considered here are not consistent with the classic picture of a conical jet starting in the vicinity of the super-massive black hole (SMBH), and may exclude a pure parabolic outflow solution, but rather suggest an intermediate solution with quasi-parabolic streams, which are frequently seen in numerical simulations. Inspection of the large opening angles near the SMBH and the range of the Lorentz factors derived from our results support our analyses. Our result suggests that the conical jet paradigm in AGNi needs to be re-examined by millimeter/sub-millimeter VLBI observations.

The narrow-line Seyfert 1 (NLS1) galaxy PMN J0948+0022 is an archetype of gamma-ray-emitting NLS1s in active galactic nuclei (AGNs). In this study, we investigate its radio structures using archival data obtained using the Karl G.~Jansky very large array (VLA) and the very long baseline array (VLBA). The new VLA images reveal an outermost radio emission separated by 9.1 arcsec. Its resolved structure and steep spectrum suggest a terminal shock in a radio lobe energized by the jet from the PMN J0948+0022 nucleus. This large-scale radio component is found at almost the same position angle as that of the pc-scale jet, indicating a stable jet direction up to ~1 Mpc. Its apparent one-sidedness implies a moderate advancing speed ($\beta>0.27$). The kinematic age of 1 x 10^7 year of the jet activity is consistent with the expected NLS1 phase of ~10^7 year in the AGN lifetime. The VLBA image reveals the jet structure at distances ranging from r=0.82 milliarcsec to 3.5 milliarcsec, corresponding to approximately 100 pc, where superluminal motions were found. The jet width profile ($\propto r^{1.12}$) and flux-density profile ($\propto r^{-1.44}$) depending on the distance $r$ along the jet suggest that the jet kinetic energy is converted into internal energy in this region. The jet is causally connected in a nearly conical streamline, which is supported by ambient pressure at 100 pc scales in the host galaxy of PMN J0948+0022.

Jets are a ubiquitous part of the accretion process, seen in a wide variety of objects ranging from active galaxies (AGN) to X-ray binary stars and even newly formed stars. AGN jets are accelerated by the supermassive black hole of their host galaxy by a coupling between the magnetic field and inflowing material. They are the source for many exciting phenomena and can profoundly influence the larger galaxy and surrounding cluster. This White Paper points out what advances can be achieved in the field by new technologies, concentrating on the zone where jets are accelerated to relativistic speeds and collimated. The ngVLA and new space VLBI missions will give higher angular resolution, sensitivity and fidelity in the radio, penetrating this zone for additional objects and allowing us to resolve fundamental questions over the physics of jet acceleration and collimation. Interferometry in other bands would allow us to probe directly flaring components. We also emphasize the need for polarimetry, which is essential to revealing the role and configuration of magnetic fields.

We study Faraday rotation in the jet of M87 inside the Bondi radius using eight Very Long Baseline Array data sets, one at 8 GHz, four at 5 GHz, and three at 2 GHz. We obtain Faraday rotation measures (RMs) measured across the bandwidth of each data set. We find that the magnitude of RM systematically decreases with increasing distance from the black hole from 5,000 to 200,000 Schwarzschild radii. The data, showing predominantly negative RM sign without significant difference of the RMs on the northern and southern jet edges, suggest that the spatial extent of the Faraday screen is much larger than the jet. We apply models of hot accretion flows, thought to be prevalent in active galactic nuclei having relatively low luminosity such as M87, and find that the decrease of RM is described well by a gas density profile $\rho \propto r^{-1}$. This behavior matches the theoretically expected signature of substantial winds, nonrelativistic un-collimated gas outflows from hot accretion flows, which is consistent with the results of various numerical simulations. The pressure profile inferred from the density profile is flat enough to collimate the jet, which can result in gradual acceleration of the jet in a magneto-hydrodynamical process. This picture is in good agreement with the observed gradual collimation and acceleration of the M87 jet inside the Bondi radius. The dominance of negative RMs suggests that jet and wind axis are misaligned such that the jet emission exposes only one side of the toroidal magnetic fields permeating the winds.

We present measurements of the wind speeds at the nightside lower clouds of Venus from observations by JAXA's mission Akatsuki during 2016, complemented with new wind measurements from ground-based observations acquired with TNG/NICS in 2012 and IRTF/SpeX in 2015 and 2017. Zonal and meridional components of the winds were measured from cloud tracking on a total of 466 Akatsuki images of Venus acquired by the camera IR2 using the 2.26-$\mathrm{\mu m}$ filter, with spatial resolutions ranging 10--80 km per pixel and covering from 2016 March 22 to October 31. More than 149,000 wind vectors were obtained with an automatic technique of template matching, and 2,947 wind vectors were inferred with the manual procedure. The meridional profiles for both components of the winds are found to be consistent with results from the Venus Express mission during 2006--2008, although stronger wind variability is found for the zonal component at equatorial latitudes where Akatsuki observations have better viewing geometry than Venus Express. The zonal winds at low latitudes also suggest a zonal variability that could be associated with solar tides or vertically propagating orographic waves. Finally, the combination of our wind measurements from TNG/NICS, IRTF/SpeX and Akatsuki images with previously published and based in data from 1978 to 2017 suggests variations of up to 30 m s$^{-1}$ in the winds at the lower clouds of the Venus nightside.

The M87 jet is extensively examined by utilizing general relativistic magnetohydrodynamic (GRMHD) simulations as well as the steady axisymmetric force-free electrodynamic (FFE) solution. Quasi-steady funnel jets are obtained in GRMHD simulations up to the scale of $\sim 100$ gravitational radius ($r_{\rm g}$) for various black hole (BH) spins. As is known, the funnel edge is approximately determined by the following equipartitions; i) the magnetic and rest-mass energy densities and ii) the gas and magnetic pressures. Our numerical results give an additional factor that they follow the outermost parabolic streamline of the FFE solution, which is anchored to the event horizon on the equatorial plane. We also identify the matter dominated, non-relativistic corona/wind play a dynamical role in shaping the funnel jet into the parabolic geometry. We confirm a quantitative overlap between the outermost parabolic streamline of the FFE jet and the edge of jet sheath in VLBI observations at $\sim 10^{1}$-$10^{5} \, r_{\rm g}$, suggesting that the M87 jet is likely powered by the spinning BH. Our GRMHD simulations also indicate a lateral stratification of the bulk acceleration (i.e., the spine-sheath structure) as well as an emergence of knotty superluminal features. The spin characterizes the location of the jet stagnation surface inside the funnel. We suggest that the limb-brightened feature could be associated with the nature of the BH-driven jet, if the Doppler beaming is a dominant factor. Our findings can be examined with (sub-)mm VLBI observations, giving a clue for the origin of the M87 jet.

This paper theoretically investigates the relations between the structure of relativistic jets and produced synchrotron images, by using a steady axisymmetric force-free jet model. We especially focus on the limb-brightened jets that are largely symmetric to the jet axes and observed in some active galactic nuclei such as M87, Mrk~501, Cyg~A, and 3C84. We find that symmetrically limb-brightened images can be produced when magnetic field lines of the jet penetrate a fast-spinning black hole as motivated by the Blandford-Znajek mechanism. On the other hand, jets with magnetic field lines that pass through a slowly spinning black hole or the Keplerian accretion disk produce highly asymmetric radio images. In addition, the edge of a counter jet tends to be luminous in the accretion-disk model even for rather small viewing angles, which may be problematic for some observed jets. We also suggest that the site of particle accelerations in relativistic jets can be constrained by fitting the radio images to observations. This kind of study focusing on the jet images far away from the central engine is complimentary to those concentrating directly on the innermost region with upcoming data of the Event Horizon Telescope.

We investigated the detailed radio structure of the jet of 1H 0323+342 using high-resolution multi-frequency Very Long Baseline Array observations. This source is known as the nearest $\gamma$-ray emitting radio-loud narrow-line Seyfert 1 (NLS1) galaxy. We discovered that the morphology of the inner jet is well characterized by a parabolic shape, indicating the jet being continuously collimated near the jet base. On the other hand, we found that the jet expands more rapidly at larger scales, resulting in a conical-like shape. The location of the "collimation break" is coincident with a bright quasi-stationary feature at 7 mas from core (corresponding to a deprojected distance of the order of $\sim$100pc), where the jet width locally contracts together with highly polarized signals, suggesting a recollimation shock. We found that the collimation region is coincident with the region where the jet speed gradually accelerates, suggesting the coexistence of the jet acceleration and collimation zone, ending up with the recollimation shock, which could be a potential site of high-energy $\gamma$-ray flares detected by the Fermi-LAT. Remarkably, these observational features of the 1H 0323+342 jet are overall very similar to those of the nearby radio galaxy M87 and HST-1 as well as some blazars, suggesting that a common jet formation mechanism might be at work. Based on the similarity of the jet profile between the two sources, we also briefly discuss the mass of the central black hole of 1H 0323+342, which is also still highly controversial on this source and NLS1s in general.

We promote the precise gamma-ray observation project Gamma-Ray Astro-Imager with Nuclear Emulsion (GRAINE), which uses balloon-borne emulsion gamma-ray telescopes. The emulsion telescope realizes observations with high angular resolution, polarization sensitivity, and large aperture area in the 0.01--100 GeV energy region. Herein, we report the data analysis of emulsion tracks and the first demonstration of gamma-ray imaging via an emulsion telescope by using the flight data from the balloon experiment performed in 2015 (GRAINE 2015). The emulsion films were scanned by the latest read-out system for a total area of 41 m$^2$ in three months, and then the gamma-ray event selection was automatically processed. Millions of electron-pair events are accumulated in the balloon-borne emulsion telescope. The emulsion telescope detected signals from a calibration source (gamma rays from the interaction of cosmic rays with an aluminum plate) with a high significance during the balloon observation and created a gamma-ray image consistent with the source size and the expected angular resolution in the energy range of 100--300 MeV. The flight performance obtained in the GRAINE 2015 experiment proves that balloon-borne emulsion telescope experiments with larger area are feasible while maintaining expected imaging performance.

Understanding the launching, acceleration, and collimation of jets powered by active galactic nuclei remains an outstanding problem in relativistic astrophysics. This is partly because observational tests of jet formation models suffer from the limited angular resolution of ground-based very long baseline interferometry that has thus far been able to probe the transverse jet structure in the acceleration and collimation zone of only two sources. Here we report radio interferometric observations of 3C 84 (NGC 1275), the central galaxy of the Perseus cluster, made with an array including the orbiting radio telescope of the RadioAstron mission. The obtained image transversely resolves the edge-brightened jet in 3C 84 only 30 microarcseconds from the core, which is ten times closer to the central engine than what has been possible in previous ground-based observations, and it allows us to measure the jet collimation profile from ~ 100 to ~10000 gravitational radii from the black hole. The previously found, almost cylindrical jet profile on scales larger than a few thousand r_g is now seen to continue at least down to a few hundred r_g from the black hole and we find a broad jet with a transverse radius larger than about 250 r_g at only 350 r_g from the core. If the bright outer jet layer is launched by the black hole ergosphere, it has to rapidly expand laterally on scales smaller than 100 r_g. If this is not the case, then this jet sheath is likely launched from the accretion disk.

Japanese Venus Climate Orbiter/AKATSUKI was proposed in 2001 with strong support by international Venus science community and approved as an ISAS (The Institute of Space and Astronautical Science) mission soon after the proposal. The mission life we expected was more than two Earth years in Venus orbit. AKATSUKI was successfully launched at 06:58:22JST on May 21, 2010, by H-IIA F17. After the separation from H-IIA, the telemetry from AKATSUKI was normally detected by DSN Goldstone station (10:00JST) and the solar cell paddles' deployment was confirmed. After a successful cruise, the malfunction happened on the propulsion system during the Venus orbit insertion (VOI) on Dec 7, 2010. The engine shut down before the planned reduction in speed to achieve. The spacecraft did not enter the Venus orbit, but entered an orbit around the Sun with a period of 203 days. Most of the fuel still had remained, but the orbital maneuvering engine was found to be broken and unusable. However, we have found an alternate way of achieving orbit by using only the reaction control system (RSC). We had adopted the alternate way for orbital maneuver and three minor maneuvers in Nov 2011 were successfully done so that AKATSUKI would meet Venus in 2015. We are considering several scenarios for VOI using only RCS.

We present Very Long Baseline Array polarimetric observations of the innermost jet of 3C$\sim$84 (NGC$\sim$1275) at 43$\sim$GHz. A significant polarized emission is detected at the hotspot of the innermost re-started jet, which is located $\sim$1 pc south from the radio core. While the previous report presented a hotspot at the southern end of the western limb, the hotspot location has been moved to the southern end of the eastern limb. Faraday rotation is detected within an entire bandwidth of the 43-GHz band. The measured rotation measure (RM) is at most (6.3$\pm$1.9)$\times10^{5}$$\sim$rad$\sim$m$^{-2}$ and might be slightly time variable on the timescale of a month by a factor of a few. Our measured RM and the RM previously reported by the CARMA and SMA observations cannot be consistently explained by the spherical accretion flow with a power-law profile. We propose that a clumpy/inhomogeneous ambient medium is responsible for the observed rotation measure. Using equipartition magnetic field, we derive the electron density of $2\times10^{4}$$\sim$cm$^{-3}$. Such an electron density is consistent with the cloud of narrow line emission region around the central engine. We also discuss the magnetic field configuration from black hole scale to pc scale and the origin of low polarization.

The Venusian atmosphere is in a state of superrotation where prevailing westward winds move much faster than the planet's rotation. Venus is covered with thick clouds that extend from about 45 to 70 km altitude, but thermal radiation emitted from the lower atmosphere and the surface on the planet's nightside escapes to space at narrow spectral windows of the near-infrared. The radiation can be used to estimate winds by tracking the silhouettes of clouds in the lower and middle cloud regions below about 57 km in altitude. Estimates of wind speeds have ranged from 50 to 70 m/s at low to mid-latitudes, either nearly constant across latitudes or with winds peaking at mid-latitudes. Here we report the detection of winds at low latitude exceeding 80 m/s using IR2 camera images from the Akatsuki orbiter taken during July and August 2016. The angular speed around the planetary rotation axis peaks near the equator, which we suggest is consistent with an equatorial jet, a feature that has not been observed previously in the Venusian atmosphere. The mechanism producing the jet remains unclear. Our observations reveal variability in the zonal flow in the lower and middle cloud region that may provide clues to the dynamics of Venus's atmospheric superrotation.

The general-relativistic magnetohydrodynamical (GRMHD) formulation for black hole-powered jets naturally gives rise to a stagnation surface, wherefrom inflows and outflows along magnetic field lines that thread the black hole event horizon originate. We derive a conservative formulation for the transport of energetic electrons which are initially injected at the stagnation surface and subsequently transported along flow streamlines. With this formulation the energy spectra evolution of the electrons along the flow in the presence of radiative and adiabatic cooling is determined. For flows regulated by synchrotron radiative losses and adiabatic cooling, the effective radio emission region is found to be finite, and geometrically it is more extended along the jet central axis. Moreover, the emission from regions adjacent to the stagnation surface is expected to be the most luminous as this is where the freshly injected energetic electrons concentrate. An observable stagnation surface is thus a strong prediction of the GRMHD jet model with the prescribed non-thermal electron injection. Future millimeter/sub-millimeter (mm/sub-mm) very-long-baseline interferometric (VLBI) observations of supermassive black hole candidates, such as the one at the center of M87, can verify this GRMHD jet model and its associated non-thermal electron injection mechanism.

This White Paper details the intentions and plans of the East Asian Very Long Baseline Interferometry (VLBI) community for pushing the frontiers of millimeter/submillimeter VLBI. To this end, we shall endeavor to actively promote coordinated efforts in the East Asia region. Our goal is to establish firm collaborations among the East Asia VLBI community in partnership with related institutes in North America and Europe and to expand existing global mm/submm VLBI arrays for (a) exploring the vicinity of black holes with an ultimate angular resolution down to 1~R$_{S}$ (Schwarzschild radius) and (b) investigating the dynamics of circumstellar gas in star-forming regions and late-type stars, and circumnuclear gas around active galactic nuclei (AGNs). In the first half of this White Paper, we highlight scientific accomplishments of the East Asia (EA) VLBI community. Various VLBI research results on M87, Sgr A*, blazars, narrow-line Seyfert~1 galaxies, and compact symmetric objects are described, and future visions of our VLBI science are briefly presented. Maser science of star formation, stellar evolution, and physics of accretion disks around AGNs are also discussed. A new vision for conducting multi-transition maser studies using mm/submm VLBI together with the Atacama Large Millimeter/Submillimeter Array (ALMA) is described. In the second half of this White Paper, we describe the EA community's vision for using mm/submm VLBI arrays in the framework or extended version of the Event Horizon Telescope (EHT) and the Global Millimeter VLBI Array (GMVA). The accomplishment of the aforementioned goal will maximize the overall scientific outcomes of mm/submm VLBI in the world.

The structure of the NGC 6251 jet on the milliarcsecond scale is investigated using images taken with the European VLBI Network and the Very Long Baseline Array. We detect a structural transition of the jet from a parabolic to a conical shape at a distance of (1-2) x 10^5 times the Schwarzschild radius from the central engine, which is close to the sphere of gravitational influence of the supermassive black hole (SMBH). We also examine the jet pressure profiles with the synchrotron minimum energy assumption to discuss the physical origin of the structural transition. The NGC 6251 jet, together with the M87 jet, suggests a fundamental process of structural transition in the jets of active galactic nuclei (AGNs). Collimated AGN jets are characterized by their external galactic medium, showing that AGN jets interplay with the SMBH and its host galaxy.

Since the ALMA North America Prototype Antenna was awarded to the Smithsonian Astrophysical Observatory (SAO), SAO and the Academia Sinica Institute of Astronomy & Astrophysics (ASIAA) are working jointly to relocate the antenna to Greenland. This paper shows the status of the antenna retrofit and the work carried out after the recommissioning and subsequent disassembly of the antenna at the VLA has taken place. The next coming months will see the start of the antenna reassembly at Thule Air Base. These activities are expected to last until the fall of 2017 when commissioning should take place. In parallel, design, fabrication and testing of the last components are taking place in Taiwan.

In the current paradigm, it is believed that the compact VLBI radio core of radio-loud AGNs represents the innermost upstream regions of relativistic outflows. These regions of AGN jets have generally been modeled by a conical outflow with roughly constant opening angle and flow speed. Nonetheless, some works suggest that a parabolic geometry would be more appropriate to fit the high energy spectral energy distribution properties and it has been recently found that, at least in some nearby radio-galaxies, the geometry of the innermost regions of the jet is parabolic. We compile here multi-frequency core sizes of archival data to investigate the typically unresolved upstream regions of the jet geometry of a sample of 56 radio-loud AGNs. Data combined from the sources considered here is not consistent with the classic picture of a conical jet starting in the vicinity of the super-massive black hole (SMBH), and may exclude a pure parabolic outflow solution, but rather suggest an intermediate solution with quasi-parabolic streams, which are frequently seen in numerical simulations. Inspection of the large opening angles near the SMBH and the range of the Lorentz factors derived from our results support our analyses. Our result suggests that the conical jet paradigm in AGNs needs to be re-examined by mm/sub-mm VLBI observations.

Nowadays there is compelling evidence for the existence of dark matter in the Universe. A general consensus has been expressed on the need for a directional sensitive detector to confirm, with a complementary approach, the candidates found in conventional searches and to finally extend their sensitivity beyond the limit of neutrino-induced background. We propose here the use of a detector based on nuclear emulsions to measure the direction of WIMP-induced nuclear recoils. The production of nuclear emulsion films with nanometric grains is established. Several measurement campaigns have demonstrated the capability of detecting sub-micrometric tracks left by low energy ions in such emulsion films. Innovative analysis technologies with fully automated optical microscopes have made it possible to achieve the track reconstruction for path lengths down to one hundred nanometers and there are good prospects to further exceed this limit. The detector concept we propose foresees the use of a bulk of nuclear emulsion films surrounded by a shield from environmental radioactivity, to be placed on an equatorial telescope in order to cancel out the effect of the Earth rotation, thus keeping the detector at a fixed orientation toward the expected direction of galactic WIMPs. We report the schedule and cost estimate for a one-kilogram mass pilot experiment, aiming at delivering the first results on the time scale of six years.

We investigate the magnetic field configuration of the M87 jet at arc-second scales by using archival polarimetric VLA data at 8, 15, 22 and 43 GHz. By stacking images over three years in order to enhance the sensitivity, we reveal, for the first time, systematic transverse gradients of the Faraday rotation measure in several knots along the jet. Combining this result with polarization properties and the dynamics of the jet, we suggest the magnetic structure in several knots at kiloparsec scales consists of a systematically wrapped, tightly wound helical configuration. Our analysis brings us a new paradigm where the M87 jet is a fundamentally current carrying system produced in the vicinity of the supermassive black hole, transferring a huge amount of the electromagnetic energy over the host galaxy scale.

We investigate the distribution of Faraday rotation measure (RM) in the M87 jet at arc-second scales by using archival polarimetric VLA data at 8, 15, 22 and 43 GHz. We resolve the structure of the RM in several knots along the jet for the first time. We derive the power spectrum in the arcsecond scale jet and find indications that the RM cannot be associated with a turbulent magnetic field with 3D Kolmogorov spectrum. Our analysis indicates that the RM probed on jet scales has a significant contribution of a Faraday screen associated with the vicinity of the jet, in contrast with that on kiloparsec scales, typically assumed to be disconnected from the jet. Comparison with previous RM analyses suggests that the magnetic fields giving rise to the RMs observed in jet scales have different properties and are well less turbulent than these observed in the lobes.

Ground-based observations at terahertz (THz) frequencies are a newly explorable area of astronomy for the next ten years. We discuss science cases for a first-generation 10-m class THz telescope, focusing on the Greenland Telescope as an example of such a facility. We propose science cases and provide quantitative estimates for each case. The largest advantage of ground-based THz telescopes is their higher angular resolution (~ 4 arcsec for a 10-m dish), as compared to space or airborne THz telescopes. Thus, high-resolution mapping is an important scientific argument. In particular, we can isolate zones of interest for Galactic and extragalactic star-forming regions. The THz windows are suitable for observations of high-excitation CO lines and [N II] 205 um lines, which are scientifically relevant tracers of star formation and stellar feedback. Those lines are the brightest lines in the THz windows, so that they are suitable for the initiation of ground-based THz observations. THz polarization of star-forming regions can also be explored since it traces the dust population contributing to the THz spectral peak. For survey-type observations, we focus on ``sub-THz'' extragalactic surveys, whose uniqueness is to detect galaxies at redshifts z ~ 1--2, where the dust emission per comoving volume is the largest in the history of the Universe. Finally we explore possibilities of flexible time scheduling, which enables us to monitor active galactic nuclei, and to target gamma-ray burst afterglows. For these objects, THz and submillimeter wavelength ranges have not yet been explored.

We report on 230 GHz (1.3 mm) VLBI observations of M87 with the Event Horizon Telescope using antennas on Mauna Kea in Hawaii, Mt. Graham in Arizona and Cedar Flat in California. For the first time, we have acquired 230 GHz VLBI interferometric phase information on M87 through measurement of closure phase on the triangle of long baselines. Most of the measured closure phases are consistent with 0$^{\circ}$ as expected by physically-motivated models for 230 GHz structure such as jet models and accretion disk models. The brightness temperature of the event-horizon-scale structure is $\sim 1 \times 10^{10}$ K derived from the compact flux density of $\sim 1$ Jy and the angular size of $\sim 40 $ $\rm \mu$as $\sim$ 5.5 $R_{{\rm s}}$, which is broadly consistent with the peak brightness of the radio cores at 1-86 GHz located within $\sim 10^2$ $R_{{\rm s}}$. Our observations occurred in the middle of an enhancement in very-high-energy (VHE) $\rm \gamma$-ray flux, presumably originating in the vicinity of the central black hole. Our measurements, combined with results of multi-wavelength observations, favor a scenario in which the VHE region has an extended size of $\sim$20-60 $R_{{\rm s}}$.

General relativistic magnetohydrodynamic (GRMHD) flows along magnetic fields threading a black hole can be divided into inflow and outflow parts, according to the result of the competition between the black hole gravity and magneto-centrifugal forces along the field line. Here we present the first self-consistent, semi-analytical solution for a cold, Poynting flux-dominated (PFD) GRMHD flow, which passes all four critical (inner and outer, Alfven and fast magnetosonic) points along a parabolic streamline. By assuming that the dominating (electromagnetic) component of the energy flux per flux tube is conserved at the surface where the inflow and outflow are separated, the outflow part of the solution can be constrained by the inflow part. The semi-analytical method can provide fiducial and complementary solutions for GRMHD simulations around the rotating black hole, given that the black hole spin, global streamline, and magnetizaion (i.e., a mass loading at the inflow/outflow separation) are prescribed. For reference, we demonstrate a self-consistent result with the work by McKinney in a quantitative level.

The ALMA North America Prototype Antenna was awarded to the Smithsonian Astrophysical Observatory (SAO) in 2011. SAO and the Academia Sinica Institute of Astronomy & Astrophysics (ASIAA), SAO's main partner for this project, are working jointly to relocate the antenna to Greenland to carry out millimeter and submillimeter VLBI observations. This paper presents the work carried out on upgrading the antenna to enable operation in the Arctic climate by the GLT Team to make this challenging project possible, with an emphasis on the unexpected telescope components that had to be either redesigned or changed. Five-years of inactivity, with the antenna laying idle in the desert of New Mexico, coupled with the extreme weather conditions of the selected site in Greenland have it necessary to significantly refurbish the antenna. We found that many components did need to be replaced, such as the antenna support cone, the azimuth bearing, the carbon fiber quadrupod, the hexapod, the HVAC, the tiltmeters, the antenna electronic enclosures housing servo and other drive components, and the cables. We selected Vertex, the original antenna manufacturer, for the main design work, which is in progress. The next coming months will see the major antenna components and subsystems shipped to a site of the US East Coast for test-fitting the major antenna components, which have been retrofitted. The following step will be to ship the components to Greenland to carry out VLBI and single dish observations. Antenna reassembly at Summit Station should take place during the summer of 2018.

A 12-m diameter radio telescope will be deployed to the Summit Station in Greenland to provide direct confirmation of a Super Massive Black Hole (SMBH) by observing its shadow image in the active galaxy M87. The telescope (Greenland Telescope: GLT) is to become one of the Very Long Baseline Interferometry (VLBI) stations at sub-millimeter (submm) regime, providing the longest baseline > 9,000 km to achieve an exceptional angular resolution of 20 micro arc sec at 350 GHz, which will enable us to resolve the shadow size of ~40 micro arc sec. The triangle with the longest baselines formed by the GLT, the Atacama Large Millimeter/submillimeter Array (ALMA) in Chile, and the Submillimeter Array (SMA) in Hawaii will play a key role for the M87 observations. We have been working on the image simulations based on realistic conditions for a better understanding of the possible observed images. In parallel, retrofitting of the telescope and the site developments are in progress. Based on three years of opacity monitoring at 225 GHz, our measurements indicate that the site is excellent for submm observations, comparable to the ALMA site. The GLT is also expected to make single-dish observations up to 1.5 THz.

We describe a new paradigm for understanding both relativistic motions and particle acceleration in the M87 jet: a magnetically dominated relativistic flow that naturally produces four relativistic magnetohydrodynamic (MHD) shocks (forward/reverse fast and slow modes). We apply this model to a set of optical super- and subluminal motions discovered by Biretta and coworkers with the {\em Hubble Space Telescope} during 1994 -- 1998. The model concept consists of ejection of a {\em single} relativistic Poynting jet, which possesses a coherent helical (poloidal + toroidal) magnetic component, at the remarkably flaring point HST-1. We are able to reproduce quantitatively proper motions of components seen in the {\em optical} observations of HST-1 with the same model we used previously to describe similar features in radio VLBI observations in 2005 -- 2006. This indicates that the quad relativistic MHD shock model can be applied generally to recurring pairs of super/subluminal knots ejected from the upstream edge of the HST-1 complex as observed from radio to optical wavelengths, with forward/reverse fast-mode MHD shocks then responsible for observed moving features. Moreover, we identify such intrinsic properties as the shock compression ratio, degree of magnetization, and magnetic obliquity and show that they are suitable to mediate diffusive shock acceleration of relativistic particles via the first-order Fermi process. We suggest that relativistic MHD shocks in Poynting-flux dominated helical jets may play a role in explaining observed emission and proper motions in many AGNs.

We present a study of sub-pc scale radio structure of the radio galaxy 3C 84/NGC 1275 based on the Very Long Baseline Array (VLBA) data at 43 GHz. We discover a limb-brightening in the "restarted" jet associated with the 2005 radio outburst. In the 1990s, the jet structure was ridge-brightening rather than limb-brightening, despite the observations being done with similar angular resolution. This indicates that the transverse jet structure has changed recently. This change in the morphology shows an interesting agreement with the $\gamma$-ray flux increase, i.e., the $\gamma$-ray flux in 1990s was at least seven times lower than the current one. One plausible explanation for the limb-brightening is the velocity structure of the jet in the context of the stratified jet, which is a successful scenario to explain the $\gamma$-ray emission in some active galactic nuclei (AGNs). If this is the case, the change in apparent transverse structure might be caused by the change in the transverse velocity structure. We argue the possibility that the transition from ridge-brightening to limb-brightening is related to the $\gamma$-ray time variability on the timescale of decades. We also discuss the collimation profile of the jet.

We present the first constraint on Faraday rotation measure (RM) at submillimeter wavelengths for the nucleus of M 87. By fitting the polarization position angles ($\chi$) observed with the SMA at four independent frequencies around $\sim$230 GHz and interpreting the change in $\chi$ as a result of \emph{external} Faraday rotation associated with accretion flow, we determine the rotation measure of the M 87 core to be between $-$7.5$\times$10$^{5}$ and 3.4$\times$10$^{5}$ rad/m$^{2}$. Assuming a density profile of the accretion flow that follows a power-law distribution and a magnetic field that is ordered, radial, and has equipartition strength, the limit on the rotation measure constrains the mass accretion rate $\dot{M}$ to be below 9.2$\times$10$^{-4}$ M$_{\odot}$~yr$^{-1}$ at a distance of 21 Schwarzchild radii from the central black hole. This value is at least two orders of magnitude smaller than the Bondi accretion rate, suggesting significant suppression of the accretion rate in the inner region of the accretion flow. Consequently, our result disfavors the classical \emph{advection dominated accretion flow} (ADAF) and prefers the \emph{adiabatic inflow-outflow solution} (ADIOS) or \emph{convection-dominated accretion flow} (CDAF) for the hot accretion flow in M 87.

The central cD galaxy of the Hydra A cluster has one of the most powerful active galactic nuclei (AGNs) in the nearby Universe (z <0.2). We report on the discovery of a dust lane in the cD galaxy using Subaru telescope. The i'-band image shows the existence of a dark band of the size of 3.6"x 0.7" (4 kpc x 0.8 kpc), which appears to be quite similar to the dust lane observed in Centaurus A. The morphology indicates that the cold disk that seen as the dust lane is almost edge-on and rotates around the AGN. Since the minor axis of the dust lane is nearly parallel to the radio jets emerging from the AGN, the disk is probably feeding its gas into the central black hole. From the absorption, we estimate the hydrogen column density of the lane is N_H=2.0 x 10^21 cm^-2, and the mass of the disk is ~8 x 10^7 M_sun. The column density is consistent with constraints obtained from Chandra X-ray observations. The age of the disk is >~ 4 x 10^7 yr. The position angle of the disk and the galaxy's photometric axis are misaligned, which may imply that the cold gas in the disk is brought via galaxy mergers. Our observations may indicate that the supply of cold gas by galaxy mergers is required for the most intensive feedback from AGNs.

The velocity field of the M87 jet from milli-arcsecond (mas) to arcsecond scales is extensively investigated together with new radio images taken by EVN observations. We detected proper motions of components located at between 160 mas from the core and the HST-1 complex for the first time. Newly derived velocity fields exhibits a systematic increase from sub-to-superluminal speed in the upstream of HST-1. If we assume that the observed velocities reflect the bulk flow, we here suggest that the M87 jet may be gradually accelerated through a distance of 10^6 times of the Schwarzschild radius of the supermassive black hole. The acceleration zone is co-spatial with the jet parabolic region, which is interpreted as the collimation zone of the jet (Asada & Nakamura 2012). The acceleration and collimation take place simultaneously, which we suggest a characteristic of magnetohydrodynamic flows. Distribution of the velocity field has a peak at HST-1, which is considered as the site of over-collimation, and shows a deceleration downstream of HST-1 where the jet is conical. Our interpretation of the velocity map in the M87 jet gives a hypothesis in AGNs that the acceleration and collimation zone of relativistic jets extends over the whole scale within the sphere of influence of the supermassive black hole.

The GLT project is deploying a new submillimeter (submm) VLBI station in Greenland. Our primary scientific goal is to image a shadow of the supermassive black hole (SMBH) of six billion solar masses in M87 at the center of the Virgo cluster of galaxies. The expected SMBH shadow size of 40-50 $\mu$as requires superbly high angular resolution, suggesting that the submm VLBI would be the only way to obtain the shadow image. The Summit station in Greenland enables us to establish baselines longer than 9,000 km with ALMA in Chile and SMA in Hawaii as well as providing a unique $u$--$v$ coverage for imaging M87. Our VLBI network will achieve a superior angular resolution of about 20 $\mu$as at 350 GHz, corresponding to $\sim2.5$ times of the Schwarzschild radius of the supermassive black hole in M87. We have been monitoring the atmospheric opacity at 230 GHz since August. 2011; we have confirmed the value on site during the winter season is comparable to the ALMA site thanks to high altitude of 3,200 m and low temperature of $-50\degr$C. We will report current status and future plan of the GLT project towards our expected first light on 2015--2016.

An international consortium is presently constructing a beamformer for the Atacama Large Millimeter/submillimeter Array (ALMA) in Chile that will be available as a facility instrument. The beamformer will aggregate the entire collecting area of the array into a single, very large aperture. The extraordinary sensitivity of phased ALMA, combined with the extremely fine angular resolution available on baselines to the Northern Hemisphere, will enable transformational new very long baseline interferometry (VLBI) observations in Bands 6 and 7 (1.3 and 0.8 mm) and provide substantial improvements to existing VLBI arrays in Bands 1 and 3 (7 and 3 mm). The ALMA beamformer will have impact on a variety of scientific topics, including accretion and outflow processes around black holes in active galactic nuclei (AGN), tests of general relativity near black holes, jet launch and collimation from AGN and microquasars, pulsar and magnetar emission processes, the chemical history of the universe and the evolution of fundamental constants across cosmic time, maser science, and astrometry.

We report new results from an HST archival program to study proper motions in the optical jet of the nearby radio galaxy M87. Using over 13 years of archival imaging, we reach accuracies below 0.1c in measuring the apparent velocities of individual knots in the jet. We confirm previous findings of speeds up to 4.5c in the inner 6" of the jet, and report new speeds for optical components in the outer part of the jet. We find evidence of significant motion transverse to the jet axis on the order of 0.6c in the inner jet features, and superluminal velocities parallel and transverse to the jet in the outer knot components, with an apparent ordering of velocity vectors possibly consistent with a helical jet pattern. Previous results suggested a global deceleration over the length of the jet in the form of decreasing maximum speeds of knot components from HST-1 outward, but our results suggest that superluminal speeds persist out to knot C, with large differentials in very nearby features all along the jet. We find significant apparent accelerations in directions parallel and transverse to the jet axis, along with evidence for stationary features in knots D, E, and I. These results are expected to place important constraints on detailed models of kpc-scale relativistic jets.

The structure and dynamics of the M87 jet from sub-milli-arcsec to arcsecond scales are continuously examined. We analysed the VLBA archival data taken at 43 and 86 GHz to measure the size of VLBI cores. Millimeter/sub-mm VLBI cores are considered as innermost jet emissions, which has been originally suggested by Blandford & K\"onigl. Those components fairly follow an extrapolated parabolic streamline in our previous study so that the jet has a single power-law structure with nearly five orders of magnitude in the distance starting from the vicinity of the supermassive black hole (SMBH), less than 10 Schwarzschild radius ($r_{\rm s}$). We further inspect the jet parabolic structure as a counterpart of the magnetohydrodynamic (MHD) nozzle in order to identify the property of a bulk acceleration. We interpret that the parabolic jet consists of Poynting-flux dominated flows, powered by large amplitude, nonlinear torsional Alfv\'en waves. We examine the non-relativistic MHD nozzle equation in a parabolic shape. The nature of trans-fast magnetosonic flow is similar to the one of transonic solution of Parker's hydrodynamic solar wind; the jet becomes super-escape as well as super-fast magnetosonic at around $\sim 10^{3}\,r_{\rm s}$, while the upstream trans-Alfv\'enic flow speed increases linearly as a function of the distance at $\sim 10^{2}$ - $10^{3}\, r_{\rm s}$. We here point out that this is the first evidence to identify these features in astrophysical jets. We propose that the M87 jet is magnetically accelerated, but thermally confined by the stratified ISM inside the sphere of gravitational influence of the SMBH potential, which may be a norm in AGN jets.

Detecting the first electron pairs with nuclear emulsion allows a precise measurement of the direction of incident gamma-rays as well as their polarization. With recent innovations in emulsion scanning, emulsion analyzing capability is becoming increasingly powerful. Presently, we are developing a balloon-borne gamma-ray telescope using nuclear emulsion. An overview and a status of our telescope is given.

The structure of the M87 jet, from milli-arcsec to arcsecond scales, is extensively investigated, utilizing the images taken with the EVN, MERLIN and VLBA. We discover that the jet maintains a parabolic streamline over a range in sizescale equal to $10^{5}$ times the Schwarzschild radius. The jet then transitions into a conical shape further downstream. This suggests that the magnetohydrodynamic jet is initially subjected to the confinement by the external gas which is dominated by the gravitational influence of the supermassive black hole. Afterwards the jet is then freely expanding with a conical shape. This geometrical transition indicates that the origin of the HST-1 complex may be a consequence of the over-collimation of the jet. Our result suggests that when even higher angular resolution is provided by a future submm VLBI experiment, we will be able to explore the origin of active galactic nuclei jets.

We carried out the intercomparison experiments between VLBI and other techniques to show the capability of VLBI time and frequency transfer by using the current geodetic VLBI technique and facilities as the summary of the experiments that we carried out since 2007. The results from the two different types of experiments show that the VLBI is more stable than GPS but is slightly noisier than two new two-way techniques (TW(DPN), ETS8(TCE)), and VLBI can measure the correct time difference as same as ETS8(TCE).

In this study, we are doing R&D for directional dark matter search with nuclear emulsion. First of all, higher resolution nuclear emulsion with fine silver halide crystals was developed in the production facility of emulsion at Nagoya university, and we confirmed that it can detect the expected nuclear recoil tracks. The readout of submicron tracks was required the new technology. We developed the expansion technique, and could readout the signal by shape analysis with optical microscopy. The two dimensional angular resolution is 36 degrees at the original track length of range from 150nm to 200nm with optical microscopy. Finally we demonstrated by using recoiled nuclei induced by 14.8MeV neutron, and confirmed the technique.Moreover, we developed the X-ray microscope system with SPring-8 as final check with higher resolution of selected candidate tracks with optical microscopy. The angular resolution was improved from 31 degrees with optical microscopy to 17degrees with X-ray microscopy at the track length of range from 150nm to 250nm. We are developing the practical system and planning for start of the test running with prototype detector.

We report a measurement of the neutrino-electron elastic scattering rate from 8B solar neutrinos based on a 123 kton-day exposure of KamLAND. The background-subtracted electron recoil rate, above a 5.5 MeV analysis threshold is 1.49+/-0.14(stat)+/-0.17(syst) events per kton-day. Interpreted as due to a pure electron flavor flux with a 8B neutrino spectrum, this corresponds to a spectrum integrated flux of 2.77+/-0.26(stat)+/-0.32(syst) x 10^6 cm^-2s^-1. The analysis threshold is driven by 208Tl present in the liquid scintillator, and the main source of systematic uncertainty is due to background from cosmogenic 11Be. The measured rate is consistent with existing measurements and with Standard Solar Model predictions which include matter enhanced neutrino oscillation.

This is the first in a series of papers that introduces a new paradigm for understanding the jet in M87: a collimated relativistic flow in which strong magnetic fields play a dominant dynamical role. Here wefocus on the flow downstream of HST-1 - an essentially stationary flaring feature that ejects trails of superluminal components. We propose that these components are quad relativistic magnetohydrodynamic shock fronts (forward/reverse fast and slow modes) in a narrow jet with a helically twisted magnetic structure. And we demonstrate the properties of such shocks with simple one-dimensional numerical simulations. Quasi-periodic ejections of similar component trails may be responsible for the M87 jet substructures observed further downstream on 100 - 1,000 pc scales. This new paradigm requires the assimilation of some new concepts into the astrophysical jet community, particularly the behavior of slow/fast-mode waves/shocks and of current-driven helical kink instabilities. However, the prospects of these ideas applying to a large number of other jet systems may make this worth the effort.

We present the case for a dark matter detector with directional sensitivity. This document was developed at the 2009 CYGNUS workshop on directional dark matter detection, and contains contributions from theorists and experimental groups in the field. We describe the need for a dark matter detector with directional sensitivity; each directional dark matter experiment presents their project's status; and we close with a feasibility study for scaling up to a one ton directional detector, which would cost around $150M.

We apply magnetohydrodynamic (MHD) modeling to the radio galaxy Hercules A for investigating the jet-driven shock, jet/lobe transition, wiggling, and magnetic field distribution associated with this source. The model consists of magnetic tower jets in a galaxy cluster environment, which has been discussed in a series of our papers. The profile of underlying ambient gas plays an important role in jet-lobe morphology. The balance between the magnetic pressure generated by axial current and the ambient gas pressure can determine the lobe radius. The jet body is confined jointly by the external pressure and gravity inside the cluster core radius R_c, while outside R_c it expands radially to form fat lobes in a steeply decreasing ambient thermal pressure gradient. The current-carrying jets are responsible for generating a strong, tightly wound helical magnetic field. This magnetic configuration will be unstable against the current-driven kink mode and it visibly grows beyond R_c where a separation between the jet forward and return currents occurs. The reversed pinch profile of global magnetic field associated with the jet and lobes produces projected B-vector distributions aligned with the jet flow and the lobe edge. AGN-driven shock powered by the expanding magnetic tower jet surrounds the jet/lobe structure and heats the ambient ICM. The lobes expand subsonically; no obvious hot spots are produced at the heads of lobes. Several key features in our MHD modeling may be qualitatively supported by the observations of Hercules A.

We present results from multifrequency polarimetry of NRAO 140 using the Very Long Baseline Array. These observations allow us to reveal the distributions of both the polarization position angle and the Faraday rotation measure (RM). These distributions are powerful tools to discern the projected and line-of-sight components of the magnetic field, respectively. We find a systematic gradient in the RM distribution, with its sign being opposite at either side of the jet with respect to the jet axis. The sign of the RM changes only with the direction of the magnetic field component along the line of sight, so this can be explained by the existence of helical magnetic components associated with the jet itself. We derive two constraints for the pitch angle of the helical magnetic field from the distributions of the RM and the projected magnetic field; the RM distribution indicates that the helical fields are tightly wound, while that of the projected magnetic field suggests they are loosely wound around the jet axis. This inconsistency may be explained if the Faraday rotator is not cospatial with theemitting region. Our results may point toward a physical picture in which an ultra-relativistic jet (spine) with a loosely wound helical magnetic field is surrounded by a sub-relativistic wind layer (sheath) with a tightly wound helical magnetic field.

We present results from magneto-hydrodynamic jet simulations in the magnetically dominated regime and find that they inflate cavities into the intracluster medium. Mock Chandra X-ray observations of our simulations show a strong resemblance to real cavities observed in X-ray images of clusters, both in terms of their morphology and thermodynamic structure. An analysis of the evolution of bubble sizes in the multi-cavity system Hydra A, as well as in a large sample of 64 cavities in 32 clusters shows that bubbles tend to expand much faster than expected in the purely adiabatic regime. Instead, we find that the bubbles follow more closely a trend predicted by our current-carrying jet models.

We argue the behavior of Poynting flux-dominated outflows from AGN in the galactic cluster systems by performing three-dimensional MHD simulations within the framework of the "magnetic tower" model. Of particular interests are the structure of MHD waves, the cylindrical radial force balance, the (de)collimation, and the stability properties of magnetic tower jets. Transition between the jet/lobe and the formation of wiggling jet by growing current-driven instability are discussed.

Stability properties of ``magnetic tower'' jets propagating in the gravitationally stratified background have been examined by performing three-dimensional magnetohydrodynamic simulations. The current-carrying, Poynting flux-dominated magnetic tower jet, which possesses a highly wound helical magnetic field, is subject to the current-driven instability (CDI). We find that, under general physical conditions including small perturbations in the initial background profiles, the propagating magnetic tower jets develop the non-axisymmetric, $m=1$ kink mode of the CDI. The kink mode grows on the local Alfv\'en crossing time scale. In addition, two types of kink modes appear in the system. At the central region where external thermal pressure confinement is strong, only the internal kink mode is excited and will grow. A large distance away from the central region where the external thermal pressure becomes low, the external kink mode is observed. As a result, the exterior of magnetic tower jets will be deformed into a large-scale wiggled structure. We also discuss extensively the different physical processes that contribute to the overall stability properties of the magnetic tower jets. Specifically, when the jet propagates in an initially unperturbed background, we find that they can survive the kink mode beyond the point predicted by the well-known Kruskal-Shafranov (K-S) criterion. The stabilization in this case comes mainly from the dynamical relaxation of magnetic twists during the propagation of magnetic towers; the magnetic pitch is reduced and the corresponding K-S critical wavelength becomes longer as the tower jet proceeds. Furthermore, we show that the pressure-driven and Kelvin-Helmholtz instabilities do not occur in the magnetic tower jets.

Based on a new approach on modeling the magnetically dominated outflows from AGNs (Li et al. 2006), we study the propagation of magnetic tower jets in gravitationally stratified atmospheres (such as a galaxy cluster environment) in large scales ($>$ tens of kpc) by performing three-dimensional magnetohydrodynamic (MHD) simulations. We present the detailed analysis of the MHD waves, the cylindrical radial force balance, and the collimation of magnetic tower jets. As magnetic energy is injected into a small central volume over a finite amount of time, the magnetic fields expand down the background density gradient, forming a collimated jet and an expanded ``lobe'' due to the gradually decreasing background density and pressure. Both the jet and lobes are magnetically dominated. In addition, the injection and expansion produce a hydrodynamic shock wave that is moving ahead of and enclosing the magnetic tower jet. This shock can eventually break the hydrostatic equilibrium in the ambient medium and cause a global gravitational contraction. This contraction produces a strong compression at the head of the magnetic tower front and helps to collimate radially to produce a slender-shaped jet. At the outer edge of the jet, the magnetic pressure is balanced by the background (modified) gas pressure, without any significant contribution from the hoop stress. On the other hand, along the central axis of the jet, hoop stress is the dominant force in shaping the central collimation of the poloidal current. The system, which possesses a highly wound helical magnetic configuration, never quite reaches a force-free equilibrium state though the evolution becomes much slower at late stages. The simulations were performed without any initial perturbations so the overall structures of the jet remain mostly axisymmetric.

We argue that the Maxwellian approximation can essentially underestimate the rates of some nuclear reactions in hot plasma under conditions very close to thermal equilibrium. This phenomenon is demonstrated explicitly on the example of reactions in self-sustained DT fusion plasma with admixture of light elements X = Li, Be, C. A kinetic analysis shows that the reactivity enhancement results from non-Maxwellian knock-on perturbations of ion distributions caused by close collisions with energetic fusion products. It is found that although the fraction of the knock-on ions is small, these particles appreciably affect the D+X and T+X reaction rates. The phenomenon discussed is likely to have general nature and can play role in other laboratory and probably astrophysical plasma processes.

Non-relativistic three-dimensional magnetohydrodynamical (MHD) simulations of Poynting flux dominated (PFD) jets are presented. Our study focuses on the propagation of strongly magnetized hypersonic, but sub-Alfv\'enic ($C^{2}_{\rm s} \ll V^{2}_{\rm jet} < V^{2}_{\rm A}$) flow and on the subsequent development of a current-driven (CD) kink instability. This instability may be responsible for the ``wiggled'' structures seen in sub-parsec scale (VLBI) jets. In the present paper, we investigate the nonlinear behavior of PFD jets in a variety of external ambient magnetized gas distributions, including those with density, pressure, and temperature gradients. Our numerical results show that the jets can develop CD distortions in the trans-Alfv\'enic flow case, even when the flow itself is still strongly magnetically dominated. An internal non-axisymmetric body mode grows on time scales of order of the Alfv\'en crossing time and distorts the structure and magnetic configuration of the jet. The kink ($m=1$) mode of the CD instability, driven by the radial component of the Lorentz force, grows faster than other higher order modes ($m>1$). In the jet frame the mode grows locally and expands radially at each axial position where the jet is unstable: the instability, therefore, does not propagate as a wave along the jet length. A naturally-occurring, external helically magnetized wind, which is (quasi-) axially current-free, surrounds the well-collimated current-carrying jet and reduces velocity shear between the jet and external medium. This stabilizes the growth of MHD Kelvin-Helmholtz surface modes in the inner jet flow.

Distributions of Faraday rotation measure (FRM) and the projected magnetic field derived by a 3-dimensional simulation of MHD jets are investigated based on our "sweeping magnetic twist model". FRM and Stokes parameters were calculated to be compared with radio observations of large scale wiggled AGN jets on kpc scales. We propose that the FRM distribution can be used to discuss the 3-dimensional structure of magnetic field around jets and the validity of existing theoretical models, together with the projected magnetic field derived from Stokes parameters. In the previous paper, we investigated the basic straight part of AGN jets by using the result of a 2-dimensional axisymmetric simulation. The derived FRM distribution has a general tendency to have a gradient across the jet axis, which is due to the toroidal component of the magnetic field generated by the rotation of the accretion disk. In this paper, we consider the wiggled structure of the AGN jets by using the result of a 3-dimensional simulation. Our numerical results show that the distributions of FRM and the projected magnetic field have a clear correlation with the large scale structure of the jet itself, namely, 3-dimensional helix. Distributions, seeing the jet from a certain direction, show a good matching with those in a part of 3C449 jet. This suggests that the jet has a helical structure and that the magnetic field (especially the toroidal component) plays an important role in the dynamics of the wiggle formation because it is due to a current-driven helical kink instability in our model.

We review recent progress in the theory of jet production, with particular emphasis on the possibility of 1) powerful jets being produced in the first few seconds after collapse of a supernova core and 2) those jets being responsible for the asymmetric explosion itself. The presently favored jet-production mechanism is an electrodynamic one, in which charged plasma is accelerated by electric fields that are generated by a rotating magnetic field anchored in the protopulsar. Recent observations of Galactic jet sources provide important clues to how {\em all} such sources may be related, both in the physical mechanism that drives the jet in the astrophysical mechanisms that create conditions conducive to jet formation. We propose a grand evolutionary scheme that attempts to unify these sources on this basis, with MHD supernovae providing the missing link. We also discuss several important issues that must be resolved before this (or another scheme) can be adopted.

Non-relativistic three-dimensional magnetohydrodynamic simulations of Poynting-flux-dominated (PFD) jets are presented. Our study focuses on the propagation of strongly magnetized hypersonic but sub-Alfv\'enic flow ($C_{\rm s}^2 << V_{\rm jet}^2 < V_{\rm A}^2$) and the development of a current-driven (CD) kink instability. This instability may be responsible for the "wiggled" structures seen in VLBI-scale AGN jets. In the present paper we investigate the nonlinear behavior of PFD jets in a variety of external ambient magnetized gas distributions, including those with density, pressure, and temperature gradients. Our numerical results show that PFD jets can develop kink distortions in the trans-Alfv\'enic flow case, even when the flow itself is still strongly magnetically dominated. In the nonlinear development of the instability, a non-axisymmetric mode grows on time scales of order the Alfv\'en crossing time (in the jet frame) and proceeds to disrupt the kinematic and magnetic structure of the jet. Because of a large scale poloidal magnetic field in the ambient medium, the growth of surface modes ({\it i.e.}, MHD Kelvin-Helmholtz instabilities) is suppressed. The CD kink mode ($m = 1$) grows faster than the other higher order modes ($m > 1$), driven in large part by the radial component of the Lorentz force.

Using the numerical data of MHD simulation for AGN jets based on our ``sweeping magnetic twist model'', we calculated the Faraday rotation measure (FRM) and the Stokes parameters to compare with observations. We propose that the FRM distribution can be used to discuss the 3-dimensional structure of magnetic field around jets, together with the projected magnetic field derived from the Stokes parameters. In the present paper, we supposed the basic straight part of AGN jet, and used the data of axisymmetric simulation. The FRM distribution we derived has a general tendency to have gradient across the jet axis, which is due to the toroidal component of the helical magnetic field generated by the rotation of the accretion disk. This kind of gradient in the FRM distribution is actually observed in some AGN jets (e.g. Asada et al. 2002), which suggests helical magnetic field around the jets and thus supports our MHD model. Following this success, we are now extending our numerical observation to the wiggled part of the jets using the data of 3-dimensional simulation based on our model in the following paper.