Researcher profile

Santiago García-Burillo

Santiago García-Burillo contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
0followers
1topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2026arXiv

Evidence of Feedback Effects in Low-luminosity Active Galactic Nuclei Revealed by JWST Spectroscopy

This letter presents an analysis of the infrared ($\sim 3-28\,μm$) spectra extracted from the nuclear ($r < 150$ pc) regions of four low-luminosity active galactic nuclei (AGN), observed by JWST NIRSpec/IFU and MIRI/MRS as an extension of the Galaxy Activity, Torus, and Outflow Survey (GATOS). We find that, compared to higher-luminosity AGN, these low-luminosity AGN exhibit distinct properties in their emission of ionized gas, polycyclic aromatic hydrocarbons (PAHs), and molecular hydrogen (H$_2$). Specifically, the low-luminosity AGN exhibit relatively weak high-ionization potential lines (e.g., [Ne V] and [O IV]), and the line ratios suggest that fast radiative shocks (with $v_{\rm s}$ of $\sim \rm 100s\,km\,s^{-1}$) are the primary excitation source of ionized gas therein. Under the low-excitation conditions of their nuclear regions, these low-luminosity AGN generally exhibit a higher fraction of PAHs with large size ($N_{\rm C} \gtrsim 200$), reflecting the preferential destruction of smaller PAH molecules by AGN feedback. Furthermore, the H$_2$ transitions in these low-luminosity AGN are not fully thermalized, with slow, plausibly jet-driven molecular shocks (with $v_{\rm s} \leq \rm 10\,km\,s^{-1}$) likely being the extra excitation source. Taken together with results from the literature, these findings indicate that feedback operates in both low- and high-luminosity AGN, albeit its impact varies with AGN luminosity. In particular, systematic variations in PAH band ratios are found across AGN, demonstrating the differing influence of feedback in AGN of varying luminosities and highlighting the potential of PAH band ratios as diagnostics for distinguishing kinetic- and radiative-mode AGN feedback.

preprint2025arXiv

Galaxy Activity, Torus and Outflow Survey (GATOS) X: Molecular gas clumpiness under the influence of AGN

The distribution of molecular gas on small scales regulates star formation and the growth of supermassive black holes in galaxy centers, yet the role of active galactic nuclei (AGN) feedback in shaping this distribution remains poorly constrained. We investigate how AGN influence the small-scale structure of molecular gas in galaxy centers, by measuring the clumpiness of CO(3 - 2) emission observed with the Atacama Large Millimeter/submillimeter Array (ALMA) in the nuclear regions (50 - 200 pc from the AGN) of 16 nearby Seyfert galaxies from the Galaxy Activity, Torus, and Outflow Survey (GATOS). To quantify clumpiness, we apply three different methods: (1) the median of the pixel-by-pixel contrast between the original and smoothed maps; (2) the ratio of the total excess flux to the total flux, after substracting the background smoothed emission; and (3) the fraction of total flux coming from clumpy regions, interpreted as the mass fraction in clumps. We find a negative correlation between molecular gas clumpiness and AGN X-ray luminosity (L_X), suggesting that higher AGN activity is associated with smoother gas distributions. All methods reveal a turnover in this relation around L_X = 10^{42} erg/s, possibly indicating a threshold above which AGN feedback becomes efficient at dispersing dense molecular structures and suppressing future star formation. Our findings provide new observational evidence that AGN feedback can smooth out dense gas structures in galaxy centers.

preprint2025arXiv

Impact of AGN and nuclear star formation on the ISM turbulence of galaxies: Insights from JWST/MIRI spectroscopy

