One of the main mechanisms that could drive mass outflows on parsec scales in AGN is thermal driving. The same X-rays that ionize and heat the plasma are also expected to make it thermally unstable. Indeed, it has been proposed that the observed clumpiness in AGN winds is caused by thermal instability (TI). While many studies employing time-dependent numerical simulations of AGN outflows have included the necessary physics for TI, none have so far managed to produce clumpiness. Here we present the first such clumpy wind simulations in 1-D and 2-D, obtained by simulating parsec scale outflows irradiated by an AGN. By combining an analysis of our extensive parameter survey with physical arguments, we show that the lack of clumps in previous numerical models can be attributed to the following three effects: (i) insufficient radiative heating or other physical processes that prevent the outflowing gas from entering the TI zone; (ii) the stabilizing effect of stretching (due to rapid radial acceleration) in cases where the gas enters the TI zone; and (iii) a flow speed effect: in circumstances where stretching is inefficient, the flow can still be so fast that it passes through the TI zone too quickly for perturbations to grow. Besides these considerations, we also find that a necessary condition to trigger TI in an outflow is for the pressure ionization parameter to decrease along a streamline once gas enters a TI zone.

The Galactic black hole X-ray binary MAXI J1820+070 had a bright outburst in 2018 when it became the second brightest X-ray source in the Sky. It was too bright for X-ray CCD instruments such as XMM-Newton and Chandra, but was well observed by photon-counting instruments such as NICER and NuSTAR. We report here on the discovery of an excess emission component during the soft state. It is best modelled with a blackbody spectrum in addition to the regular disk emission, modelled either as diskbb or kerrbb. Its temperature varies from about 0.9 to 1.1 keV which is about 30 to 80 per cent higher than the inner disc temperature of diskbb. Its flux varies between 4 and 12 percent of the disc flux. Simulations of magnetised accretion discs have predicted the possibility of excess emission associated with a non-zero torque at the Innermost Stable Circular Orbit (ISCO) about the black hole, which from other NuSTAR studies lies at about 5 gravitational radii or about 60 km (for a black hole mass is 8 M). In this case the emitting region at the ISCO has a width varying between 1.3 and 4.6 km and would encompass the start of the plunge region where matter begins to fall freely into the black hole.

We study how different opacity-temperature scalings affect the dynamical evolution of irradiated gas clouds using time-dependent, radiation-hydrodynamics (rad-HD) simulations. When clouds are optically thick, the bright side heats up and expands, accelerating the cloud via the rocket effect. Clouds that become more optically thick as they heat accelerate $\sim 35\%$ faster than clouds that become optically thin. An enhancement of $\sim 85\%$ in the acceleration can be achieved by having a broken powerlaw opacity profile, which allows the evaporating gas driving the cloud to become optically thin and not attenuate the driving radiation flux. We find that up to $\sim 2\%$ of incident radiation is re-emitted by accelerating clouds, which we estimate as the contribution of a single accelerating cloud to an emission or absorption line. Re-emission is suppressed by "bumps" in the opacity-temperature relation since these decrease the opacity of the hot, evaporating gas, primarily responsible for the re-radiation. If clouds are optically thin, they heat nearly uniformly, expand and form shocks. This triggers the Richtmyer-Meshkov instability, leading to cloud disruption and dissipation on thermal time-scales.

We use multi-frequency radiation hydrodynamics (rad-HD) to simulate radiative acceleration of a spherically symmetric stellar wind. We demonstrate the rad-HD capabilities of Athena++ for a series of test problems with multi-group radiation transfer. We then model the radiative transfer of a single spectral line through a spherically symmetric, isothermal, "CAK"-type line driven wind. We find that correctly accounting for the Doppler shift of the absorbed radiation, the force is well described by the analytic Sobolev line transfer in the supersonic parts of the solution where the flow is stationary and the effects of Abbott waves is negligible. Unlike in the analytic, steady-state solution re-radiation is important and leads to non-trivial radiation energy density and fluxes in the outer parts of the wind. We discuss a variety of applications to these multi-group methods that are currently computationally tractable.

