Search

We present the results of four magnetohydrodynamic simulations and one alpha-disk simulation of accretion disks in a compact binary system, neglecting vertical stratification and assuming a locally isothermal equation of state. We demonstrate that in the presence of net vertical field, disks that extend out to the 3:1 mean motion resonance grow eccentricity in full MHD in much the same way as in hydrodynamical disks. Hence turbulence due to the magnetorotational instability (MRI) does not impede the tidally-driven growth of eccentricity in any meaningful way. However, we find two important differences with alpha-disk theory. First, in MHD, eccentricity builds up in the inner disk with a series of episodes of radial disk breaking into two misaligned eccentric disks, separated by a region of circular orbits. Standing eccentric waves are often present in the inner eccentric disk. Second, the successful spreading of an accretion disk with MRI turbulence out to the resonant radius is nontrivial, and much harder than spreading an alpha-disk. This is due to the tendency to develop over-dense rings in which tidal torques overwhelm MRI transport and truncate the disk too early. We believe that the inability to spread the disk sufficiently was the reason why our previous attempt to excite eccentricity via the 3:1 mean motion resonance with MHD failed. Exactly how MHD disks successfully spread outward in compact binary systems is an important problem that has not yet been understood., Comment: 16 pages, 11 Figures, 1 table, accepted by ApJ

For reasons of computational constraint, most global ocean circulation models used for Earth System Modeling still rely on parameterizations of sub-grid processes, and limitations in these parameterizations affect the modeled ocean circulation and impact on predictive skill. An increasingly popular approach is to leverage machine learning approaches for parameterizations, regressing for a map between the resolved state and missing feedbacks in a fluid system as a supervised learning task. However, the learning is often performed in an `offline' fashion, without involving the underlying fluid dynamical model during the training stage. Here, we explore the `online' approach that involves the fluid dynamical model during the training stage for the learning of baroclinic turbulence and its parameterization, with reference to ocean eddy parameterization. Two online approaches are considered: a full adjoint-based online approach, related to traditional adjoint optimization approaches that require a `differentiable' dynamical model, and an approximately online approach that approximates the adjoint calculation and does not require a differentiable dynamical model. The online approaches are found to be generally more skillful and numerically stable than offline approaches. Others details relating to online training, such as window size, machine learning model set up and designs of the loss functions are detailed to aid in further explorations of the online training methodology for Earth System Modeling., Comment: 25 pages, 1 table, 8 figures

Physical quantities with long lifetimes have both theoretical significance in the study of quantum many-body systems and practical implications for quantum technologies. In this manuscript, we investigate the roles played by topological defects in the construction of quasi-conserved quantities, using as a prototypical example the Kramers-Wannier duality defect in a deformed 1d quantum transverse field Ising model. We construct the duality defect Hamiltonian in three different ways: half-chain Kramers-Wannier transformation, utilization of techniques in the Ising fusion category, and defect-modified weak integrability breaking deformation. The third method is also applicable for the study of generic integrable defects under weak integrability breaking deformations. We also work out the deformation of defect-modified higher charges in the model and study their slower decay behavior. Furthermore, we consider the corresponding duality defect twisted deformed Floquet transverse field Ising model, and investigate the stability of the isolated zero mode associated with the duality defect in the integrable Floquet Ising model, under such weak integrability breaking deformation., Comment: 20 pages, 21 figures

