Showing 1–50 of 119 results for author: Rose, S

Back from the dead: AT2019aalc as a candidate repeating TDE in an AGN

Authors: Patrik Milán Veres, Anna Franckowiak, Sjoert van Velzen, Bjoern Adebahr, Sam Taziaux, Jannis Necker, Robert Stein, Alexander Kier, Ancla Mueller, Dominik J. Bomans, Nuria Jordana-Mitjans, Marek Kowalski, Erica Hammerstein, Elena Marci-Boehncke, Simeon Reusch, Simone Garrappa, Sam Rose, Kaustav Kashyap Das

Abstract: Context. To date, three nuclear transients have been associated with high-energy neutrino events. These transients are generally thought to be powered by tidal disruptions of stars (TDEs) by massive black holes. However, AT2019aalc, hosted in a Seyfert-1 galaxy, was not yet classified due to a lack of multiwavelength observations. Interestingly, the source has re-brightened 4 years after its disco… ▽ More Context. To date, three nuclear transients have been associated with high-energy neutrino events. These transients are generally thought to be powered by tidal disruptions of stars (TDEs) by massive black holes. However, AT2019aalc, hosted in a Seyfert-1 galaxy, was not yet classified due to a lack of multiwavelength observations. Interestingly, the source has re-brightened 4 years after its discovery. Aims. We aim to classify the transient and explain the mechanism responsible for its second optical flare. Methods. We conducted a multi-wavelength monitoring program (from radio to X-rays) of AT2019aalc during its re-brightening in 2023. Results. The observations revealed a uniquely bright UV counterpart and multiple X-ray flares during the second optical flaring episode of the transient. The second flare, similarly to the first one, is also accompanied by IR dust echo emission. A long-term radio flare is found with an inverted spectrum. Optical spectroscopic observations reveal the presence of Bowen Fluorescence lines and strong high-ionization coronal lines indicating an extreme level of ionization in the system. Conclusions. The results suggest that the transient can be classified as a Bowen Fluorescence Flare (BFF), a relatively new sub-class of flaring active galactic nuclei (AGN). AT2019aalc can be also classified as an extreme coronal line emitter (ECLE). We found that, in addition to AT2019aalc, another BFF AT2021loi is spatially coincident with a high-energy neutrino event. The multi-wavelength properties of these transients suggest a possible connection between ECLEs, BFFs and TDEs in AGN. △ Less

Submitted 30 August, 2024; originally announced August 2024.

Comments: Submitted to A&A. 27 pages, 24 figures, 2 tables

Investigating the Electron Capture Supernova Candidate AT 2019abn with JWST Spectroscopy

Authors: Sam Rose, Ryan M. Lau, Jacob E. Jencson, Mansi M. Kasliwal, Kishalay De, Michael E. Ressler, Ori D. Fox, Matthew J. Hankins

Abstract: The James Webb Space Telescope (JWST) has opened up a new window to study highly reddened explosive transients. We present results from late-time (1421 days post-explosion) JWST follow-up spectroscopic observations with NIRSpec and MIRI LRS of the intermediate luminosity red transient (ILRT) AT 2019abn located in the nearby Messier 51 galaxy (8.6 Mpc). ILRTs represent a mysterious class of transie… ▽ More The James Webb Space Telescope (JWST) has opened up a new window to study highly reddened explosive transients. We present results from late-time (1421 days post-explosion) JWST follow-up spectroscopic observations with NIRSpec and MIRI LRS of the intermediate luminosity red transient (ILRT) AT 2019abn located in the nearby Messier 51 galaxy (8.6 Mpc). ILRTs represent a mysterious class of transients which exhibit peak luminosities between those of classical novae and supernovae and which are known to be highly dust obscured. Similar to the prototypical examples of this class of objects, NGC 300 2008-OT and SN 2008S, AT 2019abn has an extremely red and dusty progenitor detected only in pre-explosion Spitzer/IRAC imaging at 3.6 and 4.5 micron and not in deep optical or near-infrared HST images. We find that late time observations of AT 2019abn from NEOWISE and JWST are consistent with the late time evolution of SN 2008S. In part because they are so obscured by dust, it is unknown what produces an ILRT with hypotheses ranging from high mass stellar merger events, non-terminal stellar outbursts, or terminal supernovae explosions through electron-capture in super-AGB stars. Our JWST observations show strong mid-IR Class C PAH features at 6.3 and 8.25 micron typical of carbon-rich post-AGB sources. These features suggest the dust around AT 2019abn, either pre-existing or newly formed in the ejecta, is composed of carbonaceous grains which are not typically observed around red supergiants. However, depending on the strength and temperature of hot bottom burning, SAGBs may be expected to exhibit a carbon-rich chemistry. Thus our JWST observations are consistent with AT 2019abn having an SAGB progenitor. △ Less

Submitted 29 July, 2024; originally announced July 2024.

Comments: 12 pages, 4 figures, submitted to ApJL

Bounding elastic photon-photon scattering at $\sqrt s \approx 1\,$MeV using a laser-plasma platform

Authors: R. Watt, B. Kettle, E. Gerstmayr, B. King, A. Alejo, S. Astbury, C. Baird, S. Bohlen, M. Campbell, C. Colgan, D. Dannheim, C. Gregory, H. Harsh, P. Hatfield, J. Hinojosa, D. Hollatz, Y. Katzir, J. Morton, C. D. Murphy, A. Nurnberg, J. Osterhoff, G. Pérez-Callejo, K. Põder, P. P. Rajeev, C. Roedel , et al. (14 additional authors not shown)

Abstract: We report on a direct search for elastic photon-photon scattering using x-ray and $γ$ photons from a laser-plasma based experiment. A gamma photon beam produced by a laser wakefield accelerator provided a broadband gamma spectrum extending to above $E_γ= 200$ MeV. These were collided with a dense x-ray field produced by the emission from a laser heated germanium foil at $E_x \approx 1.4$ keV, corr… ▽ More We report on a direct search for elastic photon-photon scattering using x-ray and $γ$ photons from a laser-plasma based experiment. A gamma photon beam produced by a laser wakefield accelerator provided a broadband gamma spectrum extending to above $E_γ= 200$ MeV. These were collided with a dense x-ray field produced by the emission from a laser heated germanium foil at $E_x \approx 1.4$ keV, corresponding to an invariant mass of $\sqrt{s} = 1.22 \pm 0.22$ MeV. In these asymmetric collisions elastic scattering removes one x-ray and one high-energy $γ$ photon and outputs two lower energy $γ$ photons. No changes in the $γ$ photon spectrum were observed as a result of the collisions allowing us to place a 95% upper bound on the cross section of $1.5 \times 10^{15}\,μ$b. Although far from the QED prediction, this represents the lowest upper limit obtained so far for $\sqrt{s} \lesssim 1$ MeV. △ Less

Submitted 17 July, 2024; originally announced July 2024.

Comments: 6 pages, 10 figures

Algorithmic Changes Are Not Enough: Evaluating the Removal of Race Adjustment from the eGFR Equation

Authors: Marika M. Cusick, Glenn M. Chertow, Douglas K. Owens, Michelle Y. Williams, Sherri Rose

Abstract: Changing clinical algorithms to remove race adjustment has been proposed and implemented for multiple health conditions. Removing race adjustment from estimated glomerular filtration rate (eGFR) equations may reduce disparities in chronic kidney disease (CKD), but has not been studied in clinical practice after implementation. Here, we assessed whether implementing an eGFR equation (CKD-EPI 2021)… ▽ More Changing clinical algorithms to remove race adjustment has been proposed and implemented for multiple health conditions. Removing race adjustment from estimated glomerular filtration rate (eGFR) equations may reduce disparities in chronic kidney disease (CKD), but has not been studied in clinical practice after implementation. Here, we assessed whether implementing an eGFR equation (CKD-EPI 2021) without adjustment for Black or African American race modified quarterly rates of nephrology referrals and visits within a single healthcare system, Stanford Health Care (SHC). Our cohort study analyzed 547,194 adult patients aged 21 and older who had at least one recorded serum creatinine or serum cystatin C between January 1, 2019 and September 1, 2023. During the study period, implementation of CKD-EPI 2021 did not modify rates of quarterly nephrology referrals in those documented as Black or African American or in the overall cohort. After adjusting for capacity at SHC nephrology clinics, estimated rates of nephrology referrals and visits with CKD-EPI 2021 were 34 (95% CI 29, 39) and 188 (175, 201) per 10,000 patients documented as Black or African American. If race adjustment had not been removed, estimated rates were nearly identical: 38 (95% CI: 28, 53) and 189 (165, 218) per 10,000 patients. Changes to the eGFR equation are likely insufficient to achieve health equity in CKD care decision-making as many other structural inequities remain. △ Less

Submitted 25 April, 2024; v1 submitted 19 April, 2024; originally announced April 2024.

Comments: Accepted to Conference on Health, Inference, and Learning (CHIL) 2024

Classifying Graphemes in English Words Through the Application of a Fuzzy Inference System

Authors: Samuel Rose, Chandrasekhar Kambhampati

Abstract: In Linguistics, a grapheme is a written unit of a writing system corresponding to a phonological sound. In Natural Language Processing tasks, written language is analysed through two different mediums, word analysis, and character analysis. This paper focuses on a third approach, the analysis of graphemes. Graphemes have advantages over word and character analysis by being self-contained represent… ▽ More In Linguistics, a grapheme is a written unit of a writing system corresponding to a phonological sound. In Natural Language Processing tasks, written language is analysed through two different mediums, word analysis, and character analysis. This paper focuses on a third approach, the analysis of graphemes. Graphemes have advantages over word and character analysis by being self-contained representations of phonetic sounds. Due to the nature of splitting a word into graphemes being based on complex, non-binary rules, the application of fuzzy logic would provide a suitable medium upon which to predict the number of graphemes in a word. This paper proposes the application of a Fuzzy Inference System to split words into their graphemes. This Fuzzy Inference System results in a correct prediction of the number of graphemes in a word 50.18% of the time, with 93.51% being within a margin of +- 1 from the correct classification. Given the variety in language, graphemes are tied with pronunciation and therefore can change depending on a regional accent/dialect, the +- 1 accuracy represents the impreciseness of grapheme classification when regional variances are accounted for. To give a baseline of comparison, a second method involving a recursive IPA mapping exercise using a pronunciation dictionary was developed to allow for comparisons to be made. △ Less

Submitted 2 April, 2024; originally announced April 2024.

SN 2023zaw: an ultra-stripped, nickel-poor supernova from a low-mass progenitor

Authors: Kaustav K. Das, Christoffer Fremling, Mansi M. Kasliwal, Steve Schulze, Jesper Sollerman, Viraj Karambelkar, Sam Rose, Shreya Anand, Igor Andreoni, Marie Aubert, Sean J. Brennan, S. Bradley Cenko, Michael W. Coughlin, B. O'Connor, Kishalay De, Jim Fuller, Matthew Graham, Erica Hammerstein, Annastasia Haynie, K-Ryan Hinds, Io Kleiser, S. R. Kulkarni, Zeren Lin, Chang Liu, Ashish A. Mahabal , et al. (12 additional authors not shown)