Hercules X-1 is one of the best studied highly magnetised neutron star X-ray binaries with a wealth of archival data. We present the discovery of an ionised wind in its X-ray spectrum when the source is in the high state. The wind detection is statistically significant in most of the XMM-Newton observations, with velocities ranging from 200 to 1000 km/s. Observed features in the iron K band can be explained by both wind absorption or by a forest of iron emission lines. However, we also detect nitrogen, oxygen and neon absorption lines at the same systematic velocity in the high-resolution RGS grating spectra. The wind must be launched from the accretion disc, and could be the progenitor of the UV absorption features observed at comparable velocities, but the latter likely originate at significantly larger distances from the compact object. We find strong correlations between the ionisation level of the outflowing material and the ionising luminosity as well as the super-orbital phase. If the luminosity is driving the correlation, the wind could be launched by a combination of Compton heating and radiation pressure. If instead the super-orbital phase is the driver for the variations, the observations are likely scanning the wind at different heights above the warped accretion disc. If this is the case, we can estimate the wind mass outflow rate, corrected for the limited launching solid angle, to be roughly 70% of the mass accretion rate.

We study line driven winds for models with different radial intensity profiles: standard Shakura-Sunyaev radiating thin discs, uniform intensity discs and truncated discs where driving radiation is cutoff at some radius. We find that global outflow properties depend primarily on the total system luminosity but truncated discs can launch outflows with $\sim 2$ times higher mass flux and $\sim 50\%$ faster outflow velocity than non-truncated discs with the same total radiation flux. Streamlines interior to the truncation radius are largely unaffected and carry the same momentum flux as non-truncated models whereas those far outside the truncation radius effectively carry no outflow because the local radiation force is too weak to lift matter vertically away from the disc. Near the truncation radius the flow becomes more radial, due to the loss of pressure/radiation support from gas/radiation at larger radii. These models suggest that line driven outflows are sensitive to the geometry of the radiation field driving them, motivating the need for self-consistent disc/wind models.

We study the Rossby wave instability (RWI) in a thin accretion disc threaded by an initially toroidal magnetic field using the magnetohydrodynamics (MHD) code PLUTO. We find that for plasma $10^1 < \beta < 10^3$ , the growth rate and late-time density are suppressed, whereas for plasma $\beta > 10^3$, the magnetic field has negligible effect. The initially toroidal field is twisted inside the vortex and expelled at late times. This creates a radially directed field outside the vortex region, which can be understood via a simple kinematic model of a magnetic field in a rotating fluid and may be observable via polarised dust emission.

We study temporal variability of radiation driven winds using one dimensional, time dependent simulations and an extension of the classic theory of line driven winds developed by Castor Abbott and Klein. We drive the wind with a sinusoidally varying radiation field and find that after a relaxation time, determined by the propagation time for waves to move out of the acceleration zone of the wind, the solution settles into a periodic state. Winds driven at frequencies much higher than the dynamical frequency behave like stationary winds with time averaged radiation flux whereas winds driven at much lower frequencies oscillate between the high and low flux stationary states. Most interestingly, we find a resonance frequency near the dynamical frequency which results in velocity being enhanced or suppressed by a factor comparable to the amplitude of the flux variation. Whether the velocity is enhanced or suppressed depends on the relative phase between the radiation and the dynamical variables. These results suggest that a time-varying radiation source can induce density and velocity perturbations in the acceleration zones of line driven winds.

