Showing 1–50 of 66 results for author: Brewer, B J

Rotation of the Globular Cluster Population of the Dark Matter Deficient Galaxy NGC 1052-DF4: Implication for the Total Mass

Authors: Yuan, Li, Brendon J. Brewer, Geraint F. Lewis

Abstract: We explore the globular cluster population of NGC 1052-DF4, a dark matter deficient galaxy, using Bayesian inference to search for the presence of rotation. The existence of such a rotating component is relevant to the estimation of the mass of the galaxy, and therefore the question of whether NGC 1052-DF4 is truly deficient of dark matter, similar to NGC 1052-DF2 another galaxy in the same group.… ▽ More We explore the globular cluster population of NGC 1052-DF4, a dark matter deficient galaxy, using Bayesian inference to search for the presence of rotation. The existence of such a rotating component is relevant to the estimation of the mass of the galaxy, and therefore the question of whether NGC 1052-DF4 is truly deficient of dark matter, similar to NGC 1052-DF2 another galaxy in the same group. The rotational characteristics of seven globular clusters in NGC 1052-DF4 were investigated, finding that a non-rotating kinematic model has a higher Bayesian evidence than a rotating model, by a factor of approximately 2.5. In addition, we find that under the assumption of rotation, its amplitude must be small. This distinct lack of rotation strengthens the case that, based on its intrinsic velocity dispersion, NGC 1052-DF4 is a truly dark matter deficient galaxy. △ Less

Submitted 16 May, 2024; originally announced May 2024.

Comments: 9 pages 6 figures. Accepted for publication in PASA

Detection of the Cosmological Time Dilation of High Redshift Quasars

Authors: Geraint F. Lewis, Brendon J. Brewer

Abstract: A fundamental prediction of relativistic cosmologies is that, due to the expansion of space, observations of the distant cosmos should be time dilated and appear to run slower than events in the local universe. Whilst observations of cosmological supernovae unambiguously display the expected redshift-dependent time dilation, this has not been the case for other distant sources. Here we present the… ▽ More A fundamental prediction of relativistic cosmologies is that, due to the expansion of space, observations of the distant cosmos should be time dilated and appear to run slower than events in the local universe. Whilst observations of cosmological supernovae unambiguously display the expected redshift-dependent time dilation, this has not been the case for other distant sources. Here we present the identification of cosmic time dilation in a sample of 190 quasars monitored for over two decades in multiple wavebands by assessing various hypotheses through Bayesian analysis. This detection counters previous claims that observed quasar variability lacked the expected redshift-dependent time dilation. Hence, as well as demonstrating the claim that the lack of the redshift dependence of quasar variability represents a significant challenge to the standard cosmological model, this analysis further indicates that the properties of quasars are consistent with them being truly cosmologically distant sources. △ Less

Submitted 6 June, 2023; originally announced June 2023.

Comments: 18 pages, 3 figures. Author version of a manuscript under review at Nature Astronomy

Report number: GFL-001

What Does the Virial Coefficient of the \Hb Broad-Line Region Depend On?

Authors: Lizvette Villafaña, Peter R. Williams, Tommaso Treu, Brendon J. Brewer, Aaron J. Barth, Vivian U, Vardha N. Bennert, Hengxiao Guo, Misty C. Bentz, Gabriela Canalizo, Alexei V. Filippenko, Elinor Gates, Michael D. Joner, Matthew A. Malkan, Jong-Hak Woo, Bela Abolfathi, Thomas Bohn, K. Azalee Bostroem, Andrew Brandel, Thomas G. Brink, Sanyum Channa, Maren Cosens, Edward Donohue, Goni Halevi, Carol E. Hood , et al. (15 additional authors not shown)

Abstract: We combine our dynamical modeling black hole mass measurements from the Lick AGN Monitoring Project 2016 sample with measured cross-correlation time lags and line widths to recover individual scale factors, f, used in traditional reverberation mapping analyses. We extend our sample by including prior results from Code for AGN Reverberation and Modeling of Emission Lines (caramel) studies that have… ▽ More We combine our dynamical modeling black hole mass measurements from the Lick AGN Monitoring Project 2016 sample with measured cross-correlation time lags and line widths to recover individual scale factors, f, used in traditional reverberation mapping analyses. We extend our sample by including prior results from Code for AGN Reverberation and Modeling of Emission Lines (caramel) studies that have utilized our methods. Aiming to improve the precision of black hole mass estimates, as well as uncover any regularities in the behavior of the broad-line region (BLR), we search for correlations between f and other AGN/BLR parameters. We find (i) evidence for a correlation between the virial coefficient log10(fmean,σ) and black hole mass, (ii) marginal evidence for a similar correlation between log10(frms,σ) and black hole mass, (iii) marginal evidence for an anti-correlation of BLR disk thickness with log10(fmean,FWHM)and log10(frms,FWHM), and (iv) marginal evidence for an anti-correlation of inclination angle with log10(fmean,FWHM), log10(frms,σ), and log10(fmean,σ). Lastly, we find marginal evidence for a correlation between line-profile shape, when using the root-meansquare spectrum, log10(FWHM/σ)rms, and the virial coefficient, log10(frms,σ), and investigate how BLR properties might be related to line-profile shape using caramel models. △ Less

Submitted 13 April, 2023; originally announced April 2023.

Comments: 17 pages, 10 figures

Testing the Cosmological Principle with CatWISE Quasars: A Bayesian Analysis of the Number-Count Dipole

Authors: Lawrence Dam, Geraint F. Lewis, Brendon J. Brewer

Abstract: The Cosmological Principle, that the Universe is homogeneous and isotropic on sufficiently large scales, underpins the standard model of cosmology. However, a recent analysis of 1.36 million infrared-selected quasars has identified a significant tension in the amplitude of the number-count dipole compared to that derived from the CMB, thus challenging the Cosmological Principle. Here we present a… ▽ More The Cosmological Principle, that the Universe is homogeneous and isotropic on sufficiently large scales, underpins the standard model of cosmology. However, a recent analysis of 1.36 million infrared-selected quasars has identified a significant tension in the amplitude of the number-count dipole compared to that derived from the CMB, thus challenging the Cosmological Principle. Here we present a Bayesian analysis of the same quasar sample, testing various hypotheses using the Bayesian evidence. We find unambiguous evidence for the presence of a dipole in the distribution of quasars with a direction that is consistent with the dipole identified in the CMB. However, the amplitude of the dipole is found to be 2.7 times larger than that expected from the conventional kinematic explanation of the CMB dipole, with a statistical significance of $5.7σ$. To compare these results with theoretical expectations, we sharpen the $Λ$CDM predictions for the probability distribution of the amplitude, taking into account a number of observational and theoretical systematics. In particular, we show that the presence of the Galactic plane mask causes a considerable loss of dipole signal due to a leakage of power into higher multipoles, exacerbating the discrepancy in the amplitude. By contrast, we show using probabilistic arguments that the source evolution of quasars improves the discrepancy, but only mildly so. These results support the original findings of an anomalously large quasar dipole, independent of the statistical methodology used. △ Less

Submitted 7 August, 2023; v1 submitted 15 December, 2022; originally announced December 2022.

Comments: 15 pages, 7 figures; accepted for publication in MNRAS

Chemo-dynamical substructure in the M31 inner halo globular clusters: Further evidence for a recent accretion event

Authors: Geraint F. Lewis, Brendon J. Brewer, Dougal Mackey, Annette M. N. Ferguson, Yuan, Li, Tim Adams

Abstract: Based upon a metallicity selection, we identify a significant sub-population of the inner halo globular clusters in the Andromeda Galaxy which we name the Dulais Structure. It is distinguished as a co-rotating group of 10-20 globular clusters which appear to be kinematically distinct from, and on average more metal-poor than, the majority of the inner halo population. Intriguingly, the orbital axi… ▽ More Based upon a metallicity selection, we identify a significant sub-population of the inner halo globular clusters in the Andromeda Galaxy which we name the Dulais Structure. It is distinguished as a co-rotating group of 10-20 globular clusters which appear to be kinematically distinct from, and on average more metal-poor than, the majority of the inner halo population. Intriguingly, the orbital axis of this Dulais Structure is closely aligned with that of the younger accretion event recently identified using a sub-population of globular clusters in the outer halo of Andromeda, and this is strongly suggestive of a causal relationship between the two. If this connection is confirmed, a natural explanation for the kinematics of the globular clusters in the Dulais Structure is that they trace the accretion of a substantial progenitor (~10^11 Msun) into the halo of Andromeda during the last few billion years, that may have occurred as part of a larger group infall. △ Less

Submitted 14 November, 2022; originally announced November 2022.

Comments: 10 pages, 7 figures - accepted for publication in MNRAS

Report number: GFL-001

The Lick AGN Monitoring Project 2016: Dynamical Modeling of Velocity-Resolved H\b{eta} Lags in Luminous Seyfert Galaxies

Authors: Lizvette Villafaña, Peter R. Williams, Tommaso Treu, Brendon J. Brewer, Aaron J. Barth, Vivian U, Vardha N. Bennert, H. Alexander Vogler, Hengxiao Guo, Misty C. Bentz, Gabriela Canalizo, Alexei V. Filippenko, Elinor Gates, Frederick Hamann, Michael D. Joner, Matthew A. Malkan, Jong-Hak Woo, Bela Abolfathi, L. E. Abramson, Stephen F. Armen, Hyun-Jin Bae, Thomas Bohn, Benjamin D. Boizelle, Azalee Bostroem, Andrew Brandel , et al. (40 additional authors not shown)

Abstract: We have modeled the velocity-resolved reverberation response of the H\b{eta} broad emission line in nine Seyfert 1 galaxies from the Lick Active Galactic Nucleus (AGN) Monitioring Project 2016 sample, drawing inferences on the geometry and structure of the low-ionization broad-line region (BLR) and the mass of the central supermassive black hole. Overall, we find that the H\b{eta} BLR is generally… ▽ More We have modeled the velocity-resolved reverberation response of the H\b{eta} broad emission line in nine Seyfert 1 galaxies from the Lick Active Galactic Nucleus (AGN) Monitioring Project 2016 sample, drawing inferences on the geometry and structure of the low-ionization broad-line region (BLR) and the mass of the central supermassive black hole. Overall, we find that the H\b{eta} BLR is generally a thick disk viewed at low to moderate inclination angles. We combine our sample with prior studies and investigate line-profile shape dependence, such as log10(FWHM/σ), on BLR structure and kinematics and search for any BLR luminosity-dependent trends. We find marginal evidence for an anticorrelation between the profile shape of the broad H\b{eta} emission line and the Eddington ratio, when using the root-mean-square spectrum. However, we do not find any luminosity-dependent trends, and conclude that AGNs have diverse BLR structure and kinematics, consistent with the hypothesis of transient AGN/BLR conditions rather than systematic trends. △ Less

Submitted 28 March, 2022; originally announced March 2022.

Dynamical Modeling of the CIV Broad Line Region of the $z=2.805$ Multiply Imaged Quasar SDSS J2222+2745

Authors: Peter R. Williams, Tommaso Treu, Håkon Dahle, Stefano Valenti, Louis Abramson, Aaron J. Barth, Brendon J. Brewer, Karianne Dyrland, Michael Gladders, Keith Horne, Keren Sharon

Abstract: We present CIV BLR modeling results for the multiply imaged $z=2.805$ quasar SDSS J2222+2745. Using data covering a 5.3 year baseline after accounting for gravitational time delays, we find models that can reproduce the observed emission-line spectra and integrated CIV fluctuations. The models suggest a thick disk BLR that is inclined by $\sim$40 degrees to the observer's line of sight and with a… ▽ More We present CIV BLR modeling results for the multiply imaged $z=2.805$ quasar SDSS J2222+2745. Using data covering a 5.3 year baseline after accounting for gravitational time delays, we find models that can reproduce the observed emission-line spectra and integrated CIV fluctuations. The models suggest a thick disk BLR that is inclined by $\sim$40 degrees to the observer's line of sight and with a emissivity weighted median radius of $r_{\rm median} = 33.0^{+2.4}_{-2.1}$ light days. The kinematics are dominated by near-circular Keplerian motion with the remainder inflowing. The rest-frame lag one would measure from the models is $τ_{\rm median} = 36.4^{+1.8}_{-1.8}$ days, which is consistent with measurements based on cross-correlation. We show a possible geometry and transfer function based on the model fits and find that the model-produced velocity-resolved lags are consistent with those from cross-correlation. We measure a black hole mass of $\log_{10}(M_{\rm BH}/M_\odot) = 8.31^{+0.07}_{-0.06}$, which requires a scale factor of $\log_{10}(f_{{\rm mean},σ}) = 0.20^{+0.09}_{-0.07}$. △ Less

