Search

We use the Hubble Space Telescope to analyse the extended [O III] 5007A emission in seven bright radio-quiet type 1 quasars (QSO1s), focusing on the morphology and physical conditions of their extended Narrow-Line Regions (NLRs). We find NLRs extending 3-9 kpc, with four quasars showing roughly symmetrical structures (b/a=1.2-1.5) and three displaying asymmetric NLRs (b/a=2.4-5.6). When included with type 1 and type 2 AGNs from previous studies, the sizes of the extended [O III] regions scale with luminosity as $R[O III] \sim L[O III]^{0.5}$, consistent with photoionisation. However, when analysed separately, type 1s exhibit a steeper slope ($\gamma=0.57\pm0.05$) compared to type 2 AGNs ($\gamma=0.48\pm0.02$). We use photoionisation modeling to estimate the maximum NLRs sizes, assuming a minimum ionisation parameter of $\log(U) = -3$, an ionising luminosity based on the $L[O III]$-derived bolometric luminosity, and a minimum gas number density $n_H \sim 100\,\text{cm}^{-3}$, assuming that molecular clouds provide a reservoir for the ionised gas. The derived sizes agree well with direct measurements for a sample of type 2 quasars, but are underestimated for the current sample of QSO1s. A better agreement is obtained for the QSO1s using bolometric luminosities derived from the 5100A continuum luminosity. Radial mass profiles for the QSO1s show significant extended mass in all cases, but with less [O III]-emitting gas near the central AGN compared to QSO2s. This may suggest that the QSO1s are in a later evolutionary stage than QSO2s, further past the blow-out stage., Comment: Accepted for publication on MNRAS

We investigate the impact of spatial survey non-uniformity on the galaxy redshift distributions for forthcoming data releases of the Rubin Observatory Legacy Survey of Space and Time (LSST). Specifically, we construct a mock photometry dataset degraded by the Rubin OpSim observing conditions, and estimate photometric redshifts of the sample using a template-fitting photo-$z$ estimator, BPZ, and a machine learning method, FlexZBoost. We select the Gold sample, defined as $i<25.3$ for 10 year LSST data, with an adjusted magnitude cut for each year and divide it into five tomographic redshift bins for the weak lensing lens and source samples. We quantify the change in the number of objects, mean redshift, and width of each tomographic bin as a function of the coadd $i$-band depth for 1-year (Y1), 3-year (Y3), and 5-year (Y5) data. In particular, Y3 and Y5 have large non-uniformity due to the rolling cadence of LSST, hence provide a worst-case scenario of the impact from non-uniformity. We find that these quantities typically increase with depth, and the variation can be $10-40\%$ at extreme depth values. Using Y3 as an example, we propagate the variable depth effect to the weak lensing $3\times2$pt analysis, and assess the impact on cosmological parameters via a Fisher forecast. We find that galaxy clustering is most susceptible to variable depth, and non-uniformity needs to be mitigated below $3\%$ to recover unbiased cosmological constraints. There is little impact on galaxy-shear and shear-shear power spectra, given the expected LSST Y3 noise., Comment: 28 pages, 23 figures. Matched to the accepted version on MNRAS

Outflowing gas from supermassive black holes in the centers of active galaxies has been postulated as a major contributor to galactic evolution. To explore the interaction between narrow-line region (NLR) outflows and their host galaxies, we use Hubble Space Telescope (HST) Space Telescope Imaging Spectrograph (STIS) spectra and Wide Field Camera 3 (WFC3) images of 15 nearby (z < 0.02) active galactic nuclei (AGN) to determine the extents and geometries of their NLRs. We combine new HST WFC3 continuum and [O III] $\lambda$5007A images of 11 AGN with 4 archival AGN to match existing spectra from HST STIS. For the 6 AGN with suitable long-slit coverage of their NLRs, we use isophotal fitting of ground-based images, continuum-subtracted [O III] images, and the STIS spectra, to resolve, measure, and de-project the gas kinematics to the plane of the host galaxy disk and distinguish NLR outflows from galaxy rotation and/or kinematically disturbed gas. We find an average [O III] extent of $\sim$680pc with a correlation between gas extent and [O III] luminosity of R$_\mathrm{[O III]}$ $\propto$ L$_{\text{[O III]}}^{0.39}$. The measured extents depend strongly on the depth of the [O III] images, highlighting the importance of adopting uniform thresholds when analyzing scaling relationships. The outflows reach from 39-88% of the full NLR extents, and we find that all 6 of the AGN with STIS coverage of their entire NLRs show strong kinematic evidence for outflows, despite previous uncertainty for these AGN. This suggests that NLR outflows are ubiquitous in moderate luminosity AGN and that standard criteria for kinematic modeling are essential for identifying outflows., Comment: Accepted for publication in ApJ on July 9, 2024. The paper has 21 pages, 10 figures, and 6 tables. Version two includes minor corrections to match the journal publication. The calibrated data are available through MAST at: https://archive.stsci.edu/hlsp/nlr-agn

