Recent observations of SgrA* by the GRAVITY instrument have astrometrically tracked infrared flares (IR) at distances of $\sim 10$ gravitational radii ($r_g$). In this paper, we study a model for the flares based on 3D general relativistic magnetohydrodynamic (GRMHD) simulations of magnetically arrested accretion disks (MADs) which exhibit violent episodes of flux escape from the black hole magnetosphere. These events are attractive for flare modeling for several reasons: i) the magnetically dominant regions can resist being disrupted via magneto-rotational turbulence and shear, ii) the orientation of the magnetic field is predominantly vertical as suggested by the GRAVITY data, iii) magnetic reconnection associated with the flux eruptions could yield a self-consistent means of particle heating/acceleration during the flare events. In this analysis we track erupted flux bundles and provide distributions of sizes, energies and plasma parameter. In our simulations, the orbits tend to circularize at a range of radii from $\sim 5-40 r_g$. The magnetic energy contained within the flux bundles ranges up to $\sim10^{40}$ erg, enough to power IR and X-ray flares. We find that the motion within the magnetically supported flow is substantially sub-Keplerian, in tension with the inferred period-radius relation of the three GRAVITY flares.

The Event Horizon Telescope recently observed the first shadow of a black hole. Images like this can potentially be used to test or constrain theories of gravity and deepen the understanding in plasma physics at event horizon scales, which requires accurate parameter estimations. In this work, we present Deep Horizon, two convolutional deep neural networks that recover the physical parameters from images of black hole shadows. We investigate the effects of a limited telescope resolution and observations at higher frequencies. We trained two convolutional deep neural networks on a large image library of simulated mock data. The first network is a Bayesian deep neural regression network and is used to recover the viewing angle $i$, and position angle, mass accretion rate $\dot{M}$, electron heating prescription $R_{\rm high}$ and the black hole mass $M_{\rm BH}$. The second network is a classification network that recovers the black hole spin $a$. We find that with the current resolution of the Event Horizon Telescope, it is only possible to accurately recover a limited number of parameters of a static image, namely the mass and mass accretion rate. Since potential future space-based observing missions will operate at frequencies above 230 GHz, we also investigated the applicability of our network at a frequency of 690 GHz. The expected resolution of space-based missions is higher than the current resolution of the Event Horizon Telescope, and we show that Deep Horizon can accurately recover the parameters of simulated observations with a comparable resolution to such missions.

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.

One of the key questions in the study of relativistic jets is how magnetic reconnection occurs and whether it can effectively accelerate electrons in the jet. We performed 3D particle-in-cell (PIC) simulations of a relativistic electron-proton jet of relatively large radius that carries a helical magnetic field. We focussed our investigation on the interaction between the jet and the ambient plasma and explore how the helical magnetic field affects the excitation of kinetic instabilities such as the Weibel instability (WI), the kinetic Kelvin-Helmholtz instability (kKHI), and the mushroom instability (MI). In our simulations these kinetic instabilities are indeed excited, and particles are accelerated. At the linear stage we observe recollimation shocks near the center of the jet. As the electron-proton jet evolves into the deep nonlinear stage, the helical magnetic field becomes untangled due to reconnection-like phenomena, and electrons are repeatedly accelerated as they encounter magnetic-reconnection events in the turbulent magnetic field.

One of the main dissipation processes acting on all scales in relativistic jets is thought to be governed by magnetic reconnection. Such dissipation processes have been studied in idealized environments, such as reconnection layers, which evolve in merging islands and lead to the production of plasmoids, ultimately resulting in efficient particle acceleration. In accretion flows onto black holes, reconnection layers can be developed and destroyed rapidly during the turbulent evolution of the flow. We present a series of two-dimensional general-relativistic magnetohydrodynamic simulations of tori accreting onto rotating black holes focusing our attention on the formation and evolution of current sheets. Initially, the tori are endowed with a poloidal magnetic field having a multi-loop structure along the radial direction and with an alternating polarity. During reconnection processes, plasmoids and plasmoid chains are developed leading to a flaring activity and hence to a variable electromagnetic luminosity. We describe the methods developed to track automatically the plasmoids that are generated and ejected during the simulation, contrasting the behaviour of multi-loop initial data with that encountered in typical simulations of accreting black holes having initial dipolar field composed of one loop only. Finally, we discuss the implications that our results have on the variability to be expected in accreting supermassive black holes.

High-resolution imaging of supermassive black holes is now possible, with new applications to testing general relativity and horizon-scale accretion and relativistic jet formation processes. Over the coming decade, the EHT will propose to add new strategically placed VLBI elements operating at 1.3mm and 0.87mm wavelength. In parallel, development of next-generation backend instrumentation, coupled with high throughput correlation architectures, will boost sensitivity, allowing the new stations to be of modest collecting area while still improving imaging fidelity and angular resolution. The goal of these efforts is to move from imaging static horizon scale structure to dynamic reconstructions that capture the processes of accretion and jet launching in near real time.

In April 2019, the Event Horizon Telescope (EHT) collaboration revealed the first image of the candidate super-massive black hole (SMBH) at the centre of the giant elliptical galaxy Messier 87 (M87). This event-horizon-scale image shows a ring of glowing plasma with a dark patch at the centre, which is interpreted as the shadow of the black hole. This breakthrough result, which represents a powerful confirmation of Einstein's theory of gravity, or general relativity, was made possible by assembling a global network of radio telescopes operating at millimetre wavelengths that for the first time included the Atacama Large Millimeter/ submillimeter Array (ALMA). The addition of ALMA as an anchor station has enabled a giant leap forward by increasing the sensitivity limits of the EHT by an order of magnitude, effectively turning it into an imaging array. The published image demonstrates that it is now possible to directly study the event horizon shadows of SMBHs via electromagnetic radiation, thereby transforming this elusive frontier from a mathematical concept into an astrophysical reality. The expansion of the array over the next few years will include new stations on different continents - and eventually satellites in space. This will provide progressively sharper and higher-fidelity images of SMBH candidates, and potentially even movies of the hot plasma orbiting around SMBHs. These improvements will shed light on the processes of black hole accretion and jet formation on event-horizon scales, thereby enabling more precise tests of general relativity in the truly strong field regime.

The Event Horizon Telescope (EHT) recently produced the first horizon-scale image of a supermassive black hole. Expanding the array to include a 3-meter space telescope operating at >200 GHz enables mass measurements of many black holes, movies of black hole accretion flows, and new tests of general relativity that are impossible from the ground.

The astrophysical agenda of the 21st century requires a very sharp view of celestial objects. High angular resolution studies are essential for fundamental studies of a broad variety of astrophysical phenomena ranging from relativistic physics of black holes, their gravitational and electromagnetic imprints, violent transient processes, including those producing detectable gravitational waves, birth and evolution of planetary systems. Over the past decades, radio astronomy made huge leap in achieving ground-breaking angular resolution measured in tens of microarcseconds (one tenth of nanoradian and better). Recently a global Event Horizon Telescope (EHT) collaboration obtained first direct images of the shadow of a super-massive black hole in the nucleus of the active galaxy M87. These observations were conducted at 230 GHz. The two first generation Space Very Long Baseline Interferometry (VLBI) missions, VSOP/HALCA led by the Japan Aerospace Exploration Agency (JAXA) and RadioAstron led by the Russia Roscosmos State Corporation and Russia Academy of Sciences, achieved the highest angular resolution at frequencies from 0.3 to 22 GHz in observations conducted in the period 1997 - 2019. The next step in advancing high angular resolution radio astronomy is in combining high frequency (millimeter and sub-millimeter wavelengths) and interferometric baselines exceeding the Earth diameter. The present THEZA White Paper describes a combination which would unify technology developments in giga-/tera-hertz instrumentation and space-borne radio astronomy. The current preprint version of the THEZA White Paper is slightly re-formatted and edited comparing to the official submitted version.

Recent developments in compact object astrophysics, especially the discovery of merging neutron stars by LIGO, the imaging of the black hole in M87 by the Event Horizon Telescope (EHT) and high precision astrometry of the Galactic Center at close to the event horizon scale by the GRAVITY experiment motivate the development of numerical source models that solve the equations of general relativistic magnetohydrodynamics (GRMHD). Here we compare GRMHD solutions for the evolution of a magnetized accretion flow where turbulence is promoted by the magnetorotational instability from a set of nine GRMHD codes: Athena++, BHAC, Cosmos++, ECHO, H-AMR, iharm3D, HARM-Noble, IllinoisGRMHD and KORAL. Agreement between the codes improves as resolution increases, as measured by a consistently applied, specially developed set of code performance metrics. We conclude that the community of GRMHD codes is mature, capable, and consistent on these test problems.

Upcoming VLBI observations will resolve nearby supermassive black holes, most notably Sagittarius A* and M87, on event horizon-scales. Recent observations of Sagittarius A* with the Event Horizon Telescope have revealed horizon-scale structure. Accordingly, the detection and measurement of the back hole ``shadow" is expected to enable the existence of astrophysical black holes to be verified directly. Although the theoretical description of the shadow is straightforward, its observational appearance is largely determined by the properties of the surrounding accretion flow, which is highly turbulent. We introduce a new polarised general-relativistic radiative transfer code, \texttt{BHOSS}, which accurately solves the equations of polarised radiative transfer in arbitrary strong-gravity environments, providing physically-realistic images of astrophysical black holes on event horizon-scales, as well as also providing insight into the fundamental properties and nature of the surrounding accretion flow environment.

Worldwide very long baseline radio interferometry arrays are expected to obtain horizon-scale images of supermassive black hole candidates as well as of relativistic jets in several nearby active galactic nuclei. This motivates the development of models for magnetohydrodynamic flows in strong gravitational fields. The Black Hole Accretion Code (BHAC) intends to aid with the modelling of such sources by means of general relativistic magnetohydrodynamical (GRMHD) simulations in arbitrary stationary spacetimes. New additions were required to guarantee an accurate evolution of the magnetic field when small and large scales are captured simultaneously. We discuss the adaptive mesh refinement (AMR) techniques employed in BHAC, essential to keep several problems computationally tractable, as well as staggered-mesh-based constrained transport (CT) algorithms to preserve the divergence-free constraint of the magnetic field, including a general class of prolongation operators for face-allocated variables compatible with them. Through several standard tests, we show that the choice of divergence-control method can produce qualitative differences in simulations of scientifically relevant accretion problems. We demonstrate the ability of AMR to reduce the computational costs of accretion simulations while sufficiently resolving turbulence from the magnetorotational instability. In particular, we describe a simulation of an accreting Kerr black hole in Cartesian coordinates using AMR to follow the propagation of a relativistic jet while self-consistently including the jet engine, a problem set up-for which the new AMR implementation is particularly advantageous. The CT methods and AMR strategies discussed here are being employed in the simulations performed with BHAC used in the generation of theoretical models for the Event Horizon Telescope Collaboration.

