Search

Modelling complex deformation for soft robotics provides a guideline to understand their behaviour, leading to safe interaction with the environment. However, building a surrogate model with high accuracy and fast inference speed can be challenging for soft robotics due to the nonlinearity from complex geometry, large deformation, material nonlinearity etc. The reality gap from surrogate models also prevents their further deployment in the soft robotics domain. In this study, we proposed a physics-informed Neural Networks (PINNs) named PINN-Ray to model complex deformation for a Fin Ray soft robotic gripper, which embeds the minimum potential energy principle from elastic mechanics and additional high-fidelity experimental data into the loss function of neural network for training. This method is significant in terms of its generalisation to complex geometry and robust to data scarcity as compared to other data-driven neural networks. Furthermore, it has been extensively evaluated to model the deformation of the Fin Ray finger under external actuation. PINN-Ray demonstrates improved accuracy as compared with Finite element modelling (FEM) after applying the data assimilation scheme to treat the sim-to-real gap. Additionally, we introduced our automated framework to design, fabricate soft robotic fingers, and characterise their deformation by visual tracking, which provides a guideline for the fast prototype of soft robotics.

Isolated black holes (BHs) and neutron stars (NSs) are largely undetectable across the electromagnetic spectrum. For this reason, our only real prospect of observing these isolated compact remnants is via microlensing; a feat recently performed for the first time. However, characterisation of the microlensing events caused by BHs and NSs is still in its infancy. In this work, we perform N-body simulations to explore the frequency and physical characteristics of microlensing events across the entire sky. Our simulations find that every year we can expect $88_{-6}^{+6}$ BH, $6.8_{-1.6}^{+1.7}$ NS and $20^{+30}_{-20}$ stellar microlensing events which cause an astrometric shift larger than 2~mas. Similarly, we can expect $21_{-3}^{+3}$ BH, $18_{-3}^{+3}$ NS and $7500_{-500}^{+500}$ stellar microlensing events which cause a bump magnitude larger than 1~mag. Leveraging a more comprehensive dynamical model than prior work, we predict the fraction of microlensing events caused by BHs as a function of Einstein time to be smaller than previously thought. Comparison of our microlensing simulations to events in Gaia finds good agreement. Finally, we predict that in the combination of Gaia and GaiaNIR data there will be $14700_{-900}^{+600}$ BH and $1600_{-200}^{+300}$ NS events creating a centroid shift larger than 1~mas and $330_{-120}^{+100}$ BH and $310_{-100}^{+110}$ NS events causing bump magnitudes $> 1$. Of these, $<10$ BH and $5_{-5}^{+10}$ NS events should be detectable using current analysis techniques. These results inform future astrometric mission design, such as GaiaNIR, as they indicate that, compared to stellar events, there are fewer observable BH events than previously thought., Comment: Accepted for publication in MNRAS. 15 pages

The pyramid wavefront sensor (PyWFS) has become increasingly popular to use in adaptive optics (AO) systems due to its high sensitivity. The main drawback of the PyWFS is that it is inherently nonlinear, which means that classic linear wavefront reconstruction techniques face a significant reduction in performance at high wavefront errors, particularly when the pyramid is unmodulated. In this paper, we consider the potential use of neural networks (NNs) to replace the widely used matrix vector multiplication (MVM) control. We aim to test the hypothesis that the neural network (NN)'s ability to model nonlinearities will give it a distinct advantage over MVM control. We compare the performance of a MVM linear reconstructor against a dense NN, using daytime data acquired on the Subaru Coronagraphic Extreme Adaptive Optics system (SCExAO) instrument. In a first set of experiments, we produce wavefronts generated from 14 Zernike modes and the PyWFS responses at different modulation radii (25, 50, 75, and 100 mas). We find that the NN allows for a far more precise wavefront reconstruction at all modulations, with differences in performance increasing in the regime where the PyWFS nonlinearity becomes significant. In a second set of experiments, we generate a dataset of atmosphere-like wavefronts, and confirm that the NN outperforms the linear reconstructor. The SCExAO real-time computer software is used as baseline for the latter. These results suggest that NNs are well positioned to improve upon linear reconstructors and stand to bring about a leap forward in AO performance in the near future., Comment: Accepted in Publications of the Astronomical Society of the Pacific, September 2023

Photonic lanterns allow the decomposition of highly multimodal light into a simplified modal basis such as single-moded and/or few-moded. They are increasingly finding uses in astronomy, optics and telecommunications. Calculating propagation through a photonic lantern using traditional algorithms takes $\sim 1$ hour per simulation on a modern CPU. This paper demonstrates that neural networks can bridge the disparate opto-electronic systems, and when trained can achieve a speed-up of over 5 orders of magnitude. We show that this approach can be used to model photonic lanterns with manufacturing defects as well as successfully generalising to polychromatic data. We demonstrate two uses of these neural network models, propagating seeing through the photonic lantern as well as performing global optimisation for purposes such as photonic lantern funnels and photonic lantern nullers., Comment: 20 pages, 14 figures

