Jet opening angles of long gamma-ray bursts (lGRBs) appear to evolve in cosmic time, with lGRBs at higher redshifts being on average more narrowly beamed than those at lower redshifts. We examine the nature of this anti-correlation in the context of collimation by the progenitor stellar envelope. First, we show that the data indicate a strong correlation between gamma-ray luminosity and jet opening angle, and suggest this is a natural selection effect - only the most luminous GRBs are able to successfully launch jets with large opening angles. Then, by considering progenitor properties expected to evolve through cosmic time, we show that denser stars lead to more collimated jets; we argue that the apparent anti-correlation between opening angle and redshift can be accounted for if lGRB massive star progenitors at high redshifts have higher average density compared to those at lower redshifts. This may be viable for an evolving IMF - under the assumption that average density scales directly with mass, this relationship is consistent with the form of the IMF mass evolution suggested in the literature. The jet angle-redshift anti-correlation may also be explained if the lGRB progenitor population is dominated by massive stars at high redshift, while lower redshift lGRBs allow for a greater diversity of progenitor systems (that may fail to collimate the jet as acutely). Overall, however, we find both the jet angle-redshift anti-correlation and jet angle-luminosity correlation are consistent with the conditions of jet launch through, and collimation by, the envelope of a massive star progenitor.

MAXI J1535-571 is a Galactic black hole candidate X-ray binary that was discovered going into outburst in 2017 September. In this paper, we present comprehensive radio monitoring of this system using the Australia Telescope Compact Array (ATCA), as well as the MeerKAT radio observatory, showing the evolution of the radio jet during its outburst. Our radio observations show the early rise and subsequent quenching of the compact jet as the outburst brightened and then evolved towards the soft state. We constrain the compact jet quenching factor to be more than 3.5 orders of magnitude. We also detected and tracked (for 303 days) a discrete, relativistically-moving jet knot that was launched from the system. From the motion of the apparently superluminal knot, we constrain the jet inclination (at the time of ejection) and speed to $\leq 45^{\circ}$ and $\geq0.69$c, respectively. Extrapolating its motion back in time, our results suggest that the jet knot was ejected close in time to the transition from the hard intermediate state to soft intermediate state. The launching event also occurred contemporaneously with a short increase in X-ray count rate, a rapid drop in the strength of the X-ray variability, and a change in the type-C quasi-periodic oscillation (QPO) frequency that occurs $>$2.5 days before the first appearance of a possible type-B QPO.

Despite their importance in the star formation process, measurements of magnetic field strength in proto-planetary discs remain rare. While linear polarisation of dust and molecular lines can give insight into the magnetic field structure, only observations of the circular polarisation produced by Zeeman splitting provide a direct measurement of magnetic field strenghts. One of the most promising probes of magnetic field strengths is the paramagnetic radical CN. Here we present the first Atacama Large Millimeter/submillimeter Array (ALMA) observations of the Zeeman splitting of CN in the disc of TW~Hya. The observations indicate an excellent polarisation performance of ALMA, but fail to detect significant polarisation. An analysis of eight individual CN hyperfine components as well as a stacking analysis of the strongest (non-blended) hyperfine components yields the most stringent limits obtained so far on the magnetic field strength in a proto-planetary disc. We find that the vertical component of the magnetic field $|B_z|<0.8$~mG ($1\sigma$ limit). We also provide a $1\sigma$ toroidal field strength limit of $<30$~mG. These limits rule out some of the earlier accretion disc models, but remain consistent with the most recent detailed models with efficient advection. We detect marginal linear polarisation from the dust continuum, but the almost purely toroidal geometry of the polarisation vectors implies that his is due to radiatively aligned grains.

GX~339$-$4 is a black hole X-ray binary that is a key focus of accretion studies since it goes into outburst roughly every two-to-three years. Tracking of its radio, IR and X-ray flux during multiple outbursts reveals tight broadband correlations. The radio emission originates in a compact, self-absorbed jet, however the origin of the X-ray emission is still debated: jet base or corona? We fit 20 quasi-simultaneous radio, IR, optical and X-ray observations of GX~339$-$4 covering three separate outbursts in 2005, 2007, 2010--2011, with a composite corona + jet model, where inverse Compton emission from both regions contributes to the X-ray emission. Using a recently-proposed identifier of the X-ray variability properties known as power-spectral hue, we attempt to explain both the spectral and evolving timing characteristics, with the model. We find the X-ray spectra are best fit by inverse Compton scattering in a dominant hot corona ($kT_{\rm e}\sim$ hundreds of keV). However, radio and IR-optical constraints imply a non-negligible contribution from inverse Compton scattering off hotter electrons ($kT_{\rm e} \ge 511$~keV) in the base of the jets, ranging from a few up to $\sim50$\% of the integrated 3--100~keV flux. We also find that the physical properties of the jet show interesting correlations with the shape of the broadband X-ray variability of the source, posing intriguing suggestions for the connection between the jet and corona.