The galaxy M 87 harbors a kiloparsec-scale relativistic jet, whose origin coincides with a supermassive black hole. Observational mm-VLBI campaigns are capable of resolving the jet-launching region at the scale of the event horizon. In order to provide a context for interpreting these observations, realistic general-relativistic magnetohydrodynamical (GRMHD) models of the accretion flow are constructed. The characteristics of the observed spectral-energy distribution (SED) depend on the shape of the electrons' energy-distribution function (eDF). The dependency on the eDF is omitted in the modeling of the first Event Horizon Telescope results. In this work, we aim to model the M 87 SED from radio up to NIR/optical frequencies using a thermal-relativistic Maxwell- J\"uttner distribution, as well as a relativistic $\kappa$-distribution function. The electrons are injected based on sub-grid, particle-in-cell parametrizations for sub-relativistic reconnection. A GRMHD simulation in Cartesian-Kerr-Schild coordinates, using eight levels of adaptive mesh refinement (AMR), forms the basis of our model. To obtain spectra and images, the GRMHD data is post-processed with the ray-tracing code RAPTOR, which is capable of ray tracing through AMR GRMHD simulation data. We obtain radio spectra in both the thermal-jet and $\kappa$-jet models consistent with radio observations. Additionally, the $\kappa$-jet models also recover the NIR/optical emission. The models recover the observed source sizes and core shifts and obtain a jet power of $\approx 10^{43}$ ergs/s. In the $\kappa$-jet models, both the accretion rates and jet powers are approximately two times lower than the thermal-jet model. The frequency cut-off observed at $\nu \approx 10^{15}$ Hz is recovered when the accelerator size is $10^6$ - $10^8$ cm, this could potentially point to an upper limit for plasmoid sizes in the jet of M 87.

Radio-astronomical observations of the supermassive black-hole candidate in the galactic center will soon offer the possibility to study gravity in its strongest regimes and to test different models for these compact objects. Studies based on semi-analytic models and strong-field images of stationary plasma configurations around boson stars have stressed the difficulty to distinguish them from black holes. For the first time, we present general-relativistic magnetohydrodynamic simulations followed by consistent general-relativistic radiative-transfer calculations in a boson-star spacetime to revisit the appearance of these objects in a realistic accretion scenario. We find that the absence of an event horizon in a boson star leads to important differences in the dynamics of the accretion and results in both the formation of a small torus in the interior of the boson star and in the absence of an evacuated high-magnetization funnel in the polar regions. Synthetic reconstructed images considering realistic astronomical observing conditions show that differences in the appearance of the two compact objects are large enough to be detectable. These results, which also apply to other horizonless compact objects, strengthen confidence in the ability to determine the presence of an event horizon via radio observations and highlight the importance of self-consistent multidimensional simulations to study the compact object at the galactic center.

High-resolution Very-Long-Baseline Interferometry observations of NGC 1052 show a two sided jet with several regions of enhanced emission and a clear emission gap between the two jets.This gap shrinks with increasing frequency and vanishes around $\nu\sim43$ GHz. The observed structures are due to both the macroscopic fluid dynamics interacting with the surrounding ambient medium including an obscuring torus and the radiation microphysics. In this paper we investigate the possible physical conditions in relativistic jets of NGC 1052 by directly modelling the observed emission and spectra via state-of-the-art special-relativistic hydrodynamic (SRHD) simulations and radiative transfer calculations. To investigate the physical conditions in the relativistic jet we coupled our radiative transfer code to evolutionary algorithms and performed simultaneous modelling of the observed jet structure and the broadband radio spectrum. During the calculation of the radiation we consider both thermal and non-thermal emission. In order to compare our model to VLBI observations we take into account the sparse sampling of the u-v plane, the array properties and the imaging algorithm. We present for the first time an end-to-end pipeline for fitting numerical simulations to VLBI observations of relativistic jets taking into account the macrophysics including fluid dynamics and ambient medium configurations together with thermal/non-thermal emission and the properties of the observing array. The detailed analysis of our simulations shows that the structure and properties of the observed relativistic jets in NGC 1052 can be reconstructed by a slightly over-pressured jet ($d_k\sim1.5$) embedded in a decreasing pressure ambient medium

Particle acceleration induced by fast magnetic reconnection may help to solve current puzzles related to the interpretation of the very high energy (VHE) and neutrino missions from AGNs and compact sources in general. Our general relativistic-MHD simulations of accretion disk-corona systems reveal the growth of turbulence driven by MHD instabilities that lead to the development of fast magnetic reconnection in the corona. In addition, our simulations of relativistic MHD jets reveal the formation of several sites of fast reconnection induced by current-driven kink turbulence. The injection of thousands of test particles in these regions cause acceleration up to energies of several PeVs, thus demonstrating the ability of this process to accelerate particles and produce VHE and neutrino emission, specially in blazars. Finally, we discuss how reconnection can also explain the observed VHE luminosity-black hole mass correlation, involving hundreds of non-blazar sources like Perseus A, and black hole binaries.

Our Galactic Center, Sagittarius A* (Sgr A*), is believed to harbour a supermassive black hole (BH), as suggested by observations tracking individual orbiting stars. Upcoming sub-millimetre very-long-baseline-interferometry (VLBI) images of Sgr A* carried out by the Event-Horizon-Telescope Collaboration (EHTC) are expected to provide critical evidence for the existence of this supermassive BH. We assess our present ability to use EHTC images to determine if they correspond to a Kerr BH as predicted by Einstein's theory of general relativity (GR) or to a BH in alternative theories of gravity. To this end, we perform general-relativistic magnetohydrodynamical (GRMHD) simulations and use general-relativistic radiative transfer (GRRT) calculations to generate synthetic shadow images of a magnetised accretion flow onto a Kerr BH. In addition, and for the first time, we perform GRMHD simulations and GRRT calculations for a dilaton BH, which we take as a representative solution of an alternative theory of gravity. Adopting the VLBI configuration from the 2017 EHTC campaign, we find that it could be extremely difficult to distinguish between BHs from different theories of gravity, thus highlighting that great caution is needed when interpreting BH images as tests of GR.

To date, the most precise tests of general relativity have been achieved through pulsar timing, albeit in the weak-field regime. Since pulsars are some of the most precise and stable "clocks" in the Universe, present observational efforts are focused on detecting pulsars in the vicinity of supermassive black holes (most notably in our Galactic Centre), enabling pulsar timing to be used as an extremely precise probe of strong-field gravity. In this paper a mathematical framework to describe test-particle dynamics in general black hole spacetimes is presented, and subsequently used to study a binary system comprising a pulsar orbiting a black hole. In particular, taking into account the parameterization of a general spherically symmetric black hole metric, general analytic expressions for both the advance of the periastron and for the orbital period of a massive test particle are derived. Furthermore, these expressions are applied to four representative cases of solutions arising in both general relativity and in alternative theories of gravity. Finally, this framework is applied to the Galactic Centre $S$-stars and four distinct pulsar toy models. It is shown that by adopting a fully general-relativistic description of test-particle motion which is independent of any particular theory of gravity, observations of pulsars can help impose better constraints on alternative theories of gravity than is presently possible.

In this study we investigate jet interaction at a microscopic level in a cosmological environment, which responds to a key open question in the study of relativistic jets. Using small simulation systems during prior research, we initially studied the evolution of both electron-proton and electron-positron relativistic jets containing helical magnetic fields, by focusing on their interactions with an ambient plasma. Here, using larger jet radii, we have performed simulations of global jets containing helical magnetic fields in order to examine how helical magnetic fields affect kinetic instabilities such as the Weibel instability, the kinetic Kelvin-Helmholtz instability (kKHI) and the Mushroom instability (MI). We found that the evolution of global jets strongly depends on the size of the jet radius. For example, phase bunching of jet electrons, in particular in the electron-proton jet, is mixed with larger jet radius due to the more complicated structures of magnetic fields with excited kinetic instabilities. In our simulation study these kinetic instabilities lead to new types of instabilities in global jets. In the electron-proton jet simulation a modified recollimation occurs and jet electrons are strongly perturbed. In the electron-positron jet simulation mixed kinetic instabilities occur at early times followed by a turbulence-like structure. Simulations using much larger (and longer) systems are further required in order to thoroughly investigate the evolution of global jets containing helical magnetic fields.

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.

We present the black hole accretion code (BHAC), a new multidimensional general-relativistic magnetohydrodynamics module for the MPI-AMRVAC framework. BHAC has been designed to solve the equations of ideal general-relativistic magnetohydrodynamics in arbitrary spacetimes and exploits adaptive mesh refinement techniques with an efficient block-based approach. Several spacetimes have already been implemented and tested. We demonstrate the validity of BHAC by means of various one-, two-, and three-dimensional test problems, as well as through a close comparison with the HARM3D code in the case of a torus accreting onto a black hole. The convergence of a turbulent accretion scenario is investigated with several diagnostics and we find accretion rates and horizon-penetrating fluxes to be convergent to within a few percent when the problem is run in three dimensions. Our analysis also involves the study of the corresponding thermal synchrotron emission, which is performed by means of a new general-relativistic radiative transfer code, BHOSS. The resulting synthetic intensity maps of accretion onto black holes are found to be convergent with increasing resolution and are anticipated to play a crucial role in the interpretation of horizon-scale images resulting from upcoming radio observations of the source at the Galactic Center.

High-resolution Very-Long-Baseline Interferometry observations of active galactic nuclei have revealed asymmetric structures in the jets of radio galaxies. These asymmetric structures may be due to internal asymmetries in the jet, could be induced by the different conditions in the surrounding ambient medium including the obscuring torus, or a combination of the two. In this paper we investigate the influence of the ambient medium (including the obscuring torus) on the observed properties of jets from radio galaxies. We performed special-relativistic hydrodynamic (RHD) simulations of over-pressured and pressure-matched jets using the special-relativistic hydrodynamics code \texttt{Ratpenat}, which is based on a second-order accurate finite-volume method and an approximate Riemann solver. Using a newly developed emission code to compute the electromagnetic emission, we have investigated the influence of different ambient medium and torus configurations on the jet structure and subsequently computed the non-thermal emission produced by the jet and the thermal absorption due to the torus. To better compare the emission simulations with observations we produced synthetic radio maps, taking into account the properties of the observatory. The detailed analysis of our simulations shows that the observed asymmetries can be produced by the interaction of the jet with the ambient medium and by the absorption properties of the obscuring torus.

