Search for: All records

Stellar streams that emerge from globular clusters (GCs) are thin stellar structures spread along the orbits of progenitor clusters. Numerical modeling of these streams is essential for understanding their interaction with the host galaxy's mass distribution. Traditional methods are either computationally expensive or oversimplified, motivating us to develop a fast and accurate approach using a particle spray algorithm. By conducting a series ofN-body simulations of GCs orbiting a host galaxy, we find that the position and velocity distributions of newly escaped stream particles are consistent across various GC masses and orbital parameters. Based on these distributions, we develop a new algorithm that avoids computing the detailed internal cluster dynamics by directly drawing tracer particles from these distributions. This algorithm correctly reproduces the action space distribution of stream particles and achieves a 10% accuracy in stream morphology and velocities compared toN-body simulations. To facilitate broader use, we have implemented this algorithm in galactic dynamics codesagama,gala,galax, andgalpy.

Dark matter (DM) bars that shadow stellar bars have been previously shown to form in idealized simulations of isolated disk galaxies, but have yet to be studied in a fully cosmological context. In this work, we analyze a population of disk galaxies within the TNG50 simulation to determine the characteristics of their dark bars. We estimate bar strength and orientation using both the in-plane FourierA₂density moments and the quadrupolar coefficients of the spherical harmonic basis function expansions of the density. We additionally present two novel methods for measuring the bar pattern speed Ω_pand rotation axis orientation using these coefficients, and apply them to one sample galaxy located in a TNG50 subbox. Consistent with isolated simulations, DM bars are shorter than their stellar counterparts and are 75% weaker inA₂. DM bars dominate the shape of the inner halo potential and are apparent in the time series of quadrupolar coefficients. In our selected subbox galaxy, the stellar and dark bars remain co-aligned throughout the last 8 Gyr and have identical Ω_p. Pattern speed Ω_pevolves considerably over the last 8 Gyr, consistent with torques on the bars due to dynamical friction and gas accretion, and is seen to increase following a merger att_lb= 1.5 Gyr. Rather than remaining static in time, the bar rotation axis displays both precession and nutation possibly caused by torques outside the plane of rotation. We find that the shape of the stellar and DM mass distributions are tightly correlated with Ω_p.

Barred galaxies exhibit boxy/peanut or X-shapes (BP/X) protruding from their disks in edge-on views. Two types of BP/X morphologies exist depending on whether the X-wings meet at the center (CX) or are off-centered (OX). Orbital studies indicate that various orbital types can generate X-shaped structures. Here we provide a classification approach that identifies the specific orbit families responsible for generating OX- and CX-shaped structures. Applying this approach to three differentN-body bar models, we show that both OX and CX structures are associated with thex1 orbit family, but OX-supporting orbits possess higher angular momentum (closer tox1 orbits) than orbits in CX structures. Consequently, as the bar slows down, the contribution of higher angular momentum OX-supporting orbits decreases and that of lower angular momentum orbits increases, resulting in an evolution of the morphology from OX to CX. If the bar does not slow down, the shape of the BP/X structure and the fractions of OX/CX-supporting orbits remain substantially unchanged. Bars that do not undergo buckling but that do slow down initially show the OX structure and are dominated by high angular momentum orbits, transitioning to a CX morphology. Bars that buckle exhibit a combination of both OX- and CX-supporting orbits immediately after the buckling but become more CX dominated as their pattern speed decreases. This study demonstrates that the evolution of BP/X morphology and orbit populations strongly depends on the evolution of the bar angular momentum.

We apply the barred Schwarzschild method developed by Tahmasebzadeh et al. (2022) to a barred S0 galaxy, NGC 4371, observed by IFU instruments from the TIMER and ATLAS3D projects. We construct the gravitational potential by combining a fixed black hole mass, a spherical dark matter halo, and stellar mass distribution deprojected from 3.6 μm S$^4$G image considering an axisymmetric disc and a triaxial bar. We independently modelled kinematic data from TIMER and ATLAS3D. Both models fit the data remarkably well. We find a consistent bar pattern speed from the two sets of models with $\Omega _{\rm p} = 23.6 \pm 2.8 \, \mathrm{km \, s^{-1} \, kpc^{-1} }$ and $\Omega _{\rm p} = 22.4 \pm 3.5 \, \mathrm{km \, s^{-1} \, kpc^{-1} }$, respectively. The dimensionless bar rotation parameter is determined to be $\mathcal {R} \equiv R_{\rm cor}/R_{\rm bar}=1.88 \pm 0.37$, indicating a likely slow bar in NGC 4371. Additionally, our model predicts a high amount of dark matter within the bar region ($M_{\rm DM}/ M_{\rm total}$$\sim 0.51 \pm 0.06$), which, aligned with the predictions of cosmological simulations, indicates that fast bars are generally found in baryon-dominated discs. Based on the best-fitting model, we further decompose the galaxy into multiple 3D orbital structures, including a BP/X bar, a classical bulge, a nuclear disc, and a main disc. The BP/X bar is not perfectly included in the input 3D density model, but BP/X-supporting orbits are picked through the fitting to the kinematic data. This is the first time a real barred galaxy has been modelled utilizing the Schwarzschild method including a 3D bar.

