Source author record

Sergei Dyda

Sergei Dyda appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

9works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

9 published item(s)

preprint2022arXiv

Dynamical thermal instability in highly supersonic outflows

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.

preprint2022arXiv

Thermal instability in radiation hydrodynamics: instability mechanisms, position-dependent S-curves, and attenuation curves

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.

preprint2020arXiv

Clumpy AGN outflows due to thermal instability

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.

preprint2020arXiv

Lense-Thirring Precession of Misaligned Discs I

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.

preprint2019arXiv

Effects of Opacity Temperature Dependence on Radiatively Accelerated Clouds

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.

preprint2015arXiv

Asymmetric MHD Outflows/Jets from Accreting T Tauri Stars

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.

preprint2014arXiv

Counter-Rotating Accretion Discs

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.

preprint2014arXiv

Kelvin-Helmholtz Instability of Counter-Rotating Discs

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.

preprint2012arXiv

Vacuum Instability in Chern-Simons Gravity

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.