# Publications

## Eccentric Black Hole Mergers in Active Galactic Nuclei

ASTROPHYSICAL JOURNAL LETTERS 907 (2021) ARTN L20

H Tagawa, B Kocsis, Z Haiman, I Bartos, K Omukai, J Samsing

## The Evolution of Gas-Phase Metallicity and Resolved Abundances in Star-forming Galaxies at z ≈ 0.6 – 1.8

Monthly Notices of the Royal Astronomical Society Oxford University Press (OUP) (2020)

S Gillman, A Tiley, A Swinbank, U Dudzevičiūtė, R Sharples, I Smail, C Harrison, AJ Bunker, M Bureau, M Cirasuolo, GE Magdis, T Mendel, JP Stott

&lt;jats:title&gt;Abstract&lt;/jats:title&gt; &lt;jats:p&gt;We present an analysis of the chemical abundance properties of &#x2248;650 star-forming galaxies at z&#xA0;&#x2248;&#xA0;0.6&#xA0;&#x2013;&#xA0;1.8. Using integral-field observations from the K&#xA0;-&#xA0;band Multi-Object Spectrograph (KMOS), we quantify the [N&#x2009;ii]/H&#x3B1; emission-line ratio, a proxy for the gas-phase Oxygen abundance within the interstellar medium. We define the stellar mass&#xA0;&#x2013;&#xA0;metallicity relation at z&#xA0;&#x2248;&#xA0;0.6&#xA0;&#x2013;&#xA0;1.0 and z&#xA0;&#x2248;&#xA0;1.2&#xA0;&#x2013;&#xA0;1.8 and analyse the correlation between the scatter in the relation and fundamental galaxy properties (e.g. H&#x3B1; star-formation rate, H&#x3B1; specific star-formation rate, rotation dominance, stellar continuum half-light radius and Hubble-type morphology). We find that for a given stellar mass, more highly star-forming, larger and irregular galaxies have lower gas-phase metallicities, which may be attributable to their lower surface mass densities and the higher gas fractions of irregular systems. We measure the radial dependence of gas-phase metallicity in the galaxies, establishing a median, beam smearing-corrected, metallicity gradient of &#x394;Z/&#x394;R=&#xA0;0.002&#xA0;&#xB1;&#xA0;0.004 dex kpc&#x2212;1, indicating on average there is no significant dependence on radius. The metallicity gradient of a galaxy is independent of its rest-frame optical morphology, whilst correlating with its stellar mass and specific star-formation rate, in agreement with an inside-out model of galaxy evolution, as well as its rotation dominance. We quantify the evolution of metallicity gradients, comparing the distribution of &#x394;Z/&#x394;R in our sample with numerical simulations and observations at z&#xA0;&#x2248;&#xA0;0&#xA0;&#x2013;&#xA0;3. Galaxies in our sample exhibit flatter metallicity gradients than local star-forming galaxies, in agreement with numerical models in which stellar feedback plays a crucial role redistributing metals.&lt;/jats:p&gt;

## WISDOM project – VI. Exploring the relation between supermassive black hole mass and galaxy rotation with molecular gas

Monthly Notices of the Royal Astronomical Society Oxford University Press (OUP) 500 (2020) 1933-1952

MD Smith, M Bureau, TA Davis, M Cappellari, L Liu, K Onishi, S Iguchi, EV North, M Sarzi

<jats:title>ABSTRACT</jats:title> <jats:p>Empirical correlations between the masses of supermassive black holes (SMBHs) and properties of their host galaxies are well established. Among these is the correlation with the flat rotation velocity of each galaxy measured either at a large radius in its rotation curve or via a spatially integrated emission-line width. We propose here the use of the deprojected integrated CO emission-line width as an alternative tracer of this rotation velocity, which has already been shown useful for the Tully–Fisher (luminosity–rotation velocity) relation. We investigate the correlation between CO line widths and SMBH masses for two samples of galaxies with dynamical SMBH mass measurements, with spatially resolved and unresolved CO observations, respectively. The tightest correlation is found using the resolved sample of 25 galaxies as $\log (M_\mathrm{BH}/\mathrm{M_\odot })=(7.5\pm 0.1)+(8.5\pm 0.9)[\log (W_\mathrm{50}/\sin i \, \mathrm{km\, s}^{-1})-2.7]$, where MBH is the central SMBH mass, W50 is the full width at half-maximum of a double-horned emission-line profile, and i is the inclination of the CO disc. This relation has a total scatter of $0.6\,$ dex, comparable to those of other SMBH mass correlations, and dominated by the intrinsic scatter of $0.5\,$ dex. A tight correlation is also found between the deprojected CO line widths and the stellar velocity dispersions averaged within one effective radius. We apply our correlation to the COLD GASS sample to estimate the local SMBH mass function.</jats:p>

