Showing 1–50 of 98 results for author: Joseph, R

Magnetospheric origin of a fast radio burst constrained using scintillation

Authors: Kenzie Nimmo, Ziggy Pleunis, Paz Beniamini, Pawan Kumar, Adam E. Lanman, D. Z. Li, Robert Main, Mawson W. Sammons, Shion Andrew, Mohit Bhardwaj, Shami Chatterjee, Alice P. Curtin, Emmanuel Fonseca, B. M. Gaensler, Ronniy C. Joseph, Zarif Kader, Victoria M. Kaspi, Mattias Lazda, Calvin Leung, Kiyoshi W. Masui, Ryan Mckinven, Daniele Michilli, Ayush Pandhi, Aaron B. Pearlman, Masoud Rafiei-Ravandi , et al. (4 additional authors not shown)

Abstract: Fast radio bursts (FRBs) are micro-to-millisecond duration radio transients that originate mostly from extragalactic distances. The emission mechanism responsible for these high luminosity, short duration transients remains debated. The models are broadly grouped into two classes: physical processes that occur within close proximity to a central engine; and central engines that release energy whic… ▽ More Fast radio bursts (FRBs) are micro-to-millisecond duration radio transients that originate mostly from extragalactic distances. The emission mechanism responsible for these high luminosity, short duration transients remains debated. The models are broadly grouped into two classes: physical processes that occur within close proximity to a central engine; and central engines that release energy which moves to large radial distances and subsequently interacts with surrounding media producing radio waves. The expected emission region sizes are notably different between these two types of models. FRB emission size constraints can therefore be used to distinguish between these competing models and inform on the physics responsible. Here we present the measurement of two mutually coherent scintillation scales in the frequency spectrum of FRB 20221022A: one originating from a scattering screen located within the Milky Way, and the second originating from a scattering screen located within its host galaxy or local environment. We use the scattering media as an astrophysical lens to constrain the size of the lateral emission region, $R_{\star\mathrm{obs}} \lesssim 3\times10^{4}$ km. We find that this is inconsistent with the expected emission sizes for the large radial distance models, and is more naturally explained with an emission process that operates within or just beyond the magnetosphere of a central compact object. Recently, FRB 20221022A was found to exhibit an S-shaped polarisation angle swing, supporting a magnetospheric emission process. The scintillation results presented in this work independently support this conclusion, while highlighting scintillation as a useful tool in our understanding of FRB emission physics and progenitors. △ Less

Submitted 16 June, 2024; originally announced June 2024.

Comments: Submitted, comments welcome

Transportation and Systems Analysis Collaborations in Support of a Federal Consolidated Interim Storage Facility

Authors: Robby Joseph, Harish Gadey, Brian Craig, Lucas Vander Wal, Mark Abkowitz, Robert Claypool, Riley Cumberland, Caitlin Condon, Miriam Juckett, Steve Maheras, Gordon Petersen

Abstract: The Integrated Waste Management program under the Office of Nuclear Energy is planning for the future transportation, storage, and eventual disposal of spent nuclear fuel and high-level radioactive waste from nuclear power plant and waste custodian sites across the United States. To better enable informed decision-making regarding the back end of the nuclear fuel cycle, the IWM program has been sp… ▽ More The Integrated Waste Management program under the Office of Nuclear Energy is planning for the future transportation, storage, and eventual disposal of spent nuclear fuel and high-level radioactive waste from nuclear power plant and waste custodian sites across the United States. To better enable informed decision-making regarding the back end of the nuclear fuel cycle, the IWM program has been sponsoring the development and application of system analysis tools capable of analyzing various options for managing SNF and HLW. With these tools, integrated waste management system architecture analyses are being conducted to support the future deployment of a comprehensive nuclear waste management system that considers all major back-end aspects of the nuclear fuel cycle i.e., transportation, storage, and disposal. System analyses and assessments typically use these modeling and simulation tools to investigate implications of changes in various assumptions and parameters such as acceptance rates, receipt logic, facility capacities and capabilities, use of standardized canisters, start and stop dates of facilities, etc. △ Less

Submitted 21 March, 2024; originally announced March 2024.

Comments: Presented at the 2024 Waste Management Symposia. arXiv admin note: text overlap with arXiv:2306.12319

Implementation of Parallel Process Execution in the Next Generation System Analysis Model

Authors: Harish Gadey, Lucas Vander Wal, Robert Joseph

Abstract: The United States DOE Office of Integrated Waste Management program is planning for the transportation, storage, and eventual disposal of spent nuclear fuel and high-level radioactive waste from nuclear power plant sites across the United States. The Next Generation System Analysis Model is an agent-based simulation toolkit that is used for system level simulation and analysis of the SNF inventory… ▽ More The United States DOE Office of Integrated Waste Management program is planning for the transportation, storage, and eventual disposal of spent nuclear fuel and high-level radioactive waste from nuclear power plant sites across the United States. The Next Generation System Analysis Model is an agent-based simulation toolkit that is used for system level simulation and analysis of the SNF inventory in the United States. This tool was developed as part of a collaborative effort between Argonne National Laboratory and Oak Ridge National Laboratory. The analyst using NGSAM has the ability to define several factors like the number of storage facilities, capacity at each facility, transportation schedules, shipment rates, and other conditions. The primary purpose of NGSAM is to provide the system analyst with the tools to model the integrated waste management system and gain insights on waste management alternatives, impact of storage choices, generating cost estimates, and developing an integrated approach with emphasis on flexibility. △ Less

Submitted 14 March, 2024; originally announced March 2024.

Comments: Presented at the '2024 Waste Management Symposia'

Development and validation of an artificial intelligence model to accurately predict spinopelvic parameters

Authors: Edward S. Harake, Joseph R. Linzey, Cheng Jiang, Rushikesh S. Joshi, Mark M. Zaki, Jaes C. Jones, Siri S. Khalsa, John H. Lee, Zachary Wilseck, Jacob R. Joseph, Todd C. Hollon, Paul Park

Abstract: Objective. Achieving appropriate spinopelvic alignment has been shown to be associated with improved clinical symptoms. However, measurement of spinopelvic radiographic parameters is time-intensive and interobserver reliability is a concern. Automated measurement tools have the promise of rapid and consistent measurements, but existing tools are still limited by some degree of manual user-entry re… ▽ More Objective. Achieving appropriate spinopelvic alignment has been shown to be associated with improved clinical symptoms. However, measurement of spinopelvic radiographic parameters is time-intensive and interobserver reliability is a concern. Automated measurement tools have the promise of rapid and consistent measurements, but existing tools are still limited by some degree of manual user-entry requirements. This study presents a novel artificial intelligence (AI) tool called SpinePose that automatically predicts spinopelvic parameters with high accuracy without the need for manual entry. Methods. SpinePose was trained and validated on 761 sagittal whole-spine X-rays to predict sagittal vertical axis (SVA), pelvic tilt (PT), pelvic incidence (PI), sacral slope (SS), lumbar lordosis (LL), T1-pelvic angle (T1PA), and L1-pelvic angle (L1PA). A separate test set of 40 X-rays was labeled by 4 reviewers, including fellowship-trained spine surgeons and a fellowship-trained radiologist with neuroradiology subspecialty certification. Median errors relative to the most senior reviewer were calculated to determine model accuracy on test images. Intraclass correlation coefficients (ICC) were used to assess inter-rater reliability. Results. SpinePose exhibited the following median (interquartile range) parameter errors: SVA: 2.2(2.3)mm, p=0.93; PT: 1.3(1.2)°, p=0.48; SS: 1.7(2.2)°, p=0.64; PI: 2.2(2.1)°, p=0.24; LL: 2.6(4.0)°, p=0.89; T1PA: 1.1(0.9)°, p=0.42; and L1PA: 1.4(1.6)°, p=0.49. Model predictions also exhibited excellent reliability at all parameters (ICC: 0.91-1.0). Conclusions. SpinePose accurately predicted spinopelvic parameters with excellent reliability comparable to fellowship-trained spine surgeons and neuroradiologists. Utilization of predictive AI tools in spinal imaging can substantially aid in patient selection and surgical planning. △ Less

Submitted 8 February, 2024; originally announced February 2024.

Comments: 10 pages, 5 figures, to appear in Journal of Neurosurgery: Spine

Factor Importance Ranking and Selection using Total Indices

Authors: Chaofan Huang, V. Roshan Joseph

Abstract: Factor importance measures the impact of each feature on output prediction accuracy. Many existing works focus on the model-based importance, but an important feature in one learning algorithm may hold little significance in another model. Hence, a factor importance measure ought to characterize the feature's predictive potential without relying on a specific prediction algorithm. Such algorithm-a… ▽ More Factor importance measures the impact of each feature on output prediction accuracy. Many existing works focus on the model-based importance, but an important feature in one learning algorithm may hold little significance in another model. Hence, a factor importance measure ought to characterize the feature's predictive potential without relying on a specific prediction algorithm. Such algorithm-agnostic importance is termed as intrinsic importance in Williamson et al. (2023), but their estimator again requires model fitting. To bypass the modeling step, we present the equivalence between predictiveness potential and total Sobol' indices from global sensitivity analysis, and introduce a novel consistent estimator that can be directly estimated from noisy data. Integrating with forward selection and backward elimination gives rise to FIRST, Factor Importance Ranking and Selection using Total (Sobol') indices. Extensive simulations are provided to demonstrate the effectiveness of FIRST on regression and binary classification problems, and a clear advantage over the state-of-the-art methods. △ Less

Submitted 11 January, 2024; v1 submitted 1 January, 2024; originally announced January 2024.

Measurement of the Isolated Nuclear Two-Photon Decay in $^{72}\mathrm{Ge}$

Authors: D. Freire-Fernández, W. Korten, R. J. Chen, S. Litvinov, Yu. A. Litvinov, M. S. Sanjari, H. Weick, F. C. Akinci, H. M. Albers, M. Armstrong, A. Banerjee, K. Blaum, C. Brandau, B. A. Brown, C. G. Bruno, J. J. Carroll, X. Chen, Ch. J. Chiara, M. L. Cortes, S. F. Dellmann, I. Dillmann, D. Dmytriiev, O. Forstner, H. Geissel, J. Glorius , et al. (35 additional authors not shown)

Abstract: The nuclear two-photon or double-gamma ($2γ$) decay is a second-order electromagnetic process whereby a nucleus in an excited state emits two gamma rays simultaneously. To be able to directly measure the $2γ$ decay rate in the low-energy regime below the electron-positron pair-creation threshold, we combined the isochronous mode of a storage ring with Schottky resonant cavities. The newly develope… ▽ More The nuclear two-photon or double-gamma ($2γ$) decay is a second-order electromagnetic process whereby a nucleus in an excited state emits two gamma rays simultaneously. To be able to directly measure the $2γ$ decay rate in the low-energy regime below the electron-positron pair-creation threshold, we combined the isochronous mode of a storage ring with Schottky resonant cavities. The newly developed technique can be applied to isomers with excitation energies down to $\sim100$\,keV and half-lives as short as $\sim10$\,ms. The half-life for the $2γ$ decay of the first-excited $0^+$ state in bare $^{72}\mathrm{Ge}$ ions was determined to be $23.9\left(6\right)$\,ms, which strongly deviates from expectations. △ Less