Type Iax supernovae represent the largest class of peculiar white-dwarf supernovae. The type Iax SN~2012Z in NGC 1309 is the only white dwarf supernova with a detected progenitor system in pre-explosion observations. Deep \textit{Hubble Space Telescope} images taken before SN~2012Z show a luminous, blue source that we have interpreted as a helium-star companion (donor) to the exploding white dwarf. We present here late-time \textit{HST} observations taken $\sim$1400 days after the explosion to test this model. We find the SN light curve can empirically be fit by an exponential decay model in magnitude units. The fitted asymptotic brightness is within $10\%$ of our latest measurements and approximately twice the brightness of the pre-explosion source. The decline of the light curve is too slow to be powered by $^{56}$Co or $^{57}$Co decay: if radioactive decay is the dominate power source, it must be from longer half-life species like $^{55}$Fe. Interaction with circumstellar material may contribute to the light curve, as may shock heating of the companion star. Companion-star models underpredict the observed flux in the optical, producing most of their flux in the UV at these epochs. A radioactively-heated bound remnant, left after only a partial disruption of the white dwarf, is also capable of producing the observed excess late-time flux. Our analysis suggests that the total ejecta + remnant mass is consistent with the Chandrasekhar mass for a range of type Iax supernovae., Comment: 20 pages, 13 figures, accepted to ApJ

We present an overview of the SkyMapper optical follow-up program for gravitational-wave event triggers from the LIGO/Virgo observatories, which aims at identifying early GW170817-like kilonovae out to $\sim 200$ Mpc distance. We describe our robotic facility for rapid transient follow-up, which can target most of the sky at $\delta<+10\deg $ to a depth of $i_\mathrm{AB}\approx 20$ mag. We have implemented a new software pipeline to receive LIGO/Virgo alerts, schedule observations and examine the incoming real-time data stream for transient candidates. We adopt a real-bogus classifier using ensemble-based machine learning techniques, attaining high completeness ($\sim$98%) and purity ($\sim$91%) over our whole magnitude range. Applying further filtering to remove common image artefacts and known sources of transients, such as asteroids and variable stars, reduces the number of candidates by a factor of more than 10. We demonstrate the system performance with data obtained for GW190425, a binary neutron star merger detected during the LIGO/Virgo O3 observing campaign. In time for the LIGO/Virgo O4 run, we will have deeper reference images allowing transient detection to $i_\mathrm{AB}\approx $21 mag., Comment: 13 pages, 6 figures, Accepted for publication in PASA

Astrometric detection involves a precise measurement of stellar positions, and is widely regarded as the leading concept presently ready to find earth-mass planets in temperate orbits around nearby sun-like stars. The TOLIMAN space telescope[39] is a low-cost, agile mission concept dedicated to narrow-angle astrometric monitoring of bright binary stars. In particular the mission will be optimised to search for habitable-zone planets around Alpha Centauri AB. If the separation between these two stars can be monitored with sufficient precision, tiny perturbations due to the gravitational tug from an unseen planet can be witnessed and, given the configuration of the optical system, the scale of the shifts in the image plane are about one millionth of a pixel. Image registration at this level of precision has never been demonstrated (to our knowledge) in any setting within science. In this paper we demonstrate that a Deep Convolutional Auto-Encoder is able to retrieve such a signal from simplified simulations of the TOLIMAN data and we present the full experimental pipeline to recreate out experiments from the simulations to the signal analysis. In future works, all the more realistic sources of noise and systematic effects present in the real-world system will be injected into the simulations., Comment: Preprint version of the manuscript to appear in the Volume "Intelligent Astrophysics" of the series "Emergence, Complexity and Computation", Book eds. I. Zelinka, D. Baron, M. Brescia, Springer Nature Switzerland, ISSN: 2194-7287

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., Comment: 23 pages, 12 figures, accepted for MNRAS

The rigorous quantification of uncertainty in geophysical inversions is a challenging problem. Inversions are often ill-posed and the likelihood surface may be multi-modal; properties of any single mode become inadequate uncertainty measures, and sampling methods become inefficient for irregular posteriors or high-dimensional parameter spaces. We explore the influences of different choices made by the practitioner on the efficiency and accuracy of Bayesian geophysical inversion methods that rely on Markov chain Monte Carlo sampling to assess uncertainty, using a multi-sensor inversion of the three-dimensional structure and composition of a region in the Cooper Basin of South Australia as a case study. The inversion is performed using an updated version of the Obsidian distributed inversion software. We find that the posterior for this inversion has complex local covariance structure, hindering the efficiency of adaptive sampling methods that adjust the proposal based on the chain history. Within the context of a parallel-tempered Markov chain Monte Carlo scheme for exploring high-dimensional multi-modal posteriors, a preconditioned Crank-Nicholson proposal outperforms more conventional forms of random walk. Aspects of the problem setup, such as priors on petrophysics or on 3-D geological structure, affect the shape and separation of posterior modes, influencing sampling performance as well as the inversion results. Use of uninformative priors on sensor noise can improve inversion results by enabling optimal weighting among multiple sensors even if noise levels are uncertain. Efficiency could be further increased by using posterior gradient information within proposals, which Obsidian does not currently support, but which could be emulated using posterior surrogates., Comment: 17 pages, 5 figures, 1 table; submitted to Geoscientific Model Development

