2025-03-07 | | Total: 32
Fast radio bursts (FRBs) are millisecond-duration radio bursts with unidentified extra-galactic origin. Some FRBs exhibit mild magneto-ionic environmental variations, possibly attributed to plasma turbulence or geometric configuration variation in a binary system. Here we report an abrupt magneto-ionic environment variation of FRB 20220529, a repeating FRB from a disk galaxy at redshift 0.1839. Initially, its Faraday rotation measure (RM) was $21 \pm 96~{\rm rad~m^{-2}}$ over 17 months. In December 2023, it jumped to $1976.9~{\rm rad~m^{-2}}$, exceeding twenty times of the standard deviation of the previous RM variation, and returned to the typical values within two weeks. Such a drastic RM variation suggests a dense magnetized clump moving across the line of sight, possibly due to coronal mass ejection associated with a stellar flare. It indicates that the FRB likely has a companion star that produced the stellar flare.
Based on the iron and hydrogen similar elliptical morphologies of the type Ia supernova (SN Ia) remnant (SNR) 0509-67.5, I suggest that the ambient gas shaped the SN ejecta, that it is a remnant of an old planetary nebula, and that the explosion was more or less spherical. Adding that there is no observed stellar survivor in this SNR, I conclude that the SN Ia scenarios that best account for the explosion of this SN inside a planetary nebula (SNIP), are the lonely white dwarf scenarios, i.e., the core-degenerate (CD) and the double degenerate (DD) with merger to explosion delay (MED) time, i.e., the DD-MED scenario. Other scenarios encounter challenges, but I cannot completely rule them out. If true, the suggestion of SNR~0509-67.5 being a SNIP implies that the fraction of SNIPs might be higher than the previously estimated 50 percent of all normal SNe Ia. In the frame of the CD and the DD-MED scenarios, I attribute the observed double-shell structure in calcium to Rayleigh-Taylor instability during the explosion process: the caps of the Rayleigh-Taylor instability mushrooms form the outer shell. Instabilities during the explosion process might form clumps of some elements moving much faster than their mean velocity, possibly explaining fast-moving intermediate-mass elements.
We present simple analytic corrections to the standard blackbody fitting used for early kilonova emission. We consider a spherical, relativistically expanding shell that radiates thermally at a single temperature in its own rest frame. Due to relativistic effects, including Doppler boosting, time delay, and temperature evolution -- the observed temperature is smeared across different polar angles by approximately $\sim10\%$. While the observed spectrum remains roughly consistent with a single-temperature blackbody, neglecting relativistic effects leads to significant systematic inaccuracies: the inferred photospheric velocity and temperature are overestimated by up to $\sim50\%$ for mildly relativistic velocities. By applying our analytic corrections, these deviations are reduced to within $10\%$, even in cases where the photosphere is receding and cooling is considered. Applying our corrections to observed kilonovae (AT2017gfo and the thermal component of GRB211211A) reveals that standard blackbody fitting overestimated the inferred velocities and temperatures by $10\%-40\%$.
Blazars are interesting source candidates for astrophysical neutrino emission. Multi-messenger lepto-hadronic models based on proton-photon (p-gamma) interactions result in predictions for the neutrino spectra which are typically strongly peaked at high energies. In contrast, statistical analyses looking to associate blazars and neutrinos often assume a power-law spectral shape, putting the emphasis at lower energies. We aim to examine the impact of such spectral modelling assumptions on the associations of neutrinos with blazars. We use hierarchical_nu, a Bayesian framework for point source searches, and incorporate the theoretical predictions for neutrino spectra through a dedicated spectral model and priors on the relevant parameters. Our spectral model is based on the predictions presented in Rodrigues et al. (2024) for a selection of intermediate and high synchrotron peaked blazars that have been connected with IceCube neutrino alert events. We apply our model to the 10 years of publicly available IceCube data, focusing on the Northern hemisphere. Out of 29 source candidates, we find five sources, including TXS 0506+056, that have an association probability greater than 0.5 to at least one event. The p-gamma spectra typically lead to a lower overall number of associated events compared to the power-law case, but retain or even enhance strong associations to high-energy events. Our results demonstrate that including more information from theoretical predictions can allow for more interpretable source-neutrino connections.
Ultralight Dark Matter (ULDM) offers an alternative to cold dark matter, characterized by wavelike behavior on galactic scales. This contribution summarizes our work~\cite{Blas:2024duy} on how solitonic cores induced by ULDM affect gravitational waves (GWs), enabling their detection through next-generation observatories like LISA, Einstein Telescope and Cosmic Explorer. We show that continuous GWs from spinning neutron stars near the Galactic Center (GC) can probe solitons for \(m \lesssim 10^{-22}\, \mathrm{eV}\) and surpass current constraints on the ULDM couplings to the Standard Model for masses \(m \lesssim 10^{-20}\, \mathrm{eV}\). Additionally, GW modulation could reveal solitons in extragalactic systems, particularly in more massive halos, highlighting the broad potential of GW-based methods for uncovering ULDM properties.
Low-luminosity Active Galactic Nuclei (AGN) are believed to be surrounded by a collisionless, highly magnetized accretion flow. As a result, Particle-in-Cell simulations are the best tools to study the immediate vicinity of the event horizons of these supermassive black holes. We present a GPU-based general relativistic particle-in-cell (GRPIC) code framework called Aperture. Aperture is developed in C++, with compute kernels written in CUDA and HIP to take advantage of the massive acceleration modern GPUs enable. The code is organized in a fully modular way, allowing easy extensions to new physics problems. In this paper, we describe in detail the particle pusher, field solver, and charge-conserving current deposition algorithms employed in Aperture, and present test cases to validate their correctness. Then, we apply the code to study spark gaps and plasma injection in black hole magnetospheres. We find that the apparent location and time-evolution of the gap depend on the observer. Our results reconcile the previous conflicting findings from 1D and 2D simulations in the literature.
The detection of gravitational waves from binary black holes (BBHs) started the hunt for their pre-merger electromagnetic emission. In that respect, numerical simulations have been looking for the "smoking gun" signal that could help identify pre-merger systems. Here we study if any of the expected features of circumbinary discs, such as the periodic modulation from the orbiting "lump", could be used to identify pre-merger BBHs or if they could be easily confused with other systems. Indeed, while the timing feature associated with the "lump" seems to be present for a large part of the parameter space defined by the binary separation and mass ratio in circular binaries, it was recently proposed to form thanks to an instability occurring naturally at the edge of accretion discs around single black holes (SBH). In order to check if features of a circumbinary disc could be reproduced by a SBH system, we search for at least one SBH fit able to replicate the given synthetic observations of a circumbinary disc. We found that many of the features from a circumbinary disc can be reproduced by a SBH system with different masses, distances or inner disc positions. Interestingly, while we can always find a SBH model providing a good enough fit to the data, the presence of two variabilities, associated with the lump and the binary, or binary-lump beat, period, is a necessary condition for a wide range of BBH system parameters and should be used as a test to disqualify some BBH candidates.
Magnetars are the most magnetic objects in the Universe, serving as unique laboratories to test physics under extreme magnetic conditions that cannot be replicated on Earth. They were discovered in the late 1970s through their powerful X-ray flares, and were subsequently identified as neutron stars characterized by steady and transient emission across the radio, infrared, optical, X-ray, and gamma-ray bands. In this chapter, we summarize the current state of our experimental and theoretical knowledge on magnetars, as well as briefly discussing their relationship with supernovae, gamma-ray bursts, fast radio bursts, and the transient multi-band sky at large.
We investigate the axially symmetric accretion of low angular momentum hydrodynamic matter onto a rotating black hole. The gravitational field under consideration is assumed to be described by a pseudo-Newtonian Kerr potential. The accreting matter consists of different species defined by a relativistic equation of state with a variable adiabatic index.We construct and solve the hydrodynamical conservation equations governing such a flow, and find out the corresponding stationary integral solutions. We find that depending on the values of initial boundary conditions, accretion flow may exhibit multi-transonic behaviour, and a standing shock may form. We investigate, in minute detail, how the spin angular momentum of the black hole, as well as the composition of the accreting matter influence the dynamics of accretion flow and the astrophysics of shock formation in the aforementioned accreting black hole systems.
We present multiband observations and analysis of EP240801a, a low-energy, extremely soft gamma-ray burst (GRB) discovered on August 1, 2024 by the Einstein Probe (EP) satellite, with a weak contemporaneous signal also detected by Fermi/GBM. Optical spectroscopy of the afterglow, obtained by GTC and Keck, identified the redshift of $z = 1.6734$. EP240801a exhibits a burst duration of 148 s in X-rays and 22.3 s in gamma-rays, with X-rays leading by 80.61 s. Spectral lag analysis indicates the gamma-ray signal arrived 8.3 s earlier than the X-rays. Joint spectral fitting of EP/WXT and Fermi/GBM data yields an isotropic energy $E_{\gamma,\rm{iso}} = (5.57^{+0.54}_{-0.50})\times 10^{51}\,\rm{erg}$, a peak energy $E_{\rm{peak}} = 14.90^{+7.08}_{-4.71}\,\rm{keV}$, a fluence ratio $\rm S(25-50\,\rm{keV})/S(50-100\,\rm{keV}) = 1.67^{+0.74}_{-0.46}$, classifying EP240801a as an X-ray flash (XRF). The host-galaxy continuum spectrum, inferred using Prospector, was used to correct its contribution for the observed outburst optical data. Unusual early $R$-band behavior and EP/FXT observations suggest multiple components in the afterglow. Three models are considered: two-component jet model, forward-reverse shock model and forward-shock model with energy injection. Both three provide reasonable explanations. The two-component jet model and the energy injection model imply a relatively small initial energy and velocity of the jet in the line of sight, while the forward-reverse shock model remains typical. Under the two-component jet model, EP240801a may resemble GRB 221009A (BOAT) if the bright narrow beam is viewed on-axis. Therefore, EP240801a can be interpreted as an off-beam (narrow) jet or an intrinsically weak GRB jet. Our findings provide crucial clues for uncovering the origin of XRFs.
Swift J1727.8--1613 is a black hole X-ray binary that differs from other black hole X-ray binaries in that it has an extra hard component in addition to a reflection component. We perform spectral analysis with simultaneous Insight-HXMT, NICER and NuSTAR observations when the source was in the hard and hard intermediate states. For the presentation of the extra components, we investigate the correlation between the relxill parameters. We find that the correlation between the spectral index and the reflection fraction is consistent with MAXI J1820+070 when the jet dominates the reflection. This provides the second sample to have such a correlation during an outburst. Interestingly, when the reflection component is attributed to the corona, the spectral fit results in a small reflection fraction and the correlation between the spectral index and reflection fraction is in agreement with the overall trend built-in You et al. 2023 with a large sample of outbursts from other X-ray binaries. Hence Swift J1727.8--1613 turns out to be the first sample to bridge the MAXI J1820+070 to the majority of X-ray binaries according to the dual correlations observed between the spectral index and the reflection fraction. We speculate that a configuration of a jet plus a hot inner flow may account for this peculiar outburst behavior of Swift J1727.8--1613.
In this work, we numerically explore the dynamics of a point mass (e.g., a stellar-mass black hole) moving within a thin accretion disk of a massive object (i.e., a supermassive black hole) with three dimensional hydrodynamical simulations using \texttt{Athena++}. We are particularly interested in the regime that the Hill radius of the point mass is greater than the disk thickness, but the point mass is not sufficiently massive to open a ``gap" in a high viscosity disk. This parameter regime may be seen for stellar-mass objects moving within thin accretion disks of massive black holes, but less studied for the planet migration scenario in protoplanetary disks. We find that the disk migration may be significantly slower than the type I migration, depending on the surface density gradient of the disk. Furthermore, the circumsingle disk around the point mass plays an important role in damping the orbital eccentricity. If the gravitational interaction between the point mass and circumsingle disk material is turned off, the orbital eccentricity may be pumped to $\mathcal{O}(10^{-2})$ level. Because space-borne gravitational wave detectors such as LISA (Laser Interferometer Space Antenna) is able to measure the eccentricity of extreme mass-ratio inspirals to the level of $\mathcal{O}(10^{-5})$, this finding highlights the importance of understanding stellar-mass black hole feedback mechanisms which will modify the structures of the circumsingle disk, as they will impact LISA observables such as the eccentricity.
Relativistic jets from accreting black holes (BHs) radiate non-thermal emission which is highly variable in different time scales. Magnetic fields anchored to a rotating BH or accretion disc accelerate and collimate jets of the BH systems. Previous studies on black holes of different mass scales, including supermassive and stellar-mass black holes, only report flux quasi-periodic oscillations in radio, optical, X-ray and gamma-ray bands. No quasi-periodic variations in polarization have yet been detected in any black hole systems. Here, we report the first detection of GHz radio polarization oscillations in GRS 1915+105, which harbors a spinning stellar-mass BH with a relativistic jet. Our observations show that during the increasing phase of radio emission, linear polarization and flux exhibit similar oscillation periods of $\sim 17$ and $33$ seconds, and their variation patterns anti-correlate with each other. These rare, short-period oscillations in both polarization and flux would be important to understand instabilities and special dynamics in magnetized jets.
Supergiant Fast X-ray Transients (SFXT) are a sub-class of High Mass X-ray Binaries (HMXB) in which a compact object accretes part of the clumpy wind of the blue supergiant companion, triggering a series of brief, X-ray flares lasting a few kiloseconds. Currently, only about fifteen SFXTs are known. The EXTraS catalog provides the timing signatures of every source observed by the EPIC instrument on-board XMM-Newton. Among the most peculiar sources, in terms of variability, we selected 4XMM J181330.1-17511 (J1813). We analyzed all publicly available X-ray data pointed at the J1813 position to determine the source's duty cycle and to provide a comprehensive description of its timing and spectral behavior during its active phase. Additionally, we searched for the optical and infrared counterpart of the X-ray source in public databases and fitted its Spectral Energy Distribution (SED). The optical-to-MIR SED of J1813 is consistent with a highly-absorbed (A$_V\sim38$) B0 star at $\sim$10 kpc. During its X-ray active phase, the source is characterized by continuous $\sim$thousands seconds-long flares with peak luminosities (2-12 keV) ranging from $10^{34}$ to $4 \times 10^{35}$ erg s$^{-1}$. Its X-ray spectrum is consistent with a high-absorbed power-law model with N$_H \sim 1.8 \times 10^{23}$ cm$^{-2}$ and $\Gamma \sim 1.66$. No spectral variability was observed as a function of time or flux. J1813 is in a quiescent state $\sim$60\% of the time, with an upper-limit luminosity of $8 \times 10^{32}$ erg s$^{-1}$ (at 10 kpc), implying an observed long-term X-ray flux variability $>$500. The optical counterpart alone indicates J1813 is a HMXB. Its transient nature, duty cycle, the amplitude of observed X-ray variability, the shape and luminosity of the X-ray flares -- and the lack of known X-ray outbursts ($>10^{36}$ erg s$^{-1}$) -- strongly support the identification of J1813 as an SFXT.
In view of the great uncertainty of the equation of state (EOS) of high-density nuclear matter, establishing EOS-independent universal relations between global properties of neutron stars provides a practical way to constrain the unobservable or difficult-to-observe properties through astronomical observations. It is common to construct universal relations between EOS-dependent properties (e.g., moment of inertia, tidal deformation, etc.) or combined properties (e.g., compactness). Improving the precision of the universal relations may provide stricter constraint on the properties of neutron star. We find that in 3-dimensional space with mass and radius as the base coordinates, the points corresponding to a certain property of neutron star described by different EOSs are almost located in the same surface. Thus the universal relation between the property and the stellar mass-radius can be expressed through describing the surface. It is shown that the resulting universal relations have higher precisions. As an example, we construct high-precision universal relations for the moment of inertia, the $f$-mode frequency, and the dimensionless tidal deformation respect to the mass-radius. As the observational data of neutron star mass and radius from NICER grows in data and accuracy, these universal relations allow for more precise constraints on the unobservable or difficult-to-observe properties.
We present photometric and spectroscopic data of the type IIb supernova (SN) 2024abfo in NGC 1493 (at 11 Mpc). The ATLAS survey discovered the object just a few hours after the explosion, and observed a fast rise on the first day. Signs of the sharp shock break-out peak and the subsequent cooling phase are observed in the ultraviolet and the bluest optical bands in the first couple of days, while no peak is visible in the reddest filters. Subsequently, in analogy with normal SNe IIb, the light curve of SN 2024abfo rises again in all bands to the broad peak, with the maximum light reached around one month after the explosion. Its absolute magnitude at peak is $M_r=-16.5\pm0.1$ mag, making it a faint SN IIb. The early spectra are dominated by Balmer lines with broad P-Cygni profiles indicating ejecta velocity of 22,500 \kms. One month after the explosion, the spectra display a transition towards being He-dominated, though the H lines do not completely disappear, supporting the classification of SN 2024abfo as a relatively H-rich SN IIb. We identify the progenitor of SN 2024abfo in archival images of the Hubble Space Telescope, the Dark Energy Survey, and the XMM-Newton space telescope, in multiple optical filters. From its spectral energy distribution, the progenitor is consistent with being a white supergiant, Deneb-like star, with a photospheric temperature of 8110 K, a radius of 176 \Rsun, a luminosity of $\log(L/L_{\odot})=5.08$, having an initial mass of 15-16 \Msun. This detection supports an emerging trend of SN IIb progenitors being more luminous and hotter than SN II ones, and being primaries of massive binaries. Within the SN IIb class, fainter events such as SN 2024abfo tend to have cooler and more expanded progenitors than luminous SNe IIb.
The impact of neutrino flavor conversion on the supernova mechanism is yet to be fully understood. We present multi-energy and multi-angle solutions of the neutrino quantum kinetic equations in three flavors, taking into account neutrino advection and non-forward collisions with the background medium. Flavor evolution is explored within a spherically symmetric shell surrounding the region of neutrino decoupling in the interior of a core-collapse supernova, relying on the outputs of a spherically symmetric core-collapse supernova model with a mass of $18.6 M_\odot$. We select two representative post-bounce times: $t_{\rm pb} = 0.25$ s (no angular crossings are present and flavor conversion is triggered by slow collective effects) and $t_{\rm pb} = 1$ s (angular crossings trigger fast flavor instabilities). We find that flavor equipartition is achieved in the antineutrino sector between $\bar\nu_e$ and $\bar\nu_x = (\bar\nu_\mu + \bar\nu_\tau)/2$ for both post-bounce times. In the neutrino sector, flavor equipartition between $\nu_e$ and $\nu_x$ seems more likely at later post-bounce times, where the neutrino emission properties among different flavors tend to approach each other, but it is not a generic feature. The exponential growth of the $\nu_\mu$--$\nu_\tau$ asymmetry due to three-flavor effects is responsible for differences between the quasi-steady configurations obtained in the three-flavor solution and in the two-flavor approximation. This has consequences on the neutrino heating rate, which is generally larger when all three flavors are taken into account and can increase up to $30\%$ with respect to the case where flavor conversion is neglected.
We study particle acceleration in strongly turbulent pair plasmas using novel 3D Particle-in-Cell simulations, featuring particle injection from an external heat bath and diffusive escape. We demonstrate the formation of steady-state, nonthermal particle distributions with maximum energies reaching the Hillas limit. The steady state is characterized by the equilibration of plasma kinetic and magnetic pressures, which imposes upper limits on the acceleration rate. With growing cold plasma magnetization $\sigma_0$, nonthermal power-law spectra become harder, and the fraction of energy channeled into escaping cosmic rays increases. At $\sigma_0 \gtrsim 1$, the escaping cosmic rays amount to more than 50\% of the dissipated energy. Our method allows for kinetic studies of particle acceleration under steady-state conditions, with applications to a variety of astrophysical systems.
Binary population synthesis simulations allow detailed modelling of gravitational-wave sources from a variety of formation channels. These population models can be compared to the observed catalogue of merging binaries to infer the uncertain astrophysical input parameters describing binary formation and evolution, as well as relative rates between various formation pathways. However, it is computationally infeasible to run population synthesis simulations for all variations of uncertain input physics. We demonstrate the use of normalising flows to emulate population synthesis results and interpolate between astrophysical input parameters. Using current gravitational-wave observations of binary black holes, we use our trained normalising flows to infer branching ratios between multiple formation channels, and simultaneously infer common-envelope efficiency and natal spins across a continuous parameter range. Given our set of formation channel models, we infer the natal spin to be $0.04^{+0.04}_{-0.01}$, and the common-envelope efficiency to be $>3.7$ at 90% credibility, with the majority of underlying mergers coming from the common-envelope channel. Our framework allows us to measure population synthesis inputs where we do not have simulations, and better constrain the astrophysics underlying current gravitational-wave populations.
In classic time-dependent 1D accretion disk models, the inner radiation pressure dominated regime is viscously unstable. However, late-time observations of accretion disks formed in tidal disruption events (TDEs) do not exhibit evidence of such instabilities. The common theoretical response is to modify the viscosity parametrization, but typically used viscosity parametrization are generally ad hoc. In this study, we take a different approach, and investigate a time-dependent 1D $\alpha$-disk model in which the pressure is dominated by magnetic fields rather than photons. We compare the time evolution of thermally stable, strongly magnetized TDE disks to the simpler linear viscosity model. We find that the light curves of magnetized disks evolve as $L_{\rm UV}\propto t^{-5/6}$ for decades to centuries, and that this same evolution can be reproduced by the linear viscosity model for specific parameter choices. Additionally, we show that TDEs remain UV-bright for many years, suggesting we could possibly find fossil TDEs decades after their bursts. We estimate that ULTRASAT could detect hundreds of such events, providing an opportunity to study late-stage TDE physics and supermassive black hole (SMBH) properties. Finally, we explore the connection between TDE disks and quasi-periodic eruptions (QPEs) suggested by recent observations. One theoretical explanation involves TDE disks expanding to interact with extreme mass ratio inspirals (EMRIs), which produce X-ray flares as the EMRI passes through the disk. Our estimates indicate that magnetized TDE disks should exhibit QPEs earlier than those observed in AT2019qiz, suggesting that the QPEs may have begun before their first detection.
The recent observation of the ultra-high-energy neutrino event KM3-230213A by the KM3NeT experiment offers a compelling avenue to explore physics beyond the Standard Model. In this work, we explore a simplest possibility that this event originates from the decay of a super-heavy dark matter (SHDM). We consider a minimal extension of the type-I seesaw by including a singlet scalar and a singlet fermion DM, both being odd under a $Z_2$ symmetry. At high scale, the $Z_2$ symmetry is spontaneously broken by the vev of the $Z_2$ odd scalar, leading to i) mixing between DM and heavy right-handed neutrinos and ii) formation of domain walls (DW). In the former case, the decay of the DM to $\nu,h$ can give rise to the PeV neutrino event KM3-230213A. While, in the latter case, the disappearance of the DW can give rise to stochastic gravitational waves. We derive constraints on the DM lifetime as a function of DM mass ensuring consistency with IceCube, Auger upper limits and the observed KM3-230213A event. We found that KM3-230213A gives stringent constraint on DM mass ranging from $1.7\times10^8$ GeV to $5.5\times10^9$ GeV with lifetime $6.3\times10^{29}$ s to $3.6\times10^{29}$ s.
This paper introduces Firmamento, a new online platform designed for astronomical research, particularly for studying blazars and other multi-messenger emitters. Firmamento provides access to a wealth of astronomical data, including imaging, spectral, and timing information, along with tools for data analysis and machine learning. A key feature is the Error Region Counterpart Identifier (ERCI), which helps locate potential counterparts to gamma-ray and other high-energy sources pinpointing objects exhibiting blazar-like characteristics. Firmamento grants access to spectral data across the entire electromagnetic spectrum, offering well-populated SED generation and VHE detectability estimation. The platform also supports user-provided data, allowing researchers to incorporate their own source lists.
Numerical simulations of multidimensional astrophysical fluids present considerable challenges. However, the development of exascale computing has significantly enhanced computational capabilities, motivating the development of new codes that can take full advantage of these resources. In this article, we introduce HERACLES++, a new hydrodynamics code with high portability, optimized for exascale machines with different architectures and running efficiently both on CPUs and GPUs. The code is Eulerian and employs a Godunov finite-volume method to solve the hydrodynamics equations, which ensures accuracy in capturing shocks and discontinuities. It includes different Riemann solvers, equations of state, and gravity solvers. It works in Cartesian and spherical coordinates, either in 1-D, 2-D, or 3-D, and uses passive scalars to handle gases with several species. The code accepts a user-supplied heating or cooling term to treat a variety of astrophysical contexts. In addition to the usual series of benchmarking tests, we use HERACLES++ to simulate the propagation of a supernova shock in a red-supergiant star envelope, from minutes after core collapse until shock emergence. In 1-D, the results from HERACLES++ are in agreement with those of V1D for the same configuration. In 3-D, the Rayleigh-Taylor instability develops and modifies the 1-D picture by introducing density and composition fluctuations as well as turbulence. The focus on a wedge, rather than the full solid angle, and the ability of running HERACLES++ with a large number of GPUs allow for long-term simulations of 3-D supernova ejecta with a sub-degree resolution. Future development is to extend HERACLES++ to a radiation-hydrodynamics code.
In this work, we investigate the tidal deformability of regular black holes (RBHs). Employing different phenomenological models, we analyze their response to both test fields and gravitational perturbations, interpreting the latter within the framework of Einstein's field equations in the presence of an appropriate exotic matter distribution. Numerical and analytical methods reveal that RBHs exhibit non-trivial tidal responses, influenced by their regularization parameters and exotic matter distributions. The results obtained for test fields and gravitational perturbations are in qualitative agreement. This hints at the possibility that similar conclusions could hold if these spacetimes were interpreted as solutions of a modified gravitational action. Our findings suggest that RBHs possess distinct, though subtle, tidal signatures, which may serve as observational probes of their internal structure in gravitational wave detections.
The velocity-dependent Newtonian analogous potentials (NAPs) corresponding to general relativistic (GR) spacetimes accurately capture most of the relativistic features, including all classical tests of GR, effectively representing spacetime geometries in Newtonian terms. The NAP formulated by Tejeda \& Rosswog (TR13) for Schwarzschild spacetime has been applied to the standard thin accretion disk around a black hole (BH) as well as in the context of streamlines of noninteracting particles accreting onto a Schwarzschild BH, showing good agreement with the exact relativistic solutions. As a further application, here we explore the extent to which TR13 NAP could describe a transonic hydrodynamical spherical accretion flow in Schwarzschild spacetime within the framework of standard Newtonian hydrodynamics. Instead of obtaining a typical single "saddle-type" sonic transition, a "saddle-spiral pair" is produced, with the inner sonic point being an (unphysical) "spiral type" and the outer being a usual "saddle type." The Bondi accretion rate at outer sonic radii, however, remains consistent with that of the GR case. The primary reason for the deviation of our findings from the classical Bondi solution is likely due to the inconsistency between the Euler-type equation in the presence of velocity-dependent TR13 NAP within the standard Newtonian hydrodynamics framework, and the corresponding GR Euler equation, regardless of the fluid's energy. Our study suggests that a (modified) hydrodynamical formalism is needed to effectively implement such potentials in transonic accretion studies that align with the spirit of TR13-like NAP, while remaining consistent with the GR hydrodynamics. This could then essentially circumvent GR hydrodynamics or GR magnetohydrodynamics equations.