The Vera C. Rubin Observatory is a unique facility for survey astronomy that will soon be commissioned and begin operations. Crucial to many of its scientific goals is the achievement of sustained high image quality, limited only by the seeing at the site. This will be maintained through an Active Optics System (AOS) that controls optical element misalignments and corrects mirror figure error to minimize aberrations caused by both thermal and gravitational distortions. However, the large number of adjustment degrees of freedom available on the Rubin Observatory introduces a range of degeneracies, including many that are \textit{noise-induced} due to imperfect measurement of the wavefront errors. We present a structured methodology for identifying these degeneracies through an analysis of image noise level. We also present a novel scaling strategy based on Truncated Singular Value Decomposition (TSVD) that mitigates the degeneracy, and optimally distributes the adjustment over the available degrees of freedom. Our approach ensures the attainment of optimal image quality, while avoiding excursions around the noise-induced subspace of degeneracies, marking a significant improvement over the previous techniques adopted for Rubin, which were based on an Optimal Integral Controller (OIC). This new approach is likely to also yield significant benefits for all telescopes that incorporate large numbers of degrees of freedom of adjustment., Comment: 19 pages, 8 figures, 4 tables. Submitted to ApJ

We have characterized the ionized, neutral, and warm molecular gas kinematics in the Seyfert 1 galaxy NGC 3227 using observations from the Hubble Space Telescope Space Telescope Imaging Spectrograph, Apache Point Observatory's Kitt Peak Ohio State Multi-Object Spectrograph, Gemini-North's Near-Infrared Integral Field Spectrometer, and the Atacama Large Millimeter Array. We fit multiple Gaussians to several spatially-resolved emission lines observed with long-slit and integral-field spectroscopy and isolate the kinematics based on apparent rotational and outflowing motions. We use the kinematics to determine an orientation for the bicone along which the outflows travel, and find that the biconical structure has an inclination of $40 ^{+5}_{-4}${\deg} from our line of sight, and a half-opening angle with an inner and outer boundary of $47 ^{+6}_{-2}${\deg} and $68 ^{+1}_{-1}${\deg}, respectively. We observe ionized outflows traveling 500 km s$^{-1}$ at distances up to 7$''$ (800 pc) from the SMBH, and disturbed ionized gas up to a distance of 15$''$ (1.7 kpc). Our analysis reveals that the ionized outflows are launched from within 20 pc of the SMBH, at the same location as a bridge of cold gas across the nucleus detected in ALMA CO(2-1) observations. We measure a turnover radius where the gas starts decelerating at a distance of $26 \pm 6$ pc from the AGN. Compared to a turnover radius in the range of $31- 63$ pc from a radiative driving model, we confirm that radiative driving is the dominant acceleration mechanism for the narrow line region (NLR) outflows in NGC 3227., Comment: 26 pages, 17 figures

