Search

Ongoing deep integral-field spectroscopy (IFS) observations of disc galaxies provide opportunities for comparison with the Milky Way (MW) to understand galaxy evolution. However, such comparisons are marred by many challenges such as selection effects, differences in observations and methodology, and proper validation of full-spectrum fitting methods. In this study, we present a novel code GalCraft to address these challenges by generating mock IFS data cubes of the MW using simple stellar population models and a mock MW stellar catalogue derived from E-galaxia. We use the widely adopted full-spectrum fitting code ppxf to investigate the ability to recover kinematics and stellar populations for an edge-on mock MW IFS observation. We confirm that differences in kinematics, mean age, |$[\mathrm{M/H}]$|⁠ , and |$[\mathrm{\alpha /Fe}]$| between thin and thick discs can be distinguished. However, the age distribution is overestimated in the ranges between |$2{-}4$| and |$12{-}14$|  Gyr compared to the expected values. This is likely due to the age spacing and degeneracy of SSP templates. We find systematic offsets in the recovered kinematics due to insufficient spectral resolution and the variation of line-of-sight velocity distribution with age and |$[\mathrm{M/H}]$|⁠. With future higher resolution and multi- |$[\mathrm{\alpha /Fe}]$| simple stellar population templates, GalCraft will be useful to validate key signatures such as |$[\mathrm{\alpha /Fe}]$| – |$[\mathrm{M/H}]$| distribution at different R and |$|z|$| and potentially infer radial migration and kinematic heating efficiency to study detailed chemodynamical evolution of MW-like galaxies. [ABSTRACT FROM AUTHOR]

We use high-resolution magnetohydrodynamic simulations of isolated disc galaxies to investigate the co-evolution of magnetic fields with a self-regulated, star-forming interstellar medium (ISM). The simulations are conducted using the ramses adaptive mesh refinement code on the standard agora initial condition, with gas cooling, star formation, and feedback. We run galaxies with a variety of initial magnetic field strengths. The fields evolve and achieve approximate saturation within 500 Myr, but at different levels. The galaxies reach a quasi-steady state, with slowly declining star formation due to both gas consumption and increase in the field strength at intermediate ISM densities. We connect this behaviour to differences in the gas properties and overall structure of the galaxies. Stronger magnetic fields limit supernova bubble sizes. Different cases support the ISM using varying combinations of magnetic pressure, turbulence, and thermal energy. Initially, |$\gtrsim\!\! 1\ \mu \mathrm{ G}$| magnetic fields evolve modestly and dominate support at all radii. Conversely, initially weaker fields grow through feedback and turbulence but never dominate the support. This is reflected in the stability of the gas disc. This interplay determines the overall distribution of star formation in each case. We conclude that an initially weak field can grow to produce a realistic model of a local disc galaxy, but starting with typically assumed field strengths (⁠|$\gtrsim\!\! 1\ \mu \mathrm{ G}$|⁠) will not. [ABSTRACT FROM AUTHOR]

We use high-resolution (⁠|$\simeq$| 35pc) hydrodynamical simulations of galaxy formation to investigate the relation between gas accretion and star formation in galaxies hosted by dark matter haloes of mass |$10^{12}$| |${{\mathrm{M}}_{\odot }}$| at |$z = 2$|⁠. At high-redshift, cold-accreted gas is expected to be readily available for star formation, while gas accreted in a hot mode is expected to require a longer time to cool down before being able to form stars. Contrary to these expectations, we find that the majority of cold-accreted gas takes several hundred Myr longer to form stars than hot-accreted gas after it reaches the inner circumgalactic medium (CGM). Approximately 10 per cent of the cold-accreted gas flows rapidly through the inner CGM on to the galactic disc. The remaining 90 per cent is trapped in a turbulent accretion region that extends up to |$\sim 50$|  per cent of the virial radius, from which it takes several hundred Myr for the gas to be transported to the star-forming disc. In contrast, most hot shock-heated gas avoids this 'slow track', and accretes directly from the CGM on to the disc where stars can form. We find that shock-heating of cold gas after accretion in the inner CGM and supernova-driven outflows contribute to, but do not fully explain, the delay in star formation. These processes combined slow down the delivery of cold-accreted gas to the galactic disc and consequently limit the rate of star formation in Milky Way mass galaxies at |$z \gt 2$|⁠. [ABSTRACT FROM AUTHOR]

This study outlines a detailed investigation using CCD UBV and Gaia DR3 data sets of the two open clusters Ruprecht 1 (Rup‐1) and Ruprecht 171 (Rup‐171). Fundamental astrophysical parameters such as color excesses, photometric metallicities, ages, and isochrone distances were based on UBV‐data analyses, whereas membership probability calculations, structural and astrophysical parameters, as well as the kinematic analyses were based on Gaia DR3‐data. We identified 74 and 596 stars as the most probable cluster members with membership probabilities over 50% for Rup‐1 and Rup‐171, respectively. The color excesses E(B−V)$$ E\left(B-V\right) $$ were obtained as 0.166±0.022$$ 0.166\pm 0.022 $$ and 0.301±0.027$$ 0.301\pm 0.027 $$ mag for Rup‐1 and Rup‐171, respectively. Photometric metallicity analyses were performed by considering F‐G type main‐sequence member stars and found to be [Fe/H]=−0.09±0.16$$ \left[\mathrm{Fe}/\mathrm{H}\right]=-0.09\pm 0.16 $$ and [Fe/H]=−0.20±0.20$$ \left[\mathrm{Fe}/\mathrm{H}\right]=-0.20\pm 0.20 $$ dex for Rup‐1 and Rup‐171, respectively. Ages and distances were based on both UBV and Gaia‐data analyses; according to isochrone‐fitting these values were estimated to be t=580±60$$ t=580\pm 60 $$ Myr, d=1469±57$$ d=1469\pm 57 $$ pc for Rup‐1 and t=2700±200$$ t=2700\pm 200 $$ Myr, d=1509±69$$ d=1509\pm 69 $$ pc for Rup‐171. The present‐day mass function slope of Rup‐1 was estimated as 1.26±0.32$$ 1.26\pm 0.32 $$ and Rup‐171 as 1.53±1.49$$ 1.53\pm 1.49 $$. Galactic orbit integration analyses showed that both of the clusters might be formed outside the solar circle. [ABSTRACT FROM AUTHOR]

Early-type dwarf galaxies constitute a prevalent population in the central regions of rich groups and clusters in the local Universe. These low-luminosity and low-mass stellar systems play a fundamental role in the assembly of the luminous galaxies observed today, according to the Lambda cold dark matter hierarchical theory. The origin of early-type dwarfs has been linked to the transformation of disc galaxies interacting with the intracluster medium, especially in dense environments. However, the existence of low-luminosity early-type galaxies in low-density environments presents a challenge to this scenario. This study presents a comprehensive photometric and spectroscopic analysis of the early-type dwarf galaxy CGCG014−074 using deep Gemini GMOS (Gemini Multi-Object Spectrograph) data, focusing on its peculiarities and evolutionary implications. CGCG014−074 exhibits distinct features, including a rotating inner disc, an extended stellar formation with a quiescent phase since about 2 Gyr ago, and the presence of boxy isophotes. From the kinematic analysis, we confirm CGCG014−074 as a nucleated early-type dwarf galaxy with embedded disc. The study of its stellar population parameters using different methods provides significant insights into the galaxy's evolutionary history. These results show an old and metal-poor nucleus (⁠|${\sim}9.3$| Gyr and |$\mathrm{[Z/H]}\sim -0.84$| dex), while the stellar disc is younger (⁠|${\sim}4.4$| Gyr) with a higher metallicity (⁠|$\mathrm{[Z/H]}\sim -0.40$| dex). These distinctive features collectively position CGCG014−074 as a likely building block galaxy that has evolved passively throughout its history. [ABSTRACT FROM AUTHOR]