Abstract: We present SN 2023zaw $-$ a sub-luminous ($\mathrm{M_r} = -16.7$ mag) and rapidly-evolving supernova ($\mathrm{t_{1/2,r}} = 4.9$ days), with the lowest nickel mass ($\approx0.002$ $\mathrm{M_\odot}$) measured among all stripped-envelope supernovae discovered to date. The photospheric spectra are dominated by broad He I and Ca NIR emission lines with velocities of $\sim10\ 000 - 12\ 000$… ▽ More We present SN 2023zaw $-$ a sub-luminous ($\mathrm{M_r} = -16.7$ mag) and rapidly-evolving supernova ($\mathrm{t_{1/2,r}} = 4.9$ days), with the lowest nickel mass ($\approx0.002$ $\mathrm{M_\odot}$) measured among all stripped-envelope supernovae discovered to date. The photospheric spectra are dominated by broad He I and Ca NIR emission lines with velocities of $\sim10\ 000 - 12\ 000$ $\mathrm{km\ s^{-1}}$. The late-time spectra show prominent narrow He I emission lines at $\sim$1000$\ \mathrm{km\ s^{-1}}$, indicative of interaction with He-rich circumstellar material. SN 2023zaw is located in the spiral arm of a star-forming galaxy. We perform radiation-hydrodynamical and analytical modeling of the lightcurve by fitting with a combination of shock-cooling emission and nickel decay. The progenitor has a best-fit envelope mass of $\approx0.2$ $\mathrm{M_\odot}$ and an envelope radius of $\approx50$ $\mathrm{R_\odot}$. The extremely low nickel mass and low ejecta mass ($\approx0.5$ $\mathrm{M_\odot}$) suggest an ultra-stripped SN, which originates from a mass-losing low mass He-star (ZAMS mass $<$ 10 $\mathrm{M_\odot}$) in a close binary system. This is a channel to form double neutron star systems, whose merger is detectable with LIGO. SN 2023zaw underscores the existence of a previously undiscovered population of extremely low nickel mass ($< 0.005$ $\mathrm{M_\odot}$) stripped-envelope supernovae, which can be explored with deep and high-cadence transient surveys. △ Less

Submitted 7 August, 2024; v1 submitted 12 March, 2024; originally announced March 2024.

Comments: Accepted for publication in ApJL

Journal ref: The Astrophysical Journal Letters, July 2024, Volume 969, Issue 1, id.L11, 18 pp

Dramatic rebrightening of the type-changing stripped-envelope supernova SN 2023aew

Authors: Yashvi Sharma, Jesper Sollerman, Shrinivas R. Kulkarni, Takashi J. Moriya, Steve Schulze, Stan Barmentloo, Michael Fausnaugh, Avishay Gal-Yam, Anders Jerkstrand, Tomás Ahumada, Eric C. Bellm, Kaustav K. Das, Andrew Drake, Christoffer Fremling, Saarah Hall, K. R. Hinds, Theophile Jegou du Laz, Viraj Karambelkar, Mansi M. Kasliwal, Frank J. Masci, Adam A. Miller, Guy Nir, Daniel A. Perley, Josiah N. Purdum, Yu-Jing Qin , et al. (10 additional authors not shown)

Abstract: Multi-peaked supernovae with precursors, dramatic light-curve rebrightenings, and spectral transformation are rare, but are being discovered in increasing numbers by modern night-sky transient surveys like the Zwicky Transient Facility (ZTF). Here, we present the observations and analysis of SN 2023aew, which showed a dramatic increase in brightness following an initial luminous (-17.4 mag) and lo… ▽ More Multi-peaked supernovae with precursors, dramatic light-curve rebrightenings, and spectral transformation are rare, but are being discovered in increasing numbers by modern night-sky transient surveys like the Zwicky Transient Facility (ZTF). Here, we present the observations and analysis of SN 2023aew, which showed a dramatic increase in brightness following an initial luminous (-17.4 mag) and long (~100 days) unusual first peak (possibly precursor). SN 2023aew was classified as a Type IIb supernova during the first peak but changed its type to resemble a stripped-envelope supernova (SESN) after the marked rebrightening. We present comparisons of SN 2023aew's spectral evolution with SESN subtypes and argue that it is similar to SNe Ibc during its main peak. P-Cygni Balmer lines are present during the first peak, but vanish during the second peak's photospheric phase, before H$α$ resurfaces again during the nebular phase. The nebular lines ([O I], [Ca II], Mg I], H$α$) exhibit a double-peaked structure which hints towards a clumpy or non-spherical ejecta. We analyze the second peak in the light curve of SN 2023aew and find it to be broader than normal SESNe as well as requiring a very high $^{56}$Ni mass to power the peak luminosity. We discuss the possible origins of SN 2023aew including an eruption scenario where a part of the envelope is ejected during the first peak which also powers the second peak of the light curve through SN-CSM interaction. △ Less

Submitted 5 February, 2024; originally announced February 2024.

Comments: 22 pages, 11 figures, 5 tables

Notes on Bolyai's 'Appendix'

Authors: Steven Rose

Abstract: This paper aims to provide an explanatory edition of Bolyai's 'Appendix Demonstrating the Absolute Science of Space', first published in 1832. In this treatise Bolyai began by extending neutral (or 'absolute') geometry by deriving a number of theorems which are independent of Euclid's parallel postulate. Then, while retaining Euclid first four postulates, he explored the consequences of replacing… ▽ More This paper aims to provide an explanatory edition of Bolyai's 'Appendix Demonstrating the Absolute Science of Space', first published in 1832. In this treatise Bolyai began by extending neutral (or 'absolute') geometry by deriving a number of theorems which are independent of Euclid's parallel postulate. Then, while retaining Euclid first four postulates, he explored the consequences of replacing the parallel postulate with a different assumption, that through a given point not on a given straight line, more than than one straight line can be drawn which does not intersect the given line. On this basis Bolyai developed a non-Euclidean geometry nowadays called hyperbolic geometry. The English translation of Bolyai's text is reprinted in the paper, with explanatory notes provided at the end of each section. Some of the theorems are discussed in detail, with a particular focus on those in which Bolyai derived the trigonometric identities of the hyperbolic plane and on his astonishing quadrature of the circle. Bolyai's terse style of exposition can be quite challenging, so I have tried to fill in some of the mathematical details which he chose to omit, either through lack of space or simply because he thought them to obvious too include. △ Less

Submitted 22 July, 2023; originally announced December 2023.

Comments: 96 pages, 49 figures

MSC Class: 01A55

Collisional Shaping of Nuclear Star Cluster Density Profiles

Authors: Sanaea C. Rose, Morgan MacLeod

Abstract: A supermassive black hole (SMBH) surrounded by a dense, nuclear star cluster resides at the center of many galaxies. In this dense environment, high-velocity collisions frequently occur between stars. About $10 \%$ of the stars within the Milky Way's nuclear star cluster collide with other stars before evolving off the main-sequence. Collisions preferentially affect tightly-bound stars, which orbi… ▽ More A supermassive black hole (SMBH) surrounded by a dense, nuclear star cluster resides at the center of many galaxies. In this dense environment, high-velocity collisions frequently occur between stars. About $10 \%$ of the stars within the Milky Way's nuclear star cluster collide with other stars before evolving off the main-sequence. Collisions preferentially affect tightly-bound stars, which orbit most quickly and pass through regions of the highest stellar density. Over time, collisions therefore shape the bulk properties of the nuclear star cluster. We examine the effect of collisions on the cluster's stellar density profile. We show that collisions produce a turning point in the density profile which can be determined analytically. Varying the initial density profile and collision model, we characterize the evolution of the stellar density profile over $10$ Gyr. We find that old, initially cuspy populations exhibit a break around $0.1$ pc in their density profile, while shallow density profiles retain their initial shape outside of $0.01$ pc. The initial density profile is always preserved outside of a few tenths of parsec irrespective of initial conditions. Lastly, we comment on the implications of collisions for the luminosity and color of stars in the collisionly-shaped inner cluster. △ Less

Submitted 30 October, 2023; originally announced October 2023.

Comments: Submitted to ApJ Letters. Comments welcome

Spike-based Neuromorphic Computing for Next-Generation Computer Vision

Authors: Md Sakib Hasan, Catherine D. Schuman, Zhongyang Zhang, Tauhidur Rahman, Garrett S. Rose

Abstract: Neuromorphic Computing promises orders of magnitude improvement in energy efficiency compared to traditional von Neumann computing paradigm. The goal is to develop an adaptive, fault-tolerant, low-footprint, fast, low-energy intelligent system by learning and emulating brain functionality which can be realized through innovation in different abstraction layers including material, device, circuit,… ▽ More Neuromorphic Computing promises orders of magnitude improvement in energy efficiency compared to traditional von Neumann computing paradigm. The goal is to develop an adaptive, fault-tolerant, low-footprint, fast, low-energy intelligent system by learning and emulating brain functionality which can be realized through innovation in different abstraction layers including material, device, circuit, architecture and algorithm. As the energy consumption in complex vision tasks keep increasing exponentially due to larger data set and resource-constrained edge devices become increasingly ubiquitous, spike-based neuromorphic computing approaches can be viable alternative to deep convolutional neural network that is dominating the vision field today. In this book chapter, we introduce neuromorphic computing, outline a few representative examples from different layers of the design stack (devices, circuits and algorithms) and conclude with a few exciting applications and future research directions that seem promising for computer vision in the near future. △ Less

Submitted 16 March, 2024; v1 submitted 14 October, 2023; originally announced October 2023.

Comments: Pending to be published as a book chapter in the book 'Computer Vision: Challenges, Trends, and Opportunities' from CRC Press

Harnessing Unipolar Threshold Switches for Enhanced Rectification

Authors: Md Mazharul Islam, Shamiul Alam, Garrett S. Rose, Aly Fathy, Sumeet Kumar Gupta, Ahmedullah Aziz

Abstract: Phase transition materials (PTM) have drawn significant attention in recent years due to their abrupt threshold switching characteristics and hysteretic behavior. Augmentation of the PTM with a transistor has been shown to provide enhanced selectivity (as high as ~107 for Ag/HfO2/Pt) leading to unique circuit-level advantages. Previously, a unipolar PTM, Ag-HfO2-Pt, was reported as a replacement f… ▽ More Phase transition materials (PTM) have drawn significant attention in recent years due to their abrupt threshold switching characteristics and hysteretic behavior. Augmentation of the PTM with a transistor has been shown to provide enhanced selectivity (as high as ~107 for Ag/HfO2/Pt) leading to unique circuit-level advantages. Previously, a unipolar PTM, Ag-HfO2-Pt, was reported as a replacement for diodes due to its polaritydependent high selectivity and hysteretic properties. It was shown to achieve ~50% higher DC output compared to a diode-based design in a Cockcroft-Walton multiplier circuit. In this paper, we take a deeper dive into this design. We augment two different PTMs (unipolar Ag-HfO2-Pt and bipolar VO2) with diodeconnected MOSFETs to retain the benefits of hysteretic rectification. Our proposed hysteretic diodes (Hyperdiodes) exhibit a low forward voltage drop owing to their volatile hysteretic characteristics. However, augmenting a hysteretic PTM with a transistor brings an additional stability concern due to their complex interplay. Hence, we perform a comprehensive stability analysis for a range of threshold voltages (-0.2 < Vth < 0.8) and transistor sizes to ensure operational stability and to choose the most optimum design parameters. We then test a standalone AgHfO2-Pt and an Ag-HfO2-Pt-based Hyperdiode in two different types of voltage multipliers and report ~500 and ~20 times lower settling time, respectively. △ Less