Submitted 18 December, 2023; originally announced December 2023.

Comments: Submitted to Physical Review Letters, 8 pages, 3 figures, 2 tables, 7 equations

Rational Kriging

Authors: V. Roshan Joseph

Abstract: This article proposes a new kriging that has a rational form. It is shown that the generalized least squares estimate of the mean from rational kriging is much more well behaved than that from ordinary kriging. Parameter estimation and uncertainty quantification for rational kriging are proposed using a Gaussian process framework. Its potential applications in emulation and calibration of computer… ▽ More This article proposes a new kriging that has a rational form. It is shown that the generalized least squares estimate of the mean from rational kriging is much more well behaved than that from ordinary kriging. Parameter estimation and uncertainty quantification for rational kriging are proposed using a Gaussian process framework. Its potential applications in emulation and calibration of computer models are also discussed. △ Less

Submitted 8 December, 2023; originally announced December 2023.

Updating the first CHIME/FRB catalog of fast radio bursts with baseband data

Authors: The CHIME/FRB Collaboration, :, Mandana Amiri, Bridget C. Andersen, Shion Andrew, Kevin Bandura, Mohit Bhardwaj, P. J. Boyle, Charanjot Brar, Daniela Breitman, Tomas Cassanelli, Pragya Chawla, Amanda M. Cook, Alice P. Curtin, Matt Dobbs, Fengqiu Adam Dong, Gwendolyn Eadie, Emmanuel Fonseca, B. M. Gaensler, Utkarsh Giri, Antonio Herrera-Martin, Hans Hopkins, Adaeze L. Ibik, Ronniy C. Joseph, J. F. Kaczmarek , et al. (36 additional authors not shown)

Abstract: In 2021, a catalog of 536 fast radio bursts (FRBs) detected with the Canadian Hydrogen Intensity Mapping Experiment (CHIME) radio telescope was released by the CHIME/FRB Collaboration. This large collection of bursts, observed with a single instrument and uniform selection effects, has advanced our understanding of the FRB population. Here we update the results for 140 of these FRBs for which chan… ▽ More In 2021, a catalog of 536 fast radio bursts (FRBs) detected with the Canadian Hydrogen Intensity Mapping Experiment (CHIME) radio telescope was released by the CHIME/FRB Collaboration. This large collection of bursts, observed with a single instrument and uniform selection effects, has advanced our understanding of the FRB population. Here we update the results for 140 of these FRBs for which channelized raw voltage ('baseband') data are available. With the voltages measured by the telescope's antennas, it is possible to maximize the telescope sensitivity in any direction within the primary beam, an operation called 'beamforming'. This allows us to increase the signal-to-noise ratio (S/N) of the bursts and to localize them to sub-arcminute precision. The improved localization is also used to correct the beam response of the instrument and to measure fluxes and fluences with a ~10% uncertainty. Additionally, the time resolution is increased by three orders of magnitude relative to that in the first CHIME/FRB catalog, and, applying coherent dedispersion, burst morphologies can be studied in detail. Polarization information is also available for the full sample of 140 FRBs, providing an unprecedented dataset to study the polarization properties of the population. We release the baseband data beamformed to the most probable position of each FRB. These data are analyzed in detail in a series of accompanying papers. △ Less

Submitted 22 May, 2024; v1 submitted 31 October, 2023; originally announced November 2023.

Enhancing Sample Quality through Minimum Energy Importance Weights

Authors: Chaofan Huang, V. Roshan Joseph

Abstract: Importance sampling is a powerful tool for correcting the distributional mismatch in many statistical and machine learning problems, but in practice its performance is limited by the usage of simple proposals whose importance weights can be computed analytically. To address this limitation, Liu and Lee (2017) proposed a Black-Box Importance Sampling (BBIS) algorithm that computes the importance we… ▽ More Importance sampling is a powerful tool for correcting the distributional mismatch in many statistical and machine learning problems, but in practice its performance is limited by the usage of simple proposals whose importance weights can be computed analytically. To address this limitation, Liu and Lee (2017) proposed a Black-Box Importance Sampling (BBIS) algorithm that computes the importance weights for arbitrary simulated samples by minimizing the kernelized Stein discrepancy. However, this requires knowing the score function of the target distribution, which is not easy to compute for many Bayesian problems. Hence, in this paper we propose another novel BBIS algorithm using minimum energy design, BBIS-MED, that requires only the unnormalized density function, which can be utilized as a post-processing step to improve the quality of Markov Chain Monte Carlo samples. We demonstrate the effectiveness and wide applicability of our proposed BBIS-MED algorithm on extensive simulations and a real-world Bayesian model calibration problem where the score function cannot be derived analytically. △ Less

Submitted 31 December, 2023; v1 submitted 11 October, 2023; originally announced October 2023.

Discovering the Unknowns: A First Step

Authors: V. Roshan Joseph, William E. Lewis, Henry S. Yuchi, Kathryn A. Maupin

Abstract: This article aims at discovering the unknown variables in the system through data analysis. The main idea is to use the time of data collection as a surrogate variable and try to identify the unknown variables by modeling gradual and sudden changes in the data. We use Gaussian process modeling and a sparse representation of the sudden changes to efficiently estimate the large number of parameters… ▽ More This article aims at discovering the unknown variables in the system through data analysis. The main idea is to use the time of data collection as a surrogate variable and try to identify the unknown variables by modeling gradual and sudden changes in the data. We use Gaussian process modeling and a sparse representation of the sudden changes to efficiently estimate the large number of parameters in the proposed statistical model. The method is tested on a realistic dataset generated using a one-dimensional implementation of a Magnetized Liner Inertial Fusion (MagLIF) simulation model and encouraging results are obtained. △ Less

Submitted 10 October, 2023; originally announced October 2023.

Asset Bundling for Wind Power Forecasting

Authors: Hanyu Zhang, Mathieu Tanneau, Chaofan Huang, V. Roshan Joseph, Shangkun Wang, Pascal Van Hentenryck