Milky Way analogues (MWAs) provide an alternative insight into the various pathways that lead to the formation of disc galaxies with similar properties to the Milky Way (MW). In this study, we explore different selection techniques for identifying MWAs in the Sydney-AAO (Australian Astronomical Observatory) Multi-object Integral field spectrograph (SAMI) Galaxy Survey. We utilize a nearest neighbours method to define MWAs using four selection parameters including stellar mass (⁠|$M_{\star }$|⁠), star formation rate (⁠|$\rm SFR$|⁠), bulge-to-total ratio (⁠|$B/T$|⁠), and disc effective radius (⁠|$R_{\rm {e}}$|⁠). Based on 15 different selection combinations, we find that including |$M_{\star }$| and SFR is essential for minimizing biases in the average MWA properties as compared to the MW. Furthermore, given the MW's smaller than average size, selection combinations without |$R_{\rm {e}}$| result in MWAs being too large. Lastly, we find that |$B/T$| is the least important parameter out of the four tested parameters. Using all four selection criteria, we define the top 10 most MW-like galaxies in the GAMA and Cluster regions of the SAMI survey. These most MW-like galaxies are typically barred spirals, with kinematically cold rotating discs and reside in a wide range of environments. Surprisingly, we find no significant differences between the MWAs selected from the GAMA and Cluster regions. Our work highlights the importance of using multiple selection criteria for finding MWAs and also demonstrates potential biases in previous MWA studies. [ABSTRACT FROM AUTHOR]

Studying the nature of spiral arms is essential for understanding the formation of the intricate disc structure of the Milky Way. The European Space Agency's Gaia mission has provided revolutionary observational data that have uncovered detailed kinematical features of stars in the Milky Way. However, so far the nature of spiral arms continues to remain a mystery. Here, we present that the stellar kinematics traced by the classical Cepheids around the Perseus and Outer spiral arms in the Milky Way show strikingly different kinematical properties from each other: the radial and azimuthal velocities of Cepheids show positive and negative correlations in the Perseus and Outer arms, respectively. We also found that the dynamic spiral arms commonly seen in an N -body/hydrodynamic simulation of a Milky Way-like galaxy can naturally explain the observed kinematic trends. Furthermore, a comparison with such a simulation suggests that the Perseus arm is being disrupted, while the Outer arm is growing. Our findings suggest that two neighbouring spiral arms in distinct evolutionary phases – growing and disrupting phases – coexist in the Milky Way. [ABSTRACT FROM AUTHOR]

Warps are common vertical asymmetries that appear in the outer parts of the galactic discs, bending one part upwards and the other downwards. Many mechanisms can trigger warp formation, including tidal interactions. The interactions with satellites distort the edges of the disc and can also change the central morphology, impacting, for example, the development of a galactic bar. In mergers events, the bar can be weakened or even destroyed. In this study, we aim to compare barred and non-barred galaxy models and their susceptibility to warping. To analyse the effects of induced warps, we used N -body simulations of a barred and a non-barred central galaxy interacting with satellites of varying masses (⁠|$0.1 \times 10^{10}\, \mathrm{{\rm M}_{\odot }}$|⁠ , |$0.5 \times 10^{10}\, \mathrm{{\rm M}_{\odot }}$|⁠ , and |$1 \times 10^{10}\, \mathrm{{\rm M}_{\odot }}$|⁠) and initial orbital radii (10, 20, and 30 kpc). We also ran isolated simulations of the central galaxies for comparison. We found that the induced warps are stronger in the barred galaxy compared with the non-barred galaxy, in perturbed and isolated models. In addition, the masses of the satellites determine the level of destruction of the bar and the intensity of the induced warp. The time in which the bar will be weakened or destroyed depends on the orbital radius of the satellite. [ABSTRACT FROM AUTHOR]

We present quantizer , a new collisionless N -body code tailored for cosmological and galaxy-scale simulations. The code uses the particle-mesh method, combined with a novel adaptive mesh refinement (AMR) technique. The design of the AMR method of quantizer maximizes the efficient use of refinement zones by combining a cell-by-cell scheme with parallelization-friendly regular patches. We test the code using three astrophysical systems: a binary system of two particles without AMR to test orbital accuracy; an isolated elliptical galaxy modelled as a Plummer sphere; and a merger of two disc galaxies. Comparisons are made with the well-known particle-mesh AMR code ramses , and with our own direct particle-particle solver. The results show that the accuracy of quantizer is equal to or better than ramses , and the runtime is approximately 2.5 to 10 times faster than that code, depending on the problem. [ABSTRACT FROM AUTHOR]

We investigate how the stellar rotational support changes as a function of spatially resolved stellar population age (⁠|$D_n4000$|⁠) and relative central stellar surface density (⁠|$\Delta \Sigma _1$|⁠) for MaNGA isolated/central disc galaxies. We find that the galaxy rotational support indicator |$\lambda _{R_\mathrm{e}}$| varies smoothly as a function of |$\Delta \Sigma _1$| and |$D_n4000$|⁠. |$D_n4000$| versus |$\Delta \Sigma _1$| follows a 'J-shape', with |$\lambda _{R_\mathrm{e}}$| contributing to the scatters. In this 'J-shaped' pattern rotational support increases with central |$D_n4000$| when |$\Delta \Sigma _1$| is low but decreases with |$\Delta \Sigma _1$| when |$\Delta \Sigma _1$| is high. Restricting attention to low- |$\Delta \Sigma _1$| (i.e. large-radius) galaxies, we suggest that the trend of increasing rotational support with |$D_n4000$| for these objects is produced by a mix of two different processes, a primary trend characterized by growth in |$\lambda _{R_\mathrm{e}}$| along with mass through gas accretion, on top of which disturbance episodes are overlaid, which reduce rotational support and trigger increased star formation. An additional finding is that star-forming galaxies with low |$\Delta \Sigma _1$| have relatively larger radii than galaxies with higher |$\Delta \Sigma _1$| at fixed stellar mass. Assuming that these relative radii rankings are preserved while galaxies are star forming then implies clear evolutionary paths in central |$D_n4000$| versus |$\Delta \Sigma _1$|⁠. The paper closes with comments on the implications that these paths have for the evolution of pseudo-bulges versus classical bulges. The utility of using |$\rm D_n4000$| – |$\Delta \Sigma _1$| to study |$\lambda _{R_\mathrm{e}}$| reinforces the notion that galaxy kinematics correlate both with structure and with stellar-population state, and indicates the importance of a multidimensional description for understanding bulge and galaxy evolution. [ABSTRACT FROM AUTHOR]

The baryonic mass–velocity relation provides an important test of different galaxy dynamics models such as Lambda–cold dark matter (Λ CDM) and alternatives like Modified Newtonian Dynamics (MOND). Novel nonlinear density wave theory with a soliton solution gives an opportunity to test whether the derived rotational velocity expression is able to support the well known Tully–Fisher empirical relation between mass and rotation velocity in disk galaxies. Initial assumptions do not involve any larger dark matter halo that supports the stability of the very thin galactic disk nor any modified gravitational acceleration acting on galactic scales. It rather follows an important gravitational interaction between constituents of disk mass in the outer part of the disk via gravitational potential. Data are obtained by a fitting procedure applied on the sample of 81 rotational curves of late type spirals using expressions for the rotational velocity derived as an exact, a self-consistent solution of the nonlinear Schrodinger (NLS) equation for galactic surface mass density. The location of these selected objects in the baryonic mass–rotation velocity plane follows the relation log M b = 3.7 ± 0.2 log V f l a t + 2.7 ± 0.4 in marginal agreement with the findings in the literature. [ABSTRACT FROM AUTHOR]