Submitted 30 August, 2023; originally announced October 2023.

Dielectronic satellite emission from a solid-density Mg plasma: relationship to models of ionisation potential depression

Authors: G. Pérez-Callejo, T. Gawne, T. R. Preston, P. Hollebon, O. S. Humphries, H. -K. Chung, G. L. Dakovski, J. Krzywinski, M. P. Minitti, T. Burian, J. Chalupský, V. Hájková, L. Juha, V. Vozda, U. Zastrau, S. M. Vinko, S. J. Rose, J. S. Wark

Abstract: We report on experiments where solid-density Mg plasmas are created by heating with the focused output of the Linac Coherent Light Source x-ray free-electron-laser. We study the K-shell emission from the Helium and Lithium-like ions using Bragg crystal spectroscopy. Observation of the dielectronic satellites in Lithium-like ions confirms that the M-shell electrons appear bound for these high charg… ▽ More We report on experiments where solid-density Mg plasmas are created by heating with the focused output of the Linac Coherent Light Source x-ray free-electron-laser. We study the K-shell emission from the Helium and Lithium-like ions using Bragg crystal spectroscopy. Observation of the dielectronic satellites in Lithium-like ions confirms that the M-shell electrons appear bound for these high charge states. An analysis of the intensity of these satellites indicates that when modelled with an atomic-kinetics code, the ionisation potential depression model employed needs to produce depressions for these ions which lie between those predicted by the well known Stewart-Pyatt and Ecker-Kroll models. These results are largely consistent with recent Density Functional Theory calculations. △ Less

Submitted 4 March, 2024; v1 submitted 5 October, 2023; originally announced October 2023.

Generation of photoionized plasmas in the laboratory of relevance to accretion-powered x-ray sources using keV line radiation

Authors: D. Riley, R. L. Singh, S White, M. Charlwood, D. Bailie, C. Hyland, T. Audet, G. Sarri, B. Kettle, G. Gribakin, S. J. Rose, E. G. Hill, G. J. Ferland, R. J. R. Williams, F. P. Keenan

Abstract: We describe laboratory experiments to generate X-ray photoionized plasmas of relevance to accretion-powered X-ray sources such as neutron star binaries and quasars, with significant improvements over previous work. A key quantity is referenced, namely the photoionization parameter. This is normally meaningful in an astrophysical steady-state context, but is also commonly used in the literature as… ▽ More We describe laboratory experiments to generate X-ray photoionized plasmas of relevance to accretion-powered X-ray sources such as neutron star binaries and quasars, with significant improvements over previous work. A key quantity is referenced, namely the photoionization parameter. This is normally meaningful in an astrophysical steady-state context, but is also commonly used in the literature as a figure of merit for laboratory experiments that are, of necessity, time-dependent. We demonstrate emission-weighted values of ξ > 50 ergcm/s using laser-plasma X-ray sources, with higher results at the centre of the plasma which are in the regime of interest for several astrophysical scenarios. Comparisons of laboratory experiments with astrophysical codes are always limited, principally by the many orders of magnitude differences in time and spatial scales, but also other plasma parameters. However useful checks on performance can often be made for a limited range of parameters. For example, we show that our use of a keV line source, rather than the quasi-blackbody radiation fields normally employed in such experiments, has allowed the generation of the ratio of inner-shell to outer-shell photoionization expected from a blackbody source with ~keV spectral temperature. We compare calculations from our in-house plasma modelling code with those from Cloudy and find moderately good agreement for the time evolution of both electron temperature and average ionisation. However, a comparison of code predictions for a K-beta argon X-ray spectrum with experimental data reveals that our Cloudy simulation overestimates the intensities of more highly ionised argon species. This is not totally surprising as the Cloudy model was generated for a single set of plasma conditions, while the experimental data are spatially integrated. △ Less

Submitted 19 March, 2024; v1 submitted 13 September, 2023; originally announced September 2023.

Comments: High Energy Density Physics, in press; 21 pages, 10 figures

Reimagining Sense Amplifiers: Harnessing Phase Transition Materials for Current and Voltage Sensing

Authors: Md Mazharul Islam, Shamiul Alam, Mohammad Adnan Jahangir, Garrett S. Rose, Suman Datta, Vijaykrishnan Narayanan, Sumeet Kumar Gupta, Ahmedullah Aziz

Abstract: Energy-efficient sense amplifier (SA) circuits are essential for reliable detection of stored memory states in emerging memory systems. In this work, we present four novel sense amplifier (SA) topologies based on phase transition material (PTM) tailored for non-volatile memory applications. We utilize the abrupt switching and volatile hysteretic characteristics of PTMs which enables efficient and… ▽ More Energy-efficient sense amplifier (SA) circuits are essential for reliable detection of stored memory states in emerging memory systems. In this work, we present four novel sense amplifier (SA) topologies based on phase transition material (PTM) tailored for non-volatile memory applications. We utilize the abrupt switching and volatile hysteretic characteristics of PTMs which enables efficient and fast sensing operation in our proposed SA topologies. We provide comprehensive details of their functionality and assess how process variations impact their performance metrics. Our proposed sense amplifier topologies manifest notable performance enhancement. We achieve a ~67% reduction in sensing delay and a ~80% decrease in sensing power for current sensing. For voltage sensing, we achieve a ~75% reduction in sensing delay and a ~33% decrease in sensing power. Moreover, the proposed SA topologies exhibit improved variation robustness compared to conventional SAs. We also scrutinize the dependence of transistor mirroring window and PTM transition voltages on several device parameters to determine the optimum operating conditions and stance of tunability for each of the proposed SA topologies. △ Less

Submitted 30 August, 2023; originally announced August 2023.

Functional Specification of the RAVENS Neuroprocessor

Authors: Adam Z. Foshie, James S. Plank, Garrett S. Rose, Catherine D. Schuman

Abstract: RAVENS is a neuroprocessor that has been developed by the TENNLab research group at the University of Tennessee. Its main focus has been as a vehicle for chip design with memristive elements; however it has also been the vehicle for all-digital CMOS development, plus it has implementations on FPGA's, microcontrollers and software simulation. The software simulation is supported by the TENNLab neur… ▽ More RAVENS is a neuroprocessor that has been developed by the TENNLab research group at the University of Tennessee. Its main focus has been as a vehicle for chip design with memristive elements; however it has also been the vehicle for all-digital CMOS development, plus it has implementations on FPGA's, microcontrollers and software simulation. The software simulation is supported by the TENNLab neuromorphic software framework so that researchers may develop RAVENS solutions for a variety of neuromorphic computing applications. This document provides a functional specification of RAVENS that should apply to all implementations of the RAVENS neuroprocessor. △ Less

Submitted 27 July, 2023; originally announced July 2023.

Comments: 17 pages, 11 figures

Enhanced Read Resolution in Reconfigurable Memristive Synapses for Spiking Neural Networks

Authors: Hritom Das, Nishith N. Chakraborty, Catherine Schuman, Garrett S. Rose

Abstract: Synapse is a key element of any neuromorphic computing system which is mostly constructed with memristor devices. A memristor is a two-terminal analog memory device. Memristive synapse suffers from various challenges such as forming at high voltage, SET, RESET failure, and READ margin or resolution issue between two weights. Enhanced READ resolution is very important to make a memristive synapse f… ▽ More Synapse is a key element of any neuromorphic computing system which is mostly constructed with memristor devices. A memristor is a two-terminal analog memory device. Memristive synapse suffers from various challenges such as forming at high voltage, SET, RESET failure, and READ margin or resolution issue between two weights. Enhanced READ resolution is very important to make a memristive synapse functionally reliable. Usually, the READ resolution is very small for a memristive synapse with 4-bit data precision. This work considers a step-by-step analysis to enhance the READ current resolution for a current-controlled memristor-based synapse. An empirical model is used to characterize the HfO2-based memristive device. 1st and 2nd stage device of our proposed synapse can be scaled to enhance the READ current margin up to ~ 4.3x and ~ 21% respectively. Moreover, READ current resolution can be enhanced with run-time adaptation features such as READ voltage scaling and body biasing. The READ voltage scaling and body biasing can improve the READ current resolution by about 46% and 15% respectively. TENNLabs' neuromorphic computing framework is leveraged to evaluate the effect of READ current resolution on classification applications. Higher READ current resolution shows better accuracy than lower resolution with different percentages of read noise scenarios. △ Less

Submitted 23 June, 2023; originally announced June 2023.

An Efficient and Accurate Memristive Memory for Array-based Spiking Neural Networks

Authors: Hritom Das, Rocco D. Febbo, SNB Tushar, Nishith N. Chakraborty, Maximilian Liehr, Nathaniel Cady, Garrett S. Rose

Abstract: Memristors provide a tempting solution for weighted synapse connections in neuromorphic computing due to their size and non-volatile nature. However, memristors are unreliable in the commonly used voltage-pulse-based programming approaches and require precisely shaped pulses to avoid programming failure. In this paper, we demonstrate a current-limiting-based solution that provides a more predictab… ▽ More Memristors provide a tempting solution for weighted synapse connections in neuromorphic computing due to their size and non-volatile nature. However, memristors are unreliable in the commonly used voltage-pulse-based programming approaches and require precisely shaped pulses to avoid programming failure. In this paper, we demonstrate a current-limiting-based solution that provides a more predictable analog memory behavior when reading and writing memristive synapses. With our proposed design READ current can be optimized by about 19x compared to the 1T1R design. Moreover, our proposed design saves about 9x energy compared to the 1T1R design. Our 3T1R design also shows promising write operation which is less affected by the process variation in MOSFETs and the inherent stochastic behavior of memristors. Memristors used for testing are hafnium oxide based and were fabricated in a 65nm hybrid CMOS-memristor process. The proposed design also shows linear characteristics between the voltage applied and the resulting resistance for the writing operation. The simulation and measured data show similar patterns with respect to voltage pulse-based programming and current compliance-based programming. We further observed the impact of this behavior on neuromorphic-specific applications such as a spiking neural network △ Less

Submitted 6 September, 2023; v1 submitted 10 June, 2023; originally announced June 2023.

Journal ref: IEEE Transactions on Circuits and Systems I: Regular Papers, 10 August 2023

Tidal Disruption Events from the Combined Effects of Two-Body Relaxation and the Eccentric Kozai-Lidov Mechanism