We investigate how the formation and structure of circumplanetary disks (CPDs) varies with planet mass and protoplanetary disk aspect ratio. Using static mesh refinement and a near-isothermal equation of state, we perform a small parameter survey of hydrodynamic simulations with parameters appropriate for disk-embedded protoplanets at moderate to large orbital radii. We find that CPD formation occurs along a continuum, with ``diskiness'' increasing smoothly with planetary mass and decreasing disk aspect ratio. As expected from disk hydrostatic equilibrium arguments, the transition from envelope-dominated to disk-dominated structures is determined to first order by the ratio of the planetary Hill sphere radius to the disk scale height, but planets need to be significantly super-thermal to host classical rotationally supported CPDs. The circularization radius of inflowing gas (as a fraction of the Hill sphere radius) shows an approximately quadratic power-law scaling with the ratio of planetary mass to the thermal mass. Compared to more physically complete radiation hydrodynamic simulations, our runs almost maximize the possibility for classical CPD formation, and hence define a plausible necessary condition for CPDs. The low abundance of detected CPDs in disks where planetary companions are inferred from substructure data may be due to a combination of the large scale height of the protoplanetary disk, and a low frequency of sufficiently massive protoplanets. Unless their CPDs cool below the local protoplanetary disk temperature, most of the wide-orbit giant planet population will be embedded in quasi-spherical envelopes that are hard to detect. Disks, and satellite systems, are more likely to form around smaller orbital separation planets., Comment: Submitted; comments are welcome

Fragmentation in a gravitationally unstable accretion disk can be an important pathway for forming stellar/planetary companions. To characterize quantitatively the condition for and outcome of fragmentation under realistic thermodynamics, we perform global 3D simulations of gravitationally unstable disks at various cooling rates and cooling types, including the first global simulations of gravitational instability that employ full radiation transport. We find that fragmentation is a stochastic process, with the fragment generation rate per disk area $p_{\rm frag}$ showing an exponential dependence on the parameter $\beta\equiv\Omega_K t_{cool}$, where $\Omega_K$ is the Keplerian rotation frequency and $t_{cool}$ is the average cooling timescale. Compared to a prescribed constant $\beta$, radiative cooling in the optically thin/thick regime makes $p_{\rm frag}$ decrease slower/faster in $\beta$; the critical $\beta$ corresponding to $\sim 1$ fragment per orbit is $\approx$3, 5, 2 for constant $\beta$, optically thin, and optically thick cooling, respectively. The distribution function of the initial fragment mass is remarkably insensitive to disk thermodynamics. Regardless of cooling rate and optical depth, the typical initial fragment mass is $m_{frag} \approx 40 M_{tot}h^3$, with $M_{tot}$ being the total (star+disk) mass and $h=H/R$ being the disk aspect ratio. Applying this result to typical Class 0/I protostellar disks, we find $m_{frag}\sim 20 M_J$, suggesting that fragmentation more likely forms brown dwarfs. Given the finite width of the $m_{frag}$ distribution, forming massive planets is also possible., Comment: 21 pages, 11 figures, accepted for publication in ApJ

In this chapter, we summarize recent progress on the properties of accretion disks when the accretion rate exceeds the so-called Eddington limit based on multi-dimensional radiation magnetohydrodynamic simulations. We first summarize the classical models that are used to describe the accretion disks in the super-Eddington regime with an emphasis on the key uncertainties in these models. Then we show that radiation-driven outflows are ubiquitously found by numerical simulations of super-Eddington accretion disks. Some key physical processes on energy transport inside the disk are also identified by numerical simulations. Radiative and mechanical output as a function of mass accretion rates, black hole mass, spin, and magnetic field topology are summarized. Applications of super-Eddington accretion disks to different astrophysical systems, particularly tidal disruption events, are also discussed., Comment: 22 pages, 4 figures, invited review chapter for the book "New Frontiers in GRMHD Simulations"

Massive stars may form in or be captured into AGN disks. Recent 1D studies employing stellar-evolution codes have demonstrated the potential for rapid growth of such stars through accretion up to a few hundred $M_\odot$. We perform 3D radiation hydrodynamic simulations of moderately massive stars' envelopes, in order to determine the rate and critical radius $R_{\rm crit}$ of their accretion process in an isotropic gas-rich environment in the absence of luminosity-driven mass loss. We find that in the ``fast-diffusion" regime where characteristic radiative diffusion speed $c/\tau$ is faster than the gas sound speed $c_s$, the accretion rate is suppressed by feedback from gravitational and radiative advection energy flux, in addition to the stellar luminosity. Alternatively, in the ``slow-diffusion" regime where $c/\tau

