-
Position dependent radiation fields near accretion disks
Authors:
Kara Smith,
Daniel Proga,
Randall Dannen,
Sergei Dyda,
Tim Waters
Abstract:
In disk wind models for active galactic nuclei (AGN) outflows, high-energy radiation poses a significant problem wherein the gas can become overionized, effectively disabling what is often inferred to be the largest force acting on the gas: the radiation force due to spectral line opacity. Calculations of this radiation force depend on the magnitude of ionizing radiation, which can strongly depend…
▽ More
In disk wind models for active galactic nuclei (AGN) outflows, high-energy radiation poses a significant problem wherein the gas can become overionized, effectively disabling what is often inferred to be the largest force acting on the gas: the radiation force due to spectral line opacity. Calculations of this radiation force depend on the magnitude of ionizing radiation, which can strongly depend on the position above a disk where the radiation is anisotropic. As our first step to quantify the position and direction dependence of the radiation field, we assumed free streaming of photons and computed energy distributions of the mean intensity and components of flux as well as energy-integrated quantities such as mean photon energy. We find a significant dependence of radiation field properties on position, but this dependence is not necessarily the same for different field quantities. A key example is that the mean intensity is much softer than the radial flux at many points near the disk. Because the mean intensity largely controls ionization, this softening decreases the severity of the overionization problem. The position dependence of mean intensity implies the position dependence of gas opacity, which we illustrate by computing the radiation force a fluid element feels in an accelerating wind. We find that in a vertical accelerating flow, the force due to radiation is not parallel to the radiation flux. This misalignment is due to the force's geometric weighting by both the velocity field's directionality and the position dependence of the mean intensity.
△ Less
Submitted 24 April, 2024;
originally announced April 2024.
-
Line Driven Instabilities due to Continuum Radiation Transport in Stellar Winds
Authors:
Sergei Dyda,
Shane W. Davis
Abstract:
We study line driven stellar winds using time-dependent radiation hydrodynamics where the continuum radiation couples to the gas via either a scattering or absorption opacity and there is an additional radiation force due to spectral lines that we model in the Sobolev approximation. We find that in winds with scattering opacities, instabilties tend to be suppressed and the wind reaches a steady st…
▽ More
We study line driven stellar winds using time-dependent radiation hydrodynamics where the continuum radiation couples to the gas via either a scattering or absorption opacity and there is an additional radiation force due to spectral lines that we model in the Sobolev approximation. We find that in winds with scattering opacities, instabilties tend to be suppressed and the wind reaches a steady state. Winds with absorption opacities are unstable and remain clumpy at late times. Clumps persist because they are continually regenerated in the subcritical part of the flow. Azimuthal gradients in the radial velocity distribution cause a drop in the radial radiation force and provide a mechanism for generating clumps. These clumps form on super-Sobolev scales, but at late times become Sobolev-length sized indicating that our radiation transfer model is breaking down. Inferring the clump distribution at late times therefore requires radiation-hydrodynamic modeling below the Sobolev scale.
△ Less
Submitted 28 November, 2023;
originally announced November 2023.
-
Time-Dependent AGN Disc Winds I -- X-ray Irradiation
Authors:
Sergei Dyda,
Shane W. Davis,
Daniel Proga
Abstract:
We study AGN line driven disc winds using time-dependent radiation hydrodynamics. The key criterion for determining wind launching is the coupling strength of the UV radiation field via the spectral lines of the gas. The strength of these lines in turn relies crucially on the gas ionization state, determined by the local X-ray intensity. We consider a suite of models where the central ionizing rad…
▽ More
We study AGN line driven disc winds using time-dependent radiation hydrodynamics. The key criterion for determining wind launching is the coupling strength of the UV radiation field via the spectral lines of the gas. The strength of these lines in turn relies crucially on the gas ionization state, determined by the local X-ray intensity. We consider a suite of models where the central ionizing radiation is affected by scattering, absorption and re-emission by the intervening gas. In a pure attenuation model, the disc launches an episodic wind, as previous studies have shown. Including scattering or re-emission tends to weaken the wind, lowering the mass flux and outflow velocity and if sufficiently dominant, suppressing the outflow entirely. However, the exponential nature of radiative attenuation means only a modest, factor of a few, increase in the absorption cross section can overcome the wind suppression due to scattering and re-emission. We find mass outflow rates of $\sim 20\%$ or more of the assumed inflow rate through the disk, indicating that radiation driven winds may significantly alter the structure of the accretion flow. The winds also supply a large, time-varying column of material above the nominal constant disk scale height, which will determine the geometry of reprocessed emission from the central source. Our results suggest the need for accurate photoionization modeling, radiation transport as well as accretion disc physics, to study their effects on the AGN disc winds
△ Less
Submitted 27 October, 2023;
originally announced October 2023.
-
Vertical wind structure in an X-ray binary revealed by a precessing accretion disk
Authors:
P. Kosec,
E. Kara,
A. C. Fabian,
F. Fürst,
C. Pinto,
I. Psaradaki,
C. S. Reynolds,
D. Rogantini,
D. J. Walton,
R. Ballhausen,
C. Canizares,
S. Dyda,
R. Staubert,
J. Wilms
Abstract:
The accretion of matter onto black holes and neutron stars often leads to the launching of outflows that can greatly affect the environments surrounding the compact object. In supermassive black holes, these outflows can even be powerful enough to dictate the evolution of the entire host galaxy, and yet, to date, we do not understand how these so-called accretion disk winds are launched - whether…
▽ More
The accretion of matter onto black holes and neutron stars often leads to the launching of outflows that can greatly affect the environments surrounding the compact object. In supermassive black holes, these outflows can even be powerful enough to dictate the evolution of the entire host galaxy, and yet, to date, we do not understand how these so-called accretion disk winds are launched - whether by radiation pressure, magnetic forces, thermal irradiation, or a combination thereof. An important means of studying disk winds produced near the central compact object is through X-ray absorption line spectroscopy, which allows us to probe outflow properties along a single line of sight, but usually provides little information about the global 3D disk wind structure that is vital for understanding the launching mechanism and total wind energy budget. Here, we study Hercules X-1, a unique, nearly edge-on X-ray binary with a warped accretion disk precessing with a period of about 35 days. This disk precession results in changing sightlines towards the neutron star, through the ionized outflow. We perform time-resolved X-ray spectroscopy over the precession phase and detect a strong decrease in the wind column density by three orders of magnitude as our sightline progressively samples the wind at greater heights above the accretion disk. The wind becomes clumpier as it rises upwards and expands away from the neutron star. Modelling the warped disk shape, we create a 2D map of wind properties. This unique measurement of the vertical structure of an accretion disk wind allows direct comparisons to 3D global simulations to reveal the outflow launching mechanism.
△ Less
Submitted 11 April, 2023;
originally announced April 2023.
-
Line Driven Winds from Variable Accretion Discs
Authors:
Anthony Kirilov,
Sergei Dyda,
Christopher S. Reynolds
Abstract:
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…
▽ More
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.
△ Less
Submitted 9 January, 2023;
originally announced January 2023.
-
The long stare at Hercules X-1 -- I. Emission lines from the outer disk, the magnetosphere boundary and the accretion curtain
Authors:
P. Kosec,
E. Kara,
A. C. Fabian,
F. Furst,
C. Pinto,
I. Psaradaki,
C. S. Reynolds,
D. Rogantini,
D. J. Walton,
R. Ballhausen,
C. Canizares,
S. Dyda,
R. Staubert,
J. Wilms
Abstract:
Hercules X-1 is a nearly edge-on accreting X-ray pulsar with a warped accretion disk, precessing with a period of about 35 days. The disk precession allows for unique and changing sightlines towards the X-ray source. To investigate the accretion flow at a variety of sightlines, we obtained a large observational campaign on Her X-1 with XMM-Newton (380 ks exposure) and Chandra (50 ks exposure) for…
▽ More
Hercules X-1 is a nearly edge-on accreting X-ray pulsar with a warped accretion disk, precessing with a period of about 35 days. The disk precession allows for unique and changing sightlines towards the X-ray source. To investigate the accretion flow at a variety of sightlines, we obtained a large observational campaign on Her X-1 with XMM-Newton (380 ks exposure) and Chandra (50 ks exposure) for a significant fraction of a single disk precession cycle, resulting in one of the best datasets taken to date on a neutron star X-ray binary. Here we present the spectral analysis of the High State high-resolution grating and CCD datasets, including the extensive archival data available for this famous system. The observations reveal a complex Fe K region structure, with three emission line components of different velocity widths. Similarly, the high-resolution soft X-ray spectra reveal a number of emission lines of various widths. We correct for the uncertain gain of the EPIC-pn Timing mode spectra, and track the evolution of these spectral components with Her X-1 precession phase and observed luminosity. We find evidence for three groups of emission lines: one originates in the outer accretion disk (10^5 RG from the neutron star). The second line group plausibly originates at the boundary between the inner disk and the pulsar magnetosphere (10^3 RG). The last group is too broad to arise in the magnetically-truncated disk and instead must originate very close to the neutron star surface, likely from X-ray reflection from the accretion curtain (~10^2 RG).
△ Less
Submitted 18 August, 2022;
originally announced August 2022.
-
Thermal instability in radiation hydrodynamics: instability mechanisms, position-dependent S-curves, and attenuation curves
Authors:
Daniel Proga,
Tim Waters,
Sergei Dyda,
Zhaohuan Zhu
Abstract:
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 be…
▽ More
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.
△ Less
Submitted 13 August, 2022; v1 submitted 8 July, 2022;
originally announced July 2022.
-
Dynamical thermal instability in highly supersonic outflows
Authors:
Tim Waters,
Daniel Proga,
Randall Dannen,
Sergei Dyda
Abstract:
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 ou…
▽ More
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.
△ Less
Submitted 8 April, 2022; v1 submitted 14 November, 2021;
originally announced November 2021.
-
On Synthetic Absorption Line Profiles of Thermally Driven Winds from Active Galactic Nuclei
Authors:
Shalini Ganguly,
Daniel Proga,
Tim Waters,
Randall C. Dannen,
Sergei Dyda,
Margherita Giustini,
Timothy Kallman,
John Raymond,
Jon Miller,
Paola Rodriguez Hidalgo
Abstract:
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. Danne…
▽ More
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.
△ Less
Submitted 15 May, 2021; v1 submitted 11 March, 2021;
originally announced March 2021.
-
Lense-Thirring Precession of Misaligned Discs I
Authors:
Sergei Dyda,
Christopher S. Reynolds
Abstract:
We study Lense-Thirring precession of inviscid and viscous misaligned $α-$discs around a black hole using a gravitomagnetic term in the momentum equation. For weak misalignments, $i \lesssim 10^{\circ}$, the discs behave like rigid bodies, undergoing the full suite of classical harmonic oscillator dynamics including, weak and critically damped motion (due to viscosity), precession (due to Lense-Th…
▽ More
We study Lense-Thirring precession of inviscid and viscous misaligned $α-$discs around a black hole using a gravitomagnetic term in the momentum equation. For weak misalignments, $i \lesssim 10^{\circ}$, the discs behave like rigid bodies, undergoing the full suite of classical harmonic oscillator dynamics including, weak and critically damped motion (due to viscosity), precession (due to Lense-Thirring torque) and nutation (due to apsidal precession). For strong misalignments, $i \gtrsim 30^{\circ}$, we find sufficiently thin, $h/r \lesssim 0.05$ discs break, form a gap and the inner and outer sub-discs evolve quasi independently apart from slow mass transfer. Assuming the sound speed sets the communication speed of warps in the disc, we can estimate the breaking radius by requiring that the inner sub-disc precesses like a rigid body. We explicitly show for the first time using a grid code that an Einstein potential is needed to reproduce the analytic properties of the inner disc edge and find disc breaking. At large inclination angles we find multiple disc breaks, consistent with recent GRMHD simulations of highly inclined discs. Our results suggest that the inclusion of a gravitomagnetic term and appropriate pseudo-Newtonian potential captures the important quantitative features of misaligned discs.
△ Less
Submitted 27 August, 2020;
originally announced August 2020.
-
The Soft State of the Black Hole Transient Source MAXI J1820+070: Emission from the Edge of the Plunge Region?
Authors:
A. C. Fabian,
D. J Buisson,
P. Kosec,
C. S. Reynolds,
D. R. Wilkins,
J. A. Tomsick,
D. J. Walton,
P. Gandhi,
D. Altamirano,
Z. Arzoumanian,
E. M. Cackett,
S. Dyda,
J. A. Garcia,
K. C. Gendreau,
B. W Grefenstette,
F. A. Harrison,
J. Homan,
E. Kara,
R. M. Ludlam,
J. M. Miller,
J. F. Steiner
Abstract:
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 mode…
▽ More
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.
△ Less
Submitted 22 February, 2020;
originally announced February 2020.
-
Clumpy AGN outflows due to thermal instability
Authors:
Randall Dannen,
Daniel Proga,
Tim Waters,
Sergei Dyda
Abstract:
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 incl…
▽ More
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.
△ Less
Submitted 9 April, 2020; v1 submitted 31 December, 2019;
originally announced January 2020.
-
Effects of Opacity Temperature Dependence on Radiatively Accelerated Clouds
Authors:
Sergei Dyda,
Daniel Proga,
Christopher S Reynolds
Abstract:
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…
▽ More
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.
△ Less
Submitted 12 December, 2019;
originally announced December 2019.
-
Line Driven Acceleration using Multi-Frequency Radiation Hydrodynamics
Authors:
Sergei Dyda,
Christopher S. Reynolds,
Yan-Fei Jiang
Abstract:
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…
▽ More
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.
△ Less
Submitted 20 November, 2019;
originally announced November 2019.
-
An ionised accretion disc wind in Hercules X-1
Authors:
P. Kosec,
A. C. Fabian,
C. Pinto,
D. J. Walton,
S. Dyda,
C. S. Reynolds
Abstract:
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…
▽ More
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.
△ Less
Submitted 13 November, 2019; v1 submitted 18 October, 2019;
originally announced October 2019.
-
Effects of Radiation Field Geometry on Line Driven Disc Winds
Authors:
Sergei Dyda,
Daniel Proga
Abstract:
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…
▽ More
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.
△ Less
Submitted 7 August, 2018;
originally announced August 2018.
-
Time Dependent Radiation Driven Winds
Authors:
Sergei Dyda,
Daniel Proga
Abstract:
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,…
▽ More
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.
△ Less
Submitted 22 May, 2018;
originally announced May 2018.
-
Non-Axisymmetric Line Driven Disc Winds II - Full Velocity Gradient
Authors:
Sergei Dyda,
Daniel Proga
Abstract:
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 thei…
▽ More
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.
△ Less
Submitted 10 February, 2018;
originally announced February 2018.
-
Rossby Vortices in Thin Magnetized Accretion Discs
Authors:
Loren Matilsky,
Sergei Dyda,
Richard V. E. Lovelace,
Patrick Lii
Abstract:
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 exp…
▽ More
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.
△ Less
Submitted 29 July, 2018; v1 submitted 23 November, 2017;
originally announced November 2017.
-
Non-Axisymmetric Line Driven Disc Winds I - Disc Perturbations
Authors:
Sergei Dyda,
Daniel Proga
Abstract:
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…
▽ More
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, $φ$-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.
△ Less
Submitted 10 February, 2018; v1 submitted 21 October, 2017;
originally announced October 2017.
-
Irradiation of Astrophysical Objects - SED and Flux Effects on Thermally Driven Winds
Authors:
Sergei Dyda,
Randall Dannen,
Tim Waters,
Daniel Proga
Abstract:
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, bla…
▽ More
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.
△ Less
Submitted 3 February, 2017; v1 submitted 13 October, 2016;
originally announced October 2016.
-
Numerical MHD Codes for Modeling Astrophysical Flows
Authors:
A. V. Koldoba,
G. V. Ustyugova,
P. S. Lii,
M. L. Comins,
S. Dyda,
M. M. Romanova,
R. V. E. Lovelace
Abstract:
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 c…
▽ More
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.
△ Less
Submitted 31 July, 2015;
originally announced August 2015.
-
Magnetic Field Amplification via Protostellar Disc Dynamos
Authors:
Sergei Dyda,
Richard V. E. Lovelace,
Galina V. Ustyugova,
Alexander V. Koldoba,
Ira Wasserman
Abstract:
We numerically investigate the generation of a magnetic field in a protostellar disc via an $αΩ$-dynamo and the resulting magnetohydrodynamic (MHD) driven outflows. We find that for small values of the dimensionless dynamo parameter $α_d$ the poloidal field grows exponentially at a rate $σ\propto Ω_K \sqrt{α_d}$, before saturating to a value $\propto \sqrt{α_d}$. The dynamo excites dipole and octu…
▽ More
We numerically investigate the generation of a magnetic field in a protostellar disc via an $αΩ$-dynamo and the resulting magnetohydrodynamic (MHD) driven outflows. We find that for small values of the dimensionless dynamo parameter $α_d$ the poloidal field grows exponentially at a rate $σ\propto Ω_K \sqrt{α_d}$, before saturating to a value $\propto \sqrt{α_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 $αΩ$-dynamos may be responsible for generating magnetic fields strong enough to launch observed outflows.
△ Less
Submitted 2 March, 2018; v1 submitted 30 July, 2015;
originally announced July 2015.
-
Asymmetric MHD Outflows/Jets from Accreting T Tauri Stars
Authors:
Sergei Dyda,
Richard V. E. Lovelace,
Galina V. Ustyugova,
Patrick S. Lii,
Marina M. Romanova,
Alexander V. Koldoba
Abstract:
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…
▽ More
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 $β$ of the disc (where $β$ is the ratio of the plasma pressure to the magnetic pressure). For the low density discs with small plasma $β$ values, outflows are observed to be symmetric to within $10\%$ over timescales of hundreds of inner disc orbits. For the denser higher $β$ 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.
△ Less
Submitted 16 March, 2015; v1 submitted 4 February, 2015;
originally announced February 2015.
-
Counter-Rotating Accretion Discs
Authors:
Sergei Dyda,
Richard V. E. Lovelace,
Galina V. Ustyugova,
Marina M. Romanova,
Alexander V. Koldoba
Abstract:
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-resolutio…
▽ More
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 $α-$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.
△ Less
Submitted 29 October, 2014; v1 submitted 24 August, 2014;
originally announced August 2014.
-
Kelvin-Helmholtz Instability of Counter-Rotating Discs
Authors:
Dan Quach,
Sergei Dyda,
Richard V. E. Lovelace
Abstract:
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 ar…
▽ More
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.
△ Less
Submitted 29 October, 2014; v1 submitted 24 August, 2014;
originally announced August 2014.
-
Advection of Matter and B-Fields in Alpha-Discs
Authors:
S. Dyda,
R. V. E. Lovelace,
G. V. Ustyugova,
P. S. Lii,
M. M. Romanova,
A. V. Koldoba
Abstract:
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…
▽ More
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.
△ Less
Submitted 1 April, 2013; v1 submitted 3 December, 2012;
originally announced December 2012.
-
Vacuum Instability in Chern-Simons Gravity
Authors:
Sergei Dyda,
Eanna E. Flanagan,
Marc Kamionkowski
Abstract:
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 unstab…
▽ More
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.
△ Less
Submitted 28 November, 2012; v1 submitted 23 August, 2012;
originally announced August 2012.
-
Relativistic Poynting-Flux Jets as Transmission Lines
Authors:
R. V. E. Lovelace,
S. Dyda,
P. P. Kronberg
Abstract:
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ère 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…
▽ More
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ère 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)Ω$, 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).
△ Less
Submitted 28 February, 2012;
originally announced February 2012.
-
Cosmic Strings and Weak Gravitational Lensing
Authors:
Sergei Dyda,
Robert H. Brandenberger
Abstract:
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 sign…
▽ More
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.
△ Less
Submitted 10 October, 2007;
originally announced October 2007.