Fragmentation of filaments into dense cores is thought to be an important step in forming stars. The bar-mode instability of spherically collapsing cores found in previous linear analysis invokes a possibility of re-fragmentation of the cores due to their ellipsoidal (prolate or oblate) deformation. To investigate this possibility, here we perform three-dimensional self-gravitational hydrodynamics simulations that follow all the way from filament fragmentation to subsequent core collapse. We assume the gas is polytropic with index \gamma, which determines the stability of the bar-mode. For the case that the fragmentation of isolated hydrostatic filaments is triggered by the most unstable fragmentation mode, we find the bar mode grows as collapse proceeds if \gamma < 1.1, in agreement with the linear analysis. However, it takes more than ten orders-of-magnitude increase in the central density for the distortion to become non-linear. In addition to this fiducial case, we also study non-fiducial ones such as the fragmentation is triggered by a fragmentation mode with a longer wavelength and it occurs during radial collapse of filaments and find the distortion rapidly grows. In most of astrophysical applications, the effective polytropic index of collapsing gas exceeds 1.1 before ten orders-of-magnitude increase in the central density. Thus, supposing the fiducial case of filament fragmentation, re-fragmentation of dense cores would not be likely and their final mass would be determined when the filaments fragment.

Einstein's General Theory of Relativity (GR) successfully describes gravity. The most fundamental predictions of GR are black holes (BHs), but in spite of many convincing BH candidates in the Universe, there is no conclusive experimental proof of their existence using astronomical observations in the electromagnetic spectrum. Are BHs real astrophysical objects? Does GR hold in its most extreme limit or are alternatives needed? The prime target to address these fundamental questions is in the center of our own Galaxy, which hosts the closest and best-constrained supermassive BH candidate in the Universe, Sagittarius A* (Sgr A*). Three different types of experiments hold the promise to test GR in a strong-field regime using observations of Sgr A* with new-generation instruments. The first experiment aims to image the relativistic plasma emission which surrounds the event horizon and forms a "shadow" cast against the background, whose predicted size (~50 microarcseconds) can now be resolved by upcoming VLBI experiments at mm-waves such as the Event Horizon Telescope (EHT). The second experiment aims to monitor stars orbiting Sgr A* with the upcoming near-infrared interferometer GRAVITY at the Very Large Telescope (VLT). The third experiment aims to time a radio pulsar in tight orbit about Sgr A* using radio telescopes (including the Atacama Large Millimeter Array or ALMA). The BlackHoleCam project exploits the synergy between these three different techniques and aims to measure the main BH parameters with sufficient precision to provide fundamental tests of GR and probe the spacetime around a BH in any metric theory of gravity. Here, we review our current knowledge of the physical properties of Sgr A* as well as the current status of such experimental efforts towards imaging the event horizon, measuring stellar orbits, and timing pulsars around Sgr A*.

In order to discuss the potential impact of solar 'superflares' on space weather, we investigated statistical relations among energetic proton peak flux with energy higher than $ 10 \rm MeV$ ($F_p$), CME speed near the Sun ($V_{CME}$) obtained by {\it SOHO}/LASCO coronagraph and flare soft X-ray peak flux in 1-8\AA band ($F_{SXR}$) during 110 major solar proton events (SPEs) recorded from 1996 to 2014. The linear regression fit results in the scaling relations $V_{CME} \propto F_{SXR}^\alpha$, $F_p\propto F_{SXR}^\beta$ and $F_p\propto V_{CME}^\gamma$ with $\alpha = 0.30\pm 0.04$, $\beta = 1.19 \pm 0.08$ and $\gamma = 4.35 \pm 0.50$, respectively. On the basis of simple physical assumptions, on the other hand, we derive scaling relations expressing CME mass ($M_{CME}$), CME speed and energetic proton flux in terms of total flare energy ($E_{flare}$) as, $M_{CME}\propto E_{flare}^{2/3}$, $V_{CME}\propto E_{flare}^{1/6}$ and $F_{p}\propto E_{flare}^{5/6}\propto V_{CME}^5$, respectively. We then combine the derived scaling relations with observation, and estimated the upper limit of $V_{CME}$ and $F_p$ to be associated with possible solar superflares.

(Abridged) We here continue our effort to model the behaviour of matter when orbiting or accreting onto a generic black hole by developing a new numerical code employing advanced techniques geared solve the equations of in general-relativistic hydrodynamics. The new code employs a number of high-resolution shock-capturing Riemann-solvers and reconstruction algorithms, exploiting the enhanced accuracy and the reduced computational cost of AMR techniques. In addition, the code makes use of sophisticated ray-tracing libraries that, coupled with general-relativistic radiation-transfer calculations, allow us to compute accurately the electromagnetic emissions from such accretion flows. We validate the new code by presenting an extensive series of stationary accretion flows either in spherical or axial symmetry and performed either in 2D or 3D. In addition, we consider the highly nonlinear scenario of a recoiling black hole produced in the merger of a supermassive black hole binary interacting with the surrounding circumbinary disc. In this way we can present, for the first time, ray-traced images of the shocked fluid and the light-curve resulting from consistent general-relativistic radiation-transport calculations from this process. The work presented here lays the ground for the development of a generic computational infrastructure employing AMR techniques to deal accurately and self-consistently with accretion flows onto compact objects. In addition to the accurate handling of the matter, we provide a self-consistent electromagnetic emission from these scenarios by solving the associated radiative-transfer problem. While magnetic fields are presently excluded from our analysis, the tools presented here can have a number of applications to study accretion flows onto black holes or neutron stars.

We study the interaction of relativistic jets with their environment, using 3-dimensional relativistic particle-in-cell simulations for two cases of jet composition: (i) electron-proton ($e^{-}-p^{+}$) and (ii) electron-positron ($e^{\pm}$) plasmas containing helical magnetic fields. We have performed simulations of "global" jets containing helical magnetic fields in order to examine how helical magnetic fields affect kinetic instabilities such as the Weibel instability, the kinetic Kelvin-Helmholtz instability and the Mushroom instability. We have found that these kinetic instabilities are suppressed and new types of instabilities can grow. For the $e^{-}-p^{+}$ jet, a recollimation-like instability occurs and jet electrons are strongly perturbed, whereas for the $e^{\pm}$ jet, a recollimation-like instability occurs at early times followed by kinetic instability and the general structure is similar to a simulation without a helical magnetic field. We plan to perform further simulations using much larger systems to confirm these new findings.

Collaborative international efforts under the name of the Event Horizon Telescope project, using sub- mm very long baseline interferometry, are soon expected to provide the first images of the shadow cast by the candidate supermassive black hole in our Galactic center, Sagittarius A*. Observations of this shadow would provide direct evidence of the existence of astrophysical black holes. Although it is expected that astrophysical black holes are described by the axisymmetric Kerr solution, there also exist many other black hole solutions, both in general relativity and in other theories of gravity, which cannot presently be ruled out. To this end, we present calculations of black hole shadow images from various metric theories of gravity as described by our recent work on a general parameterisation of axisymmetric black holes [R. Konoplya, L. Rezzolla and A. Zhidenko, Phys. Rev. D 93, 064015 (2016)]. An algorithm to perform general ray-tracing calculations for any metric theory of gravity is first outlined and then employed to demonstrate that even for extremal metric deformation parameters of various black hole spacetimes, this parameterisation is both robust and rapidly convergent to the correct solution.

In the study of relativistic jets one of the key open questions is their interaction with the environment on the microscopic level. Here, we study the initial evolution of both electron$-$proton ($e^{-}-p^{+}$) and electron$-$positron ($e^{\pm}$) relativistic jets containing helical magnetic fields, focusing on their interaction with an ambient plasma. We have performed simulations of "global" jets containing helical magnetic fields in order to examine how helical magnetic fields affect kinetic instabilities such as the Weibel instability, the kinetic Kelvin-Helmholtz instability (kKHI) and the Mushroom instability (MI). In our initial simulation study these kinetic instabilities are suppressed and new types of instabilities can grow. In the $e^{-}-p^{+}$ jet simulation a recollimation-like instability occurs and jet electrons are strongly perturbed. In the $e^{\pm}$ jet simulation a recollimation-like instability occurs at early times followed by a kinetic instability and the general structure is similar to a simulation without helical magnetic field. Simulations using much larger systems are required in order to thoroughly follow the evolution of global jets containing helical magnetic fields.

Using the three-dimensional relativistic magnetohydrodynamic code RAISHIN, we investigated the influence of radial density profile on the spatial development of the current-driven kink instability along magnetized rotating, relativistic jets. For the purpose of our study, we used a non-periodic computational box, the jet flow is initially established across the computational grid, and a precessional perturbation at the inlet triggers the growth of the kink instability. We studied light as well as heavy jets with respect to the environment depending on the density profile. Different angular velocity amplitudes have been also tested. The results show the propagation of a helically kinked structure along the jet and relatively stable configuration for the lighter jets. The jets appear to be collimated by the magnetic field and the flow is accelerated due to conversion of electromagnetic into kinetic energy. We also identify regions of high current density in filamentary current sheets, indicative of magnetic reconnection, which are associated to the kink unstable regions and correlated to the decrease of the sigma parameter of the flow. We discuss the implications of our findings for Poynting-flux dominated jets in connection with magnetic reconnection process. We find that fast magnetic reconnection may be driven by the kink-instability turbulence and govern the transformation of magnetic into kinetic energy thus providing an efficient way to power and accelerate particles in AGN and gamma-ray-burst relativistic jets.

Fast Magnetic Reconnection is currently regarded as an important process also beyond the solar system, specially in magnetically dominated regions of galactic and extragalactic sources like the surrounds of black holes and relativistic jets. In this lecture we discuss briefly the theory of fast magnetic reconnection, specially when driven by turbulence which is very frequent in Astrophysical flows, and its implications for relativistic particle acceleration. Then we discuss these processes in the context of the sources above, showing recent analytical and multidimensional numerical MHD studies that indicate that fast reconnection can be a powerful process to accelerate particles to relativistic velocities, produce the associated high energy non-thermal emission, and account for efficient conversion of magnetic into kinetic energy in these flows.

One of the key open questions in the study of relativistic jets is their interaction with the environment. Here, we study the initial evolution of both electron-proton and electron-positron relativistic jets, focusing on their lateral interaction with the ambient plasma. We trace the generation and evolution of the toroidal magnetic field generated by both kinetic Kelvin-Helmholtz (kKH) and Mushroom instabilities (MI). This magnetic field collimates the jet. We show that in electron-proton jet, electrons are perpendicularly accelerated with jet collimation. The magnetic polarity switches from the clockwise to anti-clockwise in the middle of jet, as the instabilities weaken. For the electron-positron jet, we find strong mixture of electron-positron with the ambient plasma, that results in the creation of a bow shock. Merger of magnetic field current filaments generate density bumps which initiate a forward shock. The strong mixing between jet and ambient particles prevents full development of the jet on the studied scale. Our results therefore provide a direct evidence for both jet collimation and particle acceleration in the created bow shock. Differences in the magnetic field structures generated by electron-proton and electron-positron jets may contribute to observable differences in the polarized properties of emission by electrons.