We report on the results of optical, near-infrared (NIR) and mid-infrared observations of the black hole X-ray binary candidate (BHB) MAXI J1535-571 during its 2017/2018 outburst. During the first part of the outburst (MJD 58004-58012), the source shows an optical-NIR spectrum that is consistent with an optically thin synchrotron power-law from a jet. After MJD 58015, however, the source faded considerably, the drop in flux being much more evident at lower frequencies. Before the fading, we measure a de-reddened flux density of $\gtrsim$100 mJy in the mid-infrared, making MAXI J1535-571 one of the brightest mid-infrared BHBs known so far. A significant softening of the X-ray spectrum is evident contemporaneous with the infrared fade. We interpret it as due to the suppression of the jet emission, similar to the accretion-ejection coupling seen in other BHBs. However, MAXI J1535-571 did not transition smoothly to the soft state, instead showing X-ray hardness deviations, associated with infrared flaring. We also present the first mid-IR variability study of a BHB on minute timescales, with a fractional rms variability of the light curves of $\sim 15-22 \%$, which is similar to that expected from the internal shock jet model, and much higher than the optical fractional rms ($\lesssim 7 \%$). These results represent an excellent case of multi-wavelength jet spectral-timing and demonstrate how rich, multi-wavelength time-resolved data of X-ray binaries over accretion state transitions can help refining models of the disk-jet connection and jet launching in these systems.

Relativistic plasma jets are observed in many accreting black holes. According to theory, coiled magnetic fields close to the black hole accelerate and collimate the plasma, leading to a jet being launched. Isolating emission from this acceleration and collimation zone is key to measuring its size and understanding jet formation physics. But this is challenging because emission from the jet base cannot be easily disentangled from other accreting components. Here, we show that rapid optical flux variations from a Galactic black-hole binary are delayed with respect to X-rays radiated from close to the black hole by ~0.1 seconds, and that this delayed signal appears together with a brightening radio jet. The origin of these sub-second optical variations has hitherto been controversial. Not only does our work strongly support a jet origin for the optical variations, it also sets a characteristic elevation of <~10$^3$ Schwarzschild radii for the main inner optical emission zone above the black hole, constraining both internal shock and magnetohydrodynamic models. Similarities with blazars suggest that jet structure and launching physics could potentially be unified under mass-invariant models. Two of the best-studied jetted black hole binaries show very similar optical lags, so this size scale may be a defining feature of such systems.

The paradigm in which magnetic fields play a crucial role in launching/collimating outflows in many astrophysical objects continues to gain support. However, semi-analytical models including the effect of magnetic fields on the dynamics and morphology of jets are still missing due to the intrinsic difficulties in integrating the equations describing a collimated, relativistic flow in the presence of gravity. Only few solutions have been found so far, due to the highly nonlinear character of the equations together with the need to blindly search for singularities. These numerical problems prevented a full exploration of the parameter space. We present a new integration scheme to solve r-self-similar, stationary, axisymmetric magnetohydrodynamics equations describing collimated, relativistic outflows crossing smoothly all the singular points (Alfven point and modified slow/fast points). For the first time, we are able to integrate from the disk mid-plane to downstream of the modified fast point. We discuss an ensemble of jet solutions, emphasising trends and features that can be compared to observables. We present, for the first time with a semi-analytical MHD model, solutions showing counter-rotation of the jet for a substantial fraction of its extent. We find diverse jet configurations with bulk Lorentz factors up to 10 and potential sites for recollimation between $10^3-10^7$ gravitational radii. Such extended coverage of the intervals of quantities, such as magnetic-to-thermal energy ratios at the base or the heights/widths of the recollimation region, makes our solutions suitable for application to many different systems where jets are launched.