## Jeans modelling of the Milky Way’s nuclear stellar disc

Monthly Notices of the Royal Astronomical Society Oxford University Press (OUP) (2020)

MC Sormani, J Magorrian, F Nogueras-Lara, N Neumayer, R Schönrich, RS Klessen, A Mastrobuono-Battisti

&lt;jats:title&gt;Abstract&lt;/jats:title&gt; &lt;jats:p&gt;The nuclear stellar disc (NSD) is a flattened stellar structure that dominates the gravitational potential of the Milky Way at Galactocentric radii 30 &#x2272; R &#x2272; 300&#x2009;pc. In this paper, we construct axisymmetric Jeans dynamical models of the NSD based on previous photometric studies and we fit them to line-of-sight kinematic data of APOGEE and SiO maser stars. We find that (i) the NSD mass is lower but consistent with the mass independently determined from photometry by Launhardt et&#xA0;al. (2002). Our fiducial model has a mass contained within spherical radius r = 100&#x2009;pc of $M(r&amp;amp;lt;100\, {\rm pc}) = 3.9 \pm 1 \times 10^8 \, \rm M_\odot$ and a total mass of $M_{\rm NSD} = 6.9 \pm 2 \times 10^8 \, \rm M_\odot$. (ii) The NSD might be the first example of a vertically biased disc, i.e. with ratio between the vertical and radial velocity dispersion &#x3C3;z/&#x3C3;R &amp;amp;gt; 1. Observations and theoretical models of the star-forming molecular gas in the central molecular zone suggest that large vertical oscillations may be already imprinted at stellar birth. However, the finding &#x3C3;z/&#x3C3;R &amp;amp;gt; 1 depends on a drop in the velocity dispersion in the innermost few tens of parsecs, on our assumption that the NSD is axisymmetric, and that the available (extinction corrected) stellar samples broadly trace the underlying light and mass distributions, all of which need to be established by future observations and/or modelling. (iii) We provide the most accurate rotation curve to date for the innermost 500&#x2009;pc of our Galaxy.&lt;/jats:p&gt;

## SDSS-IV MaNGA: the indispensable role of bars in enhancing the central star formation of low-z galaxies

Monthly Notices of the Royal Astronomical Society Oxford University Press (OUP) (2020)

L Lin, C Li, C Du, E Wang, T Xiao, M Bureau, A Fraser-McKelvie, K Masters, L Lin, D Wake, L Hao

