Showing 1–38 of 38 results for author: Nguyen, D D

Hot Gas Outflow Properties of the Starburst Galaxy NGC 4945

Authors: Natalia Porraz Barrera, Sebastian Lopez, Laura A. Lopez, Adi Foord, Dustin D. Nguyen, Todd A. Thompson, Smita Mathur, Alberto D. Bolatto

Abstract: We analyze 330 ks of {\it Chandra} X-ray imaging and spectra of the nearby, edge-on starburst and Seyfert Type 2 galaxy NGC 4945 to measure the hot gas properties along the galactic outflows. We extract and model spectra from 15 regions extending from $-$0.55 kpc to $+$0.85 kpc above and below the galactic disk to determine the best-fit parameters and metal abundances. We find that the hot gas tem… ▽ More We analyze 330 ks of {\it Chandra} X-ray imaging and spectra of the nearby, edge-on starburst and Seyfert Type 2 galaxy NGC 4945 to measure the hot gas properties along the galactic outflows. We extract and model spectra from 15 regions extending from $-$0.55 kpc to $+$0.85 kpc above and below the galactic disk to determine the best-fit parameters and metal abundances. We find that the hot gas temperatures and number densities peak in the central regions and decrease along the outflows. These profiles are inconsistent with a spherical, adiabatically-expanding wind model, suggesting the need to include mass loading and/or a non-spherical outflow geometry. We estimate the mass outflow rate of the hot wind to be $1.6\:M_{\odot}~\rm{yr}^{-1}$. Emission from charge exchange is detected in the northern outflow, and we estimate it contributes 12\% to the emitted, broad-band ($0.5-7$ keV) X-ray flux. △ Less

Submitted 11 April, 2024; v1 submitted 13 December, 2023; originally announced December 2023.

Comments: 11 pages, 7 figures; submitted (13 December 2023), accepted ApJ (11 April 2023)

Neural ODEs as a discovery tool to characterize the structure of the hot galactic wind of M82

Authors: Dustin D. Nguyen, Yuan-Sen Ting, Todd A. Thompson, Sebastian Lopez, Laura A. Lopez

Abstract: Dynamic astrophysical phenomena are predominantly described by differential equations, yet our understanding of these systems is constrained by our incomplete grasp of non-linear physics and scarcity of comprehensive datasets. As such, advancing techniques in solving non-linear inverse problems becomes pivotal to addressing numerous outstanding questions in the field. In particular, modeling hot g… ▽ More Dynamic astrophysical phenomena are predominantly described by differential equations, yet our understanding of these systems is constrained by our incomplete grasp of non-linear physics and scarcity of comprehensive datasets. As such, advancing techniques in solving non-linear inverse problems becomes pivotal to addressing numerous outstanding questions in the field. In particular, modeling hot galactic winds is difficult because of unknown structure for various physical terms, and the lack of \textit{any} kinematic observational data. Additionally, the flow equations contain singularities that lead to numerical instability, making parameter sweeps non-trivial. We leverage differentiable programming, which enables neural networks to be embedded as individual terms within the governing coupled ordinary differential equations (ODEs), and show that this method can adeptly learn hidden physics. We robustly discern the structure of a mass-loading function which captures the physical effects of cloud destruction and entrainment into the hot superwind. Within a supervised learning framework, we formulate our loss function anchored on the astrophysical entropy ($K \propto P/ρ^{5/3}$). Our results demonstrate the efficacy of this approach, even in the absence of kinematic data $v$. We then apply these models to real Chandra X-Ray observations of starburst galaxy M82, providing the first systematic description of mass-loading within the superwind. This work further highlights neural ODEs as a useful discovery tool with mechanistic interpretability in non-linear inverse problems. We make our code public at this GitHub repository (https://github.com/dustindnguyen/2023_NeurIPS_NeuralODEs_M82). △ Less

Submitted 28 November, 2023; v1 submitted 3 November, 2023; originally announced November 2023.

Comments: 9 Pages, 2 Figures, Accepted at the NeurIPS 2023 workshop on Machine Learning and the Physical Sciences

Numerical modeling of thermal dust polarization from aligned grains in the envelope of evolved stars with updated POLARIS

Authors: Bao Truong, Thiem Hoang, Nguyen Chau Giang, Pham Ngoc Diep, Dieu D. Nguyen, Nguyen Bich Ngoc

Abstract: Magnetic fields are thought to influence the formation and evolution of evolved star envelopes. Thermal dust polarization from magnetically aligned grains is potentially a powerful tool for probing magnetic fields and dust properties in these circumstellar environments. In this paper, we present numerical modeling of thermal dust polarization from the envelope of IK Tau using the magnetically enha… ▽ More Magnetic fields are thought to influence the formation and evolution of evolved star envelopes. Thermal dust polarization from magnetically aligned grains is potentially a powerful tool for probing magnetic fields and dust properties in these circumstellar environments. In this paper, we present numerical modeling of thermal dust polarization from the envelope of IK Tau using the magnetically enhanced radiative torque (MRAT) alignment theory implemented in our updated POLARIS code. Due to the strong stellar radiation field, the minimum size required for RAT alignment of silicate grains is $\sim 0.005 - 0.05\,\rmμm$. Additionally, ordinary paramagnetic grains can achieve perfect alignment by MRAT in the inner regions of $r < 500\,\rm au$ due to stronger magnetic fields of $B\sim 10$ mG - 1G, producing thermal dust polarization degree of $\sim 10\,\%$. The polarization degree can be enhanced to $\sim 20-40\%$ for grains with embedded iron inclusions. We also find that the magnetic field geometry affects the alignment size and the resulting polarization degree due to the projection effect in the plane-of-sky. We also study the spectrum of polarized thermal dust emission and find the increased polarization degree toward $λ> 50\,\rmμm$ due to the alignment of small grains by MRAT. Furthermore, we investigate the impact of rotational disruption by RATs (RAT-D) and find the RAT-D effect cause a decrease in the dust polarization fraction. Finally, we compare our numerical results with available polarization data observed by SOFIA/HAWC+ for constraining dust properties, suggesting grains are unlikely to have embedded iron clusters and might have slightly elongated shapes. Our modeling results suggest further observational studies at far-infrared/sub-millimeter wavelengths to understand the properties of magnetic fields and dust in AGB envelopes. △ Less

Submitted 2 August, 2023; originally announced August 2023.

Comments: 27 pages, 23 figures, 1 table, to be submitted

Highly-mass-loaded hot galactic winds are unstable to cool filament formation

Authors: Dustin D. Nguyen, Todd A. Thompson, Evan E. Schneider, Ashley P. Tarrant

Abstract: When cool clouds are ram-pressure accelerated by a hot supersonic galactic wind, some of the clouds may be shredded by hydrodynamical instabilities and incorporated into the hot flow. Recent one-dimensional steady-state calculations show how cool cloud entrainment directly affects the bulk thermodynamics, kinematics, and observational characteristics of the hot gas. In particular, mass-loading dec… ▽ More When cool clouds are ram-pressure accelerated by a hot supersonic galactic wind, some of the clouds may be shredded by hydrodynamical instabilities and incorporated into the hot flow. Recent one-dimensional steady-state calculations show how cool cloud entrainment directly affects the bulk thermodynamics, kinematics, and observational characteristics of the hot gas. In particular, mass-loading decelerates the hot flow and changes its entropy. Here, we investigate the stability of planar and spherical mass-loaded hot supersonic flows using both perturbation analysis and three-dimensional time-dependent radiative hydrodynamical simulations. We show that mass-loading is stable over a broad range of parameters and that the 1D time-steady analytic solutions exactly reproduce the 3D time-dependent calculations, provided that the flow does not decelerate sufficiently to become subsonic. For higher values of the mass-loading, the flow develops a sonic point and becomes thermally unstable, rapidly cooling and forming elongated dense cometary filaments. We explore the mass-loading parameters required to reach a sonic point and the radiative formation of these filaments. For certain approximations, we can derive simple analytic criteria. In general a mass-loading rate similar to the initial mass outflow rate is required. In this sense, the destruction of small cool clouds by a hot flow may ultimately spontaneously generate fast cool filaments, as observed in starburst superwinds. Lastly, we find that the kinematics of filaments is sensitive to the slope of the mass-loading function. Filaments move faster than the surrounding wind if mass-loading is over long distances whereas filaments move slower than their surroundings if mass-loading is abrupt. △ Less

Submitted 21 July, 2023; originally announced July 2023.

Comments: 12 pages, 15 figures, submitted to MNRAS (21 July 2023)

Neural Astrophysical Wind Models

Authors: Dustin D. Nguyen

Abstract: The bulk kinematics and thermodynamics of hot supernovae-driven galactic winds is critically dependent on both the amount of swept up cool clouds and non-spherical collimated flow geometry. However, accurately parameterizing these physics is difficult because their functional forms are often unknown, and because the coupled non-linear flow equations contain singularities. We show that deep neural… ▽ More The bulk kinematics and thermodynamics of hot supernovae-driven galactic winds is critically dependent on both the amount of swept up cool clouds and non-spherical collimated flow geometry. However, accurately parameterizing these physics is difficult because their functional forms are often unknown, and because the coupled non-linear flow equations contain singularities. We show that deep neural networks embedded as individual terms in the governing coupled ordinary differential equations (ODEs) can robustly discover both of these physics, without any prior knowledge of the true function structure, as a supervised learning task. We optimize a loss function based on the Mach number, rather than the explicitly solved-for 3 conserved variables, and apply a penalty term towards near-diverging solutions. The same neural network architecture is used for learning both the hidden mass-loading and surface area expansion rates. This work further highlights the feasibility of neural ODEs as a promising discovery tool with mechanistic interpretability for non-linear inverse problems. △ Less

Submitted 25 June, 2023; v1 submitted 20 June, 2023; originally announced June 2023.

Comments: 7 Pages, 4 Figures, Accepted at the ICML 2023 Workshop on Machine Learning for Astrophysics; v2) fixed typos

Effect of the initial mass function on the dynamical SMBH mass estimate in the nucleated early-type galaxy FCC 47

Authors: Sabine Thater, Mariya Lyubenova, Katja Fahrion, Ignacio Martín-Navarro, Prashin Jethwa, Dieu D. Nguyen, Glenn van de Ven

Abstract: Supermassive black holes (SMBHs) and nuclear star clusters (NSCs) co-exist in many galaxies. While the formation history of the black hole is essentially lost, NSCs preserve their evolutionary history imprinted onto their stellar populations and kinematics. Studying SMBHs and NSCs in tandem might help us to ultimately reveal the build-up of galaxy centres. In this study, we combine large-scale VLT… ▽ More Supermassive black holes (SMBHs) and nuclear star clusters (NSCs) co-exist in many galaxies. While the formation history of the black hole is essentially lost, NSCs preserve their evolutionary history imprinted onto their stellar populations and kinematics. Studying SMBHs and NSCs in tandem might help us to ultimately reveal the build-up of galaxy centres. In this study, we combine large-scale VLT/MUSE and high-resolution adaptive-optics-assisted VLT/SINFONI observations of the early-type galaxy FCC 47 with the goal being to assess the effect of a spatially (non-)variable initial mass function (IMF) on the determination of the mass of the putative SMBH in this galaxy. We achieve this by performing DYNAMITE Schwarzschild orbit-superposition modelling of the galaxy and its NSC. In order to properly take account of the stellar mass contribution to the galaxy potential, we create mass maps using a varying stellar mass-to-light ratio derived from single stellar population models with fixed and with spatially varying IMFs. Using the two mass maps, we estimate black hole masses of $(7.1^{+0.8}_{-1.1})\times 10^7\,M_{\odot}$ and $(4.4^{+1.2}_{-2.1}) \times 10^7\,M_{\odot}$ at $3σ$ signifance, respectively. Compared to models with constant stellar-mass-to-light ratio, the black hole masses decrease by 15% and 48%, respectively. Therefore, a varying IMF, both in its functional form and spatially across the galaxy, has a non-negligible effect on the SMBH mass estimate. Furthermore, we find that the SMBH in FCC 47 has probably not grown over-massive compared to its very over-massive NSC. △ Less