Evaluating the accuracy and calibration of the redshift posteriors produced by photometric redshift (photo-z) estimators is vital for enabling precision cosmology and extragalactic astrophysics with modern wide-field photometric surveys. Evaluating photo-z posteriors on a per-galaxy basis is difficult, however, as real galaxies have a true redshift but not a true redshift posterior. We introduce PZFlow, a Python package for the probabilistic forward modeling of galaxy catalogs with normalizing flows. For catalogs simulated with PZFlow, there is a natural notion of "true" redshift posteriors that can be used for photo-z validation. We use PZFlow to simulate a photometric galaxy catalog where each galaxy has a redshift, noisy photometry, shape information, and a true redshift posterior. We also demonstrate the use of an ensemble of normalizing flows for photo-z estimation. We discuss how PZFlow will be used to validate the photo-z estimation pipeline of the Dark Energy Science Collaboration (DESC), and the wider applicability of PZFlow for statistical modeling of any tabular data., Comment: 19 pages, 13 figures, submitted to AJ

Large imaging surveys will rely on photometric redshifts (photo-z's), which are typically estimated through machine learning methods. Currently planned spectroscopic surveys will not be deep enough to produce a representative training sample for LSST, so we seek methods to improve the photo-z estimates that arise from non-representative training samples. Spectroscopic training samples for photo-z's are biased towards redder, brighter galaxies, which also tend to be at lower redshift than the typical galaxy observed by LSST, leading to poor photo-z estimates with outlier fractions nearly 4 times larger than for a representative training sample. In this paper, we apply the concept of training sample augmentation, where we augment simulated non-representative training samples with simulated galaxies possessing otherwise unrepresented features. When we select simulated galaxies with (g-z) color, i-band magnitude and redshift outside the range of the original training sample, we are able to reduce the outlier fraction of the photo-z estimates for simulated LSST data by nearly 50% and the normalized median absolute deviation (NMAD) by 56%. When compared to a fully representative training sample, augmentation can recover nearly 70% of the degradation in the outlier fraction and 80% of the degradation in NMAD. Training sample augmentation is a simple and effective way to improve training samples for photo-z's without requiring additional spectroscopic samples., Comment: 11 pages, 4 figures, published in ApJ Letters

The Vera C. Rubin Observatory will, over a period of 10 years, repeatedly survey the southern sky. To ensure that images generated by Rubin meet the quality requirements for precision science, the observatory will use an Active Optics System (AOS) to correct for alignment and mirror surface perturbations introduced by gravity and temperature gradients in the optical system. To accomplish this Rubin will use out-of-focus images from sensors located at the edge of the focal plane to learn and correct for perturbations to the wavefront. We have designed and integrated a deep learning model for wavefront estimation into the AOS pipeline. In this paper, we compare the performance of this deep learning approach to Rubin's baseline algorithm when applied to images from two different simulations of the Rubin optical system. We show the deep learning approach is faster and more accurate, achieving the atmospheric error floor both for high-quality images, and low-quality images with heavy blending and vignetting. Compared to the baseline algorithm, the deep learning model is 40x faster, the median error 2x better under ideal conditions, 5x better in the presence of vignetting by the Rubin camera, and 14x better in the presence of blending in crowded fields. In addition, the deep learning model surpasses the required optical quality in simulations of the AOS closed loop. This system promises to increase the survey area useful for precision science by up to 8%. We discuss how this system might be deployed when commissioning and operating Rubin., Comment: 24 pages, 21 figures

This study investigated multi-void formation in fused silica using high resolution Finite-Difference-Time-Domain (FDTD) simulations. Despite extensive research dedicated to understanding the mechanisms behind multi-void formation in materials, the fundamental aspects and mechanisms governing self-void array formation in dielectrics and polymers remain poorly understood. By modeling the voids as concentric spheres with densified shells and simulating the laser interaction with the voids, we showed that void array generation in fused silica is a linear mechanism. This study provides valuable insight into the mechanism behind the formation of void arrays in fused silica., Comment: arXiv admin note: substantial text overlap with arXiv:2305.02976

We present recent Very Long Baseline Array (VLBA) 5 GHz radio observations of the nearby, luminous Seyfert 2 galaxy NGC 1068 for comparison to similar VLBA observations made on 1997 April 26. By cross-correlating the positions of emitting regions across both epochs, we find that spatially-resolved extra-nuclear radio knots in this system have sub-relativistic transverse speeds (v < 0.1c). We discuss sources of the observed knots and how the radio emission relates to additional phases of gas in the central ~150 pcs of this system. We suggest that the most likely explanation for the observed emission is synchrotron radiation formed by shocked host media via interactions between AGN winds and the host environment., Comment: 13 pages, 4 figures, accepted for publication in ApJ

Let $G$ be a graph such that, whenever two vertices $x$ and $y$ of $G$ are joined by three internally disjoint paths, $x$ and $y$ are adjacent. Jamison and Mulder determined that the set of such graphs coincides with the set of graphs that can be built from cycles and complete graphs via 1-sums and parallel connections. This paper proves an analogous result for binary matroids., Comment: 12 pages, 6 figures

We present the hot molecular and warm ionised gas kinematics for 33 nearby ($0.001\lesssim z\lesssim0.056$) X-ray selected active galaxies using the H$_2 2.1218 \mu$m and Br$\gamma$ emission lines observed in the K-band with the Gemini Near-Infrared Field Spectrograph (NIFS). The observations cover the inner 0.04$-$2 kpc of each AGN at spatial resolutions of 4$-$250 pc with a velocity resolution of $\sigma_{\rm inst}\approx$20 ${\rm km s^{-1}}$. We find that 31 objects (94 per cent) present a kinematically disturbed region (KDR) seen in ionised gas, while such regions are observed in hot molecular gas for 25 galaxies (76 per cent). We interpret the KDR as being due to outflows with masses of 10$^2-$10$^7$ M$_\odot$ and 10$^0-$10$^4$ M$_\odot$ for the ionised and hot molecular gas, respectively. The ranges of mass-outflow rates ($\dot{M}_{\rm out}$) and kinetic power ($\dot{E}_{\rm K}$) of the outflows are 10$^{-3}-$10$^{1}$ M$_\odot$yr$^{-1}$ and $\sim$10$^{37}$$-$10$^{43}$ erg s$^{-1}$ for the ionised gas outflows, and 10$^{-5}$$-$10$^{-2}$ M$_\odot$ yr$^{-1}$ and 10$^{35}$$-$10$^{39}$ erg s$^{-1}$ for the hot molecular gas outflows. The median coupling efficiency in our sample is $\dot{E}_{K}/L_{\rm bol}\approx1.8\times10^{-3}$ and the estimated momentum fluxes of the outflows suggest they are produced by radiation-pressure in low-density environment, with possible contribution from shocks., Comment: 37 pages, published in MNRAS - A few typos in the text and in the label of Fg 1 were corrected in this version

We present dynamical models of the narrow line region (NLR) outflows in the nearby Seyfert galaxies Mrk 3, Mrk 78, NGC 1068, and NGC 4151 using observations from the Hubble Space Telescope and Apache Point Observatory. We employ long-slit spectroscopy to map the spatially-resolved outflow and rotational velocities of the ionized gas. We also perform surface brightness decompositions of host galaxy images to constrain the enclosed stellar mass distributions as functions of distance from the supermassive black holes (SMBHs). Assuming that the NLR gas is accelerated by AGN radiation pressure, and subsequently decelerated by the host galaxy and SMBH gravitational potentials, we derive outflow velocity profiles where the gas is launched in situ at multiple distances from the SMBH. We find a strong correlation between the turnover (from acceleration to deceleration) radii from our models, with the turnovers seen in the observed velocities and spatially-resolved mass outflow rates for the AGN with bolometric luminosities $>$ 10$^{44}$ erg sec$^{-1}$. This consistency indicates that radiation pressure is the dominant driving mechanism behind the NLR outflows in these moderate-luminosity AGN, with a force multiplier $\sim$500 yielding the best agreement between the modeled and observed turnover radii. However, in Meena2021 we found that this trend may not hold at lower luminosities, where our modeled turnover distance for NGC 4051 is much smaller than in the observed kinematics. This result may indicate that either additional force(s) are responsible for accelerating the NLR outflows in low-luminosity AGN, or higher spatial resolution observations are required to quantify their turnover radii., Comment: 36 pages, 10 Figures, 5 Tables, 7 supplementary Figures, Accepted for publication in ApJ

The form of the energy-momentum tensor when a quasimonochromatic field propagates into and through an antireflection-coated, sourceless, transparent, continuous, linear magneto-dielectric medium, initially at rest in the local frame, remains controversial. The Minkowski energy-momentum tensor is the main component of the electromagnetic conservation law. It has been known for over a century that the electromagnetic conservation law is unsound as evidenced by alternative energy-momentum tensors that have been proposed to ameliorate known physical deficiencies (violation of conservation of angular and linear momentum) and by the various material energy-momentum tensors and coupling forces that have been introduced to repair or complete the law. The extant resolution is to treat the continuum electromagnetic system as a subsystem and add a phenomenological material subsystem energy-momentum tensor. We show that the four-divergence of the total, electromagnetic plus material, energy-momentum tensor produces an energy continuity theorem in which the two non-zero terms depend on different powers of the refractive index $n$. Then the extant resolution of the Abraham-Minkowski controversy is self-inconsistent.

We present $Chandra$ ACIS imaging spectroscopy of the nucleus of the Seyfert 2 Galaxy Mrk 34. We identify spatially and spectrally resolved features in the band that includes Fe K$\alpha$, Fe XXV and Fe XXVI. These features indicate high-velocity ($\gtrsim15,000\,\rm{km\,s}^{-1}$ line-of-sight) material separated spanning $\sim$0.5 arcsec, within $\sim200$ pc of the nucleus. This outflow could have deprojected velocities $\sim12-28\times$ greater than the [O III] emitting outflows, and could potentially dominate the kinetic power in the outflow. This emission may point to the origins of the optical and X-ray winds observed at larger radii, and could indicate a link between ultra-fast outflows and AGN feedback on $\gtrsim$kpc scales., Comment: 13 pages, 7 figures. Resubmitted to ApJ with revisions

As we observe a rapidly growing number of astrophysical transients, we learn more about the diverse host galaxy environments in which they occur. Host galaxy information can be used to purify samples of cosmological Type Ia supernovae, uncover the progenitor systems of individual classes, and facilitate low-latency follow-up of rare and peculiar explosions. In this work, we develop a novel data-driven methodology to simulate the time-domain sky that includes detailed modeling of the probability density function for multiple transient classes conditioned on host galaxy magnitudes, colours, star formation rates, and masses. We have designed these simulations to optimize photometric classification and analysis in upcoming large synoptic surveys. We integrate host galaxy information into the SNANA simulation framework to construct the Simulated Catalogue of Optical Transients and Correlated Hosts (SCOTCH), a publicly-available catalogue of 5 million idealized transient light curves in LSST passbands and their host galaxy properties over the redshift range $0

The cycles of a graph give a natural cyclic ordering to their edge-sets, and these orderings are consistent in that two edges are adjacent in one cycle if and only if they are adjacent in every cycle in which they appear together. An orderable matroid is one whose set of circuits admits such a consistent ordering. In this paper, we consider the question of determining which matroids are orderable. Although we are able to answer this question for non-binary matroids, it remains open for binary matroids. We give examples to provide insight into the potential difficulty of this question in general. We also show that, by requiring that the ordering preserves the three arcs in every theta-graph restriction of a binary matroid $M$, we guarantee that $M$ is orderable if and only if $M$ is graphic., Comment: 28 pages, 19 figures, this is a corrected version

Active galactic nuclei (AGN) can launch outflows of ionized gas that may influence galaxy evolution, and quantifying their full impact requires spatially resolved measurements of the gas masses, velocities, and radial extents. We previously reported these quantities for the ionized narrow-line region (NLR) outflows in six low-redshift AGN, where the gas velocities and extents were determined from Hubble Space Telescope long-slit spectroscopy. However, calculating the gas masses required multi-component photoionization models to account for radial variations in the gas densities, which span $\sim$6 orders of magnitude. In order to simplify this method for larger samples with less spectral coverage, we compare these gas masses with those calculated from techniques in the literature. First, we use a recombination equation with three different estimates for the radial density profiles. These include constant densities, those derived from [S II], and power-law profiles based on constant values of the ionization parameter ($U$). Second, we use single-component photoionization models with power-law density profiles based on constant $U$, and allow $U$ to vary with radius based on the [O III]/H$\beta$ ratios. We find that assuming a constant density of $n_\mathrm{H} =$ 10$^2$ cm$^{-3}$ overestimates the gas masses for all six outflows, particularly at small radii where the outflow rates peak. The use of [S II] marginally matches the total gas masses, but also overestimates at small radii. Overall, single-component photoionization models where $U$ varies with radius are able to best match the gas mass and outflow rate profiles when there are insufficient emission lines to construct detailed models., Comment: Accepted for Publication in ApJ on March 11, 2022. The paper has 20 pages and 8 figures, with results tabulated in the Appendix. Version two includes minor corrections to match the journal publication

The accurate estimation of photometric redshifts is crucial to many upcoming galaxy surveys, for example the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). Almost all Rubin extragalactic and cosmological science requires accurate and precise calculation of photometric redshifts; many diverse approaches to this problem are currently in the process of being developed, validated, and tested. In this work, we use the photometric redshift code GPz to examine two realistically complex training set imperfections scenarios for machine learning based photometric redshift calculation: i) where the spectroscopic training set has a very different distribution in colour-magnitude space to the test set, and ii) where the effect of emission line confusion causes a fraction of the training spectroscopic sample to not have the true redshift. By evaluating the sensitivity of GPz to a range of increasingly severe imperfections, with a range of metrics (both of photo-z point estimates as well as posterior probability distribution functions, PDFs), we quantify the degree to which predictions get worse with higher degrees of degradation. In particular we find that there is a substantial drop-off in photo-z quality when line-confusion goes above ~1%, and sample incompleteness below a redshift of 1.5, for an experimental setup using data from the Buzzard Flock synthetic sky catalogues., Comment: 12 pages, 8 figures, accepted in PASP