The North–South asymmetry in the number density and bulk velocity of stars in the Solar neighbourhood provides valuable insights into the formation and evolution of the Milky Way disc. Our objective is to investigate the wave-like disc oscillations of mono-age stellar populations in the Solar neighbourhood using data from Gaia Data Release 3. We have selected a comprehensive sample of main-sequence turn-off stars. The ages of these stars can be accurately determined using isochrone fitting methods. Our findings indicate that the North–South density and mean vertical velocity asymmetries remain consistent across all age groups. The uniformity of perturbations across all subsamples suggests that all populations are responding to the same external influence, which likely affects them irrespective of their age. Moreover, the fact that these perturbations appear consistently implies they could be either ongoing or recent. Regarding vertical velocity dispersions, we observe that older stars exhibit larger dispersions. [ABSTRACT FROM AUTHOR]

The question of whether the Milky Way's disc extends to low metallicity has been the subject of debate for many years. We aim to address the question by employing a large sample of giant stars with radial velocities and homogeneous metallicities based on the Gaia Data Release 3 BP/RP(XP) spectra. We study the 3D velocity distribution of stars in various metallicity ranges, including the very metal-poor (VMP) regime (⁠|$\mathrm{[M/H]} < -2.0$|⁠). We find that a clear, stand-alone disc population, i.e. that with a ratio of rotational velocity to velocity dispersion |$v/\sigma > 1$|⁠ , starts to emerge only around |$\mathrm{[M/H]} \sim -1.3$|⁠ , and is not visible for |$\mathrm{[M/H]} < -1.6$|⁠. Using Gaussian mixture modelling, we show that there are two halo populations in the VMP regime: one stationary and one with a net prograde rotation of |$\sim\!\! 80\ \mathrm{km\, s}^{-1}$|⁠. In this low-metallicity range, we are able to place constraints on the contribution of a rotation-supported thick disc sub-population to a maximum of |$\sim\!\! 3$| per cent in our sample. We compare our results to previous claims of discy VMP stars in both observations and simulations and find that having a prograde halo component could explain most of these. [ABSTRACT FROM AUTHOR]

We study solutions of the collisionless Boltzmann equation (CBE) in a functional Koopman representation. This facilitates the use of linear spectral techniques characteristic of the analysis of Schrödinger-type equations. For illustrative purposes, we consider the classical phase mixing of a non-interacting distribution function in a quartic potential. Solutions are determined perturbatively relative to a harmonic oscillator. We impose a form of coarse-graining by choosing a finite-dimensional basis to represent the distribution function and time evolution operators, which sets a minimum length-scale on phase space structure. We observe a relationship between the dimension of the representation and the multiplicity of the harmonic oscillator eigenvalues. System dynamics are understood in terms of degenerate subspaces of the linear operator spectra. Each subspace is associated with a mode of the harmonic oscillator, the first two being bending and breathing structures. The quartic potential splits the degenerate eigenvalues within each subspace. This facilitates the formation of spiral structure as deformations from the harmonic oscillator modes. We ultimately argue that this construction provides a promising avenue for study of self-interacting systems experiencing phase mixing, which is an outstanding problem in the context of the Gaia DR2 vertical phase space spirals. [ABSTRACT FROM AUTHOR]

Making use of the APOGEE DR17 catalogue with high quality data for 143 509 red giant branch stars we explore the strength of different mechanisms that causes a star to radially migrate in the Milky Way stellar disc. At any position in the disc we find stars that are more metal-rich than the local interstellar medium. This is surprising and normally attributed to the migration of these stars after their formation inside their current Galactocentric radius. Such stars are prime candidates for studying the strength of different migratory processes. We specifically select two types of metal-rich stars: (i) super metal-rich stars (⁠|$\mathrm{[Fe/H]}\gt 0.2$|⁠) and (ii) stars that are more metal-rich than their local environment. For both, we explore the distribution of orbital parameters and ages as evidence of their migration history. We find that most super metal-rich stars have experienced some amount of churning as they have orbits with |$R_g\gtrsim 5 \,{\rm kpc}$|⁠. Furthermore, about half of the super metal-rich stars are on non-circular orbits (⁠|$\mathrm{ecc} \gt 0.15$|⁠) and therefore also have experienced blurring. The metallicity of young stars in our sample is generally the same as the metallicity of the interstellar medium, suggesting they have not radially migrated yet. Stars with lower metallicity than the local environment have intermediate to old ages. We further find that super metal-rich stars have approximately the same age distribution at all Galactocentric radii, which suggests that radial migration is a key mechanism responsible for the chemical compositions of stellar populations in the Milky Way. [ABSTRACT FROM AUTHOR]

The study presents a theoretical framework for understanding the role of dark matter on the stability of the galactic disc. We model the galaxy as a two-component system consisting of stars and gas in equilibrium with an external dark matter halo. We derive the equations governing the growth of perturbations and obtain a stability criterion that connects the potential of the dark matter halo and the gas fraction with the stability levels of the galaxy. We find that a two-component disc is more susceptible to the growth of gravitational instabilities than individual components, particularly as gas fractions increase. However, the external field, due to the dark matter halo, acts as a stabilizing agent and increases the net stability levels even in the presence of a cold gas component. We apply the stability criterion to models of the Milky Way, low surface brightness galaxies, and baryon-dominated cold rotating disc galaxies observed in the early universe. Our results show that the potential due to the dark matter halo plays a significant role in stabilizing nearby galaxies, such as the Milky Way, and low surface brightness galaxies, which would otherwise be prone to local gravitational instabilities. However, we find that the baryon-dominated cold disc galaxies observed in the early universe remain susceptible to the growth of local gravitational instabilities despite the stabilizing effect of the dark matter halo. [ABSTRACT FROM AUTHOR]

Gaseous, disc–halo interfaces are shaped by processes that are critical to galaxy evolution, including gas accretion and outflows. Extraplanar diffuse ionized gas (eDIG) layers are characterized by scale heights that largely exceed those predicted by their temperature, suggesting the presence of turbulent energy injection from star-formation feedback. However, the origin of this large-scale height remains uncertain. To explore the connection between eDIG and star-forming discs, we present a spatially resolved case study of a nearby pair of sub- |$L_*$|⁠ , intermediately inclined disc galaxies NGC 3511/3513. We decompose optical nebular lines observed using long-slit spectroscopy into narrow and broad velocity components. In NGC 3511, the broad component has three distinctive characteristics in comparison to the narrow component: (1) significantly higher velocity dispersions (a median |$\langle \sigma \rangle _{\text{Broad}} = 24$|  km s |$^{-1}$| compared to |$\langle \sigma \rangle _{\text{Narrow}} = 13$|  km s |$^{-1}$|⁠), (2) elevated [N  ii ] |$\lambda$| 6583/H |$\alpha$| and [S  ii ] |$\lambda$| 6716/H |$\alpha$| line ratios, and (3) a rotational velocity lag. These characteristics support the origin of the broad component in an extraplanar, gaseous disc. In NGC 3513, the broad component reveals disc–halo circulation via localized outflows at radius |$\lesssim 1$|  kpc. For NGC 3511, we test a vertical hydrostatic equilibrium model with pressure support supplied by thermal and turbulent motions. Under this assumption, the eDIG velocity dispersion corresponds to a scale height |$h_{z} \gtrsim 0.2 - 0.4$|  kpc at |$R = 3 - 5$|  kpc, a factor of a few above the thermal scale height (⁠|$h_{z} \lesssim 0.1$|  kpc). This highlights the importance of turbulent motions to the vertical structure of the gaseous, disc–halo interface. [ABSTRACT FROM AUTHOR]