&lt;jats:title&gt;Abstract&lt;/jats:title&gt; &lt;jats:p&gt;We analyse two-dimensional maps and radial profiles of EW(H&#x3B1;), EW(H&#x3B4;A), and Dn(4000) of low-redshift galaxies using integral field spectroscopy from the MaNGA survey. Out of &#x2248;1400 nearly face-on late-type galaxies with a redshift z &amp;amp;lt; 0.05, we identify 121 &#x201C;turnover&#x201D; galaxies that each have a central upturn in EW(H&#x3B1;), EW(H&#x3B4;A) and/or a central drop in Dn(4000), indicative of ongoing/recent star formation. The turnover features are found mostly in galaxies with a stellar mass above &#x223C;1010 M&#x2299; and NUV-r colour less than &#x2248;5. The majority of the turnover galaxies are barred, with a bar fraction of 89&#xB1;3%. Furthermore, for barred galaxies the radius of the central turnover region is found to tightly correlate with one third of the bar length. Comparing the observed and the inward extrapolated star formation rate surface density, we estimate that the central SFR have been enhanced by an order of magnitude. Conversely, only half of the barred galaxies in our sample have a central turnover feature, implying that the presence of a bar is not sufficient to lead to a central SF enhancement. We further examined the SF enhancement in paired galaxies, as well as the local environment, finding no relation. This implies that environment is not a driving factor for central SF enhancement in our sample. Our results reinforce both previous findings and theoretical expectation that galactic bars play a crucial role in the secular evolution of galaxies by driving gas inflow and enhancing the star formation and bulge growth in the center.&lt;/jats:p&gt;

## Detecting Kozai–Lidov Imprints on the Gravitational Waves of Intermediate-mass Black Holes in Galactic Nuclei

The Astrophysical Journal American Astronomical Society 901 (2020) 125-125

B Deme, B-M Hoang, S Naoz, B Kocsis

## Binary Intermediate-mass Black Hole Mergers in Globular Clusters

The Astrophysical Journal American Astronomical Society 899 (2020) 149-149

A Rasskazov, G Fragione, B Kocsis

## Revealing the Intermediate Mass Black Hole at the Heart of Dwarf Galaxy NGC404 with Sub-parsec Resolution ALMA Observations

Monthly Notices of the Royal Astronomical Society Oxford University Press 496 (2020) 4061-4078

M Bureau, M Cappellari, L Liu, M Smith

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 ≈ 5 × 105 M 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×105 M (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.

## K-CLASH: Strangulation and ram pressure stripping in galaxy cluster members at 0.3 &lt; z &lt; 0.6

Monthly Notices of the Royal Astronomical Society Oxford University Press (OUP) 496 (2020) 3841-3861

SP Vaughan, AL Tiley, RL Davies, LJ Prichard, SM Croom, M Bureau, JP Stott, A Bunker, M Cappellari, B Ansarinejad, MJ Jarvis

&lt;jats:title&gt;ABSTRACT&lt;/jats:title&gt; &lt;jats:p&gt;Galaxy clusters have long been theorized to quench the star formation of their members. This study uses integral-field unit observations from the K-band MultiObject Spectrograph (KMOS) &#x2013; Cluster Lensing And Supernova survey with Hubble (CLASH) survey (K-CLASH) to search for evidence of quenching in massive galaxy clusters at redshifts 0.3 &amp;amp;lt; z &amp;amp;lt; 0.6. We first construct mass-matched samples of exclusively star-forming cluster and field galaxies, then investigate the spatial extent of their H&#x2009;&#x3B1; emission and study their interstellar medium conditions using emission line ratios. The average ratio of H&#x2009;&#x3B1; half-light radius to optical half-light radius ($r_{\mathrm{e}, {\rm {H}\,\alpha }}/r_{\mathrm{e}, R_{\mathrm{c} } }$) for all galaxies is 1.14&#xA0;&#xB1;&#xA0;0.06, showing that star formation is taking place throughout stellar discs at these redshifts. However, on average, cluster galaxies have a smaller $r_{\mathrm{e}, {\rm {H}\alpha }}/r_{\mathrm{e}, R_{\mathrm{c} } }$ ratio than field galaxies: &#x2329;$r_{\mathrm{e}, {\rm {H}\alpha }}/r_{\mathrm{e}, R_{\mathrm{c} } }$&#x232A;&#xA0;= 0.96&#xA0;&#xB1;&#xA0;0.09 compared to 1.22&#xA0;&#xB1;&#xA0;0.08 (smaller at a 98&#x2009;per&#x2009;cent credibility level). These values are uncorrected for the wavelength difference between H&#x2009;&#x3B1; emission and Rc-band stellar light but implementing such a correction only reinforces our results. We also show that whilst the cluster and field samples follow indistinguishable mass&#x2013;metallicity (MZ) relations, the residuals around the MZ relation of cluster members correlate with cluster-centric distance; galaxies residing closer to the cluster centre tend to have enhanced metallicities (significant at the 2.6&#x3C3; level). Finally, in contrast to previous studies, we find no significant differences in electron number density between the cluster and field galaxies. We use simple chemical evolution models to conclude that the effects of disc strangulation and ram-pressure stripping can quantitatively explain our observations.&lt;/jats:p&gt;

## Spin Evolution of Stellar-mass Black Hole Binaries in Active Galactic Nuclei

The Astrophysical Journal American Astronomical Society 899 (2020) 26-26

H Tagawa, Z Haiman, I Bartos, B Kocsis

## K-CLASH: spatially-resolving star-forming galaxies in field and cluster environments at z ≈ 0.2-0.6

Monthly Notices of the Royal Astronomical Society Oxford University Press (2020)

AL Tiley, JP Stott, R Davies, LJ Prichard, A Bunker, M Bureau, M Cappellari, M Jarvis, A Robotham, L Cortese, S Bellstedt, B Ansarinejad

We present the KMOS-CLASH (K-CLASH) survey, a K-band Multi-Object Spectrograph (KMOS) survey of the spatially-resolved gas properties and kinematics of 191 (predominantly blue) Hα-detected galaxies at 0.2 ≲ z ≲ 0.6 in field and cluster environments. K-CLASH targets galaxies in four Cluster Lensing And Supernova survey with Hubble (CLASH) fields in the KMOS IZ-band, over 7′ radius (≈2–3 Mpc) fields-of-view. K-CLASH aims to study the transition of star-forming galaxies from turbulent, highly star-forming disc-like and peculiar systems at z ≈ 1–3, to the comparatively quiescent, ordered late-type galaxies at z ≈ 0, and to examine the role of clusters in the build-up of the red sequence since z ≈ 1. In this paper, we describe the K-CLASH survey, present the sample, and provide an overview of the K-CLASH galaxy properties. We demonstrate that our sample comprises star-forming galaxies typical of their stellar masses and epochs, residing both in field and cluster environments. We conclude K-CLASH provides an ideal sample to bridge the gap between existing large integral-field spectroscopy surveys at higher and lower redshifts. We find that star-forming K-CLASH cluster galaxies at intermediate redshifts have systematically lower stellar masses than their star-forming counterparts in the field, hinting at possible “downsizing” scenarios of galaxy growth in clusters at these epochs. We measure no difference between the star-formation rates of Hα-detected, star-forming galaxies in either environment after accounting for stellar mass, suggesting that cluster quenching occurs very rapidly during the epochs probed by K-CLASH, or that star-forming K-CLASH galaxies in clusters have only recently arrived there, with insufficient time elapsed for quenching to have occured.

## Formation and Evolution of Compact-object Binaries in AGN Disks

ASTROPHYSICAL JOURNAL American Astronomical Society 898 (2020) ARTN 25

H Tagawa, Z Haiman, B Kocsis

## The 16th Data Release of the Sloan Digital Sky Surveys: First Release from the APOGEE-2 Southern Survey and Full Release of eBOSS Spectra

ASTROPHYSICAL JOURNAL SUPPLEMENT SERIES 249 (2020) ARTN 3

## The 16th data release of the Sloan Digital Sky Surveys: first release from the APOGEE-2 Southern Survey and full release of eBOSS spectra

Astrophysical Journal Supplement American Astronomical Society 249 (2020) 3

R Ahumada, C Allende Prieto, A Almeida, M Bureau, M Cappellari, R Davies, E-M Mueller, R Smethurst, SDSS-IVC SDSS-IV Collaboration

This paper documents the 16th data release (DR16) from the Sloan Digital Sky Surveys (SDSS), the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the Southern Hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey and new data from the SPectroscopic IDentification of ERosita Survey programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

## Cosmic Evolution of Stellar-mass Black Hole Merger Rate in Active Galactic Nuclei

ASTROPHYSICAL JOURNAL American Astronomical Society 896 (2020) ARTN 138

Y Yang, I Bartos, Z Haiman, B Kocsis, S Marka, H Tagawa

## Trapped orbits and solar-neighbourhood kinematics

MONTHLY NOTICES OF THE ROYAL ASTRONOMICAL SOCIETY 495 (2020) 895-904

J Binney

## Angle-action variables for orbits trapped at a Lindblad resonance

MONTHLY NOTICES OF THE ROYAL ASTRONOMICAL SOCIETY 495 (2020) 886-894

J Binney

## Electromagnetic transients and gravitational waves from white dwarf disruptions by stellar black holes in triple systems

Monthly Notices of the Royal Astronomical Society Oxford University Press (OUP) (2020)

G Fragione, BD Metzger, R Perna, NWC Leigh, B Kocsis

Mergers of binaries comprised of compact objects can give rise to explosive transient events, heralding the birth of exotic objects which cannot be formed through single star evolution. Using a large number of direct N-body simulations, we explore the possibility that a white dwarf (WD) is dynamically driven to tidal disruption by a stellar-mass black hole (BH) as a consequence of the joint effects of gravitational wave (GW) emission and Lidov-Kozai oscillations imposed by the tidal field of a outer tertiary companion orbiting the inner BH-WD binary. We explore the sensitivity of our results to the distributions of natal kick velocities imparted to the BH and WD upon formation, adiabatic mass loss, semi-major axes and eccentricities of the triples, and stellar mass ratios. We find rates of WD-TDEs in the range $1.2\times 10^{-3}-1.4$ Gpc$^{-3}$ yr$^{-1}$ for $z\leq 0.1$, rarer than stellar TDEs in triples by a factor of $\sim 3$--$30$. The uncertainty in the TDE rates may be greatly reduced in the future using gravitational wave (GW) observations of Galactic binaries and triples with LISA. WD-TDEs may give rise to high energy X-ray or gamma-ray transients of duration similar to long gamma-ray bursts but lacking the signatures of a core-collapse supernova, while being accompanied by a supernova-like optical transient which lasts for only days. WD--BH and WD--NS binaries will also emit GWs in the LISA band before the TDE. The discovery and identification of triple-induced WD-TDE events by future time domain surveys and/or GWs could enable the study of the demographics of BHs in nearby galaxies.

## A precise benchmark for cluster scaling relations: Fundamental Plane, Mass Plane, and IMF in the Coma cluster from dynamical models

Monthly Notices of the Royal Astronomical Society Oxford University Press 494 (2020) 5619-5635

S Shetty, M Cappellari, RM McDermid, D Krajnovic, PT de Zeeuw, RL Davies, C Kobayashi

We study a sample of 148 early-type galaxies in the Coma cluster using SDSS photometry and spectra, and calibrate our results using detailed dynamical models for a subset of these galaxies, to create a precise benchmark for dynamical scaling relations in high-density environments. For these galaxies, we successfully measured global galaxy properties, modelled stellar populations, and created dynamical models, and support the results using detailed dynamical models of 16 galaxies, including the two most massive cluster galaxies, using data taken with the SAURON IFU. By design, the study provides minimal scatter in derived scaling relations due to the small uncertainty in the relative distances of galaxies compared to the cluster distance. Our results demonstrate low (≤55 per cent for 90th percentile) dark matter fractions in the inner 1Re of galaxies. Owing to the study design, we produce the tightest, to our knowledge, IMF–σe relation of galaxies, with a slope consistent with that seen in local galaxies. Leveraging our dynamical models, we transform the classical Fundamental Plane of the galaxies to the Mass Plane. We find that the coefficients of the Mass Plane are close to predictions from the virial theorem, and have significantly lower scatter compared to the Fundamental Plane. We show that Coma galaxies occupy similar locations in the (M*–Re) and (M*−σe) relations as local field galaxies but are older. This, and the fact we find only three slow rotators in the cluster, is consistent with the scenario of hierarchical galaxy formation and expectations of the kinematic morphology–density relation.

## Intermediate-mass Black Holes' Effects on Compact Object Binaries

ASTROPHYSICAL JOURNAL American Astronomical Society 892 (2020) ARTN 130

B Deme, Y Meiron, B Kocsis

Although their existence is not yet confirmed observationally, intermediate mass black holes (IMBHs) may play a key role in the dynamics of galactic nuclei. In this paper, we neglect the effect of the nuclear star cluster itself and investigate only how a small reservoir of IMBHs influences the secular dynamics of stellar-mass black hole binaries, using N-body simulations. We show that our simplifications are valid and that the IMBHs significantly enhance binary evaporation by pushing the binaries into the Hill-unstable region of parameter space, where they are separated by the SMBH's tidal field. For binaries in the S-cluster region of the Milky Way, IMBHs drive the binaries to merge in up to 1-6% of cases, assuming five IMBHs within 5 pc of mass 10,000 solar masses each. Observations of binaries in the Galactic center may strongly constrain the population of IMBHs therein.