We present bolometric light curves constructed from multi-wavelength photometry of Type Ia supernovae (SNe Ia) from the Carnegie Supernova Project and the CfA Supernova Group, using near-infrared observations to provide robust constraints on host galaxy dust extinction. This set of light curves form a well-measured reference set for comparison with theoretical models. Ejected mass and synthesized $^{56}$Ni mass are inferred for each SN Ia from its bolometric light curve using a semi-analytic Bayesian light curve model, and fitting formulae provided in terms of light curve width parameters from the SALT2 and SNooPy light curve fitters. A weak bolometric width-luminosity relation is confirmed, along with a correlation between ejected mass and the bolometric light curve width. SNe Ia likely to have sub-Chandrasekhar ejected masses belong preferentially to the broad-line and cool-photosphere spectroscopic subtypes, and have higher photospheric velocities and populate older, higher-mass host galaxies than SNe Ia consistent with Chandrasekhar-mass explosions. Two peculiar events, SN 2006bt and SN 2006ot, have normal peak luminosities but appear to have super-Chandrasekhar ejected masses., Comment: 26 pages, 14 figures; accepted to MNRAS. An online-only appendix in the MNRAS version is included as a supplemental appendix to the arXiv text; online-only tables, including bolometric light curves and MCMC inversion results, are included as ancillary files

The extraction of geological lineaments from digital satellite data is a fundamental application in remote sensing. The location of geological lineaments such as faults and dykes are of interest for a range of applications, particularly because of their association with hydrothermal mineralization. Although a wide range of applications have utilized computer vision techniques, a standard workflow for application of these techniques to mineral exploration is lacking. We present a framework for extracting geological lineaments using computer vision techniques which is a combination of edge detection and line extraction algorithms for extracting geological lineaments using optical remote sensing data. It features ancillary computer vision techniques for reducing data dimensionality, removing noise and enhancing the expression of lineaments. We test the proposed framework on Landsat 8 data of a mineral-rich portion of the Gascoyne Province in Western Australia using different dimension reduction techniques and convolutional filters. To validate the results, the extracted lineaments are compared to our manual photointerpretation and geologically mapped structures by the Geological Survey of Western Australia (GSWA). The results show that the best correlation between our extracted geological lineaments and the GSWA geological lineament map is achieved by applying a minimum noise fraction transformation and a Laplacian filter. Application of a directional filter instead shows a stronger correlation with the output of our manual photointerpretation and known sites of hydrothermal mineralization. Hence, our framework using either filter can be used for mineral prospectivity mapping in other regions where faults are exposed and observable in optical remote sensing data., Comment: 17 pages, 10 figures, 2 tables

Estimating the impact of environmental processes on vertical reef development in geological time is a very challenging task. pyReef-Core is a deterministic carbonate stratigraphic forward model designed to simulate the key biological and environmental processes that determine vertical reef accretion and assemblage changes in fossil reef drill cores. We present a Bayesian framework called Bayesreef for the estimation and uncertainty quantification of parameters in pyReef-Core that represent environmental conditions affecting the growth of coral assemblages on geological timescales. We demonstrate the existence of multimodal posterior distributions and investigate the challenges of sampling using Markov chain Monte-Carlo (MCMC) methods, which includes parallel tempering MCMC. We use synthetic reef-core to investigate fundamental issues and then apply the methodology to a selected reef-core from the Great Barrier Reef in Australia. The results show that Bayesreef accurately estimates and provides uncertainty quantification of the selected parameters that represent the environment and ecological conditions in pyReef-Core. Bayesreef provides insights into the complex posterior distributions of parameters in pyReef-Core, which provides the groundwork for future research in this area.

We present detailed ultraviolet, optical and near-infrared light curves of the Type Ia supernova (SN) 2012fr, which exploded in the Fornax cluster member NGC 1365. These precise high-cadence light curves provide a dense coverage of the flux evolution from $-$12 to $+$140 days with respect to the epoch of $B$-band maximum (\tmax). Supplementary imaging at the earliest epochs reveals an initial slow, nearly linear rise in luminosity with a duration of $\sim$2.5 days, followed by a faster rising phase that is well reproduced by an explosion model with a moderate amount of $^{56}$Ni mixing in the ejecta. From an analysis of the light curves, we conclude: $(i)$ explosion occurred $< 22$ hours before the first detection of the supernova, $(ii)$ the rise time to peak bolometric ($\lambda > 1800 $\AA) luminosity was $16.5 \pm 0.6$ days, $(iii)$ the supernova suffered little or no host-galaxy dust reddening, $(iv)$ the peak luminosity in both the optical and near-infrared was consistent with the bright end of normal Type Ia diversity, and $(v)$ $0.60 \pm 0.15 M_{\odot}$ of $^{56}$Ni was synthesized in the explosion. Despite its normal luminosity, SN 2012fr displayed unusually prevalent high-velocity \ion{Ca}{2} and \ion{Si}{2} absorption features, and a nearly constant photospheric velocity of the \ion{Si}{2} $\lambda$6355 line at $\sim$12,000 \kms\ beginning $\sim$5 days before \tmax. Other peculiarities in the early phase photometry and the spectral evolution are highlighted. SN 2012fr also adds to a growing number of Type Ia supernovae hosted by galaxies with direct Cepheid distance measurements., Comment: 88 pages, 29 figures, 12 tables