We numerically investigate the generation of a magnetic field in a protostellar disc via an $\alpha \Omega$-dynamo and the resulting magnetohydrodynamic (MHD) driven outflows. We find that for small values of the dimensionless dynamo parameter $\alpha_d$ the poloidal field grows exponentially at a rate $\sigma \propto \Omega_K \sqrt{\alpha_d}$, before saturating to a value $\propto \sqrt{\alpha_d}$. The dynamo excites dipole and octupole modes, but quadrupole modes are suppressed, because of the symmetries of the seed field. Initial seed fields too weak to launch MHD outflows are found to grow sufficiently to launch winds with observationally relevant mass fluxes of order $10^{-9} M_{\odot}/\rm{yr}$ for T Tauri stars. This suggests $\alpha \Omega$-dynamos may be responsible for generating magnetic fields strong enough to launch observed outflows.

We study mass outflows driven from accretion discs by radiation pressure due to spectral lines. To investigate non-axisymmetric effects, we use the Athena++ code and develop a new module to account for radiation pressure driving. In 2D, our new simulations are consistent with previous 2D axisymmetric solutions by Proga et al. who used the Zeus 2D code. Specifically, we find that the disc winds are time dependent, characterized by a dense stream confined to $\sim 45^{\circ}$ relative to the disc midplane and bounded on the polar side by a less dense, fast stream. Introducing a vertical, $\phi$-dependent, subsonic velocity perturbation in the disc midplane does not change the overall character of the solution but global outflow properties such as the mass, momentum and kinetic energy fluxes are altered by up to 100%. Non-axisymmetric density structures develop and persist mainly at the base of the wind. They are relatively small, and their densities can be a few times higher that the azimuthal average. The structure of the non-axisymmetric and axisymmetric solutions differ also in other ways. Perhaps most importantly from the observational point of view are the differences in the so called clumping factors, that serve as a proxy for emissivity due to two body processes. In particular, the spatially averaged clumping factor over the entire fast stream, while it is of a comparable value in both solutions, it varies about 10 times faster in the non-axisymmetric case.

We study non-axisymetric features of 3D line driven winds in the Sobolev approximation, where the optical depth is calculated using the full velocity gradient. We find that non-axisymmetric density features, so called clumps, form primarily at the base of the wind on super-Sobolev length scales. The density of clumps differs by a factor of $\sim 3$ from the azimuthal average, the magnitude of their velocity dispersion is comparable to the flow velocity and they produce $\sim 20\%$ variations in the column density. Clumps may be observable because differences in density produce enhancements in emission and absorption profiles or through their velocity dispersion which enhances line broadening.

We develop a general method for the self consistent calculation of the hydrodynamics of an astrophysical object irradiated by a radiation field with an arbitrary strength and spectral energy distribution (SED). Using the XSTAR photoionization code, we calculate heating and cooling rates as a function of gas photoionization parameter and temperature for several examples of SEDs: bremsstrahlung, blackbody, hard and soft state XRBs, Type 1 and Type 2 AGN. As an application of our method we study the hydrodynamics of 1-dimensional spherical winds heated by a uniform radiation field using the code Athena++. We find that in all cases explored a wind settles into a transonic, steady state. The wind evolves along the radiative heating equilibrium curve until adiabatic cooling effects become important and the flow departs from radiative equilibrium. If the flow is heated very rapidly, for example as in a thermally unstable regime, the corresponding column density of gas is low. Perhaps one of the most intriguing results of our work is the two stage acceleration of the wind that happens when there are two thermally unstable regions and the flux is relatively high. The efficiency with which the radiation field transfers energy to the wind is dependent on the SED of the external source, particularly the relative flux of soft X-rays. These results suggest that detailed photoionization calculations are essential not only to predict spectra but also to properly capture the flow dynamics.

We describe a Godunov-type magnetohydrodynamic (MHD) code based on the Miyoshi and Kusano (2005) solver which can be used to solve various astrophysical hydrodynamic and MHD problems. The energy equation is in the form of entropy conservation. The code has been implemented on several different coordinate systems: 2.5D axisymmetric cylindrical coordinates, 2D Cartesian coordinates, 2D plane polar coordinates, and fully 3D cylindrical coordinates. Viscosity and diffusivity are implemented in the code to control the accretion rate in the disk and the rate of penetration of the disk matter through the magnetic field lines. The code has been utilized for the numerical investigations of a number of different astrophysical problems, several examples of which are shown.