Submitted 19 March, 2021; originally announced March 2021.

Comments: 8 pages, 5 figures, submitted to ApJL

Space Telescope and Optical Reverberation Mapping Project. XII. Broad-Line Region Modeling of NGC 5548

Authors: P. R. Williams, A. Pancoast, T. Treu, B. J. Brewer, B. M. Peterson, A. J. Barth, M. A. Malkan, G. De Rosa, Keith Horne, G. A. Kriss, N. Arav, M. C. Bentz, E. M. Cackett, E. Dalla Bontà, M. Dehghanian, C. Done, G. J. Ferland, C. J. Grier, J. Kaastra, E. Kara, C. S. Kochanek, S. Mathur, M. Mehdipour, R. W. Pogge, D. Proga , et al. (133 additional authors not shown)

Abstract: We present geometric and dynamical modeling of the broad line region for the multi-wavelength reverberation mapping campaign focused on NGC 5548 in 2014. The dataset includes photometric and spectroscopic monitoring in the optical and ultraviolet, covering the H$β$, C IV, and Ly$α$ broad emission lines. We find an extended disk-like H$β$ BLR with a mixture of near-circular and outflowing gas traje… ▽ More We present geometric and dynamical modeling of the broad line region for the multi-wavelength reverberation mapping campaign focused on NGC 5548 in 2014. The dataset includes photometric and spectroscopic monitoring in the optical and ultraviolet, covering the H$β$, C IV, and Ly$α$ broad emission lines. We find an extended disk-like H$β$ BLR with a mixture of near-circular and outflowing gas trajectories, while the C IV and Ly$α$ BLRs are much less extended and resemble shell-like structures. There is clear radial structure in the BLR, with C IV and Ly$α$ emission arising at smaller radii than the H$β$ emission. Using the three lines, we make three independent black hole mass measurements, all of which are consistent. Combining these results gives a joint inference of $\log_{10}(M_{\rm BH}/M_\odot) = 7.64^{+0.21}_{-0.18}$. We examine the effect of using the $V$ band instead of the UV continuum light curve on the results and find a size difference that is consistent with the measured UV-optical time lag, but the other structural and kinematic parameters remain unchanged, suggesting that the $V$ band is a suitable proxy for the ionizing continuum when exploring the BLR structure and kinematics. Finally, we compare the H$β$ results to similar models of data obtained in 2008 when the AGN was at a lower luminosity state. We find that the size of the emitting region increased during this time period, but the geometry and black hole mass remain unchanged, which confirms that the BLR kinematics suitably gauge the gravitational field of the central black hole. △ Less

Submitted 1 October, 2020; originally announced October 2020.

Comments: 26 pages, 19 figures, 1 table, accepted for publication in ApJ

The Globular Cluster Population of NGC 1052-DF2: Evidence for Rotation

Authors: Geraint F. Lewis, Brendon J. Brewer, Zhen Wan

Abstract: Based upon the kinematics of ten globular clusters, it has recently been claimed that the ultra-diffuse galaxy, NCD 1052-DF2, lacks a significant quantity of dark matter. Dynamical analyses have generally assumed that this galaxy is pressure supported, with the relatively small velocity dispersion of the globular cluster population indicating the deficit of dark matter. However, the presence of a… ▽ More Based upon the kinematics of ten globular clusters, it has recently been claimed that the ultra-diffuse galaxy, NCD 1052-DF2, lacks a significant quantity of dark matter. Dynamical analyses have generally assumed that this galaxy is pressure supported, with the relatively small velocity dispersion of the globular cluster population indicating the deficit of dark matter. However, the presence of a significant rotation of the globular cluster population could substantially modify this conclusion. Here we present the discovery of such a signature of rotation in the kinematics of NGC 1052-DF2's globular clusters, with a velocity amplitude of $\sim12.44^{+4.40}_{-5.16}$ km/s, which, through Bayesian model comparison, represents a marginally better fit to the available kinematic data; note that this rotation is distinct from, and approximately perpendicular to, the recently identified rotation of the stellar component of NGC 1052-DF2. Assuming this truly represents an underlying rotation, it is shown that the determined mass depends upon the inclination of the rotational component and, with a moderate inclination, the resultant mass to light ratio can exceed $M/L\sim10$. △ Less

Submitted 7 October, 2019; originally announced October 2019.

Comments: 5 pages, 4 figures, accepted for publication in MNRAS Letters

Report number: GFL-001

Two major accretion epochs in M31 from two distinct populations of globular clusters

Authors: Dougal Mackey, Geraint F. Lewis, Brendon J. Brewer, Annette M. N. Ferguson, Jovan Veljanoski, Avon P. Huxor, Michelle L. M. Collins, Patrick Côté, Rodrigo A. Ibata, Mike J. Irwin, Nicolas Martin, Alan W. McConnachie, Jorge Peñarrubia, Nial Tanvir, Zhen Wan

Abstract: Large galaxies grow through the accumulation of dwarf galaxies. In principle it is possible to trace this growth history using the properties of a galaxy's stellar halo. Previous investigations of the galaxy M31 (Andromeda) have shown that outside a radius of 25 kpc the population of halo globular clusters is rotating in alignment with the stellar disk, as are more centrally located clusters. The… ▽ More Large galaxies grow through the accumulation of dwarf galaxies. In principle it is possible to trace this growth history using the properties of a galaxy's stellar halo. Previous investigations of the galaxy M31 (Andromeda) have shown that outside a radius of 25 kpc the population of halo globular clusters is rotating in alignment with the stellar disk, as are more centrally located clusters. The M31 halo also contains coherent stellar substructures, along with a smoothly distributed stellar component. Many of the globular clusters outside 25 kpc are associated with the most prominent substructures, while others are part of the smooth halo. Here we report a new analysis of the kinematics of these globular clusters. We find that the two distinct populations are rotating with perpendicular orientations. The rotation axis for the population associated with the smooth halo is aligned with the rotation axis for the plane of dwarf galaxies that encircles M31. We interpret these separate cluster populations as arising from two major accretion epochs, likely separated by billions of years. Stellar substructures from the first epoch are gone, but those from the more recent second epoch still remain. △ Less

Submitted 2 October, 2019; originally announced October 2019.

Comments: Authors' version of a Letter published in Nature on October 3rd, 2019

The SAMI Galaxy Survey: Bayesian Inference for Gas Disk Kinematics using a Hierarchical Gaussian Mixture Model

Authors: Mathew R. Varidel, Scott M. Croom, Geraint F. Lewis, Brendon J. Brewer, Enrico M. Di Teodoro, Joss Bland-Hawthorn, Julia J. Bryant, Christoph Federrath, Caroline Foster, Karl Glazebrook, Michael Goodwin, Brent Groves, Andrew M. Hopkins, Jon S. Lawrence, Ángel R. López-Sánchez, Anne M. Medling, Matt S. Owers, Samuel N. Richards, Richard Scalzo, Nicholas Scott, Sarah M. Sweet, Dan S. Taranu, Jesse van de Sande

Abstract: We present a novel Bayesian method, referred to as Blobby3D, to infer gas kinematics that mitigates the effects of beam smearing for observations using Integral Field Spectroscopy (IFS). The method is robust for regularly rotating galaxies despite substructure in the gas distribution. Modelling the gas substructure within the disk is achieved by using a hierarchical Gaussian mixture model. To acco… ▽ More We present a novel Bayesian method, referred to as Blobby3D, to infer gas kinematics that mitigates the effects of beam smearing for observations using Integral Field Spectroscopy (IFS). The method is robust for regularly rotating galaxies despite substructure in the gas distribution. Modelling the gas substructure within the disk is achieved by using a hierarchical Gaussian mixture model. To account for beam smearing effects, we construct a modelled cube that is then convolved per wavelength slice by the seeing, before calculating the likelihood function. We show that our method can model complex gas substructure including clumps and spiral arms. We also show that kinematic asymmetries can be observed after beam smearing for regularly rotating galaxies with asymmetries only introduced in the spatial distribution of the gas. We present findings for our method applied to a sample of 20 star-forming galaxies from the SAMI Galaxy Survey. We estimate the global H$α$ gas velocity dispersion for our sample to be in the range $\barσ_v \sim $[7, 30] km s$^{-1}$. The relative difference between our approach and estimates using the single Gaussian component fits per spaxel is $Δ\barσ_v / \barσ_v = - 0.29 \pm 0.18$ for the H$α$ flux-weighted mean velocity dispersion. △ Less

Submitted 10 March, 2019; v1 submitted 7 March, 2019; originally announced March 2019.

Comments: 23 pages, 12 figures, accepted for MNRAS

The Lick AGN Monitoring Project 2011: Dynamical Modeling of the Broad-Line Region

Authors: Peter R. Williams, Anna Pancoast, Tommaso Treu, Brendon J. Brewer, Aaron J. Barth, Vardha N. Bennert, Tabitha Buehler, Gabriela Canalizo, S. Bradley Cenko, Kelsey I. Clubb, Michael C. Cooper, Alexei V. Filippenko, Elinor Gates, Sebastian F. Hoenig, Michael D. Joner, Michael T. Kandrashoff, Clifton David Laney, Mariana S. Lazarova, Weidong Li, Matthew A. Malkan, Jacob Rex, Jeffrey M. Silverman, Erik Tollerud, Jonelle L. Walsh, Jong-Hak Woo

Abstract: We present models of the H$β$-emitting broad-line region (BLR) in seven Seyfert 1 galaxies from the Lick AGN (Active Galactic Nucleus) Monitoring Project 2011 sample, drawing inferences on the BLR structure and dynamics as well as the mass of the central supermassive black hole. We find that the BLR is generally a thick disk, viewed close to face-on, with preferential emission back toward the ioni… ▽ More We present models of the H$β$-emitting broad-line region (BLR) in seven Seyfert 1 galaxies from the Lick AGN (Active Galactic Nucleus) Monitoring Project 2011 sample, drawing inferences on the BLR structure and dynamics as well as the mass of the central supermassive black hole. We find that the BLR is generally a thick disk, viewed close to face-on, with preferential emission back toward the ionizing source. The dynamics in our sample range from near-circular elliptical orbits to inflowing or outflowing trajectories. We measure black hole masses of $\log_{10}(M_{\rm BH}/M_\odot) = 6.48^{+0.21}_{-0.18}$ for PG 1310$-$108, $7.50^{+0.25}_{-0.18}$ for Mrk 50, $7.46^{+0.15}_{-0.21}$ for Mrk 141, $7.58^{+0.08}_{-0.08}$ for Mrk 279, $7.11^{+0.20}_{-0.17}$ for Mrk 1511, $6.65^{+0.27}_{-0.15}$ for NGC 4593, and $6.94^{+0.14}_{-0.14}$ for Zw 229$-$015. We use these black hole mass measurements along with cross-correlation time lags and line widths to recover the scale factor $f$ used in traditional reverberation mapping measurements. Combining our results with other studies that use this modeling technique, bringing our sample size to 16, we calculate a scale factor that can be used for measuring black hole masses in other reverberation mapping campaigns. When using the root-mean-square (rms) spectrum and using the line dispersion to measure the line width, we find $\log_{10}(f_{{\rm rms},σ})_{\rm pred} = 0.57 \pm 0.19$. Finally, we search for correlations between $f$ and other AGN and BLR parameters and find marginal evidence that $f$ is correlated with $M_{\rm BH}$ and the BLR inclination angle, but no significant evidence of a correlation with the AGN luminosity or Eddington ratio. △ Less

Submitted 13 September, 2018; originally announced September 2018.

Comments: 26 pages, 14 figures. Accepted for publication in ApJ

kima: Exoplanet detection in radial velocities

Authors: J. P. Faria, N. C. Santos, P. Figueira, B. J. Brewer

Abstract: The radial-velocity (RV) method is one of the most successful in the detection of exoplanets, but is hindered by the intrinsic RV variations of the star, which can easily mimic or hide true planetary signals. kima is a package for the detection and characterization of exoplanets using RV data. It fits a sum of Keplerian curves to a timeseries of RV measurements and calculates the evidence for mode… ▽ More The radial-velocity (RV) method is one of the most successful in the detection of exoplanets, but is hindered by the intrinsic RV variations of the star, which can easily mimic or hide true planetary signals. kima is a package for the detection and characterization of exoplanets using RV data. It fits a sum of Keplerian curves to a timeseries of RV measurements and calculates the evidence for models with a fixed number Np of Keplerian signals, or after marginalising over Np. Moreover, kima can use a GP with a quasi-periodic kernel as a noise model, to deal with activity-induced signals. The hyperparameters of the GP are inferred together with the orbital parameters. The code is written in C++, but includes a helper Python package, pykima, which facilitates the analysis of the results. △ Less