Authors: Denyz Melchor, Brenna Mockler, Smadar Naoz, Sanaea Rose, Enrico Ramirez-Ruiz

Abstract: Tidal disruption events (TDEs) take place when a star ventures too close to a supermassive black hole (SMBH) and becomes ruptured. One of the leading proposed physical mechanisms often invoked in the literature involves weak two-body interactions experienced by the population of stars within the host SMBH's sphere of influence, commonly referred to as two-body relaxation. This process can alter th… ▽ More Tidal disruption events (TDEs) take place when a star ventures too close to a supermassive black hole (SMBH) and becomes ruptured. One of the leading proposed physical mechanisms often invoked in the literature involves weak two-body interactions experienced by the population of stars within the host SMBH's sphere of influence, commonly referred to as two-body relaxation. This process can alter the angular momentum of stars at large distances and place them into nearly radial orbits, thus driving them to disruption. On the other hand, gravitational perturbations from an SMBH companion via the eccentric Kozai-Lidov (EKL) mechanism have also been proposed as a promising stellar disruption channel. Here we demonstrate that the combination of EKL and two-body relaxation in SMBH binaries is imperative for building a comprehensive picture of the rates of TDEs. Here we explore how the density profile of the surrounding stellar distribution and the binary orbital parameters of an SMBH companion influence the rate of TDEs. We show that this combined channel naturally produces disruptions at a rate that is consistent with observations and also naturally forms repeated TDEs, where a bound star is partially disrupted on multiple orbits. Recent observations show stars being disrupted in short-period orbits, which is challenging to explain when these mechanisms are considered independently. However, the diffusive effect of two-body relaxation, combined with the secular nature of the eccentricity excitations from EKL, is found to drive stars on short eccentric orbits at a much higher rate. △ Less

Submitted 13 June, 2023; v1 submitted 8 June, 2023; originally announced June 2023.

Comments: Fixed affiliation

Optimizations for a Current-Controlled Memristor-based Neuromorphic Synapse Design

Authors: Hritom Das, Rocco D. Febbo, Charlie P. Rizzo, Nishith N. Chakraborty, James S. Plank, Garrett S. Rose

Abstract: The synapse is a key element of neuromorphic computing in terms of efficiency and accuracy. In this paper, an optimized current-controlled memristive synapse circuit is proposed. Our proposed synapse demonstrates reliability in the face of process variation and the inherent stochastic behavior of memristors. Up to an 82% energy optimization can be seen during the SET operation over prior work. In… ▽ More The synapse is a key element of neuromorphic computing in terms of efficiency and accuracy. In this paper, an optimized current-controlled memristive synapse circuit is proposed. Our proposed synapse demonstrates reliability in the face of process variation and the inherent stochastic behavior of memristors. Up to an 82% energy optimization can be seen during the SET operation over prior work. In addition, the READ process shows up to 54% energy savings. Our current-controlled approach also provides more reliable programming over traditional programming methods. This design is demonstrated with a 4-bit memory precision configuration. Using a spiking neural network (SNN), a neuromorphic application analysis was performed with this precision configuration. Our optimized design showed up to 82% improvement in control applications and a 2.7x improvement in classification applications compared with other design cases. △ Less

Submitted 7 September, 2023; v1 submitted 25 May, 2023; originally announced May 2023.

Journal ref: IEEE Journal on Emerging and Selected Topics in Circuits and Systems, 2023

Extended X-ray absorption spectroscopy using an ultrashort pulse laboratory-scale laser-plasma accelerator

Authors: B. Kettle, C. Colgan, E. Los, E. Gerstmayr, M. J. V. Streeter, F. Albert, S. Astbury, R. A. Baggott, N. Cavanagh, K. Falk, T. I. Hyde, O. Lundh, P. P. Rajeev, D. Riley, S. J. Rose, G. Sarri, C. Spindloe, K. Svendsen, D. R. Symes, M. Smid, A. G. R. Thomas, C. Thornton, R. Watt, S. P. D. Mangles

Abstract: Laser-driven compact particle accelerators can provide ultrashort pulses of broadband X-rays, well suited for undertaking X-ray absorption spectroscopy measurements on a femtosecond timescale. Here the Extended X-ray Absorption Fine Structure (EXAFS) features of the K-edge of a copper sample have been observed over a 250 eV window in a single shot using a laser wakefield accelerator, providing inf… ▽ More Laser-driven compact particle accelerators can provide ultrashort pulses of broadband X-rays, well suited for undertaking X-ray absorption spectroscopy measurements on a femtosecond timescale. Here the Extended X-ray Absorption Fine Structure (EXAFS) features of the K-edge of a copper sample have been observed over a 250 eV window in a single shot using a laser wakefield accelerator, providing information on both the electronic and ionic structure simultaneously. This unique capability will allow the investigation of ultrafast processes, and in particular, probing high-energy-density matter and physics far-from-equilibrium where the sample refresh rate is slow and shot number is limited. For example, states that replicate the tremendous pressures and temperatures of planetary bodies or the conditions inside nuclear fusion reactions. Using high-power lasers to pump these samples also has the advantage of being inherently synchronised to the laser-driven X-ray probe. A perspective on the additional strengths of a laboratory-based ultrafast X-ray absorption source is presented. △ Less

Submitted 1 July, 2024; v1 submitted 17 May, 2023; originally announced May 2023.

Stellar Collisions in the Galactic Center: Massive Stars, Collision Remnants, and Missing Red Giants

Authors: Sanaea C. Rose, Smadar Naoz, Re'em Sari, Itai Linial

Abstract: Like most galaxies, the Milky Way harbors a supermassive black hole (SMBH) at its center, surrounded by a nuclear star cluster. In this dense star cluster, direct collisions can occur between stars before they evolve off the main-sequence. Using a statistical approach, we characterize the outcomes of these stellar collisions within the inner parsec of the Galactic Center (GC). Close to the SMBH, w… ▽ More Like most galaxies, the Milky Way harbors a supermassive black hole (SMBH) at its center, surrounded by a nuclear star cluster. In this dense star cluster, direct collisions can occur between stars before they evolve off the main-sequence. Using a statistical approach, we characterize the outcomes of these stellar collisions within the inner parsec of the Galactic Center (GC). Close to the SMBH, where the velocity dispersion is larger than the escape speed from a Sun-like star, collisions lead to mass loss. We find that the stellar population within $0.01$ pc is halved within about a Gyr because of destructive collisions. Additionally, we predict a diffuse population of peculiar low-mass stars in the GC. These stars have been divested of their outer layers in the inner $0.01$ pc before migrating to larger distances from the SMBH. Between $0.01$ and $0.1$ pc from the SMBH, collisions can result in mergers. Our results suggest that repeated collisions between lower mass stars can produce massive ($\gtrsim 10$ M$_\odot$) stars, and there may be $\sim 100$ of them residing in this region. We provide predictions on the number of G objects, dust and gas enshrouded stellar objects, that may result from main-sequence stellar collisions. Lastly, we comment on uncertainties in our model and possible connections between stellar collisions and the missing red giants in the GC. △ Less

Submitted 20 April, 2023; originally announced April 2023.

Comments: Submitted to ApJ. Comments welcome. 24 pages, 13 figures

Operational Research: Methods and Applications

Authors: Fotios Petropoulos, Gilbert Laporte, Emel Aktas, Sibel A. Alumur, Claudia Archetti, Hayriye Ayhan, Maria Battarra, Julia A. Bennell, Jean-Marie Bourjolly, John E. Boylan, Michèle Breton, David Canca, Laurent Charlin, Bo Chen, Cihan Tugrul Cicek, Louis Anthony Cox Jr, Christine S. M. Currie, Erik Demeulemeester, Li Ding, Stephen M. Disney, Matthias Ehrgott, Martin J. Eppler, Güneş Erdoğan, Bernard Fortz, L. Alberto Franco , et al. (57 additional authors not shown)

Abstract: Throughout its history, Operational Research has evolved to include a variety of methods, models and algorithms that have been applied to a diverse and wide range of contexts. This encyclopedic article consists of two main sections: methods and applications. The first aims to summarise the up-to-date knowledge and provide an overview of the state-of-the-art methods and key developments in the vari… ▽ More Throughout its history, Operational Research has evolved to include a variety of methods, models and algorithms that have been applied to a diverse and wide range of contexts. This encyclopedic article consists of two main sections: methods and applications. The first aims to summarise the up-to-date knowledge and provide an overview of the state-of-the-art methods and key developments in the various subdomains of the field. The second offers a wide-ranging list of areas where Operational Research has been applied. The article is meant to be read in a nonlinear fashion. It should be used as a point of reference or first-port-of-call for a diverse pool of readers: academics, researchers, students, and practitioners. The entries within the methods and applications sections are presented in alphabetical order. The authors dedicate this paper to the 2023 Turkey/Syria earthquake victims. We sincerely hope that advances in OR will play a role towards minimising the pain and suffering caused by this and future catastrophes. △ Less

Submitted 13 January, 2024; v1 submitted 24 March, 2023; originally announced March 2023.

Journal ref: Journal of the Operational Research Society (2024) 75(3)

Monte Carlo modelling of the linear Breit-Wheeler process within the GEANT4 framework

Authors: R. A. Watt, S. J. Rose, B. Kettle, S. P. D. Mangles

Abstract: A linear Breit-Wheeler module for the code Geant4 has been developed. This allows signal-to-noise ratio calculations of linear Breit-Wheeler detection experiments to be performed within a single framework. The interaction between two photon sources is modelled by treating one as a static field, then photons from the second source are sampled and tracked through the field. To increase the efficienc… ▽ More A linear Breit-Wheeler module for the code Geant4 has been developed. This allows signal-to-noise ratio calculations of linear Breit-Wheeler detection experiments to be performed within a single framework. The interaction between two photon sources is modelled by treating one as a static field, then photons from the second source are sampled and tracked through the field. To increase the efficiency of the module, we have used a Gaussian process regression, which can lead to an increase in the calculation rate by a factor of up to 1000. To demonstrate the capabilities of this module, we use it to perform a parameter scan, modelling an experiment based on that recently reported by Kettle et al. [1]. We show that colliding $50\,$fs duration $γ$-rays, produced through bremsstrahlung emission of a $100\,$pC, $2\,$GeV laser wakefield accelerator beam, with a $50\,$ps X-ray field, generated by a germanium burn-through foil heated to temperatures $>\,150\,$eV, this experiment is capable of producing $>1\,$ Breit-Wheeler pair per shot. △ Less

Submitted 9 February, 2023; originally announced February 2023.

Comments: 7 pages, 8 figures

SpK: A fast atomic and microphysics code for the high-energy-density regime

Authors: A. J. Crilly, N. P. L. Niasse, A. R. Fraser, D. A. Chapman, K. M. McLean, S. J. Rose, J. P. Chittenden

