Normalized to: Brewer, B.

**59** article(s) in total. **305** co-authors, from **1** to **23** common article(s). Median position in authors list is **2,0**.

[1] oai:arXiv.org:1910.02778 [pdf] - 2038293

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

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

Submitted: **2019-10-07**

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$.

[2] oai:arXiv.org:1910.00808 [pdf] - 1972777

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

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

Submitted: **2019-10-02**

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.

[3] oai:arXiv.org:1903.03121 [pdf] - 1845968

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

Comments: 23 pages, 12 figures, accepted for MNRAS

Submitted: **2019-03-07**, last modified: **2019-03-10**

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$\alpha$ gas velocity
dispersion for our sample to be in the range $\bar{\sigma}_v \sim $[7, 30] km
s$^{-1}$. The relative difference between our approach and estimates using the
single Gaussian component fits per spaxel is $\Delta \bar{\sigma}_v /
\bar{\sigma}_v = - 0.29 \pm 0.18$ for the H$\alpha$ flux-weighted mean velocity
dispersion.

[4] oai:arXiv.org:1809.05113 [pdf] - 1771690

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

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

Submitted: **2018-09-13**

We present models of the H$\beta$-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},\sigma})_{\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.

[5] oai:arXiv.org:1806.08305 [pdf] - 1703383

kima: Exoplanet detection in radial velocities

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

Submitted: **2018-06-21**

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.

[6] oai:arXiv.org:1803.02318 [pdf] - 1663378

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

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

Submitted: **2018-03-06**

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.

[7] oai:arXiv.org:1705.02346 [pdf] - 1583164

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

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

Submitted: **2017-05-05**, last modified: **2017-09-27**

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.

[8] oai:arXiv.org:1707.03543 [pdf] - 1585827

Computing Entropies With Nested Sampling

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

Submitted: **2017-07-12**, last modified: **2017-08-15**

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.

[9] oai:arXiv.org:1606.03757 [pdf] - 1507581

DNest4: Diffusive Nested Sampling in C++ and Python

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

Submitted: **2016-06-12**, last modified: **2016-11-01**

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.

[10] oai:arXiv.org:1611.00008 [pdf] - 1532352

Red nuggets grow inside-out: evidence from gravitational lensing

Comments: 21 pages; accepted for publication in MNRAS

Submitted: **2016-10-31**

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\'ersic 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.

[11] oai:arXiv.org:1607.02190 [pdf] - 1475426

The elusive stellar halo of the Triangulum Galaxy

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

Submitted: **2016-07-07**

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 ${\mu}_V$ = 35.5 mags per square arcsec, or a total halo
luminosity of $L < 10^6L_{\odot}$.

[12] oai:arXiv.org:1603.08741 [pdf] - 1422270

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

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

Submitted: **2016-03-29**

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.

[13] oai:arXiv.org:1510.05648 [pdf] - 1396732

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

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

Submitted: **2015-10-19**, last modified: **2016-02-29**

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\,\AA\ to the $z$ band ($\sim\!9160$\,\AA). 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}\,$\lambda 1640$ and $\lambda 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 $\tau \propto \lambda^{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.

[14] oai:arXiv.org:1601.07495 [pdf] - 1374310

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

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

Submitted: **2016-01-27**, last modified: **2016-02-03**

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.

[15] oai:arXiv.org:1508.00662 [pdf] - 1290797

Trans-Dimensional Bayesian Inference for Gravitational Lens
Substructures

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

Submitted: **2015-08-04**, last modified: **2015-10-09**

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.

[16] oai:arXiv.org:1501.05251 [pdf] - 1276968

Dissecting magnetar variability with Bayesian hierarchical models

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

Submitted: **2015-01-21**, last modified: **2015-06-29**

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.

[17] oai:arXiv.org:1501.05954 [pdf] - 1243799

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

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

Submitted: **2015-01-23**, last modified: **2015-04-27**

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.

[18] oai:arXiv.org:1501.05951 [pdf] - 1243798

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

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

Submitted: **2015-01-23**, last modified: **2015-04-07**

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 \tau propto \lambda^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.

[19] oai:arXiv.org:1503.01146 [pdf] - 1231932

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

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

Submitted: **2015-03-03**

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.

[20] oai:arXiv.org:1501.06952 [pdf] - 1224212

Fast Bayesian Inference for Exoplanet Discovery in Radial Velocity Data

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

Submitted: **2015-01-27**

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 $\nu$ Oph and Gliese 581, finding moderate evidence for additional
signals in $\nu$ 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.

[21] oai:arXiv.org:1311.6475 [pdf] - 1201859

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

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

Submitted: **2013-11-25**, last modified: **2015-01-16**

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$\beta$ 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$\beta$ 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/$\sigma$ of the emission line profile.

[22] oai:arXiv.org:1411.3921 [pdf] - 921506

Inference for Trans-dimensional Bayesian Models with Diffusive Nested
Sampling

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

Submitted: **2014-11-14**, last modified: **2015-01-14**

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.

[23] oai:arXiv.org:1407.2941 [pdf] - 1215594

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

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

Submitted: **2014-07-10**, last modified: **2014-11-24**

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.

[24] oai:arXiv.org:1312.0919 [pdf] - 1202009

Hierarchical Reverberation Mapping

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

Submitted: **2013-12-03**

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 $\tau$
between variations in the continuum emission from the accretion disc and
subsequent response of the broad line region (BLR). The measurement of $\tau$
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
$\tau$. 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.

[25] oai:arXiv.org:1310.5177 [pdf] - 1180070

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

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

Submitted: **2013-10-18**

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.

[26] oai:arXiv.org:1308.0945 [pdf] - 1173203

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

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

Submitted: **2013-08-05**

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.

[27] oai:arXiv.org:1211.5805 [pdf] - 1158023

Probabilistic Catalogs for Crowded Stellar Fields

Comments: Accepted for publication in The Astronomical Journal

Submitted: **2012-11-25**, last modified: **2013-04-20**

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.

[28] oai:arXiv.org:1304.4643 [pdf] - 1166020

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

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

Submitted: **2013-04-16**

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.

[29] oai:arXiv.org:1206.4310 [pdf] - 1124238

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

Comments: 14 pages, 8 figures, accepted to MNRAS

Submitted: **2012-06-19**, last modified: **2012-10-22**

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.

[30] oai:arXiv.org:1205.3789 [pdf] - 1123462

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

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

Submitted: **2012-05-16**

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)-{\sigma}* 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.