Abstract: The growing penetration of intermittent, renewable generation in US power grids, especially wind and solar generation, results in increased operational uncertainty. In that context, accurate forecasts are critical, especially for wind generation, which exhibits large variability and is historically harder to predict. To overcome this challenge, this work proposes a novel Bundle-Predict-Reconcile (… ▽ More The growing penetration of intermittent, renewable generation in US power grids, especially wind and solar generation, results in increased operational uncertainty. In that context, accurate forecasts are critical, especially for wind generation, which exhibits large variability and is historically harder to predict. To overcome this challenge, this work proposes a novel Bundle-Predict-Reconcile (BPR) framework that integrates asset bundling, machine learning, and forecast reconciliation techniques. The BPR framework first learns an intermediate hierarchy level (the bundles), then predicts wind power at the asset, bundle, and fleet level, and finally reconciles all forecasts to ensure consistency. This approach effectively introduces an auxiliary learning task (predicting the bundle-level time series) to help the main learning tasks. The paper also introduces new asset-bundling criteria that capture the spatio-temporal dynamics of wind power time series. Extensive numerical experiments are conducted on an industry-size dataset of 283 wind farms in the MISO footprint. The experiments consider short-term and day-ahead forecasts, and evaluates a large variety of forecasting models that include weather predictions as covariates. The results demonstrate the benefits of BPR, which consistently and significantly improves forecast accuracy over baselines, especially at the fleet level. △ Less

Submitted 28 September, 2023; originally announced September 2023.

Reinforcing Security and Usability of Crypto-Wallet with Post-Quantum Cryptography and Zero-Knowledge Proof

Authors: Yathin Kethepalli, Rony Joseph, Sai Raja Vajrala, Jashwanth Vemula, Nenavath Srinivas Naik

Abstract: Crypto-wallets or digital asset wallets are a crucial aspect of managing cryptocurrencies and other digital assets such as NFTs. However, these wallets are not immune to security threats, particularly from the growing risk of quantum computing. The use of traditional public-key cryptography systems in digital asset wallets makes them vulnerable to attacks from quantum computers, which may increase… ▽ More Crypto-wallets or digital asset wallets are a crucial aspect of managing cryptocurrencies and other digital assets such as NFTs. However, these wallets are not immune to security threats, particularly from the growing risk of quantum computing. The use of traditional public-key cryptography systems in digital asset wallets makes them vulnerable to attacks from quantum computers, which may increase in the future. Moreover, current digital wallets require users to keep track of seed-phrases, which can be challenging and lead to additional security risks. To overcome these challenges, a new algorithm is proposed that uses post-quantum cryptography (PQC) and zero-knowledge proof (ZKP) to enhance the security of digital asset wallets. The research focuses on the use of the Lattice-based Threshold Secret Sharing Scheme (LTSSS), Kyber Algorithm for key generation and ZKP for wallet unlocking, providing a more secure and user-friendly alternative to seed-phrase, brain and multi-sig protocol wallets. This algorithm also includes several innovative security features such as recovery of wallets in case of downtime of the server, and the ability to rekey the private key associated with a specific username-password combination, offering improved security and usability. The incorporation of PQC and ZKP provides a robust and comprehensive framework for securing digital assets in the present and future. This research aims to address the security challenges faced by digital asset wallets and proposes practical solutions to ensure their safety in the era of quantum computing. △ Less

Submitted 29 August, 2023; v1 submitted 14 August, 2023; originally announced August 2023.

Single-molecule fluorescence multiplexing by multi-parameter spectroscopic detection of nanostructured FRET labels

Authors: Jiachong Chu, Ayesha Ejaz, Kyle M. Lin, Madeline R. Joseph, Aria E. Coraor, D. Allan Drummond, Allison H. Squires

Abstract: Multiplexed, real-time fluorescence detection at the single-molecule level is highly desirable to reveal the stoichiometry, dynamics, and interactions of individual molecular species within complex systems. However, traditionally fluorescence sensing is limited to 3-4 concurrently detected labels, due to low signal-to-noise, high spectral overlap between labels, and the need to avoid dissimilar dy… ▽ More Multiplexed, real-time fluorescence detection at the single-molecule level is highly desirable to reveal the stoichiometry, dynamics, and interactions of individual molecular species within complex systems. However, traditionally fluorescence sensing is limited to 3-4 concurrently detected labels, due to low signal-to-noise, high spectral overlap between labels, and the need to avoid dissimilar dye chemistries. We have engineered a palette of several dozen fluorescent labels, called FRETfluors, for spectroscopic multiplexing at the single-molecule level. Each FRETfluor is a compact nanostructure formed from the same three chemical building blocks (DNA, Cy3, and Cy5). The composition and dye-dye geometries create a characteristic Förster Resonance Energy Transfer (FRET) efficiency for each construct. In addition, we varied the local DNA sequence and attachment chemistry to alter the Cy3 and Cy5 emission properties and thereby shift the emission signatures of an entire series of FRET constructs to new sectors of the multi-parameter detection space. Unique spectroscopic emission of each FRETfluor is therefore conferred by a combination of FRET and this site-specific tuning of individual fluorophore photophysics. We show single-molecule identification of a set of 27 FRETfluors in a sample mixture using a subset of constructs statistically selected to minimize classification errors, measured using an Anti-Brownian ELectrokinetic (ABEL) trap which provides precise multi-parameter spectroscopic measurements. The ABEL trap also enables discrimination between FRETfluors attached to a target (here: mRNA) and unbound FRETfluors, eliminating the need for washes or removal of excess label by purification. We show single-molecule identification of a set of 27 FRETfluors in a sample mixture using a subset of constructs selected to minimize classification errors. △ Less

Submitted 25 January, 2024; v1 submitted 4 July, 2023; originally announced July 2023.

Comments: 62 pages, 6 figures, 19 Supplementary figures, 5 Supplementary tables, 9 Supplementary notes

System Analysis Modeling and Intermodal Transportation for Commercial Spent

Authors: Harish R Gadey, Mark W Nutt, Philip Jensen, Robert Howard, Robby Joseph, Lucas Vander Wal, Riley M Cumberland

Abstract: The United States Department of Energy has long term goals to develop solutions for managing the nations spent nuclear fuel and high-level waste inventory. The Integrated Waste Management program is employing system-level engineering and analysis principles to inform potential future waste management system architectures. Managing the SNF requires the use of system-level analysis software that con… ▽ More The United States Department of Energy has long term goals to develop solutions for managing the nations spent nuclear fuel and high-level waste inventory. The Integrated Waste Management program is employing system-level engineering and analysis principles to inform potential future waste management system architectures. Managing the SNF requires the use of system-level analysis software that considers waste generation, on-site (centralized storage), transportation infrastructure, and long-term disposal. The Next Generation System Analysis Model is an agent-based model that was developed to simulate the transportation and storage of SNF and HLW. NGSAM has the capability to detail the interaction and movement of individual components and groups, such as rail cars and casks. The SNF inventory from commercial nuclear reactors is currently in temporary storage at multiple locations spread across the US. Shipping of SNF from these locations relies on one of three transportation modes: rail, heavy-haul truck, or barge. Rail is the most preferred due to the size of the canisters and casks the SNF would be shipped in. Under some scenarios, a rail route might not be available to a reactor site or improving the rail infrastructure at shutdown sites might be too cost-prohibitive for utilities to opt for a direct rail transfer. Under such scenarios, using a barge or heavy haul truck to de-inventory the site and transfer the SNF to a nearby intermodal transfer site with adequate rail infrastructure where the payload could be transferred to a rail car might prove to be an attractive option. This work presents the various intermodal transportation options to move SNF from reactor sites to rail cars. Next, the operational steps in each of these modes to move the SNF from a reactor site and transfer it to a rail car is explored. The ideology, assumptions, and future steps are presented. △ Less

Submitted 4 June, 2023; originally announced June 2023.

Comments: American Nuclear Society 2022 International High-Level Radioactive Waste Management Conference (https://www.ans.org/pubs/proceedings/article-52728/)

Storage, Accumulation and Deceleration of Secondary Beams for Nuclear Astrophysics

Authors: J. Glorius, Yu. A. Litvinov, M. Aliotta, F. Amjad, B. Brückner, C. G. Bruno, R. Chen, T. Davinson, S. F. Dellmann, T. Dickel, I. Dillmann, P. Erbacher, O. Forstner, H. Geissel, C. J. Griffin, R. Grisenti, A. Gumberidze, E. Haettner, R. Hess, P. -M. Hillenbrand, C. Hornung, R. Joseph, B. Jurado, E. Kazanseva, R. Knöbel , et al. (39 additional authors not shown)

Abstract: Low-energy investigations on rare ion beams are often limited by the available intensity and purity of the ion species in focus. Here, we present the first application of a technique that combines in-flight production at relativistic energies with subsequent secondary beam storage, accumulation and finally deceleration to the energy of interest. Using the FRS and ESR facilities at GSI, this scheme… ▽ More Low-energy investigations on rare ion beams are often limited by the available intensity and purity of the ion species in focus. Here, we present the first application of a technique that combines in-flight production at relativistic energies with subsequent secondary beam storage, accumulation and finally deceleration to the energy of interest. Using the FRS and ESR facilities at GSI, this scheme was pioneered to provide a secondary beam of $^{118}$Te$^{52+}$ for the measurement of nuclear proton-capture at energies of 6 and 7 MeV/u. The technique provided stored beam intensities of about $10^6$ ions at high purity and brilliance, representing a major step towards low-energy nuclear physics studies using rare ion beams. △ Less

Submitted 30 May, 2023; v1 submitted 25 May, 2023; originally announced May 2023.

Journal ref: Nuclear Instruments and Methods in Physics Research B 541 (2023) 190-193

A Global-Local Approximation Framework for Large-Scale Gaussian Process Modeling

Authors: Akhil Vakayil, Roshan Joseph

Abstract: In this work, we propose a novel framework for large-scale Gaussian process (GP) modeling. Contrary to the global, and local approximations proposed in the literature to address the computational bottleneck with exact GP modeling, we employ a combined global-local approach in building the approximation. Our framework uses a subset-of-data approach where the subset is a union of a set of global poi… ▽ More In this work, we propose a novel framework for large-scale Gaussian process (GP) modeling. Contrary to the global, and local approximations proposed in the literature to address the computational bottleneck with exact GP modeling, we employ a combined global-local approach in building the approximation. Our framework uses a subset-of-data approach where the subset is a union of a set of global points designed to capture the global trend in the data, and a set of local points specific to a given testing location to capture the local trend around the testing location. The correlation function is also modeled as a combination of a global, and a local kernel. The performance of our framework, which we refer to as TwinGP, is on par or better than the state-of-the-art GP modeling methods at a fraction of their computational cost. △ Less

Submitted 17 May, 2023; originally announced May 2023.

Sequential Designs for Filling Output Spaces

Authors: Shangkun Wang, Adam P. Generale, Surya R. Kalidindi, V. Roshan Joseph

Abstract: Space-filling designs are commonly used in computer experiments to fill the space of inputs so that the input-output relationship can be accurately estimated. However, in certain applications such as inverse design or feature-based modeling, the aim is to fill the response or feature space. In this article, we propose a new experimental design framework that aims to fill the space of the outputs (… ▽ More Space-filling designs are commonly used in computer experiments to fill the space of inputs so that the input-output relationship can be accurately estimated. However, in certain applications such as inverse design or feature-based modeling, the aim is to fill the response or feature space. In this article, we propose a new experimental design framework that aims to fill the space of the outputs (responses or features). The design is adaptive and model-free, and therefore is expected to be robust to different kinds of modeling choices and input-output relationships. Several examples are given to show the advantages of the proposed method over the traditional input space-filling designs. △ Less

Submitted 11 May, 2023; originally announced May 2023.

Comments: 36 pages, 12 figures

Journal ref: Technometrics (2023)

The Stochastic Toolbox User's Guide -- xSPDE3: extensible software for stochastic ordinary and partial differential equations

Authors: Simon Kiesewetter, Ria R. Joseph, Peter D. Drummond

Abstract: The xSPDE toolbox treats stochastic partial and ordinary differential equations, with applications in biology, chemistry, engineering, medicine, physics and quantum technologies. It computes statistical averages, including time-step and/or sampling error estimation. xSPDE can provide higher order convergence, Fourier spectra and probability densities. The toolbox has graphical output and $χ^{2}$ s… ▽ More The xSPDE toolbox treats stochastic partial and ordinary differential equations, with applications in biology, chemistry, engineering, medicine, physics and quantum technologies. It computes statistical averages, including time-step and/or sampling error estimation. xSPDE can provide higher order convergence, Fourier spectra and probability densities. The toolbox has graphical output and $χ^{2}$ statistics, as well as weighted, projected, or forward-backward equations. It can generate input-output quantum spectra. All equations may have independent periodic, Dirichlet, and Neumann or Robin boundary conditions in any dimension, for any vector field component, and at either end of any interval. △ Less

Submitted 27 May, 2023; v1 submitted 8 March, 2023; originally announced March 2023.

Comments: User manual for xSPDE software on Github, at https://github.com/peterddrummond/xspde_matlab. Total of 160 pages with examples. This is a replacement with minor corrections and an updated reference list. Accepted, and will appear in Scipost

Delaying transition induced by a strip of distributed roughness using additional fine grit roughness

Authors: Robin Joseph, P Phani Kumar, Sourabh S Diwan

Abstract: Distributed roughness occurs on aerodynamic surfaces like wind/gas turbine blades and aircraft wings causing early boundary layer transition resulting in a reduction of the lift-to-drag ratio/power production. Despite being a recurring theme in engineering scenarios, the mechanism of boundary layer transition caused by distributed roughness is not well understood, and consequently, methods for mit… ▽ More Distributed roughness occurs on aerodynamic surfaces like wind/gas turbine blades and aircraft wings causing early boundary layer transition resulting in a reduction of the lift-to-drag ratio/power production. Despite being a recurring theme in engineering scenarios, the mechanism of boundary layer transition caused by distributed roughness is not well understood, and consequently, methods for mitigating its effects are scarce. In this work, we present a passive method for delaying boundary layer transition caused by distributed roughness (sandpaper strip) using a combination of secondary fine roughness strips placed immediately upstream and downstream of the distributed roughness. Hot-wire and PIV measurements are used to characterize the flow features and quantify transition delay. A combination of secondary roughness strips placed both upstream and downstream is shown to be most effective in delaying transition caused by the primary distributed roughness. Results suggest that the upstream roughness lifts the boundary layer reducing the effective Reynolds number of the primary roughness, while the downstream roughness reduces the strength of vortices shed from the primary roughness. A parametric study on the length and type of secondary roughness shows that smooth strips can also delay the transition and there is likely an optimal length of the secondary roughness beyond which increasing the extent of the downstream roughness has marginal effects on transition delay. Finally, after the onset of transition, there are no specific signatures in the flow corresponding to the secondary roughness which suggests that the secondary roughness can delay transition without substantially altering the transitional flow features. The results point to an adaptable and practical method for increasing the life cycle and efficiency of aerodynamic surfaces with distributed roughness. △ Less

Submitted 20 February, 2023; originally announced February 2023.

Comments: 24 pages, 19 figures

The impact of human expert visual inspection on the discovery of strong gravitational lenses

Authors: Karina Rojas, Thomas E. Collett, Daniel Ballard, Mark R. Magee, Simon Birrer, Elizabeth Buckley-Geer., James H. H. Chan, Benjamin Clément, José M. Diego, Fabrizio Gentile, Jimena González, Rémy Joseph, Jorge Mastache, Stefan Schuldt, Crescenzo Tortora, Tomás Verdugo, Aprajita Verma, Tansu Daylan, Martin Millon, Neal Jackson, Simon Dye, Alejandra Melo, Guillaume Mahler, Ricardo L. C. Ogando, Frédéric Courbin , et al. (31 additional authors not shown)

Abstract: We investigate the ability of human 'expert' classifiers to identify strong gravitational lens candidates in Dark Energy Survey like imaging. We recruited a total of 55 people that completed more than 25$\%$ of the project. During the classification task, we present to the participants 1489 images. The sample contains a variety of data including lens simulations, real lenses, non-lens examples, an… ▽ More We investigate the ability of human 'expert' classifiers to identify strong gravitational lens candidates in Dark Energy Survey like imaging. We recruited a total of 55 people that completed more than 25$\%$ of the project. During the classification task, we present to the participants 1489 images. The sample contains a variety of data including lens simulations, real lenses, non-lens examples, and unlabeled data. We find that experts are extremely good at finding bright, well-resolved Einstein rings, whilst arcs with $g$-band signal-to-noise less than $\sim$25 or Einstein radii less than $\sim$1.2 times the seeing are rarely recovered. Very few non-lenses are scored highly. There is substantial variation in the performance of individual classifiers, but they do not appear to depend on the classifier's experience, confidence or academic position. These variations can be mitigated with a team of 6 or more independent classifiers. Our results give confidence that humans are a reliable pruning step for lens candidates, providing pure and quantifiably complete samples for follow-up studies. △ Less

Submitted 25 April, 2023; v1 submitted 9 January, 2023; originally announced January 2023.

Comments: 16 pages, 20 Figures

Adaptive Exploration and Optimization of Materials Crystal Structures

Authors: Arvind Krishna, Huan Tran, Chaofan Huang, Rampi Ramprasad, V. Roshan Joseph

Abstract: A central problem of materials science is to determine whether a hypothetical material is stable without being synthesized, which is mathematically equivalent to a global optimization problem on a highly non-linear and multi-modal potential energy surface (PES). This optimization problem poses multiple outstanding challenges, including the exceedingly high dimensionality of the PES and that PES mu… ▽ More A central problem of materials science is to determine whether a hypothetical material is stable without being synthesized, which is mathematically equivalent to a global optimization problem on a highly non-linear and multi-modal potential energy surface (PES). This optimization problem poses multiple outstanding challenges, including the exceedingly high dimensionality of the PES and that PES must be constructed from a reliable, sophisticated, parameters-free, and thus, very expensive computational method, for which density functional theory (DFT) is an example. DFT is a quantum mechanics based method that can predict, among other things, the total potential energy of a given configuration of atoms. DFT, while accurate, is computationally expensive. In this work, we propose a novel expansion-exploration-exploitation framework to find the global minimum of the PES. Starting from a few atomic configurations, this ``known'' space is expanded to construct a big candidate set. The expansion begins in a non-adaptive manner, where new configurations are added without considering their potential energy. A novel feature of this step is that it tends to generate a space-filling design without the knowledge of the boundaries of the domain space. If needed, the non-adaptive expansion of the space of configurations is followed by adaptive expansion, where ``promising regions'' of the domain space (those with low energy configurations) are further expanded. Once a candidate set of configurations is obtained, it is simultaneously explored and exploited using Bayesian optimization to find the global minimum. The methodology is demonstrated using a problem of finding the most stable crystal structure of Aluminum. △ Less

Submitted 1 December, 2022; originally announced December 2022.

Journal ref: INFORMS Journal on Data Science, 2023

Wigner's Friend paradoxes: consistency with weak-contextual and weak-macroscopic realism models

Authors: Ria Joseph, Manushan Thenabadu, Channa Hatharasinghe, Jesse Fulton, Run-Yan Teh, P. D. Drummond, M. D. Reid

Abstract: Wigner's friend paradoxes highlight contradictions between measurements made by Friends inside a laboratory and superobservers outside a laboratory, who have access to an entangled state of the measurement apparatus. The contradictions lead to no-go theorems for observer-independent facts, thus challenging concepts of objectivity. Here, we examine the paradoxes from the perspective of establishing… ▽ More Wigner's friend paradoxes highlight contradictions between measurements made by Friends inside a laboratory and superobservers outside a laboratory, who have access to an entangled state of the measurement apparatus. The contradictions lead to no-go theorems for observer-independent facts, thus challenging concepts of objectivity. Here, we examine the paradoxes from the perspective of establishing consistency with macroscopic realism. We present versions of the Brukner-Wigner-friend and Frauchiger-Renner paradoxes in which the spin-$1/2$ system measured by the Friends corresponds to two macroscopically distinct states. The local unitary operations $U_θ$ that determine the measurement setting $θ$ are carried out using nonlinear interactions, thereby ensuring measurements need only distinguish between the macroscopically distinct states. The macroscopic paradoxes are perplexing, seemingly suggesting there is no objectivity in a macroscopic limit. However, we demonstrate consistency with a contextual weak form of macroscopic realism (wMR): The premise wMR asserts that the system can be considered to have a definite spin outcome $λ_θ$, at the time after the system has undergone the unitary rotation $U_θ$ to prepare it in a suitable pointer basis. We further show that the paradoxical outcomes imply failure of deterministic macroscopic local realism, and arise when there are unitary interactions $U_θ$ occurring due to a change of measurement setting at both sites, with respect to the state prepared by each Friend. In models which validate wMR, there is a breakdown of a subset of the assumptions that constitute the Bell-Locality premise. A similar interpretation involving a weak contextual form of realism exists for the original paradoxes. △ Less

Submitted 5 November, 2022; originally announced November 2022.

Uncovering a population of gravitational lens galaxies with magnified standard candle SN Zwicky

Authors: Ariel Goobar, Joel Johansson, Steve Schulze, Nikki Arendse, Ana Sagués Carracedo, Suhail Dhawan, Edvard Mörtsell, Christoffer Fremling, Lin Yan, Daniel Perley, Jesper Sollerman, Rémy Joseph, K-Ryan Hinds, William Meynardie, Igor Andreoni, Eric Bellm, Josh Bloom, Thomas E. Collett, Andrew Drake, Matthew Graham, Mansi Kasliwal, Shri Kulkarni, Cameron Lemon, Adam A. Miller, James D. Neill , et al. (13 additional authors not shown)

Abstract: Detecting gravitationally lensed supernovae is among the biggest challenges in astronomy. It involves a combination of two very rare phenomena: catching the transient signal of a stellar explosion in a distant galaxy and observing it through a nearly perfectly aligned foreground galaxy that deflects light towards the observer. High-cadence optical observations with the Zwicky Transient Facility, w… ▽ More Detecting gravitationally lensed supernovae is among the biggest challenges in astronomy. It involves a combination of two very rare phenomena: catching the transient signal of a stellar explosion in a distant galaxy and observing it through a nearly perfectly aligned foreground galaxy that deflects light towards the observer. High-cadence optical observations with the Zwicky Transient Facility, with an unparalleled large field of view, led to the detection of a multiply-imaged Type Ia supernova (SN Ia), ``SN Zwicky", a.k.a. SN 2022qmx. Magnified nearly twenty-five times, the system was found thanks to the ``standard candle" nature of SNe Ia. High-spatial resolution imaging with the Keck telescope resolved four images of the supernova with very small angular separation, corresponding to an Einstein radius of only $θ_E =0.167"$ and almost identical arrival times. The small $θ_E$ and faintness of the lensing galaxy is very unusual, highlighting the importance of supernovae to fully characterise the properties of galaxy-scale gravitational lenses, including the impact of galaxy substructures. △ Less

Submitted 14 June, 2023; v1 submitted 1 November, 2022; originally announced November 2022.

Comments: Matches published version in Nature Astronomy

Beyond Ultra-diffuse Galaxies. I. Mass--Size Outliers among the Satellites of Milky Way Analogs

Authors: Jiaxuan Li, Jenny E. Greene, Johnny P. Greco, Song Huang, Peter Melchior, Rachael Beaton, Kirsten Casey, Shany Danieli, Andy Goulding, Remy Joseph, Erin Kado-Fong, Ji Hoon Kim, Lauren A. MacArthur

Abstract: Large diffuse galaxies are hard to find, but understanding the environments where they live, their numbers, and ultimately their origins, is of intense interest and importance for galaxy formation and evolution. Using Subaru's Hyper Suprime-Cam Strategic Survey Program, we perform a systematic search for low surface brightness galaxies and present novel and effective methods for detecting and mode… ▽ More Large diffuse galaxies are hard to find, but understanding the environments where they live, their numbers, and ultimately their origins, is of intense interest and importance for galaxy formation and evolution. Using Subaru's Hyper Suprime-Cam Strategic Survey Program, we perform a systematic search for low surface brightness galaxies and present novel and effective methods for detecting and modeling them. As a case study, we surveyed 922 Milky Way analogs in the nearby Universe ($0.01 < z < 0.04$) and build a large sample of satellite galaxies that are outliers in the mass-size relation. These ``ultra-puffy'' galaxies (UPGs), defined to be $1.5σ$ above the average mass-size relation, represent the tail of the satellite size distribution. We find that each MW analog hosts $N_{\rm UPG} = 0.31\pm 0.05$ ultra-puffy galaxies on average, which is consistent with but slightly lower than the observed abundance at this halo mass in the Local Volume. We also construct a sample of ultra-diffuse galaxies (UDGs) in MW analogs and find an abundance of $N_{\rm UDG} = 0.44\pm0.05$ per host. With literature results, we confirm that the UDG abundance scales with the host halo mass following a sublinear power law. We argue that our definition for ultra-puffy galaxies, which is based on the mass-size relation, is more physically-motivated than the common definition of ultra-diffuse galaxies, which depends on surface brightness and size cuts and thus yields different surface mass density cuts for quenched and star-forming galaxies. △ Less

Submitted 11 September, 2023; v1 submitted 26 October, 2022; originally announced October 2022.

Comments: 19 pages, 7 figures, ApJ accepted. Also see the companion paper arXiv:2302.14108

Conglomerate Multi-Fidelity Gaussian Process Modeling, with Application to Heavy-Ion Collisions

Authors: Yi Ji, Henry Shaowu Yuchi, Derek Soeder, J. -F. Paquet, Steffen A. Bass, V. Roshan Joseph, C. F. Jeff Wu, Simon Mak

Abstract: In an era where scientific experimentation is often costly, multi-fidelity emulation provides a powerful tool for predictive scientific computing. While there has been notable work on multi-fidelity modeling, existing models do not incorporate an important "conglomerate" property of multi-fidelity simulators, where the accuracies of different simulator components are controlled by different fideli… ▽ More In an era where scientific experimentation is often costly, multi-fidelity emulation provides a powerful tool for predictive scientific computing. While there has been notable work on multi-fidelity modeling, existing models do not incorporate an important "conglomerate" property of multi-fidelity simulators, where the accuracies of different simulator components are controlled by different fidelity parameters. Such conglomerate simulators are widely encountered in complex nuclear physics and astrophysics applications. We thus propose a new CONglomerate multi-FIdelity Gaussian process (CONFIG) model, which embeds this conglomerate structure within a novel non-stationary covariance function. We show that the proposed CONFIG model can capture prior knowledge on the numerical convergence of conglomerate simulators, which allows for cost-efficient emulation of multi-fidelity systems. We demonstrate the improved predictive performance of CONFIG over state-of-the-art models in a suite of numerical experiments and two applications, the first for emulation of cantilever beam deflection and the second for emulating the evolution of the quark-gluon plasma, which was theorized to have filled the Universe shortly after the Big Bang. △ Less

Submitted 28 September, 2023; v1 submitted 27 September, 2022; originally announced September 2022.

Risk-Aware Control and Optimization for High-Renewable Power Grids

Authors: Neil Barry, Minas Chatzos, Wenbo Chen, Dahye Han, Chaofan Huang, Roshan Joseph, Michael Klamkin, Seonho Park, Mathieu Tanneau, Pascal Van Hentenryck, Shangkun Wang, Hanyu Zhang, Haoruo Zhao

Abstract: The transition of the electrical power grid from fossil fuels to renewable sources of energy raises fundamental challenges to the market-clearing algorithms that drive its operations. Indeed, the increased stochasticity in load and the volatility of renewable energy sources have led to significant increases in prediction errors, affecting the reliability and efficiency of existing deterministic op… ▽ More The transition of the electrical power grid from fossil fuels to renewable sources of energy raises fundamental challenges to the market-clearing algorithms that drive its operations. Indeed, the increased stochasticity in load and the volatility of renewable energy sources have led to significant increases in prediction errors, affecting the reliability and efficiency of existing deterministic optimization models. The RAMC project was initiated to investigate how to move from this deterministic setting into a risk-aware framework where uncertainty is quantified explicitly and incorporated in the market-clearing optimizations. Risk-aware market-clearing raises challenges on its own, primarily from a computational standpoint. This paper reviews how RAMC approaches risk-aware market clearing and presents some of its innovations in uncertainty quantification, optimization, and machine learning. Experimental results on real networks are presented. △ Less

Submitted 2 April, 2022; originally announced April 2022.

Optimal Ratio for Data Splitting

Authors: V. Roshan Joseph

Abstract: It is common to split a dataset into training and testing sets before fitting a statistical or machine learning model. However, there is no clear guidance on how much data should be used for training and testing. In this article we show that the optimal splitting ratio is $\sqrt{p}:1$, where $p$ is the number of parameters in a linear regression model that explains the data well. It is common to split a dataset into training and testing sets before fitting a statistical or machine learning model. However, there is no clear guidance on how much data should be used for training and testing. In this article we show that the optimal splitting ratio is $\sqrt{p}:1$, where $p$ is the number of parameters in a linear regression model that explains the data well. △ Less

Submitted 7 February, 2022; originally announced February 2022.

Journal ref: Statistical Analysis and Data Mining: The ASA Data Science Journal, 2022

Rubin-Euclid Derived Data Products: Initial Recommendations

Authors: Leanne P. Guy, Jean-Charles Cuillandre, Etienne Bachelet, Manda Banerji, Franz E. Bauer, Thomas Collett, Christopher J. Conselice, Siegfried Eggl, Annette Ferguson, Adriano Fontana, Catherine Heymans, Isobel M. Hook, Éric Aubourg, Hervé Aussel, James Bosch, Benoit Carry, Henk Hoekstra, Konrad Kuijken, Francois Lanusse, Peter Melchior, Joseph Mohr, Michele Moresco, Reiko Nakajima, Stéphane Paltani, Michael Troxel , et al. (95 additional authors not shown)

Abstract: This report is the result of a joint discussion between the Rubin and Euclid scientific communities. The work presented in this report was focused on designing and recommending an initial set of Derived Data products (DDPs) that could realize the science goals enabled by joint processing. All interested Rubin and Euclid data rights holders were invited to contribute via an online discussion forum… ▽ More This report is the result of a joint discussion between the Rubin and Euclid scientific communities. The work presented in this report was focused on designing and recommending an initial set of Derived Data products (DDPs) that could realize the science goals enabled by joint processing. All interested Rubin and Euclid data rights holders were invited to contribute via an online discussion forum and a series of virtual meetings. Strong interest in enhancing science with joint DDPs emerged from across a wide range of astrophysical domains: Solar System, the Galaxy, the Local Volume, from the nearby to the primaeval Universe, and cosmology. △ Less

Submitted 13 October, 2022; v1 submitted 11 January, 2022; originally announced January 2022.

Comments: Report of the Rubin-Euclid Derived Data Products Working Group, 78 pages, 11 figures

A Survey of Deep Learning Techniques for Dynamic Branch Prediction

Authors: Rinu Joseph

Abstract: Branch prediction is an architectural feature that speeds up the execution of branch instruction on pipeline processors and reduces the cost of branching. Recent advancements of Deep Learning (DL) in the post Moore's Law era is accelerating areas of automated chip design, low-power computer architectures, and much more. Traditional computer architecture design and algorithms could benefit from dyn… ▽ More Branch prediction is an architectural feature that speeds up the execution of branch instruction on pipeline processors and reduces the cost of branching. Recent advancements of Deep Learning (DL) in the post Moore's Law era is accelerating areas of automated chip design, low-power computer architectures, and much more. Traditional computer architecture design and algorithms could benefit from dynamic predictors based on deep learning algorithms which learns from experience by optimizing its parameters on large number of data. In this survey paper, we focus on traditional branch prediction algorithms, analyzes its limitations, and presents a literature survey of how deep learning techniques can be applied to create dynamic branch predictors capable of predicting conditional branch instructions. Prior surveys in this field focus on dynamic branch prediction techniques based on neural network perceptrons. We plan to improve the survey based on latest research in DL and advanced Machine Learning (ML) based branch predictors. △ Less

Submitted 29 December, 2021; originally announced December 2021.

Comments: Survey paper

A midpoint projection algorithm for stochastic differential equations on manifolds

Authors: Ria Rushin Joseph, Jesse van Rhijn, Peter D. Drummond

Abstract: Stochastic differential equations projected onto manifolds occur in physics, chemistry, biology, engineering, nanotechnology and optimization, with interdisciplinary applications. Intrinsic coordinate stochastic equations on the manifold are often computationally impractical, and numerical projections are useful in many cases. We show that the Stratonovich interpretation of the stochastic calculus… ▽ More Stochastic differential equations projected onto manifolds occur in physics, chemistry, biology, engineering, nanotechnology and optimization, with interdisciplinary applications. Intrinsic coordinate stochastic equations on the manifold are often computationally impractical, and numerical projections are useful in many cases. We show that the Stratonovich interpretation of the stochastic calculus is obtained using adiabatic elimination with a constraint potential. We derive intrinsic stochastic equations for spheroidal and hyperboloidal surfaces for comparison purposes, and review some earlier projection algorithms. In this paper, a combined midpoint projection algorithm is proposed that uses a midpoint projection onto a tangent space, combined with a subsequent normal projection to satisfy the constraints. Numerical examples are given for a range of manifolds, including circular, spheroidal, hyperboloidal, and catenoidal cases, as well as higher-order polynomial constraints and a ten-dimensional hypersphere. We show that in all cases the combined midpoint method has greatly reduced errors compared to methods using a combined Euler projection approach or purely tangential projection. Our technique can handle multiple constraints. This allows manifolds that embody several conserved quantities. The algorithm is accurate, simple and efficient. An order of magnitude error reduction in diffusion distance is typically found compared to the other methods, with reductions of several orders of magnitude in constraint errors. △ Less

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

Comments: Revised version with error analysis and more examples

Validating Synthetic Galaxy Catalogs for Dark Energy Science in the LSST Era

Authors: Eve Kovacs, Yao-Yuan Mao, Michel Aguena, Anita Bahmanyar, Adam Broussard, James Butler, Duncan Campbell, Chihway Chang, Shenming Fu, Katrin Heitmann, Danila Korytov, François Lanusse, Patricia Larsen, Rachel Mandelbaum, Christopher B. Morrison, Constantin Payerne, Marina Ricci, Eli Rykoff, F. Javier Sánchez, Ignacio Sevilla-Noarbe, Melanie Simet, Chun-Hao To, Vinu Vikraman, Rongpu Zhou, Camille Avestruz , et al. (14 additional authors not shown)

Abstract: Large simulation efforts are required to provide synthetic galaxy catalogs for ongoing and upcoming cosmology surveys. These extragalactic catalogs are being used for many diverse purposes covering a wide range of scientific topics. In order to be useful, they must offer realistically complex information about the galaxies they contain. Hence, it is critical to implement a rigorous validation proc… ▽ More Large simulation efforts are required to provide synthetic galaxy catalogs for ongoing and upcoming cosmology surveys. These extragalactic catalogs are being used for many diverse purposes covering a wide range of scientific topics. In order to be useful, they must offer realistically complex information about the galaxies they contain. Hence, it is critical to implement a rigorous validation procedure that ensures that the simulated galaxy properties faithfully capture observations and delivers an assessment of the level of realism attained by the catalog. We present here a suite of validation tests that have been developed by the Rubin Observatory Legacy Survey of Space and Time (LSST) Dark Energy Science Collaboration (DESC). We discuss how the inclusion of each test is driven by the scientific targets for static ground-based dark energy science and by the availability of suitable validation data. The validation criteria that are used to assess the performance of a catalog are flexible and depend on the science goals. We illustrate the utility of this suite by showing examples for the validation of cosmoDC2, the extragalactic catalog recently released for the LSST DESC second Data Challenge. △ Less

Submitted 13 January, 2022; v1 submitted 7 October, 2021; originally announced October 2021.

Comments: 46 pages, 33 figures

Epoch of Reionization Power Spectrum Limits from Murchison Widefield Array Data Targeted at EoR1 Field

Authors: M. Rahimi, B. Pindor, J. L. B. Line, N. Barry, C. M. Trott, R. L. Webster, C. H. Jordan, M. Wilensky, S. Yoshiura, A. Beardsley, J. Bowman, R. Byrne, A. Chokshi, B. J. Hazelton, K. Hasegawa, E. Howard, B. Greig, D. Jacobs, R. Joseph, M. Kolopanis, C. Lynch, B. McKinley, D. A. Mitchell, S. Murray, M. F. Morales , et al. (6 additional authors not shown)

Abstract: Current attempts to measure the 21cm Power Spectrum of neutral hydrogen during the Epoch of Reionization are limited by systematics which produce measured upper limits above both the thermal noise and the expected cosmological signal. These systematics arise from a combination of observational, instrumental, and analysis effects. In order to further understand and mitigate these effects, it is ins… ▽ More Current attempts to measure the 21cm Power Spectrum of neutral hydrogen during the Epoch of Reionization are limited by systematics which produce measured upper limits above both the thermal noise and the expected cosmological signal. These systematics arise from a combination of observational, instrumental, and analysis effects. In order to further understand and mitigate these effects, it is instructive to explore different aspects of existing datasets. One such aspect is the choice of observing field. To date, MWA EoR observations have largely focused on the EoR0 field. In this work, we present a new detailed analysis of the EoR1 field. The EoR1 field is one of the coldest regions of the Southern radio sky, but contains the very bright radio galaxy Fornax-A. The presence of this bright extended source in the primary beam of the interferometer makes the calibration and analysis of EoR1 particularly challenging. We demonstrate the effectiveness of a recently developed shapelet model of Fornax-A in improving the results from this field. We also describe and apply a series of data quality metrics which identify and remove systematically contaminated data. With substantially improved source models, upgraded analysis algorithms and enhanced data quality metrics, we determine EoR power spectrum upper limits based on analysis of the best $\sim$14-hours data observed during 2015 and 2014 at redshifts 6.5, 6.8 and 7.1, with the lowest $2σ$ upper limit at z=6.5 of $Δ^2 \leq (73.78 ~\mathrm{mK)^2}$ at $k=0.13~\mathrm{h~ Mpc^{-1}}$, improving on previous EoR1 measurement results. △ Less

Submitted 7 October, 2021; originally announced October 2021.

Comments: 18 pages, 17 figures, Accepted for publication in MNRAS

Data Twinning

Authors: Akhil Vakayil, V. Roshan Joseph

Abstract: In this work, we develop a method named Twinning, for partitioning a dataset into statistically similar twin sets. Twinning is based on SPlit, a recently proposed model-independent method for optimally splitting a dataset into training and testing sets. Twinning is orders of magnitude faster than the SPlit algorithm, which makes it applicable to Big Data problems such as data compression. Twinning… ▽ More In this work, we develop a method named Twinning, for partitioning a dataset into statistically similar twin sets. Twinning is based on SPlit, a recently proposed model-independent method for optimally splitting a dataset into training and testing sets. Twinning is orders of magnitude faster than the SPlit algorithm, which makes it applicable to Big Data problems such as data compression. Twinning can also be used for generating multiple splits of a given dataset to aid divide-and-conquer procedures and $k$-fold cross validation. △ Less

Submitted 6 October, 2021; originally announced October 2021.

Strong lens systems search in the Dark Energy Survey using Convolutional Neural Networks

Authors: K. Rojas, E. Savary, B. Clément, M. Maus, F. Courbin, C. Lemon, J. H. H. Chan, G. Vernardos, R. Joseph, R. Cañameras, A. Galan

Abstract: We performed a search for strong lens galaxy-scale systems in the first data release of the Dark Energy Survey (DES), from a color-selected parent sample of 18~745~029 Luminous Red Galaxies (LRGs). Our search was based on a Convolutional Neural Network (CNN) to grade our LRG selection with values between 0 (non-lens) and 1 (lens). Our training set was data-driven, i.e. using lensed sources taken f… ▽ More We performed a search for strong lens galaxy-scale systems in the first data release of the Dark Energy Survey (DES), from a color-selected parent sample of 18~745~029 Luminous Red Galaxies (LRGs). Our search was based on a Convolutional Neural Network (CNN) to grade our LRG selection with values between 0 (non-lens) and 1 (lens). Our training set was data-driven, i.e. using lensed sources taken from HST COSMOS images and where the light distribution of the lens plane was taken directly from DES images of our LRGs. A total of 76~582 cutouts obtained a score above 0.9. These were visually inspected and resulted in two catalogs. The first one contains 405 lens candidates, where 90 present clear lensing features and counterparts, while the others 315 require more evidence, such as higher resolution images or spectra to be conclusive. A total of 186 candidates were totally new identified in this search. The second catalog includes 539 ring galaxy candidates that will be useful to train CNNs against this type of false positives. For the 90 best lens candidates we carried out color-based deblending of the lens and source light without fitting any analytical profile to the data. The method turned out to be very efficient in the deblending, even for very compact objects and for objects with very complex morphology. Finally, from the 90 best lens candidates we selected 52 systems having one single deflector, to test an automated modeling pipeline which successfully modeled 79\% of the sample within an acceptable amount of computing time. △ Less

Submitted 31 August, 2021; originally announced September 2021.

Comments: 37 pages, 31 figures, submitted to A&A

Journal ref: A&A 668, A73 (2022)

Joint survey processing: combined resampling and convolution for galaxy modelling and deblending

Authors: Rémy Joseph, Peter Melchior, Fred Moolekamp

Abstract: We present an extension of the multi-band galaxy fitting method scarlet which allows the joint modeling of astronomical images from different instruments, by performing simultaneous resampling and convolution. We introduce a fast and formally accurate linear projection operation that maps a pixelated model at a given resolution onto an observation frame with a different point spread function and p… ▽ More We present an extension of the multi-band galaxy fitting method scarlet which allows the joint modeling of astronomical images from different instruments, by performing simultaneous resampling and convolution. We introduce a fast and formally accurate linear projection operation that maps a pixelated model at a given resolution onto an observation frame with a different point spread function and pixel scale. We test our implementation against the well-tested resampling and convolution method in galsim on simulated images mimicking observations with the Euclid space telescope and the Vera C. Rubin Observatory, and find that it reduces interpolation errors by an order of magnitude or more compared to galsim default settings. Tests with a wide range of levels of blending show more accurate galaxy models from joint modeling of Euclid and Rubin images compared to separate modeling of each survey by up to an order of magnitude. Our results demonstrate, for the first time, the feasibility and utility of performing non-parametric pixel-level data fusion of overlapping imaging surveys. All results can be reproduced with the specific versions of the codes and notebooks used in this study. △ Less

Submitted 14 July, 2021; originally announced July 2021.

Comments: Submitted to Astronomy and Computing

A new MWA limit on the 21 cm Power Spectrum at Redshifts $\sim$ 13 $-$ 17

Authors: S. Yoshiura, B. Pindor, J. L. B. Line, N. Barry, C. M. Trott, A. Beardsley, J. Bowman, R. Byrne, A. Chokshi, B. J. Hazelton, K. Hasegawa, E. Howard, B. Greig, D. Jacobs, C. H. Jordan, R. Joseph, M. Kolopanis, C. Lynch, B. McKinley, D. A. Mitchell, M. F. Morales, S. G. Murray, J. C. Pober, M. Rahimi, K. Takahashi , et al. (7 additional authors not shown)

Abstract: Observations in the lowest MWA band between $75-100$ MHz have the potential to constrain the distribution of neutral hydrogen in the intergalactic medium at redshift $\sim 13-17$. Using 15 hours of MWA data, we analyse systematics in this band such as radio-frequency interference (RFI), ionospheric and wide field effects. By updating the position of point sources, we mitigate the direction indepen… ▽ More Observations in the lowest MWA band between $75-100$ MHz have the potential to constrain the distribution of neutral hydrogen in the intergalactic medium at redshift $\sim 13-17$. Using 15 hours of MWA data, we analyse systematics in this band such as radio-frequency interference (RFI), ionospheric and wide field effects. By updating the position of point sources, we mitigate the direction independent calibration error due to ionospheric offsets. Our calibration strategy is optimized for the lowest frequency bands by reducing the number of direction dependent calibrators and taking into account radio sources within a wider field of view. We remove data polluted by systematics based on the RFI occupancy and ionospheric conditions, finally selecting 5.5 hours of the cleanest data. Using these data, we obtain two sigma upper limits on the 21 cm power spectrum in the range of $0.1\lessapprox k \lessapprox 1 ~\rm ~h~Mpc^{-1}$ and at $z$=14.2, 15.2 and 16.5, with the lowest limit being $6.3\times 10^6 ~\rm mK^2$ at $\rm k=0.14 \rm ~h~Mpc^{-1}$ and at $z=15.2$ with a possibility of a few \% of signal loss due to direction independent calibration. △ Less

Submitted 26 May, 2021; originally announced May 2021.

Comments: 17 pages, 18 figures, accepted to MNRAS

Constrained Minimum Energy Designs

Authors: Chaofan Huang, V. Roshan Joseph, Douglas M. Ray

Abstract: Space-filling designs are important in computer experiments, which are critical for building a cheap surrogate model that adequately approximates an expensive computer code. Many design construction techniques in the existing literature are only applicable for rectangular bounded space, but in real world applications, the input space can often be non-rectangular because of constraints on the input… ▽ More Space-filling designs are important in computer experiments, which are critical for building a cheap surrogate model that adequately approximates an expensive computer code. Many design construction techniques in the existing literature are only applicable for rectangular bounded space, but in real world applications, the input space can often be non-rectangular because of constraints on the input variables. One solution to generate designs in a constrained space is to first generate uniformly distributed samples in the feasible region, and then use them as the candidate set to construct the designs. Sequentially Constrained Monte Carlo (SCMC) is the state-of-the-art technique for candidate generation, but it still requires large number of constraint evaluations, which is problematic especially when the constraints are expensive to evaluate. Thus, to reduce constraint evaluations and improve efficiency, we propose the Constrained Minimum Energy Design (CoMinED) that utilizes recent advances in deterministic sampling methods. Extensive simulation results on 15 benchmark problems with dimensions ranging from 2 to 13 are provided for demonstrating the improved performance of CoMinED over the existing methods. △ Less

Submitted 24 April, 2021; originally announced April 2021.

Comments: Submitted to Statistics and Computing

Journal ref: Stat Comput 31, 80 (2021)

Interacting fermion dynamics in Majorana phase-space

Authors: Ria Rushin Joseph, Laura E C Rosales-Zárate, Peter D Drummond

Abstract: The problem of fermion dynamics is studied using the Q-function for fermions. This is a probabilistic phase-space representation, which we express using Majorana operators, so that the phase-space variable is a real antisymmetric matrix. We consider a general interaction Hamiltonian with four Majorana operators and arbitrary properties. Our model includes the Majorana Hubbard and Fermi Hubbard Ham… ▽ More The problem of fermion dynamics is studied using the Q-function for fermions. This is a probabilistic phase-space representation, which we express using Majorana operators, so that the phase-space variable is a real antisymmetric matrix. We consider a general interaction Hamiltonian with four Majorana operators and arbitrary properties. Our model includes the Majorana Hubbard and Fermi Hubbard Hamiltonians, as well as general quantum field theories of interacting fermions. Using the Majorana Q-function we derive a generalized Fokker-Planck equation, with results for the drift and diffusion terms. The diffusion term is proved to be traceless, which gives a dynamical interpretation as a forwards-backwards stochastic process. This approach leads to a model of quantum measurement in terms of an ontology with real vacuum fluctuations. △ Less

Submitted 24 April, 2021; originally announced April 2021.

Comments: 12 pages

Population Quasi-Monte Carlo

Authors: Chaofan Huang, V. Roshan Joseph, Simon Mak

Abstract: Monte Carlo methods are widely used for approximating complicated, multidimensional integrals for Bayesian inference. Population Monte Carlo (PMC) is an important class of Monte Carlo methods, which utilizes a population of proposals to generate weighted samples that approximate the target distribution. The generic PMC framework iterates over three steps: samples are simulated from a set of propos… ▽ More Monte Carlo methods are widely used for approximating complicated, multidimensional integrals for Bayesian inference. Population Monte Carlo (PMC) is an important class of Monte Carlo methods, which utilizes a population of proposals to generate weighted samples that approximate the target distribution. The generic PMC framework iterates over three steps: samples are simulated from a set of proposals, weights are assigned to such samples to correct for mismatch between the proposal and target distributions, and the proposals are then adapted via resampling from the weighted samples. When the target distribution is expensive to evaluate, the PMC has its computational limitation since the convergence rate is $\mathcal{O}(N^{-1/2})$. To address this, we propose in this paper a new Population Quasi-Monte Carlo (PQMC) framework, which integrates Quasi-Monte Carlo ideas within the sampling and adaptation steps of PMC. A key novelty in PQMC is the idea of importance support points resampling, a deterministic method for finding an "optimal" subsample from the weighted proposal samples. Moreover, within the PQMC framework, we develop an efficient covariance adaptation strategy for multivariate normal proposals. Lastly, a new set of correction weights is introduced for the weighted PMC estimator to improve the efficiency from the standard PMC estimator. We demonstrate the improved empirical convergence of PQMC over PMC in extensive numerical simulations and a friction drilling application. △ Less

Submitted 26 December, 2020; originally announced December 2020.

Comments: Submitted to Journal of Computational and Graphical Statistics

Journal ref: Journal of Computational and Graphical Statistics (2022)

SPlit: An Optimal Method for Data Splitting

Authors: V. Roshan Joseph, Akhil Vakayil

Abstract: In this article we propose an optimal method referred to as SPlit for splitting a dataset into training and testing sets. SPlit is based on the method of Support Points (SP), which was initially developed for finding the optimal representative points of a continuous distribution. We adapt SP for subsampling from a dataset using a sequential nearest neighbor algorithm. We also extend SP to deal wit… ▽ More In this article we propose an optimal method referred to as SPlit for splitting a dataset into training and testing sets. SPlit is based on the method of Support Points (SP), which was initially developed for finding the optimal representative points of a continuous distribution. We adapt SP for subsampling from a dataset using a sequential nearest neighbor algorithm. We also extend SP to deal with categorical variables so that SPlit can be applied to both regression and classification problems. The implementation of SPlit on real datasets shows substantial improvement in the worst-case testing performance for several modeling methods compared to the commonly used random splitting procedure. △ Less

Submitted 19 March, 2021; v1 submitted 20 December, 2020; originally announced December 2020.

SLITronomy: towards a fully wavelet-based strong lensing inversion technique

Authors: A. Galan, A. Peel, R. Joseph, F. Courbin, J. -L. Starck

Abstract: Strong gravitational lensing provides a wealth of astrophysical information on the baryonic and dark matter content of galaxies. It also serves as a valuable cosmological probe by allowing us to measure the Hubble constant independently of other methods. These applications all require the difficult task of inverting the lens equation and simultaneously reconstructing the mass profile of the lens a… ▽ More Strong gravitational lensing provides a wealth of astrophysical information on the baryonic and dark matter content of galaxies. It also serves as a valuable cosmological probe by allowing us to measure the Hubble constant independently of other methods. These applications all require the difficult task of inverting the lens equation and simultaneously reconstructing the mass profile of the lens along with the original light profile of the unlensed source. As there is no reason for either the lens or the source to be simple, we need methods that both invert the lens equation with a large number of degrees of freedom and also enforce a well-controlled regularisation that avoids the appearance of spurious structures. This can be beautifully accomplished by representing signals in wavelet space. Building on the Sparse Lens Inversion Technique (SLIT), in this work we present an improved sparsity-based method that describes lensed sources using wavelets and optimises over the parameters given an analytical lens mass profile. We apply our technique on simulated HST and E-ELT data, as well as on real HST images of lenses from the Sloan Lens ACS (SLACS) sample, assuming a lens model. We show that wavelets allow us to reconstruct lensed sources containing detailed substructures when using both present-day data and high-resolution images from future thirty-meter-class telescopes. Wavelets moreover provide a much more tractable solution in terms of quality and computation time compared to using a source model that combines smooth analytical profiles and shapelets. Requiring very little human interaction, our pixel-based technique fits into the effort to devise automated modelling schemes. It can be incorporated in the standard workflow of sampling analytical lens model parameters. The method, which we call SLITronomy, is freely available as a new plug-in to the modelling software Lenstronomy. △ Less

Submitted 4 December, 2020; originally announced December 2020.

Journal ref: A&A 647, A176 (2021)

The Impact of Tandem Redundant/Sky-Based Calibration in MWA Phase II Data Analysis

Authors: Zheng Zhang, Jonathan C. Pober, Wenyang Li, Bryna J. Hazelton, Miguel F. Morales, Cathryn M. Trott, Christopher H. Jordan, Ronniy C. Joseph, Adam Beardsley, Nichole Barry, Ruby Byrne, Steven J. Tingay, Aman Chokshi, Kenji Hasegawa, Daniel C. Jacobs, Adam Lanman, Jack L. B. Line, Christene Lynch, Benjamin McKinley, Daniel A. Mitchell, Steven Murray, Bart Pindor, Mahsa Rahimi, Keitaro Takahashi, Randall B. Wayth , et al. (4 additional authors not shown)

Abstract: Precise instrumental calibration is of crucial importance to 21-cm cosmology experiments. The Murchison Widefield Array's (MWA) Phase II compact configuration offers us opportunities for both redundant calibration and sky-based calibration algorithms; using the two in tandem is a potential approach to mitigate calibration errors caused by inaccurate sky models. The MWA Epoch of Reionization (EoR)… ▽ More Precise instrumental calibration is of crucial importance to 21-cm cosmology experiments. The Murchison Widefield Array's (MWA) Phase II compact configuration offers us opportunities for both redundant calibration and sky-based calibration algorithms; using the two in tandem is a potential approach to mitigate calibration errors caused by inaccurate sky models. The MWA Epoch of Reionization (EoR) experiment targets three patches of the sky (dubbed EoR0, EoR1, and EoR2) with deep observations. Previous work in \cite{Li_2018} and \cite{Wenyang_2019} studied the effect of tandem calibration on the EoR0 field and found that it yielded no significant improvement in the power spectrum over sky-based calibration alone. In this work, we apply similar techniques to the EoR1 field and find a distinct result: the improvements in the power spectrum from tandem calibration are significant. To understand this result, we analyze both the calibration solutions themselves and the effects on the power spectrum over three nights of EoR1 observations. We conclude that the presence of the bright radio galaxy Fornax A in EoR1 degrades the performance of sky-based calibration, which in turn enables redundant calibration to have a larger impact. These results suggest that redundant calibration can indeed mitigate some level of model-incompleteness error. △ Less

Submitted 19 September, 2020; originally announced September 2020.

Comments: 20 pages, 11 figures, accepted for publication in PASA

Robust Experimental Designs for Model Calibration

Authors: Arvind Krishna, V. Roshan Joseph, Shan Ba, William A. Brenneman, William R. Myers

Abstract: A computer model can be used for predicting an output only after specifying the values of some unknown physical constants known as calibration parameters. The unknown calibration parameters can be estimated from real data by conducting physical experiments. This paper presents an approach to optimally design such a physical experiment. The problem of optimally designing physical experiment, using… ▽ More A computer model can be used for predicting an output only after specifying the values of some unknown physical constants known as calibration parameters. The unknown calibration parameters can be estimated from real data by conducting physical experiments. This paper presents an approach to optimally design such a physical experiment. The problem of optimally designing physical experiment, using a computer model, is similar to the problem of finding optimal design for fitting nonlinear models. However, the problem is more challenging than the existing work on nonlinear optimal design because of the possibility of model discrepancy, that is, the computer model may not be an accurate representation of the true underlying model. Therefore, we propose an optimal design approach that is robust to potential model discrepancies. We show that our designs are better than the commonly used physical experimental designs that do not make use of the information contained in the computer model and other nonlinear optimal designs that ignore potential model discrepancies. We illustrate our approach using a toy example and a real example from industry. △ Less

Submitted 2 August, 2020; originally announced August 2020.

Comments: 25 pages, 10 figures

The Impact of Realistic Foreground and Instrument Models on 21cm Epoch of Reionization Experiments

Authors: Ainulnabilah Nasirudin, Steven Murray, Cathryn Trott, Bradley Greig, Ronniy Joseph, Chris Power

Abstract: Predictions for the ability of 21-cm interferometric experiments to discriminate Epoch of Reionization (EoR) signal models are typically limited by the simplicity of data models, whereby foreground signals and characteristics of the instrument are often simplified or neglected.To move towards more realistic scenarios, we explore the effects of applying more realistic foreground and instrument mode… ▽ More Predictions for the ability of 21-cm interferometric experiments to discriminate Epoch of Reionization (EoR) signal models are typically limited by the simplicity of data models, whereby foreground signals and characteristics of the instrument are often simplified or neglected.To move towards more realistic scenarios, we explore the effects of applying more realistic foreground and instrument models to the 21cm signal, and the ability to estimate astrophysical parameters with these additional complexities. We use a highly-optimized version of \textsc{21cmFAST}, integrated into \textsc{21cmMC}, to generate lightcones of the brightness temperature fluctuation for Bayesian parameter estimation. We include a statistical point-source foreground model and an instrument model based on the Murchison Widefield Array (MWA) scaled in observation time to have an effective sensitivity similar to the future Square Kilometre Array (SKA). We also extend the current likelihood prescription to account for the presence of beam convolution and foregrounds, the 2-Dimensional Power Spectrum (PS), and the correlation of PS modes. We use frequency bands between 150 and 180 MHz to constrain the ionizing efficiency ($ζ$), the minimum virial temperature of halos ($T_{\mathrm{vir}}$), the soft X-ray emissivity per unit Star Formation Rate (SFR) ($L_X/SFR$ ), and the X-ray energy threshold ($E_0$). We find that the inclusion of realistic foregrounds and instrumental components biases the parameter constraints due to unaccounted for cross-power between the EoR signal, foregrounds and thermal noise. This causes estimates of $ζ$ to be biased by up to $5σ$ but the estimates of $T_{vir}$, L$_X$/SFR and E$_0$ remain unaffected and are all within $1σ$. △ Less

Submitted 18 March, 2020; originally announced March 2020.

COSMOGRAIL XIX: Time delays in 18 strongly lensed quasars from 15 years of optical monitoring

Authors: M. Millon, F. Courbin, V. Bonvin, E. Paic, G. Meylan, M. Tewes, D. Sluse, P. Magain, J. H. H. Chan, A. Galan, R. Joseph, C. Lemon, O. Tihhonova, R. I. Anderson, M. Marmier, B. Chazelas, M. Lendl, A. H. M. J. Triaud, A. Wyttenbach

Abstract: We present the results of 15 years of monitoring lensed quasars, which was conducted by the COSMOGRAIL programme at the Leonhard Euler 1.2m Swiss Telescope. The decade-long light curves of 23 lensed systems are presented for the first time. We complement our data set with other monitoring data available in the literature to measure the time delays in 18 systems, among which nine reach a relative p… ▽ More We present the results of 15 years of monitoring lensed quasars, which was conducted by the COSMOGRAIL programme at the Leonhard Euler 1.2m Swiss Telescope. The decade-long light curves of 23 lensed systems are presented for the first time. We complement our data set with other monitoring data available in the literature to measure the time delays in 18 systems, among which nine reach a relative precision better than 15% for at least one time delay. To achieve this, we developed an automated version of the curve-shifting toolbox PyCS to ensure robust estimation of the time delay in the presence of microlensing, while accounting for the errors due to the imperfect representation of microlensing. We also re-analysed the previously published time delays of RX J1131$-$1231 and HE 0435$-$1223, by adding six and two new seasons of monitoring, respectively, and confirming the previous time-delay measurements. When the time delay measurement is possible, we corrected the light curves of the lensed images from their time delay and present the difference curves to highlight the microlensing signal contained in the data. To date, this is the largest sample of decade-long lens monitoring data, which is useful to measure $H_0$ and the size of quasar accretion discs with microlensing as well as to study quasar variability. △ Less

Submitted 25 August, 2020; v1 submitted 13 February, 2020; originally announced February 2020.

Comments: Published in A&A, 31 pages, 8 figures, 3 Tables

Journal ref: A&A 640, A105 (2020)

Deep multi-redshift limits on Epoch of Reionisation 21cm Power Spectra from Four Seasons of Murchison Widefield Array Observations

Authors: Cathryn M. Trott, C. H. Jordan, S. Midgley, N. Barry, B. Greig, B. Pindor, J. H. Cook, G. Sleap, S. J. Tingay, D. Ung, P. Hancock, A. Williams, J. Bowman, R. Byrne, A. Chokshi, B. J. Hazelton, K. Hasegawa, D. Jacobs, R. C. Joseph, W. Li, J. L. B Line, C. Lynch, B. McKinley, D. A. Mitchell, M. F. Morales , et al. (11 additional authors not shown)

Abstract: We compute the spherically-averaged power spectrum from four seasons of data obtained for the Epoch of Reionisation (EoR) project observed with the Murchison Widefield Array (MWA). We measure the EoR power spectrum over $k= 0.07-3.0~h$Mpc$^{-1}$ at redshifts $z=6.5-8.7$. The largest aggregation of 110 hours on EoR0 high-band (3,340 observations), yields a lowest measurement of (43~mK)$^2$ = 1.8… ▽ More We compute the spherically-averaged power spectrum from four seasons of data obtained for the Epoch of Reionisation (EoR) project observed with the Murchison Widefield Array (MWA). We measure the EoR power spectrum over $k= 0.07-3.0~h$Mpc$^{-1}$ at redshifts $z=6.5-8.7$. The largest aggregation of 110 hours on EoR0 high-band (3,340 observations), yields a lowest measurement of (43~mK)$^2$ = 1.8$\times$10$^3$ mK$^2$ at $k$=0.14~$h$Mpc$^{-1}$ and $z=6.5$ (2$σ$ thermal noise plus sample variance). Using the Real-Time System to calibrate and the CHIPS pipeline to estimate power spectra, we select the best observations from the central five pointings within the 2013--2016 observing seasons, observing three independent fields and in two frequency bands. This yields 13,591 2-minute snapshots (453 hours), based on a quality assurance metric that measures ionospheric activity. We perform another cut to remove poorly-calibrated data, based on power in the foreground-dominated and EoR-dominated regions of the two-dimensional power spectrum, reducing the set to 12,569 observations (419 hours). These data are processed in groups of 20 observations, to retain the capacity to identify poor data, and used to analyse the evolution and structure of the data over field, frequency, and data quality. We subsequently choose the cleanest 8,935 observations (298 hours of data) to form integrated power spectra over the different fields, pointings and redshift ranges. △ Less

Submitted 6 February, 2020; originally announced February 2020.

Comments: 19 pages, 29 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

Calibration and 21-cm Power Spectrum Estimation in the Presence of Antenna Beam Variations

Authors: Ronniy C. Joseph, C. M. Trott, R. B. Wayth, A. Nasirudin

Abstract: Detecting a signal from the Epoch of Reionisation (EoR) requires an exquisite understanding of galactic and extra-galactic foregrounds, low frequency radio instruments, instrumental calibration, and data analysis pipelines. In this work we build upon existing work that aims to understand the impact of calibration errors on 21-cm power spectrum (PS) measurements. It is well established that calibra… ▽ More Detecting a signal from the Epoch of Reionisation (EoR) requires an exquisite understanding of galactic and extra-galactic foregrounds, low frequency radio instruments, instrumental calibration, and data analysis pipelines. In this work we build upon existing work that aims to understand the impact of calibration errors on 21-cm power spectrum (PS) measurements. It is well established that calibration errors have the potential to inhibit EoR detections by introducing additional spectral features that mimic the structure of EoR signals. We present a straightforward way to estimate the impact of a wide variety of modelling residuals in EoR PS estimation. We apply this framework to the specific case of broken dipoles in Murchison Widefield Array (MWA) to understand its effect and estimate its impact on PS estimation. Combining an estimate of the percentage of MWA tiles that have at least one broken dipole (15%-40%) with an analytic description of beam errors induced by such dipoles, we compute the residuals of the foregrounds after calibration and source subtraction. We find that that incorrect beam modelling introduces bias in the 2D-PS on the order of $\sim 10^3\, \mathrm{mK}^2 \,h^{-3}\, \mathrm{Mpc}^{3}$. Although this is three orders of magnitude lower than current lowest limits, it is two orders of magnitude higher than the expected signal. Determining the accuracy of both current beam models and direction dependent calibration pipelines is therefore crucial in our search for an EoR signal. △ Less

Submitted 25 February, 2020; v1 submitted 29 November, 2019; originally announced November 2019.

Comments: 14 pages, 8 figures. Accepted for publication in MNRAS December 2019. Erratum submitted on 20 February 2020

First Season MWA Phase II EoR Power Spectrum Results at Redshift 7

Authors: W. Li, J. C. Pober, N. Barry, B. J. Hazelton, M. F. Morales, C. M. Trott, A. Lanman, M. Wilensky, I. Sullivan, A. P. Beardsley, T. Booler, J. D. Bowman, R. Byrne, B. Crosse, D. Emrich, T. M. O. Franzen, K. Hasegawa, L. Horsley, M. Johnston-Hollitt, D. C. Jacobs, C. H. Jordan, R. C. Joseph, T. Kaneuji, D. L. Kaplan, D. Kenney , et al. (22 additional authors not shown)

Abstract: The compact configuration of Phase II of the Murchison Widefield Array (MWA) consists of both a redundant subarray and pseudo-random baselines, offering unique opportunities to perform sky-model and redundant interferometric calibration. The highly redundant hexagonal cores give improved power spectrum sensitivity. In this paper, we present the analysis of nearly 40 hours of data targeting one of… ▽ More The compact configuration of Phase II of the Murchison Widefield Array (MWA) consists of both a redundant subarray and pseudo-random baselines, offering unique opportunities to perform sky-model and redundant interferometric calibration. The highly redundant hexagonal cores give improved power spectrum sensitivity. In this paper, we present the analysis of nearly 40 hours of data targeting one of the MWA's EoR fields observed in 2016. We use both improved analysis techniques presented in Barry et al. (2019) as well as several additional techniques developed for this work, including data quality control methods and interferometric calibration approaches. We show the EoR power spectrum limits at redshift 6.5, 6.8 and 7.1 based on our deep analysis on this 40-hour data set. These limits span a range in $k$ space of $0.18$ $h$ $\mathrm{Mpc^{-1}}$ $<k<1.6$ $h$ $\mathrm{Mpc^{-1}}$, with a lowest measurement of $Δ^2\leqslant2.39\times 10^3$ $\mathrm{mK}^2$ at $k=0.59$ $h$ $\mathrm{Mpc^{-1}}$ and $z=6.5$. △ Less

Submitted 20 December, 2019; v1 submitted 22 November, 2019; originally announced November 2019.

Comments: 26 pages, 15 figures, 3 tables. Updated acknowledgements. Published in ApJ

Journal ref: ApJ (2019) 887 141

Distributional Clustering: A distribution-preserving clustering method

Authors: Arvind Krishna, Simon Mak, Roshan Joseph

Abstract: One key use of k-means clustering is to identify cluster prototypes which can serve as representative points for a dataset. However, a drawback of using k-means cluster centers as representative points is that such points distort the distribution of the underlying data. This can be highly disadvantageous in problems where the representative points are subsequently used to gain insights on the data… ▽ More One key use of k-means clustering is to identify cluster prototypes which can serve as representative points for a dataset. However, a drawback of using k-means cluster centers as representative points is that such points distort the distribution of the underlying data. This can be highly disadvantageous in problems where the representative points are subsequently used to gain insights on the data distribution, as these points do not mimic the distribution of the data. To this end, we propose a new clustering method called "distributional clustering", which ensures cluster centers capture the distribution of the underlying data. We first prove the asymptotic convergence of the proposed cluster centers to the data generating distribution, then present an efficient algorithm for computing these cluster centers in practice. Finally, we demonstrate the effectiveness of distributional clustering on synthetic and real datasets. △ Less

Submitted 14 November, 2019; originally announced November 2019.

Comments: Submitted to Statistica Sinica

Proximal Adam: Robust Adaptive Update Scheme for Constrained Optimization

Authors: Peter Melchior, Rémy Joseph, Fred Moolekamp

Abstract: We implement the adaptive step size scheme from the optimization methods AdaGrad and Adam in a novel variant of the Proximal Gradient Method (PGM). Our algorithm, dubbed AdaProx, avoids the need for explicit computation of the Lipschitz constants or additional line searches and thus reduces per-iteration cost. In test cases for Constrained Matrix Factorization we demonstrate the advantages of AdaP… ▽ More We implement the adaptive step size scheme from the optimization methods AdaGrad and Adam in a novel variant of the Proximal Gradient Method (PGM). Our algorithm, dubbed AdaProx, avoids the need for explicit computation of the Lipschitz constants or additional line searches and thus reduces per-iteration cost. In test cases for Constrained Matrix Factorization we demonstrate the advantages of AdaProx in fidelity and performance over PGM, while still allowing for arbitrary penalty functions. The python implementation of the algorithm presented here is available as an open-source package at https://github.com/pmelchior/proxmin. △ Less

Submitted 2 July, 2020; v1 submitted 22 October, 2019; originally announced October 2019.

Comments: 7 pages, 5 figures