Submitted 21 June, 2018; originally announced June 2018.

Comments: 3 pages; published in the Journal of Open Source Software; all code available at https://github.com/j-faria/kima

Journal ref: Journal of Open Source Software, 3(26), 487, 2018

Stability of the Broad Line Region Geometry and Dynamics in Arp 151 Over Seven Years

Authors: A. Pancoast, A. J. Barth, K. Horne, T. Treu, B. J. Brewer, V. N. Bennert, G. Canalizo, E. L. Gates, W. Li, M. A. Malkan, D. Sand, T. Schmidt, S. Valenti, J. H. Woo, K. I. Clubb, M. C. Cooper, S. M. Crawford, S. F. Honig, M. D. Joner, M. T. Kandrashoff, M. Lazarova, A. M. Nierenberg, E. Romero-Colmenero, D. Son, E. Tollerud , et al. (2 additional authors not shown)

Abstract: The Seyfert 1 galaxy Arp 151 was monitored as part of three reverberation mapping campaigns spanning $2008-2015$. We present modeling of these velocity-resolved reverberation mapping datasets using a geometric and dynamical model for the broad line region (BLR). By modeling each of the three datasets independently, we infer the evolution of the BLR structure in Arp 151 over a total of seven years… ▽ More The Seyfert 1 galaxy Arp 151 was monitored as part of three reverberation mapping campaigns spanning $2008-2015$. We present modeling of these velocity-resolved reverberation mapping datasets using a geometric and dynamical model for the broad line region (BLR). By modeling each of the three datasets independently, we infer the evolution of the BLR structure in Arp 151 over a total of seven years and constrain the systematic uncertainties in non-varying parameters such as the black hole mass. We find that the BLR geometry of a thick disk viewed close to face-on is stable over this time, although the size of the BLR grows by a factor of $\sim 2$. The dynamics of the BLR are dominated by inflow and the inferred black hole mass is consistent for the three datasets, despite the increase in BLR size. Combining the inference for the three datasets yields a black hole mass and statistical uncertainty of $\log_{10}($M$_{\rm BH}/\rm{M}_{\odot})=6.82^{+0.09}_{-0.09}$ with a standard deviation in individual measurements of 0.13 dex. △ Less

Submitted 6 March, 2018; originally announced March 2018.

Comments: 21 pages, 8 figures, 2 tables. Accepted for publication in ApJ

Computing Entropies With Nested Sampling

Authors: Brendon J. Brewer

Abstract: The Shannon entropy, and related quantities such as mutual information, can be used to quantify uncertainty and relevance. However, in practice, it can be difficult to compute these quantities for arbitrary probability distributions, particularly if the probability mass functions or densities cannot be evaluated. This paper introduces a computational approach, based on Nested Sampling, to evaluate… ▽ More The Shannon entropy, and related quantities such as mutual information, can be used to quantify uncertainty and relevance. However, in practice, it can be difficult to compute these quantities for arbitrary probability distributions, particularly if the probability mass functions or densities cannot be evaluated. This paper introduces a computational approach, based on Nested Sampling, to evaluate entropies of probability distributions that can only be sampled. I demonstrate the method on three examples: a simple gaussian example where the key quantities are available analytically; (ii) an experimental design example about scheduling observations in order to measure the period of an oscillating signal; and (iii) predicting the future from the past in a heavy-tailed scenario. △ Less

Submitted 15 August, 2017; v1 submitted 12 July, 2017; originally announced July 2017.

Comments: Accepted for publication in Entropy. 21 pages, 3 figures. Software available at https://github.com/eggplantbren/InfoNest

The Structure of the Broad-Line Region In Active Galactic Nuclei. II. Dynamical Modeling of Data from the AGN10 Reverberation Mapping Campaign

Authors: C. J. Grier, A. Pancoast, A. J. Barth, M. M. Fausnaugh, B. J. Brewer, T. Treu, B. M. Peterson

Abstract: We present inferences on the geometry and kinematics of the broad-Hbeta line-emitting region in four active galactic nuclei monitored as a part of the fall 2010 reverberation mapping campaign at MDM Observatory led by the Ohio State University. From modeling the continuum variability and response in emission-line profile changes as a function of time, we infer the geometry of the Hbeta- emitting b… ▽ More We present inferences on the geometry and kinematics of the broad-Hbeta line-emitting region in four active galactic nuclei monitored as a part of the fall 2010 reverberation mapping campaign at MDM Observatory led by the Ohio State University. From modeling the continuum variability and response in emission-line profile changes as a function of time, we infer the geometry of the Hbeta- emitting broad line regions to be thick disks that are close to face-on to the observer with kinematics that are well-described by either elliptical orbits or inflowing gas. We measure the black hole mass to be log (MBH) = 7.25 (+/-0.10) for Mrk 335, 7.86 (+0.20, -0.17) for Mrk 1501, 7.84 (+0.14, -0.19) for 3C 120, and 6.92 (+0.24, -0.23) for PG 2130+099. These black hole mass measurements are not based on a particular assumed value of the virial scale factor f, allowing us to compute individual f factors for each target. Our results nearly double the number of targets that have been modeled in this manner, and investigate the properties of a more diverse sample by including previously modeled objects. We measure an average scale factor f in the entire sample to be log10(f) = 0.54 +/- 0.17 when the line dispersion is used to characterize the line width, which is consistent with values derived using the normalization of the MBH-sigma relation. We find that the scale factor f for individual targets is likely correlated with the black hole mass, inclination angle, and opening angle of the broad line region but we do not find any correlation with the luminosity. △ Less

Submitted 27 September, 2017; v1 submitted 5 May, 2017; originally announced May 2017.

Comments: 19 pages, 13 Figures. Accepted for publication in ApJ

Red nuggets grow inside-out: evidence from gravitational lensing

Authors: Lindsay Oldham, Matt Auger, Chris Fassnacht, Tommaso Treu, Brendon J. Brewer, L. V. E. Koopmans, David Lagattuta, Philip Marshall, John McKean, Simona Vegetti

Abstract: We present a new sample of strong gravitational lens systems where both the foreground lenses and background sources are early-type galaxies. Using imaging from HST/ACS and Keck/NIRC2, we model the surface brightness distributions and show that the sources form a distinct population of massive, compact galaxies at redshifts $0.4 \lesssim z \lesssim 0.7$, lying systematically below the size-mass re… ▽ More We present a new sample of strong gravitational lens systems where both the foreground lenses and background sources are early-type galaxies. Using imaging from HST/ACS and Keck/NIRC2, we model the surface brightness distributions and show that the sources form a distinct population of massive, compact galaxies at redshifts $0.4 \lesssim z \lesssim 0.7$, lying systematically below the size-mass relation of the global elliptical galaxy population at those redshifts. These may therefore represent relics of high-redshift red nuggets or their partly-evolved descendants. We exploit the magnifying effect of lensing to investigate the structural properties, stellar masses and stellar populations of these objects with a view to understanding their evolution. We model these objects parametrically and find that they generally require two Sérsic components to properly describe their light profiles, with one more spheroidal component alongside a more envelope-like component, which is slightly more extended though still compact. This is consistent with the hypothesis of the inside-out growth of these objects via minor mergers. We also find that the sources can be characterised by red-to-blue colour gradients as a function of radius which are stronger at low redshift -- indicative of ongoing accretion -- but that their environments generally appear consistent with that of the general elliptical galaxy population, contrary to recent suggestions that these objects are predominantly associated with clusters. △ Less

Submitted 31 October, 2016; originally announced November 2016.

Comments: 21 pages; accepted for publication in MNRAS

The elusive stellar halo of the Triangulum Galaxy

Authors: B. McMonigal, G. F. Lewis, B. J. Brewer, M. J. Irwin, N. F. Martin, A. W. McConnachie, R. A. Ibata, A. M. N. Ferguson, A. D. Mackey, S. C. Chapman

Abstract: The stellar halos of large galaxies represent a vital probe of the processes of galaxy evolution. They are the remnants of the initial bouts of star formation during the collapse of the proto-galactic cloud, coupled with imprint of ancient and on-going accretion events. Previously, we have reported the tentative detection of a possible, faint, extended stellar halo in the Local Group spiral, the T… ▽ More The stellar halos of large galaxies represent a vital probe of the processes of galaxy evolution. They are the remnants of the initial bouts of star formation during the collapse of the proto-galactic cloud, coupled with imprint of ancient and on-going accretion events. Previously, we have reported the tentative detection of a possible, faint, extended stellar halo in the Local Group spiral, the Triangulum Galaxy (M33). However, the presence of substructure surrounding M33 made interpretation of this feature difficult. Here, we employ the final data set from the Pan-Andromeda Archaeological Survey (PAndAS), combined with an improved calibration and a newly derived contamination model for the region to revisit this claim. With an array of new fitting algorithms, fully accounting for contamination and the substantial substructure beyond the prominent stellar disk in M33, we reanalyse the surrounds to separate the signal of the stellar halo and the outer halo substructure. Using more robust search algorithms, we do not detect a large scale smooth stellar halo and place a limit on the maximum surface brightness of such a feature of $μ_V$ = 35.5 mags per square arcsec, or a total halo luminosity of $L < 10^6L_{\odot}$. △ Less

Submitted 7 July, 2016; originally announced July 2016.

Comments: 17 pages, 23 figures, accepted for publication in MNRAS

Report number: GFL-001

DNest4: Diffusive Nested Sampling in C++ and Python

Authors: Brendon J. Brewer, Daniel Foreman-Mackey

Abstract: In probabilistic (Bayesian) inferences, we typically want to compute properties of the posterior distribution, describing knowledge of unknown quantities in the context of a particular dataset and the assumed prior information. The marginal likelihood, also known as the "evidence", is a key quantity in Bayesian model selection. The Diffusive Nested Sampling algorithm, a variant of Nested Sampling,… ▽ More In probabilistic (Bayesian) inferences, we typically want to compute properties of the posterior distribution, describing knowledge of unknown quantities in the context of a particular dataset and the assumed prior information. The marginal likelihood, also known as the "evidence", is a key quantity in Bayesian model selection. The Diffusive Nested Sampling algorithm, a variant of Nested Sampling, is a powerful tool for generating posterior samples and estimating marginal likelihoods. It is effective at solving complex problems including many where the posterior distribution is multimodal or has strong dependencies between variables. DNest4 is an open source (MIT licensed), multi-threaded implementation of this algorithm in C++11, along with associated utilities including: i) RJObject, a class template for finite mixture models, (ii) A Python package allowing basic use without C++ coding, and iii) Experimental support for models implemented in Julia. In this paper we demonstrate DNest4 usage through examples including simple Bayesian data analysis, finite mixture models, and Approximate Bayesian Computation. △ Less

Submitted 1 November, 2016; v1 submitted 12 June, 2016; originally announced June 2016.

Comments: Submitted. 33 pages, 9 figures. v2 removed a duplicated figure, v3 added a comparison to other packages, v4 fixed a few minor issues

Space Telescope and Optical Reverberation Mapping Project. IV. Anomalous behavior of the broad ultraviolet emission lines in NGC 5548

Authors: M. R. Goad, K. T. Korista, G. De Rosa, G. A. Kriss, R. Edelson, A. J. Barth, G. J. Ferland, C. S. Kochanek, H. Netzer, B. M. Peterson, M. C. Bentz, S. Bisogni, D. M. Crenshaw, K. D. Denney, J. Ely, M. M. Fausnaugh, C. J. Grier, A. Gupta, K. D. Horne, J. Kaastra, A. Pancoast, L. Pei, R. W. Pogge, A. Skielboe, D. Starkey , et al. (77 additional authors not shown)