In this work we discuss the recent criticism by Zdziarski of the maximal jet model derived in Falcke & Biermann (1995). We agree with Zdziarski that in general a jet's internal energy is not bounded by its rest-mass energy density. We describe the effects of the mistake on conclusions that have been made using the maximal jet model and show when a maximal jet is an appropriate assumption. The maximal jet model was used to derive a 1-D hydrodynamical model of jets in agnjet, a model that does multiwavelength fitting of quiescent/hard state X-ray binaries and low-luminosity active galactic nuclei. We correct algebraic mistakes made in the derivation of the 1-D Euler equation and relax the maximal jet assumption. We show that the corrections cause minor differences as long as the jet has a small opening angle and a small terminal Lorentz factor. We find that the major conclusion from the maximal jet model, the jet-disk symbiosis, can be generally applied to astrophysical jets. We also show that isothermal jets are required to match the flat radio spectra seen in low-luminosity X-ray binaries and active galactic nuclei, in agreement with other works.

The `fundamental plane of black hole accretion' (FP), a relation between the radio luminosities ($L_R$), X-ray luminosities ($L_X$), and masses ($M_{BH}$) of hard/quiescent state black hole binaries and low-luminosity active galactic nuclei, suggests some aspects of black hole accretion may be scale invariant. However, key questions still exist concerning the relationship between the inflow/outflow behaviour in the `classic' hard state and quiescence, which may impact this scaling. We show that the broadband spectra of A0620-00 and~\sgra~(the least luminous stellar mass/supermassive black holes on the FP) can be modelled simultaneously with a physically-motivated outflow-dominated model where the jet power and all distances are scaled by the black hole mass. We find we can explain the data of both A0620-00 and~\sgra~(in its non-thermal flaring state) in the context of two outflow-model scenarios: (1) a synchrotron-self-Compton dominated state in which the jet plasma reaches highly sub-equipartition conditions (for the magnetic field with respect to that of the radiating particles), and (2) a synchrotron dominated state in the fast-cooling regime in which particle acceleration occurs within the inner few gravitational radii of the black hole and plasma is close to equipartition. We show that it may be possible to further discriminate between models (1) and (2) through future monitoring of its submm/IR/X-ray emission, in particular via time lags between the variable emission in these bands.

The collimation and evolution of relativistic outflows in $\gamma$-ray bursts (GRBs) are determined by their interaction with the stellar envelope through which they travel before reaching the much larger distance where the energy is dissipated and $\gamma$-rays are produced. We consider the case of a Poynting flux dominated relativistic outflow and show that it suffers strong inverse-Compton (IC) scattering drag near the stellar surface and the jet is slowed down to sub-relativistic speed if its initial magnetization parameter ($\sigma_0$) is larger than about 10$^5$. If the temperature of the cocoon surrounding the jet were to be larger than about 10 keV, then an optically thick layer of electrons and positrons forms at the interface of the cocoon and the jet, and one might expect this pair screen to protect the interior of the jet from IC drag. However, the pair screen turns out to be ephemeral, and instead of shielding the jet it speeds up the IC drag on it. Although a high $\sigma_0$ jet might not survive its passage through the star, a fraction of its energy is converted to 1-100 MeV radiation that escapes the star and appears as a bright flash lasting for about 10 s.

We consider the radiative transfer problem in a plane-parallel slab of thermal electrons in the presence of an overcritical magnetic field. Under such conditions, the magnetic field behaves as a birefringent medium for the propagating photons, and the electromagnetic radiation is split into two polarization modes, ordinary and extraordinary, having different cross-sections. When the optical depth of the slab is large, the ordinary-mode photons are strongly Comptonized and the photon field is dominated by an isotropic component. The radiative transfer problem in strong magnetic fields presents many mathematical issues and analytical or numerical solutions can be obtained only under some given approximations. We investigate this problem both from the analytical and numerical point of view, providing a test of the previous analytical estimates and extending these results introducing numerical techniques. We consider the case of low temperature blackbody photons propagating in a sub-relativistic temperature plasma, which allows us to deal with a semi Fokker-Planck approximation of the radiative transfer equation. We report the numerical solution obtained for eigenvalues and eigenfunctions of the space operator, and the emerging Comptonization spectrum of the ordinary-mode photons for any eigenvalue of the space equation and for energies significantly less than the cyclotron energy, which is of the order of MeV for the intensity of the magnetic field here considered. We derived the specific intensity of the ordinary photons, under the approximation of large angle and large optical depth. These assumptions allow the equation to be treated using a diffusion-like approximation.

