Researcher profile

Matthew Liska

Matthew Liska contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
7works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

7 published item(s)

preprint2026arXiv

Evidence for Low Universal Equilibrium Black Hole Spin in Luminous Magnetically Arrested Disks

Relativistic collimated outflows, or jets, provide a crucial mode of active galactic nucleus feedback. Although jets extract their energy from the black hole (BH) rotation, their effect on the BH spin is poorly understood. Because the spin controls radiative and mechanical BH feedback, lack of first-principles models for spin evolution limits our ability to interpret observations, including the recent LIGO-Virgo-KAGRA spin constraints. Particularly important are luminous disks, which rapidly grow and strongly torque their BHs. Jetless and weakly magnetized standard luminous disks spin up their BHs to near-maximum spin, $a_{eq,NT}=0.998$. However, sufficient large-scale vertical magnetic flux can cause the inner disk to enter a magnetically arrested disk (MAD) state, whose jets can efficiently extract BH rotational energy and significantly spin down the BH. Lowell et al. (2024) found that nonradiative, thick MADs spin down their BHs to very low $a_{eq,MAD}^{thick}=0.07$. Their analytic model predicted that luminous, thin MADs also spin down their BHs to low $a_{eq,MAD}^{thin}\sim0.3\text{-}0.5$. To test this prediction, we perform 3D general relativistic (radiation) magnetohydrodynamic (GR(R)MHD) simulations of MADs across a wide range of BH spin ($-0.9\le{}a\le0.99$) and disk thickness ($0.03\le{}h/r\le0.3$, which corresponds to Eddington ratio, $0.35\le{}\dot{m}/\dot{m}_{Edd}\le\infty$). We find that luminous, thin MADs ($0.03\le{}h/r\le0.1$) efficiently spin down their BHs to a low universal equilibrium spin value, $a_{eq,MAD}^{thin}\approx0.3$: a maximally spinning BH ($a=1$) spins down to $a=0.5$ after accreting just $25\%$ of its initial mass. Our results follow quadratic convergence, $a_{eq,MAD}^{fit}\simeq0.3-2.7(h/r)^2\to0.3$ as $h/r\to0$, which we attribute to the aggressive cooling that renders disk thermodynamics irrelevant and magnetic forces insensitive to thermal $h/r$.

preprint2024arXiv

Gas Morphology of Milky Way-like Galaxies in the TNG50 Simulation: Signals of Twisting and Stretching

We present an in-depth analysis of gas morphologies for a sample of 25 Milky Way-like galaxies from the IllustrisTNG TNG50 simulation. We constrain the morphology of cold, warm, hot gas, and gas particles as a whole using a Local Shell Iterative Method (LSIM) and explore its observational implications by computing the hard-to-soft X-ray ratio, which ranges between $10^{-3}$-$10^{-2}$ in the inner $\sim 50 \rm kpc$ of the distribution and $10^{-5}$-$10^{-4}$ at the outer portion of the hot gas distribution. We group galaxies into three main categories: simple, stretched, and twisted. These categories are based on the radial reorientation of the principal axes of the reduced inertia tensor. We find that a vast majority ($77\%$) of the galaxies in our sample exhibit twisting patterns in their radial profiles. Additionally, we present detailed comparisons between 1) the gaseous distributions belonging to individual temperature regimes, 2) the cold gas distributions and stellar distributions, and 3) the gaseous distributions and dark matter (DM) halos. We find a strong correlation between the morphological properties of the cold gas and stellar distributions. Furthermore, we find a correlation between gaseous distributions with DM halo that increases with gas temperature, implying that we may use the warm-hot gaseous morphology as a tracer to probe the DM morphology. Finally, we show gaseous distributions exhibit significantly more prolate morphologies than the stellar distributions and DM halos, which we hypothesize is due to stellar and AGN feedback.

preprint2022arXiv

Black hole flares: ejection of accreted magnetic flux through 3D plasmoid-mediated reconnection