We use smoothed particle hydrodynamics simulations of isolated Milky Way-mass disc galaxies that include cold, interstellar gas to test subgrid prescriptions for star formation (SF). Our fiducial model combines a Schmidt law with a gravitational instability criterion, but we also test density thresholds and temperature ceilings. While SF histories are insensitive to the prescription for SF, the Kennicutt–Schmidt (KS) relations between SF rate and gas surface density can discriminate between models. We show that our fiducial model, with an SF efficiency per free-fall time of 1 per cent, agrees with spatially resolved and azimuthally averaged observed KS relations for neutral, atomic, and molecular gas. Density thresholds do not perform as well. While temperature ceilings selecting cold, molecular gas can match the data for galaxies with solar metallicity, they are unsuitable for very low-metallicity gas and hence for cosmological simulations. We argue that SF criteria should be applied at the resolution limit rather than at a fixed physical scale, which means that we should aim for numerical convergence of observables rather than of the properties of gas labelled as star-forming. Our fiducial model yields good convergence when the mass resolution is varied by nearly 4 orders of magnitude, with the exception of the spatially resolved molecular KS relation at low surface densities. For the gravitational instability criterion, we quantify the impact on the KS relations of gravitational softening, the SF efficiency, and the strength of supernova feedback, as well as of observable parameters such as the inclusion of ionized gas, the averaging scale, and the metallicity. [ABSTRACT FROM AUTHOR]

The origin of diverse kinematic morphologies observed in low-mass galaxies is unclear. In this study, we investigate the kinematic morphologies of central galaxies with stellar mass |$10^{8.5-9.0}{\rm M}_{\odot }$| at |$z=0$| in the TNG50-1 cosmological simulation. The majority of the low-mass galaxies in TNG50-1 are dispersion-dominated, consistent with observations. By tracing the evolutionary histories of simulated low-mass galaxies, we find that while most stars form in rotating cold gas discs, the orientation of the star-forming discs relative to the galaxies may evolve with cosmic time. If the cold gas disc remains aligning with the galaxy during its evolution, stars formed at different times share the same rotational direction, leading to a rotation-dominated system. On the contrary, frequent misalignment of cold gas disc would result in a dispersion-dominated system. In addition, we also find that the two-body scattering can have a non-negligible numerical heating effect on the simulated galaxy morphology, especially at central regions of galaxies and for relatively low-mass galaxies. By comparing results of simulations with different resolutions, our results suggest that the simulated morphology of galaxies is roughly reliable when their number of stellar particles exceeds about |$10^{4}$|⁠ , and bulge morphology of galaxies can not be resolved robustly at the resolution level of TNG50-1. [ABSTRACT FROM AUTHOR]

Despite the evidence that supermassive black holes (SMBHs) co-evolve with their host galaxy, and that most of the growth of these SMBHs occurs via merger-free processes, the underlying mechanisms which drive this secular co-evolution are poorly understood. We investigate the role that both strong and weak large-scale galactic bars play in mediating this relationship. Using 48 871 disc galaxies in a volume-limited sample from Galaxy Zoo DESI, we analyse the active galactic nucleus (AGN) fraction in strongly barred, weakly barred, and unbarred galaxies up to |$z = 0.1$| over a range of stellar masses and colours. After controlling for stellar mass and colour, we find that the optically selected AGN fraction is |$31.6 \pm 0.9$| per cent in strongly barred galaxies, |$23.3 \pm 0.8$| per cent in weakly barred galaxies, and |$14.2 \pm 0.6$| per cent in unbarred disc galaxies. These are highly statistically robust results, strengthening the tantalizing results in earlier works. Strongly barred galaxies have a higher fraction of AGNs than weakly barred galaxies, which in turn have a higher fraction than unbarred galaxies. Thus, while bars are not required in order to grow an SMBH in a disc galaxy, large-scale galactic bars appear to facilitate AGN fuelling, and the presence of a strong bar makes a disc galaxy more than twice as likely to host an AGN than an unbarred galaxy at all galaxy stellar masses and colours. [ABSTRACT FROM AUTHOR]

Our ability to trace the star-forming molecular gas is important to our understanding of the Universe. We can trace this gas using CO emission, converting the observed CO intensity into the H |$_2$| gas mass of the region using the CO-to-H |$_2$| conversion factor (⁠|$X_{\rm{{\small CO}}}$|⁠). In this paper, we use simulations to study the conversion factor and the molecular gas within galaxies. We analysed a suite of simulations of isolated disc galaxies, ranging from dwarfs to Milky Way-mass galaxies, that were run using the fire-2 subgrid models coupled to the chimes non-equilibrium chemistry solver. We use the non-equilibrium abundances from the simulations, and we also compare to results using abundances assuming equilibrium, which we calculate from the simulation in post-processing. Our non-equilibrium simulations are able to reproduce the relation between CO and H |$_2$| column densities, and the relation between |$X_{\rm{{\small CO}}}$| and metallicity, seen within observations of the Milky Way. We also compare to the xCOLD GASS survey, and find agreement with their data to our predicted CO luminosities at fixed star formation rate. We also find the multivariate function used by xCOLD GASS overpredicts the H |$_2$| mass for our simulations, motivating us to suggest an alternative multivariate function of our fitting, though we caution that this fitting is uncertain due to the limited range of galaxy conditions covered by our simulations. We also find that the non-equilibrium chemistry has little effect on the conversion factor (<5 per cent) for our high-mass galaxies, though still affects the H |$_2$| mass and |$L_{\rm{{\small CO}}}$| by |$\approx$| 25  per cent. [ABSTRACT FROM AUTHOR]

We use resolved spectroscopy from MaNGA to investigate the significance of both local and global properties of galaxies to the cessation of star formation at kpc scales. Quenched regions are identified from a sample of isolated disk galaxies by a single-parameter criterion D n (4000) - log EW(H α) > 1.6 − log 2 = 1.3 , and are divided into gas-rich quenched regions (GRQRs) and gas-poor quenched regions according to the surface density of cold gas (Σgas). Both types of quenched regions tend to be hosted by non-AGN galaxies with relatively high mass (M * ≳ 1010 M ⊙) and red colors (NUV − r ≳ 3), as well as low star formation rate and high central density at fixed mass. They span wide ranges in other properties including structural parameters that are similar to the parent sample, indicating that the conditions responsible for quenching in gas-rich regions are largely independent of the global properties of galaxies. We train random forest classifiers and regressors for predicting quenching in our sample with 15 local/global properties. Σ* is the most important property for quenching, especially for GRQRs. These results strongly indicate the important roles of low-mass hot evolved stars, which are numerous and long-lived in quenched regions and can provide substantial radiation pressure to support the surrounding gas against gravitational collapse. The different feature importance for quenching, as found previously by A. F. L. Bluck et al., is partly due to the different definitions of quenched regions, particularly the different requirements on EW(H α). [ABSTRACT FROM AUTHOR]

Edge-on galaxies act as the best laboratories to understand the origin of thin and thick discs in galaxies. Measurement of spatially resolved stellar population properties in such galaxies, particularly age, metallicity, and [α/Fe], are crucial to understanding the formation and evolution of disc galaxies. Such measurements are made possible from stellar population model fits to deep integral field spectroscopic (IFU) observations of resolved galaxies. We utilize archival MUSE IFU observations of the edge-on galaxy ESO 544-27 to uncover the formation history of its thin and thick discs through its stellar populations. We find the thin disc of the galaxy is dominated by an old (>9 Gyr) low [α/Fe] metal-rich stellar population. Its outer thick disc is dominated by an old (>9 Gyr) high [α/Fe] metal-rich component that should have formed with higher star formation efficiency than the Milky Way thick disc. We thus find [α/Fe] dichotomy in ESO 544-27 with its thin and thick discs dominated by low and high [α/Fe] stellar populations, respectively. However, we also find a metal-rich younger (<2 Gyr old) stellar population in ESO 544-27. The galaxy was nearly quenched until its star formation was reignited recently first in the outer and inner thick disc (∼1 Gyr ago) and then in the thin disc (∼600 Myr ago). We thus find that both the low [α/Fe] thin and high [α/Fe] thick discs of ESO 544-27 are inhabited primarily by similarly old metal-rich stellar populations, a contrast to that of other galaxies with known thin and thick disc stellar population properties. [ABSTRACT FROM AUTHOR]