We present the first polarimetric space VLBI imaging observations at 22 GHz. BL Lacertae was observed in 2013 November 10 with the RadioAstron space VLBI mission, including a ground array of 15 radio telescopes. The instrumental polarization of the space radio telescope is found to be within 9%, demonstrating the polarimetric imaging capabilities of RadioAstron at 22 GHz. Ground-space fringes were obtained up to a projected baseline distance of 7.9 Earth's diameters in length, allowing us to image the jet in BL Lacertae with a maximum angular resolution of 21 $\mu$as, the highest achieved to date. We find evidence for emission upstream of the radio core, which may correspond to a recollimation shock at about 40 $\mu$as from the jet apex, in a pattern that includes other recollimation shocks at approximately 100 $\mu$as and 250 $\mu$as from the jet apex. Polarized emission is detected in two components within the innermost 0.5 mas from the core, as well as in some knots 3 mas downstream. Faraday rotation analysis, obtained from combining RadioAstron 22 GHz and ground-based 15 GHz and 43 GHz images, shows a gradient in rotation measure and Faraday corrected polarization vector as a function of position angle with respect to the core, suggesting that the jet in BL Lacertae is threaded by a helical magnetic field. The intrinsic de-boosted brightness temperature in the unresolved core exceeds $3\!\times\!10^{12}$ K, suggesting at the very least departure from equipartition of energy between the magnetic field and radiating particles.

Percolation theory is applied to the phase-transition dynamics of domain pattern formation in segregating binary Bose--Einstein condensates in quasi-two-dimensional systems. Our finite-size-scaling analysis shows that the percolation threshold of the initial domain pattern emerging from the dynamic instability is close to 0.5 for strongly repulsive condensates. The percolation probability is universally described with a scaling function when the probability is rescaled by the characteristic domain size in the dynamic scaling regime of the phase-ordering kinetics, independent of the intercomponent interaction. It is revealed that an infinite domain wall sandwiched between percolating domains in the two condensates has an noninteger fractal dimension and keeps the scaling behavior during the dynamic scaling regime. This result seems to be in contrast to the argument that the dynamic scale invariance is violated in the presence of an infinite topological defect in numerical cosmology.

We have performed two-dimensional special-relativistic magnetohydrodynamic simulations of non-equilibrium over-pressured relativistic jets in cylindrical geometry. Multiple stationary recollimation shock and rarefaction structures are produced along the jet by the nonlinear interaction of shocks and rarefaction waves excited at the interface between the jet and the surrounding ambient medium. Although initially the jet is kinematically dominated, we have considered axial, toroidal and helical magnetic fields to investigate the effects of different magnetic-field topologies and strengths on the recollimation structures. We find that an axial field introduces a larger effective gas-pressure and leads to stronger recollimation shocks and rarefactions, resulting in larger flow variations. The jet boost grows quadratically with the initial magnetic field. On the other hand, a toroidal field leads to weaker recollimation shocks and rarefactions, modifying significantly the jet structure after the first recollimation rarefaction and shock. The jet boost decreases systematically. For a helical field, instead, the behaviour depends on the magnetic pitch, with a phenomenology that ranges between the one seen for axial and toroidal magnetic fields, respectively. In general, however, a helical magnetic field yields a more complex shock and rarefaction substructure close to the inlet that significantly modifies the jet structure. The differences in shock structure resulting from different field configurations and strengths may have observable consequences for disturbances propagating through a stationary recollimation shock.

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.

HH 212 is a nearby (400 pc) highly collimated protostellar jet powered by a Class 0 source in Orion. We have mapped the inner 80" (~ 0.16 pc) of the jet in SiO (J=8-7) and CO (J=3-2) simultaneously at ~ 0.5 resolution with the Atacama Millimeter/Submillimeter Array at unprecedented sensitivity. The jet consists of a chain of knots, bow shocks, and sinuous structures in between. As compared to that seen in our previous observations with the Submillimeter Array, it appears to be more continuous, especially in the northern part. Some of the knots are now seen associated with small bow shocks, with their bow wings curving back to the jet axis, as seen in pulsed jet simulations. Two of them are reasonably resolved, showing kinematics consistent with sideways ejection, possibly tracing the internal working surfaces formed by a temporal variation in the jet velocity. In addition, nested shells are seen in CO around the jet axis connecting to the knots and bow shocks, driven by them. The proper motion of the jet is estimated to be ~ 115+-50 km/s, comparing to our previous observations. The jet has a small semi-periodical wiggle, with a period of ~ 93 yrs. The amplitude of the wiggle first increases with the distance from the central source and then stays roughly constant. One possible origin of the wiggle could be the kink instability in a magnetized jet.

Relativistic jets in active galactic nuclei (AGN) are among the most powerful astrophysical objects discovered to date. Indeed, jetted AGN studies have been considered a prominent science case for SKA, and were included in several different chapters of the previous SKA Science Book (Carilli & Rawlings 2004). Most of the fundamental questions about the physics of relativistic jets still remain unanswered, and await high-sensitivity radio instruments such as SKA to solve them. These questions will be addressed specially through analysis of the massive data sets arising from the deep, all-sky surveys (both total and polarimetric flux) from SKA1. Wide-field very-long-baseline-interferometric survey observations involving SKA1 will serve as a unique tool for distinguishing between extragalactic relativistic jets and star forming galaxies via brightness temperature measurements. Subsequent SKA1 studies of relativistic jets at different resolutions will allow for unprecedented cosmological studies of AGN jets up to the epoch of re-ionization, enabling detailed characterization of the jet composition, magnetic field, particle populations, and plasma properties on all scales. SKA will enable us to study the dependence of jet power and star formation on other properties of the AGN system. SKA1 will enable such studies for large samples of jets, while VLBI observations involving SKA1 will provide the sensitivity for pc-scale imaging, and SKA2 (with its extraordinary sensitivity and dynamic range) will allow us for the first time to resolve and model the weakest radio structures in the most powerful radio-loud AGN.