Abstract: During an intensive Hubble Space Telescope (HST) Cosmic Origins Spectrograph (COS) UV monitoring campaign of the Seyfert~1 galaxy NGC 5548 performed from 2014 February to July, the normally highly correlated far-UV continuum and broad emission-line variations decorrelated for ~60 to 70 days, starting ~75 days after the first HST/COS observation. Following this anomalous state, the flux and variabi… ▽ More During an intensive Hubble Space Telescope (HST) Cosmic Origins Spectrograph (COS) UV monitoring campaign of the Seyfert~1 galaxy NGC 5548 performed from 2014 February to July, the normally highly correlated far-UV continuum and broad emission-line variations decorrelated for ~60 to 70 days, starting ~75 days after the first HST/COS observation. Following this anomalous state, the flux and variability of the broad emission lines returned to a more normal state. This transient behavior, characterised by significant deficits in flux and equivalent width of the strong broad UV emission lines, is the first of its kind to be unambiguously identified in an active galactic nucleus reverberation mapping campaign. The largest corresponding emission-line flux deficits occurred for the high-ionization collisionally excited lines, C IV and Si IV(+O IV]), and also He II(+O III]), while the anomaly in Ly-alpha was substantially smaller. This pattern of behavior indicates a depletion in the flux of photons with E_{\rm ph} > 54 eV, relative to those near 13.6 eV. We suggest two plausible mechanisms for the observed behavior: (i) temporary obscuration of the ionizing continuum incident upon BLR clouds by a moving veil of material lying between the inner accretion disk and inner BLR, perhaps resulting from an episodic ejection of material from the disk, or (ii) a temporary change in the intrinsic ionizing continuum spectral energy distribution resulting in a deficit of ionizing photons with energies > 54 eV, possibly due to a transient restructuring of the Comptonizing atmosphere above the disk. Current evidence appears to favor the latter explanation. △ Less

Submitted 29 March, 2016; originally announced March 2016.

Comments: 12 pages, 2 figures, accepted for publication in ApJ, Friday 25th March 2016. A movie of the anomalous emission-line behavior can be found in the ancilliary documents

Uncovering the planets and stellar activity of CoRoT-7 using only radial velocities

Authors: J. P. Faria, R. D. Haywood, B. J. Brewer, P. Figueira, M. Oshagh, A. Santerne, N. C. Santos

Abstract: Stellar activity can induce signals in the radial velocities of stars, complicating the detection of orbiting low-mass planets. We present a method to determine the number of planetary signals present in radial-velocity datasets of active stars, using only radial-velocity observations. Instead of considering separate fits with different number of planets, we use a birth-death Markov chain Monte Ca… ▽ More Stellar activity can induce signals in the radial velocities of stars, complicating the detection of orbiting low-mass planets. We present a method to determine the number of planetary signals present in radial-velocity datasets of active stars, using only radial-velocity observations. Instead of considering separate fits with different number of planets, we use a birth-death Markov chain Monte Carlo algorithm to infer the posterior distribution for the number of planets in a single run. In a natural way, the marginal distributions for the orbital parameters of all planets are also inferred. This method is applied to HARPS data of CoRoT-7. We confidently recover both CoRoT-7b and CoRoT-7c although the data show evidence for additional signals. △ Less

Submitted 3 February, 2016; v1 submitted 27 January, 2016; originally announced January 2016.

Comments: 7 pages, 5 figures, 2 tables, accepted for publication in A&A, language editing revisions

Journal ref: A&A 588, A31 (2016)

Space Telescope and Optical Reverberation Mapping Project. III. Optical Continuum Emission and Broad-Band Time Delays in NGC 5548

Authors: M. M. Fausnaugh, K. D. Denney, A. J. Barth, M. C. Bentz, M. C. Bottorff, M. T. Carini, K. V. Croxall, G. De Rosa, M. R. Goad, Keith Horne, M. D. Joner, S. Kaspi, M. Kim, S. A. Klimanov, C. S. Kochanek, D. C. Leonard, H. Netzer, B. M. Peterson, K. Schnulle, S. G. Sergeev, M. Vestergaard, W. -K. Zheng, Y. Zu, M. D. Anderson, P. Arevalo , et al. (72 additional authors not shown)

Abstract: We present ground-based optical photometric monitoring data for NGC 5548, part of an extended multi-wavelength reverberation mapping campaign. The light curves have nearly daily cadence from 2014 January to July in nine filters (\emph{BVRI} and \emph{ugriz}). Combined with ultraviolet data from the \emph{Hubble Space Telescope} and \emph{Swift}, we confirm significant time delays between the conti… ▽ More We present ground-based optical photometric monitoring data for NGC 5548, part of an extended multi-wavelength reverberation mapping campaign. The light curves have nearly daily cadence from 2014 January to July in nine filters (\emph{BVRI} and \emph{ugriz}). Combined with ultraviolet data from the \emph{Hubble Space Telescope} and \emph{Swift}, we confirm significant time delays between the continuum bands as a function of wavelength, extending the wavelength coverage from 1158\,Å to the $z$ band ($\sim\!9160$\,Å). We find that the lags at wavelengths longer than the {\it V} band are equal to or greater than the lags of high-ionization-state emission lines (such as He\,{\sc ii}\,$λ1640$ and $λ4686$), suggesting that the continuum-emitting source is of a physical size comparable to the inner broad-line region (BLR). The trend of lag with wavelength is broadly consistent with the prediction for continuum reprocessing by an accretion disk with $τ\propto λ^{4/3}$. However, the lags also imply a disk radius that is 3 times larger than the prediction from standard thin-disk theory, assuming that the bolometric luminosity is 10\% of the Eddington luminosity ($L = 0.1L_{\rm Edd}$). Using optical spectra from the Large Binocular Telescope, we estimate the bias of the interband continuum lags due to BLR emission observed in the filters. We find that the bias for filters with high levels of BLR contamination ($\sim\! 20\%$) can be important for the shortest continuum lags, and likely has a significant impact on the {\it u} and {\it U} bands owing to Balmer continuum emission. △ Less

Submitted 29 February, 2016; v1 submitted 19 October, 2015; originally announced October 2015.

Comments: 26 pages, 10 figures, accepted to ApJ. For a brief video describing the main results of this paper, please see: https://www.youtube.com/watch?v=yaYtcDvIoP0&feature=youtu.be

Trans-Dimensional Bayesian Inference for Gravitational Lens Substructures

Authors: Brendon J. Brewer, David Huijser, Geraint F. Lewis

Abstract: We introduce a Bayesian solution to the problem of inferring the density profile of strong gravitational lenses when the lens galaxy may contain multiple dark or faint substructures. The source and lens models are based on a superposition of an unknown number of non-negative basis functions (or "blobs") whose form was chosen with speed as a primary criterion. The prior distribution for the blobs'… ▽ More We introduce a Bayesian solution to the problem of inferring the density profile of strong gravitational lenses when the lens galaxy may contain multiple dark or faint substructures. The source and lens models are based on a superposition of an unknown number of non-negative basis functions (or "blobs") whose form was chosen with speed as a primary criterion. The prior distribution for the blobs' properties is specified hierarchically, so the mass function of substructures is a natural output of the method. We use reversible jump Markov Chain Monte Carlo (MCMC) within Diffusive Nested Sampling (DNS) to sample the posterior distribution and evaluate the marginal likelihood of the model, including the summation over the unknown number of blobs in the source and the lens. We demonstrate the method on two simulated data sets: one with a single substructure, and one with ten. We also apply the method to the g-band image of the "Cosmic Horseshoe" system, and find evidence for more than zero substructures. However, these have large spatial extent and probably only point to misspecifications in the model (such as the shape of the smooth lens component or the point spread function), which are difficult to guard against in full generality. △ Less

Submitted 9 October, 2015; v1 submitted 4 August, 2015; originally announced August 2015.

Comments: Accepted for publication in MNRAS. 12 pages, 15 figures

The Lick AGN Monitoring Project 2011: Spectroscopic Campaign and Emission-Line Light Curves

Authors: A. J. Barth, V. N. Bennert, G. Canalizo, A. V. Filippenko, E. L. Gates, J. E. Greene, W. Li, M. A. Malkan, A. Pancoast, D. J. Sand, D. Stern, T. Treu, J. -H. Woo, R. J. Assef, H. -J. Bae, B. J. Brewer, S. B. Cenko, K. I. Clubb, M. C. Cooper, A. M. Diamond-Stanic, K. D. Hiner, S. F. Hoenig, E. Hsiao, M. T. Kandrashoff, M. S. Lazarova , et al. (5 additional authors not shown)

Abstract: In the Spring of 2011 we carried out a 2.5 month reverberation mapping campaign using the 3 m Shane telescope at Lick Observatory, monitoring 15 low-redshift Seyfert 1 galaxies. This paper describes the observations, reductions and measurements, and data products from the spectroscopic campaign. The reduced spectra were fitted with a multicomponent model in order to isolate the contributions of va… ▽ More In the Spring of 2011 we carried out a 2.5 month reverberation mapping campaign using the 3 m Shane telescope at Lick Observatory, monitoring 15 low-redshift Seyfert 1 galaxies. This paper describes the observations, reductions and measurements, and data products from the spectroscopic campaign. The reduced spectra were fitted with a multicomponent model in order to isolate the contributions of various continuum and emission-line components. We present light curves of broad emission lines and the AGN continuum, and measurements of the broad H-beta line widths in mean and root-mean square (rms) spectra. For the most highly variable AGNs we also measured broad H-beta line widths and velocity centroids from the nightly spectra. In four AGNs exhibiting the highest variability amplitudes, we detect anticorrelations between broad H-beta width and luminosity, demonstrating that the broad-line region "breathes" on short timescales of days to weeks in response to continuum variations. We also find that broad H-beta velocity centroids can undergo substantial changes in response to continuum variations; in NGC 4593 the broad H-beta velocity shifted by ~250 km/s over a one-month duration. This reverberation-induced velocity shift effect is likely to contribute a significant source of confusion noise to binary black hole searches that use multi-epoch quasar spectroscopy to detect binary orbital motion. We also present results from simulations that examine biases that can occur in measurement of broad-line widths from rms spectra due to the contributions of continuum variations and photon-counting noise. △ Less

Submitted 3 March, 2015; originally announced March 2015.

Comments: 33 pages, 28 figures, accepted for publication in ApJ Supplement Series

Fast Bayesian Inference for Exoplanet Discovery in Radial Velocity Data

Authors: Brendon J. Brewer, Courtney P. Donovan

Abstract: Inferring the number of planets $N$ in an exoplanetary system from radial velocity (RV) data is a challenging task. Recently, it has become clear that RV data can contain periodic signals due to stellar activity, which can be difficult to distinguish from planetary signals. However, even doing the inference under a given set of simplifying assumptions (e.g. no stellar activity) can be difficult. I… ▽ More Inferring the number of planets $N$ in an exoplanetary system from radial velocity (RV) data is a challenging task. Recently, it has become clear that RV data can contain periodic signals due to stellar activity, which can be difficult to distinguish from planetary signals. However, even doing the inference under a given set of simplifying assumptions (e.g. no stellar activity) can be difficult. It is common for the posterior distribution for the planet parameters, such as orbital periods, to be multimodal and to have other awkward features. In addition, when $N$ is unknown, the marginal likelihood (or evidence) as a function of $N$ is required. Rather than doing separate runs with different trial values of $N$, we propose an alternative approach using a trans-dimensional Markov Chain Monte Carlo method within Nested Sampling. The posterior distribution for $N$ can be obtained with a single run. We apply the method to $ν$ Oph and Gliese 581, finding moderate evidence for additional signals in $ν$ Oph with periods of 36.11 $\pm$ 0.034 days, 75.58 $\pm$ 0.80 days, and 1709 $\pm$ 183 days; the posterior probability that at least one of these exists is 85%. The results also suggest Gliese 581 hosts many (7-15) "planets" (or other causes of other periodic signals), but only 4-6 have well determined periods. The analysis of both of these datasets shows phase transitions exist which are difficult to negotiate without Nested Sampling. △ Less

Submitted 27 January, 2015; originally announced January 2015.

Comments: Accepted for publication in MNRAS. 9 pages, 12 figures. Code at http://www.github.com/eggplantbren/Exoplanet

Space Telescope and Optical Reverberation Mapping Project. I. Ultraviolet Observations of the Seyfert 1 Galaxy NGC 5548 with the Cosmic Origins Spectrograph on Hubble Space Telescope

Authors: G. De Rosa, B. M. Peterson, J. Ely, G. A. Kriss, D. M. Crenshaw, Keith Horne, K. T. Korista, H. Netzer, R. W. Pogge, P. Arevalo, A. J. Barth, M. C. Bentz, W. N. Brandt, A. A. Breeveld, B. J. Brewer, E. Dalla Bonta, A. De Lorenzo-Caceres, K. D. Denney, M. Dietrich, R. Edelson, P. A. Evans, M. M. Fausnaugh, N. Gehrels, J. M. Gelbord, M. R. Goad , et al. (25 additional authors not shown)