We measure the broad impact of galaxy structure on galaxy formation by examining the ongoing star formation and integrated star formation history as revealed through the stellar masses of galaxies at z < 7 based on JWST CEERS data from the Extended Groth Strip (EGS). Using the morphological catalog of 3965 visually classified JWST galaxies from Ferreira et al. (2023), we investigate the evolution of stars, and when they form, as a function of morphological type as well as galaxies classified as passive and starburst through spectral energy distributions. Although disc galaxies dominate the structures of galaxies at z < 7, we find that these discs are in general either 'passive', or on the main sequence of star formation, and do not contain a large population of starburst galaxies. We also find no significant correlation between morphological type and the star formation rate or colours of galaxies at z < 7. In fact, we find that the morphologically classified 'spheroids' tend to be blue and are not found to be predominately passive systems at z > 1.5. We also find that the stellar mass function for disc galaxies does not evolve significantly during this time, whereas other galaxy types, such as the peculiar population, evolve dramatically, declining at lower redshifts. This indicates that massive peculiars are more common at higher redshifts. We further find that up to z ∼ 7, the specific star formation rate (sSFR) does not vary with visual morphology, but strongly depends on stellar mass and internal galaxy mass density. This demonstrates that at early epochs galaxy assembly is a mass-driven, rather than a morphologically driven process. Quenching of star formation is therefore a mass-dominated process throughout the universe's history, likely due to the presence of supermassive black holes. [ABSTRACT FROM AUTHOR]

Based on the integral field unit (IFU) data from Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, we develop a new method to select galaxies with biconical ionized structures, building a sample of 142 edge-on biconical ionized galaxies. We classify these 142 galaxies into 81 star-forming galaxies, 31 composite galaxies, and 30 AGNs (consisting of 23 Seyferts and 7 LI(N)ERs) according to the [N  ii ]-BPT diagram. The star-forming bicones have bar-like structures while AGN bicones display hourglass structures, and composite bicones exhibit transitional morphologies between them due to both black hole and star formation activities. Star-forming bicones have intense star formation activities in their central regions, and the primary driver of biconical structures is the central star formation rate surface density. The lack of difference in the strength of central black hole activities (traced by dust attenuation corrected [O  iii ]λ5007 luminosity and Eddington ratio) between Seyfert bicones and their control samples can be naturally explained as that the accretion disc and the galactic disc are not necessarily coplanar. Additionally, the biconical galaxies with central LI(N)ER-like line ratios are edge-on disc galaxies that show strong central dust attenuation. The radial gradients of H α surface brightness follow the r −2.35 relation, roughly consistent with r −2 profile, which is expected in the case of photoionization by a central point-like source. These observations indicate obscured AGNs or AGN echoes as the primary drivers of biconical structures in LI(N)ERs. [ABSTRACT FROM AUTHOR]

We measure the age–velocity relationship from the lag between ionized gas and stellar tangential speeds in ∼500 nearby disc galaxies from MaNGA in Sloan Digital Sky Survey IV (SDSS-IV). Selected galaxies are kinematically axisymmetric. Velocity lags are asymmetric drift, seen in the Milky Way's (MW) solar neighbourhood and other Local Group galaxies; their amplitude correlates with stellar population age. The trend is qualitatively consistent in rate (⁠|$\dot{\sigma }$|⁠) with a simple power-law model where σ ∝  t b that explains the dynamical phase-space stratification in the solar neighbourhood. The model is generalized based on disc dynamical times to other radii and other galaxies. We find in-plane radial stratification parameters σ0, r (dispersion of the youngest populations) in the range of 10–40 km s−1  and 0.2 < br < 0.5 for MaNGA galaxies. Overall, br increases with galaxy mass, decreases with radius for galaxies above 10.4 dex (M⊙) in stellar mass, but is ∼constant with radius at lower mass. The measurement scatter indicates the stratification model is too simple to capture the complexity seen in the data, unsurprising given the many possible astrophysical processes that may lead to stellar population dynamical stratification. None the less, the data show dynamical stratification is broadly present in the galaxy population, with systematic trends in mass and density. The amplitude of the asymmetric drift signal is larger for the MaNGA sample than the MW, and better represented in the mean by what is observed in the discs of M31 and M33. Either typical discs have higher surface-density or, more likely, are dynamically hotter (hence thicker) than the MW. [ABSTRACT FROM AUTHOR]

We analyse a sample of massive disc galaxies selected from the fourth-generation Sloan Digital Sky Survey/Mapping Nearby Galaxies at Apache Point Observatory survey to investigate how the evolution of these galaxies depends on their stellar and halo masses. We applied a semi-analytic spectral fitting approach to the data from different regions in the galaxies to derive several of their key physical properties. From the best-fitting model results, together with direct observables such as morphology, colour, and the Mg b /〈Fe〉 index ratio measured within 1 R e, we find that for central galaxies both their stellar and halo masses have a significant influence in their evolution. For a given halo mass, galaxies with higher stellar mass accumulate their stellar mass and become chemically enriched earlier than those with smaller stellar mass. Furthermore, at a given stellar mass, galaxies living in more massive haloes have longer star formation time-scales and are delayed in becoming chemically enriched. In contrast, the evolution of massive satellite galaxies is mostly determined by their stellar mass. The results indicate that both the assembled halo mass and the halo assembly history impact the evolution of central galaxies. Our spatially resolved analysis indicates that only the galaxy properties in the central region (0.0–0.5 R e) show the dependencies described above. This fact supports a halo-driven formation scenario since the galaxies' central regions are more likely to contain old stars formed along with the halo itself, keeping a memory of the halo formation process. [ABSTRACT FROM AUTHOR]

The extraplanar diffuse ionized gas (eDIG) represents ionized gases traced by optical/UV lines beyond the stellar extent of galaxies. We herein introduce a novel multislit narrow-band spectroscopy method to conduct spatially resolved spectroscopy of the eDIG around a sample of nearby edge-on disk galaxies (eDIG-CHANGES). In this paper, we introduce the project design and major scientific goals, as well as a pilot study of NGC 3556 (M108). The eDIG is detected to a vertical extent of a few kiloparsecs above the disk, comparable to the X-ray and radio images. We do not see significant vertical variation of the [N ii ]/H α line ratio. A rough examination of the pressure balance between different circumgalactic medium phases indicates the magnetic field is in a rough pressure balance with the X-ray emitting hot gas and may play an important role in the global motion of both the eDIG and the hot gas in the lower halo. At the location of an Hubble Space Telescope/Cosmic Origins Spectrograph observed UV bright background active galactic nucleus ∼29 kpc from the center of NGC 3556, the magnetic pressure is much lower than that of the hot gas and the ionized gas traced by UV absorption lines, although the extrapolation of the pressure profiles may cause some biases in this comparison. By comparing the position–velocity diagrams of the optical and CO lines, we also find the dynamics of the two gas phases are consistent with each other, with no evidence of a global inflow/outflow and a maximum rotation velocity of ∼150 km s−1. [ABSTRACT FROM AUTHOR]

