Search

We use numerical hydrodynamics simulations to study line driven winds launched from an accreting alpha-disc. Building on previous work where the driving radiation field is static, we compute a time-dependent radiation flux from the local, variable accretion rate of the disc. We find that prior to the establishment of a steady state in the disc, variations of ~ 15% in disc luminosity correlate with variations of ~ 2-3 in the mass flux of the wind. After a steady state is reached, when luminosity variations drop to ~ 3%, these correlations vanish as the variability in the mass flux is dominated by the intrinsic variability of the winds. This is especially evident in lower luminosity runs where intrinsic variability is higher due to a greater prevalence of failed winds. The changing mass flux occurs primarily due to the formation of clumps and voids near the disc atmosphere that propagate out into the low velocity part of the flow, a process that can be influenced by local variations in disc intensity. By computing the normalised standard deviation of the mass outflow, we show that the impact of luminosity variations on mass outflow is more visible at higher luminosity. However, the absolute change in mass outflow due to luminosity increases is larger for lower luminosity models due to the luminosity-mass flux scaling relation becoming steeper. We further discuss implications for CVs and AGN and observational prospects., 12 pages, 7 figures, accepted for publication in MNRAS

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., Comment: 11 pages, 5 figures

Local thermal instability can plausibly explain the formation of multiphase gas in many different astrophysical environments, but the theory is only well understood in the optically thin limit of the equations of radiation hydrodynamics (RHD). Here we lay groundwork for transitioning from this limit to a full RHD treatment assuming a gray opacity formalism. We consider a situation where the gas becomes thermally unstable due to the hardening of the radiation field when the main radiative processes are free-free cooling and Compton heating. We identify two ways in which this can happen: (i) when the Compton temperature increases with time, through a rise in either the intensity or energy of a hard X-ray component; and (ii) when attenuation reduces the flux of the thermal component so that the Compton temperature increases with depth through the slab. Both ways likely occur in the broad line region of active galactic nuclei where columns of gas can be ionization bounded. In such instances where attenuation is significant, thermal equilibrium solution curves become position-dependent and it no longer suffices to assess the stability of an irradiated column of gas at all depths using a single equilibrium curve. We demonstrate how to analyze a new equilibrium curve -- the attenuation curve -- for this purpose, and we show that by Field's instability criterion, a negative slope along this curve indicates that constant density slabs are thermally unstable whenever the gas temperature increases with depth., Comment: 10 pages, 1 figure, version slightly revised and accepted to ApJ Letters

Acceleration can change the ionization of X-ray irradiated gas to the point that the gas becomes thermally unstable. Cloud formation, the expected outcome of thermal instability (TI), will be suppressed in a dynamic flow, however, due to the stretching of fluid elements that accompanies acceleration. It is therefore unlikely that cloud formation occurs during the launching phase of a supersonic outflow. In this paper, we show that the most favorable conditions for dynamical TI in highly supersonic outflows are found at radii beyond the acceleration zone, where the growth rate of entropy modes is set by the linear theory rate for a static plasma. This finding implies that even mildly relativistic outflows can become clumpy, and we explicitly demonstrate this using hydrodynamical simulations of ultrafast outflows. We describe how the continuity and heat equations can be used to appreciate another impediment (beside mode disruption due to the stretching) to making an outflow clumpy: background flow conditions may not allow the plasma to enter a TI zone in the first place. The continuity equation reveals that both impediments are in fact tightly coupled, yet one is easy to overcome. Namely, time variability in the radiation field is found to be a robust means of placing gas in a TI zone. We further show how the ratio of the dynamical and thermal timescales enters linear theory; the heat equation reveals how this ratio depends on the two processes that tend to remove gas from a TI zone -- adiabatic cooling and heat advection., Accepted version to appear in ApJ (now 20 pages - Section 2 is somewhat extended/reorganized; Section 4 includes a discussion of the effects of cosmic rays). Simulations viewable at https://trwaters.github.io/dynamicalTI/