Abstract: We describe the first results from a six-month long reverberation-mapping experiment in the ultraviolet based on 170 observations of the Seyfert 1 galaxy NGC 5548 with the Cosmic Origins Spectrograph on the Hubble Space Telescope. Significant correlated variability is found in the continuum and broad emission lines, with amplitudes ranging from ~30% to a factor of two in the emission lines and a f… ▽ More We describe the first results from a six-month long reverberation-mapping experiment in the ultraviolet based on 170 observations of the Seyfert 1 galaxy NGC 5548 with the Cosmic Origins Spectrograph on the Hubble Space Telescope. Significant correlated variability is found in the continuum and broad emission lines, with amplitudes ranging from ~30% to a factor of two in the emission lines and a factor of three in the continuum. The variations of all the strong emission lines lag behind those of the continuum, with He II 1640 lagging behind the continuum by ~2.5 days and Lyman alpha 1215, C IV 1550, and Si IV 1400 lagging by ~5-6 days. The relationship between the continuum and emission lines is complex. In particular, during the second half of the campaign, all emission-line lags increased by a factor of 1.3-2 and differences appear in the detailed structure of the continuum and emission-line light curves. Velocity-resolved cross-correlation analysis shows coherent structure in lag versus line-of-sight velocity for the emission lines; the high-velocity wings of C IV respond to continuum variations more rapidly than the line core, probably indicating higher velocity BLR clouds at smaller distances from the central engine. The velocity-dependent response of Lyman alpha, however, is more complex and will require further analysis. △ Less

Submitted 27 April, 2015; v1 submitted 23 January, 2015; originally announced January 2015.

Comments: Accepted for publication in The Astrophysical Journal. 20 pages, 6 figures, 4 tables. See also STORM Paper II "Space telescope and optical reverberation mapping project. II. Reverberation mapping of the accretion disk with SWIFT and HST" by R. Edelson et al

Journal ref: ApJ 806 (2015) 128

Space Telescope and Optical Reverberation Mapping Project. II. Swift and HST Reverberation Mapping of the Accretion Disk of NGC 5548

Authors: R. Edelson, J. M. Gelbord, K. Horne, I. M. McHardy, B. M. Peterson, P. Arevalo, A. A. Breeveld, G. De Rosa, P. A. Evans, M. R. Goad, G. A. Kriss, W. N. Brandt, N. Gehrels, D. Grupe, J. A. Kennea, C. S. Kochanek, J. A. Nousek, I. Papadakis, M. Siegel, D. Starkey, P. Uttley, S. Vaughan, S. Young, A. J. Barth, M. C. Bentz , et al. (25 additional authors not shown)

Abstract: Recent intensive Swift monitoring of the Seyfert 1 galaxy NGC 5548 yielded 282 usable epochs over 125 days across six UV/optical bands and the X-rays. This is the densest extended AGN UV/optical continuum sampling ever obtained, with a mean sampling rate <0.5 day. Approximately daily HST UV sampling was also obtained. The UV/optical light curves show strong correlations (r_max = 0.57 - 0.90) and t… ▽ More Recent intensive Swift monitoring of the Seyfert 1 galaxy NGC 5548 yielded 282 usable epochs over 125 days across six UV/optical bands and the X-rays. This is the densest extended AGN UV/optical continuum sampling ever obtained, with a mean sampling rate <0.5 day. Approximately daily HST UV sampling was also obtained. The UV/optical light curves show strong correlations (r_max = 0.57 - 0.90) and the clearest measurement to date of interband lags. These lags are well-fit by a τpropto λ^4/3 wavelength dependence, with a normalization that indicates an unexpectedly large disk radius of 0.35 +/- 0.05 lt-day at 1367 A, assuming a simple face-on model. The U-band shows a marginally larger lag than expected from the fit and surrounding bands, which could be due to Balmer continuum emission from the broad-line region as suggested by Korista and Goad. The UV/X-ray correlation is weaker (r_max < 0.45) and less consistent over time. This indicates that while Swift is beginning to measure UV/optical lags in general agreement with accretion disk theory (although the derived size is larger than predicted), the relationship with X-ray variability is less well understood. Combining this accretion disk size estimate with those from quasar microlensing studies suggests that AGN disk sizes scale approximately linearly with central black hole mass over a wide range of masses. △ Less

Submitted 7 April, 2015; v1 submitted 23 January, 2015; originally announced January 2015.

Comments: Accepted for publication in the Astrophysical Journal. Seventeen pages, 10 figures, 6 tables. See also STORM Paper I: "Space Telescope and Optical Reverberation Mapping Project. I. Ultraviolet Observations of the Seyfert 1 Galaxy NGC 5548 with the Cosmic Origins Spectrograph on Hubble Space Telescope" by G. De Rosa et al., http://arxiv.org/abs/1501.05954

Journal ref: ApJ 806 (2015) 129

Dissecting magnetar variability with Bayesian hierarchical models

Authors: D. Huppenkothen, B. J. Brewer, D. W. Hogg, I. Murray, M. Frean, C. Elenbaas, A. L. Watts, Y. Levin, A. J. van der Horst, C. Kouveliotou

Abstract: Neutron stars are a prime laboratory for testing physical processes under conditions of strong gravity, high density, and extreme magnetic fields. Among the zoo of neutron star phenomena, magnetars stand out for their bursting behaviour, ranging from extremely bright, rare giant flares to numerous, less energetic recurrent bursts. The exact trigger and emission mechanisms for these bursts are not… ▽ More Neutron stars are a prime laboratory for testing physical processes under conditions of strong gravity, high density, and extreme magnetic fields. Among the zoo of neutron star phenomena, magnetars stand out for their bursting behaviour, ranging from extremely bright, rare giant flares to numerous, less energetic recurrent bursts. The exact trigger and emission mechanisms for these bursts are not known; favoured models involve either a crust fracture and subsequent energy release into the magnetosphere, or explosive reconnection of magnetic field lines. In the absence of a predictive model, understanding the physical processes responsible for magnetar burst variability is difficult. Here, we develop an empirical model that decomposes magnetar bursts into a superposition of small spike-like features with a simple functional form, where the number of model components is itself part of the inference problem. The cascades of spikes that we model might be formed by avalanches of reconnection, or crust rupture aftershocks. Using Markov Chain Monte Carlo (MCMC) sampling augmented with reversible jumps between models with different numbers of parameters, we characterise the posterior distributions of the model parameters and the number of components per burst. We relate these model parameters to physical quantities in the system, and show for the first time that the variability within a burst does not conform to predictions from ideas of self-organised criticality. We also examine how well the properties of the spikes fit the predictions of simplified cascade models for the different trigger mechanisms. △ Less

Submitted 29 June, 2015; v1 submitted 21 January, 2015; originally announced January 2015.

Comments: accepted for publication in The Astrophysical Journal; code available at https://bitbucket.org/dhuppenkothen/magnetron, data products at http://figshare.com/articles/SGR_J1550_5418_magnetron_data/1292424

Inference for Trans-dimensional Bayesian Models with Diffusive Nested Sampling

Authors: Brendon J. Brewer

Abstract: Many inference problems involve inferring the number $N$ of components in some region, along with their properties $\{\mathbf{x}_i\}_{i=1}^N$, from a dataset $\mathcal{D}$. A common statistical example is finite mixture modelling. In the Bayesian framework, these problems are typically solved using one of the following two methods: i) by executing a Monte Carlo algorithm (such as Nested Sampling)… ▽ More Many inference problems involve inferring the number $N$ of components in some region, along with their properties $\{\mathbf{x}_i\}_{i=1}^N$, from a dataset $\mathcal{D}$. A common statistical example is finite mixture modelling. In the Bayesian framework, these problems are typically solved using one of the following two methods: i) by executing a Monte Carlo algorithm (such as Nested Sampling) once for each possible value of $N$, and calculating the marginal likelihood or evidence as a function of $N$; or ii) by doing a single run that allows the model dimension $N$ to change (such as Markov Chain Monte Carlo with birth/death moves), and obtaining the posterior for $N$ directly. In this paper we present a general approach to this problem that uses trans-dimensional MCMC embedded within a Nested Sampling algorithm, allowing us to explore the posterior distribution and calculate the marginal likelihood (summed over $N$) even if the problem contains a phase transition or other difficult features such as multimodality. We present two example problems, finding sinusoidal signals in noisy data, and finding and measuring galaxies in a noisy astronomical image. Both of the examples demonstrate phase transitions in the relationship between the likelihood and the cumulative prior mass, highlighting the need for Nested Sampling. △ Less

Submitted 14 January, 2015; v1 submitted 14 November, 2014; originally announced November 2014.

Comments: Only published here for the time being. 17 pages, 10 figures. Software available at https://github.com/eggplantbren/RJObject

Modeling reverberation mapping data I: improved geometric and dynamical models and comparison with cross-correlation results

Authors: Anna Pancoast, Brendon J. Brewer, Tommaso Treu

Abstract: We present an improved and expanded simply parameterized phenomenological model of the broad line region (BLR) in active galactic nuclei (AGN) for modeling reverberation mapping data. By modeling reverberation mapping data directly, we can constrain the geometry and dynamics of the BLR and measure the black hole mass without relying on the normalization factor needed in the traditional analysis. F… ▽ More We present an improved and expanded simply parameterized phenomenological model of the broad line region (BLR) in active galactic nuclei (AGN) for modeling reverberation mapping data. By modeling reverberation mapping data directly, we can constrain the geometry and dynamics of the BLR and measure the black hole mass without relying on the normalization factor needed in the traditional analysis. For realistic simulated reverberation mapping datasets of high-quality, we can recover the black hole mass to $0.05-0.25$ dex uncertainty and distinguish between dynamics dominated by elliptical orbits and inflowing gas. While direct modeling of the integrated emission line light curve allows for measurement of the mean time lag, other details of the geometry of the BLR are better constrained by the full spectroscopic dataset of emission line profiles. We use this improved model of the BLR to explore possible sources of uncertainty in measurements of the time lag using cross-correlation function (CCF) analysis and in measurements of the black hole mass using the virial product. Sampling the range of geometries and dynamics in our model of the BLR suggests that the theoretical uncertainty in black hole masses measured using the virial product is on the order of 0.25 dex. These results support the use of the CCF to measure time lags and the virial product to measure black hole masses when direct modeling techniques cannot be applied, provided the uncertainties associated with the interpretation of the results are taken into account. △ Less

Submitted 24 November, 2014; v1 submitted 10 July, 2014; originally announced July 2014.

Comments: 20 pages, 14 figures, 4 tables, published in MNRAS

Journal ref: Pancoast, A., Brewer, B. J., & Treu, T. 2014, MNRAS, 445, 3055

Hierarchical Reverberation Mapping

Authors: Brendon J. Brewer, Tom M. Elliott

Abstract: Reverberation mapping (RM) is an important technique in studies of active galactic nuclei (AGN). The key idea of RM is to measure the time lag $τ$ between variations in the continuum emission from the accretion disc and subsequent response of the broad line region (BLR). The measurement of $τ$ is typically used to estimate the physical size of the BLR and is combined with other measurements to est… ▽ More Reverberation mapping (RM) is an important technique in studies of active galactic nuclei (AGN). The key idea of RM is to measure the time lag $τ$ between variations in the continuum emission from the accretion disc and subsequent response of the broad line region (BLR). The measurement of $τ$ is typically used to estimate the physical size of the BLR and is combined with other measurements to estimate the black hole mass $M_{\rm BH}$. A major difficulty with RM campaigns is the large amount of data needed to measure $τ$. Recently, Fine et al (2012) introduced a new approach to RM where the BLR light curve is sparsely sampled, but this is counteracted by observing a large sample of AGN, rather than a single system. The results are combined to infer properties of the sample of AGN. In this letter we implement this method using a hierarchical Bayesian model and contrast this with the results from the previous stacked cross-correlation technique. We find that our inferences are more precise and allow for more straightforward interpretation than the stacked cross-correlation results. △ Less

Submitted 3 December, 2013; originally announced December 2013.

Comments: Accepted for publication in MNRAS Letters. 5 pages, 5 figures. Source code for this paper is available at https://github.com/eggplantbren/RMHB

Modeling reverberation mapping data II: dynamical modeling of the Lick AGN Monitoring Project 2008 dataset

Authors: Anna Pancoast, Brendon J. Brewer, Tommaso Treu, Daeseong Park, Aaron J. Barth, Misty C. Bentz, Jong-Hak Woo