We have investigated particle acceleration and emission from shocks and shear flows associated with an unmagnetized relativistic jet plasma propagating into an unmagnetized ambient plasma. Strong electro-magnetic fields are generated in the jet shock via the filamentation (Weibel) instability. Shock field strength and structure depend on plasma composition (($e^{\pm}$ or $e^-$- $p^+$ plasmas) and Lorentz factor. In the velocity shear between jet and ambient plasmas, strong AC ($e^{\pm}$ plasmas) or DC ($e^-$- $p^+$ plasmas) magnetic fields are generated via the kinetic Kelvin-Helmholtz instability (kKHI), and the magnetic field structure also depends on the jet Lorentz factor. We have calculated, self-consistently, the radiation from electrons accelerated in shock generated magnetic fields. The spectra depend on the jet's initial Lorentz factor and temperature via the resulting particle acceleration and magnetic field generation. Our ongoing "Global" jet simulations containing shocks and velocity shears will provide us with the ability to calculate and model the complex time evolution and/or spectral structure observed from gamma-ray bursts, AGN jets, and supernova remnants.

We have investigated magnetic field generation in velocity shears via the kinetic Kelvin-Helmholtz instability (kKHI) using a relativistic plasma jet core and stationary plasma sheath. Our three-dimensional particle-in-cell simulations consider plasma jet cores with Lorentz factors of 1.5, 5, and 15 for both electron-proton and electron-positron plasmas. For electron-proton plasmas we find generation of strong large-scale DC currents and magnetic fields which extend over the entire shear-surface and reach thicknesses of a few tens of electron skin depths. For electron-positron plasmas we find generation of alternating currents and magnetic fields. Jet and sheath plasmas are accelerated across the shear surface in the strong magnetic fields generated by the kKHI. The mixing of jet and sheath plasmas generates transverse structure similar to that produced by the Weibel instability.

We have investigated the influence of velocity shear and a radial density profile on the spatial development of the current driven kink instability along helically magnetized relativistic jets via three-dimensional relativistic magnetohydrodynamic simulations. In this study, we use a non-periodic computational box, the jet flow is initially established across the computational grid, and a precessional perturbation at the inlet triggers growth of the kink instability. If the velocity shear radius is located inside the characteristic radius of the helical magnetic field, a static non-propagating current driven kink is excited as the perturbation propagates down the jet. Temporal growth disrupts the initial flow across the computational grid not too far from the inlet. On the other hand, if the velocity shear radius is outside the characteristic radius of the helical magnetic field, the kink is advected with the flow and grows spatially down the jet. In this case flow is maintained to much larger distances from the inlet. The effect of different radial density profiles is more subtle. When the density increases with radius, the kink appears to saturate by the end of the simulation without apparent disruption of the helical twist. This behavior suggests that relativistic jets consisting of a tenuous spine surrounded by a denser medium with a velocity shear radius outside the radius of maximum toroidal magnetic field have a relatively stable configuration.

We have investigated via two-dimensional relativistic MHD simulations the long-term evolution of turbulence created by a relativistic shock propagating through an inhomogeneous medium. In the postshock region, magnetic field is strongly amplified by turbulent motions triggered by preshock density inhomogeneities. Using a long-simulation box we have followed the magnetic-field amplification until it is fully developed and saturated. The turbulent velocity is sub-relativistic even for a strong shock. Magnetic-field amplification is controled by the turbulent motion and saturation occurs when the magnetic energy is comparable to the turbulent kinetic energy. Magnetic-field amplification and saturation depend on the initial strength and direction of the magnetic field in the preshock medium, and on the shock strength. If the initial magnetic field is perpendicular to the shock normal, the magnetic field is first compressed at the shock and then can be amplified by turbulent motion in the postshock region. Saturation occurs when the magnetic energy becomes comparable to the turbulent kinetic energy in the postshock region. If the initial magnetic field in the preshock medium is strong, the postshock region becomes turbulent but significant field amplification does not occur. If the magnetic energy after shock compression is larger than the turbulent kinetic energy in the postshock region, significant field amplification does not occur. We discuss possible applications of our results to gamma-ray bursts and active galactic nuclei.

We have investigated generation of magnetic fields associated with velocity shear between an unmagnetized relativistic jet and an unmagnetized sheath plasma. We have examined the strong magnetic fields generated by kinetic shear (Kelvin-Helmholtz) instabilities. Compared to the previous studies using counter-streaming performed by Alves et al. (2012), the structure of KKHI of our jet-sheath configuration is slightly different even for the global evolution of the strong transverse magnetic field. In our simulations the major components of growing modes are the electric field $E_{\rm z}$ and the magnetic field $B_{\rm y}$. After the $B_{\rm y}$ component is excited, an induced electric field $E_{\rm x}$ becomes significant. However, other field components remain small. We find that the structure and growth rate of KKHI with mass ratios $m_{\rm i}/m_{\rm e} = 1836$ and $m_{\rm i}/m_{\rm e} = 20$ are similar. In our simulations saturation in the nonlinear stage is not as clear as in counter-streaming cases. The growth rate for a mildly-relativistic jet case ($\gamma_{\rm j} = 1.5$) is larger than for a relativistic jet case ($\gamma_{\rm j} = 15$).

We have investigated particle acceleration and shock structure associated with an unmagnetized relativistic jet propagating into an unmagnetized plasma. Strong magnetic fields generated in the trailing jet shock lead to transverse deflection and acceleration of the electrons. We have self-consistently calculated the radiation from the electrons accelerated in the turbulent magnetic fields. We find that the synthetic spectra depend on the bulk Lorentz factor of the jet, the jet temperature, and the strength of the magnetic fields generated in the shock. We have also begun study of electron acceleration in the strong magnetic fields generated by kinetic shear (Kelvin-Helmholtz) instabilities. Our calculated spectra should lead to a better understanding of the complex time evolution and/or spectral structure from gamma-ray bursts, relativistic jets, and supernova remnants.

We have investigated the role of the equation of state in resistive relativistic magnetohydrodynamics using a newly developed resistive relativistic magnetohydrodynamic code. A number of numerical tests in one-dimension and multi-dimensions are carried out in order to check the robustness and accuracy of the new code. The code passes all the tests in situations involving both small and large uniform conductivities. Equations of state which closely approximate the single-component perfect relativistic gas are introduced. Results from selected numerical tests using different equations of state are compared. The main conclusion is that the choice of the equation of state as well as the value of the electric conductivity can result in considerable dynamical differences in simulations involving shocks, instabilities, and magnetic reconnection.

We investigate the extraction of the rotational energy of a black hole under different accreting environment. When the accretion rate is moderate, the accretion disk consists of an outer thin disk and an inner advection-dominated accretion flow. In such a combined disk, the outer thin disk can sustain a magnetic field with moderate strength at the event horizon, leading to the formation of relativistic jets with moderate luminosity and speed via the magnetohrodynamic Penrose process. When the accretion rate increases enough, on the other hand, the disk becomes geometrically thin near the horizon. In this slim disk, the denser plasmas can sustain a stronger magnetic field than that in a combined disk, leading to the formation of jets with greater luminosity and speed via the Blandford-Znajek processs. It is discussed that the former jets are associated with the Fanaroff-Riley (FR) I galaxies and the latter with FR II galaxies.

We have investigated the influence of jet rotation and differential motion on the linear and nonlinear development of the current-driven (CD) kink instability of force-free helical magnetic equilibria via three-dimensional relativistic magnetohydrodynamic simulations. In this study, we follow the temporal development within a periodic computational box. Displacement of the initial helical magnetic field leads to the growth of the CD kink instability. We find that, in accord with linear stability theory, the development of the instability depends on the lateral distribution of the poloidal magnetic field. If the poloidal field significantly decreases outwards from the axis, the initial small perturbations grow strongly, and if multiple wavelengths are excited non-linear interaction eventually disrupts the initial cylindrical configuration. When the profile of the poloidal field is shallow, the instability develops slowly and eventually saturates. We briefly discuss implications of our findings for Poynting dominated jets.

lasma instabilities excited in collisionless shocks are responsible for particle acceleration. We have investigated the particle acceleration and shock structure associated with an unmagnetized relativistic electron-positron jet propagating into an unmagnetized electron-positron plasma. Cold jet electrons are thermalized and slowed while the ambient electrons are swept up to create a partially developed hydrodynamic-like shock structure. In the leading shock, electron density increases by a factor of about 3.5 in the simulation frame. Strong electromagnetic fields are generated in the trailing shock and provide an emission site. These magnetic fields contribute to the electron's transverse deflection behind the shock. Our initial results of a jet-ambient interaction with anti-parallel magnetic fields show pile-up of magnetic fields at the colliding shock, which may lead to reconnection and associated particle acceleration. We will investigate the radiation in transient stage as a possible generation mechanism of precursors of prompt emission. In our simulations we calculate the radiation from electrons in the shock region. The detailed properties of this radiation are important for understanding the complex time evolution and spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.

We present the properties of an extensive sample of molecular clouds in the Large Magellanic Cloud (LMC) mapped at 11 pc resolution in the CO(1-0) line. We identify clouds as regions of connected CO emission, and find that the distributions of cloud sizes, fluxes and masses are sensitive to the choice of decomposition parameters. In all cases, however, the luminosity function of CO clouds is steeper than dN/dL \propto L^{-2}, suggesting that a substantial fraction of mass is in low-mass clouds. A correlation between size and linewidth, while apparent for the largest emission structures, breaks down when those structures are decomposed into smaller structures. We argue that the correlation between virial mass and CO luminosity is the result of comparing two covariant quantities, with the correlation appearing tighter on larger scales where a size-linewidth relation holds. The virial parameter (the ratio of a cloud's kinetic to self-gravitational energy) shows a wide range of values and exhibits no clear trends with the CO luminosity or the likelihood of hosting young stellar object (YSO) candidates, casting further doubt on the assumption of virialization for molecular clouds in the LMC. Higher CO luminosity increases the likelihood of a cloud harboring a YSO candidate, and more luminous YSOs are more likely to be coincident with detectable CO emission, confirming the close link between giant molecular clouds and massive star formation.

We have investigated the influence of a velocity shear surface on the linear and non-linear development of the CD kink instability of force-free helical magnetic equilibria in 3D. In this study we follow the temporal development within a periodic computational box and concentrate on flows that are sub-Alfvenic on the cylindrical jet's axis. Displacement of the initial force-free helical magnetic field leads to the growth of CD kink instability. We find that helically distorted density structure propagates along the jet with speed and flow structure dependent on the radius of the velocity shear surface relative to the characteristic radius of the helically twisted force-free magnetic field. At small velocity shear surface radius the plasma flows through the kink with minimal kink propagation speed. The kink propagation speed increases as the velocity shear radius increases and the kink becomes more embedded in the plasma flow. A decreasing magnetic pitch profile and faster flow enhance the influence of velocity shear. Simulations show continuous transverse growth in the nonlinear phase of the instability. The growth rate of the CD kink instability and the nonlinear behavior also depend on the velocity shear surface radius and flow speed, and the magnetic pitch radial profile. Larger velocity shear radius leads to slower linear growth, makes a later transition to the nonlinear stage, and with larger maximum amplitude than occur for a static plasma column. However, when the velocity shear radius is much greater than the characteristic radius of the helical magnetic field, linear and non-linear development can be similar to the development of a static plasma column.

Using our new 3-D relativistic electromagnetic particle (REMP) code parallelized with MPI, we investigated long-term particle acceleration associated with a relativistic electron-positron jet propagating in an unmagnetized ambient electron-positron plasma. The simulations were performed using a much longer simulation system than our previous simulations in order to investigate the full nonlinear stage of the Weibel instability and its particle acceleration mechanism. Cold jet electrons are thermalized and ambient electrons are accelerated in the resulting shocks. Acceleration of ambient electrons leads to a maximum ambient electron density three times larger than the original value. Behind the bow shock in the jet shock strong electromagnetic fields are generated. These fields may lead to time dependent afterglow emission. We calculated radiation from electrons propagating in a uniform parallel magnetic field to verify the technique. We also used the new technique to calculate emission from electrons based on simulations with a small system. We obtained spectra which are consistent with those generated from electrons propagating in turbulent magnetic fields with red noise. This turbulent magnetic field is similar to the magnetic field generated at an early nonlinear stage of the Weibel instability. A fully developed shock within a larger system generates a jitter/synchrotron spectrum.

In order to precisely determine temperature and density of molecular gas in the Large Magellanic Cloud, we made observations of optically thin $^{13}$CO($J=3-2$) transition by using the ASTE 10m telescope toward 9 peaks where $^{12}$CO($J=3-2$) clumps were previously detected with the same telescope. The molecular clumps include those in giant molecular cloud (GMC) Types I (with no signs of massive star formation), II (with HII regions only), and III (with HII regions and young star clusters). We detected $^{13}$CO($J=3-2$) emission toward all the peaks and found that their intensities are 3 -- 12 times lower than those of $^{12}$CO($J=3-2$). We determined the intensity ratios of $^{12}$CO($J=3-2$) to $^{13}$CO($J=3-2$), $R^{12/13}_{3-2}$, and $^{13}$CO($J=3-2$) to $^{13}$CO($J=1-0$), $R^{13}_{3-2/1-0}$, at 45$\arcsec$ resolution. These ratios were used for radiative transfer calculations in order to estimate temperature and density of the clumps. The parameters of these clumps range kinetic temperature $T\mathrm{_{kin}}$ = 15 -- 200 K, and molecular hydrogen gas density $n(\mathrm{H_2})$ = 8$\times 10^2$ -- 7$\times 10^3$ cm$^{-3}$. We confirmed that the higher density clumps show higher kinetic temperature and that the lower density clumps lower kinetic temperature at a better accuracy than in the previous work. The kinetic temperature and density increase generally from a Type I GMC to a Type III GMC. We interpret that this difference reflects an evolutionary trend of star formation in molecular clumps. The $R^{13}_{3-2/1-0}$ and kinetic temperature of the clumps are well correlated with H$\alpha$ flux, suggesting that the heating of molecular gas $n(\mathrm{H_2})$ = $10^3$ -- $10^4$ cm$^{-3}$ can be explained by stellar FUV photons.

We have investigated the relaxation of a hydrostatic hot plasma column containing toroidal magnetic field by the Current-Driven (CD) kink instability as a model of pulsar wind nebulae. In our simulations the CD kink instability is excited by a small initial velocity perturbation and develops turbulent structure inside the hot plasma column. We demonstrate that, as envisioned by Begelman, the hoop stress declines and the initial gas pressure excess near the axis decreases. The magnetization parameter \sigma, the ratio of the Poynting to the kinetic energy flux, declines from an initial value of 0.3 to about 0.01 when the CD kink instability saturates. Our simulations demonstrate that axisymmetric models strongly overestimate the elongation of the pulsar wind nebulae. Therefore, the previous requirement for an extremely low pulsar wind magnetization can be abandoned. The observed structure of the pulsar wind nebulae do not contradict the natural assumption that the magnetic energy flux still remains a good fraction of the total energy flux after dissipation of alternating fields.

We perform two-dimensional relativistic magnetohydrodynamic simulations of a mildly relativistic shock propagating through an inhomogeneous medium. We show that the postshock region becomes turbulent owing to preshock density inhomogeneity, and the magnetic field is strongly amplified due to the stretching and folding of field lines in the turbulent velocity field. The amplified magnetic field evolves into a filamentary structure in two-dimensional simulations. The magnetic energy spectrum is flatter than the Kolmogorov spectrum and indicates that a so-called small-scale dynamo is occurring in the postshock region. We also find that the amount of magnetic-field amplification depends on the direction of the mean preshock magnetic field, and the time scale of magnetic-field growth depends on the shock strength.

The Magellanic Mopra Assessment (MAGMA) is a high angular resolution CO mapping survey of giant molecular clouds (GMCs) in the Large and Small Magellanic Clouds using the Mopra Telescope. Here we report on the basic physical properties of 125 GMCs in the LMC that have been surveyed to date. The observed clouds exhibit scaling relations that are similar to those determined for Galactic GMCs, although LMC clouds have narrower linewidths and lower CO luminosities than Galactic clouds of a similar size. The average mass surface density of the LMC clouds is 50 Msol/pc2, approximately half that of GMCs in the inner Milky Way. We compare the properties of GMCs with and without signs of massive star formation, finding that non-star-forming GMCs have lower peak CO brightness than star-forming GMCs. We compare the properties of GMCs with estimates for local interstellar conditions: specifically, we investigate the HI column density, radiation field, stellar mass surface density and the external pressure. Very few cloud properties demonstrate a clear dependence on the environment; the exceptions are significant positive correlations between i) the HI column density and the GMC velocity dispersion, ii) the stellar mass surface density and the average peak CO brightness, and iii) the stellar mass surface density and the CO surface brightness. The molecular mass surface density of GMCs without signs of massive star formation shows no dependence on the local radiation field, which is inconsistent with the photoionization-regulated star formation theory proposed by McKee (1989). We find some evidence that the mass surface density of the MAGMA clouds increases with the interstellar pressure, as proposed by Elmegreen (1989), but the detailed predictions of this model are not fulfilled once estimates for the local radiation field, metallicity and GMC envelope mass are taken into account.