Abstract: SpK is part of the numerical codebase at Imperial College London used to model high energy density physics (HEDP) experiments. SpK is an efficient atomic and microphysics code used to perform detailed configuration accounting calculations of electronic and ionic stage populations, opacities and emissivities for use in post-processing and radiation hydrodynamics simulations. This is done using scre… ▽ More SpK is part of the numerical codebase at Imperial College London used to model high energy density physics (HEDP) experiments. SpK is an efficient atomic and microphysics code used to perform detailed configuration accounting calculations of electronic and ionic stage populations, opacities and emissivities for use in post-processing and radiation hydrodynamics simulations. This is done using screened hydrogenic atomic data supplemented by the NIST energy level database. An extended Saha model solves for chemical equilibrium with extensions for non-ideal physics, such as ionisation potential depression, and non thermal equilibrium corrections. A tree-heap (treap) data structure is used to store spectral data, such as opacity, which is dynamic thus allowing easy insertion of points around spectral lines without a-priori knowledge of the ion stage populations. Results from SpK are compared to other codes and descriptions of radiation transport solutions which use SpK data are given. The treap data structure and SpK's computational efficiency allows inline post-processing of 3D hydrodynamics simulations with a dynamically evolving spectrum stored in a treap. △ Less

Submitted 29 November, 2022; originally announced November 2022.

Disclosure of a Neuromorphic Starter Kit

Authors: James S. Plank, Bryson Gullett, Adam Z. Foshie, Garrett S. Rose, Catherine D. Schuman

Abstract: This paper presents a Neuromorphic Starter Kit, which has been designed to help a variety of research groups perform research, exploration and real-world demonstrations of brain-based, neuromorphic processors and hardware environments. A prototype kit has been built and tested. We explain the motivation behind the kit, its design and composition, and a prototype physical demonstration. This paper presents a Neuromorphic Starter Kit, which has been designed to help a variety of research groups perform research, exploration and real-world demonstrations of brain-based, neuromorphic processors and hardware environments. A prototype kit has been built and tested. We explain the motivation behind the kit, its design and composition, and a prototype physical demonstration. △ Less

Submitted 8 November, 2022; originally announced November 2022.

Comments: 4 pages, 3 figures

The Impact of Initial-Final Mass Relations on Black Hole Microlensing

Authors: Sam Rose, Casey Y. Lam, Jessica R. Lu, Michael Medford, Matthew W. Hosek Jr., Natasha S. Abrams, Emily Ramey, Sergiy S. Vasylyev

Abstract: Uncertainty in the initial-final mass relation (IFMR) has long been a problem in understanding the final stages of massive star evolution. One of the major challenges of constraining the IFMR is the difficulty of measuring the mass of non-luminous remnant objects (i.e. neutron stars and black holes). Gravitational wave detectors have opened the possibility of finding large numbers of compact objec… ▽ More Uncertainty in the initial-final mass relation (IFMR) has long been a problem in understanding the final stages of massive star evolution. One of the major challenges of constraining the IFMR is the difficulty of measuring the mass of non-luminous remnant objects (i.e. neutron stars and black holes). Gravitational wave detectors have opened the possibility of finding large numbers of compact objects in other galaxies, but all in merging binary systems. Gravitational lensing experiments using astrometry and photometry are capable of finding compact objects, both isolated and in binaries, in the Milky Way. In this work we improve the PopSyCLE microlensing simulation code in order to explore the possibility of constraining the IFMR using the Milky Way microlensing population. We predict that the Roman Space Telescope's microlensing survey will likely be able to distinguish different IFMRs based on the differences at the long end of the Einstein crossing time distribution and the small end of the microlensing parallax distribution, assuming the small ($π_E \lesssim 0.02$) microlensing parallaxes characteristic of black hole lenses are able to be measured accurately. We emphasize that future microlensing surveys need to be capable of characterizing events with small microlensing parallaxes in order to place the most meaningful constraints on the IFMR. △ Less

Submitted 7 November, 2022; originally announced November 2022.

Comments: 24 pages, 17 figures Accepted to ApJ

The Case for RISP: A Reduced Instruction Spiking Processor

Authors: James S. Plank, ChaoHui Zheng, Bryson Gullett, Nicholas Skuda, Charles Rizzo, Catherine D. Schuman, Garrett S. Rose

Abstract: In this paper, we introduce RISP, a reduced instruction spiking processor. While most spiking neuroprocessors are based on the brain, or notions from the brain, we present the case for a spiking processor that simplifies rather than complicates. As such, it features discrete integration cycles, configurable leak, and little else. We present the computing model of RISP and highlight the benefits of… ▽ More In this paper, we introduce RISP, a reduced instruction spiking processor. While most spiking neuroprocessors are based on the brain, or notions from the brain, we present the case for a spiking processor that simplifies rather than complicates. As such, it features discrete integration cycles, configurable leak, and little else. We present the computing model of RISP and highlight the benefits of its simplicity. We demonstrate how it aids in developing hand built neural networks for simple computational tasks, detail how it may be employed to simplify neural networks built with more complicated machine learning techniques, and demonstrate how it performs similarly to other spiking neurprocessors. △ Less

Submitted 28 June, 2022; originally announced June 2022.

Comments: 5 pages, 5 figures

Tidal disruption events from eccentric orbits and lessons learned from the noteworthy ASASSN-14ko

Authors: Chang Liu, Brenna Mockler, Enrico Ramirez-Ruiz, Ricardo Yarza, Jamie A. P. Law-Smith, Smadar Naoz, Denyz Melchor, Sanaea Rose

Abstract: Stars grazing supermassive black holes (SMBHs) on bound orbits may survive tidal disruption, causing periodic flares. Inspired by the recent discovery of the periodic nuclear transient ASASSN-14ko, a promising candidate for a repeating tidal disruption event (TDE), we study the tidal deformation of stars approaching SMBHs on eccentric orbits. With both analytical and hydrodynamics methods, we show… ▽ More Stars grazing supermassive black holes (SMBHs) on bound orbits may survive tidal disruption, causing periodic flares. Inspired by the recent discovery of the periodic nuclear transient ASASSN-14ko, a promising candidate for a repeating tidal disruption event (TDE), we study the tidal deformation of stars approaching SMBHs on eccentric orbits. With both analytical and hydrodynamics methods, we show the overall tidal deformation of a star is similar to that in a parabolic orbit provided that the eccentricity is above a critical value. This allows one to make use of existing simulation libraries from parabolic encounters to calculate the mass fallback rate in eccentric TDEs. We find the flare structures of eccentric TDEs show a complicated dependence on both the SMBH mass and the orbital period. For stars orbiting SMBHs with relatively short periods, we predict significantly shorter-lived duration flares than those in parabolic TDEs, which can be used to predict repeating events if the mass of the SMBH can be independently measured. Using an adiabatic mass loss model, we study the flare evolution over multiple passages, and show the evolved stars can survive many more passages than main sequence stars. We apply this theoretical framework to the repeating TDE candidate ASASSN-14ko and suggest that its recurrent flares originate from a moderately massive ($M\gtrsim 1 \mathrm{M_\odot}$), extended ($\approx$ a few $\mathrm{R_\odot}$), evolved star on a grazing, bound orbit around the SMBH. Future hydrodynamics simulations of multiple tidal interactions will enable realistic models on the individual flare structure and the evolution over multiple flares. △ Less

Submitted 24 January, 2023; v1 submitted 27 June, 2022; originally announced June 2022.

Comments: 19 pages, 9 figures. Updated to accepted version (ApJ)

Journal ref: ApJ, 944, 184 (2023)

L-shell X-ray conversion yields for laser-irradiated tin and silver foils

Authors: R. L. Singh, S. White, M. Charlwood, F. P. Keenan, C. Hyland, D. Bailie, T. Audet, G. Sarri, S. J. Rose, J. Morton, C. Baird, C. Spindloe, D. Riley

Abstract: We have employed the VULCAN laser facility to generate a laser plasma X-ray source for use in photoionisation experiments. A nanosecond laser pulse with an intensity of order ${10}^{15}$ W{cm}$^{-2}$ was used to irradiate thin Ag or Sn foil targets coated onto a parylene substrate, and the L-shell emission in the $3.3-4.4$ keV range was recorded for both the laser-irradiated and non-irradiated sid… ▽ More We have employed the VULCAN laser facility to generate a laser plasma X-ray source for use in photoionisation experiments. A nanosecond laser pulse with an intensity of order ${10}^{15}$ W{cm}$^{-2}$ was used to irradiate thin Ag or Sn foil targets coated onto a parylene substrate, and the L-shell emission in the $3.3-4.4$ keV range was recorded for both the laser-irradiated and non-irradiated sides. Both the experimental and simulation results show higher laser to X-ray conversion yields for Ag compared with Sn, with our simulations indicating yields approximately a factor of two higher than found in the experiments. Although detailed angular data were not available experimentally, the simulations indicate that the emission is quite isotropic on the laser-irradiated side, but shows close to a cosine variation on the non-irradiated side of the target as seen experimentally in previous work. △ Less

Submitted 23 April, 2022; originally announced April 2022.

Comments: Number of papges: 14, Number of figures: 12

MDsrv -- visual sharing and analysis of molecular dynamics simulations

Authors: Michelle Kampfrath, René Staritzbichler, Guillermo Pérez Hernández, Alexander S. Rose, Johanna K. S. Tiemann, Gerik Scheuermann, Daniel Wiegreffe, Peter W. Hildebrand

Abstract: Molecular dynamics simulation is a proven technique for computing and visualizing the time-resolved motion of macromolecules at atomic resolution. The MDsrv is a tool that streams MD trajectories and displays them interactively in web browsers without requiring advanced skills, facilitating interactive exploration and collaborative visual analysis. We have now enhanced the MDsrv to further simplif… ▽ More Molecular dynamics simulation is a proven technique for computing and visualizing the time-resolved motion of macromolecules at atomic resolution. The MDsrv is a tool that streams MD trajectories and displays them interactively in web browsers without requiring advanced skills, facilitating interactive exploration and collaborative visual analysis. We have now enhanced the MDsrv to further simplify the upload and sharing of MD trajectories and improve their online viewing and analysis. With the new instance, the MDsrv simplifies the creation of sessions, which allows the exchange of MD trajectories with preset representations and perspectives. An important innovation is that the MDsrv can now access and visualize trajectories from remote datasets, which greatly expands its applicability and use, as the data no longer needs to be accessible on a local server. In addition, initial analyses such as sequence or structure alignments, distance measurements, or RMSD calculations have been implemented, which optionally support visual analysis. Finally, the MDsrv now offers a faster and more efficient visualization of even large trajectories. △ Less

Submitted 4 May, 2022; v1 submitted 25 March, 2022; originally announced March 2022.

Comments: 9 pages, 3 figures

Investigating radiatively driven, magnetised plasmas with a university scale pulsed-power generator

Authors: Jack W. D. Halliday, Aidan Crilly, Jeremy Chittenden, Roberto C. Mancini, Stefano Merlini, Steven Rose, Danny R. Russell, Lee G. Suttle, Vicente Valenzuela-Villaseca, Simon N. Bland, Sergey V. Lebedev