[31] oai:arXiv.org:1201.1692 [pdf] - 1092799

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

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

Submitted: **2012-01-09**, last modified: **2012-03-23**

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]

[32] oai:arXiv.org:1201.1677 [pdf] - 1092795

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

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

Submitted: **2012-01-08**, last modified: **2012-03-05**

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.

[33] oai:arXiv.org:1111.0061 [pdf] - 1091318

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

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

Submitted: **2011-10-31**

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.

[34] oai:arXiv.org:1104.5663 [pdf] - 1076307

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

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

Submitted: **2011-04-29**, last modified: **2011-07-11**

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.

[35] oai:arXiv.org:1104.4794 [pdf] - 1076204

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

Comments: Accepted for publication in ApJ Letters

Submitted: **2011-04-25**

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.

[36] oai:arXiv.org:1101.4952 [pdf] - 1051689

Geometric and Dynamical Models of Reverberation Mapping Data

Comments: 16 pages, 16 Figures, Accepted to ApJ

Submitted: **2011-01-25**, last modified: **2011-04-04**

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.

[37] oai:arXiv.org:1101.1622 [pdf] - 955917

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

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

Submitted: **2011-01-08**, last modified: **2011-03-14**

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}\Delta_{V/2}=5.9(+0.9,-0.5). This is higher than predicted for
uncontracted dark matter haloes in LCDM cosmologies, log_{10}\Delta_{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.

[38] oai:arXiv.org:1010.5298 [pdf] - 1041532

Modelling of the Complex CASSOWARY/SLUGS Gravitational Lenses

Comments: Accepted for publication in MNRAS

Submitted: **2010-10-25**, last modified: **2010-11-23**

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.

[39] oai:arXiv.org:1011.0993 [pdf] - 1041729

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

Comments: 5 pages, 4 figures; accepted to MNRAS

Submitted: **2010-11-03**

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.

[40] oai:arXiv.org:0912.2380 [pdf] - 480760

Diffusive Nested Sampling

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

Submitted: **2009-12-11**, last modified: **2010-07-26**

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.

[41] oai:arXiv.org:1006.3992 [pdf] - 1033208

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

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

Submitted: **2010-06-21**

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.

[42] oai:arXiv.org:1002.0625 [pdf] - 388537

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

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

Submitted: **2010-02-02**

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.

[43] oai:arXiv.org:0906.5609 [pdf] - 25823

Entropic Priors and Bayesian Model Selection

Comments: Presented at MaxEnt 2009, the 29th International Workshop on Bayesian
Inference and Maximum Entropy Methods in Science and Engineering (July 5-10,
2009, Oxford, Mississippi, USA)

Submitted: **2009-06-30**, last modified: **2009-08-28**

We demonstrate that the principle of maximum relative entropy (ME), used
judiciously, can ease the specification of priors in model selection problems.
The resulting effect is that models that make sharp predictions are
disfavoured, weakening the usual Bayesian "Occam's Razor". This is illustrated
with a simple example involving what Jaynes called a "sure thing" hypothesis.
Jaynes' resolution of the situation involved introducing a large number of
alternative "sure thing" hypotheses that were possible before we observed the
data. However, in more complex situations, it may not be possible to explicitly
enumerate large numbers of alternatives. The entropic priors formalism produces
the desired result without modifying the hypothesis space or requiring explicit
enumeration of alternatives; all that is required is a good model for the prior
predictive distribution for the data. This idea is illustrated with a simple
rigged-lottery example, and we outline how this idea may help to resolve a
recent debate amongst cosmologists: is dark energy a cosmological constant, or
has it evolved with time in some way? And how shall we decide, when the data
are in?

[44] oai:arXiv.org:0903.3384 [pdf] - 22540

Low-Energy Astrophysics: Stimulating the Reduction of Energy Consumption
in the Next Decade

Comments: "State of the Profession" white paper submitted to the 2010 Astronomy
& Astrophysics Decadal Survey. Add your thoughts to the wiki at
http://low-energy-astro.physics.ucsb.edu

Submitted: **2009-03-19**

In this paper we address the consumption of energy by astronomers while
performing their professional duties. Although we find that astronomy uses a
negligible fraction of the US energy budget, the rate at which energy is
consumed by an average astronomer is similar to that of a typical high-flying
businessperson. We review some of the ways in which astronomers are already
acting to reduce their energy consumption. In the coming decades, all citizens
will have to reduce their energy consumption to conserve fossil fuel reserves
and to help avert a potentially catastrophic change in the Earth's climate. The
challenges are the same for astronomers as they are for everyone: decreasing
the distances we travel and investing in energy-efficient infrastructure. The
high profile of astronomy in the media, and the great public interest in our
field, can play a role in promoting energy-awareness to the wider population.
Our specific recommendations are therefore to 1) reduce travel when possible,
through efficient meeting organization, and by investing in high-bandwidth
video conference facilities and virtual-world software, 2) create
energy-efficient observatories, computing centers and workplaces, powered by
sustainable energy resources, and 3) actively publicize these pursuits.