Submitted 26 April, 2023; originally announced April 2023.

Comments: 23 pages 19 Figures, accepted for publication in A&A

Journal ref: A&A 675, A18 (2023)

B-fields And dust in interstelLar fiLAments using Dust POLarization (BALLAD-POL): I. The massive filament G11.11-0.12 observed by SOFIA/HAWC+

Authors: Nguyen Bich Ngoc, Pham Ngoc Diep, Thiem Hoang, Le Ngoc Tram, Nguyen Chau Giang, Ngan Le, Thuong Duc Hoang, Nguyen Thi Phuong, Dieu D. Nguyen, Nguyen Minh Khang, Gia Bao Truong Le

Abstract: We report the first measurement of polarized thermal dust emission toward the entire early and massive Infrared Dark Cloud G11.11$-$0.12 taken by the polarimeter SOFIA/HAWC+ at 214 $μm$ wavelength. Magnetic fields (B-fields) obtained from the polarized emission tend to be perpendicular to the filament's spine. We produce a map of B-field strengths for the center region of the filament. The strengt… ▽ More We report the first measurement of polarized thermal dust emission toward the entire early and massive Infrared Dark Cloud G11.11$-$0.12 taken by the polarimeter SOFIA/HAWC+ at 214 $μm$ wavelength. Magnetic fields (B-fields) obtained from the polarized emission tend to be perpendicular to the filament's spine. We produce a map of B-field strengths for the center region of the filament. The strengths vary in the range of 100-600 $μ\rm{G}$ and are strongest along the filament's spine. The central region is sub-Alfvénic and mostly sub-critical meaning that B-fields dominate over turbulence and are strong enough to resist gravitational collapse. The alignment and properties of dust grains are studied in the filament using the RAdiative Torque (RAT) theory. We find the decrease of polarization degree $P$ with emission intensity $I$, i.e., depolarization effect, of the form $P\propto I^{-α}$ with $α\sim$0.8-0.9, implying a significant loss of grain alignment in the filament's spine. The depolarization can be explained by the decrease in RAT alignment efficiency toward the denser regions with lower dust temperature, and cannot be explained by the B-field tangling. We study the effect of the enhanced magnetic relaxation by embedded iron inclusions on RAT alignment and find that the high polarization fraction $P\sim$20-30\% in the outer layer of the filament is potential evidence for the enhanced RAT alignment by magnetic relaxation. This is the first time this effect is evaluated in a filament. Based on the polarization fraction and RAT alignment theory, we find evidence for grain growth in the filament. △ Less

Submitted 8 June, 2023; v1 submitted 21 February, 2023; originally announced February 2023.

Simulating supermassive black hole mass measurements for a sample of ultra massive galaxies using ELT/HARMONI high spatial resolution integral-field stellar kinematics

Authors: Dieu D. Nguyen, Michele Cappellari, Miguel Pereira-Santaella

Abstract: As the earliest relics of star formation episodes of the Universe, the most massive galaxies are the key to our understanding of the stellar population, cosmic structure, and SMBH evolution. However, the details of their formation histories remain uncertain. We address these problems by planning a large survey sample of 101 ultramassive galaxies ($z\le0.3$, $|δ+24^{\circ}|<45$°, $|b|>8$°), includi… ▽ More As the earliest relics of star formation episodes of the Universe, the most massive galaxies are the key to our understanding of the stellar population, cosmic structure, and SMBH evolution. However, the details of their formation histories remain uncertain. We address these problems by planning a large survey sample of 101 ultramassive galaxies ($z\le0.3$, $|δ+24^{\circ}|<45$°, $|b|>8$°), including 76\% ellipticals, 17\% lenticulars, and 7\% spirals brighter than $M_K\le-27$~mag (stellar mass $2\times10^{12}\lesssim M_\star\lesssim5\times10^{12}$~\Msun) with ELT/HARMONI. Our sample comprises diverse galaxy environments ranging from isolated to dense-cluster galaxies. The primary goals of the project are to (1) explore the stellar dynamics inside galaxy nuclei and weigh SMBHs, (2) constrain the black hole scaling relations at the highest mass, and (3) probe the late-time assembly of these most massive galaxies through the stellar population and kinematical gradients. We describe the survey, discuss the distinct demographics and environmental properties of the sample, and simulate their HARMONI $I_z$, $I_z+J$, and $H+K$-band observations via combining the inferred stellar-mass models from Pan-STARRS observations, an assumed synthetic spectrum of stars, and SMBHs with masses estimated based on different black hole scaling relations. Our simulations produce excellent state-of-art IFS and stellar kinematics ($ΔV_{\rm rms}\lesssim1.5$\%) in a relatively short exposure time. We use these stellar kinematics in combination with the JAM to reconstruct the SMBH mass and its error using the MCMC simulation. Thus, these simulations and modelings can be benchmarks to evaluate the instrument models and pipelines dedicated to HARMONI to exploit the unprecedented capabilities of ELT. △ Less

Submitted 17 September, 2023; v1 submitted 20 February, 2023; originally announced February 2023.

Comments: Accepted to MNRAS

Dynamics of hot galactic winds launched from spherically-stratified starburst cores

Authors: Dustin D. Nguyen, Todd A. Thompson, Evan E. Schneider, Sebastian Lopez, Laura A. Lopez

Abstract: The analytic galactic wind model derived by Chevalier and Clegg in 1985 (CC85) assumes $\textit{uniform}$ energy and mass-injection within the starburst galaxy nucleus. However, the structure of nuclear star clusters, bulges, and star-forming knots are non-uniform. We generalize to cases with spherically-symmetric energy/mass injection that scale as $r^{-Δ}$ within the starburst volume $R$, provid… ▽ More The analytic galactic wind model derived by Chevalier and Clegg in 1985 (CC85) assumes $\textit{uniform}$ energy and mass-injection within the starburst galaxy nucleus. However, the structure of nuclear star clusters, bulges, and star-forming knots are non-uniform. We generalize to cases with spherically-symmetric energy/mass injection that scale as $r^{-Δ}$ within the starburst volume $R$, providing solutions for $Δ= 0$, 1/2, 1, 3/2, and 2. In marked contrast with the CC85 model ($Δ=0$), which predicts zero velocity at the center, for a singular isothermal sphere profile ($Δ=2$), we find that the flow maintains a $\textit{constant}$ Mach number of $\mathcal{M}=\sqrt{3/5} \simeq 0.77$ throughout the volume. The fast interior flow can be written as $v_{r < R} = (\dot{E}_T/3\dot{M}_T)^{1/2} \simeq 0.41 \, v_\infty$, where $v_\infty$ is the asymptotic velocity, and $\dot{E}_T$ and $\dot{M}_T$ are the total energy and mass injection rates. For $v_\infty \simeq 2000 \, \mathrm{km \, s^{-1}}$, $v_{r<R} \simeq 820 \, \mathrm{km\, s^{-1}}$ throughout the wind-driving region. The temperature and density profiles of the non-uniform models may be important for interpreting spatially-resolved maps of starburst nuclei. We compute velocity resolved spectra to contrast the $Δ=0$ (CC85) and $Δ=2$ models. Next generation X-ray space telescopes such as XRISM may assess these kinematic predictions. △ Less

Submitted 8 November, 2022; v1 submitted 13 October, 2022; originally announced October 2022.

Comments: submitted (13 October 2022), accepted MNRAS Letters (08 November 2022)

X-ray and radio observations of central black holes in nearby low-mass early-type galaxies: Preliminary evidence for low Eddington fractions

Authors: R. Urquhart, L. I. McDermott, J. Strader, A. C. Seth, L. Chomiuk, N. Neumayer, D. D. Nguyen, E. Tremou

Abstract: We present new radio and X-ray observations of two nearby ($< 4$ Mpc) low-mass early-type galaxies with dynamically-confirmed central black holes: NGC 5102 and NGC 205. NGC 5102 shows a weak nuclear X-ray source and has no core radio emission. However, for the first time we demonstrate that it shows luminous extended radio continuum emission in low-resolution, low-frequency ($< 3$ GHz) data, consi… ▽ More We present new radio and X-ray observations of two nearby ($< 4$ Mpc) low-mass early-type galaxies with dynamically-confirmed central black holes: NGC 5102 and NGC 205. NGC 5102 shows a weak nuclear X-ray source and has no core radio emission. However, for the first time we demonstrate that it shows luminous extended radio continuum emission in low-resolution, low-frequency ($< 3$ GHz) data, consistent with jet lobes on scales $\gtrsim 100$ pc formed from past accretion and jet activity. By contrast, in new, extremely deep, strictly-simultaneous Very Large Array and Chandra observations, no radio or X-ray emission is detected from the black hole in NGC 205. We consider these measurements and upper limits in the context of the few other low-mass early-type galaxies with dynamically-confirmed black holes, and show that the mean ratio of bolometric to Eddington luminosity in this sample is only $\textrm{log} \, (L_{\rm bol}/L_{\rm Edd}) = -6.57\pm0.50$. These Eddington ratios are lower than typical in a comparison sample of more massive early-type galaxies, though this conclusion is quite tentative due to our small sample of low-mass galaxies and potential biases in the comparison sample. This preliminary result is in mild tension with previous work using less sensitive observations of more distant galaxies, which predict higher X-ray luminosities than we observe for low-mass galaxies. If it is confirmed that central black holes in low-mass galaxies typically have low Eddington ratios, this presents a challenge to measuring the occupation fraction of central black holes with standard optical emission line, X-ray, or radio surveys. △ Less

Submitted 29 September, 2022; originally announced September 2022.

Comments: 13 pages, 4 figures, 3 tables. Accepted for publication in ApJ

X-ray Properties of NGC 253's Starburst-Driven Outflow

Authors: Sebastian Lopez, Laura A. Lopez, Dustin D. Nguyen, Todd A. Thompson, Smita Mathur, Alberto D. Bolatto, Neven Vulic, Amy Sardone