We present here the first results from a high-resolution survey of the 12CO(J=1-0) emission across the northern part of the poorly-enriched Small Magellanic Cloud, made with the ATNF Mopra telescope. Three molecular complexes detected in the lower resolution NANTEN survey are mapped with a beam FWHM of ~42'', to sensitivities of approximately 210 mK per 0.9 km/s channel, resolving each complex into 4-7 small clouds of masses in the range of Mvir~10^3-4 Msol and with radii no larger than 16 pc. The northern SMC CO clouds follow similar empirical relationships to the southern SMC population, yet they appear relatively under-luminous for their size, suggesting that the star-forming environment in the SMC is not homogeneous. Our data also suggests that the CO cloud population has little or no extended CO envelope on scales<30 pc, further evidence that the weak CO component in the north SMC is being disassociated by penetrating UV radiation. The new high-resolution data provide evidence for a variable correlation of the CO integrated brightness with integrated HI and um emission; in particular CO is often, but not always, found coincident with peaks of um emission, verifying the need for matching-resolution um and HI data for a complete assessment of the SMC H2 mass.

Using our new 3-D relativistic electromagnetic particle (REMP) code parallelized with MPI, we investigated long-term particle acceleration associated with a relativistic electron-positron jet propagating in an unmagnetized ambient electron-positron plasma. We have also performed simulations with electron-ion jets. The simulations were performed using a much longer simulation system than our previous simulations in order to investigate the full nonlinear stage of the Weibel instability for electron-positron jets and its particle acceleration mechanism. Cold jet electrons are thermalized and ambient electrons are accelerated in the resulting shocks for both cases. Acceleration of ambient electrons leads to a maximum ambient electron density three times larger than the original value for pair plasmas. Behind the bow shock in the jet shock strong electromagnetic fields are generated. These fields may lead to time dependent afterglow emission. We calculated radiation from electrons propagating in a uniform parallel magnetic field to verify the technique. We also used the new technique to calculate emission from electrons based on simulations with a small system with two different cases for Lorentz factors (15 and 100). We obtained spectra which are consistent with those generated from electrons propagating in turbulent magnetic fields with red noise. This turbulent magnetic field is similar to the magnetic field generated at an early nonlinear stage of the Weibel instability.

New 12CO J=4-3 and 13CO J=3-2 observations of the N159 region in the Large Magellanic Cloud have been made. The 12CO J=4-3 distribution is separated into three clumps. These new measurements toward the three clumps are used in coupled calculations of molecular rotational excitation and line radiation transfer, along with other transitions of the 12CO as well as the isotope transitions of 13CO. The temperatures and densities are determined to be ~70-80K and ~3x10^3 cm-3 in N159W and N159E and ~30K and ~1.6x10^3 cm-3 in N159S. These results are compared with the star formation activity. The N159E clump is associated with embedded cluster(s) as observed at 24 micron and the derived high temperature is explained as due to the heating by these sources. The N159E clump is likely responsible for a dark lane in a large HII region by the dust extinction. The N159W clump is associated with embedded clusters mainly toward the eastern edge of the clump only. These clusters show offsets of 20"-40" from the 12CO J=4-3 peak and are probably responsible for heating indicated by the derived high temperature. The N159W clump exhibits no sign of star formation toward the 12CO J=4-3 peak position and its western region. We suggest that the N159W peak represents a pre-star-cluster core of ~105M_sol which deserves further detailed studies. Note that recent star formation took place between N159W and N159E as indicated by several star clusters and HII regions, while the natal molecular gas toward the stars have already been dissipated by the ionization and stellar winds of the OB stars. The N159S clump shows little sign of star formation as is consistent with the lower temperature and somewhat lower density. The N159S clump is also a candidate for future star formation.

We compare the CO J =(1-0) and HI emission in the Large Magellanic Cloud (LMC) in three dimensions, i.e. including a velocity axis in addition to the two spatial axes, with the aim of elucidating the physical connection between giant molecular clouds (GMCs) and their surrounding HI gas. The CO J =1-0 dataset is from the second NANTEN CO survey and the HI dataset is from the merged Australia Telescope Compact Array (ATCA) and Parkes Telescope surveys. The major findings of our analysis are: 1) GMCs are associated with an envelope of HI emission, 2) in GMCs [average CO intensity] is proportional to [average HI intensity]^[1.1+-0.1] and 3) the HI intensity tends to increase with the star formation activity within GMCs, from Type I to Type III. An analysis of the HI envelopes associated with GMCs shows that their average linewidth is 14 km s-1 and the mean density in the envelope is 10 cm-3. We argue that the HI envelopes are gravitationally bound by GMCs. These findings are consistent with a continual increase in the mass of GMCs via HI accretion at an accretion rate of 0.05 Msun/yr over a time scale of 10 Myr. The growth of GMCs is terminated via dissipative ionization and/or stellar-wind disruption in the final stage of GMC evolution.

We studied star formation activities in the molecular clouds in the Large Magellanic Cloud. We have utilized the second catalog of 272 molecular clouds obtained by NANTEN to compare the cloud distribution with signatures of massive star formation including stellar clusters, and optical and radio HII regions. We find that the molecular clouds are classified into three types according to the activities of massive star formation; Type I shows no signature of massive star formation, Type II is associated with relatively small HII region(s) and Type III with both HII region(s) and young stellar cluster(s). The radio continuum sources were used to confirm that Type I GMCs do not host optically hidden HII regions. These signatures of massive star formation show a good spatial correlation with the molecular clouds in a sense they are located within ~100 pc of the molecular clouds. Among possible ideas to explain the GMC Types, we favor that the Types indicate an evolutionary sequence; i.e., the youngest phase is Type I, followed by Type II and the last phase is Type III, where the most active star formation takes place leading to cloud dispersal. The number of the three types of GMCs should be proportional to the time scale of each evolutionary stage if a steady state of massive star and cluster formation is a good approximation. By adopting the time scale of the youngest stellar clusters, 10 Myrs, we roughly estimate the timescales of Types I, II and III to be 6 Myrs, 13 Myrs and 7 Myrs, respectively, corresponding to a lifetime of 20-30 Myrs for the GMCs with a mass above the completeness limit, 5 x 10^4 Msun.

Plasma instabilities (e.g., Buneman, Weibel and other two-stream instabilities) excited in collisionless shocks are responsible for particle (electron, positron, and ion) acceleration. Using a new 3-D relativistic particle-in-cell code, we have investigated the particle acceleration and shock structure associated with an unmagnetized relativistic electron-positron jet propagating into an unmagnetized electron-positron plasma. The simulation has been performed using a long simulation system in order to study the nonlinear stages of the Weibel instability, the particle acceleration mechanism, and the shock structure. Cold jet electrons are thermalized and slowed while the ambient electrons are swept up to create a partially developed hydrodynamic (HD) like shock structure. In the leading shock, electron density increases by a factor of 3.5 in the simulation frame. Strong electromagnetic fields are generated in the trailing shock and provide an emission site. We discuss the possible implication of our simulation results within the AGN and GRB context.

We have investigated the development of current-driven (CD) kink instability through three-dimensional relativistic MHD simulations. A static force-free equilibrium helical magnetic configuration is considered in order to study the influence of the initial configuration on the linear and nonlinear evolution of the instability. We found that the initial configuration is strongly distorted but not disrupted by the kink instability. The instability develops as predicted by linear theory. In the non-linear regime the kink amplitude continues to increase up to the terminal simulation time, albeit at different rates, for all but one simulation. The growth rate and nonlinear evolution of the CD kink instability depends moderately on the density profile and strongly on the magnetic pitch profile. The growth rate of the kink mode is reduced in the linear regime by an increase in the magnetic pitch with radius and the non-linear regime is reached at a later time than for constant helical pitch. On the other hand, the growth rate of the kink mode is increased in the linear regime by a decrease in the magnetic pitch with radius and reaches the non-linear regime sooner than the case with constant magnetic pitch. Kink amplitude growth in the non-linear regime for decreasing magnetic pitch leads to a slender helically twisted column wrapped by magnetic field. On the other hand, kink amplitude growth in the non-linear regime nearly ceases for increasing magnetic pitch.