We present the first data release (DR1) of the SkyMapper Southern Survey, a hemispheric survey carried out with the SkyMapper Telescope at Siding Spring Observatory in Australia. Here, we present the survey strategy, data processing, catalogue construction and database schema. The DR1 dataset includes over 66,000 images from the Shallow Survey component, covering an area of 17,200 deg$^2$ in all six SkyMapper passbands $uvgriz$, while the full area covered by any passband exceeds 20,000 deg$^2$. The catalogues contain over 285 million unique astrophysical objects, complete to roughly 18 mag in all bands. We compare our $griz$ point-source photometry with PanSTARRS1 DR1 and note an RMS scatter of 2%. The internal reproducibility of SkyMapper photometry is on the order of 1%. Astrometric precision is better than 0.2 arcsec based on comparison with Gaia DR1. We describe the end-user database, through which data are presented to the world community, and provide some illustrative science queries., Comment: 31 pages, 19 figures, 10 tables, PASA, accepted

Current synoptic sky surveys monitor large areas of the sky to find variable and transient astronomical sources. As the number of detections per night at a single telescope easily exceeds several thousand, current detection pipelines make intensive use of machine learning algorithms to classify the detected objects and to filter out the most interesting candidates. A number of upcoming surveys will produce up to three orders of magnitude more data, which renders high-precision classification systems essential to reduce the manual and, hence, expensive vetting by human experts. We present an approach based on convolutional neural networks to discriminate between true astrophysical sources and artefacts in reference-subtracted optical images. We show that relatively simple networks are already competitive with state-of-the-art systems and that their quality can further be improved via slightly deeper networks and additional preprocessing steps -- eventually yielding models outperforming state-of-the-art systems. In particular, our best model correctly classifies about 97.3% of all 'real' and 99.7% of all 'bogus' instances on a test set containing 1,942 'bogus' and 227 'real' instances in total. Furthermore, the networks considered in this work can also successfully classify these objects at hand without relying on difference images, which might pave the way for future detection pipelines not containing image subtraction steps at all., Comment: 15 pages, 15 figures, 3 tables

The SkyMapper 1.3 m telescope at Siding Spring Observatory has now begun regular operations. Alongside the Southern Sky Survey, a comprehensive digital survey of the entire southern sky, SkyMapper will carry out a search for supernovae and other transients. The search strategy, covering a total footprint area of ~2000 deg2 with a cadence of $\leq 5$ days, is optimised for discovery and follow-up of low-redshift type Ia supernovae to constrain cosmic expansion and peculiar velocities. We describe the search operations and infrastructure, including a parallelised software pipeline to discover variable objects in difference imaging; simulations of the performance of the survey over its lifetime; public access to discovered transients; and some first results from the Science Verification data., Comment: 13 pages, 11 figures; submitted to PASA

This paper presents the first major data release and survey description for the ANU WiFeS SuperNovA Program (AWSNAP). AWSNAP is an ongoing supernova spectroscopy campaign utilising the Wide Field Spectrograph (WiFeS) on the Australian National University (ANU) 2.3m telescope. The first and primary data release of this program (AWSNAP-DR1) releases 357 spectra of 175 unique objects collected over 82 equivalent full nights of observing from July 2012 to August 2015. These spectra have been made publicly available via the WISeREP supernova spectroscopy repository. We analyse the AWSNAP sample of Type Ia supernova spectra, including measurements of narrow sodium absorption features afforded by the high spectral resolution of the WiFeS instrument. In some cases we were able to use the integral-field nature of the WiFeS instrument to measure the rotation velocity of the SN host galaxy near the SN location in order to obtain precision sodium absorption velocities. We also present an extensive time series of SN 2012dn, including a near-nebular spectrum which both confirms its "super-Chandrasekhar" status and enables measurement of the sub-solar host metallicity at the SN site., Comment: Submitted to Publications of the Astronomical Society of Australia (PASA). Spectra publicly released via WISeREP at http://wiserep.weizmann.ac.il/