The ability to fast reset a qubit state is crucial for quantum information processing. However, to actively reset a qubit requires engineering a pathway to interact with a dissipative bath, which often comes with the cost of reduced qubit protection from the environment. Here, we present a novel multi-purpose architecture that enables fast reset and protection of superconducting qubits during control and readout. In our design, two on-chip diplexers are connected by two transmission lines. The high-pass branch provides a flat passband for convenient allocation of readout resonators above the qubit frequencies, which is preferred for reducing measurement-induced state transitions. In the low-pass branch, we leverage a standing-wave mode below the maximum qubit frequency for a rapid reset. The qubits are located in the common stopband to inhibit dissipation during coherent operations. We demonstrate resetting a transmon qubit from its first excited state to the ground state in 100 ns, achieving a residual population of 2.7%, mostly limited by the thermal effect. The reset time may be further shortened to 27 ns by exploiting the coherent population inversion effect. We further extend the technique to resetting the qubit from its second excited state. Our approach promises scalable implementation of fast reset and qubit protection during control and readout, adding to the toolbox of dissipation engineering., Comment: 7 pages, 4 figures

Superconducting qubits, a promising candidate for universal quantum computing, currently face limitations in chip size due to reproducibility, wiring complexity, and packaging modes. Distributed quantum modules offer a viable strategy for constructing larger quantum information processing systems, though universal quantum gate operations between remote qubits have yet to be realized. Here, we demonstrate high-fidelity perfect entanglers between two remote superconducting quantum devices over 30 cm distance, leveraging the standing-wave modes in the coaxial cable connecting them. We achieve cross-entropy benchmarking (XEB) fidelities of $(99.15 \pm 0.02)\%$ and $(98.04 \pm 0.04)\%$ for CNOT and CZ gates, respectively, which are more efficient and universal than existing state transfer or feedback-based protocols. This advancement significantly enhances the feasibility of universal distributed quantum information processing, essential for the future development of large-scale quantum systems.

We use Athena++ to study the hydrodynamics of repeated star-accretion disk collisions close to supermassive black holes, and discuss their implications for the origin of quasi-periodic eruptions (QPEs) and other repeating nuclear transients. We quantify the impact of the collisions on the stellar structure, the amount of stripped stellar debris, and the debris' orbital properties. We provide simple fitting functions for the stellar mass-loss per collision; the mass-loss is much larger after repeated collisions due to the dilute stellar atmosphere shock-heated in earlier collisions. The lifetime of the QPE-emitting phase set by stellar mass-loss in star-disk collision models for QPEs is thus at most ~100 years; it is shortest for eRO-QPE2, of order a few decades. The mass of the stripped stellar debris per collision and its orbital properties imply that currently observed QPEs are not powered by direct star-disk collisions but rather by collisions between the stellar debris liberated in previous collisions and the accretion disk (`circularization shocks'). We discuss how the hydrodynamics of this interaction can explain the diverse timing properties of QPEs including the regular timing of GSN 069 and eRO-QPE2 and the large flare-to-flare timing variations observed in eRO-QPE1. QPEs with recurrence times of many days, if observed, may have more regular timing., Comment: 21 pages, 11 figures, submitted to ApJ

Image deraining aims to improve the visibility of images damaged by rainy conditions, targeting the removal of degradation elements such as rain streaks, raindrops, and rain accumulation. While numerous single image deraining methods have shown promising results in image enhancement within the spatial domain, real-world rain degradation often causes uneven damage across an image's entire frequency spectrum, posing challenges for these methods in enhancing different frequency components. In this paper, we introduce a novel end-to-end Adaptive Frequency Enhancement Network (AFENet) specifically for single image deraining that adaptively enhances images across various frequencies. We employ convolutions of different scales to adaptively decompose image frequency bands, introduce a feature enhancement module to boost the features of different frequency components and present a novel interaction module for interchanging and merging information from various frequency branches. Simultaneously, we propose a feature aggregation module that efficiently and adaptively fuses features from different frequency bands, facilitating enhancements across the entire frequency spectrum. This approach empowers the deraining network to eliminate diverse and complex rainy patterns and to reconstruct image details accurately. Extensive experiments on both real and synthetic scenes demonstrate that our method not only achieves visually appealing enhancement results but also surpasses existing methods in performance., Comment: 8pages

