Showing 1–50 of 86 results for author: Ward, S

AGN-driven outflows in clumpy media: multiphase structure and scaling relations

Authors: Samuel Ruthven Ward, Tiago Costa, Chris M. Harrison, Vincenzo Mainieri

Abstract: Small-scale winds driven from accretion discs surrounding active galactic nuclei (AGN) are expected to launch kpc-scale outflows into their host galaxies. However, the ways in which the structure of the interstellar medium (ISM) affects the multiphase content and impact of the outflow remains uncertain. We present a series of numerical experiments featuring a realistic small-scale AGN wind with ve… ▽ More Small-scale winds driven from accretion discs surrounding active galactic nuclei (AGN) are expected to launch kpc-scale outflows into their host galaxies. However, the ways in which the structure of the interstellar medium (ISM) affects the multiphase content and impact of the outflow remains uncertain. We present a series of numerical experiments featuring a realistic small-scale AGN wind with velocity $5\times 10^3-10^4\ \rm{km/s}$ interacting with an isolated galaxy disc with a manually-controlled clumpy ISM, followed at sub-pc resolution. Our simulations are performed with AREPO and probe a wide range of AGN luminosities ($L=10^{43-47}\ \rm{erg/s}$) and ISM substructures. In homogeneous discs, the AGN wind sweeps up an outflowing, cooling shell, where the emerging cold phase dominates the mass and kinetic energy budgets, reaching a momentum flux $\dot{p} \approx 7\ L/c$. However, when the ISM is clumpy, outflow properties are profoundly different. They contain small, long-lived ($> 5\ \rm{Myr}$), cold ($T<10^{4.5}\ \rm{K}$) cloudlets entrained in the faster, hot outflow phase, which are only present in the outflow if radiative cooling is included in the simulation. While the cold phase dominates the mass of the outflow, most of the kinetic luminosity is now carried by a tenuous, hot phase with $T > 10^7 \ \rm K$. While the hot phases reaches momentum fluxes $\dot{p} \approx (1 - 5)\ L/c$, energy-driven bubbles couple to the cold phase inefficiently, producing modest momentum fluxes $\dot{p} < L/c$ in the fast-outflowing cold gas. These low momentum fluxes could lead to the outflows being misclassified as momentum-driven using common observational diagnostics. We also show predictions for scaling relations between outflow properties and AGN luminosity and discuss the challenges in constraining outflow driving mechanisms and kinetic coupling efficiencies using observed quantities. △ Less

Submitted 24 July, 2024; originally announced July 2024.

Comments: Accepted to MNRAS. 23 pages, 14 plots, 1 table

Sensor Response-Time Reduction using Long-Short Term Memory Network Forecasting

Authors: Simon J. Ward, Muhamed Baljevic, Sharon M. Weiss

Abstract: The response time of a biosensor is a crucial metric in safety-critical applications such as medical diagnostics where an earlier diagnosis can markedly improve patient outcomes. However, the speed at which a biosensor reaches a final equilibrium state can be limited by poor mass transport and long molecular diffusion times that increase the time it takes target molecules to reach the active sensi… ▽ More The response time of a biosensor is a crucial metric in safety-critical applications such as medical diagnostics where an earlier diagnosis can markedly improve patient outcomes. However, the speed at which a biosensor reaches a final equilibrium state can be limited by poor mass transport and long molecular diffusion times that increase the time it takes target molecules to reach the active sensing region of a biosensor. While optimization of system and sensor design can promote molecules reaching the sensing element faster, a simpler and complementary approach for response time reduction that is widely applicable across all sensor platforms is to use time-series forecasting to predict the ultimate steady-state sensor response. In this work, we show that ensembles of long short-term memory (LSTM) networks can accurately predict equilibrium biosensor response from a small quantity of initial time-dependent biosensor measurements, allowing for significant reduction in response time by a mean and median factor of improvement of 18.6 and 5.1, respectively. The ensemble of models also provides simultaneous estimation of uncertainty, which is vital to provide confidence in the predictions and subsequent safety-related decisions that are made. This approach is demonstrated on real-time experimental data collected by exposing porous silicon biosensors to buffered protein solutions using a multi-channel fluidic cell that enables the automated measurement of 100 porous silicon biosensors in parallel. The dramatic improvement in sensor response time achieved using LSTM network ensembles and associated uncertainty quantification opens the door to trustworthy and faster responding biosensors, enabling more rapid medical diagnostics for improved patient outcomes and healthcare access, as well as quicker identification of toxins in food and the environment. △ Less

Submitted 26 April, 2024; originally announced April 2024.

Comments: 9 pages, 3 figures

Scalable hierarchical BayeSN inference: Investigating dependence of SN Ia host galaxy dust properties on stellar mass and redshift

Authors: Matthew Grayling, Stephen Thorp, Kaisey S. Mandel, Suhail Dhawan, Ana Sofia M. Uzsoy, Benjamin M. Boyd, Erin E. Hayes, Sam M. Ward

Abstract: We apply the hierarchical probabilistic SED model BayeSN to analyse a sample of 475 SNe Ia (0.015 < z < 0.4) from Foundation, DES3YR and PS1MD to investigate the properties of dust in their host galaxies. We jointly infer the dust law $R_V$ population distributions at the SED level in high- and low-mass galaxies simultaneously with dust-independent, intrinsic differences. We find an intrinsic mass… ▽ More We apply the hierarchical probabilistic SED model BayeSN to analyse a sample of 475 SNe Ia (0.015 < z < 0.4) from Foundation, DES3YR and PS1MD to investigate the properties of dust in their host galaxies. We jointly infer the dust law $R_V$ population distributions at the SED level in high- and low-mass galaxies simultaneously with dust-independent, intrinsic differences. We find an intrinsic mass step of $-0.049\pm0.016$ mag, at a significance of 3.1$σ$, when allowing for a constant intrinsic, achromatic magnitude offset. We additionally apply a model allowing for time- and wavelength-dependent intrinsic differences between SNe Ia in different mass bins, finding $\sim$2$σ$ differences in magnitude and colour around peak and 4.5$σ$ differences at later times. These intrinsic differences are inferred simultaneously with a difference in population mean $R_V$ of $\sim$2$σ$ significance, demonstrating that both intrinsic and extrinsic differences may play a role in causing the host galaxy mass step. We also consider a model which allows the mean of the $R_V$ distribution to linearly evolve with redshift but find no evidence for any evolution - we infer the gradient of this relation $η_R = -0.38\pm0.70$. In addition, we discuss in brief a new, GPU-accelerated Python implementation of BayeSN suitable for application to large surveys which is publicly available and can be used for future cosmological analyses; this code can be found here: https://github.com/bayesn/bayesn. △ Less

Submitted 29 April, 2024; v1 submitted 16 January, 2024; originally announced January 2024.

Comments: 24 pages, 8 figures, 3 tables. Accepted for publication in MNRAS. BayeSN code available at https://github.com/bayesn/bayesn

At the end of cosmic noon: Short gas depletion times in unobscured quasars at $z \sim$ 1

Authors: M. Frias Castillo, M. Rybak, J. Hodge, P. van der Werf, L. J. Abbo, F. J. Ballieux, S. Ward, C. Harrison, G. Calistro Rivera, J. P. McKean, H. R. Stacey

Abstract: Unobscured quasars (QSOs) are predicted to be the final stage in the evolutionary sequence from gas-rich mergers to gas-depleted, quenched galaxies. Studies of this population, however, find a high incidence of far-infrared-luminous sources -suggesting significant dust-obscured star formation-but direct observations of the cold molecular gas fuelling this star formation are still necessary. We pre… ▽ More Unobscured quasars (QSOs) are predicted to be the final stage in the evolutionary sequence from gas-rich mergers to gas-depleted, quenched galaxies. Studies of this population, however, find a high incidence of far-infrared-luminous sources -suggesting significant dust-obscured star formation-but direct observations of the cold molecular gas fuelling this star formation are still necessary. We present a NOEMA study of CO(2-1) emission, tracing the cold molecular gas, in ten lensed z=1-1.5 unobscured QSOs. We detected CO(2-1) in seven of our targets, four of which also show continuum emission (λ_rest = 1.3mm). After subtracting the foreground galaxy contribution to the photometry, spectral energy distribution fitting yielded stellar masses of 10^9-11 M_\odot, with star formation rates of 25-160 M_\odot yr^-1 for the host galaxies. These QSOs have lower $L'_\mathrm{CO}$ than star-forming galaxies with the same L_IR, and show depletion times spanning a large range (50-900 Myr), but with a median of just 90 Myr. We find molecular gas masses in the range 2-40 x 10^9(alpha_CO/4) M_\odot, which suggest gas fractions above ~50% for most of the targets. Despite the presence of an unobscured QSO, the host galaxies are able to retain significant amounts of cold gas. However, with a median depletion time of ~90 Myr, the intense burst of star formation taking place in these targets will quickly deplete their molecular gas reservoirs in the absence of gas replenishment, resulting in a quiescent host galaxy. The non-detected QSOs are three of the four radio-loud QSOs in the sample, and their properties indicate that they are likely already transitioning into quiescence. Recent cosmological simulations tend to overestimate the depletion times expected for these z~1 QSO-host galaxies, which is likely linked to their difficulty producing starbursts across the general high-redshift galaxy population. △ Less

Submitted 12 January, 2024; v1 submitted 9 January, 2024; originally announced January 2024.

Comments: 20 pages. Accepted for publication in A&A