The light curves of Type Ia supernovae (SNe Ia) are powered by the radioactive decay of $^{56}$Ni to $^{56}$Co at early times, and the decay of $^{56}$Co to $^{56}$Fe from ~60 days after explosion. We examine the evolution of the [Co III] 5892 A emission complex during the nebular phase for SNe Ia with multiple nebular spectra and show that the line flux follows the square of the mass of $^{56}$Co as a function of time. This result indicates both efficient local energy deposition from positrons produced in $^{56}$Co decay, and long-term stability of the ionization state of the nebula. We compile 77 nebular spectra of 25 SN Ia from the literature and present 17 new nebular spectra of 7 SNe Ia, including SN2014J. From these we measure the flux in the [Co III] 5892 A line and remove its well-behaved time dependence to infer the initial mass of $^{56}$Ni ($M_{Ni}$) produced in the explosion. We then examine $^{56}$Ni yields for different SN Ia ejected masses ($M_{ej}$ - calculated using the relation between light curve width and ejected mass) and find the $^{56}$Ni masses of SNe Ia fall into two regimes: for narrow light curves (low stretch s~0.7-0.9), $M_{Ni}$ is clustered near $M_{Ni}$ ~ 0.4$M_\odot$ and shows a shallow increase as $M_{ej}$ increases from ~1-1.4$M_\odot$; at high stretch, $M_{ej}$ clusters at the Chandrasekhar mass (1.4$M_\odot$) while $M_{Ni}$ spans a broad range from 0.6-1.2$M_\odot$. This could constitute evidence for two distinct SN Ia explosion mechanisms., Comment: 16 pages, 12 figures (main text), plus data tables in appendix. Spectra released on WISeREP. Submitted to MNRAS, comments welcome

Among Type Ia supernovae (SNe~Ia) exist a class of overluminous objects whose ejecta mass is inferred to be larger than the canonical Chandrasekhar mass. We present and discuss the UV/optical photometric light curves, colors, absolute magnitudes, and spectra of three candidate Super-Chandrasekhar mass SNe--2009dc, 2011aa, and 2012dn--observed with the Swift Ultraviolet/Optical Telescope. The light curves are at the broad end for SNe Ia, with the light curves of SN~2011aa being amongst the broadest ever observed. We find all three to have very blue colors which may provide a means of excluding these overluminous SNe from cosmological analysis, though there is some overlap with the bluest of "normal" SNe Ia. All three are overluminous in their UV absolute magnitudes compared to normal and broad SNe Ia, but SNe 2011aa and 2012dn are not optically overluminous compared to normal SNe Ia. The integrated luminosity curves of SNe 2011aa and 2012dn in the UVOT range (1600-6000 Angstroms) are only half as bright as SN~2009dc, implying a smaller 56Ni yield. While not enough to strongly affect the bolometric flux, the early time mid-UV flux makes a significant contribution at early times. The strong spectral features in the mid-UV spectra of SNe 2009dc and 2012dn suggest a higher temperature and lower opacity to be the cause of the UV excess rather than a hot, smooth blackbody from shock interaction. Further work is needed to determine the ejecta and 56Ni masses of SNe 2011aa and 2012dn and fully explain their high UV luminosities., Comment: 12 pages, 8 figures Accepted for publication in the Astrophysical Journal Data available on the Swift supernova website http://swift.gsfc.nasa.gov/docs/swift/sne/swift_sn.html

We present a first systematic comparison of superluminous Type Ia supernovae (SNe Ia) at late epochs, including previously unpublished photometric and spectroscopic observations of SN 2007if, SN 2009dc and SNF20080723-012. Photometrically, the objects of our sample show a diverse late-time behaviour, some of them fading quite rapidly after a light-curve break at ~150-200d. The latter is likely the result of flux redistribution into the infrared, possibly caused by dust formation, rather than a true bolometric effect. Nebular spectra of superluminous SNe Ia are characterised by weak or absent [Fe III] emission, pointing at a low ejecta ionisation state as a result of high densities. To constrain the ejecta and 56Ni masses of superluminous SNe Ia, we compare the observed bolometric light curve of SN 2009dc with synthetic model light curves, focusing on the radioactive tail after ~60d. Models with enough 56Ni to explain the light-curve peak by radioactive decay, and at the same time sufficient mass to keep the ejecta velocities low, fail to reproduce the observed light-curve tail of SN 2009dc because of too much gamma-ray trapping. We instead propose a model with ~1 solar mass of 56Ni and ~2 solar masses of ejecta, which may be interpreted as the explosion of a Chandrasekhar-mass white dwarf (WD) enshrouded by 0.6-0.7 solar masses of C/O-rich material, as it could result from a merger of two massive C/O WDs. This model reproduces the late light curve of SN 2009dc well. A flux deficit at peak may be compensated by light from the interaction of the ejecta with the surrounding material., Comment: 16 pages, 10 figures, accepted for publication in MNRAS; references updated and typos corrected in v2