To date, the dayside thermal structure of ultrahot Jupiters (UHJs) is generally considered to be inverted, but their nightside thermal structure has been less explored. Here we explore the impact of nightside thermal emission on high-resolution infrared transmission spectroscopy, which should not be neglected, especially for UHJs. We present a general equation for the high-resolution transmission spectrum that includes planetary nightside thermal emission. This provides a new way to infer the thermal structure of the planetary nightside with high-resolution transmission spectroscopy. Using the cross-correlation technique, we find evidence for the presence of an H$_2$O emission signature on the UHJ WASP-33b during the transit, indicating an inverted temperature structure on its nightside. Such a result suggests a stronger heat transport through the circulation than currently expected. An alternative explanation is that the rotating visible hemisphere during transit leads to the potential contribution of the limb and dayside atmospheres to the detected emission signature. In the future, the combination of high-resolution full-phase curve spectroscopic observations and general circulation models will hopefully solve this puzzle and provide a complete picture of the three-dimensional nature of the chemistry, circulation, and thermal structure of UHJs., Comment: 17 pages, 9 figures, 1 table. Accepted for publication in ApJL

In this paper, the discrete-time modified algebraic Riccati equation (MARE) is solved when the system model is completely unavailable. To achieve this, firstly a brand new iterative method based on the standard discrete-time algebraic Riccati equation (DARE) and its input weighting matrix is proposed to solve the MARE. For the single-input case, the iteration can be initialized by an arbitrary positive input weighting if and only if the MARE has a stabilizing solution; nevertheless a pre-given input weighting matrix of a sufficiently large magnitude is used to perform the iteration for the multi-input case when the characteristic parameter belongs to a specified subset. Benefit from the developed specific iteration structure, the Q-learning (QL) algorithm can be employed to subtly solve the MARE where only the system input/output data is used thus the system model is not required. Finally, a numerical simulation example is given to verify the effectiveness of the theoretical results and the algorithm.

The rising part of a tidal disruption event light curve provides unique insight into early emission and the onset of accretion. Various mechanisms are proposed to explain the pre-peak emission, including shocks from debris interaction and reprocessing of disk emission. We study the pre-peak emission and its influence on the gas circularization by a series of gray radiation hydrodynamic simulations with varying black hole mass. We find that given a super-Eddington fallback rate of 10\dot{M}_{Edd}, the stream-stream collision can occur multiple times and drive strong outflows of up to 9\dot{M}_{Edd}. By dispersing gas to \gtrsim 100rs, the outflow can delay gas circularization and leads to sub-Eddington accretion rates during the first few stream-stream collisions. The stream-stream collision shock and circularization shock can sustain a luminosity of ~10^{44}erg/s for days. The luminosity is generally sub-Eddington and shows a weak correlation with accretion rate at early time. The outflow is optically thick, yielding a reprocessing layer with a size of ~10^{14} cm and photospheric temperature of ~4\times10^{4}K., Comment: 21 Pages, 19 figures, submitted to ApJ, comments welcome