We present a review of the Supergiant Fast X-ray Transients (SFXT) Project, a systematic investigation of the properties of SFXTs with a strategy that combines Swift monitoring programs with outburst follow-up observations. This strategy has quickly tripled the available sets of broad-band data of SFXT outbursts, and gathered a wealth of out-of-outburst data, which have led us to a broad-band spectral characterization, an assessment of the fraction of the time these sources spend in each phase, and their duty cycle of inactivity. We present some new observational results obtained through our outburst follow-ups, as fitting examples of the exceptional capabilities of Swift in catching bright flares and monitor them panchromatically.

Supergiant fast X-ray transients (SFXTs) are a class of high-mass X-ray binaries with possible counterparts in the high energy gamma rays. The Swift SFXT Project has conducted a systematic investigation of the properties of SFTXs on timescales ranging from minutes to years and in several intensity states (from bright flares, to intermediate intensity states, and down to almost quiescence). We also performed broad-band spectroscopy of outbursts, and intensity-selected spectroscopy outside of outbursts. We demonstrated that while the brightest phase of the outburst only lasts a few hours, further activity is observed at lower fluxes for a remarkably longer time, up to weeks. Furthermore, we assessed the fraction of the time these sources spend in each phase, and their duty cycle of inactivity. We present the most recent results from our investigation. The spectroscopic and, most importantly, timing properties of SFXTs we have uncovered with Swift will serve as a guide in search for the high energy emission from these enigmatic objects.

We report on the results from observations of the most recent outbursts of XTE J1739-302 and IGR J17544-2619, which are considered to be the prototypes of the supergiant fast X-ray transient (SFXT) class. They triggered the Swift/BAT on 2011 February 22 and March 24, respectively, and each time a prompt Swift slew allowed us to obtain the rich broad-band data we present. The XRT light curves show the descending portion of very bright flares that reached luminosities of ~2x10^{36} and ~5x10^{36} erg/s, respectively. The broad-band spectra, when fit with the usual phenomenological models adopted for accreting neutron stars, yield values of both high energy cut-off and e-folding energy consistent with those obtained from previously reported outbursts from these sources. In the context of more physical models, the spectra of both sources can be well fitted either with a two-blackbody model, or with a single unsaturated Comptonization model. In the latter case, the model can be either a classical static Comptonization model, such as COMPTT, or the recently developed COMPMAG model, which includes thermal and bulk Comptonization for cylindrical accretion onto a magnetized neutron star. We discuss the possible accretion scenarios derived by the different models, and we also emphasize the fact that the electron density derived from the Comptonization models, in the regions where the X-ray spectrum presumably forms, is lower than that estimated using the continuity equation at the magnetospheric radius and the source X-ray luminosity, and we give some possible explanations.

Predicting the emerging X-ray spectra in several astrophysical objects is of great importance, in particular when the observational data are compared with theoretical models. To this aim, we have developed an algorithm solving the radiative transfer equation in the Fokker-Planck approximation when both thermal and bulk Comptonization take place. The algorithm is essentially a relaxation method, where stable solutions are obtained when the system has reached its steady-state equilibrium. We obtained the solution of the radiative transfer equation in the two-dimensional domain defined by the photon energy E and optical depth of the system tau using finite-differences for the partial derivatives, and imposing specific boundary conditions for the solutions. We treated the case of cylindrical accretion onto a magnetized neutron star. We considered a blackbody seed spectrum of photons with exponential distribution across the accretion column and for an accretion where the velocity reaches its maximum at the stellar surface and at the top of the accretion column, respectively. In both cases higher values of the electron temperature and of the optical depth tau produce flatter and harder spectra. Other parameters contributing to the spectral formation are the steepness of the vertical velocity profile, the albedo at the star surface, and the radius of the accretion column. The latter parameter modifies the emerging spectra in a specular way for the two assumed accretion profiles. The algorithm has been implemented in the XSPEC package for X-ray spectral fitting and is specifically dedicated to the physical framework of accretion at the polar cap of a neutron star with a high magnetic field (> 10^{12} G), which is expected to be typical of accreting systems such as X-ray pulsars and supergiant fast X-ray transients.