We put constraints on the properties of the progenitors of peculiar calcium-rich transients using the distribution of locations within their host galaxies. We confirm that this class of transients do not follow the galaxy stellar mass profile and are more likely to be found in remote locations of their apparent hosts. We test the hypothesis that these transients are from low metallicity progenitors by comparing their spatial distributions with the predictions of self-consistent cosmological simulations that include star formation and chemical enrichment. We find that while metal-poor stars and our transient sample show a consistent preference for large offsets, metallicity alone cannot explain the extreme cases. Invoking a lower age limit on the progenitor helps to improve the match, indicating these events may result from a very old metal-poor population. We also investigate the radial distribution of globular cluster systems, and show that they too are consistent with the class of calcium-rich transients. Because photometric upper limits exist for globular clusters for some members of the class, a production mechanism related to the dense environment of globular clusters is not favoured for the calcium-rich events. However the methods developed in this paper may be used in the future to constrain the effects of low metallicity on radially distant core-collapse events or help establish a correlation with globular clusters for other classes of peculiar explosions., Comment: 8 pages, 4 figures, accepted for publication in MNRAS

Accurate measurement of galaxy cluster masses is an essential component not only in studies of cluster physics, but also for probes of cosmology. However, different mass measurement techniques frequently yield discrepant results. The SDSS MaxBCG catalog's mass-richness relation has previously been constrained using weak lensing shear, Sunyaev-Zeldovich (SZ), and X-ray measurements. The mass normalization of the clusters as measured by weak lensing shear is >~25% higher than that measured using SZ and X-ray methods, a difference much larger than the stated measurement errors in the analyses. We constrain the mass-richness relation of the MaxBCG galaxy cluster catalog by measuring the gravitational lensing magnification of type I quasars in the background of the clusters. The magnification is determined using the quasars' variability and the correlation between quasars' variability amplitude and intrinsic luminosity. The mass-richness relation determined through magnification is in agreement with that measured using shear, confirming that the lensing strength of the clusters implies a high mass normalization, and that the discrepancy with other methods is not due to a shear-related systematic measurement error. We study the dependence of the measured mass normalization on the cluster halo orientation. As expected, line-of-sight clusters yield a higher normalization; however, this minority of haloes does not significantly bias the average mass-richness relation of the catalog., Comment: 9 pages. Accepted for publication in ApJ

We apply the reconstruction technique to the clustering of galaxies from the SDSS DR7 LRG sample, sharpening the baryon acoustic oscillation (BAO) feature and achieving a 1.9% measurement of the distance to z=0.35. This is the first application of reconstruction of the BAO feature in a galaxy redshift survey. We update the reconstruction algorithm of Eisenstein et al, 2007 to account for the effects of survey geometry as well as redshift-space distortions and validate it on 160 LasDamas simulations. We demonstrate that reconstruction sharpens the BAO feature in the angle averaged galaxy correlation function, reducing the nonlinear smoothing scale \Sigma_nl from 8.1 Mpc/h to 4.4 Mpc/h. Reconstruction also significantly reduces the effects of redshift-space distortions at the BAO scale, isotropizing the correlation function. This sharpened BAO feature yields an unbiased distance estimate (< 0.2%) and reduces the scatter from 3.3% to 2.1%. We demonstrate the robustness of these results to the various reconstruction parameters, including the smoothing scale, the galaxy bias and the linear growth rate. Applying this reconstruction algorithm to the SDSS LRG DR7 sample improves the significance of the BAO feature in these data from 3.3 sigma for the unreconstructed correlation function, to 4.2 sigma after reconstruction. We estimate a relative distance scale D_V/r_s to z=0.35 of 8.88+/-0.17, where r_s is the sound horizon and D_V = (D_A^2/H)^{1/3} is a combination of the angular diameter distance D_A and Hubble parameter H. Assuming a sound horizon of 154.25 Mpc, this translates into a distance measurement D_V (z=0.35) = 1.356+/-0.025 Gpc. We find that reconstruction reduces the distance error in the DR7 sample from 3.5% to 1.9%, equivalent to a survey with three times the volume of SDSS., Comment: 15 pages, 15 figures, MNRAS submitted. This is the first of a series of three papers. See http://www.astro.yale.edu/padmanabhan/dr7recon/ for associated data files and figures

We present ultraviolet (UV) spectroscopy and photometry of four Type Ia supernovae (SNe 2004dt, 2004ef, 2005M, and 2005cf) obtained with the UV prism of the Advanced Camera for Surveys on the Hubble Space Telescope. This dataset provides unique spectral time series down to 2000 Angstrom. Significant diversity is seen in the near maximum-light spectra (~ 2000--3500 Angstrom) for this small sample. The corresponding photometric data, together with archival data from Swift Ultraviolet/Optical Telescope observations, provide further evidence of increased dispersion in the UV emission with respect to the optical. The peak luminosities measured in uvw1/F250W are found to correlate with the B-band light-curve shape parameter dm15(B), but with much larger scatter relative to the correlation in the broad-band B band (e.g., ~0.4 mag versus ~0.2 mag for those with 0.8 < dm15 < 1.7 mag). SN 2004dt is found as an outlier of this correlation (at > 3 sigma), being brighter than normal SNe Ia such as SN 2005cf by ~0.9 mag and ~2.0 mag in the uvw1/F250W and uvm2/F220W filters, respectively. We show that different progenitor metallicity or line-expansion velocities alone cannot explain such a large discrepancy. Viewing-angle effects, such as due to an asymmetric explosion, may have a significant influence on the flux emitted in the UV region. Detailed modeling is needed to disentangle and quantify the above effects., Comment: 17 pages, 13 figures, accepted by ApJ