The warm absorbers observed in more than half of all nearby active galactic nuclei (AGN) are tracers of ionized outflows located at parsec scale distances from the central engine. If the smallest inferred ionization parameters correspond to plasma at a few $10^4$~K, then the gas undergoes a transition from being bound to unbound provided it is further heated to $\sim 10^6$~K at larger radii. Dannen et al. recently discovered that under these circumstances, thermally driven wind solutions are unsteady and even show very dense clumps due to thermal instability. To explore the observational consequences of these new wind solutions, we compute line profiles based on the one-dimensional simulations of Dannen et al. We show how the line profiles from even a simple steady state wind solution depend on the ionization energy (IE) of absorbing ions, which is a reflection of the wind ionization stratification. To organize the diversity of the line shapes, we group them into four categories: weak Gaussians, saturated boxy profiles with and without an extended blue wing, and broad weak profiles. The lines with profiles in the last two categories are produced by ions with the highest IE that probe the fastest regions. Their maximum blueshifts agree with the highest flow velocities in thermally unstable models, both steady state and clumpy versions. In contrast, the maximum blueshifts of the most high IE lines in thermally stable models can be less than half of the actual solution velocities. Clumpy solutions can additionally imprint distinguishable absorption troughs at widely separated velocities., Accepted for publication in ApJ, 20 pages, 15 figures, 3 tables, videos and FITS files of line profiles available on: http://www.physics.unlv.edu/astro/clumpywindsims-lps.html

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., Accepted to MNRAS. 23 pages, 16 figures and 7 tables

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., 7 pages, 3 figures

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., Comment: 11 pages, 7 figures

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., 10 pages, 8 figures

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., 12 pages, 9 figures, 2 tables

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 < ��< 10^3$ , the growth rate and late-time density are suppressed, whereas for plasma $��> 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., 10 pages, 10 figures

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., Comment: 11 pages, 13 figures, 1 table

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

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., Comment: 23 pages, 17 figures, submitted to New Astronomy

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., Comment: 16 pages, 17 figures, v2 has additional discussion regarding observations

A brief review of the origin of jets from disc-accreting rotating magnetized stars is given. In most models, the interior of the disc is characterized by a turbulent viscosity and magnetic diffusivity (‘alpha’ discs) whereas the coronal region outside the disc is treated using ideal magnetohydrodynamics (MHD). Extensive MHD simulations have established the occurrence of long-lasting outflows in the case of both slowly and rapidly rotating stars. (1) Slowly rotating stars exhibit a new type of outflow, conical winds. Conical winds are generated when stellar magnetic flux is bunched up by the inward motion of the accretion disc. Near their region of origin, the winds have a thin conical shell shape with half opening angle of ∼30∘. At large distances, their toroidal magnetic field collimates the outflow forming current carrying, matter dominated jets. These winds are predominantly magnetically and not centrifugally driven. About 10-30% of the disc matter from the inner disc is launched in the conical wind. Conical winds may be responsible for episodic as well as long lasting outflows in different types of stars. (2) Rapidly rotating stars in the ‘propeller regime’ exhibit twocomponent outflows. One component is similar to the matter dominated conical wind, where a large fraction of the disc matter may be ejected in this regime. The second component is a high-velocity, low-density magnetically dominated axial jet where matter flows along the open polar field lines of the star. The axial jet has a mass flux of about 10% that of the conical wind, but its energy flux, due to the Poynting flux, can be as large as for the conical wind. The jet’s magnetically dominated angular momentum flux causes the star to spin down rapidly. Propeller-driven outflows may be responsible for protostellar jets and their rapid spin-down.When the artificial requirement of symmetry about the equatorial plane is dropped, the conical winds are found to come alternately from one side of the disc and then the other, even for the case where the stellar magnetic field is a centered axisymmetric dipole.Recent MHD simulations of disc accretion to rotating stars in the propeller regime have been done with no turbulent viscosity and no diffusivity. The strong turbulence observed is due to the magneto-rotational instability. This turbulence drives accretion in the disc and leads to episodic conical winds and jets.

Observations of galaxies and models of accreting systems point to the occurrence of counter-rotating discs where the inner part of the disc ($r

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., 8 pages, 2 figures; final published version