We present a new method to infer the 3D luminosity distributions of edge-on barred galaxies with boxy-peanut/X (BP/X) shaped structures from their 2D surface brightness distributions. Our method relies on forward modelling of newly introduced parametric 3D density distributions for the BP/X bar, disc and other components using an existing image fitting software package (imfit). We validate our method using an N-body simulation of a barred disc galaxy with a moderately strong BP/X shape. For fixed orientation angles, the derived 3D BP/X-shaped density distribution is shown to yield a gravitational potential that is accurate to at least 5 per cent and forces that are accurate to at least 15 per cent, with average errors being $\sim 1.5~{{\ \rm per\ cent}}$ for both. When additional quantities of interest, such as the orientation of the bar to the line of sight, its pattern speed, and the stellar mass-to-light ratio are unknown they can be recovered to high accuracy by providing the parametric density distribution to the Schwarzschild modelling code FORSTAND. We also explore the ability of our models to recover the mass of the central supermassive black hole. This method is the first to be able to accurately recover both the orientation of the bar to the line of sight and its pattern speed when the disc is perfectly edge-on.

UsingN-body simulations, we explore the effects of growing a supermassive black hole (SMBH) prior to or during the formation of a stellar bar. Keeping the final mass and growth rate of the SMBH fixed, we show that if it is introduced before or while the bar is still growing, the SMBH does not cause a decrease in bar amplitude. Rather, in most cases, it is strengthened. In addition, an early-growing SMBH always either decreases the buckling amplitude, delays buckling, or both. This weakening of buckling is caused by an increase in the disk vertical velocity dispersion at radii well beyond the nominal black hole sphere of influence. While we find considerable stochasticity and sensitivity to initial conditions, the only case where the SMBH causes a decrease in bar amplitude is when it is introduced after the bar has attained a steady state. In this case, we confirm previous findings that the decrease in bar strength is a result of scattering of bar-supporting orbits with small pericenter radii. By heating the inner disk both radially and vertically, an early-growing SMBH increases the fraction of stars that can be captured by the inner Lindblad resonance (ILR) and the vertical ILR, thereby strengthening both the bar and the boxy-peanut-shaped bulge. Using orbital frequency analysis of star particles, we show that when an SMBH is introduced early and the bar forms around it, the bar is populated by different families of regular bar-supporting orbits than when the bar forms without an SMBH.

We describe the results of a new reverberation mapping program focused on the nearby Seyfert galaxy NGC 3227. Photometric and spectroscopic monitoring was carried out from 2022 December to 2023 June with the Las Cumbres Observatory network of telescopes. We detected time delays in several optical broad emission lines, with Hβhaving the longest delay at${\tau }_{\mathrm{cent}}={4.0}_{-0.9}^{+0.9}$days and Heiihaving the shortest delay with${\tau }_{\mathrm{cent}}={0.9}_{-0.8}^{+1.1}$days. We also detect velocity-resolved behavior of the Hβemission line, with different line-of-sight velocities corresponding to different observed time delays. Combining the integrated Hβtime delay with the width of the variable component of the emission line and a standard scale factor suggests a black hole mass of$M_{\mathrm{BH}}={1.1}_{-0.3}^{+0.2}\times {10}^7$M_⊙. Modeling of the full velocity-resolved response of the Hβemission line with the phenomenological codeCARAMELfinds a similar mass of$M_{\mathrm{BH}}={1.2}_{-0.7}^{+1.5}\times {10}^7$M_⊙and suggests that the Hβ-emitting broad-line region (BLR) may be represented by a biconical or flared disk structure that we are viewing at an inclination angle ofθ_i≈ 33° and with gas motions that are dominated by rotation. The new photoionization-based BLR modeling toolBELMACfinds general agreement with the observations when assuming the best-fitCARAMELresults; however,BELMACprefers a thick-disk geometry and kinematics that are equally composed of rotation and inflow. Both codes infer a radially extended and flattened BLR that is not outflowing.