We introduce a technique to measure gravitational lensing magnification using the variability of type I quasars. Quasars' variability amplitudes and luminosities are tightly correlated, on average. Magnification due to gravitational lensing increases the quasars' apparent luminosity, while leaving the variability amplitude unchanged. Therefore, the mean magnification of an ensemble of quasars can be measured through the mean shift in the variability-luminosity relation. As a proof of principle, we use this technique to measure the magnification of quasars spectroscopically identified in the Sloan Digital Sky Survey, due to gravitational lensing by galaxy clusters in the SDSS MaxBCG catalog. The Palomar-QUEST Variability Survey, reduced using the DeepSky pipeline, provides variability data for the sources. We measure the average quasar magnification as a function of scaled distance (r/R200) from the nearest cluster; our measurements are consistent with expectations assuming NFW cluster profiles, particularly after accounting for the known uncertainty in the clusters' centers. Variability-based lensing measurements are a valuable complement to shape-based techniques because their systematic errors are very different, and also because the variability measurements are amenable to photometric errors of a few percent and to depths seen in current wide-field surveys. Given the data volume expected from current and upcoming surveys, this new technique has the potential to be competitive with weak lensing shear measurements of large scale structure., Comment: Accepted for publication in ApJ

We identify 3,113 highly variable objects in 7,200 square degrees of the Palomar-QUEST Survey, which each varied by more than 0.4 magnitudes simultaneously in two broadband optical filters on timescales from hours to roughly 3.5 years. The primary goal of the selection is to find blazars by their well-known violent optical variability. Because most known blazars have been found in radio and/or X-ray wavelengths, a sample discovered through optical variability may have very different selection effects, elucidating the range of behavior possible in these systems. A set of blazars selected in this unusual manner will improve our understanding of the physics behind this extremely variable and diverse class of AGN. The object positions, variability statistics, and color information are available using the Palomar-QUEST CasJobs server. The time domain is just beginning to be explored over large sky areas; we do not know exactly what a violently variable sample will hold. About 20% of the sample has been classified in the literature; over 70% of those objects are known or likely AGN. The remainder largely consists of a variety of variable stars, including a number of RR Lyrae and cataclysmic variables., Comment: 22 pages (preprint format), 2 figures. Accepted for publication in ApJ. References updated

We study the ensemble optical variability of 276 FSRQs and 86 BL Lacs in the Palomar-QUEST Survey with the goal of searching for common fluctuation properties, examining the range of behavior across the sample, and characterizing the appearance of blazars in such a survey so that future work can more easily identify such objects. The survey, which covers 15,000 square degrees multiple times over 3.5 years, allows for the first ensemble blazar study of this scale. Variability amplitude distributions are shown for the FSRQ and BL Lac samples for numerous time lags, and also studied through structure function analyses. Individual blazars show a wide range of variability amplitudes, timescales, and duty cycles. Of the best sampled objects, 35% are seen to vary by more than 0.4 magnitudes; for these, the fraction of measurements contributing to the high amplitude variability ranges constantly from about 5% to 80%. Blazar variability has some similarities to that of type I quasars but includes larger amplitude fluctuations on all timescales. FSRQ variability amplitudes are particularly similar to those of QSOs on timescales of several months, suggesting significant contributions from the accretion disk to the variable flux at these timescales. Optical variability amplitudes are correlated with the maximum apparent velocities of the radio jet for the subset of FSRQs with MOJAVE VLBA measurements, implying that the optically variable flux's strength is typically related to that of the radio emission. We also study CRATES radio-selected FSRQ candidates, which show similar variability characteristics to known FSRQs; this suggests a high purity for the CRATES sample., Comment: 29 pages, 12 figures. Accepted for publication in ApJ

The ensemble variability properties of nearly 23,000 quasars are studied using the Palomar-QUEST Survey. The survey has covered 15,000 square degrees multiple times over 3.5 years using 7 optical filters, and has been calibrated specifically for variability work. Palomar-QUEST allows for the study of rare objects using multiple epochs of consistently calibrated, homogeneous data, obviating the common problem of generating comparable measurements from disparate datasets. A power law fit to the quasar structure function versus time yields an index of 0.432 +/- 0.024 for our best measured sample. We see the commonly reported anticorrelation between average optical variability amplitude and optical luminosity, and measure the logarithmic decrease in variability amplitude to scale as the logarithm of the luminosity times 0.205 +/- 0.002. Black hole mass is positively correlated with variability amplitude over three orders of magnitude in mass. Quasar variability amplitude is seen to decrease with Eddington ratio as a step function with a transition around Eddington ratio of 0.5. The higher variability at low Eddington ratios is due to excess power at timescales shorter than roughly 300 days. X-ray and radio measurements exist for subsets of the quasar sample. We observe an anticorrelation between optical variability amplitude and X-ray luminosity. No significant correlation is seen between average optical variability properties and radio luminosity. The timescales of quasar fluctuations are suggestive of accretion disk instabilities. The relationships seen between variability, Eddington ratio, and radio and X-ray emission are discussed in terms of a possible link between the behavior of quasars and black hole X-ray binaries., Comment: 25 pages, 13 figures, accepted for publication in ApJ