Abstract: We present first results from a novel experimental platform which is able to access physics relevant to topics including indirect-drive magnetised ICF; laser energy deposition; various topics in atomic physics; and laboratory astrophysics (for example the penetration of B-fields into HED plasmas). This platform uses the X-Rays from a wire array Z-Pinch to irradiate a silicon target, producing an o… ▽ More We present first results from a novel experimental platform which is able to access physics relevant to topics including indirect-drive magnetised ICF; laser energy deposition; various topics in atomic physics; and laboratory astrophysics (for example the penetration of B-fields into HED plasmas). This platform uses the X-Rays from a wire array Z-Pinch to irradiate a silicon target, producing an outflow of ablated plasma. The ablated plasma expands into ambient, dynamically significant B-fields (~5 T) which are supported by the current flowing through the Z-Pinch. The outflows have a well-defined (quasi-1D) morphology, enabling the study of fundamental processes typically only available in more complex, integrated schemes. Experiments were fielded on the MAGPIE pulsed-power generator (1.4 MA, 240 ns rise time). On this machine a wire array Z-Pinch produces an X-Ray pulse carrying a total energy of ~15 kJ over ~30 ns. This equates to an average brightness temperature of around 10 eV on-target. △ Less

Submitted 22 March, 2022; originally announced March 2022.

Capturing Actionable Dynamics with Structured Latent Ordinary Differential Equations

Authors: Paidamoyo Chapfuwa, Sherri Rose, Lawrence Carin, Edward Meeds, Ricardo Henao

Abstract: End-to-end learning of dynamical systems with black-box models, such as neural ordinary differential equations (ODEs), provides a flexible framework for learning dynamics from data without prescribing a mathematical model for the dynamics. Unfortunately, this flexibility comes at the cost of understanding the dynamical system, for which ODEs are used ubiquitously. Further, experimental data are co… ▽ More End-to-end learning of dynamical systems with black-box models, such as neural ordinary differential equations (ODEs), provides a flexible framework for learning dynamics from data without prescribing a mathematical model for the dynamics. Unfortunately, this flexibility comes at the cost of understanding the dynamical system, for which ODEs are used ubiquitously. Further, experimental data are collected under various conditions (inputs), such as treatments, or grouped in some way, such as part of sub-populations. Understanding the effects of these system inputs on system outputs is crucial to have any meaningful model of a dynamical system. To that end, we propose a structured latent ODE model that explicitly captures system input variations within its latent representation. Building on a static latent variable specification, our model learns (independent) stochastic factors of variation for each input to the system, thus separating the effects of the system inputs in the latent space. This approach provides actionable modeling through the controlled generation of time-series data for novel input combinations (or perturbations). Additionally, we propose a flexible approach for quantifying uncertainties, leveraging a quantile regression formulation. Results on challenging biological datasets show consistent improvements over competitive baselines in the controlled generation of observational data and inference of biologically meaningful system inputs. △ Less

Submitted 16 June, 2022; v1 submitted 25 February, 2022; originally announced February 2022.

Comments: Accepted for the 38th Conference on Uncertainty in Artificial Intelligence (UAI 2022). Github code can be found at https://github.com/paidamoyo/structured_latent_ODEs

The Combined Effects of Two-Body Relaxation Processes and the Eccentric Kozai-Lidov Mechanism on the EMRI Rate

Authors: Smadar Naoz, Sanaea C. Rose, Erez Michaely, Denyz Melchor, Enrico Ramirez-Ruiz, Brenna Mockler, Jeremy D. Schnittman