We analyze the conditions for detection of CO(1-0) emission in the Large Magellanic Cloud (LMC), using the recently completed second NANTEN CO survey. In particular, we investigate correlations between CO integrated intensity and HI integrated intensity, peak brightness temperature, and line width at a resolution of 2.6' (~40 pc). We find that significant HI column density and peak brightness temperature are necessary but not sufficient conditions for CO detection, with many regions of strong HI emission not associated with molecular clouds. The large scatter in CO intensities for a given HI intensity persists even when averaging on scales of >200 pc, indicating that the scatter is not solely due to local conversion of HI into H_2 near GMCs. We focus on two possibilities to account for this scatter: either there exist spatial variations in the I(CO) to N(H_2) conversion factor, or a significant fraction of the atomic gas is not involved in molecular cloud formation. A weak tendency for CO emission to be suppressed for large HI linewidths supports the second hypothesis, insofar as large linewidths may be indicative of warm HI, and calls into question the likelihood of forming molecular clouds from colliding HI flows. We also find that the ratio of molecular to atomic gas shows no significant correlation (or anti-correlation) with the stellar surface density, though a correlation with midplane hydrostatic pressure P_h is found when the data are binned in P_h. The latter correlation largely reflects the increasing likelihood of CO detection at high HI column density.

Using our new 3-D relativistic electromagnetic particle (REMP) code parallelized with MPI, we have investigated long-term particle acceleration associated with an relativistic electron-positron jet propagating in an unmagnetized ambient electron-positron plasma. The simulations have been performed using a much longer simulation system than our previous simulations in order to investigate the full nonlinear stage of the Weibel instability and its particle acceleration mechanism. Cold jet electrons are thermalized and ambient electrons are accelerated in the resulting shocks. The acceleration of ambient electrons leads to a maximum ambient electron density three times larger than the original value. Behind the bow shock in the jet shock strong electromagnetic fields are generated. These fields may lead to the afterglow emission. We have calculated the time evolution of the spectrum from two electrons propagating in a uniform parallel magnetic field to verify the technique.

We solve the Riemann problem for the deceleration of an arbitrarily magnetized relativistic flow injected into a static unmagnetized medium in one dimension. We find that for the same initial Lorentz factor, the reverse shock becomes progressively weaker with increasing magnetization \sigma (the Poynting-to kinetic energy flux ratio), and the shock becomes a rarefaction wave when \sigma exceeds a critical value, \sigma_c, defined by the balance between the magnetic pressure in the flow and the thermal pressure in the forward shock. In the rarefaction wave regime, we find that the rarefied region is accelerated to a Lorentz factor that is significantly larger than the initial value. This acceleration mechanism is due to the strong magnetic pressure in the flow. We discuss the implications of these results for models of gamma-ray bursts and active galactic nuclei.

Nonthermal radiation observed from astrophysical systems containing relativistic jets and shocks, e.g., gamma-ray bursts (GRBs), active galactic nuclei (AGNs), and microquasars commonly exhibit power-law emission spectra. Recent PIC simulations of relativistic electron-ion (or electron-positron) jets injected into a stationary medium show that particle acceleration occurs within the downstream jet. In collisionless, relativistic shocks, particle (electron, positron, and ion) acceleration is due to plasma waves and their associated instabilities (e.g., the Weibel (filamentation) instability) created in the shock region. The simulations show that the Weibel instability is responsible for generating and amplifying highly non-uniform, small-scale magnetic fields. These fields contribute to the electron's transverse deflection behind the jet head. The resulting "jitter" radiation from deflected electrons has different properties compared to synchrotron radiation, which assumes a uniform magnetic field. Jitter radiation may be important for understanding the complex time evolution and/or spectra in gamma-ray bursts, relativistic jets in general, and supernova remnants.

Nonthermal radiation observed from astrophysical systems containing relativistic jets and shocks, e.g., gamma-ray bursts (GRBs), active galactic nuclei (AGNs), and Galactic microquasar systems usually have power-law emission spectra. Recent PIC simulations of relativistic electron-ion (electron-positron) jets injected into a stationary medium show that particle acceleration occurs within the downstream jet. In the presence of relativistic jets, instabilities such as the Buneman instability, other two-streaming instability, and the Weibel (filamentation) instability create collisionless shocks, which are responsible for particle (electron, positron, and ion) acceleration. The simulation results show that the Weibel instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields. These magnetic fields contribute to the electron's transverse deflection behind the jet head. The ``jitter'' radiation from deflected electrons in small-scale magnetic fields has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation, a case of diffusive synchrotron radiation, may be important to understand the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.

The second survey of the molecular clouds in 12CO (J = 1-0) was carried out in the Large Magellanic Cloud by NANTEN. The sensitivity of this survey is twice as high as that of the previous NANTEN survey, leading to a detection of molecular clouds with M_CO > 2 x 10^4 M_sun. We identified 272 molecular clouds, 230 of which are detected at three or more observed positions. We derived the physical properties, such as size, line width, virial mass, of the 164 GMCs which have an extent more than the beam size of NANTEN in both the major and minor axes. The CO luminosity and virial mass of the clouds show a good correlation of M_VIR propto L_CO^{1.1 +- 0.1} with a Spearman rank correlation of 0.8 suggesting that the clouds are in nearly virial equilibrium. Assuming the clouds are in virial equilibrium, we derived an X_CO-factor to be ~ 7 x 10^20 cm^-2 (K km s^-1)^-1. The mass spectrum of the clouds is fitted well by a power law of N_cloud(>M_CO) proportional to M_CO^{-0.75 +- 0.06} above the completeness limit of 5 x 10^4 M_sun. The slope of the mass spectrum becomes steeper if we fit only the massive clouds; e.g., N_cloud (>M_CO) is proportional to M_CO^{-1.2 +- 0.2} for M_CO > 3 x 10^5 M_sun.

Nonthermal radiation observed from astrophysical systems containing relativistic jets and shocks e.g. gamma-ray bursts (GRBs) active galactic nuclei (AGNs) and microquasars commonly exhibit power-law emission spectra. Recent PIC simulations of relativistic electron-ion (or electron-positron) jets injected into a stationary medium show that particle acceleration occurs within the downstream jet. In collisionless relativistic shocks particle (electron, positron and ion) acceleration is due to plasma waves and their associated instabilities (e.g. the Weibel (filamentation) instability) created in the shock region. The simulations show that the Weibel instability is responsible for generating and amplifying highly non-uniform small-scale magnetic fields. These fields contribute to the electron's transverse deflection behind the jet head. The resulting ``jitter'' radiation from deflected electrons has different properties compared to synchrotron radiation which assumes a uniform magnetic field. Jitter radiation may be important for understanding the complex time evolution and/or spectra in gamma-ray bursts, relativistic jets in general and supernova remnants.

Star formation at earlier cosmological times takes place in an interstellar medium with low metallicity. The Large Magellanic Cloud (LMC) is ideally suited to study star formation in such an environment. The physical and chemical state of the ISM in a star forming environment can be constrained by observations of submm and FIR spectral lines of the main carbon carrying species, CO, CI and CII, which originate in the surface layers of molecular clouds illuminated by the UV radiation of the newly formed, young stars. We present high-angular resolution sub-millimeter observations in the N159W region in the LMC obtained with the NANTEN2 telescope of the 12CO J = 4-3, J = 7-6, and 13CO J = 4-3 rotational and [CI] 3P1-3P0 and 3P2-3P1 fine-structure transitions. The 13CO J =4-3 and [CI] 3P2-3P1 transitions are detected for the first time in the LMC. We derive the physical and chemical properties of the low-metallicity molecular gas using an escape probability code and a self-consistent solution of the chemistry and thermal balance of the gas in the framework of a clumpy cloud PDR model. The separate excitation analysis of the submm CO lines and the carbon fine structure lines shows that the emitting gas in the N159W region has temperatures of about 80 K and densities of about 10^4 cm^-3. The estimated C to CO abundance ratio close to unity is substantially higher than in dense massive star-forming regions in the Milky Way. The analysis of all observed lines together, including the [CII] line intensity reported in the literature, in the context of a clumpy cloud PDR model constrains the UV intensity to about \chi ~220 and an average density of the clump ensemble of about 10^5 cm^-3, thus confirming the presence of high density material in the LMC N159W region.

Nonthermal radiation observed from astrophysical systems containing relativistic jets and shocks, e.g., active galactic nuclei (AGNs), gamma-ray bursts (GRBs), and Galactic microquasar systems usually have power-law emission spectra. Recent PIC simulations using injected relativistic electron-ion (electro-positron) jets show that acceleration occurs within the downstream jet. Shock acceleration is a ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., the Buneman instability, other two-streaming instability, and the Weibel instability) created in the shocks are responsible for particle (electron, positron, and ion) acceleration. The simulation results show that the Weibel instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields. These magnetic fields contribute to the electron's transverse deflection behind the jet head. The ``jitter'' radiation from deflected electrons has different properties than synchrotron radiation which assumes a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.

We have carried out sub-mm 12CO(J=3-2) observations of 6 giant molecular clouds (GMCs) in the Large Magellanic Cloud (LMC) with the ASTE 10m sub-mm telescope at a spatial resolution of 5 pc and very high sensitivity. We have identified 32 molecular clumps in the GMCs and revealed significant details of the warm and dense molecular gas with n(H2) $\sim$ 10$^{3-5}$ cm$^{-3}$ and Tkin $\sim$ 60 K. These data are combined with 12CO(J=1-0) and 13CO(J=1-0) results and compared with LVG calculations. We found that the ratio of 12CO(J=3-2) to 12CO(J=1-0) emission is sensitive to and is well correlated with the local Halpha flux. We interpret that differences of clump propeties represent an evolutionary sequence of GMCs in terms of density increase leading to star formation.Type I and II GMCs (starless GMCs and GMCs with HII regions only, respectively) are at the young phase of star formation where density does not yet become high enough to show active star formation and Type III GMCs (GMCs with HII regions and young star clusters) represents the later phase where the average density is increased and the GMCs are forming massive stars. The high kinetic temperature correlated with \Halpha flux suggests that FUV heating is dominant in the molecular gas of the LMC.

We performed relativistic magnetohydrodynamic simulations of the hydrodynamic boosting mechanism for relativistic jets explored by Aloy & Rezzolla (2006) using the RAISHIN code. Simulation results show that the presence of a magnetic field changes the properties of the shock interface between the tenuous, overpressured jet ($V^z_j$) flowing tangentially to a dense external medium. Magnetic fields can lead to more efficient acceleration of the jet, in comparison to the pure-hydrodynamic case. A ``poloidal'' magnetic field ($B^z$), tangent to the interface and parallel to the jet flow, produces both a stronger outward moving shock and a stronger inward moving rarefaction wave. This leads to a large velocity component normal to the interface in addition to acceleration tangent to the interface, and the jet is thus accelerated to larger Lorentz factors than those obtained in the pure-hydrodynamic case. Likewise, a strong ``toroidal'' magnetic field ($B^y$), tangent to the interface but perpendicular to the jet flow, also leads to stronger acceleration tangent to the shock interface relative to the pure-hydrodynamic case. Overall, the acceleration efficiency in the ``poloidal'' case is less than that of the ``toroidal'' case but both geometries still result in higher Lorentz factors than the pure-hydrodynamic case. Thus, the presence and relative orientation of a magnetic field in relativistic jets can significant modify the hydrodynamic boost mechanism studied by Aloy & Rezzolla (2006).