Active galactic nuclei (AGN), star formation (SF), and galaxy interactions can drive turbulence in the gas of the ISM, which in turn plays a role in the SF within galaxies. The impact on molecular gas is of particular importance, as it serves as the primary fuel for SF. Our goal is to investigate the origin of turbulence and the emission of molecular gas, as well as low- and intermediate-ionization gas, in the inner few kpc of both AGN hosts and SF galaxies. We use JWST MIRI/MRS observations of a sample consisting of 54 galaxies at z<0.1. We present fluxes of the H2 S(5)6.9091, [Ar II]6.9853, [FeII]5.3403, and [Ar III]8.9914 lines, along with velocity dispersion from W80. For galaxies with coronal emission, [Mg V]5.6098 is also included. Line ratios are compared to photoionization and shock models to explore the origin of the gas emission. AGNs exhibit broader emission lines than SFGs, with the largest velocity dispersions observed in radio-strong (RS) AGNs. H2 gas is less turbulent compared to ionized gas, while coronal gas presents higher velocity dispersions. The W80 values for the ionized gas exhibits a decrease from the nucleus out to radii of approximately 0.5--1 kpc, followed by an outward increase up to 2-3 kpc. In contrast, the H2 line widths generally display increasing profiles with distance from the center. Correlations W80 and line ratios such as H2 S(5)/[ArII] and [FeII]/[ArII] indicate that the most turbulent gas is associated with shocks, enhancing H2 and [FeII] emissions. We speculate that these shocked gas regions are produced by AGN outflows and jet-cloud interactions in AGN-dominated sources, while in SFGs, they may be created by stellar winds and mergers. This shock-induced gas heating may be an important mechanism of AGN (or stellar) feedback, preventing the gas from cooling and forming new stars.

preprint2025arXiv

NGC 6860, Mrk 915, and MCG -01-24-012. II. Inflowing and outflowing cold molecular gas and the connection with ionized gas in Seyfert galaxies

We present a study of the cold molecular gas kinematics in the inner ~ 4-7 kpc (projected sizes) of three nearby Seyfert galaxies, with AGN luminosities of ~ 10$^{44}$ erg/s, using observations of the CO(2-1) emission line, obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) at ~ 0.5-0.8$^{\prime\prime}$ (~ 150-400 pc) spatial resolutions. After modeling the CO profiles with multiple Gaussian components, we detected regions with double-peak profiles that exhibit kinematics distinct from the dominant rotational motion. In NGC 6860, a molecular outflow surrounding the bipolar emission of the [O III] ionized gas is observed extending up to $R_{out}$ ~ 560 pc from the nucleus. There is evidence of molecular inflows along the stellar bar, although an alternative scenario, involving a decoupled rotation in a circumnuclear disk (CND) can also explain the observed kinematics. Mrk 915 shows double-peak CO profiles along one of its spiral arms. Due to its ambiguous disk orientation, part of the CO emission can be interpreted as a molecular gas inflow or an outflow reaching $R_{out}$ ~ 2.8 kpc. MCG -01-24-012 has double-peak profiles associated with a CND, perpendicular to the [O III] bipolar emission. The CO in the CND is rotating while outflowing within $R_{out}$ ~ 3 kpc, with the disturbances possibly being caused by the passage of the ionized gas outflow. Overall, the mass inflow rates are larger than the accretion rate needed to produce the observed luminosities, suggesting that only a fraction of the inflowing gas ends up feeding the central black holes. Although we found signatures of AGN feedback on the cold molecular phase, the mass outflow rates of ~ 0.09-3 M$_\odot$/yr indicate an overall weak impact at these AGN luminosities. Nonetheless, we may be witnessing the start of the depletion and ejection of the molecular gas reservoir that has accumulated over time.

preprint2025arXiv

NOEMA$^{\rm 3D}$: A first kpc resolution study of a $z\sim1.5$ main sequence barred galaxy channeling gas into a growing bulge

We present a very deep CO(3-2) observation of a massive, gas-rich, main sequence, barred spiral galaxy at $z\approx1.52$. Our data were taken with the IRAM-NOEMA interferometer for a 12-antenna equivalent on-source integration time of $\sim$ 50 hours. We fit the major axis kinematics using forward modelling of a rotating disk, and then subtract the two-dimensional beam convolved best-fit model revealing signatures of planar non-circular motions in the residuals. The inferred in-plane radial velocities are remarkably large, of the order of $\approx60$ km/s. Direct comparisons with a high-resolution, simulated, gas-rich, barred galaxy, obtained with the moving mesh code AREPO and the TNG sub-grid model, show that the observed non-circular gas flows can be explained as radial flows driven by the central bar, with an inferred net inflow rate of the order of the SFR. Given the recent evidence for a higher-than-expected fraction of barred disk galaxies at cosmic noon, our results suggest that rapid gas inflows due to bars could be important evolutionary drivers for the dominant population of star-forming galaxies at the peak epoch of star and galaxy formation.