[45] oai:arXiv.org:0902.3907 [pdf] - 21660

Gaussian Process Modelling of Asteroseismic Data

Comments: Accepted for publication in MNRAS

Submitted: **2009-02-23**

The measured properties of stellar oscillations can provide powerful
constraints on the internal structure and composition of stars. To begin this
process, oscillation frequencies must be extracted from the observational data,
typically time series of the star's brightness or radial velocity. In this
paper, a probabilistic model is introduced for inferring the frequencies and
amplitudes of stellar oscillation modes from data, assuming that there is some
periodic character to the oscillations, but that they may not be exactly
sinusoidal. Effectively we fit damped oscillations to the time series, and
hence the mode lifetime is also recovered. While this approach is
computationally demanding for large time series (> 1500 points), it should at
least allow improved analysis of observations of solar-like oscillations in
subgiant and red giant stars, as well as sparse observations of semiregular
stars, where the number of points in the time series is often low. The method
is demonstrated on simulated data and then applied to radial velocity
measurements of the red giant star xi Hydrae, yielding a mode lifetime between
0.41 and 2.65 days with 95% posterior probability. The large frequency
separation between modes is ambiguous, however we argue that the most plausible
value is 6.3 microHz, based on the radial velocity data and the star's position
in the HR diagram.

[46] oai:arXiv.org:0809.0939 [pdf] - 16010