CHANG-ES (Continuum Halos in Nearby Galaxies—an EVLA Survey) is an ambitious project to target 35 nearby disk galaxies that are edge-on to the line of sight. The orientation permits both the disk and halo regions to be studied. The observations were initially at 1.5 GHz (L-band) and 6.0 GHz (C-band) in a variety of VLA array configurations, and in all four Stokes parameters, which allowed for spatially resolved images in total intensity plus polarization. The inclusion of polarization is unique to an edge-on galaxy survey and reveals the galaxies' halo magnetic fields. This paper will summarize the results to date, some of which are new phenomena, never seen prior to CHANG-ES. For example, we see that 'X-type' fields, as well as rotation measure reversals, are common features of spiral galaxies. Further observations at 3.0 GHz (S-band) as well as future scientific opportunities will also be described. [ABSTRACT FROM AUTHOR]

Barred galaxies constitute about two-thirds of observed disc galaxies. Bars affect not only the mass distribution of gas and stars but also that of the dark matter. An elongation of the inner dark matter halo is known as the halo bar. We aim to characterize the structure of the halo bars, with the goal of correlating them with the properties of the stellar bars. We use a suite of simulated galaxies with various bar strengths, including gas and star formation. We quantify the strengths, shapes, and densities of these simulated stellar bars. We carry out numerical experiments with frozen and analytic potentials in order to understand the role played by a live responsive stellar bar. We find that the halo bar generally follows the trends of the disc bar. The strengths of the halo and stellar bars are tightly correlated. Stronger bars induce a slight increase in dark matter density within the inner halo. Numerical experiments show that a non-responsive frozen stellar bar would be capable of inducing a dark matter bar, but it would be weaker than the live case by a factor of roughly two. [ABSTRACT FROM AUTHOR]