Abstract: We analyze image and spectral data from $\approx$365~ks of observations from the {\it Chandra} X-ray Observatory of the nearby, edge-on starburst galaxy NGC 253 to constrain properties of the hot phase of the outflow. We focus our analysis on the $-$1.1 to $+$0.63 kpc region of the outflow and define several regions for spectral extraction where we determine best-fit temperatures and metal abundan… ▽ More We analyze image and spectral data from $\approx$365~ks of observations from the {\it Chandra} X-ray Observatory of the nearby, edge-on starburst galaxy NGC 253 to constrain properties of the hot phase of the outflow. We focus our analysis on the $-$1.1 to $+$0.63 kpc region of the outflow and define several regions for spectral extraction where we determine best-fit temperatures and metal abundances. We find that the temperatures and electron densities peak in the central $\sim$250 pc region of the outflow and decrease with distance. These temperature and density profiles are in disagreement with an adiabatic spherically expanding starburst wind model and suggest the presence of additional physics such as mass loading and non-spherical outflow geometry. Our derived temperatures and densities yield few-Myr cooling times in the nuclear region, which may imply that the hot gas can undergo bulk radiative cooling as it escapes along the minor axis. Our metal abundances of O, Ne, Mg, Si, S, and Fe all peak in the central region and decrease with distance along the outflow, with the exception of Ne which maintains a flat distribution. The metal abundances indicate significant dilution outside of the starburst region. We also find estimates on the mass outflow rates which are $2.8\:M_{\odot}/\rm{yr}$ in the northern outflow and $3.2\:M_{\odot}/\rm{yr}$ in the southern outflow. Additionally, we detect emission from charge exchange and find it has a significant contribution ($20-42$\%) to the total broad-band ($0.5-7$ keV) X-ray emission in the central and southern regions of the outflow. △ Less

Submitted 19 January, 2023; v1 submitted 19 September, 2022; originally announced September 2022.

Comments: 15 pages, 9 figures; submitted (19 September 2022), accepted ApJ (22 November 2022)

WISDOM Project XII. Clump properties and turbulence regulated by clump-clump collisions in the dwarf galaxy NGC404

Authors: Lijie Liu, Martin Bureau, Guang-Xing Li, Timothy A. Davis, Dieu D. Nguyen, Fu-Heng Liang, Woorak Choi, Mark R. Smith, Satoru Iguchi

Abstract: We present a study of molecular structures (clumps and clouds) in the dwarf galaxy NGC 404 using high-resolution (0.86x0.51 pc^2) Atacama Large Millimeter/sub-millimeter Array ^{12}CO(2-1) observations. We find two distinct regions in NGC 404: a gravitationally-stable central region (Toomre parameter Q=3-30) and a gravitationally-unstable molecular ring (Q<=1). The molecular structures in the cent… ▽ More We present a study of molecular structures (clumps and clouds) in the dwarf galaxy NGC 404 using high-resolution (0.86x0.51 pc^2) Atacama Large Millimeter/sub-millimeter Array ^{12}CO(2-1) observations. We find two distinct regions in NGC 404: a gravitationally-stable central region (Toomre parameter Q=3-30) and a gravitationally-unstable molecular ring (Q<=1). The molecular structures in the central region have a steeper size -- linewidth relation and larger virial parameters than those in the molecular ring, suggesting gas is more turbulent in the former. In the molecular ring, clumps exhibit a shallower mass -- size relation and larger virial parameters than clouds, implying density structures and dynamics are regulated by different physical mechanisms at different spatial scales. We construct an analytical model of clump-clump collisions to explain the results in the molecular ring. We propose that clump-clump collisions are driven by gravitational instabilities coupled with galactic shear, that lead to a population of clumps whose accumulation lengths (i.e. average separations) are approximately equal to their tidal radii. Our model-predicted clump masses and sizes (and mass -- size relation) and turbulence energy injection rates (and size -- linewidth relation) match the observations in the molecular ring very well, suggesting clump-clump collisions is the main mechanism regulating clump properties and gas turbulence in that region. As expected, our collision model does not apply to the central region, where turbulence is likely driven by clump migration. △ Less

Submitted 10 August, 2022; originally announced August 2022.

Comments: Accepted for publication on MNRAS

Galactic winds and bubbles from nuclear starburst rings

Authors: Dustin D. Nguyen, Todd A. Thompson

Abstract: Galactic outflows from local starburst galaxies typically exhibit a layered geometry, with cool $10^4\,$K flow sheathing a hotter $10^7\,$K, cylindrically-collimated, X-ray emitting plasma. Here, we argue that winds driven by energy-injection in a ring-like geometry can produce this distinctive large-scale multi-phase morphology. The ring configuration is motivated by the observation that massive… ▽ More Galactic outflows from local starburst galaxies typically exhibit a layered geometry, with cool $10^4\,$K flow sheathing a hotter $10^7\,$K, cylindrically-collimated, X-ray emitting plasma. Here, we argue that winds driven by energy-injection in a ring-like geometry can produce this distinctive large-scale multi-phase morphology. The ring configuration is motivated by the observation that massive young star clusters are often distributed in a ring at the host galaxy's inner Lindblad resonance, where larger-scale spiral arm structure terminates. We present parameterized three-dimensional radiative hydrodynamical simulations that follow the emergence and dynamics of energy-driven hot winds from starburst rings. In this Letter, we show that the flow shocks on itself within the inner ring hole, maintaining high $10^7$\,K temperatures, whilst flows that emerge from the wind-driving ring unobstructed can undergo rapid bulk cooling down to $10^4\,$K, producing a fast hot bi-conical outflow enclosed by a sheath of cooler nearly co-moving material without ram-pressure acceleration. The hot flow is collimated along the ring axis, even in the absence of pressure confinement from a galactic disk or magnetic fields. In the early stages of expansion, the emerging wind forms a bubble-like shape reminiscent of the Milky Way's eROSITA and Fermi bubbles and can reach velocities usually associated with AGN-driven winds. We discuss the physics of the ring configuration, the conditions for radiative bulk cooling, and the implications for future X-ray observations. △ Less

Submitted 4 August, 2022; v1 submitted 26 May, 2022; originally announced May 2022.

Comments: 19 pages, 8 figures, submitted (may 25, 2022), accepted ApJL (aug 03, 2022)

Extremely Buried Nucleus of IRAS 17208$-$0014 Observed at Sub-Millimeter and Near-Infrared Wavelengths

Authors: Shunsuke Baba, Masatoshi Imanishi, Takuma Izumi, Taiki Kawamuro, Dieu D. Nguyen, Takao Nakagawa, Naoki Isobe, Shusuke Onishi, Kosei Matsumoto

Abstract: The ultraluminous infrared galaxy IRAS 17208$-$0014 is a late-stage merger that hosts a buried active galactic nucleus (AGN). To investigate its nuclear structure, we performed high spatial resolution ($\sim0.\!\!^{\prime\prime}04\sim32\,\mathrm{pc}$) Atacama Large Millimeter/submillimeter Array (ALMA) observations in Band 9 ($\sim$450\,\micron\ or $\sim$660\,GHz), along with near-infrared AKARI s… ▽ More The ultraluminous infrared galaxy IRAS 17208$-$0014 is a late-stage merger that hosts a buried active galactic nucleus (AGN). To investigate its nuclear structure, we performed high spatial resolution ($\sim0.\!\!^{\prime\prime}04\sim32\,\mathrm{pc}$) Atacama Large Millimeter/submillimeter Array (ALMA) observations in Band 9 ($\sim$450\,\micron\ or $\sim$660\,GHz), along with near-infrared AKARI spectroscopy in 2.5--5.0\,\micron. The Band 9 dust continuum peaks at the AGN location, and toward this position CO($J$=6--5) and CS($J$=14--13) are detected in absorption. Comparison with non-local thermal equilibrium calculations indicates that, within the central beam ($r\sim20\,\mathrm{pc}$), there exists a concentrated component that is dense ($10^7\,\mathrm{cm}^{-2}$) and warm ($>$200\,K) and has a large column density ($N_\mathrm{H_2}>10^{23}\,\mathrm{cm}^{-2}$). The AKARI spectrum shows deep and broad CO ro-vibrational absorption at 4.67\,\micron. Its band profile is well reproduced with a similarly dense and large column but hotter ($\sim$1000\,K) gas. The region observed through absorption in the near-infrared is highly likely in the nuclear direction, as in the sub-millimeter, but with a narrower beam including a region closer to the nucleus. The central component is considered to possess a hot structure where vibrationally excited HCN emission originates. The most plausible heating source for the gas is X-rays from the AGN. The AKARI spectrum does not show other AGN signs in 2.5--4\,\micron, but this absence may be usual for AGNs buried in a hot mid-infrared core. Besides, based on our ALMA observations, we relate various nuclear structures of IRAS 17208$-$0014 that have been proposed in the literature. △ Less

Submitted 22 February, 2022; originally announced February 2022.

Comments: 28 pages, 16 figures, accepted for publication in ApJ

Cross-checking SMBH mass estimates in NGC 6958 -- I: Stellar dynamics from adaptive optics-assisted MUSE observations

Authors: Sabine Thater, Davor Krajnović, Peter M. Weilbacher, Dieu D. Nguyen, Martin Bureau, Michele Cappellari, Timothy A. Davis, Satoru Iguchi, Richard McDermid, Kyoko Onishi, Marc Sarzi, Glenn van de Ven

Abstract: Supermassive black hole masses (MBH) can dynamically be estimated with various methods and using different kinematic tracers. Different methods have only been cross-checked for a small number of galaxies and often show discrepancies. To understand these discrepancies, detailed cross-comparisons of additional galaxies are needed. We present the first part of our cross-comparison between stellar- an… ▽ More Supermassive black hole masses (MBH) can dynamically be estimated with various methods and using different kinematic tracers. Different methods have only been cross-checked for a small number of galaxies and often show discrepancies. To understand these discrepancies, detailed cross-comparisons of additional galaxies are needed. We present the first part of our cross-comparison between stellar- and gas-based MBH estimates in the nearby fast-rotating early-type galaxy NGC 6958. The measurements presented here are based on ground-layer adaptive optics-assisted Multi-Unit Spectroscopic Explorer (MUSE) science verification data at around 0.6 arcsec spatial resolution. The spatial resolution is a key ingredient for the measurement and we provide a Gaussian parametrisation of the adaptive optics-assisted point spread function (PSF) for various wavelengths. From the MUSE data, we extracted the stellar kinematics and constructed dynamical models. Using an axisymmetric Schwarzschild technique, we measured an MBH of (3.6+2.7-2.4)\times 10^8 Msun at 3σsignificance taking kinematical and dynamical systematics (e.g.,radially-varying mass-to-light ratio) into account. We also added a dark halo, but our data does not allow to constrain the dark matter fraction. Adding dark matter with an abundance matching prior results in a 25 per cent more massive black hole. Jeans anisotropic models return MBH of (4.6+2.5-2.7) \times 10^8 Msun and (8.6+0.8-0.8) \times 10^8 Msun at 3σconfidence for spherical and cylindrical alignment of the velocity ellipsoid, respectively. In a follow-up study, we will compare the stellar-based MBH with those from cold and warm gas tracers, which will provide additional constraints for the MBH for NGC 6958, and insights into assumptions that lead to potential systematic uncertainty. △ Less

Submitted 2 November, 2021; originally announced November 2021.

Comments: Accepted for publication in MNRAS, 21 pages, 18 Figures (Including Appendix). Supplementary information is added at the end of the file (3 pages)

Modeling extinction and reddening effects by circumstellar dust in the Betelgeuse envelope in the presence of radiative torque disruption

Authors: Bao Truong, Le Ngoc Tram, Thiem Hoang, Nguyen Chau Giang, Pham Ngoc Diep, Dieu D. Nguyen, Nguyen Thi Phuong, Thuong D. Hoang, Nguyen Bich Ngoc, Nguyen Fuda, Hien Phan, Tuan Van Bui