Abstract: We present dynamical modeling of the broad line region (BLR) for a sample of five Seyfert 1 galaxies using reverberation mapping data taken by the Lick AGN Monitoring Project in 2008. By modeling the AGN continuum light curve and H$β$ line profiles directly we are able to constrain the geometry and kinematics of the BLR and make a measurement of the black hole mass that does not depend upon the vi… ▽ More We present dynamical modeling of the broad line region (BLR) for a sample of five Seyfert 1 galaxies using reverberation mapping data taken by the Lick AGN Monitoring Project in 2008. By modeling the AGN continuum light curve and H$β$ line profiles directly we are able to constrain the geometry and kinematics of the BLR and make a measurement of the black hole mass that does not depend upon the virial factor, $f$, needed in traditional reverberation mapping analysis. We find that the geometry of the BLR is generally a thick disk viewed close to face-on. While the H$β$ emission is found to come preferentially from the far side of the BLR, the mean size of the BLR is consistent with the lags measured with cross-correlation analysis. The BLR kinematics are found to be consistent with either inflowing motions or elliptical orbits, often with some combination of the two. We measure black hole masses of $\log_{10}(M_{\rm\,BH}/M_\odot)=6.62^{+0.10}_{-0.13}$ for Arp 151, $7.42^{+0.26}_{-0.27}$ for Mrk 1310, $7.51^{+0.23}_{-0.14}$ for NGC 5548, $6.42^{+0.24}_{-0.18}$ for NGC 6814, and $6.99^{+0.32}_{-0.25}$ for SBS 1116+583A. The $f$ factors measured individually for each AGN are found to correlate with inclination angle, although not with $M_{\rm\,BH}$, $L_{5100}$, or FWHM/$σ$ of the emission line profile. △ Less

Submitted 16 January, 2015; v1 submitted 25 November, 2013; originally announced November 2013.

Comments: 21 pages, 24 figures, 3 tables, Accepted for publication in MNRAS, corrected masses for NGC 5548 and NGC 6814 in the abstract

The SWELLS Survey. VI. hierarchical inference of the initial mass functions of bulges and discs

Authors: Brendon J. Brewer, Philip J. Marshall, Matthew W. Auger, Tommaso Treu, Aaron A. Dutton, Matteo Barnabè

Abstract: The long-standing assumption that the stellar initial mass function (IMF) is universal has recently been challenged by a number of observations. Several studies have shown that a "heavy" IMF (e.g., with a Salpeter-like abundance of low mass stars and thus normalisation) is preferred for massive early-type galaxies, while this IMF is inconsistent with the properties of less massive, later-type gala… ▽ More The long-standing assumption that the stellar initial mass function (IMF) is universal has recently been challenged by a number of observations. Several studies have shown that a "heavy" IMF (e.g., with a Salpeter-like abundance of low mass stars and thus normalisation) is preferred for massive early-type galaxies, while this IMF is inconsistent with the properties of less massive, later-type galaxies. These discoveries motivate the hypothesis that the IMF may vary (possibly very slightly) across galaxies and across components of individual galaxies (e.g. bulges vs discs). In this paper we use a sample of 19 late-type strong gravitational lenses from the SWELLS survey to investigate the IMFs of the bulges and discs in late-type galaxies. We perform a joint analysis of the galaxies' total masses (constrained by strong gravitational lensing) and stellar masses (constrained by optical and near-infrared colours in the context of a stellar population synthesis [SPS] model, up to an IMF normalisation parameter). Using minimal assumptions apart from the physical constraint that the total stellar mass within any aperture must be less than the total mass within the aperture, we find that the bulges of the galaxies cannot have IMFs heavier (i.e. implying high mass per unit luminosity) than Salpeter, while the disc IMFs are not well constrained by this data set. We also discuss the necessity for hierarchical modelling when combining incomplete information about multiple astronomical objects. This modelling approach allows us to place upper limits on the size of any departures from universality. More data, including spatially resolved kinematics (as in paper V) and stellar population diagnostics over a range of bulge and disc masses, are needed to robustly quantify how the IMF varies within galaxies. △ Less

Submitted 18 October, 2013; originally announced October 2013.

Comments: Accepted for publication in MNRAS. 15 pages, 8 figures. Code available at https://github.com/eggplantbren/SWELLS_Hierarchical

A transdimensional Bayesian method to infer the star formation history of resolved stellar populations

Authors: J. J. Walmswell, J. J. Eldridge, B. J. Brewer, C. A. Tout

Abstract: We propose a new method to infer the star formation histories of resolved stellar populations. With photometry one may plot observed stars on a colour-magnitude diagram (CMD) and then compare with synthetic CMDs representing different star formation histories. This has been accomplished hitherto by parametrising the model star formation history as a histogram, usually with the bin widths set by fi… ▽ More We propose a new method to infer the star formation histories of resolved stellar populations. With photometry one may plot observed stars on a colour-magnitude diagram (CMD) and then compare with synthetic CMDs representing different star formation histories. This has been accomplished hitherto by parametrising the model star formation history as a histogram, usually with the bin widths set by fixed increases in the logarithm of time. A best fit is then found with maximum likelihood methods and we consider the different means by which a likelihood can be calculated. We then apply Bayesian methods by parametrising the star formation history as an unknown number of Gaussian bursts with unknown parameters. This parametrisation automatically provides a smooth function of time. A Reversal Jump Markov Chain Monte Carlo method is then used to find both the most appropriate number of Gaussians, thus avoiding avoid overfitting, and the posterior probability distribution of the star formation rate. We apply our method to artificial populations and to observed data. We discuss the other advantages of the method: direct comparison of different parametrisations and the ability to calculate the probability that a given star is from a given Gaussian. This allows the investigation of possible sub-populations. △ Less

Submitted 5 August, 2013; originally announced August 2013.

Comments: Accepted for publication in MNRAS, 1 August 2013. 19 pages, 32 figures

The Lick AGN Monitoring Project 2011: Fe II Reverberation from the Outer Broad-Line Region

Authors: Aaron J. Barth, Anna Pancoast, Vardha N. Bennert, Brendon J. Brewer, Gabriela Canalizo, Alexei V. Filippenko, Elinor L. Gates, Jenny E. Greene, Weidong Li, Matthew A. Malkan, David J. Sand, Daniel Stern, Tommaso Treu, Jong-Hak Woo, Roberto J. Assef, Hyun-Jin Bae, Tabitha Buehler, S. Bradley Cenko, Kelsey I. Clubb, Michael C. Cooper, Aleksandar M. Diamond-Stanic, Sebastian F. Hoenig, Michael D. Joner, C. David Laney, Mariana S. Lazarova , et al. (4 additional authors not shown)

Abstract: The prominent broad Fe II emission blends in the spectra of active galactic nuclei have been shown to vary in response to continuum variations, but past attempts to measure the reverberation lag time of the optical Fe II lines have met with only limited success. Here we report the detection of Fe II reverberation in two Seyfert 1 galaxies, NGC 4593 and Mrk 1511, based on data from a program carrie… ▽ More The prominent broad Fe II emission blends in the spectra of active galactic nuclei have been shown to vary in response to continuum variations, but past attempts to measure the reverberation lag time of the optical Fe II lines have met with only limited success. Here we report the detection of Fe II reverberation in two Seyfert 1 galaxies, NGC 4593 and Mrk 1511, based on data from a program carried out at Lick Observatory in Spring 2011. Light curves for emission lines including H-beta and Fe II were measured by applying a fitting routine to decompose the spectra into several continuum and emission-line components, and we use cross-correlation techniques to determine the reverberation lags of the emission lines relative to V-band light curves. In both cases the measured lag (t_cen) of Fe II is longer than that of H-beta, although the inferred lags are somewhat sensitive to the choice of Fe II template used in the fit. For spectral decompositions done using the Fe II template of Veron-Cetty et al. (2004), we find t_cen(Fe II)/t_cen(H-beta) = 1.9+-0.6 in NGC 4593 and 1.5+-0.3 in Mrk 1511. The detection of highly correlated variations between Fe II and continuum emission demonstrates that the Fe II emission in these galaxies originates in photoionized gas, located predominantly in the outer portion of the broad-line region. △ Less

Submitted 16 April, 2013; originally announced April 2013.

Comments: 14 pages, 6 figures. Accepted for publication in ApJ

Probabilistic Catalogs for Crowded Stellar Fields

Authors: Brendon J. Brewer, Daniel Foreman-Mackey, David W. Hogg

Abstract: We present and implement a probabilistic (Bayesian) method for producing catalogs from images of stellar fields. The method is capable of inferring the number of sources N in the image and can also handle the challenges introduced by noise, overlapping sources, and an unknown point spread function (PSF). The luminosity function of the stars can also be inferred even when the precise luminosity of… ▽ More We present and implement a probabilistic (Bayesian) method for producing catalogs from images of stellar fields. The method is capable of inferring the number of sources N in the image and can also handle the challenges introduced by noise, overlapping sources, and an unknown point spread function (PSF). The luminosity function of the stars can also be inferred even when the precise luminosity of each star is uncertain, via the use of a hierarchical Bayesian model. The computational feasibility of the method is demonstrated on two simulated images with different numbers of stars. We find that our method successfully recovers the input parameter values along with principled uncertainties even when the field is crowded. We also compare our results with those obtained from the SExtractor software. While the two approaches largely agree about the fluxes of the bright stars, the Bayesian approach provides more accurate inferences about the faint stars and the number of stars, particularly in the crowded case. △ Less

Submitted 20 April, 2013; v1 submitted 25 November, 2012; originally announced November 2012.

Comments: Accepted for publication in The Astronomical Journal

The SWELLS survey - V. A Salpeter stellar initial mass function in the bulges of massive spiral galaxies

Authors: A. A. Dutton, T. Treu, B. J. Brewer, P. J. Marshall, M. W. Auger, M. Barnabe, D. C. Koo, A. S. Bolton, L. V. E. Koopmans

Abstract: Recent work has suggested that the stellar initial mass function (IMF) is not universal, but rather is correlated with galaxy stellar mass, stellar velocity dispersion, or morphological type. In this paper, we investigate variations of the IMF within individual galaxies. For this purpose, we use strong lensing and gas kinematics to measure independently the normalisation of the IMF of the bulge an… ▽ More Recent work has suggested that the stellar initial mass function (IMF) is not universal, but rather is correlated with galaxy stellar mass, stellar velocity dispersion, or morphological type. In this paper, we investigate variations of the IMF within individual galaxies. For this purpose, we use strong lensing and gas kinematics to measure independently the normalisation of the IMF of the bulge and disk components of a sample of 5 massive spiral galaxies with substantial bulge components taken from the SWELLS survey. We find that the stellar mass of the bulges are tightly constrained by the lensing and kinematic data. A comparison with masses based on stellar population synthesis models fitted to optical and near infrared photometry favors a Salpeter-like normalisation of the IMF. Conversely, the disk masses are less well constrained due to degeneracies with the dark matter halo, but are consistent with Milky Way type IMFs in agreement with previous studies. The disks are submaximal at 2.2 disk scale lengths, but due to the contribution of the bulges, the galaxies are baryon dominated at 2.2 disk scale lengths. Globally, our inferred IMF normalisation is consistent with that found for early-type galaxies of comparable stellar mass (> 10^11 M_sun). Our results suggest a non-universal IMF within the different components of spiral galaxies, adding to the well-known differences in stellar populations between disks and bulges. △ Less

Submitted 22 October, 2012; v1 submitted 19 June, 2012; originally announced June 2012.

Comments: 14 pages, 8 figures, accepted to MNRAS

The Lick AGN Monitoring Project 2011: Dynamical Modeling of the Broad Line Region in Mrk 50

Authors: Anna Pancoast, Brendon J. Brewer, Tommaso Treu, Aaron J. Barth, Vardha N. Bennert, Gabriela Canalizo, Alexei V. Filippenko, Elinor L. Gates, Jenny E. Greene, Weidong Li, Matthew A. Malkan, David J. Sand, Daniel Stern, Jong-Hak Woo, Roberto J. Assef, Hyun-Jin Bae, Tabitha Buehler, S. Bradley Cenko, Kelsey I. Clubb, Michael C. Cooper, Aleksandar M. Diamond-Stanic, Kyle D. Hiner, Sebastian F. Hoenig, Michael D. Joner, Michael T. Kandrashoff , et al. (11 additional authors not shown)