preprint2022arXiv

ALMA imaging of the cold molecular and dusty disk in the type 2 active nucleus of the Circinus galaxy

We aim to shed light on the physical properties and kinematics of the molecular material in the nucleus of one of the closest type 2 active galaxies. To this end, we obtained high angular resolution ALMA observations of the nucleus of the Circinus galaxy. The observations map the emission at 350GHz and 690GHz with spatial resolutions of ~3.8pc and ~2.2pc, respectively. The continuum emission traces cold ($T\lesssim100$K) dust in a circumnuclear disk with spiral arms on scales of 25pc, plus a marginally resolved nuclear emission peak. The latter is not extended in polar direction as claimed based on earlier ALMA observations. A significant amount (of the order of 40%) of the 350GHz emission is not related to dust, but most likely free-free emission instead. We detect CO(3-2) and CO(6-5) as well as HCO$^+$(4-3), HCN(4-3), and CS(4-3). The CO emission is extended, showing a spiral pattern, similar to the extended dust emission. Towards the nucleus, CO is excited to higher transitions and its emission is self-absorbed, leading to an apparent hole in the CO(3-2) but not the CO(6-5) emission. On the other hand, the high gas density tracers HCO$^+$, HCN, and CS show a strong, yet unresolved (($\lesssim4$pc) concentration of the emission at the nucleus, pointing at a very small &#39;torus&#39;. The kinematics are dominated by rotation and point at a geometrically thin disk down to the resolution limit of our observations. In contrast to several other AGNs, no HCN enhancement is found towards the nucleus. The Circinus nucleus is therefore composed of at least two distinct components: (1) an optically thin, warm outflow of ionised gas containing clouds of dust; and (2) a cold molecular and dusty disk. These findings support the most recent radiative transfer calculations of the obscuring structures in AGNs, which find a similar two-component structure. (Abridged)

preprint2022arXiv

The complex infrared dust continuum emission of NGC1068: ground-based N- and Q-band spectroscopy and new radiative transfer models

Thanks to ground-based infrared and sub-mm observations the study of the dusty torus of nearby AGN has greatly advanced in the last years. With the aim of further investigating the nuclear mid-infrared emission of the archetypal Seyfert 2 galaxy NGC1068, here we present a fitting to the N- and Q-band Michelle/Gemini spectra. We initially test several available SED libraries, including a smooth, clumpy and two phase torus models, and a clumpy disk plus wind model. Although, the smooth torus model describe the spectra of NGC1068 if we allow to vary some model parameters among the two spectral bands. Motivated by this result, we produced new SEDs using the radiative transfer code SKIRT. We use two concentric tori that allow us to test a more complex geometry. We test different values for the inner and outer radii, half opening angle, radial and polar exponent of the power-law density profile, opacity, and viewing angle. Furthermore, we also test the dust grains size and different optical and calorimetric properties of silicate grains. The best fitting model consists of two concentric components with outer radii of 1.8 and 28 pc, respectively. We find that the size and the optical and calorimetric properties of graphite and silicate grains in the dust structure are key to reproduce the spectra of NGC1068. We conclude that the dust in NGC1068 reaches different scales, where the highest contribution to the mid-infrared is given by a central and compact component. A less dense and extended component is present, which can be either part of the same torus (conforming a flared disk) or can represent the emission of a polar dust component, as already suggested from interferometric observations.

preprint2021arXiv

The Diverse Molecular Gas Content of Massive Galaxies Undergoing Quenching at z~1

We present a detailed study of the molecular gas content and stellar population properties of three massive galaxies at 1 < z < 1.3 that are in different stages of quenching. The galaxies were selected to have a quiescent optical/near-infrared spectral energy distribution and a relatively bright emission at 24 micron, and show remarkably diverse properties. CO emission from each of the three galaxies is detected in deep NOEMA observations, allowing us to derive molecular gas fractions Mgas/Mstar of 13-23%. We also reconstruct the star formation histories by fitting models to the observed photometry and optical spectroscopy, finding evidence for recent rejuvenation in one object, slow quenching in another, and rapid quenching in the third system. To better constrain the quenching mechanism we explore the depletion times for our sample and other similar samples at z~0.7 from the literature. We find that the depletion times are highly dependent on the method adopted to measure the star formation rate: using the UV+IR luminosity we obtain depletion times about 6 times shorter than those derived using dust-corrected [OII] emission. When adopting the star formation rates from spectral fitting, which are arguably more robust, we find that recently quenched galaxies and star-forming galaxies have similar depletion times, while older quiescent systems have longer depletion times. These results offer new, important constraints for physical models of galaxy quenching.

preprint2021arXiv

Vibrationally excited HC3N emission in NGC1068: Tracing the recent star formation in the starburst ring

Using ALMA data, we have studied the HC$_3$N and continuum emission in the starburst ring (SB ring) and the circumnuclear disc (CND) of the SB/AGN composite galaxy NGC 1068. We have detected emission from vibrationally excited HC$_3$N (HC$_3$N*) only towards one star-forming region of the SB ring. Remarkably, HC$_3$N* was not detected towards the CND despite its large HC$_3$N $v=0$ column density. From LTE and non-LTE modelling of HC$_3$N*, we obtained a dust temperature of $ T_\text{dust} \sim 250$ K and a density of $n_{\text{H}_2}=6\times10^5$ cm$^{-3}$ for this star-forming region. The estimated IR luminosity of $5.8\times10^8$ L$_\odot$ is typical of proto-Super Star Clusters (proto-SSC) observed in the SB galaxy NGC 253. We use the continuum emissions at $147$ GHz and $350$ GHz, along with CO and Pa $α$, to estimate the ages of other $14$ SSCs in the SB ring. We find the youngest SSCs to be associated with the region connecting the nuclear bar with the SB ring, supporting the scenario of sequential star formation. For the CND, our analysis yields $T_\text{dust} \leqslant 100$ K and $n_{\text{H}_2}\sim(3-6)\times10^5$cm$^{-3}$. The very different dust temperatures found for the CND and the proto-SSC indicates that, while the dust in the proto-SSC is being efficiently heated from the inside by the radiation from massive proto-stars, the CND is being heated externally by the AGN, which in the IR optically thin case can only heat the dust to $56$ K. We discuss the implications of the non-detection of HC$_3$N* near the luminous AGN in NGC 1068 on the interpretation of the HC$_3$N* emission observed in the SB/AGN composite galaxies NGC 4418 and Arp 220.

preprint2020arXiv

A proto-pseudobulge in ESO 320-G030 fed by a massive molecular inflow driven by a nuclear bar

Galaxies with nuclear bars are believed to efficiently drive gas inward, generating a nuclear starburst and possibly an active galactic nucleus (AGN). We confirm this scenario for the isolated, double-barred, luminous infrared galaxy ESO 320-G030 based on an analysis of Herschel and ALMA spectroscopic observations. Herschel/PACS and SPIRE observations of ESO 320-G030 show absorption/emission in 18 lines of H2O, which we combine with the ALMA H2O 423-330 448 GHz line (Eupper~400 K) and continuum images to study the nuclear region. Radiative transfer models indicate that 3 nuclear components are required to account for the H2O and continuum data. An envelope, with R~130-150 pc, T_dust~50 K, and N_H2~2x10^{23} cm^{-2}, surrounds a nuclear disk with R~40 pc and tau_100um~1.5-3 (N_H2~2x10^{24} cm^{-2}) and an extremely compact (R~12 pc), warm (~100 K), and buried (tau_100um>5, N_H2>~5x10^{24} cm^{-2}) core component. The three nuclear components account for 70% of the galaxy L_IR (SFR~16-18 Msun yr^{-1}). The nucleus is fed by a molecular inflow observed in CO 2-1 with ALMA, which is associated with the nuclear bar. With decreasing radius (r=450-225 pc), the mass inflow rate increases up to ~20 Msun yr^{-1}, which is similar to the nuclear SFR, indicating that the starburst is sustained by the inflow. At lower r, the inflow is best probed by the far-infrared OH ground-state doublets, with an estimated inflow rate of ~30 Msun yr^{-1}. The short timescale of ~20 Myr for nuclear gas replenishment indicates quick secular evolution, and indicates that we are witnessing an intermediate stage (<100 Myr) proto-pseudobulge fed by a massive inflow that is driven by a strong nuclear bar. We also apply the H2O model to the Herschel far-infrared spectroscopic observations of H2^{18}O, OH, $^{18}OH, OH+, H2O^+, H3O^+, NH, NH2, NH3, CH, CH^+, ^{13}CH^+, HF, SH, and C3, and estimate their abundances.

preprint2020arXiv

ALMA polarimetry measures magnetically aligned dust grains in the torus of NGC 1068

The obscuring structure surrounding active galactic nuclei (AGN) can be explained as a dust and gas flow cycle that fundamentally connects the AGN with their host galaxies. This structure is believed to be associated with dusty winds driven by radiation pressure. However, the role of magnetic fields, which are invoked in almost all models for accretion onto a supermassive black hole and outflows, is not thoroughly studied. Here we report the first detection of polarized thermal emission by means of magnetically aligned dust grains in the dusty torus of NGC 1068 using ALMA Cycle 4 polarimetric dust continuum observations ($0.07&#34;$, $4.2$ pc; 348.5 GHz, $860$ $μ$m). The polarized torus has an asymmetric variation across the equatorial axis with a peak polarization of $3.7\pm0.5$\% and position angle of $109\pm2^{\circ}$ (B-vector) at $\sim8$ pc east from the core. We compute synthetic polarimetric observations of magnetically aligned dust grains assuming a toroidal magnetic field and homogeneous grain alignment. We conclude that the measured 860 $μ$m continuum polarization arises from magnetically aligned dust grains in an optically thin region of the torus. The asymmetric polarization across the equatorial axis of the torus arises from 1) an inhomogeneous optical depth, and 2) a variation of the velocity dispersion, i.e. variation of the magnetic field turbulence at sub-pc scales, from the eastern to the western region of the torus. These observations and modeling constrain the torus properties beyond spectral energy distribution results. This study strongly supports that magnetic fields up to a few pc contribute to the accretion flow onto the active nuclei.

preprint2020arXiv

Plateau de Bure High-z Blue-Sequence Survey 2 (PHIBSS2): Search for Secondary Sources, CO Luminosity Functions in the Field, and the Evolution of Molecular Gas Density through Cosmic Time

We report on the results of a search for serendipitous sources in CO emission in 110 cubes targeting CO(2-1), CO(3-2), and CO(6-5) at z ~ 1-2 from the second Plateau de Bure High-z Blue-Sequence Survey (PHIBSS2). The PHIBSS2 observations were part of a 4-year legacy program at the IRAM Plateau de Bure Interferometer aimed at studying early galaxy evolution from the perspective of molecular gas reservoirs. We present a catalog of 67 candidate secondary sources from this search, with 45 out of the 110 data cubes showing sources in addition to the primary target that appear to be field detections, unrelated to the central sources. This catalog includes the redshifts, line widths, fluxes, as well as an estimation of their reliability based on their false positive probability. We perform a search in the 3D-HST/CANDELS catalogs for the secondary CO detections and tentatively find that ~64% of these have optical counterparts, which we use to constrain their redshifts. Finally, we use our catalog of candidate CO detections to derive the CO(2-1), CO(3-2), CO(4-3), CO(5-4), and CO(6-5) luminosity functions over a range of redshifts, as well as the molecular gas mass density evolution. Despite the different methodology, these results are in very good agreement with previous observational constraints derived from blind searches in deep fields. They provide an example of the type of &#34;deep field&#34; science that can be carried out with targeted observations.