Applications of Bayesian Probability Theory in Astrophysics

Comments: Paper chapters removed, but are available from arxiv. Alternatively
the full thesis is available at http://web.maths.unsw.edu.au/~brewer/

Submitted: **2008-09-04**

Bayesian Inference is a powerful approach to data analysis that is based
almost entirely on probability theory. In this approach, probabilities model
{\it uncertainty} rather than randomness or variability. This thesis is
composed of a series of papers that have been published in various astronomical
journals during the years 2005-2008. The unifying thread running through the
papers is the use of Bayesian Inference to solve underdetermined inverse
problems in astrophysics. Firstly, a methodology is developed to solve a
question in gravitational lens inversion - using the observed images of
gravitational lens systems to reconstruct the undistorted source profile and
the mass profile of the lensing galaxy. A similar technique is also applied to
the task of inferring the number and frequency of modes of oscillation of a
star from the time series observations that are used in the field of
asteroseismology. For these complex problems, many of the required calculations
cannot be done analytically, and so Markov Chain Monte Carlo algorithms have
been used. Finally, probabilistic reasoning is applied to a controversial
question in astrobiology: does the fact that life formed quite soon after the
Earth constitute evidence that the formation of life is quite probable, given
the right macroscopic conditions?

[47] oai:arXiv.org:0807.4969 [pdf] - 15009

The Implications of the Early Formation of Life on Earth

Comments: Submitted, but seems to have fallen into a black hole since

Submitted: **2008-07-30**

One of the most interesting unsolved questions in science today is the
question of life on other planets. At the present time it is safe to say that
we do not have much of an idea as to whether life is common or exceedingly rare
in the universe, and this will probably not be solved for certain unless
definitive evidence of extraterrestrial life is found in the future. Our
presence on Earth is just as consistent with the hypothesis that life is
extremely rare as it is with the hypothesis that it is common, since if there
was only one planet with intelligent life, we would find ourselves on it.
However, we have more information than this, such as the the surprisingly short
length of time it took for life to arise on Earth. Previous authors have
analysed this information, concluding that it is evidence that the probability
of abiogenesis is moderate ($>$ 13% with 95% probability) and cannot be
extremely small. In this paper I use simple probabilistic model to show that
this conclusion was based more on an unintentional assumption than on the data.
While the early formation of life on Earth provides some evidence in the
direction of life being common, it is far from conclusive, and in particular
does not rule out the possibility that abiogenesis has only occurred once in
the history of the universe.

[48] oai:arXiv.org:0807.2145 [pdf] - 14487

Unlensing HST Observations of the Einstein Ring 1RXS J1131-1231: A
Bayesian Analysis

Comments: MNRAS accepted

Submitted: **2008-07-14**

We present a source and lens reconstruction for the optical Einstein ring
gravitational lens system RXS J1131-1231. We resolve detail in the source,
which is the host galaxy of a $z=0.658$ quasar, down to a resolution of 0.045
arc seconds (this is the size of the smallest conclusively resolved structures,
rather than the pixel scale), using a Bayesian technique with a realistic model
for the prior information. The source reconstruction reveals a substantial
amount of complex structure in the host galaxy, which is $\sim$ 8 kpc in extent
and contains several bright compact substructures, with the quasar source
residing in one of these bright substructures. Additionally, we recover the
mass distribution of the lensing galaxy, assuming a simply-parameterised model,
using information from both the quasar images and the extended images. This
allows a direct comparison of the amount of information about the lens that is
provided by the quasar images in comparison to the extended images. In this
system, we find that the extended images provide significantly more information
about the lens than the quasar images alone, especially if we do not include
prior constraints on the central position of the lens.

[49] oai:arXiv.org:0806.4616 [pdf] - 208882

A Molecular Einstein Ring at z=4.12: Imaging the Dynamics of a Quasar
Host Galaxy Through a Cosmic Lens

Comments: 9 pages, 7 figures, to appear in ApJ (accepted June 27, 2008)

Submitted: **2008-06-27**