Quasar Feedback Survey: molecular gas affected by central outflows and by ~10 kpc radio lobes reveal dual feedback effects in `radio quiet' quasars

Authors: A. Girdhar, C. M. Harrison, V. Mainieri, R. Fernández Aranda, D. M. Alexander, F. Arrigoni Battaia, M. Bianchin, G. Calistro Rivera, C. Circosta, T. Costa, A. C. Edge, E. P. Farina, D. Kakkad, P. Kharb, S. J. Molyneux, D. Mukherjee, A. Njeri, Silpa S., G. Venturi, S. R. Ward

Abstract: We present a study of molecular gas, traced via CO (3-2) from ALMA data, of four z< 0.2, `radio quiet', type 2 quasars (log [L(bol)/(erg/s)] = 45.3 - 46.2; log [L(1.4 GHz)/(W/Hz)] = 23.7 - 24.3). Targets were selected to have extended radio lobes (>= 10 kpc), and compact, moderate-power jets (1 - 10 kpc; log [Pjet/(erg/s)]= 43.2 - 43.7). All targets show evidence of central molecular outflows, or… ▽ More We present a study of molecular gas, traced via CO (3-2) from ALMA data, of four z< 0.2, `radio quiet', type 2 quasars (log [L(bol)/(erg/s)] = 45.3 - 46.2; log [L(1.4 GHz)/(W/Hz)] = 23.7 - 24.3). Targets were selected to have extended radio lobes (>= 10 kpc), and compact, moderate-power jets (1 - 10 kpc; log [Pjet/(erg/s)]= 43.2 - 43.7). All targets show evidence of central molecular outflows, or injected turbulence, within the gas disks (traced via high-velocity wing components in CO emission-line profiles). The inferred velocities (Vout = 250 - 440 km/s) and spatial scales (0.6 - 1.6 kpc), are consistent with those of other samples of luminous low-redshift AGN. In two targets, we observe extended molecular gas structures beyond the central disks, containing 9 - 53 % of the total molecular gas mass. These structures tend to be elongated, extending from the core, and wrap-around (or along) the radio lobes. Their properties are similar to the molecular gas filaments observed around radio lobes of, mostly `radio loud', Brightest Cluster Galaxies. They have: projected distances of 5 - 13 kpc; bulk velocities of 100 - 340 km/s; velocity dispersion of 30 - 130 km/s; inferred mass outflow rates of 4 - 20 Msolar/yr; and estimated kinetic powers of log [Ekin/(erg/s)]= 40.3 - 41.7. Our observations are consistent with simulations that suggest moderate-power jets can have a direct (but modest) impact on molecular gas on small scales, through direct jet-cloud interactions. Then, on larger scales, jet-cocoons can push gas aside. Both processes could contribute to the long-term regulation of star formation. △ Less

Submitted 8 December, 2023; v1 submitted 6 November, 2023; originally announced November 2023.

Comments: Accepted for publication in MNRAS; references updated; typos corrected

The Quasar Feedback Survey: characterising CO excitation in quasar host galaxies

Authors: S. J. Molyneux, G. Calistro Rivera, C. De Breuck, C. M. Harrison, V. Mainieri, A. Lundgren, D. Kakkad, C. Circosta, A. Girdhar, T. Costa, J. R. Mullaney, P. Kharb, F. Arrigoni Battaia, E. P. Farina, D. M. Alexander, S. R. Ward, Silpa S., R. Smit

Abstract: We present a comprehensive study of the molecular gas properties of 17 Type 2 quasars at $z <$ 0.2 from the Quasar Feedback Survey (L$_{[OIII]}$ > $10^{42.1}$ $\rm ergs^{-1}$), selected by their high [OIII] luminosities and displaying a large diversity of radio jet properties, but dominated by LIRG-like galaxies. With these data, we are able to investigate the impact of AGN and AGN feedback mechan… ▽ More We present a comprehensive study of the molecular gas properties of 17 Type 2 quasars at $z <$ 0.2 from the Quasar Feedback Survey (L$_{[OIII]}$ > $10^{42.1}$ $\rm ergs^{-1}$), selected by their high [OIII] luminosities and displaying a large diversity of radio jet properties, but dominated by LIRG-like galaxies. With these data, we are able to investigate the impact of AGN and AGN feedback mechanisms on the global molecular interstellar medium. Using APEX and ALMA ACA observations, we measure the total molecular gas content using the CO(1-0) emission and homogeneously sample the CO spectral line energy distributions (SLEDs), observing CO transitions (J$_{up}$ = 1, 2, 3, 6, 7). We observe high $r_{21}$ ratios (r$_{21}$ = L'$_{CO(2-1)}$/L'$_{CO(1-0)}$) with a median $r_{21}$ = 1.06, similar to local (U)LIRGs (with $r_{21}$ $\sim$ 1) and higher than normal star-forming galaxies (with r$_{21}$ $\sim$ 0.65). Despite the high $r_{21}$ values, for the 7 targets with the required data we find low excitation in CO(6-5) & CO(7-6) ($r_{61}$ and $r_{62}$ < 0.6 in all but one target), unlike high redshift quasars in the literature, which are far more luminous and show higher line ratios. The ionised gas traced by [OIII] exhibit systematically higher velocities than the molecular gas traced by CO. We conclude that any effects of quasar feedback (e.g. via outflows and radio jets) do not have a significant instantaneous impact on the global molecular gas content and excitation and we suggest that it only occurs on more localised scales. △ Less

Submitted 27 October, 2023; v1 submitted 16 October, 2023; originally announced October 2023.

Comments: 32 pages (20 in the main body of the paper and 12 in the appendix), 28 figures (10 in main body of paper and 18 in appendix) Accepted for publication in MNRAS. Data available at https://doi.org/10.25405/data.ncl.24312502

Bird-Snack: Bayesian Inference of dust law $R_V$ Distributions using SN Ia Apparent Colours at peaK

Authors: Sam M. Ward, Suhail Dhawan, Kaisey S. Mandel, Matthew Grayling, Stephen Thorp

Abstract: To reduce systematic uncertainties in Type Ia supernova (SN Ia) cosmology, the host galaxy dust law shape parameter, $R_V$, must be accurately constrained. We thus develop a computationally-inexpensive pipeline, Bird-Snack, to rapidly infer dust population distributions from optical-near infrared SN colours at peak brightness, and determine which analysis choices significantly impact the populatio… ▽ More To reduce systematic uncertainties in Type Ia supernova (SN Ia) cosmology, the host galaxy dust law shape parameter, $R_V$, must be accurately constrained. We thus develop a computationally-inexpensive pipeline, Bird-Snack, to rapidly infer dust population distributions from optical-near infrared SN colours at peak brightness, and determine which analysis choices significantly impact the population mean $R_V$ inference, $μ_{R_V}$. Our pipeline uses a 2D Gaussian process to measure peak $BVriJH$ apparent magnitudes from SN light curves, and a hierarchical Bayesian model to simultaneously constrain population distributions of intrinsic and dust components. Fitting a low-to-moderate-reddening sample of 65 low-redshift SNe yields $μ_{R_V}=2.61^{+0.38}_{-0.35}$, with $68\%(95\%)$ posterior upper bounds on the population dispersion, $σ_{R_V}<0.92(1.96)$. This result is robust to various analysis choices, including: the model for intrinsic colour variations, fitting the shape hyperparameter of a gamma dust extinction distribution, and cutting the sample based on the availability of data near peak. However, these choices may be important if statistical uncertainties are reduced. With larger near-future optical and near-infrared SN samples, Bird-Snack can be used to better constrain dust distributions, and investigate potential correlations with host galaxy properties. Bird-Snack is publicly available; the modular infrastructure facilitates rapid exploration of custom analysis choices, and quick fits to simulated datasets, for better interpretation of real-data inferences. △ Less

Submitted 11 October, 2023; originally announced October 2023.

Comments: Accepted for publication in MNRAS. 21 pages, 12 figures, 14 tables

Keck Infrared Transient Survey I: Survey Description and Data Release 1

Authors: S. Tinyanont, R. J. Foley, K. Taggart, K. W. Davis, N. LeBaron, J. E. Andrews, M. J. Bustamante-Rosell, Y. Camacho-Neves, R. Chornock, D. A. Coulter, L. Galbany, S. W. Jha, C. D. Kilpatrick, L. A. Kwok, C. Larison, J. R. Pierel, M. R. Siebert, G. Aldering, K. Auchettl, J. S. Bloom, S. Dhawan, A. V. Filippenko, K. D. French, A. Gagliano, M. Grayling , et al. (13 additional authors not shown)

Abstract: We present the Keck Infrared Transient Survey (KITS), a NASA Key Strategic Mission Support program to obtain near-infrared (NIR) spectra of astrophysical transients of all types, and its first data release, consisting of 105 NIR spectra of 50 transients. Such a data set is essential as we enter a new era of IR astronomy with the James Webb Space Telescope (JWST) and the upcoming Nancy Grace Roman… ▽ More We present the Keck Infrared Transient Survey (KITS), a NASA Key Strategic Mission Support program to obtain near-infrared (NIR) spectra of astrophysical transients of all types, and its first data release, consisting of 105 NIR spectra of 50 transients. Such a data set is essential as we enter a new era of IR astronomy with the James Webb Space Telescope (JWST) and the upcoming Nancy Grace Roman Space Telescope (Roman). NIR spectral templates will be essential to search JWST images for stellar explosions of the first stars and to plan an effective Roma} SN Ia cosmology survey, both key science objectives for mission success. Between 2022 February and 2023 July, we systematically obtained 274 NIR spectra of 146 astronomical transients, representing a significant increase in the number of available NIR spectra in the literature. The first data release includes data from the 2022A semester. We systematically observed three samples: a flux-limited sample that includes all transients $<$17 mag in a red optical band (usually ZTF r or ATLAS o bands); a volume-limited sample including all transients within redshift $z < 0.01$ ($D \approx 50$ Mpc); and an SN Ia sample targeting objects at phases and light-curve parameters that had scant existing NIR data in the literature. The flux-limited sample is 39% complete (60% excluding SNe Ia), while the volume-limited sample is 54% complete and is 79% complete to $z = 0.005$. All completeness numbers will rise with the inclusion of data from other telescopes in future data releases. Transient classes observed include common Type Ia and core-collapse supernovae, tidal disruption events (TDEs), luminous red novae, and the newly categorized hydrogen-free/helium-poor interacting Type Icn supernovae. We describe our observing procedures and data reduction using Pypeit, which requires minimal human interaction to ensure reproducibility. △ Less

Submitted 13 September, 2023; originally announced September 2023.

Establishing the impact of luminous AGN with multi-wavelength observations and simulations

Authors: C. M. Harrison, A. Girdhar, S. R. Ward

Abstract: Cosmological simulations fail to reproduce realistic galaxy populations without energy injection from active galactic nuclei (AGN) into the interstellar medium (ISM) and circumgalactic medium (CGM); a process called `AGN feedback'. Consequently, observational work searches for evidence that luminous AGN impact their host galaxies. Here, we review some of this work. Multi-phase AGN outflows are com… ▽ More Cosmological simulations fail to reproduce realistic galaxy populations without energy injection from active galactic nuclei (AGN) into the interstellar medium (ISM) and circumgalactic medium (CGM); a process called `AGN feedback'. Consequently, observational work searches for evidence that luminous AGN impact their host galaxies. Here, we review some of this work. Multi-phase AGN outflows are common, some with potential for significant impact. Additionally, multiple feedback channels can be observed simultaneously; e.g., radio jets from `radio quiet' quasars can inject turbulence on ISM scales, and displace CGM-scale molecular gas. However, caution must be taken comparing outflows to simulations (e.g., kinetic coupling efficiencies) to infer feedback potential, due to a lack of comparable predictions. Furthermore, some work claims limited evidence for feedback because AGN live in gas-rich, star-forming galaxies. However, simulations do not predict instantaneous, global impact on molecular gas or star formation. The impact is expected to be cumulative, over multiple episodes. △ Less

Submitted 7 July, 2023; originally announced July 2023.

Comments: Accepted for publication in IAU Symposium 378 Conference Proceedings "Black Hole Winds at all Scales"

Quantum Decay of an Optical Soliton

Authors: Stuart Ward, Rouzbeh Allahverdi, Arash Mafi

Abstract: Optical solitons are known to be classically stable objects which are robust to perturbations. In this work, we show that due to quantum mechanical effects, an optical soliton that is initially in a classical soliton coherent state will shed photons into the continuum and hence decay. The standard formulation of the quantized soliton uses the linearized version of the quantum nonlinear Schrodinger… ▽ More Optical solitons are known to be classically stable objects which are robust to perturbations. In this work, we show that due to quantum mechanical effects, an optical soliton that is initially in a classical soliton coherent state will shed photons into the continuum and hence decay. The standard formulation of the quantized soliton uses the linearized version of the quantum nonlinear Schrodinger equation in the background of the classical soliton, and the quantized soliton remains stable in this approximation. We show that if higher-order interaction terms are taken into account, the soliton is no longer stable, and its photon number decreases quadratically as a function of the number of soliton cycles. We compute the power spectrum for the continuum radiation and find a narrow band that is localized about the initial soliton momentum with a cut-off that is inversely proportional to the initial soliton width. △ Less

Submitted 8 February, 2023; v1 submitted 1 February, 2023; originally announced February 2023.

The DEHVILS Survey Overview and Initial Data Release: High-Quality Near-Infrared Type Ia Supernova Light Curves at Low Redshift

Authors: Erik R. Peterson, David O. Jones, Daniel Scolnic, Bruno O. Sánchez, Aaron Do, Adam G. Riess, Sam M. Ward, Arianna Dwomoh, Thomas de Jaeger, Saurabh W. Jha, Kaisey S. Mandel, Justin D. R. Pierel, Brodie Popovic, Benjamin M. Rose, David Rubin, Benjamin J. Shappee, Stephen Thorp, John L. Tonry, R. Brent Tully, Maria Vincenzi

Abstract: While the sample of optical Type Ia Supernova (SN Ia) light curves (LCs) usable for cosmological parameter measurements surpasses 2000, the sample of published, cosmologically viable near-infrared (NIR) SN Ia LCs, which have been shown to be good "standard candles," is still $\lesssim$ 200. Here, we present high-quality NIR LCs for 83 SNe Ia ranging from $0.002 < z < 0.09$ as a part of the Dark En… ▽ More While the sample of optical Type Ia Supernova (SN Ia) light curves (LCs) usable for cosmological parameter measurements surpasses 2000, the sample of published, cosmologically viable near-infrared (NIR) SN Ia LCs, which have been shown to be good "standard candles," is still $\lesssim$ 200. Here, we present high-quality NIR LCs for 83 SNe Ia ranging from $0.002 < z < 0.09$ as a part of the Dark Energy, H$_0$, and peculiar Velocities using Infrared Light from Supernovae (DEHVILS) survey. Observations are taken using UKIRT's WFCAM, where the median depth of the images is 20.7, 20.1, and 19.3 mag (Vega) for $Y$, $J$, and $H$-bands, respectively. The median number of epochs per SN Ia is 18 for all three bands ($YJH$) combined and 6 for each band individually. We fit 47 SN Ia LCs that pass strict quality cuts using three LC models, SALT3, SNooPy, and BayeSN and find scatter on the Hubble diagram to be comparable to or better than scatter from optical-only fits in the literature. Fitting NIR-only LCs, we obtain standard deviations ranging from 0.128-0.135 mag. Additionally, we present a refined calibration method for transforming 2MASS magnitudes to WFCAM magnitudes using HST CALSPEC stars that results in a 0.03 mag shift in the WFCAM $Y$-band magnitudes. △ Less

Submitted 10 April, 2023; v1 submitted 27 January, 2023; originally announced January 2023.

Comments: 24 pages, 9 figures. Accepted by MNRAS

A BayeSN Distance Ladder: $H_0$ from a consistent modelling of Type Ia supernovae from the optical to the near infrared

Authors: Suhail Dhawan, Stephen Thorp, Kaisey S. Mandel, Sam M. Ward, Gautham Narayan, Saurabh W. Jha, Thaisen Chant

Abstract: The local distance ladder estimate of the Hubble constant ($H_0$) is important in cosmology, given the recent tension with the early universe inference. We estimate $H_0$ from the Type Ia supernova (SN~Ia) distance ladder, inferring SN~Ia distances with the hierarchical Bayesian SED model, BayeSN. This method has a notable advantage of being able to continuously model the optical and near-infrared… ▽ More The local distance ladder estimate of the Hubble constant ($H_0$) is important in cosmology, given the recent tension with the early universe inference. We estimate $H_0$ from the Type Ia supernova (SN~Ia) distance ladder, inferring SN~Ia distances with the hierarchical Bayesian SED model, BayeSN. This method has a notable advantage of being able to continuously model the optical and near-infrared (NIR) SN~Ia light curves simultaneously. We use two independent distance indicators, Cepheids or the tip of the red giant branch (TRGB), to calibrate a Hubble-flow sample of 67 SNe~Ia with optical and NIR data. We estimate $H_0 = 74.82 \pm 0.97$ (stat) $\pm\, 0.84$ (sys) km\,s$^{-1}$\,Mpc$^{-1}$ when using the calibration with Cepheid distances to 37 host galaxies of 41 SNe~Ia, and $70.92 \pm 1.14$ (stat) $\pm\,1.49$ (sys) km\,s$^{-1}$\,Mpc$^{-1}$ when using the calibration with TRGB distances to 15 host galaxies of 18 SNe~Ia. For both methods, we find a low intrinsic scatter $σ_{\rm int} \lesssim 0.1$ mag. We test various selection criteria and do not find significant shifts in the estimate of $H_0$. Simultaneous modelling of the optical and NIR yields up to $\sim$15\% reduction in $H_0$ uncertainty compared to the equivalent optical-only cases. With improvements expected in other rungs of the distance ladder, leveraging joint optical-NIR SN~Ia data can be critical to reducing the $H_0$ error budget. △ Less

Submitted 2 August, 2023; v1 submitted 14 November, 2022; originally announced November 2022.

Comments: 10 pages, 8 figures, published in MNRAS

The Young Supernova Experiment Data Release 1 (YSE DR1): Light Curves and Photometric Classification of 1975 Supernovae

Authors: P. D. Aleo, K. Malanchev, S. Sharief, D. O. Jones, G. Narayan, R. J. Foley, V. A. Villar, C. R. Angus, V. F. Baldassare, M. J. Bustamante-Rosell, D. Chatterjee, C. Cold, D. A. Coulter, K. W. Davis, S. Dhawan, M. R. Drout, A. Engel, K. D. French, A. Gagliano, C. Gall, J. Hjorth, M. E. Huber, W. V. Jacobson-Galán, C. D. Kilpatrick, D. Langeroodi , et al. (58 additional authors not shown)

Abstract: We present the Young Supernova Experiment Data Release 1 (YSE DR1), comprised of processed multi-color Pan-STARRS1 (PS1) griz and Zwicky Transient Facility (ZTF) gr photometry of 1975 transients with host-galaxy associations, redshifts, spectroscopic/photometric classifications, and additional data products from 2019 November 24 to 2021 December 20. YSE DR1 spans discoveries and observations from… ▽ More We present the Young Supernova Experiment Data Release 1 (YSE DR1), comprised of processed multi-color Pan-STARRS1 (PS1) griz and Zwicky Transient Facility (ZTF) gr photometry of 1975 transients with host-galaxy associations, redshifts, spectroscopic/photometric classifications, and additional data products from 2019 November 24 to 2021 December 20. YSE DR1 spans discoveries and observations from young and fast-rising supernovae (SNe) to transients that persist for over a year, with a redshift distribution reaching z~0.5. We present relative SN rates from YSE's magnitude- and volume-limited surveys, which are consistent with previously published values within estimated uncertainties for untargeted surveys. We combine YSE and ZTF data, and create multi-survey SN simulations to train the ParSNIP and SuperRAENN photometric classification algorithms; when validating our ParSNIP classifier on 472 spectroscopically classified YSE DR1 SNe, we achieve 82% accuracy across three SN classes (SNe Ia, II, Ib/Ic) and 90% accuracy across two SN classes (SNe Ia, core-collapse SNe). Our classifier performs particularly well on SNe Ia, with high (>90%) individual completeness and purity, which will help build an anchor photometric SNe Ia sample for cosmology. We then use our photometric classifier to characterize our photometric sample of 1483 SNe, labeling 1048 (~71%) SNe Ia, 339 (~23%) SNe II, and 96 (~6%) SNe Ib/Ic. YSE DR1 provides a training ground for building discovery, anomaly detection, and classification algorithms, performing cosmological analyses, understanding the nature of red and rare transients, exploring tidal disruption events and nuclear variability, and preparing for the forthcoming Vera C. Rubin Observatory Legacy Survey of Space and Time. △ Less

Submitted 21 February, 2023; v1 submitted 14 November, 2022; originally announced November 2022.

Comments: Accepted to ApJS; 64 pages; 35 figures; 10 tables

Relative intrinsic scatter in hierarchical Type Ia supernova siblings analyses: Application to SNe 2021hpr, 1997bq & 2008fv in NGC 3147

Authors: Sam M. Ward, Stephen Thorp, Kaisey S. Mandel, Suhail Dhawan, David O. Jones, Kirsty Taggart, Ryan J. Foley, Gautham Narayan, Kenneth C. Chambers, David A. Coulter, Kyle W. Davis, Thomas de Boer, Kaylee de Soto, Nicholas Earl, Alex Gagliano, Hua Gao, Jens Hjorth, Mark E. Huber, Luca Izzo, Danial Langeroodi, Eugene A. Magnier, Peter McGill, Armin Rest, César Rojas-Bravo, Radosław Wojtak

Abstract: We present Young Supernova Experiment $grizy$ photometry of SN 2021hpr, the third Type Ia supernova sibling to explode in the Cepheid calibrator galaxy, NGC 3147. Siblings are useful for improving SN-host distance estimates, and investigating the contributions towards the SN Ia intrinsic scatter (post-standardisation residual scatter in distance estimates). We thus develop a principled Bayesian fr… ▽ More We present Young Supernova Experiment $grizy$ photometry of SN 2021hpr, the third Type Ia supernova sibling to explode in the Cepheid calibrator galaxy, NGC 3147. Siblings are useful for improving SN-host distance estimates, and investigating the contributions towards the SN Ia intrinsic scatter (post-standardisation residual scatter in distance estimates). We thus develop a principled Bayesian framework for analyzing SN Ia siblings. At its core is the cosmology-independent relative intrinsic scatter parameter, $σ_{Rel}$: the dispersion of siblings distance estimates relative to one another within a galaxy. It quantifies the contribution towards the total intrinsic scatter, $σ_0$, from within-galaxy variations about the siblings' common properties. It also affects the combined-distance uncertainty. We present analytic formulae for computing a $σ_{Rel}$-posterior from individual siblings distances (estimated using any SN-model). Applying a newly trained BayeSN model, we fit the light curves of each sibling in NGC 3147 individually, to yield consistent distance estimates. However, the wide $σ_{Rel}$-posterior means $σ_{Rel}\approxσ_0$ is not ruled out. We thus combine the distances by marginalizing over $σ_{Rel}$ with an informative prior: $σ_{Rel}\sim U(0,σ_0)$. Simultaneously fitting the trio's light curves improves constraints on distance, and each sibling's individual dust parameters, compared to individual fits. Higher correlation also tightens dust parameter constraints. Therefore, $σ_{Rel}$-marginalization yields robust estimates of siblings distances for cosmology, and dust parameters for siblings-host correlation studies. Incorporating NGC 3147's Cepheid-distance yields $H_0=78.4\pm 6.5\,$km/s/Mpc. Our work motivates analyses of homogeneous siblings samples, to constrain $σ_{Rel}$, and its SN-model dependence. △ Less

Submitted 1 September, 2023; v1 submitted 21 September, 2022; originally announced September 2022.

Comments: Accepted for publication in ApJ. 22 pages, 11 figures, 8 tables. Single-galaxy siblings analysis

Cosmological simulations predict that AGN preferentially live in gas-rich, star-forming galaxies despite effective feedback

Authors: Samuel Ruthven Ward, Chris Harrison, Tiago Costa, Vincenzo Mainieri

Abstract: Negative feedback from active galactic nuclei (AGN) is the leading mechanism for the quenching of massive galaxies in the vast majority of modern galaxy evolution models. However, direct observational evidence that AGN feedback causes quenching on a population scale is lacking. Studies have shown that luminous AGN are preferentially located in gas-rich and star-forming galaxies, an observation tha… ▽ More Negative feedback from active galactic nuclei (AGN) is the leading mechanism for the quenching of massive galaxies in the vast majority of modern galaxy evolution models. However, direct observational evidence that AGN feedback causes quenching on a population scale is lacking. Studies have shown that luminous AGN are preferentially located in gas-rich and star-forming galaxies, an observation that has sometimes been suggested to be in tension with a negative AGN feedback picture. We investigate three of the current cosmological simulations (IllustrisTNG, EAGLE and SIMBA) along with post-processed models for molecular hydrogen gas masses and perform similar tests to those used by observers. We find that the simulations predict: (i) no strong negative trends between AGN luminosity and molecular gas fraction or sSFR; (ii) both high-luminosity ($L_{bol}>10^{44}$ erg/s) and high-Eddington ratio (>1%) AGN are preferentially located in galaxies with high molecular gas fractions and sSFR; and (iii) that the gas-depleted and quenched fractions of AGN host galaxies are lower than a control sample of non-active galaxies. These three findings are in qualitative agreement with observational samples at $z=0$ and $z=2$ and show that such results are not in tension with the presence of strong AGN feedback, which all simulations we employ require to produce realistic massive galaxies. However, we also find quantifiable differences between predictions from the simulations, which could allow us to observationally test the different subgrid feedback models. △ Less

Submitted 28 April, 2022; originally announced April 2022.

Comments: Accepted for publication in MNRAS. 22 pages, 12 figures, 2 tables

Capture Agent Free Biosensing using Porous Silicon Arrays and Machine Learning

Authors: Simon J. Ward, Tengfei Cao, Xiang Zhou, Catie Chang, Sharon M. Weiss

Abstract: Biosensors are an essential tool for medical diagnostics, environmental monitoring and food safety. Typically, biosensors are designed to detect specific analytes through functionalization with the appropriate capture agents. However, the use of capture agents limits the number of analytes that can be simultaneously detected and reduces the robustness of the biosensor. In this work, we report a ve… ▽ More Biosensors are an essential tool for medical diagnostics, environmental monitoring and food safety. Typically, biosensors are designed to detect specific analytes through functionalization with the appropriate capture agents. However, the use of capture agents limits the number of analytes that can be simultaneously detected and reduces the robustness of the biosensor. In this work, we report a versatile, capture agent free biosensor platform based on an array of porous silicon (PSi) thin films, which has the potential to robustly detect a wide variety of analytes based on their physical and chemical properties in the nanoscale porous media. The ability of this system to reproducibly classify, quantify, and discriminate three proteins is demonstrated to concentrations down to at least 0.02g/L (between 300nM and 450nM) by utilizing PSi array elements with a unique combination of pore size and buffer pH, employing linear discriminant analysis for dimensionality reduction, and using support vector machines as a classifier. This approach represents a significant step towards a low cost, simple and robust biosensor platform that is able to detect a vast range of biomolecules. △ Less

Submitted 22 January, 2022; originally announced January 2022.

Comments: 15 pages, 3 figures, 2 tables

Journal ref: Biosensors 13 (2023) 1-12

Quasar Feedback Survey: Multi-phase outflows, turbulence and evidence for feedback caused by low power radio jets inclined into the galaxy disk

Authors: A. Girdhar, C. M. Harrison, V. Mainieri, A. Bittner, T. Costa, P. Kharb, D. Mukherjee, F. Arrigoni Battaia, D. M. Alexander, G. Calistro Rivera, C. Circosta, C. De Breuck, A. C. Edge, E. P. Farina, D. Kakkad, G. B. Lansbury, S. J. Molyneux, J. R. Mullaney, Silpa S., A. P. Thomson, S. R. Ward

Abstract: We present a study of a luminous, z=0.15, type-2 quasar (log [L([OIII])/(erg/s)]=42.8) from the Quasar Feedback Survey. It is classified as 'radio-quiet' (log [L(1.4 GHz)/(W/Hz)]=23.8); however, radio imaging reveals ~1 kpc low-power jets (log [Pjet/(erg/s)]=44) inclined into the plane of the galaxy disk. We combine MUSE and ALMA observations to map stellar kinematics and ionised and molecular gas… ▽ More We present a study of a luminous, z=0.15, type-2 quasar (log [L([OIII])/(erg/s)]=42.8) from the Quasar Feedback Survey. It is classified as 'radio-quiet' (log [L(1.4 GHz)/(W/Hz)]=23.8); however, radio imaging reveals ~1 kpc low-power jets (log [Pjet/(erg/s)]=44) inclined into the plane of the galaxy disk. We combine MUSE and ALMA observations to map stellar kinematics and ionised and molecular gas properties. The jets are seen to drive galaxy-wide bi-conical turbulent outflows, reaching W80 = 1000-1300 km/s, in the ionised phase (traced via optical emission-lines), which also have increased electron densities compared to the quiescent gas. The turbulent gas is driven perpendicular to the jet axis and is escaping along the galaxy minor axis, reaching 7.5 kpc on both sides. Traced via CO(3-2) emission, the turbulent material in molecular gas phase is one-third as spatially extended and has 3 times lower velocity-dispersion as compared to ionised gas. The jets are seen to be strongly interacting with the interstellar medium (ISM) through enhanced ionised emission and disturbed/depleted molecular gas at the jet termini. We see further evidence for jet-induced feedback through significantly higher stellar velocity-dispersion aligned, and co-spatial with, the jet axis (<5 deg). We discuss possible negative and positive feedback scenarios arising due to the interaction of the low-power jets with the ISM in the context of recent jet-ISM interaction simulations, which qualitatively agree with our observations. We discuss how jet-induced feedback could be an important feedback mechanism even in bolometrically luminous 'radio-quiet' quasars. △ Less

Submitted 11 January, 2022; v1 submitted 6 January, 2022; originally announced January 2022.

Comments: Accepted for publication in MNRAS. Main manuscript has 21 pages with 8 figures. Supplementary material is available for download under "Ancillary files" or by downloading the source file listed under "Other formats"

Randomness and Frustration in a S = 1/2 Square-Lattice Heisenberg Antiferromagnet

Authors: Ellen Fogh, Otto Mustonen, Peter Babkevich, Vamshi M. Katukuri, Helen C. Walker, Lucile Mangin-Thro, Maarit Karppinen, Simon Ward, Bruce Normand, Henrik M. Rønnow

Abstract: We explore the interplay between randomness and magnetic frustration in the series of $S = 1/2$ Heisenberg square-lattice compounds Sr$_2$CuTe$_{1-x}$W$_x$O$_6$. Substituting W for Te alters the magnetic interactions dramatically, from strongly nearest-neighbor to next-nearest-neighbor antiferromagnetic coupling. We perform neutron scattering measurements to probe the magnetic ground state and exc… ▽ More We explore the interplay between randomness and magnetic frustration in the series of $S = 1/2$ Heisenberg square-lattice compounds Sr$_2$CuTe$_{1-x}$W$_x$O$_6$. Substituting W for Te alters the magnetic interactions dramatically, from strongly nearest-neighbor to next-nearest-neighbor antiferromagnetic coupling. We perform neutron scattering measurements to probe the magnetic ground state and excitations over a range of $x$. We propose a bond-disorder model that reproduces ground states with only short-ranged spin correlations in the mixed compounds. The calculated neutron diffraction patterns and powder spectra agree well with the measured data and allow detailed predictions for future measurements. We conclude that quenched randomness plays the major role in defining the physics of Sr$_2$CuTe$_{1-x}$W$_x$O$_6$ with frustration being less significant. △ Less

Submitted 20 May, 2022; v1 submitted 6 December, 2021; originally announced December 2021.

Journal ref: Physical Review B 105, 184410 (2022)

Borrowing from Similar Code: A Deep Learning NLP-Based Approach for Log Statement Automation

Authors: Sina Gholamian, Paul A. S. Ward

Abstract: Software developers embed logging statements inside the source code as an imperative duty in modern software development as log files are necessary for tracking down runtime system issues and troubleshooting system management tasks. However, the current logging process is mostly manual, and thus, proper placement and content of logging statements remain as challenges. To overcome these challenges,… ▽ More Software developers embed logging statements inside the source code as an imperative duty in modern software development as log files are necessary for tracking down runtime system issues and troubleshooting system management tasks. However, the current logging process is mostly manual, and thus, proper placement and content of logging statements remain as challenges. To overcome these challenges, methods that aim to automate log placement and predict its content, i.e., 'where and what to log', are of high interest. Thus, we focus on predicting the location (i.e., where) and description (i.e., what) for log statements by utilizing source code clones and natural language processing (NLP), as these approaches provide additional context and advantage for log prediction. Specifically, we guide our research with three research questions (RQs): (RQ1) how similar code snippets, i.e., code clones, can be leveraged for log statements prediction? (RQ2) how the approach can be extended to automate log statements' descriptions? and (RQ3) how effective the proposed methods are for log location and description prediction? To pursue our RQs, we perform an experimental study on seven open-source Java projects. We introduce an updated and improved log-aware code-clone detection method to predict the location of logging statements (RQ1). Then, we incorporate natural language processing (NLP) and deep learning methods to automate the log statements' description prediction (RQ2). Our analysis shows that our hybrid NLP and code-clone detection approach (NLP CC'd) outperforms conventional clone detectors in finding log statement locations on average by 15.60% and achieves 40.86% higher performance on BLEU and ROUGE scores for predicting the description of logging statements when compared to prior research (RQ3). △ Less

Submitted 2 December, 2021; originally announced December 2021.

Comments: 27 pages

A Reference Survey for Supernova Cosmology with the Nancy Grace Roman Space Telescope

Authors: B. M. Rose, C. Baltay, R. Hounsell, P. Macias, D. Rubin, D. Scolnic, G. Aldering, R. Bohlin, M. Dai, S. E. Deustua, R. J. Foley, A. Fruchter, L. Galbany, S. W. Jha, D. O. Jones, B. A. Joshi, P. L. Kelly, R. Kessler, R. P. Kirshner, K. S. Mandel, S. Perlmutter, J. Pierel, H. Qu, D. Rabinowitz, A. Rest , et al. (11 additional authors not shown)

Abstract: This note presents an initial survey design for the Nancy Grace Roman High-latitude Time Domain Survey. This is not meant to be a final or exhaustive list of all the survey strategy choices, but instead presents a viable path towards achieving the desired precision and accuracy of dark energy measurements using Type Ia supernovae (SNe Ia). We describe a survey strategy that use six filters (RZYJH… ▽ More This note presents an initial survey design for the Nancy Grace Roman High-latitude Time Domain Survey. This is not meant to be a final or exhaustive list of all the survey strategy choices, but instead presents a viable path towards achieving the desired precision and accuracy of dark energy measurements using Type Ia supernovae (SNe Ia). We describe a survey strategy that use six filters (RZYJH and F) and the prism on the Roman Wide Field Instrument. This survey has two tiers, one "wide" which targets SNe Ia at redshifts up to 1 and one "deep" targeting redshifts up to 1.7; for each, four filters are used (with Y and J used in both tiers). We propose one field each in the north and south continuous viewing zones, and expect to obtain high-quality distances of $\sim$12,000 SNe Ia with $\sim$5,000 at z > 1. We propose a wide-tier area of $\sim$19 deg$^2$ and a deep tier of $\sim$5 deg$^2$. Exposure times range from 100 s to 900 s for imaging and 900 s to 3600 s for the prism. These exposure times would reach $\sim$25.5 mag and $\sim$26.5 mag for the wide and deep tiers respectively, with deep co-add stacks reaching $\sim$28 mag and $\sim$29 mag. The total survey spans two years, with a total allocation time of six months, and a cadence of $\sim$5 days. △ Less

Submitted 4 November, 2021; originally announced November 2021.

Comments: A report to NASA from the Roman Supernova Science Investigation Teams

A Comprehensive Survey of Logging in Software: From Logging Statements Automation to Log Mining and Analysis

Authors: Sina Gholamian, Paul A. S. Ward

Abstract: Logs are widely used to record runtime information of software systems, such as the timestamp and the importance of an event, the unique ID of the source of the log, and a part of the state of a task's execution. The rich information of logs enables system developers (and operators) to monitor the runtime behaviors of their systems and further track down system problems and perform analysis on log… ▽ More Logs are widely used to record runtime information of software systems, such as the timestamp and the importance of an event, the unique ID of the source of the log, and a part of the state of a task's execution. The rich information of logs enables system developers (and operators) to monitor the runtime behaviors of their systems and further track down system problems and perform analysis on log data in production settings. However, the prior research on utilizing logs is scattered and that limits the ability of new researchers in this field to quickly get to the speed and hampers currently active researchers to advance this field further. Therefore, this paper surveys and provides a systematic literature review and mapping of the contemporary logging practices and log statements' mining and monitoring techniques and their applications such as in system failure detection and diagnosis. We study a large number of conference and journal papers that appeared on top-level peer-reviewed venues. Additionally, we draw high-level trends of ongoing research and categorize publications into subdivisions. In the end, and based on our holistic observations during this survey, we provide a set of challenges and opportunities that will lead the researchers in academia and industry in moving the field forward. △ Less

Submitted 1 January, 2022; v1 submitted 24 October, 2021; originally announced October 2021.

Comments: Under second revision at Transactions on Software Engineering -- Second Revision

Transaction Placement in Sharded Blockchains

Authors: Liuyang Ren, Paul A. S. Ward

Abstract: While many researchers adopt a sharding approach to design scaling blockchains, few works have studied the transaction placement problem incurred by sharding protocols. The widely-used hashing placement algorithm renders an overwhelming portion of transactions as cross-shard. In this paper, we analyze the high cost of cross-shard transactions and reveal that most Bitcoin transactions have simple d… ▽ More While many researchers adopt a sharding approach to design scaling blockchains, few works have studied the transaction placement problem incurred by sharding protocols. The widely-used hashing placement algorithm renders an overwhelming portion of transactions as cross-shard. In this paper, we analyze the high cost of cross-shard transactions and reveal that most Bitcoin transactions have simple dependencies and can become single-shard under a placement algorithm taking transaction dependencies into account. In addition, we perform a case study of OptChain, which is the state-of-the-art transaction placement algorithm for sharded blockchains, and find a shortcoming of it. A simple fix is proposed, and our evaluation results demonstrate that the proposed fix effectively helps OptChain overcome the shortcoming and significantly improve the system performance under a special workload. The authors of OptChain made some revisions to the algorithm description after noticing our work. Their updated algorithm does not exhibit the shortcoming under the workloads employed by this paper. △ Less

Submitted 9 June, 2022; v1 submitted 15 September, 2021; originally announced September 2021.

What Distributed Systems Say: A Study of Seven Spark Application Logs

Authors: Sina Gholamian, Paul A. S. Ward

Abstract: Execution logs are a crucial medium as they record runtime information of software systems. Although extensive logs are helpful to provide valuable details to identify the root cause in postmortem analysis in case of a failure, this may also incur performance overhead and storage cost. Therefore, in this research, we present the result of our experimental study on seven Spark benchmarks to illustr… ▽ More Execution logs are a crucial medium as they record runtime information of software systems. Although extensive logs are helpful to provide valuable details to identify the root cause in postmortem analysis in case of a failure, this may also incur performance overhead and storage cost. Therefore, in this research, we present the result of our experimental study on seven Spark benchmarks to illustrate the impact of different logging verbosity levels on the execution time and storage cost of distributed software systems. We also evaluate the log effectiveness and the information gain values, and study the changes in performance and the generated logs for each benchmark with various types of distributed system failures. Our research draws insightful findings for developers and practitioners on how to set up and utilize their distributed systems to benefit from the execution logs. △ Less

Submitted 18 August, 2021; originally announced August 2021.

Comments: 2021 The 40th International Symposium on Reliable Distributed Systems (SRDS)

A Supersymmetric Approach to the Problem of Micro-bending Attenuation in Optical Waveguides

Authors: Stuart Ward, Rouzbeh Allahverdi, Arash Mafi

Abstract: Micro-bending is a well-known source of loss in optical waveguides. By treating the micro-bending as a stochastic process, the problem of loss mitigation can be modeled in terms of a Fokker-Planck equation. Given an initial refractive index profile, and taking micro-bending into account, we develop a formalism to derive a new refractive index profile which potentially results in less loss. Our for… ▽ More Micro-bending is a well-known source of loss in optical waveguides. By treating the micro-bending as a stochastic process, the problem of loss mitigation can be modeled in terms of a Fokker-Planck equation. Given an initial refractive index profile, and taking micro-bending into account, we develop a formalism to derive a new refractive index profile which potentially results in less loss. Our formalism is based on applying the techniques of Supersymmetric Quantum Mechanics to a Fokker-Planck equation that is associated with a particular refractive index profile. We derive a non-linear differential equation, whose solutions determine whether an index profile can undergo a supersymmetric transformation that results in less loss. As an explicit example, we consider a monomial index profile. We show that there exists a range of values for the monomial exponent which results in the new index profile having less loss. △ Less

Submitted 3 August, 2021; originally announced August 2021.

Signal Processing Techniques to Reduce the Limit of Detection for Thin Film Biosensors

Authors: Simon J. Ward, Rabeb Layouni, Sofia Arshavsky-Graham, Ester Segal, Sharon M. Weiss

Abstract: The ultimate detection limit of optical biosensors is often limited by various noise sources, including those introduced by the optical measurement setup. While sophisticated modifications to instrumentation may reduce noise, a simpler approach that can benefit all sensor platforms is the application of signal processing to minimize the deleterious effects of noise. In this work, we show that appl… ▽ More The ultimate detection limit of optical biosensors is often limited by various noise sources, including those introduced by the optical measurement setup. While sophisticated modifications to instrumentation may reduce noise, a simpler approach that can benefit all sensor platforms is the application of signal processing to minimize the deleterious effects of noise. In this work, we show that applying complex Morlet wavelet convolution to Fabry-Pérot interference fringes characteristic of thin film reflectometric biosensors effectively filters out white noise and low frequency reflectance variations. Subsequent calculation of an average difference in phase between the filtered analyte and reference signals enables a significant reduction in the limit of detection (LOD) enabling closer competition with current state-of-the-art techniques. This method is applied on experimental data sets of thin film porous silicon sensors (PSi) in buffered solution and complex media obtained from two different laboratories. The demonstrated improvement in LOD achieved using wavelet convolution and average phase difference paves the way for PSi optical biosensors to operate with clinically relevant detection limits for medical diagnostics, environmental monitoring, and food safety. △ Less

Submitted 12 March, 2021; originally announced March 2021.

Comments: 11 pages, 3 Figures, 2 Tables

Journal ref: ACS Sensors 6, 2967-2978 (2021)

The Quasar Feedback Survey: Discovering hidden Radio-AGN and their connection to the host galaxy ionised gas

Authors: M. E. Jarvis, C. M. Harrison, V. Mainieri, D. M. Alexander, F. Arrigoni Battaia, G. Calistro Rivera, C. Circosta, T. Costa, C. De Breuck, A. C. Edge, A. Girdhar, D. Kakkad, P. Kharb, G. B. Lansbury, S. J. Molyneux, D. Mukherjee, J. R. Mullaney, E. P. Farina, Silpa S., A. P. Thomson, S. R. Ward

Abstract: We present the first results from the Quasar Feedback Survey, a sample of 42 z<0.2, [O III] luminous AGN (L[O III]>10^42.1 ergs/s) with moderate radio luminosities (i.e. L(1.4GHz)>10^23.4 W/Hz; median L(1.4GHz)=5.9x10^23 W/Hz). Using high spatial resolution (~0.3-1 arcsec), 1.5-6 GHz radio images from the Very Large Array, we find that 67 percent of the sample have spatially extended radio feature… ▽ More We present the first results from the Quasar Feedback Survey, a sample of 42 z<0.2, [O III] luminous AGN (L[O III]>10^42.1 ergs/s) with moderate radio luminosities (i.e. L(1.4GHz)>10^23.4 W/Hz; median L(1.4GHz)=5.9x10^23 W/Hz). Using high spatial resolution (~0.3-1 arcsec), 1.5-6 GHz radio images from the Very Large Array, we find that 67 percent of the sample have spatially extended radio features, on ~1-60 kpc scales. The radio sizes and morphologies suggest that these may be lower radio luminosity versions of compact, radio-loud AGN. By combining the radio-to-infrared excess parameter, spectral index, radio morphology and brightness temperature, we find radio emission in at least 57 percent of the sample that is associated with AGN-related processes (e.g. jets, quasar-driven winds or coronal emission). This is despite only 9.5-21 percent being classified as radio-loud using traditional criteria. The origin of the radio emission in the remainder of the sample is unclear. We find that both the established anti-correlation between radio size and the width of the [O III] line, and the known trend for the most [O III] luminous AGN to be associated with spatially-extended radio emission, also hold for our sample of moderate radio luminosity quasars. These observations add to the growing evidence of a connection between the radio emission and ionised gas in quasar host galaxies. This work lays the foundation for deeper investigations into the drivers and impact of feedback in this unique sample. △ Less

Submitted 26 February, 2021; originally announced March 2021.

Comments: Accepted for publication in MNRAS. Data products from this paper and the survey pilot papers are available through our website: https://blogs.ncl.ac.uk/quasarfeedbacksurvey/. The extensive supplementary material (containing additional figures and information on individual targets) is available for download under "Ancillary files" or by downloading the source file listed under "Other formats"

Testing the Consistency of Dust Laws in SN Ia Host Galaxies: A BayeSN Examination of Foundation DR1

Authors: Stephen Thorp, Kaisey S. Mandel, David O. Jones, Sam M. Ward, Gautham Narayan

Abstract: We apply BayeSN, our new hierarchical Bayesian model for the SEDs of Type Ia supernovae (SNe Ia), to analyse the $griz$ light curves of 157 nearby SNe Ia ($0.015<z<0.08$) from the public Foundation DR1 dataset. We train a new version of BayeSN, continuous from 0.35--0.95 $μ$m, which we use to model the properties of SNe Ia in the rest-frame $z$-band, study the properties of dust in their host gala… ▽ More We apply BayeSN, our new hierarchical Bayesian model for the SEDs of Type Ia supernovae (SNe Ia), to analyse the $griz$ light curves of 157 nearby SNe Ia ($0.015<z<0.08$) from the public Foundation DR1 dataset. We train a new version of BayeSN, continuous from 0.35--0.95 $μ$m, which we use to model the properties of SNe Ia in the rest-frame $z$-band, study the properties of dust in their host galaxies, and construct a Hubble diagram of SN Ia distances determined from full $griz$ light curves. Our $griz$ Hubble diagram has a low total RMS of 0.13 mag using BayeSN, compared to 0.16 mag using SALT2. Additionally, we test the consistency of the dust law $R_V$ between low- and high-mass host galaxies by using our model to fit the full time- and wavelength-dependent SEDs of SNe Ia up to moderate reddening (peak apparent $B-V \lesssim 0.3$). Splitting the population at the median host mass, we find $R_V=2.84\pm0.31$ in low-mass hosts, and $R_V=2.58\pm0.23$ in high-mass hosts, both consistent with the global value of $R_V=2.61\pm0.21$ that we estimate for the full sample. For all choices of mass split we consider, $R_V$ is consistent across the step within $\lesssim1.2σ$. Modelling population distributions of dust laws in low- and high-mass hosts, we find that both subsamples are highly consistent with the full sample's population mean $μ(R_V) = 2.70\pm0.25$ with a 95% upper bound on the population $σ(R_V) < 0.61$. The $R_V$ population means are consistent within $\lesssim1.2σ$. We find that simultaneous fitting of host-mass-dependent dust properties within our hierarchical model does not account for the conventional mass step. △ Less

Submitted 29 October, 2021; v1 submitted 10 February, 2021; originally announced February 2021.

Comments: 22 pages, 14 figures. Accepted for publication in MNRAS

Journal ref: MNRAS 508, 4310-4331 (2021)

Predicting student performance using data from an auto-grading system

Authors: Huanyi Chen, Paul A. S. Ward

Abstract: As online auto-grading systems appear, information obtained from those systems can potentially enable researchers to create predictive models to predict student behaviour and performances. In the University of Waterloo, the ECE 150 (Fundamentals of Programming) Instructional Team wants to get an insight into how to allocate the limited teaching resources better to achieve improved educational outc… ▽ More As online auto-grading systems appear, information obtained from those systems can potentially enable researchers to create predictive models to predict student behaviour and performances. In the University of Waterloo, the ECE 150 (Fundamentals of Programming) Instructional Team wants to get an insight into how to allocate the limited teaching resources better to achieve improved educational outcomes. Currently, the Instructional Team allocates tutoring time in a reactive basis. They help students "as-requested". This approach serves those students with the wherewithal to request help; however, many of the students who are struggling do not reach out for assistance. Therefore, we, as the Research Team, want to explore if we can determine students which need help by looking into the data from our auto-grading system, Marmoset. In this paper, we conducted experiments building decision-tree and linear-regression models with various features extracted from the Marmoset auto-grading system, including passing rate, testcase outcomes, number of submissions and submission time intervals (the time interval between the student's first reasonable submission and the deadline). For each feature, we interpreted the result at the confusion matrix level. Specifically for poor-performance students, we show that the linear-regression model using submission time intervals performs the best among all models in terms of Precision and F-Measure. We also show that for students who are misclassified into poor-performance students, they have the lowest actual grades in the linear-regression model among all models. In addition, we show that for the midterm, the submission time interval of the last assignment before the midterm predicts the midterm performance the most. However, for the final exam, the midterm performance contributes the most on the final exam performance. △ Less

Submitted 1 February, 2021; originally announced February 2021.

Re-examining the Role of Nuclear Fusion in a Renewables-Based Energy Mix

Authors: T. E. G. Nicholas, T. P. Davis, F. Federici, J. E. Leland, B. S. Patel, C. Vincent, S. H. Ward

Abstract: Fusion energy is often regarded as a long-term solution to the world's energy needs. However, even after solving the critical research challenges, engineering and materials science will still impose significant constraints on the characteristics of a fusion power plant. Meanwhile, the global energy grid must transition to low-carbon sources by 2050 to prevent the worst effects of climate change. W… ▽ More Fusion energy is often regarded as a long-term solution to the world's energy needs. However, even after solving the critical research challenges, engineering and materials science will still impose significant constraints on the characteristics of a fusion power plant. Meanwhile, the global energy grid must transition to low-carbon sources by 2050 to prevent the worst effects of climate change. We review three factors affecting fusion's future trajectory: (1) the significant drop in the price of renewable energy, (2) the intermittency of renewable sources and implications for future energy grids, and (3) the recent proposition of intermediate-level nuclear waste as a product of fusion. Within the scenario assumed by our premises, we find that while there remains a clear motivation to develop fusion power plants, this motivation is likely weakened by the time they become available. We also conclude that most current fusion reactor designs do not take these factors into account and, to increase market penetration, fusion research should consider relaxed nuclear waste design criteria, raw material availability constraints and load-following designs with pulsed operation. △ Less

Submitted 14 January, 2021; originally announced January 2021.

Comments: 31 pages, 3 figures, submitted to Energy Policy

Journal ref: Energy Policy, Volume 149, February 2021, 112043

Supersymmetric Analysis of Stochastic Micro-Bending in Optical Waveguides

Authors: Stuart Ward, Rouzbeh Allahverdi, Arash Mafi

Abstract: Micro-bending attenuation in an optical waveguide can be modeled by a Fokker-Planck equation. It is shown that a supersymmetric transformation applied to the Fokker-Planck equation is equivalent to a change in the refractive index profile, resulting in a larger or smaller attenuation. For a broad class of monomial index profiles, it is always possible to obtain an index profile with a larger micro… ▽ More Micro-bending attenuation in an optical waveguide can be modeled by a Fokker-Planck equation. It is shown that a supersymmetric transformation applied to the Fokker-Planck equation is equivalent to a change in the refractive index profile, resulting in a larger or smaller attenuation. For a broad class of monomial index profiles, it is always possible to obtain an index profile with a larger micro-bending attenuation using a supersymmetric transformation. However, obtaining a smaller attenuation is not always possible and is restricted to a subset of index profiles. △ Less

Submitted 6 August, 2021; v1 submitted 24 September, 2020; originally announced September 2020.

Recent Advances in Wearable Sensors with Application in Rehabilitation Motion Analysis

Authors: Shayan Shokri, Shane Ward, Pierre-Amaury M. Anton, Paolo Siffredi, Guglielmo Papetti

Abstract: The increase in world elderly population has significantly underlined the need for continuous health care measurement, specifically in rehabilitation monitoring. The new technologies has enabled people to have in home healthcare services, meanwhile, motion analysis methods are widely used for human activity monitoring as a remote healthcare service. Wearable sensors have indicated promising result… ▽ More The increase in world elderly population has significantly underlined the need for continuous health care measurement, specifically in rehabilitation monitoring. The new technologies has enabled people to have in home healthcare services, meanwhile, motion analysis methods are widely used for human activity monitoring as a remote healthcare service. Wearable sensors have indicated promising results both in convenience and technical performance. These sensors are extensively used in human motion analysis and advancement of wireless communications has intensively contributed to this field. Exploiting wireless technology and wearable sensors contributes to more effective help in emergency cases and has significantly decreased the hospitalization time. This paper reviews the most recent advances in wearable sensors used in motion analysis, specifically in the field of rehabilitation. Firstly, common wearable sensor technologies are introduced and then wearable sensors deploying Carbon Nano Tubes (CNT) are specifically reviewed. The next section is dedicated to sensor fusion in which possibility and performance of integration of new technologies are reviewed. This technique has been widely exploited to bring forth certainty in clinical results. Lastly, the challenges and future possibilities for advancement in motion analysis sensors is discussed. △ Less

Submitted 17 October, 2020; v1 submitted 13 September, 2020; originally announced September 2020.

Comments: 9 pages

Infinite-Parameter ADHM Transform

Authors: R. S. Ward

Abstract: The Atiyah-Drinfeld-Hitchin-Manin (ADHM) transform and its various generalizations are examples of non-linear integral transforms between finite-dimensional moduli spaces. This note describes a natural infinite-dimansional generalization, where the transform becomes a map from boundary data to a family of solutions of the self-duality equations in a domain. The Atiyah-Drinfeld-Hitchin-Manin (ADHM) transform and its various generalizations are examples of non-linear integral transforms between finite-dimensional moduli spaces. This note describes a natural infinite-dimansional generalization, where the transform becomes a map from boundary data to a family of solutions of the self-duality equations in a domain. △ Less

Submitted 26 July, 2020; v1 submitted 2 June, 2020; originally announced June 2020.

Comments: 13 pages, for special issue of Quarterly Journal of Mathematics dedicated to Michael Atiyah; v2 minor amendments, accepted version

Testing for complete spatial randomness on three dimensional bounded convex shapes

Authors: Scott Ward, Edward A. K. Cohen, Niall Adams

Abstract: There is currently a gap in theory for point patterns that lie on the surface of objects, with researchers focusing on patterns that lie in a Euclidean space, typically planar and spatial data. Methodology for planar and spatial data thus relies on Euclidean geometry and is therefore inappropriate for analysis of point patterns observed in non-Euclidean spaces. Recently, there has been extensions… ▽ More There is currently a gap in theory for point patterns that lie on the surface of objects, with researchers focusing on patterns that lie in a Euclidean space, typically planar and spatial data. Methodology for planar and spatial data thus relies on Euclidean geometry and is therefore inappropriate for analysis of point patterns observed in non-Euclidean spaces. Recently, there has been extensions to the analysis of point patterns on a sphere, however, many other shapes are left unexplored. This is in part due to the challenge of defining the notion of stationarity for a point process existing on such a space due to the lack of rotational and translational isometries. Here, we construct functional summary statistics for Poisson processes defined on convex shapes in three dimensions. Using the Mapping Theorem, a Poisson process can be transformed from any convex shape to a Poisson process on the unit sphere which has rotational symmetries that allow for functional summary statistics to be constructed. We present the first and second order properties of such summary statistics and demonstrate how they can be used to test whether an observed pattern exhibits complete spatial randomness or spatial preference on the original convex space. A study of the Type I and II errors of our test statistics are explored through simulations on ellipsoids of varying dimensions. △ Less

Submitted 10 February, 2020; originally announced February 2020.

MSC Class: 62H11

The FAIR Funder pilot programme to make it easy for funders to require and for grantees to produce FAIR Data

Authors: P. Wittenburg, H. Pergl Sustkova, A. Montesanti, S. M. Bloemers, S. H. de Waard, M. A. Musen, J. B. Graybeal, K. M. Hettne, A. Jacobsen, R. Pergl, R. W. W. Hooft, C. Staiger, C. W. G. van Gelder, S. L. Knijnenburg, A. C. van Arkel, B. Meerman, M. D. Wilkinson, S-A Sansone, P. Rocca-Serra, P. McQuilton, A. N. Gonzalez-Beltran, G. J. C. Aben, P. Henning, S. Alencar, C. Ribeiro , et al. (35 additional authors not shown)

Abstract: There is a growing acknowledgement in the scientific community of the importance of making experimental data machine findable, accessible, interoperable, and reusable (FAIR). Recognizing that high quality metadata are essential to make datasets FAIR, members of the GO FAIR Initiative and the Research Data Alliance (RDA) have initiated a series of workshops to encourage the creation of Metadata for… ▽ More There is a growing acknowledgement in the scientific community of the importance of making experimental data machine findable, accessible, interoperable, and reusable (FAIR). Recognizing that high quality metadata are essential to make datasets FAIR, members of the GO FAIR Initiative and the Research Data Alliance (RDA) have initiated a series of workshops to encourage the creation of Metadata for Machines (M4M), enabling any self-identified stakeholder to define and promote the reuse of standardized, comprehensive machine-actionable metadata. The funders of scientific research recognize that they have an important role to play in ensuring that experimental results are FAIR, and that high quality metadata and careful planning for FAIR data stewardship are central to these goals. We describe the outcome of a recent M4M workshop that has led to a pilot programme involving two national science funders, the Health Research Board of Ireland (HRB) and the Netherlands Organisation for Health Research and Development (ZonMW). These funding organizations will explore new technologies to define at the time that a request for proposals is issued the minimal set of machine-actionable metadata that they would like investigators to use to annotate their datasets, to enable investigators to create such metadata to help make their data FAIR, and to develop data-stewardship plans that ensure that experimental data will be managed appropriately abiding by the FAIR principles. The FAIR Funders design envisions a data-management workflow having seven essential stages, where solution providers are openly invited to participate. The initial pilot programme will launch using existing computer-based tools of those who attended the M4M Workshop. △ Less

Submitted 6 March, 2019; v1 submitted 26 February, 2019; originally announced February 2019.

Comments: This is a pre-print of the FAIR Funders pilot, an outcome of the first Metadata for Machines workshop, see: https://www.go-fair.org/resources/go-fair-workshop-series/metadata-for-machines-workshops/. Corresponding author: E. A Schultes, ORCID 0000-0001-8888-635X

Thermal control of spin excitations in the coupled Ising-chain material RbCoCl3

Authors: M. Mena, N. Hänni, S. Ward, E. Hirtenlechner, R. Bewley, C. Hubig, U. Schollwöck, B. Normand, K. W. Krämer, D. F. McMorrow, Ch. Rüegg

Abstract: We have used neutron spectroscopy to investigate the spin dynamics of the quantum (S = 1/2) antiferromagnetic Ising chains in RbCoCl3. The structure and magnetic interactions in this material conspire to produce two magnetic phase transitions at low temperatures, presenting an ideal opportunity for thermal control of the chain environment. The high-resolution spectra we measure of two-domain-wall… ▽ More We have used neutron spectroscopy to investigate the spin dynamics of the quantum (S = 1/2) antiferromagnetic Ising chains in RbCoCl3. The structure and magnetic interactions in this material conspire to produce two magnetic phase transitions at low temperatures, presenting an ideal opportunity for thermal control of the chain environment. The high-resolution spectra we measure of two-domain-wall excitations therefore characterize precisely both the continuum response of isolated chains and the "Zeeman-ladder" bound states of chains in three different effective staggered fields in one and the same material. We apply an extended Matsubara formalism to obtain a quantitative description of the entire dataset, Monte Carlo simulations to interpret the magnetic order, and finite-temperature DMRG calculations to fit the spectral features of all three phases. △ Less

Submitted 30 June, 2020; v1 submitted 17 November, 2018; originally announced November 2018.

Comments: 6 pages, 4 figures. Supplemental material included

Journal ref: Phys. Rev. Lett. 124, 257201 (2020)

Quantum magnetism in molecular spin ladders probed with muon-spin spectroscopy

Authors: T. Lancaster, F. Xiao, B. M. Huddart, R. C. Williams, F. L. Pratt, S. J. Blundell, S. J. Clark, R. Scheuermann, T. Goko, S. Ward, J. L. Manson, Ch. Rüegg, K. W. Krämer

Abstract: We present the results of muon-spin spectroscopy ($μ^{+}$SR) measurements on the molecular spin ladder system (Hpip)$_{2}$CuBr$_{4(1-x)}$Cl$_{4x}$, [Hpip=(C$_{5}$H$_{12}$N)]. Using transverse field $μ^{+}$SR we are able to identify characteristic behaviour in each of the regions of the phase diagram of the $x=0$ strong-rung spin ladder system (Hpip)$_{2}$CuBr$_4$. Comparison of our results to thos… ▽ More We present the results of muon-spin spectroscopy ($μ^{+}$SR) measurements on the molecular spin ladder system (Hpip)$_{2}$CuBr$_{4(1-x)}$Cl$_{4x}$, [Hpip=(C$_{5}$H$_{12}$N)]. Using transverse field $μ^{+}$SR we are able to identify characteristic behaviour in each of the regions of the phase diagram of the $x=0$ strong-rung spin ladder system (Hpip)$_{2}$CuBr$_4$. Comparison of our results to those of the dimer-based molecular magnet Cu(pyz)(gly)(ClO$_{4}$) shows several common features. We locate the crossovers in partially disordered (Hpip)$_{2}$CuBr$_{4(1-x)}$Cl$_{4x}$ ($x=0.05$), where a region of behaviour intermediate between quantum disordered and Luttinger liquid-like is identified. Our interpretation of the results incorporates an analysis of the probable muon stopping states in (Hpip)$_{2}$CuBr$_4$ based on density functional calculations and suggests how the muon plus its local distortion can lead to a local probe unit with good sensitivity to the magnetic state. Using longitudinal field $μ^{+}$SR we compare the dynamic response of the $x=1$ strong-rung material (Hpip)$_{2}$CuCl$_{4}$ to that of the strong-leg material (C$_{7}$H$_{10}$N)$_{2}$CuBr$_{4}$ (known as DIMPY) and demonstrate that our results are in agreement with predictions based on interacting fermionic quasiparticle excitations in these materials. △ Less

Submitted 25 June, 2018; originally announced June 2018.

Comments: 16 pages, 13 figures

Hopf solitons on compact manifolds

Authors: R. S. Ward

Abstract: Hopf solitons in the Skyrme-Faddeev system on $R^3$ typically have a complicated structure, in particular when the Hopf number Q is large. By contrast, if we work on a compact 3-manifold M, and the energy functional consists only of the Skyrme term (the strong-coupling limit), then the picture simplifies. There is a topological lower bound $E\geq Q$ on the energy, and the local minima of E can loo… ▽ More Hopf solitons in the Skyrme-Faddeev system on $R^3$ typically have a complicated structure, in particular when the Hopf number Q is large. By contrast, if we work on a compact 3-manifold M, and the energy functional consists only of the Skyrme term (the strong-coupling limit), then the picture simplifies. There is a topological lower bound $E\geq Q$ on the energy, and the local minima of E can look simple even for large Q. The aim here is to describe and investigate some of these solutions, when M is $S^3$, $T^3$ or $S^2 \times S^1$. In addition, we review the more elementary baby-Skyrme system, with M being $S^2$ or $T^2$. △ Less

Submitted 2 February, 2018; originally announced February 2018.

Comments: 12 pages, 2 figures, to be published in JMP

Smaller desert dust cooling effect estimated from analysis of dust size and abundance

Authors: Jasper F. Kok, David A. Ridley, Qing Zhou, Ron L. Miller, Chun Zhao, Colette L. Heald, Daniel S. Ward, Samuel Albani, Karsten Haustein

Abstract: Desert dust aerosols affect Earth's global energy balance through direct interactions with radiation, and through indirect interactions with clouds and ecosystems. But the magnitudes of these effects are so uncertain that it remains unclear whether atmospheric dust has a net warming or cooling effect on global climate. Consequently, it is still uncertain whether large changes in atmospheric dust l… ▽ More Desert dust aerosols affect Earth's global energy balance through direct interactions with radiation, and through indirect interactions with clouds and ecosystems. But the magnitudes of these effects are so uncertain that it remains unclear whether atmospheric dust has a net warming or cooling effect on global climate. Consequently, it is still uncertain whether large changes in atmospheric dust loading over the past century have slowed or accelerated anthropogenic climate change, or what the effects of potential future changes in dust loading will be. Here we present an analysis of the size and abundance of dust aerosols to constrain the direct radiative effect of dust. Using observational data on dust abundance, in situ measurements of dust optical properties and size distribution, and climate and atmospheric chemical transport model simulations of dust lifetime, we find that the dust found in the atmosphere is substantially coarser than represented in current global climate models. Since coarse dust warms climate, the global dust direct radiative effect is likely to be less cooling than the ~-0.4 W/m2 estimated by models in a current global aerosol model ensemble. Instead, we constrain the dust direct radiative effect to a range between -0.48 and +0.20 W/m2, which includes the possibility that dust causes a net warming of the planet. △ Less

Submitted 20 October, 2017; originally announced October 2017.

Comments: Published as article in Nature Geoscience

Journal ref: Nature Geoscience, 10, 274-278

The Final Chapter In The Saga Of YIG

Authors: A. J. Princep, R. A. Ewings, S. Ward, S. Tóth, C. Dubs, D. Prabhakaran, A. T. Boothroyd

Abstract: The magnetic insulator Yttrium Iron Garnet can be grown with exceptional quality, has a ferrimagnetic transition temperature of nearly 600 K, and is used in microwave and spintronic devices that can operate at room temperature. The most accurate prior measurements of the magnon spectrum date back nearly 40 years, but cover only 3 of the lowest energy modes out of 20 distinct magnon branches. Here… ▽ More The magnetic insulator Yttrium Iron Garnet can be grown with exceptional quality, has a ferrimagnetic transition temperature of nearly 600 K, and is used in microwave and spintronic devices that can operate at room temperature. The most accurate prior measurements of the magnon spectrum date back nearly 40 years, but cover only 3 of the lowest energy modes out of 20 distinct magnon branches. Here we have used time-of-flight inelastic neutron scattering to measure the full magnon spectrum throughout the Brillouin zone. We find that the existing model of the excitation spectrum, well known from an earlier work titled "The Saga of YIG", fails to describe the optical magnon modes. Using a very general spin Hamiltonian, we show that the magnetic interactions are both longer-ranged and more complex than was previously understood. The results provide the basis for accurate microscopic models of the finite temperature magnetic properties of Yttrium Iron Garnet, necessary for next-generation electronic devices. △ Less

Submitted 18 May, 2017; originally announced May 2017.

Comments: 10 pages, 3 figures, 4 supplementary figures, 1 table, 1 supplementary table

Journal ref: npj Quantum Materials, volume 2, Article number: 63 (2017)

Bound states and field-polarized Haldane modes in a quantum spin ladder

Authors: S. Ward, P. Bouillot, C. Kollath, T. Giamarchi, K. P. Schmidt, B. Normand, K. W. Krämer, D. Biner, R. Bewley, T. Guidi, M. Boehm, D. F. McMorrow, Ch. Rüegg

Abstract: The challenge of one-dimensional systems is to understand their physics beyond the level of known elementary excitations. By high-resolution neutron spectroscopy in a quantum spin ladder material, we probe the leading multiparticle excitation by characterizing the two-magnon bound state at zero field. By applying high magnetic fields, we create and select the singlet (longitudinal) and triplet (tr… ▽ More The challenge of one-dimensional systems is to understand their physics beyond the level of known elementary excitations. By high-resolution neutron spectroscopy in a quantum spin ladder material, we probe the leading multiparticle excitation by characterizing the two-magnon bound state at zero field. By applying high magnetic fields, we create and select the singlet (longitudinal) and triplet (transverse) excitations of the fully spin-polarized ladder, which have not been observed previously and are close analogs of the modes anticipated in a polarized Haldane chain. Theoretical modelling of the dynamical response demonstrates our complete quantitative understanding of these states. △ Less

Submitted 1 September, 2016; originally announced September 2016.

Comments: 6 pages, 3 figures plus supplementary material 7 pages 5 figures

Journal ref: Phys. Rev. Lett. 118, 177202 (2017)

Integrable (2k)-Dimensional Hitchin Equations

Authors: R. S. Ward

Abstract: This letter describes a completely-integrable system of Yang-Mills-Higgs equations which generalizes the Hitchin equations on a Riemann surface to arbitrary k-dimensional complex manifolds. The system arises as a dimensional reduction of a set of integrable Yang-Mills equations in 4k real dimensions. Our integrable system implies other generalizations such as the Simpson equations and the non-abel… ▽ More This letter describes a completely-integrable system of Yang-Mills-Higgs equations which generalizes the Hitchin equations on a Riemann surface to arbitrary k-dimensional complex manifolds. The system arises as a dimensional reduction of a set of integrable Yang-Mills equations in 4k real dimensions. Our integrable system implies other generalizations such as the Simpson equations and the non-abelian Seiberg-Witten equations. Some simple solutions in the k=2 case are described. △ Less

Submitted 25 April, 2016; originally announced April 2016.

Comments: 11 pages, 1 figure

Symmetric Instantons and Discrete Hitchin Equations

Authors: R. S. Ward

Abstract: Self-dual Yang-Mills instantons on $R^4$ correspond to algebraic ADHM data. The ADHM equations for $S^1$-symmetric instantons give a one-dimensional integrable lattice system, which may be viewed as an discretization of the Nahm equations. In this note, we see that generalized ADHM data for $T^2$-symmetric instantons gives an integrable two-dimensional lattice system, which may be viewed as a disc… ▽ More Self-dual Yang-Mills instantons on $R^4$ correspond to algebraic ADHM data. The ADHM equations for $S^1$-symmetric instantons give a one-dimensional integrable lattice system, which may be viewed as an discretization of the Nahm equations. In this note, we see that generalized ADHM data for $T^2$-symmetric instantons gives an integrable two-dimensional lattice system, which may be viewed as a discrete version of the Hitchin equations. △ Less

Submitted 11 February, 2016; v1 submitted 30 September, 2015; originally announced September 2015.

Comments: 9 pages, 1 figure, simplified version, accepted for publication in Journal of Integrable Systems

Report number: DCPT-15/59

Geometry of Solutions of Hitchin Equations on R^2

Authors: R. S. Ward

Abstract: We study smooth SU(2) solutions of the Hitchin equations on R^2, with the determinant of the complex Higgs field being a polynomial of degree n. When n>=3, there are moduli spaces of solutions, in the sense that the natural L^2 metric is well-defined on a subset of the parameter space. We examine rotationally-symmetric solutions for n=1 and n=2, and then focus on the n=3 case, elucidating the modu… ▽ More We study smooth SU(2) solutions of the Hitchin equations on R^2, with the determinant of the complex Higgs field being a polynomial of degree n. When n>=3, there are moduli spaces of solutions, in the sense that the natural L^2 metric is well-defined on a subset of the parameter space. We examine rotationally-symmetric solutions for n=1 and n=2, and then focus on the n=3 case, elucidating the moduli and describing the asymptotic geometry as well as the geometry of two totally-geodesic surfaces. △ Less

Submitted 4 January, 2016; v1 submitted 22 April, 2015; originally announced April 2015.

Comments: 13 pages, various minor changes, added references, accepted for publication in Nonlinearity

Report number: DCPT-15/21

Building the next pyramid

Authors: Joseph West, Greg Gallagher, Kevin Waters, Stephen Ward, Tia Ward

Abstract: The results of experimental tests of a novel method for moving large (pyramid construction size) stone blocks by rolling them are presented. The method is implemented by tying 12 identical rods of appropriately chosen radius to the faces of the block forming a rough dodecagon prism. Experiments using a 1,000 kg block show that it can be moved across level open ground with a dynamic coefficient of… ▽ More The results of experimental tests of a novel method for moving large (pyramid construction size) stone blocks by rolling them are presented. The method is implemented by tying 12 identical rods of appropriately chosen radius to the faces of the block forming a rough dodecagon prism. Experiments using a 1,000 kg block show that it can be moved across level open ground with a dynamic coefficient of friction of less than 0.06. This value is a factor of five lower than that obtained for dragging the block, and the best values reported for dragging by others, at 0.3. the results are more dramatic than those obtained on smaller scale experiments on a 29.6 kg block, also reported here. For full scale pyramid blocks, the wooden "rods" woudl need to be posts of order 30 cm in diameter, similar in size to those used as masts on ships in the Nile. △ Less

Submitted 25 February, 2015; originally announced February 2015.

Comments: 12 pages, 3 figures. arXiv admin note: substantial text overlap with arXiv:1408.3603

Academic Cloud Computing Research: Five Pitfalls and Five Opportunities

Authors: Adam Barker, Blesson Varghese, Jonathan Stuart Ward, Ian Sommerville

Abstract: This discussion paper argues that there are five fundamental pitfalls, which can restrict academics from conducting cloud computing research at the infrastructure level, which is currently where the vast majority of academic research lies. Instead academics should be conducting higher risk research, in order to gain understanding and open up entirely new areas. We call for a renewed mindset and… ▽ More This discussion paper argues that there are five fundamental pitfalls, which can restrict academics from conducting cloud computing research at the infrastructure level, which is currently where the vast majority of academic research lies. Instead academics should be conducting higher risk research, in order to gain understanding and open up entirely new areas. We call for a renewed mindset and argue that academic research should focus less upon physical infrastructure and embrace the abstractions provided by clouds through five opportunities: user driven research, new programming models, PaaS environments, and improved tools to support elasticity and large-scale debugging. The objective of this paper is to foster discussion, and to define a roadmap forward, which will allow academia to make longer-term impacts to the cloud computing community. △ Less

Submitted 19 June, 2014; originally announced June 2014.

Comments: Accepted and presented at the 6th USENIX Workshop on Hot Topics in Cloud Computing (HotCloud'14)

Dynamics of monopole walls

Authors: R. Maldonado, R. S. Ward

Abstract: The moduli space of centred Bogomolny-Prasad-Sommmerfield 2-monopole fields is a 4-dimensional manifold M with a natural metric, and the geodesics on M correspond to slow-motion monopole dynamics. The best-known case is that of monopoles on R^3, where M is the Atiyah-Hitchin space. More recently, the case of monopoles periodic in one direction (monopole chains) was studied a few years ago. Our aim… ▽ More The moduli space of centred Bogomolny-Prasad-Sommmerfield 2-monopole fields is a 4-dimensional manifold M with a natural metric, and the geodesics on M correspond to slow-motion monopole dynamics. The best-known case is that of monopoles on R^3, where M is the Atiyah-Hitchin space. More recently, the case of monopoles periodic in one direction (monopole chains) was studied a few years ago. Our aim in this note is to investigate M for doubly-periodic fields, which may be visualized as monopole walls. We identify some of the geodesics on M as fixed-point sets of discrete symmetries, and interpret these in terms of monopole scattering and bound orbits, concentrating on novel features that arise as a consequence of the periodicity. △ Less

Submitted 19 May, 2014; originally announced May 2014.

Comments: 12 pages, 3 figures

Report number: DCPT-14/19

Journal ref: Phys. Lett. B 734 (2014) 328

Geometry of Periodic Monopoles

Authors: Rafael Maldonado, R S Ward

Abstract: BPS monopoles on $\mathbb{R}^2\times S^1$ correspond, via the generalized Nahm transform, to certain solutions of the Hitchin equations on the cylinder $\mathbb{R}\times S^1$. The moduli space M of two monopoles with their centre-of-mass fixed is a 4-dimensional manifold with a natural hyperkähler metric, and its geodesics correspond to slow-motion monopole scattering. The purpose of this paper is… ▽ More BPS monopoles on $\mathbb{R}^2\times S^1$ correspond, via the generalized Nahm transform, to certain solutions of the Hitchin equations on the cylinder $\mathbb{R}\times S^1$. The moduli space M of two monopoles with their centre-of-mass fixed is a 4-dimensional manifold with a natural hyperkähler metric, and its geodesics correspond to slow-motion monopole scattering. The purpose of this paper is to study the geometry of M in terms of the Nahm/Hitchin data, i.e. in terms of structures on $\mathbb{R}\times S^1$. In particular, we identify the moduli, derive the asymptotic metric on M, and discuss several geodesic surfaces and geodesics on M. The latter include novel examples of monopole dynamics. △ Less

Submitted 11 July, 2014; v1 submitted 26 September, 2013; originally announced September 2013.

Comments: 17 pages, 4 figures

Report number: DCPT-13/35

Journal ref: Phys. Rev. D 88, 125013 (2013)

Undefined By Data: A Survey of Big Data Definitions

Authors: Jonathan Stuart Ward, Adam Barker

Abstract: The term big data has become ubiquitous. Owing to a shared origin between academia, industry and the media there is no single unified definition, and various stakeholders provide diverse and often contradictory definitions. The lack of a consistent definition introduces ambiguity and hampers discourse relating to big data. This short paper attempts to collate the various definitions which have gai… ▽ More The term big data has become ubiquitous. Owing to a shared origin between academia, industry and the media there is no single unified definition, and various stakeholders provide diverse and often contradictory definitions. The lack of a consistent definition introduces ambiguity and hampers discourse relating to big data. This short paper attempts to collate the various definitions which have gained some degree of traction and to furnish a clear and concise definition of an otherwise ambiguous term. △ Less

Submitted 20 September, 2013; originally announced September 2013.

Comments: Big data definition paper, 2 pages

Spin ladders and quantum simulators for Luttinger liquids

Authors: S. Ward, P. Bouillot, H. Ryll, K. Kiefer, K. W. Krämer, Ch. Rüegg, C. Kollath, T. Giamarchi

Abstract: Magnetic insulators have proven to be usable as quantum simulators for itinerant interacting quantum systems. In particular the compound (C$_{5}$H$_{12}$N)$_{2}$CuBr$_{4}$ (short (Hpip)$_{2}$CuBr$_{4}$) was shown to be a remarkable realization of a Tomonaga-Luttinger liquid (TLL) and allowed to quantitatively test the TLL theory. Substitution weakly disorders this class of compounds and allows thu… ▽ More Magnetic insulators have proven to be usable as quantum simulators for itinerant interacting quantum systems. In particular the compound (C$_{5}$H$_{12}$N)$_{2}$CuBr$_{4}$ (short (Hpip)$_{2}$CuBr$_{4}$) was shown to be a remarkable realization of a Tomonaga-Luttinger liquid (TLL) and allowed to quantitatively test the TLL theory. Substitution weakly disorders this class of compounds and allows thus to use them to tackle questions pertaining to the effect of disorder in TLL as well, such as the formation of the Bose glass. As a first step in this direction we present in this paper a study of the properties of the related (Hpip)$_{2}$CuCl$_{4}$ compound. We determine the exchange couplings and compute the temperature and magnetic field dependence of the specific heat, using a finite temperature Density Matrix Renormalization group (DMRG) procedure. Comparison with the measured specific heat at zero magnetic field confirms the exchange parameters and Hamiltonian for the (Hpip)$_{2}$CuCl$_{4}$ compound, giving the basis needed to start studying the disorder effects. △ Less

Submitted 12 August, 2013; originally announced August 2013.

Comments: Part of the special issue "Physics in one dimension", A van Houselt, J Schäfer, H J W Zandvliet and R Claessen ed

Journal ref: J. Phys.: Condens. Matter 25 014004 (2013)

A Cloud Computing Survey: Developments and Future Trends in Infrastructure as a Service Computing

Authors: Jonathan Stuart Ward, Adam Barker

Abstract: Cloud computing is a recent paradigm based around the notion of delivery of resources via a service model over the Internet. Despite being a new paradigm of computation, cloud computing owes its origins to a number of previous paradigms. The term cloud computing is well defined and no longer merits rigorous taxonomies to furnish a definition. Instead this survey paper considers the past, present a… ▽ More Cloud computing is a recent paradigm based around the notion of delivery of resources via a service model over the Internet. Despite being a new paradigm of computation, cloud computing owes its origins to a number of previous paradigms. The term cloud computing is well defined and no longer merits rigorous taxonomies to furnish a definition. Instead this survey paper considers the past, present and future of cloud computing. As an evolution of previous paradigms, we consider the predecessors to cloud computing and what significance they still hold to cloud services. Additionally we examine the technologies which comprise cloud computing and how the challenges and future developments of these technologies will influence the field. Finally we examine the challenges that limit the growth, application and development of cloud computing and suggest directions required to overcome these challenges in order to further the success of cloud computing. △ Less

Submitted 6 June, 2013; originally announced June 2013.