The disks of Active Galactic Nuclei (AGNs) are expected to be populated by numerous stars, either formed in the outer regions of the disk via gravitational instability, or captured from the nearby nuclear star cluster. Regardless of their formation mechanism, these stars experience altered evolutionary paths, mostly shaped by the accretion of dense disk material. In this study, through the comparison of different timescales, we chart the evolutionary outcomes of these AGN stars as a function of disk radius and across a range of supermassive black hole (SMBH) masses, spanning from $10^6$ to $10^9 \rm M_\odot$, for two popular AGN disk models. We find that, in the outer regions of the disk, stars evolve similarly to those in the interstellar medium, but in the inner and denser regions accretion quickly turns low-mass stars into massive stars, and their fate depends on just how quickly they accrete. If accretion occurs at a faster rate than nuclear burning, they can reach a quasi-steady `immortal' state. If stars accrete faster than they can thermally adjust, runaway accretion occurs, potentially preventing a quasi-steady state and altering the disk structure. During the AGN lifetime, in the regions of the disk that produce massive stars, supernovae (SNe) and Gamma-Ray Bursts (GRBs) may occur within the disk over a wide range of optical depths and ambient densities. Subsequently, in the final phase of the AGN, as the disk becomes depleted, formerly immortal stars will be unable to replenish their fuel, leading to additional SNe and GRBs.
We present a survey of 1D kinetic particle-in-cell simulations of quasi-parallel non-relativistic shocks to identify the environments favorable for electron acceleration. We explore an unprecedented range of shock speeds $v_{\rm sh}\approx 0.067-0.267\,c$, Alfvén Mach numbers $\mathcal{M}_{\rm A} = 5-40$, sonic Mach numbers $\mathcal{M}_{\rm s} = 5-160$, as well as the proton-to-electron mass ratios $m_{\rm i}/m_{\rm e}=16-1836$. We find that high Alfvén Mach number shocks can channel a large fraction of their kinetic energy into nonthermal particles, self-sustaining magnetic turbulence and acceleration to larger and larger energies. The fraction of injected particles is $\lesssim 0.5\%$ for electrons and $\approx 1\%$ for protons, and the corresponding energy efficiencies are $\lesssim 2\%$ and $\approx 10\%$, respectively. The extent of the nonthermal tail is sensitive to the Alfvén Mach number; when $\mathcal{M}_{\rm A}\lesssim 10$, the nonthermal electron distribution exhibits minimal growth beyond the average momentum of the downstream thermal protons, independently of the proton-to-electron mass ratio. Acceleration is slow for shocks with low sonic Mach numbers, yet nonthermal electrons still achieve momenta exceeding the downstream thermal proton momentum when the shock Alfvén Mach number is large enough. We provide simulation-based parametrizations of the transition from thermal to nonthermal distribution in the downstream (found at a momentum around $p_{\rm i,e}/m_{\rm i}v_{\rm sh} \approx 3\sqrt{m_{\rm i,e}/m_{\rm i}}$), as well as the ratio of nonthermal electron to proton number density. The results are applicable to many different environments and are important for modeling shock-powered nonthermal radiation.
Hot subdwarf B (sdB) stars are stripped helium-burning stars that are often found in close binaries, where they experience strong tidal interactions. The dissipation of tidally excited gravity waves alter their rotational evolution throughout the sdB lifetime. While many sdB binaries have well-measured rotational and orbital frequencies, there have been few theoretical efforts to accurately calculate the tidal torque produced by gravity waves. In this work, we directly calculate the tidal excitation of internal gravity waves in realistic sdB stellar models, and integrate the coupled spin-orbit evolution of sdB binaries. We find that for canonical sdB ($M_\mathrm{sdB}=0.47\,M_\odot$) binaries, the transitional orbital period below which they could reach tidal synchronization in the sdB lifetime is $\sim \! 0.2\;\mathrm{days}$, with weak dependence on the companion masses. For low-mass sdBs ($M_\mathrm{sdB}=0.37\,M_\odot$) formed from more massive progenitor stars, the transitional orbital period becomes $\sim \! 0.15\;\mathrm{days}$. These values are very similar to the tidal synchronization boundary ($\sim \! 0.2\;\mathrm{days}$) evident from observations. We discuss the dependence of tidal torques on stellar radii, and we make predictions for the rapidly rotating white dwarfs formed from synchronized sdB binaries.
The production of complex astronomical data is accelerating, especially with newer telescopes producing ever more large-scale surveys. The increased quantity, complexity, and variety of astronomical data demand a parallel increase in skill and sophistication in developing, deciding, and deploying statistical methods. Understanding limitations and appreciating nuances in statistical and machine learning methods and the reasoning behind them is essential for improving data-analytic proficiency and acumen. Aiming to facilitate such improvement in astronomy, we delineate cautionary tales in statistics via six maxims, with examples drawn from the astronomical literature. Inspired by the significant quality improvement in business and manufacturing processes by the routine adoption of Six Sigma, we hope the routine reflection on these Six Maxims will improve the quality of both data analysis and scientific findings in astronomy.
In order to improve the asteroseismic modelling efforts for red-giant stars, the numerical computation of theoretical individual oscillation modes for evolved red-giant models has to be made feasible. We aim to derive a method for circumventing the computational cost of computing oscillation spectra for models of red-giant stars with an average large frequency separation $\Delta\nu<15$ $\mu$Hz, thereby allowing for asteroseismic investigations of giants utilising individual frequencies. The proposed Truncated Scanning Method serves as a novel method detailing how the observable individual $l=0,1,2$ frequencies of red giants may be computed on realistic timescales through so-called model truncation. By carefully removing the innermost region of the stellar models, the g-mode influence on the oscillation spectra may be avoided, allowing estimation of the observable regions from the resulting pure p-mode oscillations. The appropriate observable frequency regions are subsequently scanned for the complete and un-truncated stellar model. The observable regions are determined by considering the limitations on observability from the internal mode coupling and damping, yielding consistent frequency spectra obtained at a much reduced computational cost. The Truncated Scanning Method proves the feasibility of obtaining the individual frequencies of red-giant models for a wide range of applications and research, demonstrating an improved computational efficiency by a factor of 10 or better. This means that the inclusion of $l=1,2$ individual frequencies is now a possibility in future asteroseismic modelling efforts of red-giant stars. Further potential avenues for improvements to this method are outlined for future pursuits.
this https URL . Submitted, comments welcome!