A new general relativistic magnetohydrodynamics (GRMHD) code ``RAISHIN'' used to simulate jet generation by rotating and non-rotating black holes with a geometrically thin Keplarian accretion disk finds that the jet develops a spine-sheath structure in the rotating black hole case. Spine-sheath structure and strong magnetic fields significantly modify the Kelvin-Helmholtz (KH) velocity shear driven instability. The RAISHIN code has been used in its relativistic magnetohydrodynamic (RMHD) configuration to study the effects of strong magnetic fields and weakly relativistic sheath motion, c/2, on the KH instability associated with a relativistic, Lorentz factor equal 2.5, jet spine-sheath interaction. In the simulations sound speeds up to c/1.7 and Alfven wave speeds up to 0.56 c are considered. Numerical simulation results are compared to theoretical predictions from a new normal mode analysis of the RMHD equations. Increased stability of a weakly magnetized system resulting from c/2 sheath speeds and stabilization of a strongly magnetized system resulting from c/2 sheath speeds is found.

Numerical simulations of weakly magnetized and strongly magnetized relativistic jets embedded in a weakly magnetized and strongly magnetized stationary or weakly relativistic (v = c/2) sheath have been performed. A magnetic field parallel to the flow is used in these simulations performed by the new GRMHD numerical code RAISHIN used in its RMHD configuration. In the numerical simulations the Lorentz factor $\gamma = 2.5$ jet is precessed to break the initial equilibrium configuration. In the simulations sound speeds are $\lesssim c/\sqrt 3$ in the weakly magnetized simulations and $\lesssim 0.3c$ in the strongly magnetized simulations. The Alfven wave speed is $\lesssim 0.07c$ in the weakly magnetized simulations and $\lesssim 0.56c$ in the strongly magnetized simulations. The results of the numerical simulations are compared to theoretical predictions from a normal mode analysis of the linearized relativistic magnetohydrodynamic (RMHD) equations capable of describing a uniform axially magnetized cylindrical relativistic jet embedded in a uniform axially magnetized relativistically moving sheath. The theoretical dispersion relation allows investigation of effects associated with maximum possible sound speeds, Alfven wave speeds near light speed and relativistic sheath speeds. The prediction of increased stability of the weakly magnetized system resulting from c/2 sheath speeds and the stabilization of the strongly magnetized system resulting from c/2 sheath speeds is verified by the numerical simulation results.

We investigate jet formation in black-hole systems using 3-D General Relativistic Particle-In-Cell (GRPIC) and 3-D GRMHD simulations. GRPIC simulations, which allow charge separations in a collisionless plasma, do not need to invoke the frozen condition as in GRMHD simulations. 3-D GRPIC simulations show that jets are launched from Kerr black holes as in 3-D GRMHD simulations, but jet formation in the two cases may not be identical. Comparative study of black hole systems with GRPIC and GRMHD simulations with the inclusion of radiate transfer will further clarify the mechanisms that drive the evolution of disk-jet systems.

We have performed numerical simulations of weakly and strongly magnetized relativistic jets embedded in a weakly and strongly magnetized stationary or mildly relativistic (0.5c) sheath using the RAISHIN code. In the numerical simulations a jet with Lorentz factor \gamma=2.5 is precessed to break the initial equilibrium configuration. Results of the numerical simulations are compared to theoretical predictions from a normal mode analysis of the linearized RMHD equations describing a uniform axially magnetized cylindrical relativistic jet embedded in a uniform axially magnetized moving sheath. The prediction of increased stability of a weakly-magnetized system with mildly relativistic sheath flow to Kelvin-Helmholtz instabilities and the stabilization of a strongly-magnetized system with mildly relativistic sheath flow is confirmed by the numerical simulations.

We have performed several simulations of black hole systems (non-rotating, black hole spin parameter a=0.0 and rapidly rotating, a=0.95) with a geometrically thin Keplerian disk using the newly developed RAISHIN code. The simulation results show the formation of jets driven by the Lorentz force and the gas pressure gradient. The jets have mildly relativistic speed (\sim 0.4c). The matter is continuously supplied from the accretion disk and the jet propagates outward until each applicable terminal simulation time (non-rotating: t/\tau_S = 275 and rotating: t/tau_S = 200, \tau_S \equiv r_S/c). It appears that a rotating black hole creates an additional, faster, and more collimated matter-dominated inner outflow (\sim 0.5c) formed and accelerated by the twisted magnetic field resulting from frame-dragging in the black hole ergosphere. This is the first known simulation confirming the formation of an inner magnetically-driven, matter-dominated jet by the frame-dragging effect from a black hole co-rotating with a thin Keplerian disk threaded by a vertical magnetic field. This result indicates that jet kinematic structure depends on black hole rotation and on the initial magnetic field configuration and strength.

We have developed a new three-dimensional general relativistic magnetohydrodynamic (GRMHD) code, RAISHIN, using a conservative, high resolution shock-capturing scheme. The numerical fluxes are calculated using the Harten, Lax, & van Leer (HLL) approximate Riemann solver scheme. The flux-interpolated, constrained transport scheme is used to maintain a divergence-free magnetic field. In order to examine the numerical accuracy and the numerical efficiency, the code uses four different reconstruction methods: piecewise linear methods with Minmod and MC slope-limiter function, convex essentially non-oscillatory (CENO) method, and piecewise parabolic method (PPM) using multistep TVD Runge-Kutta time advance methods with second and third-order time accuracy. We describe code performance on an extensive set of test problems in both special and general relativity. Our new GRMHD code has proven to be accurate in second order and has successfully passed with all tests performed, including highly relativistic and magnetized cases in both special and general relativity.

We have performed a fully three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulation of jet formation from a thin accretion disk around a Schwarzschild black hole with a free-falling corona. The initial simulation results show that a bipolar jet (velocity nearly 0.3c) is created as shown by previous two-dimensional axisymmetric simulations with mirror symmetry at the equator. The 3-D simulation ran over one hundred light-crossing time units which is considerably longer than the previous simulations. We show that the jet is initially formed as predicted due in part to magnetic pressure from the twisting the initially uniform magnetic field and from gas pressure associated with shock formation in the region around r = 3r_S. At later times, the accretion disk becomes thick and the jet fades resulting in a wind that is ejected from the surface of the thickened (torus-like) disk. It should be noted that no streaming matter from a donor is included at the outer boundary in the simulation (an isolated black hole not binary black hole). The wind flows outwards with a wider angle than the initial jet. The widening of the jet is consistent with the outward moving torsional Alfven waves (TAWs). This evolution of disk-jet coupling suggests that the jet fades with a thickened accretion disk due to the lack of streaming material from an accompanying star.

Weibel instability created in collisionless shocks is responsible for particle (electron, positron, and ion) acceleration. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic electron-ion jet fronts propagating into an ambient plasma without initial magnetic fields with a longer simulation system in order to investigate nonlinear stage of the Weibel instability and its acceleration mechanism. The current channels generated by the Weibel instability induce the radial electric fields. The z component of the Poynting vector (E x B) become positive in the large region along the jet propagation direction. This leads to the acceleration of jet electrons along the jet. In particular the E x B drift with the large scale current channel generated by the ion Weibel instability accelerate electrons effectively in both parallel and perpendicular directions.

We have performed 2.5-dimensional general relativistic magnetohydrodynamic (MHD) simulations of collapsars including a rotating black hole. This paper is an extension of our previous paper (Mizuno et al. 2004). The current calculation focuses on the effect of black hole rotation using general relativistic MHD with simplified microphysics. Initially, we assume that the core collapse is failed in this star. A few solar mass rotating black hole is inserted by hand into the calculation. We consider two cases, the co-rotating case and counter-rotating case with respect to the black hole rotation. Although the counter-rotating case may be unrealistic for collapsar, we perform as a maximally dragging case of magnetic field. The simulation results show the formation of a disk-like structure and the generation of a jet-like outflow near the central black hole. The jet-like outflow propagates outwardly with the twisted magnetic field and becomes collimated. We have found that the jets are generated and accelerated mainly by the magnetic field. The total jet velocity in the rotating black hole case is comparable to that of the non-rotating black hole case (Mizuno et al. 2004), 0.3c. When the rotation of the black hole is faster, the magnetic field is twisted strongly owing to the frame-dragging effect. The magnetic energy stored by the twisting magnetic field is converted to kinetic energy of the jet directly rather than propagating as an Alfven wave. Thus, as the rotation of the black hole becomes faster, the poloidal velocity of the jet becomes faster. In the rapidly rotating black hole case the jet-like outflow can be produced by the frame dragging effect only through the twisting of magnetic field even if there is no stellar rotation.

We have performed 2.5-dimensional general relativistic magnetohydrodynamic (MHD) simulations of the gravitational collapse of a magnetized rotating massive star as a model of gamma ray bursts (GRBs). The current calculation focuses on general relativistic MHD with simplified microphysics (we ignore neutrino cooling, physical equation of state, and photodisintegration). Initially, we assume that the core collapse has failed in this star. A few $M_{\odot}$ black hole is inserted by hand into the calculation. The simulations presented in the paper follow the accretion of gas into a black hole that is assumed to have formed before the calculation begins.The simulation results show the formation of a disk-like structure and the generation of a jetlike outflow inside the shock wave launched at the core bounce. We have found that the jet is accelerated by the magnetic pressure and the centrifugal force and is collimated by the pinching force of the toroidal magnetic field amplified by the rotation and the effect of geometry of the poloidal magnetic field. The maximum velocity of the jet is mildly relativistic (0.3c). The velocity of the jet becomes larger as the initial rotational velocity of stellar matter gets faster. On the other hand, the dependence on the initial magnetic field strength is a bit more complicated: the velocity of the jet increases with the initial field strength in the weak field regime, then is saturated at some intermediate field strength, and decreases beyond the critical field strength. These results are related to the stored magnetic energy determined by the balance between the propagation time of the Alfven wave and the rotation time of the disk (or twisting time).

We have performed 2.5-dimensional general relativistic magnetohydrodynamic (MHD) simulations of the gravitational collapse of a magnetized rotating massive star as a model of gamma ray bursts (GRBs). This simulation showed the formation of a disk-like structure and the generation of a jet-like outflow inside the shock wave launched at the core bounce. We have found the jet is accelerated by the magnetic pressure and the centrifugal force and is collimated by the pinching force of the toroidal magnetic field amplified by the rotation and the effect of geometry of the poloidal magnetic field. The maximum velocity of the jet is mildly relativistic ($\sim$ 0.3 c).