We use Cloudy photoionisation models to predict the flux profiles for optical/IR emission lines that trace the footprint of X-ray gas, such as [Fe X] 6375A and [Si X] 1.43$\mu$m. These are a subset of coronal lines, from ions with ionisation potential $\geq$ that of O VII, i.e., 138eV. The footprint lines are formed in gas over the same range in ionisation state as the H and He-like of O and Ne ions, which are also the source of X-ray emission lines. The footprint lines can be detected with optical and IR telescopes, such as the Hubble Space Telescope/STIS and James Webb Space Telescope/NIRSpec, and can potentially be used to measure the kinematics of the extended X-ray emission gas. As a test case, we use the footprints to quantify the properties of the X-ray outflow in the Seyfert 1 galaxy NGC 4151. To confirm the accuracy of our method, we compare our model predictions to the measured flux from archival STIS spectra and previous ground-based studies, and the results are in good agreement. We also use our X-ray footprint method to predict the mass profile for the X-ray emission-line gas in NGC 4151 and derive a total spatially-integrated X-ray mass of $7.8(\pm 2.1) \times 10^{5}~M_{\odot}$, in comparison to $5.4(\pm 1.1) \times 10^{5}~M_{\odot}$ measured from a Chandra X-ray analysis. Our results indicate that high-ionisation footprint emission lines in the optical and near-infrared can be used to accurately trace the kinematics and physical conditions of AGN ionised, X-ray emission-line gas.

The cogirth, $g^\ast(M)$, of a matroid $M$ is the size of a smallest cocircuit of $M$. Finding the cogirth of a graphic matroid can be done in polynomial time, but Vardy showed in 1997 that it is NP-hard to find the cogirth of a binary matroid. In this paper, we show that $g^\ast(M)\leq \frac{1}{2}\vert E(M)\vert$ when $M$ is binary, unless $M$ simplifies to a projective geometry. We also show that, when equality holds, $M$ simplifies to a Bose-Burton geometry, that is, a matroid of the form $PG(r-1,2)-PG(k-1,2)$. These results extend to matroids representable over arbitrary finite fields., Comment: 8 pages