Abstract: Circumstellar dust is formed and evolved within the envelope of evolved stars, including Asymptotic Giant Branch (AGB) and Red Supergiant (RSG). The extinction of stellar light by circumstellar dust is vital for interpreting RSG/AGB observations and determining high-mass RSG progenitors of core-collapse supernovae. Nevertheless, circumstellar dust properties are not well understood. Modern underst… ▽ More Circumstellar dust is formed and evolved within the envelope of evolved stars, including Asymptotic Giant Branch (AGB) and Red Supergiant (RSG). The extinction of stellar light by circumstellar dust is vital for interpreting RSG/AGB observations and determining high-mass RSG progenitors of core-collapse supernovae. Nevertheless, circumstellar dust properties are not well understood. Modern understanding of dust evolution suggests that intense stellar radiation can radically change the dust properties across the circumstellar envelope through the RAdiative Torque Disruption (RAT-D) mechanism. In this paper, we study the impacts of RAT-D on the grain size distribution (GSD) of circumstellar dust and model its effects on photometric observations of $α$ Orionis (Betelgeuse). Due to the RAT-D effects, large grains formed in the dust formation zone are disrupted into smaller species of size $a < 0.5\,\rmμm$. Using the GSD constrained by the RAT-D effects, we model the visual extinction of background stars and Betelgeuse. We find that the extinction decreases at near-UV, optical, and infrared wavelengths while increasing at far-UV wavelengths. The resulting flux potentially reproduces the observation from the near-UV to near-IR range. Our results can be used to explain dust extinction and photometric observations toward other RSG/AGB stars. △ Less

Submitted 29 July, 2022; v1 submitted 22 October, 2021; originally announced October 2021.

Comments: 21 pages, 18 figures, 2 tables, accepted to ApJ

The MBHBM$^{\star}$ Project -- II. Molecular Gas Kinematics in the Lenticular Galaxy NGC 3593 Reveal a Supermassive Black Hole

Authors: Dieu D. Nguyen, Martin Bureau, Sabine Thater, Kristina Nyland, Mark den Brok, Michelle Cappellari, Timothy A. Davis, Jenny E. Greene, Nadine Neumayer, Masatoshi Imanishi, Takuma Izumi, Taiki Kawamuro, Shunsuke Baba, Phuong M. Nguyen, Satoru Iguchi, Takafumi Tsukui, Lam N. T., Than Ho