Abstract: We present dynamical modeling of the broad line region (BLR) in the Seyfert 1 galaxy Mrk 50 using reverberation mapping data taken as part of the Lick AGN Monitoring Project (LAMP) 2011. We model the reverberation mapping data directly, constraining the geometry and kinematics of the BLR, as well as deriving a black hole mass estimate that does not depend on a normalizing factor or virial coeffici… ▽ More We present dynamical modeling of the broad line region (BLR) in the Seyfert 1 galaxy Mrk 50 using reverberation mapping data taken as part of the Lick AGN Monitoring Project (LAMP) 2011. We model the reverberation mapping data directly, constraining the geometry and kinematics of the BLR, as well as deriving a black hole mass estimate that does not depend on a normalizing factor or virial coefficient. We find that the geometry of the BLR in Mrk 50 is a nearly face-on thick disk, with a mean radius of 9.6(+1.2,-0.9) light days, a width of the BLR of 6.9(+1.2,-1.1) light days, and a disk opening angle of 25\pm10 degrees above the plane. We also constrain the inclination angle to be 9(+7,-5) degrees, close to face-on. Finally, the black hole mass of Mrk 50 is inferred to be log10(M(BH)/Msun) = 7.57(+0.44,-0.27). By comparison to the virial black hole mass estimate from traditional reverberation mapping analysis, we find the normalizing constant (virial coefficient) to be log10(f) = 0.78(+0.44,-0.27), consistent with the commonly adopted mean value of 0.74 based on aligning the M(BH)-σ* relation for AGN and quiescent galaxies. While our dynamical model includes the possibility of a net inflow or outflow in the BLR, we cannot distinguish between these two scenarios. △ Less

Submitted 16 May, 2012; originally announced May 2012.

Comments: Accepted for publication in ApJ. 8 pages, 6 figures

The SWELLS survey. IV. Precision measurements of the stellar and dark matter distributions in a spiral lens galaxy

Authors: Matteo Barnabè, Aaron A. Dutton, Philip J. Marshall, Matthew W. Auger, Brendon J. Brewer, Tommaso Treu, Adam S. Bolton, David C. Koo, Léon V. E. Koopmans

Abstract: We construct a fully self-consistent mass model for the lens galaxy J2141 at z=0.14, and use it to improve on previous studies by modelling its gravitational lensing effect, gas rotation curve and stellar kinematics simultaneously. We adopt a very flexible axisymmetric mass model constituted by a generalized NFW dark matter halo and a stellar mass distribution obtained by deprojecting the MGE fit… ▽ More We construct a fully self-consistent mass model for the lens galaxy J2141 at z=0.14, and use it to improve on previous studies by modelling its gravitational lensing effect, gas rotation curve and stellar kinematics simultaneously. We adopt a very flexible axisymmetric mass model constituted by a generalized NFW dark matter halo and a stellar mass distribution obtained by deprojecting the MGE fit to the high-resolution K'-band LGSAO imaging data of the galaxy, with the (spatially constant) M/L ratio as a free parameter. We model the stellar kinematics by solving the anisotropic Jeans equations. We find that the inner logarithmic slope of the dark halo is weakly constrained (gamma = 0.82^{+0.65}_{-0.54}), and consistent with an unmodified NFW profile. We infer the galaxy to have (i) a dark matter fraction within 2.2 disk radii of 0.28^{+0.15}_{-0.10}, independent of the galaxy stellar population, implying a maximal disk for J2141; (ii) an apparently uncontracted dark matter halo, with concentration c_{-2} = 7.7_{-2.5}^{+4.2} and virial velocity v_{vir} = 242_{-39}^{+44} km/s, consistent with LCDM predictions; (iii) a slightly oblate halo (q_h = 0.75^{+0.27}_{-0.16}), consistent with predictions from baryon-affected models. Comparing the stellar mass inferred from the combined analysis (log_{10} Mstar/Msun = 11.12_{-0.09}^{+0.05}) with that inferred from SPS modelling of the galaxies colours, and accounting for a cold gas fraction of 20+/-10%, we determine a preference for a Chabrier IMF over Salpeter IMF by a Bayes factor of 5.7 (substantial evidence). We infer a value beta_{z} = 1 - sigma^2_{z}/sigma^2_{R} = 0.43_{-0.11}^{+0.08} for the orbital anisotropy parameter in the meridional plane, in agreement with most studies of local disk galaxies, and ruling out at 99% CL that the dynamics of this system can be described by a two-integral distribution function. [Abridged] △ Less

Submitted 23 March, 2012; v1 submitted 9 January, 2012; originally announced January 2012.

Comments: Accepted for publication in MNRAS. 17 pages, 9 figures

The SWELLS survey. III. Disfavouring "heavy" initial mass functions for spiral lens galaxies

Authors: Brendon J. Brewer, Aaron A. Dutton, Tommaso Treu, Matthew W. Auger, Philip J. Marshall, Matteo Barnabè, Adam S. Bolton, David C. Koo, Léon V. E. Koopmans

Abstract: We present gravitational lens models for 20 strong gravitational lens systems observed as part of the Sloan WFC Edge-on Late-type Lens Survey (SWELLS) project. Fifteen of the lenses are taken from paper I while five are newly discovered systems. The systems are galaxy-galaxy lenses where the foreground deflector has an inclined disc, with a wide range of morphological types, from late-type spiral… ▽ More We present gravitational lens models for 20 strong gravitational lens systems observed as part of the Sloan WFC Edge-on Late-type Lens Survey (SWELLS) project. Fifteen of the lenses are taken from paper I while five are newly discovered systems. The systems are galaxy-galaxy lenses where the foreground deflector has an inclined disc, with a wide range of morphological types, from late-type spiral to lenticular. For each system, we compare the total mass inside the critical curve inferred from gravitational lens modelling to the stellar mass inferred from stellar population synthesis (SPS) models, computing the stellar mass fraction f* = M(SPS)/M(lens). We find that, for the lower mass SWELLS systems, adoption of a Salpeter stellar initial mass function (IMF) leads to estimates of f* that exceed 1. This is unphysical, and provides strong evidence against the Salpeter IMF being valid for these systems. Taking the lower mass end of the SWELLS sample sigma(SIE) < 230 km/s, we find that the IMF is lighter (in terms of stellar mass-to-light ratio) than Salpeter with 98% probability, and consistent with the Chabrier IMF and IMFs between the two. This result is consistent with previous studies of spiral galaxies based on independent techniques. In combination with recent studies of massive early-type galaxies that have favoured a heavier Salpeter-like IMF, this result strengthens the evidence against a universal stellar IMF. △ Less

Submitted 5 March, 2012; v1 submitted 8 January, 2012; originally announced January 2012.

Comments: Accepted for publication in MNRAS. Some changes (none major) to address the referee's comments. 18 pages, 8 figures

The Lick AGN Monitoring Project 2011: Reverberation Mapping of Markarian 50

Authors: A. J. Barth, A. Pancoast, S. J. Thorman, V. N. Bennert, D. J. Sand, W. Li, G. Canalizo, A. V. Filippenko, E. L. Gates, J. E. Greene, M. A. Malkan, D. Stern, T. Treu, J. -H. Woo, R. J. Assef, H. -J. Bae, B. J. Brewer, T. Buehler, S. B. Cenko, K. I. Clubb, M. C. Cooper, A. M. Diamond-Stanic, K. D. Hiner, S. F. Hoenig, M. D. Joner , et al. (24 additional authors not shown)

Abstract: The Lick AGN Monitoring Project 2011 observing campaign was carried out over the course of 11 weeks in Spring 2011. Here we present the first results from this program, a measurement of the broad-line reverberation lag in the Seyfert 1 galaxy Mrk 50. Combining our data with supplemental observations obtained prior to the start of the main observing campaign, our dataset covers a total duration of… ▽ More The Lick AGN Monitoring Project 2011 observing campaign was carried out over the course of 11 weeks in Spring 2011. Here we present the first results from this program, a measurement of the broad-line reverberation lag in the Seyfert 1 galaxy Mrk 50. Combining our data with supplemental observations obtained prior to the start of the main observing campaign, our dataset covers a total duration of 4.5 months. During this time, Mrk 50 was highly variable, exhibiting a maximum variability amplitude of a factor of 4 in the U-band continuum and a factor of 2 in the H-beta line. Using standard cross-correlation techniques, we find that H-beta and H-gamma lag the V-band continuum by tau_cen = 10.64(-0.93,+0.82) and 8.43(-1.28,+1.30) days, respectively, while the lag of He II 4686 is unresolved. The H-beta line exhibits a symmetric velocity-resolved reverberation signature with shorter lags in the high-velocity wings than in the line core, consistent with an origin in a broad-line region dominated by orbital motion rather than infall or outflow. Assuming a virial normalization factor of f=5.25, the virial estimate of the black hole mass is (3.2+-0.5)*10^7 solar masses. These observations demonstrate that Mrk 50 is among the most promising nearby active galaxies for detailed investigations of broad-line region structure and dynamics. △ Less

Submitted 31 October, 2011; originally announced November 2011.

Comments: Accepted for publication in ApJ Letters. 6 pages, 4 figures

The SWELLS Survey. I. A large spectroscopically selected sample of edge-on late-type lens galaxies

Authors: T. Treu, A. A. Dutton, M. W. Auger, P. J. Marshall, A. S. Bolton, B. J. Brewer, D. Koo, L. V. E. Koopmans

Abstract: The relative contribution of baryons and dark matter to the inner regions of spiral galaxies provides critical clues to their formation and evolution, but it is generally difficult to determine. For spiral galaxies that are strong gravitational lenses, however, the combination of lensing and kinematic observations can be used to break the disk-halo degeneracy. In turn, such data constrain fundamen… ▽ More The relative contribution of baryons and dark matter to the inner regions of spiral galaxies provides critical clues to their formation and evolution, but it is generally difficult to determine. For spiral galaxies that are strong gravitational lenses, however, the combination of lensing and kinematic observations can be used to break the disk-halo degeneracy. In turn, such data constrain fundamental parameters such as i) the mass density profile slope and axis ratio of the dark matter halo, and by comparison with dark matter-only numerical simulations the modifications imposed by baryons; ii) the mass in stars and therefore the overall star formation efficiency, and the amount of feedback; iii) by comparison with stellar population synthesis models, the normalization of the stellar initial mass function. In this first paper of a series, we present a sample of 16 secure, 1 probable, and 6 possible strong lensing spiral galaxies, for which multi-band high-resolution images and rotation curves were obtained using the Hubble Space Telescope and Keck-II Telescope as part of the Sloan WFC Edge-on Late-type Lens Survey (SWELLS). The sample includes 8 newly discovered secure systems. [abridged] We find that the SWELLS sample of secure lenses spans a broad range of morphologies (from lenticular to late-type spiral), spectral types (quantified by Halpha emission), and bulge to total stellar mass ratio (0.22-0.85), while being limited to M_*>10^{10.5} M_sun. The SWELLS sample is thus well-suited for exploring the relationship between dark and luminous matter in a broad range of galaxies. We find that the deflector galaxies obey the same size-mass relation as that of a comparison sample of elongated non-lens galaxies selected from the SDSS survey. We conclude that the SWELLS sample is consistent with being representative of the overall population of high-mass high-inclination disky galaxies. △ Less

Submitted 11 July, 2011; v1 submitted 29 April, 2011; originally announced April 2011.

Comments: 21 pages, 6 figures, MNRAS, in press

The Mass of the Black Hole in Arp 151 from Bayesian Modeling of Reverberation Mapping Data

Authors: Brendon J. Brewer, Tommaso Treu, Anna Pancoast, Aaron J. Barth, Vardha N. Bennert, Misty C. Bentz, Alexei V. Filippenko, Jenny E. Greene, Matthew A. Malkan, Jong-Hak Woo

Abstract: Supermassive black holes are believed to be ubiquitous at the centers of galaxies. Measuring their masses is extremely challenging yet essential for understanding their role in the formation and evolution of cosmic structure. We present a direct measurement of the mass of a black hole in an active galactic nucleus (Arp 151) based on the motion of the gas responsible for the broad emission lines. B… ▽ More Supermassive black holes are believed to be ubiquitous at the centers of galaxies. Measuring their masses is extremely challenging yet essential for understanding their role in the formation and evolution of cosmic structure. We present a direct measurement of the mass of a black hole in an active galactic nucleus (Arp 151) based on the motion of the gas responsible for the broad emission lines. By analyzing and modeling spectroscopic and photometric time series, we find that the gas is well described by a disk or torus with an average radius of 3.99 +- 1.25 light days and an opening angle of 68.9 (+21.4, -17.2) degrees, viewed at an inclination angle of 67.8 +- 7.8 degrees (that is, closer to face-on than edge-on). The black hole mass is inferred to be 10^(6.51 +- 0.28) solar masses. The method is fully general and can be used to determine the masses of black holes at arbitrary distances, enabling studies of their evolution over cosmic time. △ Less

