Pulsar wind nebulae (PWN) are fascinating systems and archetypal sources for high-energy astrophysics in general. Due to their vicinity, brightness, to the fact that they shine at multi-wavelengths, and especially to their long-living emission at gamma rays, modelling their properties is particularly important for the correct interpretation of the visible Galaxy. A complication in this respect is the variety of properties and morphologies they show at different ages. Here, we discuss the differences among the evolutionary phases of PWN, how they have been modeled in the past and what progresses have been recently made. We approach the discussion from a phenomenological, theoretical (especially numerical) and observational point of view, with particular attention to the most recent results and open questions about the physics of such intriguing sources.

Therecently discovered massive binary system Apep is the most powerful synchrotron emitter among the known Galactic colliding-wind binaries. This makes this particular system of great interest to investigate stellar winds and the non-thermal processes associated with their shocks. This source was detected at various radio bands, and in addition the wind-collision region was resolved by means of very-long baseline interferometric observations. We use a non-thermal emission model for colliding-wind binaries to derive physical properties of this system. The observed morphology in the resolved maps allows us to estimate the system projection angle on the sky to be $\psi \approx 85^\circ$ . The observed radio flux densities also allow us to characterise both the intrinsic synchrotron spectrum of the source and its modifications due to free–free absorption in the stellar winds at low frequencies; from this, we derive mass–loss rates of the stars of $\dot{M}_\mathrm{WN} \approx 4\times10^{-5}\;\mathrm{M}_\odot\,\mathrm{yr}^{-1}$ and $\dot{M}_\mathrm{WC} \approx 2.9\times10^{-5}\;\mathrm{M}_\odot\,\mathrm{yr}^{-1}$ . Finally, the broadband spectral energy distribution is calculated for different combinations of the remaining free parameters, namely the intensity of the magnetic field and the injected power in non-thermal particles. We show that the degeneracy of these two parameters can be solved with observations in the high-energy domain, most likely in the hard X-rays but also possibly in $\gamma$ -rays under favourable conditions.

We present a broadband radio study of the transient jets ejected from the black hole candidate X-ray binary MAXI J1535–571, which underwent a prolonged outburst beginning on 2017 September 2. We monitored MAXI J1535–571 with the Murchison Widefield Array (MWA) at frequencies from 119 to 186 MHz over six epochs from 2017 September 20 to 2017 October 14. The source was quasi-simultaneously observed over the frequency range 0.84–19 GHz by UTMOST (the Upgraded Molonglo Observatory Synthesis Telescope) the Australian Square Kilometre Array Pathfinder (ASKAP), the Australia Telescope Compact Array (ATCA), and the Australian Long Baseline Array (LBA). Using the LBA observations from 2017 September 23, we measured the source size to be $34\pm1$ mas. During the brightest radio flare on 2017 September 21, the source was detected down to 119 MHz by the MWA, and the radio spectrum indicates a turnover between 250 and 500 MHz, which is most likely due to synchrotron self-absorption (SSA). By fitting the radio spectrum with a SSA model and using the LBA size measurement, we determined various physical parameters of the jet knot (identified in ATCA data), including the jet opening angle ( $\phi_{\rm op} = 4.5\pm1.2^{\circ}$ ) and the magnetic field strength ( $B_{\rm s} = 104^{+80}_{-78}$ mG). Our fitted magnetic field strength agrees reasonably well with that inferred from the standard equipartition approach, suggesting the jet knot to be close to equipartition. Our study highlights the capabilities of the Australian suite of radio telescopes to jointly probe radio jets in black hole X-ray binaries via simultaneous observations over a broad frequency range, and with differing angular resolutions. This suite allows us to determine the physical properties of X-ray binary jets. Finally, our study emphasises the potential contributions that can be made by the low-frequency part of the Square Kilometre Array (SKA-Low) in the study of black hole X-ray binaries.

To study the dynamics of relativistic flows in astrophysical objects such as radio jets, we have developed a new special relativistic hydrodynamic (RHD) code based on the weighted essentially non-oscillatory (WENO) scheme, a high-order finite difference scheme. The code includes different WENO versions, and high-order time integration methods such as the 4th-order accurate Runge-Kutta (RK4) and strong stability preserving RK (SSPRK), as well as the equations of state (EOSs) that closely approximate the EOS of the single-component perfect gas in relativistic regime. Additionally, it is optimized for the reproduction of complex structures in multi-dimensional flows, and implements a modification of eigenvalues for the acoustic modes to effectively control carbuncle instability. As the first application of the code, we have simulated ultra-relativistic jets of FR II radio galaxies, and studied the nonlinear flow structures, such as shocks, velocity shear, and turbulence, through large-scale.

Following the discovery of quantum phenomena at laboratory scale (Couder & Fort 2006), de Broglie pilot wave theory (De Broglie 1962) has been revived under a hydrodynamic guise (Bush 2015). Theoretically, it boils down to solving the transport equations for the energy and linear momentum densities of a postulated fundamental fluid in terms of classical wave equations, which inherently are Lorentz-invariant and scale-invariant. Instead of the conventional harmonic solutions, for astronomical and gravitational problems the novel solutions for the homogeneous wave equation in spherical coordinates are more suitable (Munera et al. 1995, Munera & Guzman 1997, and Munera 2000). Two groups of solutions are particularly relevant: (a) The inherently-quantized helicoidal solutions that may be applicable to describe spiral galaxies, and (b) The non-harmonic solutions with time (t) and distance (r) entangled in the single variable q = Ct/r (C is the two-way local electromagnetic speed). When these functions are plotted against 1/q they manifestly depict quantum effects in the near field, and Newtonian-like gravity in the far-field. The near-field predicts quantized effects similar to ring structures and to Titius-Bode structures, both in our own solar system and in exoplanets, the correlation between predicted and observed structures being typically larger than 99 per cent. In the far-field, some non-harmonic functions have a rate of decrement with distance slower than inverse-square thus explaining the flat rotation rate of galaxies. Additional implications for Trojan orbits, and quantized effects in photon deflection were also noted.