We have developed the world's first canopy height map of the distribution area of world-level giant trees. This mapping is crucial for discovering more individual and community world-level giant trees, and for analyzing and quantifying the effectiveness of biodiversity conservation measures in the Yarlung Tsangpo Grand Canyon (YTGC) National Nature Reserve. We proposed a method to map the canopy height of the primeval forest within the world-level giant tree distribution area by using a spaceborne LiDAR fusion satellite imagery (Global Ecosystem Dynamics Investigation (GEDI), ICESat-2, and Sentinel-2) driven deep learning modeling. And we customized a pyramid receptive fields depth separable CNN (PRFXception). PRFXception, a CNN architecture specifically customized for mapping primeval forest canopy height to infer the canopy height at the footprint level of GEDI and ICESat-2 from Sentinel-2 optical imagery with a 10-meter spatial resolution. We conducted a field survey of 227 permanent plots using a stratified sampling method and measured several giant trees using UAV-LS. The predicted canopy height was compared with ICESat-2 and GEDI validation data (RMSE =7.56 m, MAE=6.07 m, ME=-0.98 m, R^2=0.58 m), UAV-LS point clouds (RMSE =5.75 m, MAE =3.72 m, ME = 0.82 m, R^2= 0.65 m), and ground survey data (RMSE = 6.75 m, MAE = 5.56 m, ME= 2.14 m, R^2=0.60 m). We mapped the potential distribution map of world-level giant trees and discovered two previously undetected giant tree communities with an 89% probability of having trees 80-100 m tall, potentially taller than Asia's tallest tree. This paper provides scientific evidence confirming southeastern Tibet--northwestern Yunnan as the fourth global distribution center of world-level giant trees initiatives and promoting the inclusion of the YTGC giant tree distribution area within the scope of China's national park conservation.

Nonlinear Compton scattering driven by ultraintense lasers presents a promising avenue for enhancing the photon energy, brilliance, and setup compactness of $\gamma$-ray sources. However, a significant challenge lies in achieving a high polarization degree with commonly generated unpolarized electrons, thus addressing a longstanding puzzle in the field. Here we investigate the polarization dynamics of photons emitted by an unpolarized electron beam interacting with a counter-propagating ultraintense laser pulse numerically, and propose a novel method to generate highly polarized $\gamma$ rays via nonlinear Compton scattering with the aid of vacuum dichroism effect. Our simulations reveal that high-brilliance $\gamma$ rays with polarization beyond 90\% are feasible in a single-shot interaction, rivaling the highest achieved by any $\gamma$-ray sources to date, based on a developed Monte Carlo method incorporating polarization-resolved tree processes of nonlinear Compton scattering and Breit-Wheeler pair production and one-loop vacuum polarization. This generation method showcases an extraordinary ultra-high polarization degree and a user-friendly all-optical experimental setup, while harnessing the high photon energy and brilliance characteristic of nonlinear Compton scattering sources, thus making it of great potential for experimental applications.

Turbulence is crucial for protoplanetary disk dynamics, and Vertical Shear Instability (VSI) is a promising mechanism in outer disk regions to generate turbulence. We use Athena++ radiation module to study VSI in full and transition disks, accounting for radiation transport and stellar irradiation. We find that the thermal structure and cooling timescale significantly influence VSI behavior. The inner rim location and radial optical depth affect disk kinematics. Compared with previous vertically-isothermal simulations, our full disk and transition disks with small cavities have a superheated atmosphere and cool midplane with long cooling timescales, which suppresses the corrugation mode and the associated meridional circulation. This temperature structure also produces a strong vertical shear at $\mathrm{\tau_*}$ = 1, producing an outgoing flow layer at $\tau_* < 1$ on top of an ingoing flow layer at $\tau_* \sim 1$. The midplane becomes less turbulent, while the surface becomes more turbulent with effective $\alpha$ reaching $\sim10^{-2}$ at $\tau_* \lesssim$1. This large surface stress drives significant surface accretion, producing substructures. Using temperature and cooling time measured/estimated from radiation-hydro simulations, we demonstrate that less computationally-intensive simulations incorporating simple orbital cooling can almost reproduce radiation-hydro results. By generating synthetic images, we find that substructures are more pronounced in disks with larger cavities. The higher velocity dispersion at the gap edge could also slow particle settling. Both properties are consistent with recent Near-IR and ALMA observations. Our simulations predict that regions with significant temperature changes are accompanied by significant velocity changes, which can be tested by ALMA kinematics/chemistry observations., Comment: 37 pages, 22 figures, accepted for publication in ApJ

