Jet-torus connection in radio galaxies: Relativistic hydrodynamics and synthetic emission

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.Comment: 14 pages, 17 figures, submitted to A&

Flares in the Galactic Centre – I. Orbiting flux tubes in magnetically arrested black hole accretion discs

Recent observations of Sgr A* by the GRAVITY instrument have astrometrically tracked infrared (IR) flares at distances of ∼10 gravitational radii (rg). In this paper, we study a model for the flares based on 3D general relativistic magnetohydrodynamic (GRMHD) simulations of magnetically arrested accretion discs (MADs) that exhibit violent episodes of flux escape from the black hole magnetosphere. These events are attractive for flare modelling for several reasons: (i) the magnetically dominant regions can resist being disrupted via magnetorotational turbulence and shear; (ii) the orientation of the magnetic field is predominantly vertical as suggested by the GRAVITY data; and (iii) the 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 ∼5 to 40rg⁠. The magnetic energy contained within the flux bundles ranges up to ∼1040erg⁠, 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

Flares in the Galactic Centre - I:Orbiting flux tubes in magnetically arrested black hole accretion discs

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.Comment: accepted for publication by MNRAS, 18-Jan-202

Star Forming Dense Cloud Cores in the TeV {\gamma}-ray SNR RX J1713.7-3946

RX J1713.7-3946 is one of the TeV {\gamma}-ray supernova remnants (SNRs) emitting synchrotron X rays. The SNR is associated with molecular gas located at ~1 kpc. We made new molecular observations toward the dense cloud cores, peaks A, C and D, in the SNR in the 12CO(J=2-1) and 13CO(J=2-1) transitions at angular resolution of 90". The most intense core in 13CO, peak C, was also mapped in the 12CO(J=4-3) transition at angular resolution of 38". Peak C shows strong signs of active star formation including bipolar outflow and a far-infrared protostellar source and has a steep gradient with a r^{-2.2$\pm$0.4} variation in the average density within radius r. Peak C and the other dense cloud cores are rim-brightened in synchrotron X rays, suggesting that the dense cloud cores are embedded within or on the outer boundary of the SNR shell. This confirms the earlier suggestion that the X rays are physically associated with the molecular gas (Fukui et al. 2003). We present a scenario where the densest molecular core, peak C, survived against the blast wave and is now embedded within the SNR. Numerical simulations of the shock-cloud interaction indicate that a dense clump can indeed survive shock erosion, since shock propagation speed is stalled in the dense clump. Additionally, the shock-cloud interaction induces turbulence and magnetic field amplification around the dense clump that may facilitate particle acceleration in the lower-density inter-clump space leading to the enhanced synchrotron X rays around dense cores.Comment: 22 pages, 7 figures, to accepted in The Astrophysical Journal. A full color version with higher resolution figures is available at http://www.a.phys.nagoya-u.ac.jp/~sano/ApJ10/ms_sano.pd

State of Annual Paid Leave–Doctors’ Working Conditions

The article's abstract is not available. &nbsp

Oligarchic planetesimal accretion and giant planet formation II

The equation of state calculated by Saumon and collaborators has been adopted in most core-accretion simulations of giant-planet formation performed to date. Since some minor errors have been found in their original paper, we present revised simulations of giant-planet formation that considers a corrected equation of state. We employ the same code as Fortier and collaborators in repeating our previous simulations of the formation of Jupiter. Although the general conclusions of Fortier and collaborators remain valid, we obtain significantly lower core masses and shorter formation times in all cases considered. The minor errors in the previously published equation of state have been shown to affect directly the adiabatic gradient and the specific heat, causing an overestimation of both the core masses and formation times.Comment: 4 pages, 2 figures, Accepted for publication in Astronomy and Astrophysic

Dynamical structure factors of S=1/2S=1/2 two-leg spin ladder systems

We investigate dynamical properties of $S=1/2$ two-leg spin ladder systems. In a strong coupling region, an isolated mode appears in the lowest excited states, while in a weak coupling region, an isolated mode is reduced and the lowest excited states become a lower bound of the excitation continuum. We find in the system with equal intrachain and interchain couplings that due to a cyclic four-spin interaction, the distribution of the weights for the dynamical structure factor and characteristics of the lowest excited states are strongly influenced. The dynamical properties of two systems proposed for ${\rm SrCu_2O_3}$ are also discussed.Comment: 5 pages, 6 figure

Enhanced reaction kinetics in biological cells

The cell cytoskeleton is a striking example of "active" medium driven out-of-equilibrium by ATP hydrolysis. Such activity has been shown recently to have a spectacular impact on the mechanical and rheological properties of the cellular medium, as well as on its transport properties : a generic tracer particle freely diffuses as in a standard equilibrium medium, but also intermittently binds with random interaction times to motor proteins, which perform active ballistic excursions along cytoskeletal filaments. Here, we propose for the first time an analytical model of transport limited reactions in active media, and show quantitatively how active transport can enhance reactivity for large enough tracers like vesicles. We derive analytically the average interaction time with motor proteins which optimizes the reaction rate, and reveal remarkable universal features of the optimal configuration. We discuss why active transport may be beneficial in various biological examples: cell cytoskeleton, membranes and lamellipodia, and tubular structures like axons.Comment: 10 pages, 2 figure

The role of the initial surface density profiles of the disc on giant planet formation: comparing with observations

In order to explain the main characteristics of the observed population of extrasolar planets and the giant planets in the Solar System, we need to get a clear understanding of which are the initial conditions that allowed their formation. To this end we develop a semi-analytical model for computing planetary systems formation based on the core instability model for the gas accretion of the embryos and the oligarchic growth regime for the accretion of the solid cores. With this model we explore not only different initial discs profiles motivated by similarity solutions for viscous accretion discs, but we also consider different initial conditions to generate a variety of planetary systems assuming a large range of discs masses and sizes according to the last results in protoplanetary discs observations. We form a large population of planetary systems in order to explore the effects in the formation of assuming different discs and also the effects of type I and II regimes of planetary migration, which were found to play fundamental role in reproducing the distribution of observed exoplanets. Our results show that the observed population of exoplanets and the giant planets in the Solar System are well represented when considering a surface density profile with a power law in the inner part characterized by an exponent of -1, which represents a softer profile when compared with the case most similar to the MMSN model case.Comment: 14 pages, 12 figures, MNRAS, 412, 211

The role of the initial surface density profiles of the disc on giant planet formation: comparing with observations

In order to explain the main characteristics of the observed population of extrasolar planets and the giant planets in the Solar System, we need to get a clear understanding of which are the initial conditions that allowed their formation. To this end we develop a semi-analytical model for computing planetary systems formation based on the core instability model for the gas accretion of the embryos and the oligarchic growth regime for the accretion of the solid cores. With this model we explore not only different initial discs profiles motivated by similarity solutions for viscous accretion discs, but we also consider different initial conditions to generate a variety of planetary systems assuming a large range of discs masses and sizes according to the last results in protoplanetary discs observations. We form a large population of planetary systems in order to explore the effects in the formation of assuming different discs and also the effects of type I and II regimes of planetary migration, which were found to play fundamental role in reproducing the distribution of observed exoplanets. Our results show that the observed population of exoplanets and the giant planets in the Solar System are well represented when considering a surface density profile with a power law in the inner part characterized by an exponent of -1, which represents a softer profile when compared with the case most similar to the MMSN model case.Comment: 14 pages, 12 figures, MNRAS, 412, 211