The emission of steady compact jets observed in the hard spectral state of X-ray binaries is likely to be powered by internal shocks caused by fluctuations of the outflow velocity. The dynamics of the internal shocks and the resulting spectral energy distribution (SED) of the jet is very sensitive to the shape of the Power Spectral Density (PSD) of the fluctuations of the jet Lorentz factor. It turns out that Lorentz factor fluctuations injected at the base of the jet with a flicker noise power spectrum (i.e. P(f) ∝1/f) naturally produce the canonical flat SED observed from radio to IR band in X-ray binary systems in the hard state. This model also predicts a strong, wavelength dependent, variability that resembles the observed one. In particular, strong sub-second variability is predicted in the infrared and optical bands. The assumed fluctuations of the jet Lorentz factor are likely to be triggered by the variability of the accretion flow which is best traced by the X-ray emission. In the case of GX339-4 for which high quality and simultaneous multi-wavelength data are available, we performed simulations assuming that the fluctuation of the jet Lorentz factor have the same PSD as the observed X-ray PSD. The synthetic SED calculated under this assumption provides a remarkable match to the observed radio to IR SED. In this case the model also produces strong mid-infrared spectral variability that is comparable to that reported in this source.

We present an analysis of the parsec-scale jet structure of the quasar 4C+21.35 with a resolution of 0.1 milliarcseconds based on 63 epochs of Very Long Baseline Array observations at 43 GHz from 2007 June to 2014 May along with the Fermi LAT γ-ray light curve and multi-frequency optical photometric and polarimetric data. We find that the innermost jet of the quasar consists of a very compact core of size ~0.03 mas, as well as feature A1 located 0.16 ± 0.03 mas from the core. The distance of A1 remains fairly stable, but its position angle with respect to the core changes from -10 to +10 deg. We detect 4 superluminal knots in the inner jet with apparent speeds ranging from 10c to 20c. The first two components appeared in the jet during the high γ-ray state of the quasar from mid-2010 to early 2011, while the fourth knot appears to be connected with the γ-ray active state in late 2013 - early 2014. The first knot can be associated with the dramatic VHE flare in 2010 June and possesses an extreme Doppler factor ~60. We find that maxima in the γ-ray light curve coincide with epochs of interaction between the moving knots and the core and feature A1. This suggests that the core and A1 are recollimation shocks where γ-ray flares occur. The Chandra 0.5-6 keV image reveals the existence of X-ray emission in the kiloparsec scale jet of the quasar that can be explained via inverse Compton scattering off the cosmic microwave background by relativistic electrons if no deceleration occurs between the parsec- and kiloparsec-scale jets.

Many theoretical models of Active Galactic Nuclei (AGN) predict that the X-ray corona, lying above the black hole, constitutes the base of the X-ray jet. Thus, by studying the exact geometry of the close black hole environment, we can pinpoint the launching site of the jet. Detection of negative X-ray reverberation time delays (i.e. soft band X-ray variations lagging behind the corresponding hard band X-ray variations) can yield significant information about the geometrical properties of the AGN, such as the location of the X-ray source, as well as the physical properties of the the black hole, such as its mass and spin. In the frame-work of the lamp-post geometry, I present the first systematic X-ray time-lag modelling results of an ensemble of 12 AGN, using a fully general relativistic (GR) ray tracing approach for the estimation of the systems' response functions. By combing these state-of-the art GR response models with statistically innovative fitting routines, I derive the geometrical layout of the close BH environment for each source, unveiling the position of the AGN jet-base.

We numerically investigate the thermally unstable accretion discs around spinning black holes with different spins. We adopted an additional evolutionary viscosity equation to replace the standard alpha-prescription based on the results of two MHD simulations. We find an interesting oscillation when accretion switches to slim disc mode. The oscillation arises from the sonic point of accretion flow and propagates outwards. We mimic the bolometric light-curve and find a series of harmonics on its power spectrum. The frequency ratio of those harmonics is a regular integer series. The lowest frequency of the harmonics is identical to the prediction of trapped p-mode in QPO theory.

We briefly summarize the method of simulating Sgr A* polarized sub-mm spectra from the accretion flow and fitting the observed spectrum. The dynamical flow model is based on three-dimensional general relativistic magneto hydrodynamic simulations. Fully self-consistent radiative transfer of polarized cyclo-synchrotron emission is performed. We compile a mean sub-mm spectrum of Sgr A* and fit it with the mean simulated spectra. We estimate the ranges of inclination angle θ=42°–75°, mass accretion rate Ṁ=(1.4-7.0)×10−8M⊙year−1, and electron temperature Te=(3–4)×1010K at 6M. We discuss multiple caveats in dynamical modeling, which must be resolved to make further progress.