We present high-resolution (0.3") Very Large Array (VLA) imaging of the
molecular gas in the host galaxy of the high redshift quasar PSS J2322+1944
(z=4.12). These observations confirm that the molecular gas (CO) in the host
galaxy of this quasar is lensed into a full Einstein ring, and reveal the
internal dynamics of the molecular gas in this system. The ring has a diameter
of ~1.5", and thus is sampled over ~20 resolution elements by our observations.
Through a model-based lens inversion, we recover the velocity gradient of the
molecular reservoir in the quasar host galaxy of PSS J2322+1944. The Einstein
ring lens configuration enables us to zoom in on the emission and to resolve
scales down to ~1 kpc. From the model-reconstructed source, we find that the
molecular gas is distributed on a scale of 5 kpc, and has a total mass of
M(H2)=1.7 x 10^10 M_sun. A basic estimate of the dynamical mass gives M_dyn =
4.4 x 10^10 (sin i)^-2 M_sun, that is, only ~2.5 times the molecular gas mass,
and ~30 times the black hole mass (assuming that the dynamical structure is
highly inclined). The lens configuration also allows us to tie the optical
emission to the molecular gas emission, which suggests that the active galactic
nucleus (AGN) does reside within, but not close to the center of the molecular
reservoir. Together with the (at least partially) disturbed structure of the
CO, this suggests that the system is interacting. Such an interaction, possibly
caused by a major `wet' merger, may be responsible for both feeding the quasar
and fueling the massive starburst of 680 M_sun/yr in this system, in agreement
with recently suggested scenarios of quasar activity and galaxy assembly in the
early universe.

[50] oai:arXiv.org:0706.0795 [pdf] - 1927

Solar-like oscillations in the metal-poor subgiant nu Indi: II. Acoustic
spectrum and mode lifetime

Comments: 5 pages, 7 figures, A&A accepted

Submitted: **2007-06-06**

Convection in stars excites resonant acoustic waves which depend on the sound
speed inside the star, which in turn depends on properties of the stellar
interior. Therefore, asteroseismology is an unrivaled method to probe the
internal structure of a star. We made a seismic study of the metal-poor
subgiant star nu Indi with the goal of constraining its interior structure. Our
study is based on a time series of 1201 radial velocity measurements spread
over 14 nights obtained from two sites, Siding Spring Observatory in Australia
and ESO La Silla Observatory in Chile. The power spectrum of the high precision
velocity time series clearly presents several identifiable peaks between 200
and 500 uHz showing regularity with a large and small spacing of 25.14 +- 0.09
uHz and 2.96 +- 0.22 uHz at 330 uHz. Thirteen individual modes have been
identified with amplitudes in the range 53 to 173 cm/s. The mode damping time
is estimated to be about 16 days (1-sigma range between 9 and 50 days),
substantially longer than in other stars like the Sun, the alpha Cen system or
the giant xi Hya.

[51] oai:arXiv.org:astro-ph/0703747 [pdf] - 316990

Solar-like oscillations in the G2 subgiant beta Hydri from dual-site
observations

Comments: 13 pages, 14 figures, accepted by ApJ. Fixed minor typo (ref to Fig
14)

Submitted: **2007-03-29**, last modified: **2007-05-14**

We have observed oscillations in the nearby G2 subgiant star beta Hyi using
high-precision velocity observations obtained over more than a week with the
HARPS and UCLES spectrographs. The oscillation frequencies show a regular comb
structure, as expected for solar-like oscillations, but with several l=1 modes
being strongly affected by avoided crossings. The data, combined with those we
obtained five years earlier, allow us to identify 28 oscillation modes. By
scaling the large frequency separation from the Sun, we measure the mean
density of beta Hyi to an accuracy of 0.6%. The amplitudes of the oscillations
are about 2.5 times solar and the mode lifetime is 2.3 d. A detailed comparison
of the mixed l=1 modes with theoretical models should allow a precise estimate
of the age of the star.

[52] oai:arXiv.org:astro-ph/0609270 [pdf] - 84862

The evolution of host mass and black hole mass in QSOs from the 2dF QSO
Redshift Survey

Comments: 11 pages, 4 figures, accepted for publication in MNRAS

Submitted: **2006-09-10**

We investigate the relation between the mass of super-massive black holes
(Mbh) in QSOs and the mass of the dark matter halos hosting them (Mdh). We
measure the widths of broad emission lines (Mgii lambda 2798, Civ lambda 1549)
from QSO composite spectra as a function of redshift. These widths are then
used to determine virial black hole mass estimates.
We compare our virial black hole mass estimates to dark matter halo masses
for QSO hosts derived by Croom et al. (2005) based on measurements of QSO
clustering. This enables us to trace the Mbh-Mdh relation over the redshift
range z=0.5 to 2.5. We calculate the mean zero-point of the Mbh-Mdh relation to
be Mbh=10^(8.4+/-0.2)Msun for an Mdh=10^(12.5)Msun. These data are then
compared with several models connecting Mbh and Mdh as well as recent
hydrodynamical simulations of galaxy evolution. We note that the flux limited
nature of QSO samples can cause a Malmquist-type bias in the measured
zero-point of the Mbh-Mdh relation. The magnitude of this bias depends on the
scatter in the Mbh-Mdh relation, and we reevaluate the zero-point assuming
three published values for this scatter.
(abridged)

[53] oai:arXiv.org:astro-ph/0608571 [pdf] - 84448

Bayesian Inference from Observations of Solar-Like Oscillations

Comments: Accepted for publication in ApJ

Submitted: **2006-08-26**

Stellar oscillations can provide a wealth of information about a star, which
can be extracted from observed time series of the star's brightness or radial
velocity. In this paper we address the question of how to extract as much
information as possible from such a dataset. We have developed a Markov Chain
Monte Carlo (MCMC) code that is able to infer the number of oscillation
frequencies present in the signal and their values (with corresponding
uncertainties), without having to fit the amplitudes and phases. Gaps in the
data do not have any serious consequences for this method; in cases where
severe aliasing exists, any ambiguity in the frequency determinations will be
reflected in the results. It also allows us to infer parameters of the
frequency pattern, such as the large separation Delta nu. We have previously
applied this method to the star nu Indi (Bedding et al 2006), and here we
describe the method fully and apply it to simulated datasets, showing that the
code is able to give correct results even when some of the model assumptions
are violated. In particular, the non-sinusoidal nature of the individual
oscillation modes due to stochastic excitation and damping has no major impact
on the usefulness of our approach.

[54] oai:arXiv.org:astro-ph/0606714 [pdf] - 83154

The Einstein Ring 0047-2808 Revisited: A Bayesian Inversion

Comments: Accepted for publication in ApJ

Submitted: **2006-06-29**

In a previous paper, we outlined a new Bayesian method for inferring the
properties of extended gravitational lenses, given data in the form of resolved
images. This method holds the most promise for optimally extracting information
from the observed image, whilst providing reliable uncertainties in all
parameters. Here, we apply the method to the well studied optical Einstein ring
0047-2808. Our results are in broad agreement with previous studies, showing
that the density profile of the lensing galaxy is aligned within a few degrees
of the light profile, and suggesting that the source galaxy (at redshift 3.6)
is a binary system, although its size is only of order 1-2 kpc. We also find
that the mass of the elliptical lensing galaxy enclosed by the image is
(2.91$\pm$0.01)$\times10^{11}$ M$_{\sun}$. Our method is able to achieve
improved resolution for the source reconstructions, although we also find that
some of the uncertainties are greater than has been found in previous analyses,
due to the inclusion of extra pixels and a more general lens model.

[55] oai:arXiv.org:astro-ph/0604453 [pdf] - 81522

Solar-like oscillations in the metal-poor subgiant nu Indi: constraining
the mass and age using asteroseismology

Comments: accpted for publication in ApJ

Submitted: **2006-04-20**

Asteroseismology is a powerful method for determining fundamental properties
of stars. We report the first application to a metal-poor object, namely the
subgiant star nu Ind. We measured precise velocities from two sites, allowing
us to detect oscillations and infer a large frequency separation of Delta_nu =
24.25 +/- 0.25 microHz. Combining this value with the location of the star in
the H-R diagram and comparing with standard evolutionary models, we were able
to place constraints on the stellar parameters. In particular, our results
indicate that nu Ind has a low mass (0.85 +/- 0.04 M_sun) and is at least 9 Gyr
old.

[56] oai:arXiv.org:astro-ph/0509863 [pdf] - 76444

Strong Gravitational Lens Inversion: A Bayesian Approach

Comments: Accepted for publication in ApJ

Submitted: **2005-09-28**

If an extended source, such as a galaxy, is gravitationally lensed by a
massive object in the foreground, the lensing distorts the observed image. It
is straightforward to simulate what the observed image would be for a
particular lens and source combination. In practice, one observes the lensed
image on the sky, but blurred by atmospheric and telescopic effects and also
contaminated with noise. The question that then arises is, given this
incomplete data, what combinations of lens mass distribution and source surface
brightness profile could plausibly have produced this image? This is a classic
example of an inverse problem, and the method for solving it is given by the
framework of Bayesian inference. In this paper we demonstrate the application
of Bayesian inference to the problem of gravitational lens reconstruction, and
illustrate the use of Markov Chain Monte Carlo simulations which can be used
when the analytical calculations become too difficult. Previous methods for
performing gravitational lens inversion are seen in a new light, as special
cases of the general approach presented in this paper. Thus, we are able to
answer, at least in principle, lingering questions about the uncertainties in
the reconstructed source and lens parameters, taking into account all of the
data and any prior information we may have.

[57] oai:arXiv.org:astro-ph/0507471 [pdf] - 74640

The light curve of the semiregular variable L2 Puppis: II. Evidence for
solar-like excitation of the oscillations

Comments: 7 pages, 7 figures, MNRAS in press

Submitted: **2005-07-20**

We analyse visual observations of the pulsations of the red giant variable L2
Pup. The data cover 77 years between 1927 and 2005, thus providing an extensive
empirical base for characterizing properties of the oscillations. The power
spectrum of the light curve shows a single mode resolved into multiple peaks
under a narrow envelope. We argue that this results from stochastic excitation,
as seen in solar oscillations. The random fluctuations in phase also support
this idea. A comparison with X Cam, a true Mira star with the same pulsation
period, and W Cyg, a true semiregular star, illustrates the basic differences
in phase behaviours. The Mira shows very stable phase, consistent with
excitation by the kappa-mechanism, whereas W Cyg shows large phase fluctuations
that imply stochastic excitation. We find L2 Pup to be intermediate, implying
that both mechanisms play a role in its pulsation. Finally, we also checked the
presence of low-dimensional chaos and could safely exclude it.

[58] oai:arXiv.org:astro-ph/0501202 [pdf] - 70314

When Darwin Met Einstein: Gravitational Lens Inversion with Genetic
Algorithms

Comments: 9 pages, to appear in PASA

Submitted: **2005-01-11**

Gravitational lensing can magnify a distant source, revealing structural
detail which is normally unresolvable. Recovering this detail through an
inversion of the influence of gravitational lensing, however, requires
optimisation of not only lens parameters, but also of the surface brightness
distribution of the source. This paper outlines a new approach to this
inversion, utilising genetic algorithms to reconstruct the source profile. In
this initial study, the effects of image degradation due to instrumental and
atmospheric effects are neglected and it is assumed that the lens model is
accurately known, but the genetic algorithm approach can be incorporated into
more general optimisation techniques, allowing the optimisation of both the
parameters for a lensing model and the surface brightness of the source.

[59] oai:arXiv.org:astro-ph/0410138 [pdf] - 67988

Probing Sub-parsec Structure in the Lyman Alpha Forest with
Gravitational Microlensing

Comments: 8 pages, 11 figures, to appear in MNRAS (note resolution of some
figures reduced due to size limitations)

Submitted: **2004-10-06**

We present the results of microlens ray-tracing simulations showing the
effect of absorbing material between a source quasar and a lensing galaxy in a
gravitational lens system. We find that, in addition to brightness fluctuations
due to microlensing, the strength of the absorption line relative to the
continuum varies with time, with the properties of the variations depending on
the structure of the absorbing material. We conclude that such variations will
be measurable via UV spectroscopy of image A of the gravitationally lensed
quasar Q2237+0305 if the Lyman Alpha clouds between the quasar and the lensing
galaxy possess structure on scales smaller than $\sim 0.1$ pc. The time scale
for the variations is on the order of order years to decades, although very
short term variability can occur. While the Lyman alpha lines may not be
accessible at all wavelengths, this approach is applicable to any absorption
system, including metal lines.