In the companion paper we began the task of systematically studying the detection of planets in wide orbits ($a > 1.5 R_E$) via microlensing surveys. In this paper we continue, focusing on repeating events. We find that, if all planetary systems are similar to our own Solar System, reasonable extensions of the present observing strategies would allow us to detect 3-6 repeating events per year along the direction to the Bulge. Indeed, if planetary systems with multiple planets are common, then future monitoring programs which lead to the discovery of thousands of stellar-lens events will likely discover events in which several different planets within a single system serve as lenses, with light curves exhibiting multiple repetitions. In this paper we discuss observing strategies to maximize the discovery of all wide-orbit planet-lens events. We also compare the likely detection rates of planets in wide orbits to those of planets located in the zone for resonant lensing. We find that, depending on the values of the planet masses and stellar radii of the lensed sources (which determine whether or not finite source size is important), and also on the sensitivity of the photometry used by observers, the detection of planets in wide orbits may be the primary route to the discovery of planets via microlensing. We also discuss how the combination of resonant and wide-orbit events can help us to learn about the distribution of planetary system properties (S 6.1). In addition, by determining the fraction of short-duration events due to planets, we indirectly derive information about the fraction of all short-duration events that may be due to low-mass MACHOs (S 6.2)., Comment: 51 pages, 7 figures. To be published in the Astrophysical Journal, 20 February 1999. This completes the introduction to the discovery of planets in wide orbits begun in astro-ph/9808075, also to appear in ApJ on 20 February 1999

We propose and evaluate the feasibility of a new strategy to search for planets via microlensing. This new strategy is designed to detect planets in "wide" orbits, i.e., with orbital separation, $a$ greater than $\sim 1.5 R_E$. Planets in wide orbits may provide the dominant channel for the microlensing discovery of planets, particularly low-mass (e.g., Earth-mass) planets. This paper concentrates on events in which a single planet serves as a lens, leading to an isolated event of short duration. We point out that a distribution of events due to lensing by stars with wide-orbit planets is necessarily accompanied by a distribution of shorter- duration events. The fraction of events in the latter distribution is proportional to the average value of $\sqrt{q}$, where $q$ is the ratio between \pl and stellar masses. The position of the peak or peaks also provides a measure of the mass ratios typical of planetary systems. We study detection strategies that can optimize our ability to discover isolated short-duration events due to lensing by planets, and find that monitoring employing sensitive photometry is particularly useful. If planetary systems similar to our own are common, even modest changes in detection strategy should lead to the discovery of a few isolated events of short duration every year. We therefore also address the issue of the contamination due to stellar populations of any microlensing signal due to low-mass MACHOs. We describe how, even for isolated events of short duration, it will be possible to test the hypothesis that the lens was a planet instead of a low-mass MACHO, if the central star of the planetary system contributes a measurable fraction of the baseline flux., Comment: 37 pages, 6 figure. To be published in the Astrophysical Journal. This is part one of a series of papers on microlensing by planetary systems containing wide-orbit planets; the series represents a reorganization and extension of astro-ph/9711013

We propose and evaluate the feasibility of a new strategy to search for planets via microlensing observations. This new strategy is designed to detect planets in "wide" orbits, i.e., with orbital separation, a, greater than ~1.5 R_E. Planets in wide orbits may provide the dominant channel for the discovery of planets via microlensing, particularly low-mass (e.g., Earth-mass) planets. Because the ongoing microlensing observations and extensions of them should be able to discover planets in wide orbits, we provide a foundation for the search through detailed calculations and simulations that quantify the expected results and compare the relative benefits of various search strategies. If planetary systems similar to our own or to some of the known extra-solar systems are common, then the predicted detection rates of wide-orbit events are high, generally in the range 2-10% of the present detection rate for apparently single events by stars. The expected high rates should allow the microlensing observing teams to either place significant limits on the presence of planetary systems in the Galactic Bulge, or begin to probe the population in detail within the next few years. We also address the issues of (1) whether planets discovered via microlensing are likely to harbor life, (2) the feasibility of follow-up observations to learn more about planet microlenses, and (3) the contamination due to stellar populations of any microlensing signal due to low-mass MACHOs., Comment: 72 pages, 16 figures, 4 tables, incorporated into a single file. Submitted to the Astrophysical Journal, email: rdistefano@cfa.harvard.edu