Many barred galaxies exhibit upturns (shoulders) in their bar-major-axis density profile. Simulation studies have suggested that shoulders are supported by loopedx₁orbits, occur in growing bars, and can appear after bar buckling. We investigate the orbital support and evolution of shoulders via frequency analyses of orbits in simulations. We confirm that looped orbits are shoulder-supporting, and can remain so, to a lesser extent, after being vertically thickened. We show that looped orbits appear at the resonance ( Ω_φ− Ω_P)/Ω_R= 1/2 (analogous to the classical inner Lindblad resonance, and here called ILR) with vertical-to-radial frequency ratios 1 ≲ Ω_z/Ω_R≲ 3/2 (verticallywarmorbits).Coolorbits at the ILR (those with Ω_z/Ω_R> 3/2) are vertically thin and have no loops, contributing negligibly to shoulders. As bars slow and thicken, either secularly or by buckling, they populate warm orbits at the ILR. Further thickening carries these orbits toward crossing the vertical ILR [vILR, ( Ω_φ− Ω_P)/Ω_z= 1/2], where they convert in-plane motion to vertical motion, become chaotic, kinematically hotter, and less shoulder-supporting. Hence, persistent shoulders require bars to trap new stars, consistent with the need for a growing bar. Since buckling speeds up trapping on warm orbits at the ILR, it can be followed by shoulder formation, as seen in simulations. This sequence supports the recent observational finding that shoulders likely precede the emergence of BP-bulges. The python module for the frequency analysis,naif, is made available.

We present a new constraint on the mass of the black hole in the active S0 galaxy NGC 5273. Due to the proximity of the galaxy at 16.6 ± 2.1 Mpc, we were able to resolve and extract the bulk motions of stars near the central black hole using adaptive-optics-assisted observations with the Gemini Near-infrared Integral Field Spectrograph, as well as constrain the large-scale kinematics using archival Spectroscopic Areal Unit for Research and Optical Nebulae spectroscopy. High-resolution Hubble Space Telescope imaging allowed us to generate a surface-brightness decomposition, determine approximate mass-to-light ratios for the bulge and disk, and obtain an estimate for the disk inclination. We constructed an extensive library of dynamical models using the Schwarzschild orbit-superposition code FORSTAND, exploring a range of disk and bulge shapes, halo masses, etc. We determined a black hole mass ofM_•= [0.5–2] × 10⁷M_⊙, where the low side of the range is in agreement with the reverberation mapping measurement ofM_•= [4.7 ± 1.6] × 10⁶M_⊙. NGC 5273 is one of the few nearby galaxies that hosts a broad-lined active galactic nucleus, allowing a crucial comparison of black hole masses derived from independent mass-measurement techniques.

We present the results of a new reverberation mapping campaign for the broad-line active galactic nucleus (AGN) in the edge-on spiral IC 4329A. Monitoring of the optical continuum withV-band photometry and broad emission-line flux variability with moderate-resolution spectroscopy allowed emission-line light curves to be measured for Hβ, Hγ, and Heiiλ4686. We find a time delay of${16.3}_{-2.3}^{+2.6}$days for Hβ, a similar time delay of${16.0}_{-2.6}^{+4.8}$days for Hγ, and an unresolved time delay of$-{0.6}_{-3.9}^{+3.9}$days for Heii. The time delay for Hβis consistent with the predicted value from the relationship between AGN luminosity and broad-line region radius, after correction for the ∼2.4 mag of intrinsic extinction at 5100 Å. Combining the measured time delay for Hβwith the broad emission-line width and an adopted value of 〈f〉 = 4.8, we find a central supermassive black hole mass of$M_{\mathrm{BH}}={6.8}_{-1.1}^{+1.2}\times {10}^7$M_⊙. Velocity-resolved time delays were measured across the broad Hβemission-line profile and may be consistent with an “M”-like shape. Modeling of the full reverberation response of Hβwas able to provide only modest constraints on some parameters, but does exhibit agreement with the black hole mass and average time delay. The models also suggest that the AGN structure is misaligned by a large amount from the edge-on galaxy disk. This is consistent with expectations from the unified model of AGNs, in which broad emission lines are expected to be visible only for AGNs that are viewed at relatively face-on inclinations.