Magnetic reconnection can power bright, rapid flares originating from the inner magnetosphere of accreting black holes. We conduct extremely high resolution ($5376\times2304\times2304$ cells) general-relativistic magnetohydrodynamics simulations, capturing plasmoid-mediated reconnection in a 3D magnetically arrested disk for the first time. We show that an equatorial, plasmoid-unstable current sheet forms in a transient, non-axisymmetric, low-density magnetosphere within the inner few Schwarzschild radii. Magnetic flux bundles escape from the event horizon through reconnection at the universal plasmoid-mediated rate in this current sheet. The reconnection feeds on the highly-magnetized plasma in the jets and heats the plasma that ends up trapped in flux bundles to temperatures proportional to the jet's magnetization. The escaped flux bundles can complete a full orbit as low-density hot spots, consistent with Sgr A$^{*}$ observations by the GRAVITY interferometer. Reconnection near the horizon produces sufficiently energetic plasma to explain flares from accreting black holes, such as the TeV emission observed from M87. The drop in mass accretion rate during the flare, and the resulting low-density magnetosphere make it easier for very high energy photons produced by reconnection-accelerated particles to escape. The extreme resolution results in a converged plasmoid-mediated reconnection rate that directly determines the timescales and properties of the flare.

preprint2022arXiv

Black hole to photosphere: 3D GRMHD simulations of collapsars reveal wobbling and hybrid composition jets

Long-duration $γ$-ray bursts (GRBs) accompany the collapse of massive stars and carry information about the central engine. However, no 3D models have been able to follow these jets from their birth by a black-hole (BH) to the photosphere. We present the first such 3D general-relativity magnetohydrodynamic simulations, which span over 6 orders of magnitude in space and time. The collapsing stellar envelope forms an accretion disk, which drags inwardly the magnetic flux that accumulates around the BH, becomes dynamically important and launches bipolar jets. The jets reach the photosphere at $\sim10^{12}$ cm with an opening angle $θ_j\sim6^\circ$ and a Lorentz factor $Γ_j\lesssim 30$, unbinding $\gtrsim90\%$ of the star. We find that (i) the disk-jet system spontaneously develops misalignment relative to the BH rotational axis. As a result, the jet wobbles with an angle $θ_t\sim12^\circ$, which can naturally explain quiescent times in GRB lightcurves. The effective opening angle for detection $θ_j+θ_t$ suggests that the intrinsic GRB rate is lower by an order of magnitude than standard estimates. This suggests that successful GRBs may be rarer than currently thought and emerge in only $\sim 0.1\%$ of supernovae Ib/c, implying that jets are either not launched or choked inside most supernova Ib/c progenitors. (ii) The magnetic energy in the jet decreases due to mixing with the star, resulting in jets with a hybrid composition of magnetic and thermal components at the photosphere, where $\sim 10\%$ of the gas maintains magnetization $σ\gtrsim 0.1$. This indicates that both a photospheric component and reconnection may play a role in the prompt emission.

preprint2022arXiv

Bridging Bondi and Event Horizon Scales: 3D GRMHD Simulations Reveal X-Shaped Radio Galaxy Morphology

X-shaped radio galaxies (XRGs) produce misaligned X-shaped jet pairs and make up $\lesssim10\%$ of radio galaxies. XRGs are thought to emerge in galaxies featuring a binary supermassive black hole ($\rm SMBH$), $\rm SMBH$ merger, or large-scale ambient medium asymmetry. We demonstrate that XRG morphology can naturally form without such special, preexisting conditions. Our 3D general-relativistic magnetohydrodynamic (GRMHD) simulation for the first time follows magnetized rotating gas from outside the $\rm SMBH$ sphere of influence of radius $R_{\rm B}$ to the $\rm SMBH$ of gravitational radius $R_{\rm g}$, at the largest scale separation $R_{\rm B}/R_{\rm g} = 10^3$ to date. Initially, our axisymmetric system of constant-density hot gas contains weak vertical magnetic field and rotates in an equatorial plane of a rapidly spinning $\rm SMBH$. We seed the gas with small-scale $2\%$-level pressure perturbations. Infalling gas forms an accretion disk, and the $\rm SMBH$ launches relativistically-magnetized collimated jets reaching well outside $R_{\rm B}$. Under the pressure of the infalling gas, the jets intermittently turn on and off, erratically wobble, and inflate pairs of cavities in different directions, resembling an X-shaped jet morphology. Synthetic X-ray images reveal multiple pairs of jet-powered shocks and cavities. Large-scale magnetic flux accumulates on the $\rm SMBH$, becomes dynamically important, and leads to a magnetically arrested disk state. The $\rm SMBH$ accretes at $2\%$ of the Bondi rate ($\dot{M}\simeq2.4\times10^{-3}M_{\odot}\,{\rm yr}^{-1}$ for M87*), and launches twin jets at $η=150\%$ efficiency. These jets are powerful enough ($P_{\rm jets}\simeq2\times10^{44}\,{\rm erg\,s}^{-1}$) to escape along the spin axis and end the short-lived jets state whose transient nature can account for the rarity of XRGs