Submitted 25 April, 2011; originally announced April 2011.

Comments: Accepted for publication in ApJ Letters

Geometric and Dynamical Models of Reverberation Mapping Data

Authors: Anna Pancoast, Brendon J. Brewer, Tommaso Treu

Abstract: We present a general method to analyze reverberation mapping data that provides both estimates for the black hole mass and for the geometry and dynamics of the broad line region (BLR) in active galactic nuclei (AGN). Our method directly infers the spatial and velocity distribution of the BLR from the data, allowing us to easily derive a velocity-resolved transfer function and allowing for a self-c… ▽ More We present a general method to analyze reverberation mapping data that provides both estimates for the black hole mass and for the geometry and dynamics of the broad line region (BLR) in active galactic nuclei (AGN). Our method directly infers the spatial and velocity distribution of the BLR from the data, allowing us to easily derive a velocity-resolved transfer function and allowing for a self-consistent estimate of the black hole mass without a virial coefficient. We obtain estimates and reasonable uncertainties of the BLR model parameters by implementing a Markov Chain Monte Carlo algorithm using the formalism of Bayesian probability theory. We use Gaussian Processes to interpolate the the continuum light curve data and create mock light curves that can be fitted to the data. We test our method by creating simulated reverberation mapping data-sets with known true parameter values and by trying to recover these parameter values using our models. We are able to recover the parameters with realistic uncertainties that depend upon the variability of the AGN and the quality of the reverberation mapping campaign. With a geometry model we can recover the mean radius of the BLR to within ~0.1dex random uncertainty for simulated data with an integrated line flux uncertainty of 1.5%, while with a dynamical model we can recover the black hole mass and the mean radius to within ~0.05dex random uncertainty, for simulated data with a line profile average signal to noise ratio of 4 per spectral pixel. These uncertainties do not include modeling errors, which are likely to be present in the analysis of real data, and should therefore be considered as lower limits to the accuracy of the method. △ Less

Submitted 4 April, 2011; v1 submitted 25 January, 2011; originally announced January 2011.

Comments: 16 pages, 16 Figures, Accepted to ApJ

The SWELLS survey. II. Breaking the disk-halo degeneracy in the spiral galaxy gravitational lens SDSS J2141-0001

Authors: A. A. Dutton, B. J. Brewer, P. J. Marshall, M. W. Auger, T. Treu, D. C. Koo, A. S. Bolton, B. P. Holden, L. V. E. Koopmans

Abstract: The degeneracy among the disk, bulge and halo contributions to galaxy rotation curves prevents an understanding of the distribution of baryons and dark matter in disk galaxies. In an attempt to break this degeneracy, we present an analysis of the spiral galaxy strong gravitational lens SDSS J2141-0001, discovered as part of the SLACS survey. We present new Hubble Space Telescope multicolor imaging… ▽ More The degeneracy among the disk, bulge and halo contributions to galaxy rotation curves prevents an understanding of the distribution of baryons and dark matter in disk galaxies. In an attempt to break this degeneracy, we present an analysis of the spiral galaxy strong gravitational lens SDSS J2141-0001, discovered as part of the SLACS survey. We present new Hubble Space Telescope multicolor imaging, gas and stellar kinematics data derived from long-slit spectroscopy, and K-band LGS adaptive optics imaging, both from the Keck telescopes. We model the galaxy as a sum of concentric axisymmetric bulge, disk and halo components and infer the contribution of each component, using information from gravitational lensing and gas kinematics. This analysis yields a best-fitting total (disk plus bulge) stellar mass of log_{10}(Mstar/Msun) = 10.99(+0.11,-0.25). The photometric data combined with stellar population synthesis models yield log_{10}(Mstar/Msun) = 10.97\pm0.07, and 11.21\pm0.07 for the Chabrier and Salpeter IMFs, respectively. Accounting for the expected gas fraction of \simeq 20% reduces the lensing plus kinematics stellar mass by 0.10\pm0.05 dex, resulting in a Bayes factor of 11.9 in favor of a Chabrier IMF. The dark matter halo is roughly spherical, with minor to major axis ratio q_{halo}=0.91(+0.15,-0.13). The dark matter halo has a maximum circular velocity of V_{max}=276(+17,-18) km/s, and a central density parameter of log_{10}Δ_{V/2}=5.9(+0.9,-0.5). This is higher than predicted for uncontracted dark matter haloes in LCDM cosmologies, log_{10}Δ_{V/2}=5.2, suggesting that either the halo has contracted in response to galaxy formation, or that the halo has a higher than average concentration. At 2.2 disk scale lengths the dark matter fraction is f_{DM}=0.55(+0.20,-0.15), suggesting that SDSS J2141-0001 is sub-maximal. △ Less

Submitted 14 March, 2011; v1 submitted 8 January, 2011; originally announced January 2011.

Comments: 24 pages, 20 figures, accepted to MNRAS, minor changes

Journal ref: 2011 MNRAS.417.1621D

A Compact Early-type Galaxy at z = 0.6 Under a Magnifying Lens: Evidence For Inside-out Growth

Authors: M. W. Auger, T. Treu, B. J. Brewer, P. J. Marshall

Abstract: We use Keck laser guide star adaptive optics imaging and exploit the magnifying effects of strong gravitational lensing (the effective resolution is FWHM ~ 200 pc) to investigate the sub-kpc scale of an intermediate-redshift (z = 0.63) massive early-type galaxy being lensed by a foreground early-type galaxy; we dub this class of strong gravitational lens systems EELs, e.g., early-type/early-type l… ▽ More We use Keck laser guide star adaptive optics imaging and exploit the magnifying effects of strong gravitational lensing (the effective resolution is FWHM ~ 200 pc) to investigate the sub-kpc scale of an intermediate-redshift (z = 0.63) massive early-type galaxy being lensed by a foreground early-type galaxy; we dub this class of strong gravitational lens systems EELs, e.g., early-type/early-type lenses. We find that the background source is massive (M* = 10^{10.9} M_sun) and compact (r_e = 1.1 kpc), and a two-component fit is required to model accurately the surface brightness distribution, including an extended low-surface-brightness component. This extended component may arise from the evolution of higher-redshift `red nuggets' or may already be in place at z ~ 2 but is unobservable due to cosmological surface brightness dimming. △ Less

Submitted 3 November, 2010; originally announced November 2010.

Comments: 5 pages, 4 figures; accepted to MNRAS

Modelling of the Complex CASSOWARY/SLUGS Gravitational Lenses

Authors: Brendon J. Brewer, Geraint F. Lewis, Vasily Belokurov, Michael J. Irwin, Terry J. Bridges, N. Wyn Evans

Abstract: We present the first high-resolution images of CSWA 31, a gravitational lens system observed as part of the SLUGS (Sloan Lenses Unravelled by Gemini Studies) program. These systems exhibit complex image structure with the potential to strongly constrain the mass distribution of the massive lens galaxies, as well as the complex morphology of the sources. In this paper, we describe the strategy used… ▽ More We present the first high-resolution images of CSWA 31, a gravitational lens system observed as part of the SLUGS (Sloan Lenses Unravelled by Gemini Studies) program. These systems exhibit complex image structure with the potential to strongly constrain the mass distribution of the massive lens galaxies, as well as the complex morphology of the sources. In this paper, we describe the strategy used to reconstruct the unlensed source profile and the lens galaxy mass profiles. We introduce a prior distribution over multi-wavelength sources that is realistic as a representation of our knowledge about the surface brightness profiles of galaxies and groups of galaxies. To carry out the inference computationally, we use Diffusive Nested Sampling, an efficient variant of Nested Sampling that uses Markov Chain Monte Carlo (MCMC) to sample the complex posterior distributions and compute the normalising constant. We demonstrate the efficacy of this approach with the reconstruction of the group-group gravitational lens system CSWA 31, finding the source to be composed of five merging spiral galaxies magnified by a factor of 13. △ Less

Submitted 23 November, 2010; v1 submitted 25 October, 2010; originally announced October 2010.

Comments: Accepted for publication in MNRAS

A comparison of Bayesian and Fourier methods for frequency determination in asteroseismology

Authors: Timothy R. White, Brendon J. Brewer, Timothy R. Bedding, Dennis Stello, Hans Kjeldsen

Abstract: Bayesian methods are becoming more widely used in asteroseismic analysis. In particular, they are being used to determine oscillation frequencies, which are also commonly found by Fourier analysis. It is important to establish whether the Bayesian methods provide an improvement on Fourier methods. We compare, using simulated data, the standard iterative sine-wave fitting method against a Markov Ch… ▽ More Bayesian methods are becoming more widely used in asteroseismic analysis. In particular, they are being used to determine oscillation frequencies, which are also commonly found by Fourier analysis. It is important to establish whether the Bayesian methods provide an improvement on Fourier methods. We compare, using simulated data, the standard iterative sine-wave fitting method against a Markov Chain Monte Carlo (MCMC) code that has been introduced to infer purely the frequencies of oscillation modes (Brewer et al. 2007). A uniform prior probability distribution function is used for the MCMC method. We find the methods do equally well at determining the correct oscillation frequencies, although the Bayesian method is able to highlight the possibility of a misidentification due to aliasing, which can be useful. In general, we suggest that the least computationally intensive method is preferable. △ Less

Submitted 21 June, 2010; originally announced June 2010.

Comments: 11 pages, 8 figures, accepted for publication in Communications in Asterosesimology

AAOmega Observations of 47 Tucanae: Evidence for a Past Merger?

Authors: Richard R. Lane, Brendon J. Brewer, László L. Kiss, Geraint F. Lewis, Rodrigo A. Ibata, Arnaud Siebert, Timothy R. Bedding, Péter Székely, Gyula M. Szabó

Abstract: The globular cluster 47 Tucanae is well studied but it has many characteristics that are unexplained, including a significant rise in the velocity dispersion profile at large radii, indicating the exciting possibility of two distinct kinematic populations. In this Letter we employ a Bayesian approach to the analysis of the largest available spectral dataset of 47 Tucanae to determine whether thi… ▽ More The globular cluster 47 Tucanae is well studied but it has many characteristics that are unexplained, including a significant rise in the velocity dispersion profile at large radii, indicating the exciting possibility of two distinct kinematic populations. In this Letter we employ a Bayesian approach to the analysis of the largest available spectral dataset of 47 Tucanae to determine whether this apparently two-component population is real. Assuming the two models were equally likely before taking the data into account, we find that the evidence favours the two-component population model by a factor of ~3x10^7. Several possible explanations for this result are explored, namely the evaporation of low-mass stars, a hierarchical merger, extant remnants of two initially segregated populations, and multiple star formation epochs. We find the most compelling explanation for the two-component velocity distribution is that 47 Tuc formed as two separate populations arising from the same proto-cluster cloud which merged <7.3 +/- 1.5 Gyr ago. This may also explain the extreme rotation, low mass-to-light ratio and mixed stellar populations of this cluster. △ Less

Submitted 2 February, 2010; originally announced February 2010.

Comments: 14 pages, 2 figures, 1 table. Accepted for publication in ApJ Letters

Journal ref: Astrophysical Journal Letters (2010), 711, L122-L126

Diffusive Nested Sampling

Authors: Brendon J. Brewer, Livia B. Pártay, Gábor Csányi

Abstract: We introduce a general Monte Carlo method based on Nested Sampling (NS), for sampling complex probability distributions and estimating the normalising constant. The method uses one or more particles, which explore a mixture of nested probability distributions, each successive distribution occupying ~e^-1 times the enclosed prior mass of the previous distribution. While NS technically requires inde… ▽ More We introduce a general Monte Carlo method based on Nested Sampling (NS), for sampling complex probability distributions and estimating the normalising constant. The method uses one or more particles, which explore a mixture of nested probability distributions, each successive distribution occupying ~e^-1 times the enclosed prior mass of the previous distribution. While NS technically requires independent generation of particles, Markov Chain Monte Carlo (MCMC) exploration fits naturally into this technique. We illustrate the new method on a test problem and find that it can achieve four times the accuracy of classic MCMC-based Nested Sampling, for the same computational effort; equivalent to a factor of 16 speedup. An additional benefit is that more samples and a more accurate evidence value can be obtained simply by continuing the run for longer, as in standard MCMC. △ Less

Submitted 26 July, 2010; v1 submitted 11 December, 2009; originally announced December 2009.

Comments: Accepted for publication in Statistics and Computing. C++ code available at http://lindor.physics.ucsb.edu/DNest

Journal ref: Statistics and Computing, 2011, 21, 4, 649-656