Noise in quantum information processing poses a significant obstacle to achieving precise results. Quantum error mitigation techniques are crucial for improving the accuracy of experimental expectation values in this process. In the experiments, it is commonly observed that some measured events violate certain principles, such as symmetry constraints. These events can be considered junk information and should be discarded in a post-selection process. In this work, we introduce a quantum error mitigation method named Self-Trained Quantum Noise Filter (SQNF), which leverages the junk information to differentiate errors from the experimental population distributions, thereby aiming to approximate the error-free distribution. Our numerical results demonstrate that the proposed method can significantly reduce the infidelity of population distributions compared to the traditional post-selection method. Notably, the infidelity reduction is achieved without additional experimental resource consumption. Our method is scalable and applicable to multi-qubit computing systems.

Correlated errors can significantly impact the quantum error correction, which challenges the assumption that errors occur in different qubits independently in both space and time. Superconducting qubits have been found to suffer correlated errors across multiple qubits, which could be attributable to ionizing radiations and cosmic rays. Nevertheless, the direct evidence and a quantitative understanding of this relationship are currently lacking. In this work, we propose to continuously monitor multi-qubit simultaneous charge-parity jumps to detect correlated errors and find that occur more frequently than multi-qubit simultaneous bit flips. Then, we propose to position two cosmic-ray muon detectors directly beneath the sample box in a dilution refrigerator and successfully observe the correlated errors in a superconducting qubit array triggered by muons. By introducing a lead shielding layer on the refrigerator, we also reveal that the majority of other correlated errors are primarily induced by gamma rays. Furthermore, we find the superconducting film with a higher recombination rate of quasiparticles used in the qubits is helpful in reducing the duration of correlated errors. Our results provide experimental evidence of the impact of gamma rays and muons on superconducting quantum computation and offer practical insights into mitigation strategies for quantum error correction. In addition, we observe the average occurrence rate of muon-induced correlated errors in our processor is approximately 0.40 min$^{-1}$cm$^{-2}$, which is comparable to the muon event rate detected by the muon detector with 0.506 min$^{-1}$cm$^{-2}$. This demonstrates the potential applications of superconducting qubit arrays as low-energy threshold sensors in the field of high-energy physics., Comment: 7 pages and 5 figures for the main text, 20 pages and 20 figures for the supplementary materials

We use local stratified shearing-box simulations with magnetic field-aligned thermal conduction to study an idealized model of the coupling between a cold, radiatively efficient accretion disc, and an overlying, hot, two-temperature corona. Evaporation of a cold disc by conduction from the hot corona has been proposed as a means of mediating the soft-to-hard state transitions observed in X-ray binary systems. We model the coronal plasma in our local disc patch as an MHD fluid subject to both free-streaming ion conduction and a parameterized cooling function that captures the collisional transfer of energy from hot ions to colder, rapidly cooling leptons. In all of our models, independent of the initial net vertical magnetic flux (NF) threading the disc, we find no evidence of disc evaporation. The ion heat flux into the disc is radiated away before conduction can heat the disc's surface layers. When an initial NF is present, steady-state temperature, density, and outflow velocities in our model coronae are unaffected by conduction. Instead of facilitating disc evaporation, thermal conduction is more likely to feed the disc with plasma condensing out of the corona, particularly in flows without NF. Our work indicates that uncertainties in the amount of NF threading the disc hold far greater influence over whether or not the disc will evaporate into a radiatively inefficient accretion flow compared to thermal conduction. We speculate that a change in net flux mediates disc truncation/evaporation., Comment: 14 pages, 7 figures, 1 table, 2 appendices, submitted to MNRAS