preprint2022arXiv

On the jet-ejecta interaction in 3D GRMHD simulations of binary neutron star merger aftermath

Short $γ$-ray burst (sGRB) jets form in the aftermath of a neutron star merger, drill through disk winds and dynamical ejecta, and extend over four to five orders of magnitude in distance before breaking out of the ejecta. We present the first 3D general-relativistic magnetohydrodynamic sGRB simulations to span this enormous scale separation. They feature three possible outcomes: jet+cocoon, cocoon, and neither. Typical sGRB jets break out of the dynamical ejecta if (i) the bound ejecta's isotropic equivalent mass along the pole at the time of the BH formation is $ \lesssim10^{-4}~{\rm M_{\odot}} $, setting a limit on the delay time between the merger and BH formation, otherwise, the jets perish inside the ejecta and leave the jet-inflated cocoon to power a low-luminosity sGRB; (ii) the post-merger remnant disk contains strong large-scale vertical magnetic field, $\gtrsim10^{15}$ G; and (iii) if the jets are weak ($\lesssim10^{50}$ erg), the ejecta's isotropic equivalent mass along the pole must be small ($\lesssim10^{-2}~{\rm M_{\odot}}$). Generally, the jet structure is shaped by the early interaction with disk winds rather than the dynamical ejecta. As long as our jets break out of the ejecta, they retain a significant magnetization ($\lesssim1$), suggesting that magnetic reconnection is a fundamental property of sGRB emission. The angular structure of the outflow isotropic equivalent energy after breakout consistently features a flat core followed by a steep power-law distribution (slope $\gtrsim3$), similar to hydrodynamic jets. In the cocoon-only outcome, the dynamical ejecta broadens the outflow angular distribution and flattens it (slope $\sim1.5$).

preprint2021arXiv

Black hole to breakout: 3D GRMHD simulations of collapsar jets reveal a wide range of transients

We present a suite of the first 3D GRMHD collapsar simulations, which extend from the self-consistent jet launching by an accreting Kerr black hole (BH) to the breakout from the star. We identify three types of outflows, depending on the angular momentum, $ l $, of the collapsing material and the magnetic field, $ B $, on the BH horizon: (i) Subrelativistic outflow (low $ l $ and high $ B $), (ii) Stationary accretion shock instability (SASI; high $ l $ and low $ B $), (iii) Relativistic jets (high $ l $ and high $ B $). In the absence of jets, free-fall of the stellar envelope provides a good estimate for the BH accretion rate. Jets can substantially suppress the accretion rate, and their duration can be limited by the magnetization profile in the star. We find that progenitors with large (steep) inner density power-law indices ($ \gtrsim 2 $), face extreme challenges as gamma-ray burst (GRB) progenitors due to excessive luminosity, global time evolution in the lightcurve throughout the burst and short breakout times, inconsistent with observations. Our results suggest that the wide variety of observed explosion appearances (supernova/supernova+GRB/low-luminosity GRBs) and the characteristics of the emitting relativistic outflows (luminosity and duration) can be naturally explained by the differences in the progenitor structure. Our simulations reveal several important jet features: (i) strong magnetic dissipation inside the star, resulting in weakly magnetized jets by breakout that may have significant photospheric emission and (ii) spontaneous emergence of tilted accretion disk-jet flows, even in the absence of any tilt in the progenitor.