Abstract: Gravitational wave (GW) emissions from extreme-mass-ratio inspirals (EMRIs) are promising sources for low-frequency GW-detectors. They result from a compact object, such as a stellar-mass black-hole (BH), captured by a supermassive black hole (SMBH). Several physical processes have been proposed to form EMRIs. In particular, weak two-body interactions over a long time scale (i.e., relaxation proce… ▽ More Gravitational wave (GW) emissions from extreme-mass-ratio inspirals (EMRIs) are promising sources for low-frequency GW-detectors. They result from a compact object, such as a stellar-mass black-hole (BH), captured by a supermassive black hole (SMBH). Several physical processes have been proposed to form EMRIs. In particular, weak two-body interactions over a long time scale (i.e., relaxation processes) have been proposed as a likely mechanism to drive the BH orbit to high eccentricity. Consequently, it is captured by the SMBH and becomes an EMRI. Here we demonstrate that EMRIs are naturally formed in SMBH binaries. Gravitational perturbations from an SMBH companion, known as the eccentric Kozai-Lidov (EKL) mechanism, combined with relaxation processes, yield a significantly more enhanced rate than any of these processes operating alone. Since EKL is sensitive to the orbital configuration, two-body relaxation can alter the orbital parameters, rendering the system in a more EKL-favorable regime. As SMBH binaries are expected to be prevalent in the Universe, this process predicts a substantially high EMRI rate. △ Less

Submitted 21 March, 2022; v1 submitted 24 February, 2022; originally announced February 2022.

Comments: 11 pages, 6 figures, accepted to ApJ-Letters

An isolated mass gap black hole or neutron star detected with astrometric microlensing

Authors: Casey Y. Lam, Jessica R. Lu, Andrzej Udalski, Ian Bond, David P. Bennett, Jan Skowron, Przemek Mroz, Radek Poleski, Takahiro Sumi, Michal K. Szymanski, Szymon Kozlowski, Pawel Pietrukowicz, Igor Soszynski, Krzysztof Ulaczyk, Lukasz Wyrzykowski, Shota Miyazaki, Daisuke Suzuki, Naoki Koshimoto, Nicholas J. Rattenbury, Matthew W. Hosek Jr., Fumio Abe, Richard Barry, Aparna Bhattacharya, Akihiko Fukui, Hirosane Fujii , et al. (20 additional authors not shown)

Abstract: We present the analysis of five black hole candidates identified from gravitational microlensing surveys. Hubble Space Telescope astrometric data and densely sampled lightcurves from ground-based microlensing surveys are fit with a single-source, single-lens microlensing model in order to measure the mass and luminosity of each lens and determine if it is a black hole. One of the five targets (OGL… ▽ More We present the analysis of five black hole candidates identified from gravitational microlensing surveys. Hubble Space Telescope astrometric data and densely sampled lightcurves from ground-based microlensing surveys are fit with a single-source, single-lens microlensing model in order to measure the mass and luminosity of each lens and determine if it is a black hole. One of the five targets (OGLE-2011-BLG-0462/MOA-2011-BLG-191 or OB110462 for short) shows a significant $>1$ mas coherent astrometric shift, little to no lens flux, and has an inferred lens mass of 1.6 - 4.4 $M_\odot$. This makes OB110462 the first definitive discovery of a compact object through astrometric microlensing and it is most likely either a neutron star or a low-mass black hole. This compact object lens is relatively nearby (0.70-1.92 kpc) and has a slow transverse motion of $<$30 km/s. OB110462 shows significant tension between models well-fit to photometry vs. astrometry, making it currently difficult to distinguish between a neutron star and a black hole. Additional observations and modeling with more complex system geometries, such as binary sources are needed to resolve the puzzling nature of this object. For the remaining four candidates, the lens masses are $<2 M_\odot$ and they are unlikely to be black holes; two of the four are likely white dwarfs or neutron stars. We compare the full sample of five candidates to theoretical expectations on the number of black holes in the Milky Way ($\sim 10^8$) and find reasonable agreement given the small sample size. △ Less

Submitted 31 May, 2022; v1 submitted 3 February, 2022; originally announced February 2022.

Comments: Accepted to ApJ Letters, with corresponding ApJ Supplement. 10 page Letter (6 figures, 2 tables) + 51 page Supplement (27 figures, 20 tables, 9 appendices). Some minor updates from the refereeing process, but no change to main conclusions

The Formation of Intermediate Mass Black Holes in Galactic Nuclei

Authors: Sanaea C. Rose, Smadar Naoz, Re'em Sari, Itai Linial

Abstract: Most stellar evolution models predict that black holes (BHs) should not exist above approximately $50-70$ M$_\odot$, the lower limit of the pair-instability mass gap. However, recent LIGO/Virgo detections indicate the existence of BHs with masses at and above this threshold. We suggest that massive BHs, including intermediate mass black holes (IMBHs), can form in galactic nuclei through collisions… ▽ More Most stellar evolution models predict that black holes (BHs) should not exist above approximately $50-70$ M$_\odot$, the lower limit of the pair-instability mass gap. However, recent LIGO/Virgo detections indicate the existence of BHs with masses at and above this threshold. We suggest that massive BHs, including intermediate mass black holes (IMBHs), can form in galactic nuclei through collisions between stellar-mass black holes and the surrounding main-sequence stars. Considering dynamical processes such as collisions, mass segregation, and relaxation, we find that this channel can be quite efficient, forming IMBHs as massive as $10^4$ M$_\odot$. This upper limit assumes that (1) the BHs accrete a substantial fraction of the stellar mass captured during each collision and (2) that the rate at which new stars are introduced into the region near the SMBH is high enough to offset depletion by stellar disruptions and star-star collisions. We discuss deviations from these key assumptions in the text. Our results suggest that BHs in the pair-instability mass gap and IMBHs may be ubiquitous in galactic centers. This formation channel has implications for observations. Collisions between stars and BHs can produce electromagnetic signatures, for example, from x-ray binaries and tidal disruption events. Additionally, formed through this channel, both black holes in the mass gap and IMBHs can merge with the supermassive black hole at the center of a galactic nucleus through gravitational waves. These gravitational wave events are extreme and intermediate mass ratio inspirals (EMRIs and IMRIs, respectively). △ Less

Submitted 6 July, 2022; v1 submitted 31 December, 2021; originally announced January 2022.

The data-driven future of high energy density physics

Authors: Peter W. Hatfield, Jim A. Gaffney, Gemma J. Anderson, Suzanne Ali, Luca Antonelli, Suzan Başeğmez du Pree, Jonathan Citrin, Marta Fajardo, Patrick Knapp, Brendan Kettle, Bogdan Kustowski, Michael J. MacDonald, Derek Mariscal, Madison E. Martin, Taisuke Nagayama, Charlotte A. J. Palmer, J. Luc Peterson, Steven Rose, J J Ruby, Carl Shneider, Matt J. V. Streeter, Will Trickey, Ben Williams

Abstract: The study of plasma physics under conditions of extreme temperatures, densities and electromagnetic field strengths is significant for our understanding of astrophysics, nuclear fusion and fundamental physics. These extreme physical systems are strongly non-linear and very difficult to understand theoretically or optimize experimentally. Here, we argue that machine learning models and data-driven… ▽ More The study of plasma physics under conditions of extreme temperatures, densities and electromagnetic field strengths is significant for our understanding of astrophysics, nuclear fusion and fundamental physics. These extreme physical systems are strongly non-linear and very difficult to understand theoretically or optimize experimentally. Here, we argue that machine learning models and data-driven methods are in the process of reshaping our exploration of these extreme systems that have hitherto proven far too non-linear for human researchers. From a fundamental perspective, our understanding can be helped by the way in which machine learning models can rapidly discover complex interactions in large data sets. From a practical point of view, the newest generation of extreme physics facilities can perform experiments multiple times a second (as opposed to ~daily), moving away from human-based control towards automatic control based on real-time interpretation of diagnostic data and updates of the physics model. To make the most of these emerging opportunities, we advance proposals for the community in terms of research design, training, best practices, and support for synthetic diagnostics and data analysis. △ Less

Submitted 22 November, 2021; originally announced November 2021.

Comments: 14 pages, 4 figures. This work was the result of a meeting at the Lorentz Center, University of Leiden, 13th-17th January 2020. This is a preprint of Hatfield et al., Nature, 593, 7859, 351-361 (2021) https://www.nature.com/articles/s41586-021-03382-w

Journal ref: Nature, 593, 7859, 351-361, 2021

Why Machine Learning Cannot Ignore Maximum Likelihood Estimation

Authors: Mark J. van der Laan, Sherri Rose

Abstract: The growth of machine learning as a field has been accelerating with increasing interest and publications across fields, including statistics, but predominantly in computer science. How can we parse this vast literature for developments that exemplify the necessary rigor? How many of these manuscripts incorporate foundational theory to allow for statistical inference? Which advances have the great… ▽ More The growth of machine learning as a field has been accelerating with increasing interest and publications across fields, including statistics, but predominantly in computer science. How can we parse this vast literature for developments that exemplify the necessary rigor? How many of these manuscripts incorporate foundational theory to allow for statistical inference? Which advances have the greatest potential for impact in practice? One could posit many answers to these queries. Here, we assert that one essential idea is for machine learning to integrate maximum likelihood for estimation of functional parameters, such as prediction functions and conditional densities. △ Less

Submitted 22 October, 2021; originally announced October 2021.

Comments: 30 pages. Forthcoming as a chapter in the Handbook of Matching and Weighting in Causal Inference

Raining Rocks: An analytical formulation for collision timescales in planetary systems

Authors: Santiago Torres, Smadar Naoz, Gongjie Li, Sanaea C. Rose

Abstract: The dynamical interaction of minor bodies (such as comets or asteroids) with planets plays an essential role in the planetary system's architecture and evolution. As a result of these interactions, structures like the Kuiper belt and the Oort cloud can be created. In particular, the collision of minor bodies with planets can drastically change the planet's internal and orbital evolution. We presen… ▽ More The dynamical interaction of minor bodies (such as comets or asteroids) with planets plays an essential role in the planetary system's architecture and evolution. As a result of these interactions, structures like the Kuiper belt and the Oort cloud can be created. In particular, the collision of minor bodies with planets can drastically change the planet's internal and orbital evolution. We present an analytic formulation to determine the collision timescale for a minor body to impact a planet for arbitrary geometry. By comparing with a suite of detailed N-body simulations and an analytical method for collision timescales in the solar system, we confirmed the accuracy of our analytic formulation. As a proof of concept, we focused on the collision rate of minor bodies randomly distributed around a Jupiter-like planet, emulating a Kuiper belt-like disk. We show that our analytical method yields in good agreement with the numerical simulations. The formalism presented here thus provides a succinct and accurate alternative to numerical calculations. △ Less

Submitted 13 June, 2022; v1 submitted 5 October, 2021; originally announced October 2021.

Comments: 10 pages, 7 figures

Conditional Cross-Design Synthesis Estimators for Generalizability in Medicaid

Authors: Irina Degtiar, Tim Layton, Jacob Wallace, Sherri Rose

Abstract: While much of the causal inference literature has focused on addressing internal validity biases, both internal and external validity are necessary for unbiased estimates in a target population of interest. However, few generalizability approaches exist for estimating causal quantities in a target population when the target population is not well-represented by a randomized study but is reflected… ▽ More While much of the causal inference literature has focused on addressing internal validity biases, both internal and external validity are necessary for unbiased estimates in a target population of interest. However, few generalizability approaches exist for estimating causal quantities in a target population when the target population is not well-represented by a randomized study but is reflected when additionally incorporating observational data. To generalize to a target population represented by a union of these data, we propose a class of novel conditional cross-design synthesis estimators that combine randomized and observational data, while addressing their respective biases. The estimators include outcome regression, propensity weighting, and double robust approaches. All use the covariate overlap between the randomized and observational data to remove potential unmeasured confounding bias. We apply these methods to estimate the causal effect of managed care plans on health care spending among Medicaid beneficiaries in New York City. △ Less

Submitted 27 September, 2021; originally announced September 2021.

Comments: 25 pages, 4 figures; supplement of 31 pages, 12 figures, and 4 tables

MSC Class: 62G05 (Primary); 62P25 (Secondary)

Journal ref: Biometrics (2023)

Uncertainty in Lung Cancer Stage for Outcome Estimation via Set-Valued Classification

Authors: Savannah Bergquist, Gabriel Brooks, Mary Beth Landrum, Nancy Keating, Sherri Rose

Abstract: Difficulty in identifying cancer stage in health care claims data has limited oncology quality of care and health outcomes research. We fit prediction algorithms for classifying lung cancer stage into three classes (stages I/II, stage III, and stage IV) using claims data, and then demonstrate a method for incorporating the classification uncertainty in outcomes estimation. Leveraging set-valued cl… ▽ More Difficulty in identifying cancer stage in health care claims data has limited oncology quality of care and health outcomes research. We fit prediction algorithms for classifying lung cancer stage into three classes (stages I/II, stage III, and stage IV) using claims data, and then demonstrate a method for incorporating the classification uncertainty in outcomes estimation. Leveraging set-valued classification and split conformal inference, we show how a fixed algorithm developed in one cohort of data may be deployed in another, while rigorously accounting for uncertainty from the initial classification step. We demonstrate this process using SEER cancer registry data linked with Medicare claims data. △ Less

Submitted 2 July, 2021; originally announced July 2021.

Comments: Code available at: https://github.com/sl-bergquist/cancer_classification

Journal ref: Statistics in Medicine (2022)

A laser-plasma platform for photon-photon physics

Authors: B. Kettle, D. Hollatz, E. Gerstmayr, G. M. Samarin, A. Alejo, S. Astbury, C. Baird, S. Bohlen, M. Campbell, C. Colgan, D. Dannheim, C. Gregory, H. Harsh, P. Hatfield, J. Hinojosa, Y. Katzir, J. Morton, C. D. Murphy, A. Nurnberg, J. Osterhoff, G. Pérez-Callejo, K. Poder, P. P. Rajeev, C. Roedel, F. Roeder , et al. (13 additional authors not shown)

Abstract: We describe a laser-plasma platform for photon-photon collision experiments to measure fundamental quantum electrodynamic processes such as the linear Breit-Wheeler process with real photons. The platform has been developed using the Gemini laser facility at the Rutherford Appleton Laboratory. A laser wakefield accelerator and a bremsstrahlung convertor are used to generate a collimated beam of ph… ▽ More We describe a laser-plasma platform for photon-photon collision experiments to measure fundamental quantum electrodynamic processes such as the linear Breit-Wheeler process with real photons. The platform has been developed using the Gemini laser facility at the Rutherford Appleton Laboratory. A laser wakefield accelerator and a bremsstrahlung convertor are used to generate a collimated beam of photons with energies of hundreds of MeV, that collide with keV x-ray photons generated by a laser heated plasma target. To detect the pairs generated by the photon-photon collisions, a magnetic transport system has been developed which directs the pairs onto scintillation-based and hybrid silicon pixel single particle detectors. We present commissioning results from an experimental campaign using this laser-plasma platform for photon-photon physics, demonstrating successful generation of both photon sources, characterisation of the magnetic transport system and calibration of the single particle detectors, and discuss the feasibility of this platform for the observation of the Breit-Wheeler process. The design of the platform will also serve as the basis for the investigation of strong-field quantum electrodynamic processes such as the nonlinear Breit-Wheeler and the Trident process, or eventually, photon-photon scattering. △ Less

Submitted 5 July, 2021; v1 submitted 29 June, 2021; originally announced June 2021.

Comments: 28 pages, 14 figures

Journal ref: New Journal of Physics, 23 (2021) 115006

Identifying Undercompensated Groups Defined By Multiple Attributes in Risk Adjustment

Authors: Anna Zink, Sherri Rose

Abstract: Risk adjustment in health care aims to redistribute payments to insurers based on costs. However, risk adjustment formulas are known to underestimate costs for some groups of patients. This undercompensation makes these groups unprofitable to insurers and creates incentives for insurers to discriminate. We develop a machine learning method for "group importance" to identify unprofitable groups def… ▽ More Risk adjustment in health care aims to redistribute payments to insurers based on costs. However, risk adjustment formulas are known to underestimate costs for some groups of patients. This undercompensation makes these groups unprofitable to insurers and creates incentives for insurers to discriminate. We develop a machine learning method for "group importance" to identify unprofitable groups defined by multiple attributes, improving on the arbitrary nature of existing evaluations. This procedure was designed to evaluate the risk adjustment formulas used in the U.S. health insurance Marketplaces as well as Medicare. We find that a number of previously unidentified groups with multiple chronic conditions are undercompensated in the Marketplaces risk adjustment formula, while groups without chronic conditions tend to be overcompensated in the Marketplaces. The magnitude of undercompensation when defining groups with multiple attributes is larger than with single attributes. No complex groups were found to be consistently under- or overcompensated in the Medicare risk adjustment formula. Our work provides policy makers with new information on potential targets of discrimination in the health care system and a path towards more equitable health coverage. △ Less

Submitted 26 July, 2021; v1 submitted 18 May, 2021; originally announced May 2021.

Journal ref: BMJ Health & Care Informatics (2021)

Considerations Across Three Cultures: Parametric Regressions, Interpretable Algorithms, and Complex Algorithms

Authors: Ani Eloyan, Sherri Rose

Abstract: We consider an extension of Leo Breiman's thesis from "Statistical Modeling: The Two Cultures" to include a bifurcation of algorithmic modeling, focusing on parametric regressions, interpretable algorithms, and complex (possibly explainable) algorithms. We consider an extension of Leo Breiman's thesis from "Statistical Modeling: The Two Cultures" to include a bifurcation of algorithmic modeling, focusing on parametric regressions, interpretable algorithms, and complex (possibly explainable) algorithms. △ Less

Submitted 13 April, 2021; originally announced April 2021.

Comments: 7 pages, forthcoming in Observational Studies

Journal ref: Observational Studies (2021); 7(1):191-196. https://muse.jhu.edu/article/799734

Ariel: Enabling planetary science across light-years

Authors: Giovanna Tinetti, Paul Eccleston, Carole Haswell, Pierre-Olivier Lagage, Jérémy Leconte, Theresa Lüftinger, Giusi Micela, Michel Min, Göran Pilbratt, Ludovic Puig, Mark Swain, Leonardo Testi, Diego Turrini, Bart Vandenbussche, Maria Rosa Zapatero Osorio, Anna Aret, Jean-Philippe Beaulieu, Lars Buchhave, Martin Ferus, Matt Griffin, Manuel Guedel, Paul Hartogh, Pedro Machado, Giuseppe Malaguti, Enric Pallé , et al. (293 additional authors not shown)

Abstract: Ariel, the Atmospheric Remote-sensing Infrared Exoplanet Large-survey, was adopted as the fourth medium-class mission in ESA's Cosmic Vision programme to be launched in 2029. During its 4-year mission, Ariel will study what exoplanets are made of, how they formed and how they evolve, by surveying a diverse sample of about 1000 extrasolar planets, simultaneously in visible and infrared wavelengths.… ▽ More Ariel, the Atmospheric Remote-sensing Infrared Exoplanet Large-survey, was adopted as the fourth medium-class mission in ESA's Cosmic Vision programme to be launched in 2029. During its 4-year mission, Ariel will study what exoplanets are made of, how they formed and how they evolve, by surveying a diverse sample of about 1000 extrasolar planets, simultaneously in visible and infrared wavelengths. It is the first mission dedicated to measuring the chemical composition and thermal structures of hundreds of transiting exoplanets, enabling planetary science far beyond the boundaries of the Solar System. The payload consists of an off-axis Cassegrain telescope (primary mirror 1100 mm x 730 mm ellipse) and two separate instruments (FGS and AIRS) covering simultaneously 0.5-7.8 micron spectral range. The satellite is best placed into an L2 orbit to maximise the thermal stability and the field of regard. The payload module is passively cooled via a series of V-Groove radiators; the detectors for the AIRS are the only items that require active cooling via an active Ne JT cooler. The Ariel payload is developed by a consortium of more than 50 institutes from 16 ESA countries, which include the UK, France, Italy, Belgium, Poland, Spain, Austria, Denmark, Ireland, Portugal, Czech Republic, Hungary, the Netherlands, Sweden, Norway, Estonia, and a NASA contribution. △ Less

Submitted 10 April, 2021; originally announced April 2021.

Comments: Ariel Definition Study Report, 147 pages. Reviewed by ESA Science Advisory Structure in November 2020. Original document available at: https://www.cosmos.esa.int/documents/1783156/3267291/Ariel_RedBook_Nov2020.pdf/

Report number: ESA/SCI(2020)1

Temperature Equilibration Due to Charge State Fluctuations in Dense Plasmas

Authors: R. A. Baggott, S. J. Rose, S. P. D. Mangles

Abstract: The charge states of ions in dense plasmas fluctuate due to collisional ionization and recombination. Here we show how, by modifying the ion interaction potential, these fluctuations can mediate energy exchange between the plasma electrons and ions. Moreover, we develop a theory for this novel electron-ion energy transfer mechanism. Calculations using a random walk approach for the fluctuations su… ▽ More The charge states of ions in dense plasmas fluctuate due to collisional ionization and recombination. Here we show how, by modifying the ion interaction potential, these fluctuations can mediate energy exchange between the plasma electrons and ions. Moreover, we develop a theory for this novel electron-ion energy transfer mechanism. Calculations using a random walk approach for the fluctuations suggest that the energy exchange rate from charge state fluctuations could be comparable to direct electron-ion collisions. This mechanism is, however, predicted to exhibit a complex dependence on the temperature and ionization state of the plasma, which could contribute to our understanding of significant variation in experimental measurements of equilibration times. △ Less

Submitted 14 July, 2021; v1 submitted 5 March, 2021; originally announced March 2021.

Journal ref: Phys. Rev. Lett. 127, 035002 (2021)

A Review of Generalizability and Transportability

Authors: Irina Degtiar, Sherri Rose

Abstract: When assessing causal effects, determining the target population to which the results are intended to generalize is a critical decision. Randomized and observational studies each have strengths and limitations for estimating causal effects in a target population. Estimates from randomized data may have internal validity but are often not representative of the target population. Observational data… ▽ More When assessing causal effects, determining the target population to which the results are intended to generalize is a critical decision. Randomized and observational studies each have strengths and limitations for estimating causal effects in a target population. Estimates from randomized data may have internal validity but are often not representative of the target population. Observational data may better reflect the target population, and hence be more likely to have external validity, but are subject to potential bias due to unmeasured confounding. While much of the causal inference literature has focused on addressing internal validity bias, both internal and external validity are necessary for unbiased estimates in a target population. This paper presents a framework for addressing external validity bias, including a synthesis of approaches for generalizability and transportability, the assumptions they require, as well as tests for the heterogeneity of treatment effects and differences between study and target populations. △ Less

Submitted 23 February, 2021; originally announced February 2021.

Comments: 30 pages, 3 figures

MSC Class: 62-02

Journal ref: Annual Review of Statistics and Its Application (2023)

Design of an Enhanced Reconfigurable Chaotic Oscillator using G4FET-NDR Based Discrete Map

Authors: Md Sakib Hasan, Aysha S. Shanta, Partha Sarathi Paul, Maisha Sadia, Md Badruddoja Majumder, Garrett S. Rose

Abstract: In this paper, a novel chaotic map is introduced usinga voltage controlled negative differential resistance (NDR) circuitcomposed of ann-channel and ap-channel silicon-on-insulator(SOI) four-gate transistor (G4FET). The multiple gates of theG4FET are leveraged to create a discrete chaotic map with threebifurcation parameters. The three tunable parameters are thegain of a transimpedance amplifier (… ▽ More In this paper, a novel chaotic map is introduced usinga voltage controlled negative differential resistance (NDR) circuitcomposed of ann-channel and ap-channel silicon-on-insulator(SOI) four-gate transistor (G4FET). The multiple gates of theG4FET are leveraged to create a discrete chaotic map with threebifurcation parameters. The three tunable parameters are thegain of a transimpedance amplifier (TIA), top-gate voltage ofn-channel G4FET, and top-gate voltage ofp-channel G4FET. Twomethods are proposed for building chaotic oscillators using thisdiscrete map. The effect of altering bifurcation parameters onchaotic operation is illustrated using bifurcation diagrams andLyapunov exponent. A design methodology for building flexibleand reconfigurable logic gate is outlined and the consequentenhancement in functionality space caused by the existence ofthree independent bifurcation parameters is demonstrated andcompared with previous work. △ Less

Submitted 1 January, 2021; originally announced January 2021.

Comments: This paper is accepted and presented in '14TH IEEE DALLAS CIRCUITS AND SYSTEMS CONFERENCE', November 16, 2020. The name of the paper can be found in conference schedule: https://engineering.utdallas.edu/DCAS/schedule.html

Designing an Interactive Visualization System for Monitoring Participant Compliance in a Large-Scale, Longitudinal Study

Authors: Poorna Talkad Sukumar, Thomas Breideband, Gonzalo Martinez, Megan Caruso, Sierra Rose, Cooper Steputis, Sidney D'Mello, Gloria Mark, Aaron Striegel

Abstract: Frequent monitoring of participant compliance is necessary when conducting large-scale, longitudinal studies to ensure that the collected data is of sufficiently high quality. While the need for achieving high compliance has been underscored and there are discussions on incentives and factors affecting compliance, little is shared about the actual processes and tools used for monitoring compliance… ▽ More Frequent monitoring of participant compliance is necessary when conducting large-scale, longitudinal studies to ensure that the collected data is of sufficiently high quality. While the need for achieving high compliance has been underscored and there are discussions on incentives and factors affecting compliance, little is shared about the actual processes and tools used for monitoring compliance in such studies. Monitoring participant compliance with respect to multi-modal data can be a tedious process, especially if there are only a few personnel involved. In this case study, we describe the iterative design of an interactive visualization system we developed for monitoring compliance and refined based on changing requirements in an ongoing study. We find that the visualization system, leveraging the digital medium, both facilitates the exploratory tasks of monitoring participant compliance and supports asynchronous collaboration among non-co-located researchers. Our documented requirements for checking participant compliance as well as the design of the visualization system can help inform the compliance-monitoring process in future studies. △ Less

Submitted 22 December, 2020; originally announced December 2020.

Comments: 8 pages, 4 figures

Ethical Machine Learning in Health Care

Authors: Irene Y. Chen, Emma Pierson, Sherri Rose, Shalmali Joshi, Kadija Ferryman, Marzyeh Ghassemi

Abstract: The use of machine learning (ML) in health care raises numerous ethical concerns, especially as models can amplify existing health inequities. Here, we outline ethical considerations for equitable ML in the advancement of health care. Specifically, we frame ethics of ML in health care through the lens of social justice. We describe ongoing efforts and outline challenges in a proposed pipeline of e… ▽ More The use of machine learning (ML) in health care raises numerous ethical concerns, especially as models can amplify existing health inequities. Here, we outline ethical considerations for equitable ML in the advancement of health care. Specifically, we frame ethics of ML in health care through the lens of social justice. We describe ongoing efforts and outline challenges in a proposed pipeline of ethical ML in health, ranging from problem selection to post-deployment considerations. We close by summarizing recommendations to address these challenges. △ Less

Submitted 7 October, 2020; v1 submitted 22 September, 2020; originally announced September 2020.

Comments: Annual Reviews in Biomedical Data Science 2021

On Socially Distant Neighbors: Using Binaries to Constrain the Density of Objects in the Galactic Center

Authors: Sanaea C. Rose, Smadar Naoz, Abhimat K. Gautam, Andrea M. Ghez, Tuan Do, Devin Chu, Eric Becklin

Abstract: Stars often reside in binary configurations. The nuclear star cluster surrounding the supermassive black hole (SMBH) in the Galactic Center (GC) is expected to include a binary population. In this dense environment, a binary frequently encounters and interacts with neighboring stars. These interactions vary from small perturbations to violent collisions. In the former case, weak gravitational inte… ▽ More Stars often reside in binary configurations. The nuclear star cluster surrounding the supermassive black hole (SMBH) in the Galactic Center (GC) is expected to include a binary population. In this dense environment, a binary frequently encounters and interacts with neighboring stars. These interactions vary from small perturbations to violent collisions. In the former case, weak gravitational interactions unbind a soft binary over the evaporation timescale, which depends on the binary properties as well as the density of surrounding objects and velocity dispersion. Similarly, collisions can also unbind a binary, and the collision rate depends on the density. Thus, the detection of a binary with known properties can constrain the density profile in the GC with implications for the number of compact objects, which are otherwise challenging to detect. We estimate the density necessary to unbind a binary within its lifetime for an orbit of arbitrary eccentricity about the SMBH. We find that the eccentricity has a minimal impact on the density constraint. In this proof of concept, we demonstrate that this procedure can probe the density in the GC using hypothetical young and old binaries as examples. Similarly, a known density profile provides constraints on the binary orbital separation. Our results highlight the need to consider multiple dynamical processes in tandem. In certain cases, often closer to the SMBH, the collision timescale rather than the evaporation timescale gives the more stringent density constraint, while other binaries farther from the SMBH provide unreliable density constraints because they migrate inwards due to mass segregation. △ Less

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

Comments: Accepted to ApJ, 20 pages, 11 figures