Quasar feedback is a key ingredient in shaping galaxy evolution. A rare population of extremely red quasars (ERQs) at $z=2-3$ are often associated with high-velocity [O III]$\lambda5008$ outflows and may represent sites of strong feedback. In this paper, we present an X-ray study of 50 ERQs to investigate the link between the X-ray and outflow properties of these intriguing objects. Using hardness ratio analysis, we confirm that the ERQs are heavily obscured systems with gas column density reaching $N_\mathrm{H}=10^{23-24}\,\mathrm{cm^{-2}}$. We identify 20 X-ray-non-detected ERQs at high mid-infrared luminosities of $\nu L_\mathrm{\nu,6\mu m}\gtrsim3\times10^{46}\,\mathrm{erg\,s^{-1}}$. By stacking the X-ray observations, we find that the non-detected ERQs are on average underluminous in X-rays by a factor of $\sim10$ for their mid-infrared luminosities. We consider such X-ray weakness to be due to both heavy gas absorption and intrinsic factors. Moreover, we find that the X-ray-weak sources also display higher-velocity outflows. One option to explain this trend is that weaker X-rays facilitate more vigorous line-driven winds, which then accelerate the \oiii-emitting gas to kpc-scales. Alternatively, super-Eddington accretion could also lead to intrinsic X-ray weakness and more powerful continuum-driven outflow., Comment: 24 pages, 10 figures, Accepted to ApJ

We investigate the dynamical interaction between cosmic rays (CRs) and the multiphase interstellar medium (ISM) using numerical magnetohydrodynamic (MHD) simulations with a two-moment CR solver and TIGRESS simulations of star-forming galactic disks. We previously studied transport of CRs within TIGRESS outputs using a "post-processing" approach, and we now assess the effects of the MHD backreaction to CR pressure. We confirm our previous conclusion that there are three quite different regimes of CR transport in multiphase ISM gas, while also finding that simulations with "live MHD" predict a smoother CR pressure distribution. The CR pressure near the midplane is comparable to other pressure components in the gas, but the scale height of CRs is far larger. Next, with a goal of understanding the role of CRs in driving galactic outflows, we conduct a set of controlled simulations of the extraplanar region above $z=500$ pc, with imposed boundary conditions flowing from the midplane into this region. We explore a range of thermal and kinematic properties for the injected thermal gas, encompassing both hot, fast-moving outflows, and cooler, slower-moving outflows. The boundary conditions for CR energy density and flux are scaled from the supernova rate in the underlying TIGRESS model. Our simulations reveal that CRs efficiently accelerate extra-planar material if the latter is mostly warm/warm-hot gas, in which CRs stream at the Alfv\'en speed and the effective sound speed increases as density decreases. In contrast, CRs have very little effect on fast, hot outflows where the Alfv\'en speed is small, even when the injected CR momentum flux exceeds the injected MHD momentum flux., Comment: Resubmitted to ApJ after minor revisions

It is well established that supermassive black hole (SMBH) feedback is crucial for regulating the evolution of massive, if not all, galaxies. However, modelling the interplay between SMBHs and their host galaxies is challenging due to the vast dynamic range. Previous simulations have utilized simple subgrid models for SMBH accretion, while recent advancements track the properties of the unresolved accretion disc, usually based on the thin $\alpha$-disc model. However, this neglects accretion in the radiatively inefficient regime, expected to occur through a thick disc for a significant portion of an SMBH's lifetime. To address this, we present a novel 'unified' accretion disc model for SMBHs, harnessing results from the analytical advection-dominated inflow-outflow solution (ADIOS) model and state-of-the-art GR(R)MHD simulations. Going from low to high Eddington ratios, our model transitions from an ADIOS flow to a thin $\alpha$-disc via a truncated disc, incorporating self-consistently SMBH spin evolution due to Lense-Thirring precession. Utilizing the moving mesh code AREPO, we perform simulations of single and binary SMBHs within gaseous discs to validate our model and assess its impact. The disc state significantly affects observable luminosities, and we predict markedly different electromagnetic counterparts in SMBH binaries. Crucially, the assumed disc model shapes SMBH spin magnitudes and orientations, parameters that gravitational wave observatories like LISA and IPTA are poised to constrain. Our simulations emphasize the importance of accurately modelling SMBH accretion discs and spin evolution, as they modulate the available accretion power, profoundly shaping the interaction between SMBHs and their host galaxies., Comment: Accepted for publication in MNRAS; 29 pages, 17 figures, 1 appendix; doi:10.1093/mnras/stae1422