Observations of jets from young stellar objects reveal the asymmetric outflows from some sources. A large set of $2.5$D MHD simulations has been carried out for axisymmetric viscous/diffusive disc accretion to rotating magnetized stars for the purpose of assessing the conditions where the outflows or jets are asymmetric relative to the equatorial plane. We consider initial magnetic fields that are symmetric about the equatorial plane and consist of a radially distributed field threading the disc (disc-field) and a stellar dipole field.({\bf 1}). For pure disc-fields the symmetry or asymmetry of the outflows is affected by the midplane plasma $\beta$ of the disc (where $\beta$ is the ratio of the plasma pressure to the magnetic pressure). For the low density discs with small plasma $\beta$ values, outflows are observed to be symmetric to within $10\%$ over timescales of hundreds of inner disc orbits. For the denser higher $\beta$ discs, the coupling of the upper and lower coronal plasmas is broken, and quasi-periodic field motion in the two hemispheres becomes different. This asymmetry leads to asymmetric episodic outflows. ({\bf 2.}) Accreting stars with a stellar dipole field and no disc-field exhibit episodic, two component outflows - a magnetospheric wind and an inner disc wind from somewhat larger radial distances. Both are characterized by similar velocity profiles but the magnetospheric wind has densities $\gtrsim 10$ times that of the disc wind. ({\bf 3}.)Adding a disc-field parallel to the stellar dipole field acts to enhance the magnetospheric winds but suppress the disc wind. ({\bf 4}.) In contrast, adding a disc-field which is anti-parallel to the stellar dipole field in the disc acts to suppress the magnetospheric and disc winds. Our simulations reproduce some key features of observations of asymmetric outflows of T Tauri stars.

Counter-rotating discs can arise from the accretion of a counter-rotating gas cloud onto the surface of an existing co-rotating disc or from the counter-rotating gas moving radially inward to the outer edge of an existing disc. At the interface, the two components mix to produce gas or plasma with zero net angular momentum which tends to free-fall towards the disc center. We discuss high-resolution axisymmetric hydrodynamic simulations of a viscous counter-rotating disc for cases where the two components are vertically separated and radially separated. The viscosity is described by an isotropic $\alpha-$viscosity including all terms in the viscous stress tensor. For the vertically separated components a shear layer forms between them. The middle of this layer free-falls to the disk center. The accretion rates are increased by factors $\sim 10^2-10^4$ over that of a conventional disc rotating in one direction with the same viscosity. The vertical width of the shear layer and the accretion rate are strongly dependent on the viscosity and the mass fraction of the counter-rotating gas. In the case of radially separated components where the inner disc co-rotates and the outer disc rotates in the opposite direction, a gap between the two components opens and closes quasi-periodically. The accretion rates are $\gtrsim 25$ times larger than those for a disc rotating in one direction with the same viscosity.

Observations of galaxies and models of accreting systems point to the occurrence of counter-rotating discs where the inner part of the disc ($r<r_0$) is co-rotating and the outer part is counter-rotating. This work analyzes the linear stability of radially separated co- and counter-rotating thin discs. The strong instability found is the supersonic Kelvin-Helmholtz instability. The growth rates are of the order of or larger than the angular rotation rate at the interface. The instability is absent if there is no vertical dependence of the perturbation. That is, the instability is essentially three-dimensional. The nonlinear evolution of the instability is predicted to lead to a mixing of the two components, strong heating of the mixed gas, and vertical expansion of the gas, and annihilation of the angular momenta of the two components. As a result the heated gas will free-fall towards the disc's center over the surface of the inner disc.