Abstract: As part of the Measuring Black Holes in Below Milky Way-mass (M$^\star$) galaxies (MBHBM$^\star$) Project, we present a dynamical measurement of the supermassive black hole (SMBH) mass in the nearby lenticular galaxy NGC 3593, using cold molecular gas $^{12}$CO(2-1) emission observed at an angular resolution of $\approx0''.3$ ($\approx10$ pc) with the Atacama Large Millimeter/submillimeter Array (… ▽ More As part of the Measuring Black Holes in Below Milky Way-mass (M$^\star$) galaxies (MBHBM$^\star$) Project, we present a dynamical measurement of the supermassive black hole (SMBH) mass in the nearby lenticular galaxy NGC 3593, using cold molecular gas $^{12}$CO(2-1) emission observed at an angular resolution of $\approx0''.3$ ($\approx10$ pc) with the Atacama Large Millimeter/submillimeter Array (ALMA). Our ALMA observations reveal a circumnuclear molecular gas disc (CND) elongated along the galaxy major axis and rotating around the SMBH. Using dynamical modelling, the molecular gas kinematics allow us to infer a SMBH mass $M_{\rm BH}=2.40_{-1.05}^{+1.87}\times10^6$ M$_\odot$ (only statistical uncertainties at the $3σ$ level). We also detect a massive core of cold molecular gas (CMC) of mass $M_{\rm CMC}=(5.4\pm1.2)\times10^6$ M$_\odot$ and effective (half-mass) radius $r_{\rm CMC,e}=11.2\pm2.8$ pc, co-spatial with a nuclear star cluster (NSC) of mass $M_{\rm NSC}=(1.67\pm0.48)\times10^7$ M$_\odot$ and effective radius $r_{\rm NSC,e}=5.0\pm1.0$~pc (or $0''.15\pm0''.03$). The mass profiles of the CMC and NSC are well described by Sérsic functions with indices $1-1.4$. Our $M_{\rm BH}$ and $M_{\rm NSC}$ estimates for NGC 3593 agree well with the recently compiled $M_{\rm BH}$-$M_{\rm NSC}$ scaling relation. Although the $M_{\rm NSC}$ uncertainty is twice the inferred $M_{\rm BH}$, the rapid central rise of the rotation velocities of the CND (as the radius decreases) clearly suggests a SMBH. Indeed, our dynamical models show that even if $M_{\rm NSC}$ is at the upper end of its allowed range, the evidence for a black hole does not vanish, but remains with a lower limit of $M_{\rm BH}>3\times10^5$ M$_\odot$. △ Less

Submitted 16 October, 2021; originally announced October 2021.

Comments: 27 pages, 18 Figures, 8 Tables. Accepted to MNRAS

Spatially Resolving the Star Formation Histories of Three Nearby Nuclear Star Clusters

Authors: Christian H. Hannah, Anil C. Seth, Dieu D. Nguyen, Antoine Dumont, Nikolay Kacharov, Nadine Neumayer, Mark den Brok

Abstract: The formation of nuclear star clusters (NSCs) remains an open question. In this work, we use spatially resolved HST/STIS spectroscopic observations of three nearby NSCs (hosted by NGC 5102, NGC 5206, and NGC 205) to constrain their formation histories by exploring radial variations of the stellar populations within each cluster. Utilizing full-spectrum fitting, we find substantial age and metallic… ▽ More The formation of nuclear star clusters (NSCs) remains an open question. In this work, we use spatially resolved HST/STIS spectroscopic observations of three nearby NSCs (hosted by NGC 5102, NGC 5206, and NGC 205) to constrain their formation histories by exploring radial variations of the stellar populations within each cluster. Utilizing full-spectrum fitting, we find substantial age and metallicity gradients within the central 0."9 (16~pc) of the NSC in NGC 5102 where populations near the center are young/metal-rich (age ~400~Myr and [M/H] ~-0.4) and become older/metal-poor at larger radii (mean age ~1~Gyr and mean [M/H] ~-1.6 in the radial range [0."3, 0."9]). This behavior suggests that the young/metal-rich population at the center was formed from a period of in situ formation, while the older/metal-poor populations were likely formed by inspiraled globular clusters. The two broad populations observed in the NGC 5102 NSC (young/metal-rich and old/metal-poor) appear to be linked to the transition between the two morphological components of the NSC derived from the surface-brightness profile in Nguyen et al. (2018)(arXiv:1711.04314). The radial ranges explored in NGC 5206 and NGC 205 were much smaller due to poor data quality; in NGC 5206 we find a similar metallicity gradient to NGC 5102 (but with much lower significance), while the data for NGC 205 is too poor to reach any conclusions. Overall, this data highlights the links between the morphological and stellar population complexity of NSCs and their formation mechanisms. △ Less

Submitted 24 September, 2021; originally announced September 2021.

Comments: 18 pages, 12 figures, to be published in The Astronomical Journal

Hard X-ray Irradiation Potentially Drives Negative AGN Feedback by Altering Molecular Gas Properties

Authors: Taiki Kawamuro, Claudio Ricci, Takuma Izumi, Masatoshi Imanishi, Shunsuke Baba, Dieu D. Nguyen, Kyoko Onishi

Abstract: To investigate the role of active galactic nucleus (AGN) X-ray irradiation on the interstellar medium (ISM), we systematically analyzed Chandra and ALMA CO($J$=2-1) data for 26 ultra-hard X-ray ($>$ 10 keV) selected AGNs at redshifts below 0.05. While Chandra unveils the distribution of X-ray-irradiated gas via Fe-K$α$ emission, the CO($J$=2-1) observations reveal that of cold molecular gas. At hi… ▽ More To investigate the role of active galactic nucleus (AGN) X-ray irradiation on the interstellar medium (ISM), we systematically analyzed Chandra and ALMA CO($J$=2-1) data for 26 ultra-hard X-ray ($>$ 10 keV) selected AGNs at redshifts below 0.05. While Chandra unveils the distribution of X-ray-irradiated gas via Fe-K$α$ emission, the CO($J$=2-1) observations reveal that of cold molecular gas. At high resolutions $\lesssim$ 1 arcsec, we derive Fe-K$α$ and CO($J$=2-1) maps for the nuclear 2 arcsec region, and for the external annular region of 2 arcsec-4 arcsec, where 2 arcsec is $\sim$ 100-600 pc for most of our AGNs. First, focusing on the external regions, we find the Fe-K$α$ emission for six AGNs above 2$σ$. Their large equivalent widths ($\gtrsim$ 1 keV) suggest a fluorescent process as their origin. Moreover, by comparing 6-7 keV/3-6 keV ratio, as a proxy of Fe-K$α$, and CO($J$=2-1) images for three AGNs with the highest significant Fe-K$α$ detections, we find a possible spatial separation. These suggest the presence of X-ray-irradiated ISM and the change in the ISM properties. Next, examining the nuclear regions, we find that (1) The 20-50 keV luminosity increases with the CO($J$=2-1) luminosity. (2) The ratio of CO($J$=2-1)-to-HCN($J$=1-0) luminosities increases with 20-50 keV luminosity, suggesting a decrease in the dense gas fraction with X-ray luminosity. (3) The Fe-K$α$-to-X-ray continuum luminosity ratio decreases with the molecular gas mass. This may be explained by a negative AGN feedback scenario: the mass accretion rate increases with gas mass, and simultaneously, the AGN evaporates a portion of the gas, which possibly affects star formation. △ Less

Submitted 20 September, 2021; originally announced September 2021.

Comments: 50 pages (27 pages for the main body), 26 figures, accepted for publication in ApJS

Studying magnetic fields and dust in M17 using polarized thermal dust emission observed by SOFIA/HAWC+

Authors: Thuong Duc Hoang, Nguyen Bich Ngoc, Pham Ngoc Diep, Le Ngoc Tram, Thiem Hoang, Wanggi Lim, Dieu D. Nguyen, Ngan Le, Nguyen Thi Phuong, Nguyen Fuda, Tuan Van Bui, Kate Pattle, Gia Bao Truong Le, Hien Phan, Nguyen Chau Giang

Abstract: We report the highest spatial resolution measurement of magnetic fields in M17 using thermal dust polarization taken by SOFIA/HAWC+ centered at 154 $μ$m wavelength. Using the Davis-Chandrasekhar-Fermi method, we found the presence of strong magnetic fields of $980 \pm 230\;μ$G and $1665 \pm 885\;μ$G in lower-density (M17-N) and higher-density (M17-S) regions, respectively. The magnetic field morph… ▽ More We report the highest spatial resolution measurement of magnetic fields in M17 using thermal dust polarization taken by SOFIA/HAWC+ centered at 154 $μ$m wavelength. Using the Davis-Chandrasekhar-Fermi method, we found the presence of strong magnetic fields of $980 \pm 230\;μ$G and $1665 \pm 885\;μ$G in lower-density (M17-N) and higher-density (M17-S) regions, respectively. The magnetic field morphology in M17-N possibly mimics the fields in gravitational collapse molecular cores while in M17-S the fields run perpendicular to the matter structure and display a pillar and an asymmetric hourglass shape. The mean values of the magnetic field strength are used to determine the Alfvénic Mach numbers ($\mathcal{M_A}$) of M17-N and M17-S which turn out to be sub-Alfvénic, or magnetic fields dominate turbulence. We calculate the mass-to-flux ratio, $λ$, and obtain $λ=0.07$ for M17-N and $0.28$ for M17-S. The sub-critical values of $λ$ are in agreement with the lack of massive stars formed in M17. To study dust physics, we analyze the relationship between the dust polarization fraction, $p$, and the thermal emission intensity, $I$, gas column density, $N({\rm H_2})$, and dust temperature, $T_{\rm d}$. The polarization fraction decreases with intensity as $I^{-α}$ with $α= 0.51$. The polarization fraction also decreases with increasing $N(\rm H_{2})$, which can be explained by the decrease of grain alignment by radiative torques (RATs) toward denser regions with a weaker radiation field and/or tangling of magnetic fields. The polarization fraction tends to increase with $T_{\rm d}$ first and then decreases when $T_ {\rm d} > 50$ K. The latter feature seen in the M17-N, where the gas density changes slowly with $T_{d}$, is consistent with the RAT disruption effect. △ Less

Submitted 12 November, 2021; v1 submitted 23 August, 2021; originally announced August 2021.

Comments: Submitted to Apj

Report number: 929:27 (21pp)

Journal ref: The Astrophysical Journal 2022 April 10

Mass-Loading and Non-Spherical Divergence in Hot Galactic Winds: Implications for X-ray Observations

Authors: Dustin D. Nguyen, Todd A. Thompson

Abstract: Cool clouds are expected to be destroyed and incorporated into hot supernova-driven galactic winds. The mass-loading of a wind by the cool medium modifies the bulk velocity, temperature, density, entropy, and abundance profiles of the hot phase relative to an un-mass-loaded outflow. We provide general equations and limits for this physics that can be used to infer the rate of cool gas entrainment… ▽ More Cool clouds are expected to be destroyed and incorporated into hot supernova-driven galactic winds. The mass-loading of a wind by the cool medium modifies the bulk velocity, temperature, density, entropy, and abundance profiles of the hot phase relative to an un-mass-loaded outflow. We provide general equations and limits for this physics that can be used to infer the rate of cool gas entrainment from X-ray observations, accounting for non-spherical expansion. In general, mass-loading flattens the density and temperature profiles, decreases the velocity and increases the entropy if the Mach number is above a critical value. We first apply this model to a recent high-resolution galactic outflow simulation where the mass-loading can be directly inferred. We show that the temperature, entropy, and composition profiles are well-matched, providing evidence that this physics sets the bulk hot gas profiles. We then model the diffuse X-ray emission from the local starburst M82. The non-spherical (more cylindrical) outflow geometry is directly taken from the observed X-ray surface brightness profile. These models imply a total mass-loading rate that is about equal to that injected in the starburst, $\simeq 10$ M$_\odot$ yr$^{-1}$, and they predict an asymptotic hot wind velocity of $\sim 1000\,{\rm km \ s^{-1}}$ that is $\sim1.5-2$ times smaller than previous predictions. We also show how the observed entropy profile can be used to constrain the outflow velocity, making predictions for future missions like XRISM. We argue that the observed X-ray limb-brightening may be explained by mass-loading at the outflow's edges. △ Less

Submitted 5 October, 2021; v1 submitted 16 August, 2021; originally announced August 2021.

Comments: 17 pages, 12 figures, submitted to MNRAS (Accepted Oct 2021)

Resolved nuclear kinematics link the formation and growth of nuclear star clusters with the evolution of their early and late-type hosts

Authors: Francesca Pinna, Nadine Neumayer, Anil Seth, Eric Emsellem, Dieu D. Nguyen, Torsten Boeker, Michele Cappellari, Richard M. McDermid, Karina Voggel, C. Jakob Walcher

Abstract: We present parsec-scale kinematics of eleven nearby galactic nuclei, derived from adaptive-optics assisted integral-field spectroscopy at (near-infrared) CO band-head wavelengths. We focus our analysis on the balance between ordered rotation and random motions, which can provide insights into the dominant formation mechanism of nuclear star clusters (NSCs). We divide our target sample into late- a… ▽ More We present parsec-scale kinematics of eleven nearby galactic nuclei, derived from adaptive-optics assisted integral-field spectroscopy at (near-infrared) CO band-head wavelengths. We focus our analysis on the balance between ordered rotation and random motions, which can provide insights into the dominant formation mechanism of nuclear star clusters (NSCs). We divide our target sample into late- and early-type galaxies, and discuss the nuclear kinematics of the two sub-samples, aiming at probing any link between NSC formation and host galaxy evolution. The results suggest that the dominant formation mechanism of NSCs is indeed affected by the different evolutionary paths of their hosts across the Hubble sequence. More specifically, nuclear regions in late-type galaxies are on average more rotation dominated, and the formation of nuclear stellar structures is potentially linked to the presence of gas funnelled to the center. Early-type galaxies, in contrast, tend to display slower-rotating NSCs with lower ellipticity. However, some exceptions suggest that in specific cases, early-type hosts can form NSCs in a way similar to spirals. △ Less

Submitted 26 October, 2021; v1 submitted 19 July, 2021; originally announced July 2021.

Comments: 46 pages, 21 figures, 3 tables, published in ApJ

Journal ref: ApJ 921 (2021) 8

Black hole mass measurement using ALMA observations of [CI] and CO emissions in the Seyfert 1 galaxy NGC7469

Authors: Dieu D. Nguyen, Takuma Izumi, Sabine Thater, Masatoshi Imanishi, Taiki Kawamuro, Shunsuke Baba, Suzuka Nakano, Jean L. Turner, Kotaro Kohno, Satoki Matsushita, Sergio Martin, David S. Meier, Phuong M. Nguyen, Lam T. Nguyen

Abstract: We present a supermassive black hole (SMBH) mass measurement in the Seyfert 1 galaxy NGC7469 using Atacama Large Millimeter/submillimeter Array (ALMA) observations of the atomic-${\rm [CI]}$(1-0) and molecular-$^{12}$CO(1-0) emission lines at the spatial resolution of $\approx0.3$" (or $\approx$ 100 pc). These emissions reveal that NGC7469 hosts a circumnuclear gas disc (CND) with a ring-like stru… ▽ More We present a supermassive black hole (SMBH) mass measurement in the Seyfert 1 galaxy NGC7469 using Atacama Large Millimeter/submillimeter Array (ALMA) observations of the atomic-${\rm [CI]}$(1-0) and molecular-$^{12}$CO(1-0) emission lines at the spatial resolution of $\approx0.3$" (or $\approx$ 100 pc). These emissions reveal that NGC7469 hosts a circumnuclear gas disc (CND) with a ring-like structure and a two-arm/bi-symmetric spiral pattern within it, surrounded by a starbursting ring. The CND has a relatively low $σ/V\approx0.35$ ($r\sim0.5$") and $\sim0.19$ ($r>0.5"$), suggesting that the gas is dynamically settled and suitable for dynamically deriving the mass of its central source. As is expected from X-ray dominated region (XDR) effects that dramatically increase an atomic carbon abundance by dissociating CO molecules, we suggest that the atomic [CI](1-0) emission is a better probe of SMBH masses than CO emission in AGNs. Our dynamical model using the ${\rm [CI]}$(1-0) kinematics yields a $M_{\rm BH}=1.78^{+2.69}_{-1.10}\times10^7$M$_\odot$ and $M/L_{\rm F547M}=2.25^{+0.40}_{-0.43}$ (M$_\odot$/L$_\odot$). The model using the CO(1-0) kinematics also gives a consistent $M_{\rm BH}$ with a larger uncertainty, up to an order of magnitude, i.e.\ $M_{\rm BH}=1.60^{+11.52}_{-1.45}\times10^7$M$_\odot$. This newly dynamical $M_{\rm BH}$ is $\approx$ 2 times higher than the mass determined from the reverberation mapped (RM) method using emissions arising in the unresolved broad-line region (BLR). Given this new $M_{\rm BH}$, we are able to constrain the specific RM dimensionless scaling factor of $f=7.2^{+4.2}_{-3.4}$ for the AGN BLR in NGC7469. The gas within the unresolved BLR thus has a Keplerian virial velocity component and the inclination of $i\approx11.0^\circ$$_{-2.5}^{+2.2}$, confirming its face-on orientation in a Seyfert 1 AGN by assuming a geometrically thin BLR model. △ Less

Submitted 8 April, 2021; originally announced April 2021.

Comments: 22 pages, 16 figures, 7 tables. Accepted for publication on MNRAS

ALMA 0.02"-resolution observations reveal HCN-abundance-enhanced counter-rotating and outflowing dense molecular gas at the NGC 1068 nucleus

Authors: Masatoshi Imanishi, Dieu D. Nguyen, Keiichi Wada, Yoshiaki Hagiwara, Satoru Iguchi, Takuma Izumi, Nozomu Kawakatu, Kouichiro Nakanishi, Kyoko Onishi

Abstract: We present ALMA ~0.02"-resolution observations of the nucleus of the nearby (~14 Mpc) type-2 AGN NGC 1068 at HCN/HCO+/HNC J=3-2 lines, as well as at their 13C isotopologue and vibrationally excited lines, to scrutinize the morphological/dynamical/chemical/physical properties of dense molecular gas in the putative dusty molecular torus around a mass-accreting supermassive black hole. We confirm alm… ▽ More We present ALMA ~0.02"-resolution observations of the nucleus of the nearby (~14 Mpc) type-2 AGN NGC 1068 at HCN/HCO+/HNC J=3-2 lines, as well as at their 13C isotopologue and vibrationally excited lines, to scrutinize the morphological/dynamical/chemical/physical properties of dense molecular gas in the putative dusty molecular torus around a mass-accreting supermassive black hole. We confirm almost east-west-oriented dense molecular gas emission both morphologically and dynamically, which we regard as coming from the torus. Bright emission is compact (<3 pc), and low-surface-brightness emission extends out to 5-7 pc. These dense molecular gas properties are not symmetric between the eastern and western torus. The HCN J=3-2 emission is stronger than the HCO+ J=3-2 emission within the ~7 pc torus region, with an estimated dense molecular mass of (0.4-1.0)x10^6Msun. We interpret that HCN abundance is enhanced in the torus. We detect signatures of outflowing dense molecular gas and a vibrationally excited HCN J=3-2 line. Finally, we find that in the innermost (<1 pc) part of the torus, the dense molecular line rotation velocity, relative to the systemic velocity, is the opposite of that in the outer (>2 pc) part, in both the eastern and western torus. We prefer a scenario of counter-rotating dense molecular gas with innermost almost-Keplerian-rotation and outer slowly rotating (far below Keplerian) components. Our high-spatial-resolution dense molecular line data reveal that torus properties of NGC 1068 are much more complicated than the simple axi-symmetrically rotating torus picture in the classical AGN unification paradigm. △ Less

Submitted 18 August, 2020; originally announced August 2020.

Comments: 45 pages, 20 Figures. Accepted for publication in ApJ

Revealing the intermediate-mass black hole at the heart of the dwarf galaxy NGC 404 with sub-parsec resolution ALMA observations

Authors: Timothy A. Davis, Dieu D. Nguyen, Anil C. Seth, Jenny E. Greene, Kristina Nyland, Aaron J. Barth, Martin Bureau, Michele Cappellari, Mark den Brok, Satoru Iguchi, Federico Lelli, Lijie Liu, Nadine Neumayer, Eve V. North, Kyoko Onishi, Marc Sarzi, Mark D. Smith, Thomas G. Williams

Abstract: We estimate the mass of the intermediate-mass black hole at the heart of the dwarf elliptical galaxy NGC 404 using Atacama Large Millimeter/submillimeter Array (ALMA) observations of the molecular interstellar medium at an unprecedented linear resolution of ~0.5 pc, in combination with existing stellar kinematic information. These ALMA observations reveal a central disc/torus of molecular gas clea… ▽ More We estimate the mass of the intermediate-mass black hole at the heart of the dwarf elliptical galaxy NGC 404 using Atacama Large Millimeter/submillimeter Array (ALMA) observations of the molecular interstellar medium at an unprecedented linear resolution of ~0.5 pc, in combination with existing stellar kinematic information. These ALMA observations reveal a central disc/torus of molecular gas clearly rotating around the black hole. This disc is surrounded by a morphologically and kinematically complex flocculent distribution of molecular clouds, that we resolve in detail. Continuum emission is detected from the central parts of NGC 404, likely arising from the Rayleigh-Jeans tail of emission from dust around the nucleus, and potentially from dusty massive star-forming clumps at discrete locations in the disc. Several dynamical measurements of the black hole mass in this system have been made in the past, but they do not agree. We show here that both the observed molecular gas and stellar kinematics independently require a ~5x10$^5$ Msun black hole once we include the contribution of the molecular gas to the potential. Our best estimate comes from the high-resolution molecular gas kinematics, suggesting the black hole mass of this system is 5.5$^{+4.1}_{-3.8}\times$10$^5$ Msun (at the 99% confidence level), in good agreement with our revised stellar kinematic measurement and broadly consistent with extrapolations from the black hole mass - velocity dispersion and black hole mass - bulge mass relations. This highlights the need to accurately determine the mass and distribution of each dynamically important component around intermediate-mass black holes when attempting to estimate their masses. △ Less

Submitted 10 July, 2020; originally announced July 2020.

Comments: 18 pages, 15 figures, MNRAS accepted

ALMA Observations of Multiple-CO and C Lines Toward the Active Galactic Nucleus of NGC 7469: X-Ray-dominated Region Caught in the Act

Authors: Takuma Izumi, Dieu D. Nguyen, Masatoshi Imanishi, Taiki Kawamuro, Shunsuke Baba, Suzuka Nakano, Kotaro Kohno, Satoki Matsushita, David S. Meier, Jean L. Turner, Tomonari Michiyama, Nanase Harada, Sergio Martín, Kouichiro Nakanishi, Shuro Takano, Tommy Wiklind, Naomasa Nakai, Pei-Ying Hsieh

Abstract: We used the Atacama Large Millimeter/submillimeter Array (ALMA) to map $^{12}$CO($J$ = 1-0), $^{12}$CO($J$ = 2-1), $^{12}$CO($J$ = 3-2), $^{13}$CO($J$ = 2-1), and [CI]($^3P_1$-$^3P_0$) emission lines around the type 1 active galactic nucleus (AGN) of NGC 7469 ($z = 0.0164$) at $\sim 100$ pc resolutions. The CO lines are bright both in the circumnuclear disk (central $\sim 300$ pc) and the surround… ▽ More We used the Atacama Large Millimeter/submillimeter Array (ALMA) to map $^{12}$CO($J$ = 1-0), $^{12}$CO($J$ = 2-1), $^{12}$CO($J$ = 3-2), $^{13}$CO($J$ = 2-1), and [CI]($^3P_1$-$^3P_0$) emission lines around the type 1 active galactic nucleus (AGN) of NGC 7469 ($z = 0.0164$) at $\sim 100$ pc resolutions. The CO lines are bright both in the circumnuclear disk (central $\sim 300$ pc) and the surrounding starburst (SB) ring ($\sim 1$ kpc diameter), with two bright peaks on either side of the AGN. By contrast, the [CI]($^3P_1$-$^3P_0$) line is strongly peaked on the AGN. Consequently, the brightness temperature ratio of [CI]($^3P_1$-$^3P_0$) to $^{13}$CO(2-1) is $\sim 20$ at the AGN, as compared to $\sim 2$ in the SB ring. Our local thermodynamic equilibrium (LTE) and non-LTE models indicate that the enhanced line ratios (or CI enhancement) are due to an elevated C$^0$/CO abundance ratio ($\sim 3-10$) and temperature ($\sim 100-500$ K) around the AGN as compared to the SB ring (abundance ratio $\sim 1$, temperature $\lesssim 100$ K), which accords with the picture of the X-ray-dominated Region (XDR). Based on dynamical modelings, we also provide CO(1-0)-to- and [CI]($^3P_1$-$^3P_0$)-to-molecular mass conversion factors at the central $\sim 100$ pc of this AGN as $α_{\rm CO} = 4.1$ and $α_{\rm CI} = 4.4~M_\odot$ (K km s$^{-1}$ pc$^2$)$^{-1}$, respectively. Our results suggest that the CI enhancement is potentially a good marker of AGNs that could be used in a new submillimeter diagnostic method toward dusty environments. △ Less

Submitted 3 July, 2020; v1 submitted 16 June, 2020; originally announced June 2020.

Comments: 25 pages, 12 figures, 8 tables. Accepted for publication in ApJ. Minor updates in a replacement

Temperature and Metallicity Gradients in the Hot Gas Outflows of M82

Authors: Laura A. Lopez, Smita Mathur, Dustin D. Nguyen, Todd A. Thompson, Grace M. Olivier

Abstract: We utilize deep Chandra X-ray Observatory imaging and spectra of M82, the prototype of a starbursting galaxy with a multiphase wind, to map the hot plasma properties along the minor axis of the galaxy. We extract spectra from 11 regions up to 2.5 kpc from the starbursting midplane and model the data as a multi-temperature, optically thin thermal plasma with contributions from a non-thermal (power-… ▽ More We utilize deep Chandra X-ray Observatory imaging and spectra of M82, the prototype of a starbursting galaxy with a multiphase wind, to map the hot plasma properties along the minor axis of the galaxy. We extract spectra from 11 regions up to 2.5 kpc from the starbursting midplane and model the data as a multi-temperature, optically thin thermal plasma with contributions from a non-thermal (power-law) component and from charge exchange (CX). We examine the gradients in best-fit parameters, including the intrinsic column density, plasma temperature, metal abundances, and number density of the hot gas as a function of distance from the M82 nucleus. We find that the temperatures and number densities of the warm-hot and hot plasma peak at the starbursting ridge and decreases along the minor axis. The temperature and density profiles are inconsistent with spherical adiabatic expansion of a super-heated wind and suggest mass loading and mixing of the hot phase with colder material. Non-thermal emission is detected in all of the regions considered, and CX comprises 8-25% of the total absorption-corrected, broad-band (0.5-7 keV) X-ray flux. We show that the abundances of O, Ne, Mg, and Fe are roughly constant across the regions considered, while Si and S peak within 500 pc of the central starburst. These findings support a direct connection between the M82 superwind and the warm-hot, metal-rich circumgalactic medium (CGM). △ Less

Submitted 10 October, 2020; v1 submitted 15 June, 2020; originally announced June 2020.

Comments: 15 pages, 8 figures, ApJ in press

AGN X-ray irradiation of CO gas in NGC 2110 revealed by $Chandra$ and ALMA

Authors: Taiki Kawamuro, Takuma Izumi, Kyoko Onishi, Masatoshi Imanishi, Dieu D. Nguyen, Shunsuke Baba

Abstract: We report spatial distributions of the Fe-K$α$ line at 6.4 keV and the CO($J$ = 2--1) line at 230.538 GHz in NGC 2110, which are respectively revealed by $Chandra$ and ALMA at $\approx$ 0.5 arcsec. A $Chandra$ 6.2--6.5 keV-to-3.0--6.0 keV image suggests that the Fe-K$α$ emission extends preferentially in a northwest-to-southeast direction out to $\sim$ 3 arcsec, or 500 pc, on each side. Spatially-… ▽ More We report spatial distributions of the Fe-K$α$ line at 6.4 keV and the CO($J$ = 2--1) line at 230.538 GHz in NGC 2110, which are respectively revealed by $Chandra$ and ALMA at $\approx$ 0.5 arcsec. A $Chandra$ 6.2--6.5 keV-to-3.0--6.0 keV image suggests that the Fe-K$α$ emission extends preferentially in a northwest-to-southeast direction out to $\sim$ 3 arcsec, or 500 pc, on each side. Spatially-resolved spectral analyses support this by finding significant Fe-K$α$ emission lines only in northwest and southeast regions. Moreover, their equivalent widths are found $\sim$ 1.5 keV, indicative for the fluorescence by nuclear X-ray irradiation as the physical origin. By contrast, CO($J$ = 2--1) emission is weak therein. For quantitative discussion, we derive ionization parameters by following an X-ray dominated region (XDR) model. We then find them high enough to interpret the weakness as the result of X-ray dissociation of CO and/or H$_2$. Another possibility also remains that CO molecules follow a super-thermal distribution, resulting in brighter emission in higher-$J$ lines. Further follow-up observations are encouraged to draw a conclusion on what predominantly changes the inter-stellar matter properties, and whether the X-ray irradiation eventually affects the surrounding star formation as an AGN feedback. △ Less

Submitted 20 April, 2020; originally announced April 2020.

Comments: 15 pages, 11 figures, 3 tables. Accepted for publication in ApJ

Testing the robustness of black hole mass measurements with ALMA and MUSE

Authors: Sabine Thater, Davor Krajnović, Dieu D. Nguyen, Satoru Iguchi, Peter M. Weilbacher

Abstract: We present our ongoing work of using two independent tracers to estimate the supermassive black hole mass in the nearby early-type galaxy NGC 6958; namely integrated stellar and molecular gas kinematics. We used data from the Atacama Large Millimeter/submillimeter Array (ALMA), and the adaptive-optics assisted Multi-Unit Spectroscopic Explorer (MUSE) and constructed state-of-the-art dynamical mode… ▽ More We present our ongoing work of using two independent tracers to estimate the supermassive black hole mass in the nearby early-type galaxy NGC 6958; namely integrated stellar and molecular gas kinematics. We used data from the Atacama Large Millimeter/submillimeter Array (ALMA), and the adaptive-optics assisted Multi-Unit Spectroscopic Explorer (MUSE) and constructed state-of-the-art dynamical models. The different methods provide black hole masses of $(2.89\pm 2.05) \times 10^8M_{\odot}$ from stellar kinematics and $(1.35\pm 0.09) \times 10^8M_{\odot}$ from molecular gas kinematics which are consistent within their $3σ$ uncertainties. Compared to recent M$_{\rm BH}$ - $σ_{\rm e}$ scaling relations, we derive a slightly over-massive black hole. Our results also confirm previous findings that gas-based methods tend to provide lower black hole masses than stellar-based methods. More black hole mass measurements and an extensive analysis of the method-dependent systematics are needed in the future to understand this noticeable discrepancy. △ Less

Submitted 26 November, 2019; originally announced November 2019.

Comments: 5 pages, 3 figures, Proceeding of IAU Symposium 353, Galactic Dynamics in the Era of Large Surveys, ed. M. Valluri & J. A. Sellwood, Cambridge Univ. Press, in press; submitted September, 2019

The MBHBM* Project -- I: Measurement of the Central Black Hole Mass in Dwarf Galaxy NGC 3504 Using Molecular Gas Kinematics

Authors: Dieu D. Nguyen, Mark den Brok, Anil C. Seth, Satoru Iguchi, Jenny E. Greene, Timothy Davis, Masatoshi Imanishi, Takuma, Izumi, Michelle Cappellari, Nadine Neumayer, Kristina Nyland, Takafumi Tsukui, Kouichiro Nakanishi, Phuong M. Nguyen, Quang L. Nguyen, Sabine Thater, Martin Bureau, Kyoko Onishi, Karina T. Voggel, Ngan M. Le, Trung V. Dinh

Abstract: We present the first measurement of the mass of a supermassive black hole (SMBH) in the nearby double-barred spiral galaxy NGC 3504 as part of the Measuring Black Holes Below the Milky Way ($M_{\star}$) mass galaxies (MBHBM$_{\star}$) Project. Our analysis is based on Atacama Large Millimeter/submillimeter Array (ALMA) Cycle-5 observations of the ${\rm ^{12}CO(2-1)}$ emission line. NGC 3504 has a… ▽ More We present the first measurement of the mass of a supermassive black hole (SMBH) in the nearby double-barred spiral galaxy NGC 3504 as part of the Measuring Black Holes Below the Milky Way ($M_{\star}$) mass galaxies (MBHBM$_{\star}$) Project. Our analysis is based on Atacama Large Millimeter/submillimeter Array (ALMA) Cycle-5 observations of the ${\rm ^{12}CO(2-1)}$ emission line. NGC 3504 has a circumnuclear gas disk (CND), which has a relatively high-velocity dispersion of 30 \kms. Our dynamical models of the CND yield a \Mbh~of $M_{\rm BH}=1.02^{+0.18}_{-0.15}\times10^7$\Msun and a mass-to-light ratio in $H$-band of \ml$_{\rm H}=0.66^{+1.44}_{-0.65}$ (\Msun/\Lsun). This black hole (BH) mass is consistent with BH--galaxy scaling relations. We also detect a central deficiency in the ${\rm ^{12}CO(2-1)}$ integrated intensity map with a diameter of 2.7 pc at the putative position of the SMBH. However, this hole is filled by a dense gas tracer ${\rm CS(5-4)}$ that peaks at the galaxy center found in one of the three low-velocity-resolution continuum spectral correlators. The ${\rm CS(5-4)}$ line has the same kinematics with the ${\rm ^{12}CO(2-1)}$ line within the CND, suggesting that it is also an alternative transition for measuring the central \Mbh~in NGC 3504 probably more accurately than the current commonly used of ${\rm ^{12}CO(2-1)}$ due to its centralization. △ Less

Submitted 11 February, 2019; originally announced February 2019.

Comments: 25 pages, 15 Figures, 6 Tables, Submitted to ApJ

Improved Dynamical Constraints on the Masses of the Central Black Holes in Nearby Low-mass Early-type Galactic Nuclei And the First Black Hole Determination for NGC 205

Authors: Dieu D. Nguyen, Anil C. Seth, Nadine Neumayer, Satoru Iguchi, Michele Cappellari, Jay Strader, Laura Chomiuk, Evangelia Tremou, Fabio Pacucci, Kouichiro Nakanishi, Arash Bahramian, Phuong M. Nguyen, Mark den Brok, Christopher Ahn, Karina T. Voggel, Nikolay Kacharov, Takafumi Tsukui, Cuc K. Ly, Antoine Dumont, Renuka Pechetti

Abstract: We improve the dynamical black hole (BH) mass estimates in three nearby low-mass early-type galaxies--NGC 205, NGC 5102, and NGC 5206. We use new \hst/STIS spectroscopy to fit the star formation histories of the nuclei in these galaxies, and use these measurements to create local color--mass-to-light ratio (\ml) relations. We then create new mass models from \hst~imaging and combined with adaptive… ▽ More We improve the dynamical black hole (BH) mass estimates in three nearby low-mass early-type galaxies--NGC 205, NGC 5102, and NGC 5206. We use new \hst/STIS spectroscopy to fit the star formation histories of the nuclei in these galaxies, and use these measurements to create local color--mass-to-light ratio (\ml) relations. We then create new mass models from \hst~imaging and combined with adaptive optics kinematics, we use Jeans dynamical models to constrain their BH masses. The masses of the central BHs in NGC 5102 and NGC 5206 are both below one million solar masses and are consistent with our previous estimates, $9.12_{-1.53}^{+1.84}\times10^5$\Msun~and $6.31_{-2.74}^{+1.06}\times10^5$\Msun~(3$σ$ errors), respectively. However, for NGC 205, the improved models suggest the presence of a BH for the first time, with a best-fit mass of $6.8_{-6.7}^{+95.6}\times10^3$\Msun~(3$σ$ errors). This is the least massive central BH mass in a galaxy detected using any method. We discuss the possible systematic errors of this measurement in detail. Using this BH mass, the existing upper limits of both X-ray, and radio emissions in the nucleus of NGC 205 suggest an accretion rate $\lesssim$$10^{-5}$ of the Eddington rate. We also discuss the color--\mleff~relations in our nuclei and find that the slopes of these vary significantly between nuclei. Nuclei with significant young stellar populations have steeper color--\mleff~relations than some previously published galaxy color--\mleff~relations. △ Less

Submitted 14 January, 2019; originally announced January 2019.

Comments: 31 pages, 19 figures, 6 tables, Accepted to ApJ

The Black Hole in the Most Massive Ultracompact Dwarf Galaxy M59-UCD3

Authors: Christopher P. Ahn, Anil C. Seth, Michele Cappellari, Davor Krajnović, Jay Strader, Karina T. Voggel, Jonelle L. Walsh, Arash Bahramian, Holger Baumgardt, Jean Brodie, Igor Chilingarian, Laura Chomiuk, Mark den Brok, Matthias Frank, Michael Hilker, Richard M. McDermid, Steffen Mieske, Nadine Neumayer, Dieu D. Nguyen, Renuka Pechetti, Aaron J. Romanowsky, Lee Spitler

Abstract: We examine the internal properties of the most massive ultracompact dwarf galaxy (UCD), M59-UCD3, by combining adaptive optics assisted near-IR integral field spectroscopy from Gemini/NIFS, and Hubble Space Telescope (HST) imaging. We use the multi-band HST imaging to create a mass model that suggests and accounts for the presence of multiple stellar populations and structural components. We combi… ▽ More We examine the internal properties of the most massive ultracompact dwarf galaxy (UCD), M59-UCD3, by combining adaptive optics assisted near-IR integral field spectroscopy from Gemini/NIFS, and Hubble Space Telescope (HST) imaging. We use the multi-band HST imaging to create a mass model that suggests and accounts for the presence of multiple stellar populations and structural components. We combine these mass models with kinematics measurements from Gemini/NIFS to find a best-fit stellar mass-to-light ratio ($M/L$) and black hole (BH) mass using Jeans Anisotropic Models (JAM), axisymmetric Schwarzschild models, and triaxial Schwarzschild models. The best fit parameters in the JAM and axisymmetric Schwarzschild models have black holes between 2.5 and 5.9 million solar masses. The triaxial Schwarzschild models point toward a similar BH mass, but show a minimum $χ^2$ at a BH mass of $\sim 0$. Models with a BH in all three techniques provide better fits to the central $V_{rms}$ profiles, and thus we estimate the BH mass to be $4.2^{+2.1}_{-1.7} \times 10^{6}$ M$_\odot$ (estimated 1$σ$ uncertainties). We also present deep radio imaging of M59-UCD3 and two other UCDs in Virgo with dynamical BH mass measurements, and compare these to X-ray measurements to check for consistency with the fundamental plane of BH accretion. We detect faint radio emission in M59cO, but find only upper limits for M60-UCD1 and M59-UCD3 despite X-ray detections in both these sources. The BH mass and nuclear light profile of M59-UCD3 suggests it is the tidally stripped remnant of a $\sim$10$^{9-10}$ M$_\odot$ galaxy. △ Less

Submitted 6 April, 2018; originally announced April 2018.

Comments: 17 pages, 14 figures, 5 tables

Upper limits on the presence of central massive black holes in two ultra-compact dwarf galaxies in Centaurus A

Authors: Karina T. Voggel, Anil C. Seth, Nadine Neumayer, Steffen Mieske, Igor Chilingarian, Christopher Ahn, Holger Baumgardt, Michael Hilker, Dieu D. Nguyen, Aaron J. Romanowsky, Jonelle L. Walsh, Mark den Brok, Jay Strader

Abstract: The recent discovery of massive black holes (BHs) in the centers of high-mass ultra compact dwarf galaxies (UCDs) suggests that at least some are the stripped nuclear star clusters of dwarf galaxies. We present the first study that investigates whether such massive BHs, and therefore stripped nuclei, also exist in low-mass ($M<10^{7}M_{\odot}$) UCDs. We constrain the BH masses of two UCDs located… ▽ More The recent discovery of massive black holes (BHs) in the centers of high-mass ultra compact dwarf galaxies (UCDs) suggests that at least some are the stripped nuclear star clusters of dwarf galaxies. We present the first study that investigates whether such massive BHs, and therefore stripped nuclei, also exist in low-mass ($M<10^{7}M_{\odot}$) UCDs. We constrain the BH masses of two UCDs located in Centaurus A (UCD320 and UCD330) using Jeans modeling of the resolved stellar kinematics from adaptive optics VLT/SINFONI data. No massive BHs are found in either UCD. We find a $3\,σ$ upper limit on the central BH mass in UCD\,330 of $M_{\bullet}<1.0\times10^{5}M_{\odot}$, which corresponds to 1.7\% of the total mass. This excludes a high mass fraction BH and would only allow a low-mass BHs similar to those claimed to be detected in Local Group GCs. For UCD320, poorer data quality results in a less constraining $3\,σ$ upper limit of $M_{\bullet}<1\times10^{6}M_{\odot}$, which is equal to 37.7\% of the total mass. The dynamical $M/L$ of UCD320 and UCD330 are not inflated compared to predictions from stellar population models. The non-detection of BHs in these low-mass UCDs is consistent with the idea that elevated dynamical $M/L$s do indicate the presence of a substantial BH. Despite not detecting massive BHs, these systems could still be stripped nuclei. The strong rotation ($v/σ$ of 0.3 to 0.4) in both UCDs and the two-component light profile in UCD330 support the idea that these UCDs may be stripped nuclei of low-mass galaxies where the BH occupation fraction is not yet known. △ Less

Submitted 26 March, 2018; originally announced March 2018.

Comments: Accepted for publication in ApJ, 21 pages, 10 Figures

Improved dynamical constraints on the mass of the central black hole in NGC 404

Authors: Dieu D. Nguyen

Abstract: We determine the dynamical black hole mass in NGC404 including modeling of the nuclear stellar populations. We combine \hst/STIS spectroscopy with WFC3 images to create a local color-\ml~relation derived from stellar population modeling of the STIS data. We then use this to create a mass model for the nuclear region. We use Jeans modeling to fit this mass model to adaptive optics stellar kinematic… ▽ More We determine the dynamical black hole mass in NGC404 including modeling of the nuclear stellar populations. We combine \hst/STIS spectroscopy with WFC3 images to create a local color-\ml~relation derived from stellar population modeling of the STIS data. We then use this to create a mass model for the nuclear region. We use Jeans modeling to fit this mass model to adaptive optics stellar kinematic observations from Gemini/NIFS. From our stellar dynamical modeling, we find a 3$σ$ upper limit on the black hole mass of $1.5\times10^5$\Msun. △ Less

Submitted 6 December, 2017; originally announced December 2017.

Comments: Proceedings of the Star Formation in Different Environments, ICISE, Quy Nhon, Vietnam, 2016, (eds. D. Johnstone, T. Hoang, F. Nakamura, Q. Nguyen-Luong, and J. T. Tranh Van)

Report number: SFDE/2016/68

Nearby Early-Type Galactic Nuclei at High Resolution: Dynamical Black Hole and Nuclear Star Cluster Mass Measurements

Authors: Dieu D. Nguyen, Anil C. Seth, Nadine Neumayer, Sebastian Kamann, Karina T. Voggel, Michele Cappellari, Arianna Picotti, Phuong M. Nguyen, Torsten Böker, Victor Debattista, Nelson Caldwell, Richard McDermid, Bastian Nathan, Christopher E. Ahn, Renuka Pechetti

Abstract: We present a detailed study of the nuclear star clusters (NSCs) and massive black holes (BHs) of four of the nearest low-mass early-type galaxies: M32, NGC205, NGC5012, and NGC5206. We measure dynamical masses of both the BHs and NSCs in these galaxies using Gemini/NIFS or VLT/SINFONI stellar kinematics, Hubble Space Telescope (HST) imaging, and Jeans Anisotropic Models. We detect massive BHs in M… ▽ More We present a detailed study of the nuclear star clusters (NSCs) and massive black holes (BHs) of four of the nearest low-mass early-type galaxies: M32, NGC205, NGC5012, and NGC5206. We measure dynamical masses of both the BHs and NSCs in these galaxies using Gemini/NIFS or VLT/SINFONI stellar kinematics, Hubble Space Telescope (HST) imaging, and Jeans Anisotropic Models. We detect massive BHs in M32, NGC5102, and NGC5206, while in NGC205, we find only an upper limit. These BH mass estimates are consistent with previous measurements in M32 and NGC205, while those in NGC5102 and NGC5206 are estimated for the first time, and both found to be $<$$10^6~M_{\odot}$. This adds to just a handful of galaxies with dynamically measured sub-million $M_{\odot}$ central BHs. Combining these BH detections with our recent work on NGC404's BH, we find that 80\% (4/5) of nearby, low-mass ($10^9-10^{10}~M_{\odot}$; $σ_{\star}\sim20-70$ km/s) early-type galaxies host BHs. Such a high occupation fraction suggests the BH seeds formed in the early epoch of cosmic assembly likely resulted in abundant seeds, favoring a low-mass seed mechanism of the remnants, most likely from the first generation of massive stars. We find dynamical masses of the NSCs ranging from $2-80\times10^6~M_{\odot}$ and compare these masses to scaling relations for NSCs based primarily on photometric mass estimates. Color gradients suggest younger stellar populations lie at the centers of the NSCs in three of the four galaxies (NGC205, NGC5102, and NGC5206), while the morphology of two are complex and are best-fit with multiple morphological components (NGC5102 and NGC5206). The NSC kinematics show they are rotating, especially in M32 and NGC5102 ($V/σ_{\star}\sim0.7$). △ Less

Submitted 12 April, 2018; v1 submitted 12 November, 2017; originally announced November 2017.

Comments: 30 pages, 13 figures, 12 tables. Accepted to ApJ

A Multi-wavelength Study of the Turbulent Central Engine of the Low-mass AGN hosted by NGC404

Authors: Kristina Nyland, Timothy A. Davis, Dieu D. Nguyen, Anil Seth, Joan M. Wrobel, Atish Kamble, Mark Lacy, Katherine Alatalo, Margarita Karovska, W. Peter Maksym, Dipanjan Mukherjee, Lisa M. Young

Abstract: The nearby dwarf galaxy NGC404 harbors a low-luminosity active galactic nucleus (AGN) powered by the lowest-mass (< 150,000 solar-masses) central massive black hole (MBH) with a dynamical mass constraint currently known, thus providing a rare low-redshift analog to the MBH "seeds" that formed in the early Universe. Here, we present new imaging of the nucleus of NGC404 at 12-18 GHz with the Karl G.… ▽ More The nearby dwarf galaxy NGC404 harbors a low-luminosity active galactic nucleus (AGN) powered by the lowest-mass (< 150,000 solar-masses) central massive black hole (MBH) with a dynamical mass constraint currently known, thus providing a rare low-redshift analog to the MBH "seeds" that formed in the early Universe. Here, we present new imaging of the nucleus of NGC404 at 12-18 GHz with the Karl G. Jansky Very Large Array (VLA) and observations of the CO(2-1) line with the Atacama Large Millimeter/Submillimeter Array (ALMA). For the first time, we have successfully resolved the nuclear radio emission, revealing a centrally peaked, extended source spanning 17 pc. Combined with previous VLA observations, our new data place a tight constraint on the radio spectral index and indicate an optically-thin synchrotron origin for the emission. The peak of the resolved radio source coincides with the dynamical center of NGC404, the center of a rotating disk of molecular gas, and the position of a compact, hard X-ray source. We also present evidence for shocks in the NGC404 nucleus from archival narrowband HST imaging, Chandra X-ray data, and Spitzer mid-infrared spectroscopy, and discuss possible origins for the shock excitation. Given the morphology, location, and steep spectral index of the resolved radio source, as well as constraints on nuclear star formation from the ALMA CO(2-1) data, we find the most likely scenario for the origin of the radio source in the center of NGC404 to be a radio outflow associated with a confined jet driven by the active nucleus. △ Less

Submitted 7 July, 2017; originally announced July 2017.

Comments: Accepted to ApJ - 19 pages, 9 figures

Improved Dynamical Constraints on the Mass of the central Black Hole in NGC 404

Authors: Dieu D. Nguyen, Anil C. Seth, Mark den Brok, Nadine Neumayer, Michele Cappellari, Aaron J. Barth, Nelson Caldwell, Benjamin F. Williams, Breanna Binder

Abstract: We explore the nucleus of the nearby 10$^9M_{\odot}$~early-type galaxy (ETGs), NGC~404, using Hubble Space Telescope (HST)/STIS spectroscopy and WFC3 imaging. We first present evidence for nuclear variability in UV, optical, and infrared filters over a time period of 15~years. This variability adds to the already substantial evidence for an accreting black hole at the center of NGC~404. We then re… ▽ More We explore the nucleus of the nearby 10$^9M_{\odot}$~early-type galaxy (ETGs), NGC~404, using Hubble Space Telescope (HST)/STIS spectroscopy and WFC3 imaging. We first present evidence for nuclear variability in UV, optical, and infrared filters over a time period of 15~years. This variability adds to the already substantial evidence for an accreting black hole at the center of NGC~404. We then redetermine the dynamical black hole mass in NGC~404 including modeling of the nuclear stellar populations. We combine HST/STIS spectroscopy with WFC3 images to create a local color-$M/L$~relation derived from stellar population modeling of the STIS data. We then use this to create a mass model for the nuclear region. We use Jeans modeling to fit this mass model to adaptive optics (AO) stellar kinematic observations from Gemini/NIFS. From our stellar dynamical modeling, we find a 3$σ$ upper limit on the black hole mass of $1.5\times10^5M_{\odot}$. Given the accretion evidence for a black hole, this upper limit makes NGC~404 the lowest mass central black hole with dynamical mass constraints. We find that the kinematics of H$_2$ emission line gas show evidence for non-gravitational motions preventing the use of gas dynamical modeling to constrain the black hole mass. Our stellar population modeling also reveals that the central, counter-rotating region of the nuclear cluster is dominated by $\sim$1~Gyr old populations. △ Less

Submitted 3 November, 2016; v1 submitted 28 October, 2016; originally announced October 2016.

Comments: 23 papes, 16 figures, 6 tables, Resubmitted after accounting for referee comments

Extended Structure and Fate of the Nucleus in Henize 2-10

Authors: Dieu D. Nguyen, Anil C. Seth, Amy E. Reines, Mark den Brok, David Sand, Brian McLeod

Abstract: We investigate the structure and nuclear region of the black hole (BH) hosting galaxy Henize2-10. Surface brightness (SB) profiles are analyzed using Magellan/Megacam $g$- and $r$-band images. Excluding the central starburst, we find a best-fit two component Sérsic profile with $n_{\rm in} \sim 0.6$, $r_{\text{eff,in}} \sim$ 260~pc, and $n_{\rm out}\sim 1.8$, $r_{\text{eff, out}}\sim$ 1 kpc. Integ… ▽ More We investigate the structure and nuclear region of the black hole (BH) hosting galaxy Henize2-10. Surface brightness (SB) profiles are analyzed using Magellan/Megacam $g$- and $r$-band images. Excluding the central starburst, we find a best-fit two component Sérsic profile with $n_{\rm in} \sim 0.6$, $r_{\text{eff,in}} \sim$ 260~pc, and $n_{\rm out}\sim 1.8$, $r_{\text{eff, out}}\sim$ 1 kpc. Integrating out to our outer most data point (100 arc sec $\sim$ 4.3 kpc), we calculate $M_g=-19.2$ and $M_r=-19.8$. The corresponding enclosed stellar mass is $M_{\star}\sim(10\pm3)\times10^9$ M$_\odot$, $\sim3\times$ larger than previous estimates. Apart from the central $\lesssim$500 pc, with blue colors and an irregular morphology, the galaxy appears to be an early-type system. The outer color is quite red, $(g-r)_0=0.75$, suggesting a dominant old population. We study the nuclear region of the galaxy using archival Gemini/NIFS $K$-band adaptive optics spectroscopy and {\it Hubble Space Telescope} imaging. We place an upper limit on the BH mass of $\sim10^7M_{\odot}$ from the NIFS data, consistent with that from the $M_{\rm BH}$-radio-X-ray fundamental plane. No coronal lines are seen, but a Br$γ$ source is located at the position of the BH with a luminosity consistent with the X-ray emission. The starburst at the center of Henize 2-10 has led to the formation of several super star clusters, which are within $\sim$100 pc of the BH. We examine the fate of the nucleus by estimating the dynamical masses and dynamical friction timescales of the clusters. The most massive clusters ($\sim 10^6 M_\odot$) have $τ_{\rm dyn} \lesssim 200$ Myr, and thus Henize 2-10 may represent a rare snapshot of nuclear star cluster formation around a pre-existing massive BH. △ Less

Submitted 21 August, 2014; v1 submitted 19 August, 2014; originally announced August 2014.

Comments: 15 pages, 9 figures, ApJ accepted