Empirical studies of the relationship between baryonic matter in galaxies and the gravitational potential of their host haloes are important to constrain our theoretical framework for galaxy formation and evolution. One such relation, between the atomic hydrogen (H  i) mass of central galaxies (M H  i ,c) and the total mass of their host haloes (M halo), has attracted significant interest in the last few years. In this work, we use the extended GALEX Arecibo SDSS Survey to examine the scatter of the H  i– halo mass relation for a representative sample of central galaxies. Our findings reveal a flat median relation at |$\mathrm{ log}_{10}(M_{\rm {H\,{\small I},c}}/\rm {M}_{\odot }) \approx 9.40$|⁠ , across |$11.1 \lt \mathrm{ log}_{10}(M_{\rm {halo}}/{{\rm M}_{\odot }}) \lt 14.1$|⁠. This flat relation stems from the statistical dominance of star-forming, disc galaxies at low M halo in combination with the increasing prevalence of passive, high stellar concentration systems at higher M halo. The scatter of this relation and the stellar specific angular momentum of centrals have a strong link (Spearman's rank correlation coefficient ≥0.5). Comparisons with simulations suggest that the kinematic state of host haloes may be primarily driving this scatter. Our findings highlight that the H  i –halo mass parameter space is too complex to be completely represented by simple median or average relations and we show that tensions with previous works are most likely due to selection biases. We recommend that future observational studies, and their comparisons with theoretical models, bin central galaxies also by their secondary properties to enable a statistically robust understanding of the processes regulating the cold gas content within central galaxies of dark matter haloes. [ABSTRACT FROM AUTHOR]

The presence of a stellar bar in a disc galaxy indicates that the galaxy hosts in its main part a dynamically settled disc and that bar-driven processes are taking place in shaping its evolution. Studying the cosmic evolution of the bar fraction in disc galaxies is therefore essential to understand galaxy evolution in general. Previous studies have found, using the Hubble Space Telescope (HST), that the bar fraction significantly declines from the local Universe to redshifts near one. Using the first four pointings from the JWST Cosmic Evolution Early Release Science Survey and the initial public observations for the Public Release Imaging for Extragalactic Research, we extend the studies of the bar fraction in disc galaxies to redshifts 1 ≤ |$z$| ≤ 3, that is, for the first time beyond redshift two. We only use galaxies that are also present in the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey on the Extended Groth Strip and Ultra Deep Survey HST observations. An optimized sample of 368 close-to-face-on galaxies is visually classified to find the fraction of bars in disc galaxies in two redshift bins: 1 ≤ |$z$| ≤ 2 and 2 < |$z$| ≤ 3. The bar fraction decreases from |$\approx 17.8^{+ 5.1}_{- 4.8}$|  per cent to |$\approx 13.8^{+ 6.5}_{- 5.8}$|  per cent (from the lower to the higher redshift bin), but is about twice the bar fraction found using bluer HST filters. Our results show that bar-driven evolution might commence at early cosmic times and that dynamically settled discs are already present at a lookback time of ∼11 Gyr. [ABSTRACT FROM AUTHOR]

Observed accretion rates onto the Milky Way and other local spirals fall short of that required to sustain star formation for cosmological timescales. A potential avenue for this unseen accretion is a rotating inflow in the volume-filling hot phase (⁠|$\sim 10^6\, {\rm K}$|⁠) of the circumgalactic medium (CGM), as suggested by some cosmological simulations. Using hydrodynamic simulations and a new analytic solution valid in the slow-rotation limit, we show that a hot inflow spins up as it approaches the galaxy, while remaining hot, subsonic, and quasi-spherical. Within the radius of angular momentum support (⁠|$\sim 15\, {\rm kpc}$| for the Milky Way) the hot flow flattens into a disc geometry and then cools from |$\sim 10^6$| to |$\sim 10^4\, {\rm K}$| at the disc–halo interface. Cooling affects all hot gas, rather than just a subset of individual gas clouds, implying that accretion via hot inflows does not rely on local thermal instability in contrast with 'precipitation' models for galaxy accretion. Prior to cooling and accretion the inflow completes ≈ t cool/ t ff radians of rotation, where t cool/ t ff is the cooling time to free-fall time ratio in hot gas immediately outside the galaxy. The ratio t cool/ t ff may thus govern the development of turbulence and enhancement of magnetic fields in gas accreting onto low-redshift spirals. We show that if rotating hot inflows are common in Milky-Way-size disc galaxies, as predicted, then signatures of the expected hot gas rotation profile should be observable with X-ray telescopes and fast radio burst surveys. [ABSTRACT FROM AUTHOR]

The radial acceleration relation (RAR) connects the total gravitational acceleration of a galaxy at a given radius, a tot(r), with that accounted for by baryons at the same radius, a bar(r). The shape and tightness of the RAR for rotationally-supported galaxies have characteristics in line with MOdified Newtonian Dynamics (MOND) and can also arise within the cosmological constant + cold dark matter (ΛCDM) paradigm. We use zoom simulations of 20 galaxies with stellar masses of M ⋆ ≃ 107–11 M⊙ to study the RAR in the FIRE-2 simulations. We highlight the existence of simulated galaxies with non-monotonic RAR tracks that 'hook' down from the average relation. These hooks are challenging to explain in Modified Inertia theories of MOND, but naturally arise in all of our ΛCDM-simulated galaxies that are dark-matter dominated at small radii and have feedback-induced cores in their dark matter haloes. We show, analytically and numerically, that downward hooks are expected in such cored haloes because they have non-monotonic acceleration profiles. We also extend the relation to accelerations below those traced by disc galaxy rotation curves. In this regime, our simulations exhibit 'bends' off of the MOND-inspired extrapolation of the RAR, which, at large radii, approach a tot ≈ a bar/ f b, where f b is the cosmic baryon fraction. Future efforts to search for these hooks and bends in real galaxies will provide interesting tests for MOND and ΛCDM. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

It is well established that stellar discs are destabilized by sharp features in their phase space, driving recurrent spiral modes. We explore the extent to which surface-density breaks in disc galaxies – which represent sharp changes in the gradient of the disc density – drive new spiral modes. We employ linear perturbation theory to investigate how disc breaks alter the eigenmode spectrum of an otherwise pure exponential disc. We find that the presence of a density break gives rise to a set of new, vigorously growing, modes. For a given multiplicity, these edge modes occur in pairs, with closely separated resonances between each pair. The growth rate of edge modes decreases when the break is weakened or moved outward to lower-density regions of the disc. Both down- and up-bending profiles excite edge modes, whose origin can be best understood via the gravitational torques they exert on the underlying disc. When the profile is down-bending (Type II) the faster growing mode is the inner one while in the up-bending (Type III) case the outer mode is faster growing. In both cases, the faster growing mode has a corotation almost coincident with the break. We show that the torques of the edge modes tend to smoothen the break. [ABSTRACT FROM AUTHOR]

We present N -body simulations of the process of bulge formation in disc galaxies due to inward migration of massive stellar clumps. The process is accompanied by dark halo heating, with a quasi-isothermal core replacing the initial central density cusp, transforming an initially dark matter dominated central region into a baryon dominated one. The characteristics of the clumps are chosen to be compatible with low redshift observations of stellar clumps in DYNAMO-HST galaxies, which may be relatively long lived in terms of being robust against internal starburst-instigated disruption. We thus test for disruption due to tidal stripping using different clump internal radial profiles; Plummer, Hernquist, and Jaffe, in ascending order of see per central density profile. Our calculations predict that in order for clump migration to be effective in building galactic bulges and dark halo cores, steeply increasing central clump profiles, or a less massive or less concentrated haloes, are preferred. The dependence on such factors may contribute to the diversity in observed total mass distributions and resulting rotation curves in galaxies. When the process is most efficient, a 'bulge-halo conspiracy', with a singular isothermal total density akin to that observed bright galaxies, results. [ABSTRACT FROM AUTHOR]

Disk galaxies have a typical dimension of a few tens of kiloparsecs. It follows from the theory of general relativity that any signal originating from the galactic center will be noticed at the outskirts of the galaxy only tens of thousands of years later. This retardation effect, however, is absent in modelling used to calculate rotation curves throughout entire galaxies and their external gas. The considerable differences between Newtonian predictions and observed velocities are currently removed either by assuming dark matter or by suggesting various modifications to the laws of gravity, MOND being a long standing alternative to Newtonian gravity. In previous papers we have shown that by applying general relativity in a rigorous fashion, without neglecting retardation, one can explain the rotational velocities of galactic matter without modifying gravity or adding dark matter. Moreover, it was shown that dark matter effects, as they appear in gravitational lensing, the Tully-Fisher relation, and mass estimations based on the virial theorem could also be explained as retarded-gravity effects. It must be noted, however, that the proposed theory relies on the existence of a mass flow (of a changing rate) into the galaxy; a requirement that was not directly observed. In the original paper on the subject only one galaxy (M33) was analysed in detail. This was later amended with a published study of eleven galaxies. Here we give a more comprehensive retardation analysis of 143 galaxies of different types from the SPARC Galaxy collection. We show that in most cases we obtain very accurate fits to the data. [ABSTRACT FROM AUTHOR]

It is widely believed that resonant orbits play an important role in formation and evolution of bars and large-scale spirals in galaxy discs. These resonant orbits have been studied in a number of specific potentials, often with an imposed bar component. In this paper, I show that families of resonant (e.g. two-dimensional x 1) orbits of differing eccentricities can be excited at a common pattern speed, in a variety of axisymmetric potentials. These families only exist over finite ranges of frequency in most of these potentials. Populations of such resonant eccentric orbits (REOs) can provide the backbone of both bars and spirals. At each frequency in the allowed range, there is a maximum eccentricity, beyond which the REOs generically become quasi-stable (or 'sticky'), then unstable (or chaotic), as the eccentricity increases, at values that depend on the potential and the orbit frequency. Sticky and chaotic orbits have been extensively studied recently with invariant/unstable manifolds in a variety of phase planes; however, it is found that studying them as a function of eccentricity and pattern speed provides a particularly useful framework for classifying them and their stability transitions. The characteristics of these orbit families depend on the galaxy potential and the pattern speed, and as backbones of bars and spirals can help understand a number of observed or predicted regularities. These include: the size and speed of bars in different potentials, the range of pattern speeds and wind-up rates in spirals within galaxy discs, and constraints wave growth. [ABSTRACT FROM AUTHOR]

Spiral density waves can arise in galactic discs as linear instabilities of the underlying stellar distribution function. Such an instability grows exponentially in amplitude at some fixed growth rate β before saturating non-linearly. However, the mechanisms behind saturation, and the resulting saturated spiral amplitude, have received little attention. Here, we argue that one important saturation mechanism is the non-linear trapping of stars near the spiral's corotation resonance. Under this mechanism, we show analytically that an m -armed spiral instability will saturate when the libration frequency of resonantly trapped orbits reaches |$\omega _\mathrm{lib} \sim \mathrm{a\, \, few}\times m^{1/2} \beta$|⁠. For a galaxy with a flat rotation curve, this implies a maximum relative spiral surface density |$\vert \delta \Sigma /\Sigma _0\vert \sim \mathrm{a\, \, few} \times (\beta /\Omega _\mathrm{p})^2 \cot \alpha$|⁠ , where Ωp is the spiral pattern speed and α is its pitch angle. This result is in reasonable agreement with recent N -body simulations, and suggests that spirals driven by internally-generated instabilities reach relative amplitudes of at most a few tens of per cent; higher amplitude spirals, like in M51 and NGC 1300, are likely caused by very strong bars and/or tidal perturbations. [ABSTRACT FROM AUTHOR]

Recent puzzling observations, such as the H 0 tension, large-scale anisotropies, and massive disk galaxies at high redshifts, have been challenging the standard cosmological model. While one possible explanation is that the standard model is incomplete, other theories are based on the contention that the redshift model as a distance indicator might be biased. These theories can explain the recent observations, but they are challenged by the absence of a direct empirical reproducible observation that the redshift model can indeed be inconsistent. Here, I describe a simple experiment that shows that the spectra of galaxies depend on their rotational velocity relative to the rotational velocity of the Milky Way. Moreover, it shows that the redshift of galaxies that rotate in the opposite direction relative to the Milky Way is significantly smaller compared with the redshift of galaxies that rotate in the same direction relative to the Milky Way (p < 0.006). Three different datasets were used independently, each one was prepared in a different manner, and all of them showed similar redshift bias. A fourth dataset of galaxies from the Southern Galactic pole was also analyzed and shows similar results. All four datasets are publicly available. While a maximum average z difference of ∼0.012 observed with galaxies of relatively low redshift (z < 0.25) is not extreme, the bias is consistent and canpotentially lead to explanations to puzzling observations such as the H 0 tension. [ABSTRACT FROM AUTHOR]

Theoretical models of spiral arms suggest that the spiral arms provoke a vertical bulk motion in disc stars. By analysing the breathing motion, a coherent asymmetric vertical motion around the mid-plane of the Milky Way disc, with Gaia DR3, we found that a compressing breathing motion presents along the Local arm. On the other hand, with an N -body simulation of an isolated Milky Way-like disc galaxy, we found that the transient and dynamic spiral arms induce compressing breathing motions when the arms are in the growth phase, while the expanding breathing motion appears in the disruption phase. The observed clear alignment of the compressing breathing motion with the Local arm is similar to what is seen in the growth phase of the simulated spiral arms. Hence, we suggest that the Local arm's compressing breathing motion can be explained by the Local arm being in the growth phase of a transient and dynamic spiral arm. We also identified the tentative signatures of the expanding breathing motion associated with the Perseus arm and also the Outer arm coinciding with the compressing breathing motion. This may infer that the Perseus and Outer arms are in the disruption and growth phases, respectively. [ABSTRACT FROM AUTHOR]

Magnetic fields are ubiquitous in the Universe. Recently, cosmological simulations of galaxies have successfully begun to incorporate magnetic fields and their evolution in galaxies and their haloes. However, so far they have mostly focused on Milky Way-like galaxies. Here, we analyse a sample of high-resolution cosmological zoom simulations of disc galaxies in haloes with mass |${M}_\rm {200c}$| from |$10^{10}$| to |$10^{13}\, \rm {M}_\odot$|⁠ , simulated with the Auriga galaxy formation model. We show that with sufficient numerical resolution the magnetic field amplification and saturation is converged. The magnetic field strength reaches equipartition with turbulent energy density for galaxies in haloes with |${M}_\rm {200c}\gtrsim 10^{11.5}\, \mathrm{M_\odot }$|⁠. For galaxies in less massive haloes, the magnetic field strength saturates at a fraction of equipartition that decreases with decreasing halo mass. For our lowest mass haloes, the magnetic field saturates significantly below 10  per cent of equipartition. We quantify the resolution we need to obtain converged magnetic field strengths and discuss our resolution requirements also in the context of the IllustrisTNG cosmological box simulations. We show that, at z  = 0, rotation-dominated galaxies in our sample exhibit for the most part an ordered large-scale magnetic field, with fewer field reversals in more massive galaxies. Finally, we compare the magnetic fields in our cosmological galaxies at z  = 0 with simulations of isolated galaxies in a collapsing halo set-up. Our results pave the way for detailed studies of cosmic rays and other physical processes in similar cosmological galaxy simulations that crucially depend on the strength and structure of magnetic fields. [ABSTRACT FROM AUTHOR]

The merger history of a galaxy is thought to be one of the major factors determining its internal dynamics, with galaxies having undergone different types or mergers (e.g. dry, minor, or major mergers) predicted to show different dynamical properties. We study the instantaneous orbital distribution of galaxies in the eagle simulation, colouring the orbits of the stellar particles by their stellar age, in order to understand whether stars form in particular orbits (e.g. in a thin or thick disc). We first show that eagle reproduces well the observed stellar mass fractions in different stellar orbital families as a function of stellar mass and spin parameter at z  = 0. We find that the youngest stars reside in a thin disc component that can extend to the very inner regions of galaxies, and that older stars have warmer orbits, with the oldest ones showing orbits consistent with both hot and counter-rotating classifications, which is consistent with the trend found in the Milky Way and other disc galaxies. We also show that counter-rotating orbits trace galaxy mergers – in particular dry mergers, and that in the absence of mergers, counter-rotating orbits can also be born from highly misaligned gas accretion that leads to star formation. [ABSTRACT FROM AUTHOR]

We investigate the relationship between the baryonic angular momentum and mass for a sample of 36 isolated disc galaxies with resolved neutral hydrogen (H  i) kinematics and infrared Wide-field Infrared Survey Explorer photometry drawn from – and representative in terms of morphologies, stellar masses, and H  i -to-star fraction of – the carefully constructed Analysis of the interstellar Medium in Isolated GAlaxies (AMIGA) sample of isolated galaxies. Similarly to previous studies performed on non-isolated galaxies, we find that the relation is well described by a power law |$j_{\rm bar} \propto M_{\rm bar}^\alpha$|⁠. We also find a slope of α = 0.54 ± 0.08 for the AMIGA galaxies, in line with previous studies in the literature; however, we find that the specific angular momenta of the AMIGA galaxies are on average higher than those of non-isolated galaxies in the literature. This is consistent with theories stipulating that environmental processes involving galaxy–galaxy interaction are able to impact the angular momentum content of galaxies. However, no correlation was found between the angular momentum and the degree of isolation, suggesting that there may exist a threshold local number density beyond which the effects of the environment on the angular momentum become important. [ABSTRACT FROM AUTHOR]

Large-scale cosmological galaxy formation simulations typically prevent gas in the interstellar medium (ISM) from cooling below |$\approx 10^4\, \mathrm{K}$|⁠. This has been motivated by the inability to resolve the Jeans mass in molecular gas (⁠|$\ll 10^5\, \mathrm{M}_{\odot }$|⁠) which would result in undesired artificial clumping. We show that the classical Jeans criteria derived for Newtonian gravity are not applicable in the simulated ISM if the spacing of resolution elements representing the dense ISM is below the gravitational force softening length and gravity is therefore softened and not Newtonian. We re-derive the Jeans criteria for softened gravity in Lagrangian codes and use them to analyse gravitational instabilities at and below the hydrodynamical resolution limit for simulations with adaptive and constant gravitational softening lengths. In addition, we define criteria for which a numerical runaway collapse of dense gas clumps can occur caused by oversmoothing of the hydrodynamical properties relative to the gravitational force resolution. This effect is illustrated using simulations of isolated disc galaxies with the smoothed particle hydrodynamics code swift. We also demonstrate how to avoid the formation of artificial clumps in gas and stars by adjusting the gravitational and hydrodynamical force resolutions. [ABSTRACT FROM AUTHOR]

Ring galaxies are rare in the Universe. A head-on or off-centre collision between an intruder galaxy and a disc galaxy leads to a collisional ring galaxy (CRG) when the intruder-to-target mass ratio (ITMR) is over 0.1 in Newtonian dynamics. Within the framework of Milgromian dynamics, the strong equivalence principle is violated due to the presence of an external field. When two galaxies collide towards each other, the dynamical mass of the phantom dark halo in a galaxy is suppressed by the external field induced by the other galaxy. As a consequence of such suppression, the gravitational perturbation for the target galaxy introduced by a collision is weakened. In this case, a minor collision may not be capable of generating a CRG. In this work, we address this issue and perform a series of numerical simulations of collisions by tuning the values of ITMR. We find that the critical ITMR is 0.5 in MOND, which is much larger than that in Newtonian dynamics. The observed massive ring galaxies, such as Arp 147, can be effectively interpreted by CRGs in Milgrom's Modified Newtonian Dynamics (MOND). This interpretation does not necessitate the presence of dark matter halos for either the target or intruder galaxies. Moreover, for a larger inclination angle or a larger impact parameter, the off-centred ring structure is fainter. The larger critical ITMR indicates that it is harder to form a CRG in MOND. To account for the observed ring structures of the NGC 922-like galaxies in MOND, it is necessary to invoke other scenarios than a dry minor collision. [ABSTRACT FROM AUTHOR]

Stellar kinematics provides a window into the gravitational field, and therefore into the distribution of all mass, including dark matter. Deep Potential is a method for determining the gravitational potential from a snapshot of stellar positions in phase space, using mathematical tools borrowed from deep learning to model the distribution function and solve the Collisionless Boltzmann equation. In this work, we extend the Deep Potential method to rotating systems, and then demonstrate that it can accurately recover the gravitational potential, density distribution, and pattern speed of a simulated barred disc galaxy, using only a frozen snapshot of the stellar velocities. We demonstrate that we are able to recover the bar pattern speed to within |$15 \,\rm {per\, cent}$| in our simulated galaxy using stars in a 4 kpc subvolume centred on a Solar-like position, and to within |$20 \,\rm{per\,cent}$| in a 2 kpc subvolume. In addition, by subtracting the mock 'observed' stellar density from the recovered total density, we are able to infer the radial profile of the dark matter density in our simulated galaxy. This extension of Deep Potential is an important step in allowing its application to the Milky Way, which has rotating features, such as a central bar and spiral arms, and may moreover provide a new method of determining the pattern speed of the Milky Way bar. [ABSTRACT FROM AUTHOR]