We have carried out and analyzed a set of axisymmetric MHD simulations of the evolution of a turbulent/diffusive accretion disc around an initially unmagnetized star. The disc is initially threaded by a weak magnetic field where the magnetic pressure is significantly less than the kinetic pressure in the disc. The viscosity and magnetic diffusivity are modelled by two "alpha" parameters, while the coronal region above the disc is treated using ideal MHD. The initial magnetic field is taken to consist of three poloidal field loops threading the disc. The motivation for this study is to understand the advection of disc matter and magnetic field by the turbulent/diffusive disc. At early times the innermost field loop twists and its field lines become open. The twisting of the opened field lines leads to the formation of both an inner collimated, magnetically-dominated jet, and at larger distances from the axis a matter dominated uncollimated wind. For later times, the strength of the magnetic field decreases owing to field reconnection and annihilation in the disc. For the early times, we have derived from the simulations both the matter accretion speed in the disc and the accretion speed of the magnetic field. We show that the derived matter accretion speed agrees approximately with the predictions of a model where the accretion speed is the sum of two terms, one due to the disc's viscosity (which gives a radial outflow of angular momentum in the disc), and a second due to the twisted magnetic field at the disc's surface (which gives a vertical outflow of angular momentum). For early times we find that the magnetic contribution is roughly twice the viscous contribution for the case where the alpha parameters are both equal to 0.1. At later times the magnetic contribution to the matter speed becomes small compared to the viscous contribution.

We explore perturbations about a Friedmann-Robertson-Walker background in Chern-Simons gravity. At large momenta one of the two circularly polarized tensor modes becomes ghostlike. We argue that nevertheless the theory does not exhibit classical runaway solutions, except possibly in the relativistic nonlinear regime. However, the ghost modes cause the vacuum state to be quantum mechanically unstable, with a decay rate that is naively infinite. The decay rate can be made finite only if one interprets the theory as an effective quantum field theory valid up to some momentum cutoff, which violates Lorentz invariance. By demanding that the energy density in photons created by vacuum decay over the lifetime of the Universe not violate observational bounds, we derive strong constraints on the two dimensional parameter space of the theory, consisting of the cutoff and the Chern-Simons mass.

Recent radio emission, polarization, and Faraday rotation maps of the radio jet of the galaxy 3C 303 have shown that one knot of this jet has a {\it galactic}-scale electric current of $\sim 3\times 10^{18}$ Amp\`ere flowing along the jet axis (Kronberg et al. 2011). We develop the theory of relativistic Poynting-flux jets which are modeled as a transmission line carrying a DC current $I_0$, having a potential drop $V_0$, and a definite impedance ${\cal Z}_0 =90(u_z/c)\Omega$, where $u_z$ is the bulk velocity of the jet plasma. The electromagnetic energy flow in the jet is ${\cal Z}_0 I_0^2$. The observed current in 3C 303 can be used to calculate the electromagnetic energy flow in this magnetically dominated jet. Time-dependent but not necessarily small perturbations of a Poynting-flux jet - possibly triggered by a gas cloud penetrating the jet - are described by "telegrapher's equations," which predict the propagation speed of disturbances and the effective wave impedance ${\cal Z}$. The disturbance of a Poynting jet by the cloud gives rise to localized dissipation in the jet which may explain the enhanced synchrotron radiation in the knots of the 3C 303 jet, and in the apparently stationary knot HST-1 in the jet from the nucleus of the galaxy M87 (Biretta et al. 1999).

We study the deflection of light in the background of a "wiggly" cosmic string, and investigate whether it is possible to detect cosmic strings by means of weak gravitational lensing. For straight strings without small-scale structure there are no signals. In the case of strings with small-scale structure leading to a local gravitational attractive force towards the string, there is a small signal, namely a preferential elliptical distortion of the shape of background galaxies in the direction corresponding to the projection of the string onto the sky. The signal can be statistically distinguished from the signal produced by a linear distribution of black holes by employing an ellipticity axis distribution statistic.