Showing 1–50 of 2,302 results for author: Smith, R

An infrared census of R Coronae Borealis Stars II -- Spectroscopic classifications and implications for the rate of low-mass white dwarf mergers

Authors: Viraj R. Karambelkar, Mansi M. Kasliwal, Patrick Tisserand, Shreya Anand, Michael C. B. Ashley, Lars Bildsten, Geoffrey C. Clayton, Courtney C. Crawford, Kishalay De, Nicholas Earley, Matthew J. Hankins, Xander Hall, Astrid Lamberts, Ryan M. Lau, Dan McKenna, Anna Moore, Eran O. Ofek, Roger M. Smith, Roberto Soria, Jamie Soon, Tony Travouillon

Abstract: We present results from a systematic infrared (IR) census of R Coronae Borealis (RCB) stars in the Milky Way, using data from the Palomar Gattini IR (PGIR) survey. R Coronae Borealis stars are dusty, erratic variable stars presumably formed from the merger of a He-core and a CO-core white dwarf (WD). PGIR is a 30 cm $J$-band telescope with a 25 deg$^{2}$ camera that surveys 18000 deg$^{2}$ of the… ▽ More We present results from a systematic infrared (IR) census of R Coronae Borealis (RCB) stars in the Milky Way, using data from the Palomar Gattini IR (PGIR) survey. R Coronae Borealis stars are dusty, erratic variable stars presumably formed from the merger of a He-core and a CO-core white dwarf (WD). PGIR is a 30 cm $J$-band telescope with a 25 deg$^{2}$ camera that surveys 18000 deg$^{2}$ of the northern sky ($δ>-28^{o}$) at a cadence of 2 days. Using PGIR J-band lightcurves for $\sim$60 million stars together with mid-IR colors from WISE, we selected a sample of 530 candidate RCB stars. We obtained near-IR spectra for these candidates and identified 53 RCB stars in our sample. Accounting for our selection criteria, we find that there are a total of $\approx350^{+150}_{-100}$ RCB stars in the Milky Way. Assuming typical RCB lifetimes, this corresponds to an RCB formation rate of 0.8 - 5 $\times$ 10$^{-3}$ yr$^{-1}$, consistent with observational and theoretical estimates of the He-CO WD merger rate. We searched for quasi-periodic pulsations in the PGIR lightcurves of RCB stars and present pulsation periods for 16 RCB stars. We also examined high-cadenced TESS lightcurves for RCB and the chemically similar, but dustless hydrogen-deficient carbon (dLHdC) stars. We find that dLHdC stars show variations on timescales shorter than RCB stars, suggesting that they may have lower masses than RCB stars. Finally, we identified 3 new spectroscopically confirmed and 12 candidate Galactic DY Per type stars - believed to be colder cousins of RCB stars - doubling the sample of Galactic DY Per type stars. △ Less

Submitted 11 July, 2024; originally announced July 2024.

Comments: accepted for publication in PASP

Unraveling the role of merger histories in the population of Insitu stars: linking IllustrisTNG cosmological simulation to H3 survey

Authors: Razieh Emami, Lars Hernquist, Randall Smith, James F. Steiner, Grant Tremblay, Douglas Finkbeiner, Mark Vogelsberger, Josh Grindlay, Federico Marinacci, Kung-Yi Su, Cecilia Garraffo, Yuan-Sen Ting, Phillip A. Cargile, Rebecca L. Davies, Chloë E. Benton, Yijia Li, Letizia Bugiani, Amir H. Khoram, Sownak Bose

Abstract: We undertake a comprehensive investigation into the distribution of insitu stars within Milky Way-like galaxies, leveraging TNG50 simulations and comparing their predictions with data from the H3 survey. Our analysis reveals that 28% of galaxies demonstrate reasonable agreement with H3, while only 12% exhibit excellent alignment in their profiles, regardless of the specific spatial cut employed to… ▽ More We undertake a comprehensive investigation into the distribution of insitu stars within Milky Way-like galaxies, leveraging TNG50 simulations and comparing their predictions with data from the H3 survey. Our analysis reveals that 28% of galaxies demonstrate reasonable agreement with H3, while only 12% exhibit excellent alignment in their profiles, regardless of the specific spatial cut employed to define insitu stars. To uncover the underlying factors contributing to deviations between TNG50 and H3 distributions, we scrutinize correlation coefficients among internal drivers(e.g., virial radius, star formation rate [SFR]) and merger-related parameters (such as the effective mass-ratio, mean distance, average redshift, total number of mergers, average spin-ratio and maximum spin alignment between merging galaxies). Notably, we identify significant correlations between deviations from observational data and key parameters such as the median slope of virial radius, mean SFR values, and the rate of SFR change across different redshift scans. Furthermore, positive correlations emerge between deviations from observational data and parameters related to galaxy mergers. We validate these correlations using the Random Forest Regression method. Our findings underscore the invaluable insights provided by the H3 survey in unraveling the cosmic history of galaxies akin to the Milky Way, thereby advancing our understanding of galactic evolution and shedding light on the formation and evolution of Milky Way-like galaxies in cosmological simulations. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Comments: 20 pages, 13 figures

The large-scale structure around the Fornax-Eridanus Complex

Authors: Maria Angela Raj, Petra Awad, Reynier F. Peletier, Rory Smith, Ulrike Kuchner, Rien van de Weygaert, Noam I. Libeskind, Marco Canducci, Peter Tino, Kerstin Bunte

Abstract: Our objectives are to map the filamentary network around the Fornax-Eridanus Complex and probe the influence of the local environment on galaxy morphology. We employ the novel machine-learning tool, 1-DREAM (1-Dimensional, Recovery, Extraction, and Analysis of Manifolds) to detect and model filaments around the Fornax cluster. We then use the morphology-density relation of galaxies to examine the… ▽ More Our objectives are to map the filamentary network around the Fornax-Eridanus Complex and probe the influence of the local environment on galaxy morphology. We employ the novel machine-learning tool, 1-DREAM (1-Dimensional, Recovery, Extraction, and Analysis of Manifolds) to detect and model filaments around the Fornax cluster. We then use the morphology-density relation of galaxies to examine the variation in the galaxies' morphology with respect to their distance from the central axis of the detected filaments. We detect 27 filaments that vary in length and galaxy-number density around the Fornax-Eridanus Complex. These filaments showcase a variety of environments; some filaments encompass groups/clusters, while others are only inhabited by galaxies in pristine filamentary environments. We also reveal a well-known structure -- the Fornax Wall, that passes through the Dorado group, Fornax cluster, and Eridanus supergroup. Regarding the morphology of galaxies, we find that early-type galaxies (ETGs) populate high-density filaments and high-density regions of the Fornax Wall. Furthermore, the fraction of ETGs decreases as the distance to the filament spine increases. Of the total galaxy population in filaments, ~7% are ETGs and ~24% are late-type galaxies (LTGs) located in pristine environments of filaments, while ~27% are ETGs and ~42% are LTGs in groups/clusters within filaments. This study reveals the Cosmic Web around the Fornax Cluster and asserts that filamentary environments are heterogeneous in nature. When investigating the role of the environment on galaxy morphology, it is essential to consider both, the local number-density and a galaxy's proximity to the filament spine. Within this framework, we ascribe the observed morphological segregation in the Fornax Wall to pre-processing of galaxies within groups embedded in it. △ Less

Submitted 3 July, 2024; originally announced July 2024.

Comments: Accepted for publication in A&A. 21 pages with 15 figures

Pop Stacks with a Bypass

Authors: Lapo Cioni, Luca Ferrari, Rebecca Smith

Abstract: We consider sorting procedures for permutations making use of pop stacks with a bypass operation, and explore the combinatorial properties of the associated algorithms. We consider sorting procedures for permutations making use of pop stacks with a bypass operation, and explore the combinatorial properties of the associated algorithms. △ Less

Submitted 24 June, 2024; originally announced June 2024.

Comments: In Proceedings GASCom 2024, arXiv:2406.14588

ACM Class: G.2.1

Journal ref: EPTCS 403, 2024, pp. 73-78

Non-stabilizerness in kinetically-constrained Rydberg atom arrays

Authors: Ryan Smith, Zlatko Papić, Andrew Hallam

Abstract: Non-stabilizer states are a fundamental resource for universal quantum computation. However,despite broad significance in quantum computing, the emergence of "many-body" non-stabilizerness in interacting quantum systems remains poorly understood due to its analytical intractability. Here we show that Rydberg atom arrays provide a natural reservoir of non-stabilizerness that extends beyond single q… ▽ More Non-stabilizer states are a fundamental resource for universal quantum computation. However,despite broad significance in quantum computing, the emergence of "many-body" non-stabilizerness in interacting quantum systems remains poorly understood due to its analytical intractability. Here we show that Rydberg atom arrays provide a natural reservoir of non-stabilizerness that extends beyond single qubits and arises from quantum correlations engendered by the Rydberg blockade. We demonstrate that this non-stabilizerness can be experimentally accessed in two complementary ways, either by performing quench dynamics or via adiabatic ground state preparation. Using the analytical framework based on matrix product states, we explain the origin of Rydberg nonstabilizerness via a quantum circuit decomposition of the wave function. △ Less

Submitted 20 June, 2024; originally announced June 2024.

Comments: 15 pages, 7 figures

The Blue Multi Unit Spectroscopic Explorer (BlueMUSE) on the VLT: science drivers and overview of instrument design

Authors: Johan Richard, Rémi Giroud, Florence Laurent, Davor Krajnović, Alexandre Jeanneau, Roland Bacon, Manuel Abreu, Angela Adamo, Ricardo Araujo, Nicolas Bouché, Jarle Brinchmann, Zhemin Cai, Norberto Castro, Ariadna Calcines, Diane Chapuis, Adélaïde Claeyssens, Luca Cortese, Emanuele Daddi, Christopher Davison, Michael Goodwin, Robert Harris, Matthew Hayes, Mathilde Jauzac, Andreas Kelz, Jean-Paul Kneib , et al. (24 additional authors not shown)

Abstract: BlueMUSE is a blue-optimised, medium spectral resolution, panoramic integral field spectrograph under development for the Very Large Telescope (VLT). With an optimised transmission down to 350 nm, spectral resolution of R$\sim$3500 on average across the wavelength range, and a large FoV (1 arcmin$^2$), BlueMUSE will open up a new range of galactic and extragalactic science cases facilitated by its… ▽ More BlueMUSE is a blue-optimised, medium spectral resolution, panoramic integral field spectrograph under development for the Very Large Telescope (VLT). With an optimised transmission down to 350 nm, spectral resolution of R$\sim$3500 on average across the wavelength range, and a large FoV (1 arcmin$^2$), BlueMUSE will open up a new range of galactic and extragalactic science cases facilitated by its specific capabilities. The BlueMUSE consortium includes 9 institutes located in 7 countries and is led by the Centre de Recherche Astrophysique de Lyon (CRAL). The BlueMUSE project development is currently in Phase A, with an expected first light at the VLT in 2031. We introduce here the Top Level Requirements (TLRs) derived from the main science cases, and then present an overview of the BlueMUSE system and its subsystems fulfilling these TLRs. We specifically emphasize the tradeoffs that are made and the key distinctions compared to the MUSE instrument, upon which the system architecture is built. △ Less

Submitted 19 June, 2024; originally announced June 2024.

Comments: 20 pages, 10 figures, proceedings of the SPIE astronomical telescopes and instrumentation conference, Yokohama, 16-21 June

The Design, Implementation, and Performance of the LZ Calibration Systems

Authors: J. Aalbers, D. S. Akerib, A. K. Al Musalhi, F. Alder, C. S. Amarasinghe, A. Ames, T. J. Anderson, N. Angelides, H. M. Araújo, J. E. Armstrong, M. Arthurs, A. Baker, S. Balashov, J. Bang, E. E. Barillier, J. W. Bargemann, K. Beattie, T. Benson, A. Bhatti, A. Biekert, T. P. Biesiadzinski, H. J. Birch, E. Bishop, G. M. Blockinger, B. Boxer , et al. (179 additional authors not shown)

Abstract: LUX-ZEPLIN (LZ) is a tonne-scale experiment searching for direct dark matter interactions and other rare events. It is located at the Sanford Underground Research Facility (SURF) in Lead, South Dakota, USA. The core of the LZ detector is a dual-phase xenon time projection chamber (TPC), designed with the primary goal of detecting Weakly Interacting Massive Particles (WIMPs) via their induced low e… ▽ More LUX-ZEPLIN (LZ) is a tonne-scale experiment searching for direct dark matter interactions and other rare events. It is located at the Sanford Underground Research Facility (SURF) in Lead, South Dakota, USA. The core of the LZ detector is a dual-phase xenon time projection chamber (TPC), designed with the primary goal of detecting Weakly Interacting Massive Particles (WIMPs) via their induced low energy nuclear recoils. Surrounding the TPC, two veto detectors immersed in an ultra-pure water tank enable reducing background events to enhance the discovery potential. Intricate calibration systems are purposely designed to precisely understand the responses of these three detector volumes to various types of particle interactions and to demonstrate LZ's ability to discriminate between signals and backgrounds. In this paper, we present a comprehensive discussion of the key features, requirements, and performance of the LZ calibration systems, which play a crucial role in enabling LZ's WIMP-search and its broad science program. The thorough description of these calibration systems, with an emphasis on their novel aspects, is valuable for future calibration efforts in direct dark matter and other rare-event search experiments. △ Less

Submitted 20 June, 2024; v1 submitted 2 May, 2024; originally announced June 2024.

ZTF SN Ia DR2: Exploring SN Ia properties in the vicinity of under-dense environments

Authors: M. Aubert, P. Rosnet, B. Popovic, F. Ruppin, M. Smith, M. Rigault, G. Dimitriadis, A. Goobar, J. Johansson, C. Barjou-Delayre, U. Burgaz, B. Carreres, F. Feinstein, D. Fouchez, L. Galbany, M. Ginolin, T. de Jaeger, M. M. Kasliwal, Y. -L. Kim, L. Lacroix, F. J. Masci, T. E. Müller-Bravo, B. Racine, C. Ravoux, N. Regnault , et al. (7 additional authors not shown)

Abstract: The unprecedented statistics of detected Type Ia supernovae (SNe Ia) brought by the Zwicky Transient Facility enables us to probe the impact of the Large-Scale Structure on the properties of these objects. The goal of this paper is to explore the possible impact of the under-dense part of the large-scale structure on the intrinsic SALT2 light curve properties of SNe Ia and uncover possible biases… ▽ More The unprecedented statistics of detected Type Ia supernovae (SNe Ia) brought by the Zwicky Transient Facility enables us to probe the impact of the Large-Scale Structure on the properties of these objects. The goal of this paper is to explore the possible impact of the under-dense part of the large-scale structure on the intrinsic SALT2 light curve properties of SNe Ia and uncover possible biases in SN Ia analyses. With a volume-limited selection of ZTF-Cosmo-DR2 Type Ia supernovae overlapping with the SDSS-DR7 survey footprint, we investigate the distribution of their properties with regard to voids detected in the SDSS-DR7 galaxy sample. We further use Voronoi volumes as proxy for local density environments within the large-scale structure. We find a moderate dependency of the stretch toward the localisation around the void centre and none when considering colour. The local Voronoi volumes mostly affect the fraction of low/high stretch supernovae. With the current statistics available, we consider that the impact of high or low local density environment can be considered as a proxy for the colour of the host galaxy. Under-dense environments should not cause any biases in supernova analyses. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: 10 pages, 8 figures. Submitted to A&A

Negative thermal expansion in transition-metal dicyanides: the hidden role of the underlying diamondoid framework

Authors: Quentin Guéroult, Jonathan Bulled, Henry Patteson, Chloe Coates, Ronald Smith, Helen Playford, David Keen, Andrew Goodwin

Abstract: The transition-metal dicyanides M(CN)$_2$ (M = Zn, Cd) are amongst the most important negative thermal expansion (NTE) materials known, favoured for the magnitude, isotropy, and thermal persistence of the NTE behaviour they show. The conventional picture of the NTE mechanism in this family is one of correlated rotations and translations of M(C/N)$_4$ polyhedra acting to draw the diamondoid network… ▽ More The transition-metal dicyanides M(CN)$_2$ (M = Zn, Cd) are amongst the most important negative thermal expansion (NTE) materials known, favoured for the magnitude, isotropy, and thermal persistence of the NTE behaviour they show. The conventional picture of the NTE mechanism in this family is one of correlated rotations and translations of M(C/N)$_4$ polyhedra acting to draw the diamondoid network of M--CN--M linkages in on itself. An implication of this mechanism is increased transverse vibrational motion of C and N atoms relative to the isotropic displacements of M atoms, which act as anchors. Here, we use a combination of neutron total scattering measurements and \emph{ab initio} calculations to reassess the vibrational behaviour of the M(CN)$_2$ family. We find that M, C, and N atoms all exhibit similar degrees of local thermal motion, such that the cyanide linkages behave as pseudo-springs connecting M$\ldots$M pairs. This interpretation leads us to uncover a `hidden' dispersion in the M(CN)$_2$ phonon dispersions, closely related to that of diamond and silicon themselves. By virtue of this mapping, a simple geometric model based on the competing energy scales of network stretching and flexing -- long applied to interpret NTE modes in C and Si -- turns out to capture the key NTE physics of M(CN)$_2$, especially at low temperatures. Our study highlights the potential insight gained by coarse-graining the complex lattice dynamics of framework materials in terms of what we call `framework modes' -- the correlated distortions of the underlying network structure itself. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Where to place a mosquito trap for West Nile Virus surveillance?

Authors: Anwesha Chakravarti, Bo Li, Dan Bartlett, Patrick Irwin, Rebecca Smith

Abstract: The rapid spread of West Nile Virus (WNV) is a growing concern. With no vaccines or specific medications available, prevention through mosquito control is the only solution to curb the spread. Mosquito traps, used to detect viral presence in mosquito populations, are essential tools for WNV surveillance. But how do we decide where to place a mosquito trap? And what makes a good trap location, anyw… ▽ More The rapid spread of West Nile Virus (WNV) is a growing concern. With no vaccines or specific medications available, prevention through mosquito control is the only solution to curb the spread. Mosquito traps, used to detect viral presence in mosquito populations, are essential tools for WNV surveillance. But how do we decide where to place a mosquito trap? And what makes a good trap location, anyway? We present a robust statistical approach to determine a mosquito trap's ability to predict human WNV cases in the Chicago metropolitan area and its suburbs. We then use this value to detect the landscape, demographic, and socioeconomic factors associated with a mosquito trap's predictive ability. This approach enables resource-limited mosquito control programs to identify better trap locations while reducing trap numbers to increase trap-based surveillance efficiency. The approach can also be applied to a wide range of different environmental surveillance programs. △ Less

Submitted 10 June, 2024; originally announced June 2024.

Comments: 22 pages, 9 figures

Equivalence of Polarizability and Circuit Models for Waveguide-Fed Metamaterial Elements

Authors: David R. Smith, Mohsen Sazegar, Insang Yoo

Abstract: A common variant of a metasurface antenna consists of an array of metamaterial elements coupled to a waveguide feed. The guided wave excites the metamaterial elements, coupling energy from the waveguide mode to radiation. Under appropriate conditions, each sub-wavelength metamaterial element can be modeled as a polarizable dipole, with the polarizability determined by an extraction procedure from… ▽ More A common variant of a metasurface antenna consists of an array of metamaterial elements coupled to a waveguide feed. The guided wave excites the metamaterial elements, coupling energy from the waveguide mode to radiation. Under appropriate conditions, each sub-wavelength metamaterial element can be modeled as a polarizable dipole, with the polarizability determined by an extraction procedure from the computed or measured waveguide scattering parameters. Here we establish the equivalence of this polarizability description of a metamaterial element with an equivalent circuit model, providing an additional tool for metasurface design that offers significant insight and a path towards efficiently modeling very large apertures. With this equivalence established, more complicated external circuits that include lumped elements and devices such as diodes and transistors can be integrated into the metamaterial element, which can then be transformed into an equivalent polarizability for modeling in the coupled dipole framework. We derive appropriate circuit models for several basic metamaterial elements, which provide direct relationships between the equivalent circuit parameters of an element and its effective polarizability. These expressions are confirmed using scattering parameters for several example structures obtained via full-wave simulations. △ Less

Submitted 10 June, 2024; originally announced June 2024.

Detection of New Galaxy Candidates at $z\ >$ 11 in the JADES Field Using JWST NIRCam

Authors: Priyanka Chakraborty, Arnab Sarkar, Scott Wolk, Benjamin Schneider, Nancy Brickhouse, Kenneth Lanzetta, Adam Foster, Randall Smith

Abstract: We report the detection of seven new galaxy candidates with redshift $z$ $>$ 11 within the JWST Advanced Deep Extragalactic Survey (JADES) GOODS-S and GOODS-N fields. These new candidates are detected through meticulous analysis of NIRCam photometry in eight filters spanning a wavelength range of 0.8-5.0 $μ$m. Photometric redshifts of these galaxy candidates are independently measured utilizing sp… ▽ More We report the detection of seven new galaxy candidates with redshift $z$ $>$ 11 within the JWST Advanced Deep Extragalactic Survey (JADES) GOODS-S and GOODS-N fields. These new candidates are detected through meticulous analysis of NIRCam photometry in eight filters spanning a wavelength range of 0.8-5.0 $μ$m. Photometric redshifts of these galaxy candidates are independently measured utilizing spectral energy distribution (SED) fitting techniques using \texttt{EAZY} and \texttt{BAGPIPES} codes, followed by visual scrutiny. Two of these galaxy candidates are located in GOODS-S field, while the remaining five galaxies are located in GOODS-N field. Our analysis reveals that the stellar masses of these galaxies typically range from log $M_{\ast}$/$M_{\odot}$ = 7.75--8.75. Futhermore, these galaxies are typically young with their mass-weighted ages spanning from 80 to 240 Myr. Their specific star formation rates (sSFR), quantified as $\log (\text{sSFR}/\text{Gyr}$), are measured to vary between $\sim 0.95$ to 1.46. These new galaxy candidates offer a robust sample for probing the physical properties of galaxies within the first few hundred Myr of the history of the Universe. We also analyze the relationship between star formation rate (SFR) and stellar mass ($M_\ast$) within our sample. Using linear regression, our analysis yields a slope of $0.71 \pm 0.12$, which we then compare with results from previous studies. Continued investigation through spectroscopic analysis using JWST/NIRSpec is needed to spectroscopically confirm these high-redshift galaxy candidates and investigate further into their physical properties. We plan to follow up on these candidates with future NIRSpec observations. △ Less

Submitted 16 June, 2024; v1 submitted 7 June, 2024; originally announced June 2024.

Comments: 12 pages, 10 figures, Submitted to MNRAS

Probing the Scalar WIMP-Pion Coupling with the first LUX-ZEPLIN data

Authors: J. Aalbers, D. S. Akerib, A. K. Al Musalhi, F. Alder, C. S. Amarasinghe, A. Ames, T. J. Anderson, N. Angelides, H. M. Araújo, J. E. Armstrong, M. Arthurs, A. Baker, S. Balashov, J. Bang, E. E. Barillier, J. W. Bargemann, K. Beattie, T. Benson, A. Bhatti, A. Biekert, T. P. Biesiadzinski, H. J. Birch, E. J. Bishop, G. M. Blockinger, B. Boxer , et al. (178 additional authors not shown)

Abstract: Weakly interacting massive particles (WIMPs) may interact with a virtual pion that is exchanged between nucleons. This interaction channel is important to consider in models where the spin-independent isoscalar channel is suppressed. Using data from the first science run of the LUX-ZEPLIN dark matter experiment, containing 60 live days of data in a 5.5~tonne fiducial mass of liquid xenon, we repor… ▽ More Weakly interacting massive particles (WIMPs) may interact with a virtual pion that is exchanged between nucleons. This interaction channel is important to consider in models where the spin-independent isoscalar channel is suppressed. Using data from the first science run of the LUX-ZEPLIN dark matter experiment, containing 60 live days of data in a 5.5~tonne fiducial mass of liquid xenon, we report the results on a search for WIMP-pion interactions. We observe no significant excess and set an upper limit of $1.5\times10^{-46}$~cm$^2$ at a 90\% confidence level for a WIMP mass of 33~GeV/c$^2$ for this interaction. △ Less

Submitted 4 June, 2024; originally announced June 2024.

The 2024 Brain Tumor Segmentation (BraTS) Challenge: Glioma Segmentation on Post-treatment MRI

Authors: Maria Correia de Verdier, Rachit Saluja, Louis Gagnon, Dominic LaBella, Ujjwall Baid, Nourel Hoda Tahon, Martha Foltyn-Dumitru, Jikai Zhang, Maram Alafif, Saif Baig, Ken Chang, Gennaro D'Anna, Lisa Deptula, Diviya Gupta, Muhammad Ammar Haider, Ali Hussain, Michael Iv, Marinos Kontzialis, Paul Manning, Farzan Moodi, Teresa Nunes, Aaron Simon, Nico Sollmann, David Vu, Maruf Adewole , et al. (60 additional authors not shown)

Abstract: Gliomas are the most common malignant primary brain tumors in adults and one of the deadliest types of cancer. There are many challenges in treatment and monitoring due to the genetic diversity and high intrinsic heterogeneity in appearance, shape, histology, and treatment response. Treatments include surgery, radiation, and systemic therapies, with magnetic resonance imaging (MRI) playing a key r… ▽ More Gliomas are the most common malignant primary brain tumors in adults and one of the deadliest types of cancer. There are many challenges in treatment and monitoring due to the genetic diversity and high intrinsic heterogeneity in appearance, shape, histology, and treatment response. Treatments include surgery, radiation, and systemic therapies, with magnetic resonance imaging (MRI) playing a key role in treatment planning and post-treatment longitudinal assessment. The 2024 Brain Tumor Segmentation (BraTS) challenge on post-treatment glioma MRI will provide a community standard and benchmark for state-of-the-art automated segmentation models based on the largest expert-annotated post-treatment glioma MRI dataset. Challenge competitors will develop automated segmentation models to predict four distinct tumor sub-regions consisting of enhancing tissue (ET), surrounding non-enhancing T2/fluid-attenuated inversion recovery (FLAIR) hyperintensity (SNFH), non-enhancing tumor core (NETC), and resection cavity (RC). Models will be evaluated on separate validation and test datasets using standardized performance metrics utilized across the BraTS 2024 cluster of challenges, including lesion-wise Dice Similarity Coefficient and Hausdorff Distance. Models developed during this challenge will advance the field of automated MRI segmentation and contribute to their integration into clinical practice, ultimately enhancing patient care. △ Less

Submitted 28 May, 2024; originally announced May 2024.

Comments: 10 pages, 4 figures, 1 table

Asymptotic Theory for Estimation of the Husler-Reiss Distribution via Block Maxima Method

Authors: Hank Flury, Jan Hannig, Richard Smith

Abstract: The Hüsler-Reiss distribution describes the limit of the pointwise maxima of a bivariate normal distribution. This distribution is defined by a single parameter, $λ$. We provide asymptotic theory for maximum likelihood estimation of $λ$ under a block maxima approach. Our work assumes independent and identically distributed bivariate normal random variables, grouped into blocks where the block size… ▽ More The Hüsler-Reiss distribution describes the limit of the pointwise maxima of a bivariate normal distribution. This distribution is defined by a single parameter, $λ$. We provide asymptotic theory for maximum likelihood estimation of $λ$ under a block maxima approach. Our work assumes independent and identically distributed bivariate normal random variables, grouped into blocks where the block size and number of blocks increase simultaneously. With these assumptions our results provide conditions for the asymptotic normality of the Maximum Likelihood Estimator (MLE). We characterize the bias of the MLE, provide conditions under which this bias is asymptotically negligible, and discuss how to choose the block size to minimize a bias-variance trade-off. The proofs are an extension of previous results for choosing the block size in the estimation of univariate extreme value distributions (Dombry and Ferreria 2019), providing a potential basis for extensions to multivariate cases where both the marginal and dependence parameters are unknown. The proofs rely on the Argmax Theorem applied to a localized loglikelihood function, combined with a Lindeberg-Feller Central Limit Theorem argument to establish asymptotic normality. Possible applications of the method include composite likelihood estimation in Brown-Resnick processes, where it is known that the bivariate distributions are of Hüsler-Reiss form. △ Less

Submitted 24 May, 2024; originally announced May 2024.

The Data Acquisition System of the LZ Dark Matter Detector: FADR

Authors: J. Aalbers, D. S. Akerib, A. K. Al Musalhi, F. Alder, C. S. Amarasinghe, A. Ames, T. J. Anderson, N. Angelides, H. M. Araújo, J. E. Armstrong, M. Arthurs, A. Baker, S. Balashov, J. Bang, E. E. Barillier, J. W. Bargemann, K. Beattie, T. Benson, A. Bhatti, A. Biekert, T. P. Biesiadzinski, H. J. Birch, E. Bishop, G. M. Blockinger, B. Boxer , et al. (190 additional authors not shown)

Abstract: The Data Acquisition System (DAQ) for the LUX-ZEPLIN (LZ) dark matter detector is described. The signals from 745 PMTs, distributed across three subsystems, are sampled with 100-MHz 32-channel digitizers (DDC-32s). A basic waveform analysis is carried out on the on-board Field Programmable Gate Arrays (FPGAs) to extract information about the observed scintillation and electroluminescence signals.… ▽ More The Data Acquisition System (DAQ) for the LUX-ZEPLIN (LZ) dark matter detector is described. The signals from 745 PMTs, distributed across three subsystems, are sampled with 100-MHz 32-channel digitizers (DDC-32s). A basic waveform analysis is carried out on the on-board Field Programmable Gate Arrays (FPGAs) to extract information about the observed scintillation and electroluminescence signals. This information is used to determine if the digitized waveforms should be preserved for offline analysis. The system is designed around the Kintex-7 FPGA. In addition to digitizing the PMT signals and providing basic event selection in real time, the flexibility provided by the use of FPGAs allows us to monitor the performance of the detector and the DAQ in parallel to normal data acquisition. The hardware and software/firmware of this FPGA-based Architecture for Data acquisition and Realtime monitoring (FADR) are discussed and performance measurements are described. △ Less

Submitted 23 May, 2024; originally announced May 2024.

Comments: 18 pages, 24 figures

Is the EJRA proportionate and therefore justified? A critical review of the EJRA policy at Cambridge

Authors: Oliver Linton, Raghavendra Rau, Patrick Baert, Peter Bossaerts, Jon Crowcroft, G. R. Evans, Paul Ewart, Nick Gay, Paul Kattuman, Stefan Scholtes, Hamid Sabourian, Richard J. Smith

Abstract: This paper critically evaluates the HESA (Higher Education Statistics Agency) Data Report for the Employer Justified Retirement Age (EJRA) Review Group at the University of Cambridge (\cite{CambridgeHESA2024}), identifying significant methodological flaws and misinterpretations. Our analysis reveals issues such as unclear application of data filters, inconsistent variable treatment, and erroneous… ▽ More This paper critically evaluates the HESA (Higher Education Statistics Agency) Data Report for the Employer Justified Retirement Age (EJRA) Review Group at the University of Cambridge (\cite{CambridgeHESA2024}), identifying significant methodological flaws and misinterpretations. Our analysis reveals issues such as unclear application of data filters, inconsistent variable treatment, and erroneous statistical conclusions. The Report suggests that the EJRA increased job creation rates at Cambridge, but we show Cambridge consistently had lower job creation rates for Established Academic Careers compared to other Russell Group universities, both before and after EJRA implementation in 2011, with no evidence for a significant change in this deficit post implementation. This suggests that EJRA is not a significant factor driving job creation rates. Since other universities without an EJRA exhibit higher job creation rates, this suggests job creation can be sustained without such a policy. We conclude that the EJRA did not achieve its intended goal of increasing opportunities for young academics and may have exacerbated existing disparities compared to other leading universities. We recommend EJRA be abolished at Cambridge since it does not meet its aims and could be viewed as unlawful age discrimination.\newline \it{This version is a revision reflecting some of the comments made by members of the university EJRA review group in public discussion, see \cite{Holmes24}.} △ Less

Submitted 12 June, 2024; v1 submitted 23 May, 2024; originally announced May 2024.

Comments: This is a revision of our earlier paper

Euclid. I. Overview of the Euclid mission

Authors: Euclid Collaboration, Y. Mellier, Abdurro'uf, J. A. Acevedo Barroso, A. Achúcarro, J. Adamek, R. Adam, G. E. Addison, N. Aghanim, M. Aguena, V. Ajani, Y. Akrami, A. Al-Bahlawan, A. Alavi, I. S. Albuquerque, G. Alestas, G. Alguero, A. Allaoui, S. W. Allen, V. Allevato, A. V. Alonso-Tetilla, B. Altieri, A. Alvarez-Candal, A. Amara, L. Amendola , et al. (1086 additional authors not shown)

Abstract: The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14… ▽ More The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue`Euclid on Sky'

Sub-relativistic Outflow and Hours-Timescale Large-amplitude X-ray Dips during Super-Eddington Accretion onto a Low-mass Massive Black Hole in the Tidal Disruption Event AT2022lri

Authors: Yuhan Yao, Muryel Guolo, Francesco Tombesi, Ruancun Li, Suvi Gezari, Javier A. García, Lixin Dai, Ryan Chornock, Wenbin Lu, S. R. Kulkarni, Keith C. Gendreau, Dheeraj R. Pasham, S. Bradley Cenko, Erin Kara, Raffaella Margutti, Yukta Ajay, Thomas Wevers, Tom M. Kwan, Igor Andreoni, Joshua S. Bloom, Andrew J. Drake, Matthew J. Graham, Erica Hammerstein, Russ R. Laher, Natalie LeBaron , et al. (10 additional authors not shown)

Abstract: We present the tidal disruption event (TDE) AT2022lri, hosted in a nearby ($\approx\!144$ Mpc) quiescent galaxy with a low-mass massive black hole ($10^4\,M_\odot < M_{\rm BH} < 10^6\,M_\odot$). AT2022lri belongs to the TDE-H+He subtype. More than 1 Ms of X-ray data were collected with NICER, Swift, and XMM-Newton from 187 d to 672 d after peak. The X-ray luminosity gradually declined from… ▽ More We present the tidal disruption event (TDE) AT2022lri, hosted in a nearby ($\approx\!144$ Mpc) quiescent galaxy with a low-mass massive black hole ($10^4\,M_\odot < M_{\rm BH} < 10^6\,M_\odot$). AT2022lri belongs to the TDE-H+He subtype. More than 1 Ms of X-ray data were collected with NICER, Swift, and XMM-Newton from 187 d to 672 d after peak. The X-ray luminosity gradually declined from $1.5\times 10^{44}\,{\rm erg\,s^{-1}}$ to $1.5\times 10^{43}\,{\rm erg\,s^{-1}}$ and remains much above the UV and optical luminosity, consistent with a super-Eddington accretion flow viewed face-on. Sporadic strong X-ray dips atop a long-term decline are observed, with variability timescale of $\approx\!0.5$ hr--1 d and amplitude of $\approx\!2$--8. When fitted with simple continuum models, the X-ray spectrum is dominated by a thermal disk component with inner temperature going from $\sim\! 146$ eV to $\sim\! 86$ eV. However, there are residual features that peak around 1 keV, which, in some cases, cannot be reproduced by a single broad emission line. We analyzed a subset of time-resolved spectra with two physically motivated models describing either a scenario where ionized absorbers contribute extra absorption and emission lines or where disk reflection plays an important role. Both models provide good and statistically comparable fits, show that the X-ray dips are correlated with drops in the inner disk temperature, and require the existence of sub-relativistic (0.1--0.3$c$) ionized outflows. We propose that the disk temperature fluctuation stems from episodic drops of the mass accretion rate triggered by magnetic instabilities or/and wobbling of the inner accretion disk along the black hole's spin axis. △ Less

Submitted 18 May, 2024; originally announced May 2024.

Comments: 35 pages, 20 figures, submitted

Self-consistent modelling of the Milky Way structure using live potentials

Authors: Eva Durán-Camacho, Ana Duarte-Cabral, Alex R. Pettitt, Robin G. Treß, Paul C. Clark, Ralf S. Klessen, Kamran R. J. Bogue, Rowan J. Smith, Mattia C. Sormani

Abstract: To advance our understanding of the evolution of the interstellar medium (ISM) of our Galaxy, numerical models of Milky Way (MW) type galaxies are widely used. However, most models only vaguely resemble the MW (e.g. in total mass), and often use imposed analytic potentials (which cannot evolve dynamically). This poses a problem in asserting their applicability for the interpretation of observation… ▽ More To advance our understanding of the evolution of the interstellar medium (ISM) of our Galaxy, numerical models of Milky Way (MW) type galaxies are widely used. However, most models only vaguely resemble the MW (e.g. in total mass), and often use imposed analytic potentials (which cannot evolve dynamically). This poses a problem in asserting their applicability for the interpretation of observations of our own Galaxy. The goal of this work is to identify a numerical model that is not only a MW-type galaxy, but one that can mimic some of the main observed structures of our Galaxy, using dynamically evolving potentials, so that it can be used as a base model to study the ISM cycle in a galaxy like our own. This paper introduces a suite of 15 MW-type galaxy models developed using the {\sc arepo} numerical code, that are compared to Galactic observations of $^{12}$CO and \ion{H}{I} emission via longitude-velocity plots, from where we extract and compare the skeletons of major galactic features and the terminal gas velocities. We found that our best-fitting model to the overall structure, also reproduces some of the more specific observed features of the MW, including a bar with a pattern speed of $30.0 \pm 0.2$ km\,s$^{-1}$\,kpc$^{-1}$, a bar half-length of $3.2 \pm 0.8$\,kpc. Our model shows large streaming motions around spiral arms, and strong radial motions well beyond the inner bar. This model highlights the complex motions of a dynamic MW-type galaxy and has the potential to offer valuable insight into how our Galaxy regulates the ISM and star formation. △ Less

Submitted 11 June, 2024; v1 submitted 15 May, 2024; originally announced May 2024.

Comments: Accepted for publication in MNRAS. 24 pages, 23 figures, 3 tables

Formation channels of the diffuse lights in the groups and clusters over time

Authors: Kyungwon Chun, Jihye Shin, Jongwan Ko, Rory Smith, Jaewon Yoo

Abstract: We explore the formation of the intragroup light (IGL) and intracluster light (ICL), representing diffuse lights within groups and clusters, since $z=1.5$. For this, we perform multi-resolution cosmological N-body simulations using the ``galaxy replacement technique" (GRT) and identify the progenitors in which the diffuse light stars existed when they fell into the groups or clusters. Our findings… ▽ More We explore the formation of the intragroup light (IGL) and intracluster light (ICL), representing diffuse lights within groups and clusters, since $z=1.5$. For this, we perform multi-resolution cosmological N-body simulations using the ``galaxy replacement technique" (GRT) and identify the progenitors in which the diffuse light stars existed when they fell into the groups or clusters. Our findings reveal that typical progenitors contributing to diffuse lights enter the host halo with the massive galaxies containing a stellar mass of $10 < \log M_{\rm{gal}}~[M_{\odot}]< 11$, regardless of the mass or dynamical state of the host halos at $z=0$. In cases where the host halos are dynamically unrelaxed or more massive, diffuse lights from massive progenitors with $\log M_{\rm{gal}}~[M_{\odot}]> 11$ are more prominent, with over half of them already pre-processed before entering the host halo. Additionally, we find that the main formation mechanism of diffuse lights is the stripping process of satellites, and a substantial fraction ($40-45\%$) of diffuse light stars is linked to the merger tree of the BCG. Remarkably, all trends persist for groups and clusters at higher redshifts. The fraction of diffuse lights in the host halos with a similar mass decreases as the redshift increases, but they are already substantial at $z=1.5$ ($\sim10\%$). However, it's crucial to acknowledge that detection limits related to the observable radius and faint-end surface brightness may obscure numerous diffuse light stars and even alter the main formation channel of diffuse lights. △ Less

Submitted 13 May, 2024; originally announced May 2024.

Comments: 18 pages, 13 figures. Accepted for publication in The Astrophysical Journal, May 10th, 2024

Large deviations in statistics of the local time and occupation time for a run and tumble particle

Authors: Soheli Mukherjee, Pierre Le Doussal, Naftali R. Smith

Abstract: We investigate the statistics of the local time $\mathcal{T} = \int_0^T δ(x(t)) dt$ that a run and tumble particle (RTP) $x(t)$ in one dimension spends at the origin, with or without an external drift. By relating the local time to the number of times the RTP crosses the origin, we find that the local time distribution $P(\mathcal{T})$ satisfies the large deviation principle… ▽ More We investigate the statistics of the local time $\mathcal{T} = \int_0^T δ(x(t)) dt$ that a run and tumble particle (RTP) $x(t)$ in one dimension spends at the origin, with or without an external drift. By relating the local time to the number of times the RTP crosses the origin, we find that the local time distribution $P(\mathcal{T})$ satisfies the large deviation principle $P(\mathcal{T}) \sim \, e^{-T \, I(\mathcal{T} / T)} $ in the large observation time limit $T \to \infty$. Remarkably, we find that in presence of drift the rate function $I(ρ)$ is nonanalytic: We interpret its singularity as dynamical phase transitions of first order. We then extend these results by studying the statistics of the amount of time $\mathcal{R}$ that the RTP spends inside a finite interval (i.e., the occupation time), with qualitatively similar results. In particular, this yields the long-time decay rate of the probability $P(\mathcal{R} = T)$ that the particle does not exit the interval up to time $T$. We find that the conditional endpoint distribution exhibits an interesting change of behavior from unimodal to bimodal as a function of the size of the interval. To study the occupation time statistics, we extend the Donsker-Varadhan large-deviation formalism to the case of RTPs, for general dynamical observables and possibly in the presence of an external potential. △ Less

Submitted 11 May, 2024; originally announced May 2024.

Comments: 17 pages, 5 figures

New observational recipes for measuring dynamical state of galaxy clusters

Authors: Hyowon Kim, Rory Smith, Jongwan Ko, Jong-Ho Shinn, Kyungwon Chun, Jihye Shin, Jaewon Yoo

Abstract: During cluster assembly, a cluster's virialization process leaves behind signatures that can provide information on its dynamical state. However, no clear consensus yet exists on the best way to achieve this. Therefore, we attempt to derive improved recipes for classifying the dynamical state of clusters in observations using cosmological simulations. Cluster halo mass and their subhalos' mass are… ▽ More During cluster assembly, a cluster's virialization process leaves behind signatures that can provide information on its dynamical state. However, no clear consensus yet exists on the best way to achieve this. Therefore, we attempt to derive improved recipes for classifying the dynamical state of clusters in observations using cosmological simulations. Cluster halo mass and their subhalos' mass are used to $ 10^{14}M_{\odot} h^{-1}$ and $10^{10}M_{\odot} h^{-1}$ to calculate five independent dynamical state indicators. We experiment with recipes by combining two to four indicators for detecting specific merger stages like recent and ancient mergers. These recipes are made by plotting merging clusters and a control sample of relaxed clusters in multiple indicators parameter space, and then applying a rotation matrix method to derive the best way to separate mergers from the control sample. The success of the recipe is quantified using the success rate and the overlap percentage of the merger and control histograms along the newly rotated $x$-axis. This provides us with recipes using different numbers of combined indicators and for different merger stage. Among the recipes, the stellar mass gap and center offset are the first and second most dominant of the indicators, and using more indicators improves the effectiveness of the recipe. When applied to observations, our results show good agreement with literature values of cluster dynamical state. △ Less

Submitted 10 May, 2024; originally announced May 2024.

Comments: 17pages, 9 figures, 3 tables, ApJ Accepted

Do spiral arms enhance star formation efficiency?

Authors: Miguel Querejeta, Adam K. Leroy, Sharon E. Meidt, Eva Schinnerer, Francesco Belfiore, Eric Emsellem, Ralf S. Klessen, Jiayi Sun, Mattia Sormani, Ivana Bešlic, Yixian Cao, Mélanie Chevance, Dario Colombo, Daniel A. Dale, Santiago García-Burillo, Simon C. O. Glover, Kathryn Grasha, Brent Groves, Eric. W. Koch, Lukas Neumann, Hsi-An Pan, Ismael Pessa, Jérôme Pety, Francesca Pinna, Lise Ramambason , et al. (10 additional authors not shown)

Abstract: Spiral arms are some of the most spectacular features in disc galaxies, and also present in our own Milky Way. It has been argued that star formation should proceed more efficiently in spiral arms as a result of gas compression. Yet, observational studies have so far yielded contradictory results. Here we examine arm/interarm surface density contrasts at ~100 pc resolution in 28 spiral galaxies fr… ▽ More Spiral arms are some of the most spectacular features in disc galaxies, and also present in our own Milky Way. It has been argued that star formation should proceed more efficiently in spiral arms as a result of gas compression. Yet, observational studies have so far yielded contradictory results. Here we examine arm/interarm surface density contrasts at ~100 pc resolution in 28 spiral galaxies from the PHANGS survey. We find that the arm/interarm contrast in stellar mass surface density (Sigma_*) is very modest, typically a few tens of percent. This is much smaller than the contrasts measured for molecular gas (Sigma_mol) or star formation rate (Sigma_SFR) surface density, which typically reach a factor of ~2-3. Yet, Sigma_mol and Sigma_SFR contrasts show a significant correlation with the enhancement in Sigma_*, suggesting that the small stellar contrast largely dictates the stronger accumulation of gas and star formation. All these contrasts increase for grand-design spirals compared to multi-armed and flocculent systems (and for galaxies with high stellar mass). The median star formation efficiency (SFE) of the molecular gas is 16% higher in spiral arms than in interarm regions, with a large scatter, and the contrast increases significantly (median SFE contrast 2.34) for regions of particularly enhanced stellar contrast (Sigma_* contrast >1.97). The molecular-to-atomic gas ratio (Sigma_mol/Sigma_atom) is higher in spiral arms, pointing to a transformation of atomic to molecular gas. In conclusion, the boost in the star formation efficiency of molecular gas in spiral arms is generally modest or absent, except for locations with exceptionally large stellar contrasts. (abridged) △ Less

Submitted 8 May, 2024; originally announced May 2024.

Comments: 26 pages, 16 figures. Accepted for publication in A&A

Identifying and exploiting alpha in linear asset pricing models with strong, semi-strong, and latent factors

Authors: Ron P. Smith, M. Hashem Pesaran

Abstract: The risk premia of traded factors are the sum of factor means and a parameter vector we denote by φ which is identified from the cross section regression of alpha of individual securities on the vector of factor loadings. If phi is non-zero one can construct "phi-portfolios" which exploit the systematic components of non-zero alpha. We show that for known values of betas and when phi is non-zero t… ▽ More The risk premia of traded factors are the sum of factor means and a parameter vector we denote by φ which is identified from the cross section regression of alpha of individual securities on the vector of factor loadings. If phi is non-zero one can construct "phi-portfolios" which exploit the systematic components of non-zero alpha. We show that for known values of betas and when phi is non-zero there exist phi-portfolios that dominate mean-variance portfolios. The paper then proposes a two-step bias corrected estimator of phi and derives its asymptotic distribution allowing for idiosyncratic pricing errors, weak missing factors, and weak error cross-sectional dependence. Small sample results from extensive Monte Carlo experiments show that the proposed estimator has the correct size with good power properties. The paper also provides an empirical application to a large number of U.S. securities with risk factors selected from a large number of potential risk factors according to their strength and constructs phi-portfolios and compares their Sharpe ratios to mean variance and S&P 500 portfolio. △ Less

Submitted 21 May, 2024; v1 submitted 3 May, 2024; originally announced May 2024.

Comments: The paper is composed of the main paper with mathematical appendices, plus two online supplements

Diffuse scattering from dynamically compressed single-crystal zirconium following the pressure-induced $α\toω$ phase transition

Authors: P. G. Heighway, S. Singh, M. G. Gorman, D. McGonegle, J. H. Eggert, R. F. Smith

Abstract: The prototypical $α\toω$ phase transition in zirconium is an ideal test-bed for our understanding of polymorphism under extreme loading conditions. After half a century of study, a consensus had emerged that the transition is realized via one of two distinct displacive mechanisms, depending on the nature of the compression path. However, recent dynamic-compression experiments equipped with in situ… ▽ More The prototypical $α\toω$ phase transition in zirconium is an ideal test-bed for our understanding of polymorphism under extreme loading conditions. After half a century of study, a consensus had emerged that the transition is realized via one of two distinct displacive mechanisms, depending on the nature of the compression path. However, recent dynamic-compression experiments equipped with in situ diffraction diagnostics performed in the past few years have revealed new transition mechanisms, demonstrating that our understanding of the underlying atomistic dynamics and transition kinetics is in fact far from complete. We present classical molecular dynamics simulations of the $α\toω$ phase transition in single-crystal zirconium shock-compressed along the [0001] axis using a machine-learning-class potential. The transition is predicted to proceed primarily via a modified version of the two-stage Usikov-Zilberstein mechanism, whereby the high-pressure $ω$-phase heterogeneously nucleates at boundaries between grains of an intermediate $β$-phase. We further observe the fomentation of atomistic disorder at the junctions between $β$ grains, leading to the formation of highly defective interstitial material between the $ω$ grains. We directly compare synthetic x-ray diffraction patterns generated from our simulations with those obtained using femtosecond diffraction in recent dynamic-compression experiments, and show that the simulations produce the same unique, anisotropic diffuse scattering signal unlike any previously seen from an elemental metal. Our simulations suggest that the diffuse signal arises from a combination of thermal diffuse scattering, nanoparticle-like scattering from residual kinetically stabilized $α$ and $β$ grains, and scattering from interstitial defective structures. △ Less

Submitted 29 April, 2024; originally announced April 2024.

Comments: 20 pages, 17 figures in main article; 5 pages, 5 figures in supplementary material. This article will be submitted to Physical Review B

Constraints On Covariant WIMP-Nucleon Effective Field Theory Interactions from the First Science Run of the LUX-ZEPLIN Experiment

Authors: J. Aalbers, D. S. Akerib, A. K. Al Musalhi, F. Alder, C. S. Amarasinghe, A. Ames, T. J. Anderson, N. Angelides, H. M. Araújo, J. E. Armstrong, M. Arthurs, A. Baker, S. Balashov, J. Bang, E. E. Barillier, J. W. Bargemann, K. Beattie, T. Benson, A. Bhatti, A. Biekert, T. P. Biesiadzinski, H. J. Birch, E. J. Bishop, G. M. Blockinger, B. Boxer , et al. (179 additional authors not shown)

Abstract: The first science run of the LUX-ZEPLIN (LZ) experiment, a dual-phase xenon time project chamber operating in the Sanford Underground Research Facility in South Dakota, USA, has reported leading limits on spin-independent WIMP-nucleon interactions and interactions described from a non-relativistic effective field theory (NREFT). Using the same 5.5~t fiducial mass and 60 live days of exposure we re… ▽ More The first science run of the LUX-ZEPLIN (LZ) experiment, a dual-phase xenon time project chamber operating in the Sanford Underground Research Facility in South Dakota, USA, has reported leading limits on spin-independent WIMP-nucleon interactions and interactions described from a non-relativistic effective field theory (NREFT). Using the same 5.5~t fiducial mass and 60 live days of exposure we report on the results of a relativistic extension to the NREFT. We present constraints on couplings from covariant interactions arising from the coupling of vector, axial currents, and electric dipole moments of the nucleon to the magnetic and electric dipole moments of the WIMP which cannot be described by recasting previous results described by an NREFT. Using a profile-likelihood ratio analysis, in an energy region between 0~keV$_\text{nr}$ to 270~keV$_\text{nr}$, we report 90% confidence level exclusion limits on the coupling strength of five interactions in both the isoscalar and isovector bases. △ Less

Submitted 26 April, 2024; originally announced April 2024.

Comments: 7 pages, 4 figures

Are early-type galaxies quenched by present-day environment? A study of dwarfs in the Fornax Cluster

Authors: Romero-Gómez, J., Reynier F. Peletier, J. A. L. Aguerri, R. Smith

Abstract: Galaxies undergo processes throughout their lifetimes that ultimately lead to the expulsion of the gas and the cessation of the star-forming activity. This phenomenon commonly known as quenching, can be caused by environmental processes. For this we use the results of Romero-Gómez et al. (2024), who analyzed galaxies from the SAMI-Fornax and ATLAS$^{3D}$ survey. Using t$_{90}$ as an approximation… ▽ More Galaxies undergo processes throughout their lifetimes that ultimately lead to the expulsion of the gas and the cessation of the star-forming activity. This phenomenon commonly known as quenching, can be caused by environmental processes. For this we use the results of Romero-Gómez et al. (2024), who analyzed galaxies from the SAMI-Fornax and ATLAS$^{3D}$ survey. Using t$_{90}$ as an approximation for the quenching time and comparing it with the infall time derived from phase-space models, we determine the probability of the quenching being produced by the local environment of galaxies. Our results reveal a relation between galaxy mass and quenching probability. Down to M$_{\star}$ $\sim$10$^{10}$ M$_{\odot}$, galaxies exhibit almost zero probability of quenching, suggesting their independence from environmental effects. As we move into the mass regime of dwarf galaxies, the probability increases with decreasing mass, highlighting their sensitivity to environmental quenching. For the dwarfs, 10$^{7}$ - 10$^{9}$ M$_{\odot}$, 36$\pm$9% of our observational data are consistent with this hypothesis, challenging the idea that the present-day cluster, Fornax, is the primary driver of quenching in the low mass galaxies. We compare these results with cosmological simulations, selecting galaxies under similar conditions to our observational sample. The simulated sample shows lower quenching probabilities as we move down in mass, only 5$\pm$1% of galaxies meet the quenching criteria. This discrepancy between observations and simulations underlines that modelling quenching is still in its infancy. In general, the number of observed galaxies quenched by their environment is lower than expected, which suggests that pre-processing plays a larger role in galaxy evolution. Ultimately, our results highlight the need for higher-quality simulations and refinement of galaxy formation and evolution models. △ Less

Submitted 23 April, 2024; originally announced April 2024.

Comments: Accepted for publication in A&A. 10 pages, 4 figures

On de Sitter vacua in O(d,d) invariant cosmology

Authors: Yang Liu, Antonio Padilla, Paul M. Saffin, Robert G. C. Smith

Abstract: We perform a thorough analysis of de Sitter vacua in O(d,d) invariant cosmologies. Starting with a homogeneous and isotropic framework we examine conditions for the existence of such vacua, non-perturbative in α' in both the string frame and the Einstein frame. We elucidate the nature of the instability in the string frame vacuum. For the Einstein frame, we demonstrate that the de Sitter vacuum ca… ▽ More We perform a thorough analysis of de Sitter vacua in O(d,d) invariant cosmologies. Starting with a homogeneous and isotropic framework we examine conditions for the existence of such vacua, non-perturbative in α' in both the string frame and the Einstein frame. We elucidate the nature of the instability in the string frame vacuum. For the Einstein frame, we demonstrate that the de Sitter vacuum cannot be eternal. We then extend our analysis to include Bianchi I universes where the O(d,d) symmetry includes scale factor exchange as well as scale factor duality. We show how the theory can be extended to the anisotropic case so that it admits de Sitter vacua, noting the crucial role played by the O(d,d) symmetry in satisfying any additional constraints. △ Less

Submitted 7 May, 2024; v1 submitted 23 April, 2024; originally announced April 2024.

Comments: Paper updated with references added, most notably that of Rost and Nunez which includes details of additional dS solutions

SPHEREx: NASA's Near-Infrared Spectrophotmetric All-Sky Survey

Authors: Brendan P. Crill, Michael Werner, Rachel Akeson, Matthew Ashby, Lindsey Bleem, James J. Bock, Sean Bryan, Jill Burnham, Joyce Byunh, Tzu-Ching Chang, Yi-Kuan Chiang, Walter Cook, Asantha Cooray, Andrew Davis, Olivier Doré, C. Darren Dowell, Gregory Dubois-Felsmann, Tim Eifler, Andreas Faisst, Salman Habib, Chen Heinrich, Katrin Heitmann, Grigory Heaton, Christopher Hirata, Viktor Hristov , et al. (29 additional authors not shown)

Abstract: SPHEREx, the Spectro-Photometer for the History of the Universe, Epoch of Reionization, and ices Explorer, is a NASA MIDEX mission planned for launch in 2024. SPHEREx will carry out the first all-sky spectral survey at wavelengths between 0.75 micron and 5 micron with spectral resolving power ~40 between 0.75 and 3.8 micron and ~120 between 3.8 and 5 micron At the end of its two-year mission, SPHE… ▽ More SPHEREx, the Spectro-Photometer for the History of the Universe, Epoch of Reionization, and ices Explorer, is a NASA MIDEX mission planned for launch in 2024. SPHEREx will carry out the first all-sky spectral survey at wavelengths between 0.75 micron and 5 micron with spectral resolving power ~40 between 0.75 and 3.8 micron and ~120 between 3.8 and 5 micron At the end of its two-year mission, SPHEREx will provide 0.75-to-5 micron spectra of each 6.2"x6.2" pixel on the sky - 14 billion spectra in all. This paper updates an earlier description of SPHEREx presenting changes made during the mission's Preliminary Design Phase, including a discussion of instrument integration and test and a summary of the data processing, analysis, and distribution plans. △ Less

Submitted 16 April, 2024; originally announced April 2024.

Journal ref: Proceedings Volume 11443, Space Telescopes and Instrumentation 2020: Optical, Infrared, and Millimeter Wave; 114430I (2020)

Peak Time-Windowed Risk Estimation of Stochastic Processes

Authors: Jared Miller, Niklas Schmid, Matteo Tacchi, Didier Henrion, Roy S. Smith

Abstract: This paper develops a method to upper-bound extreme-values of time-windowed risks for stochastic processes. Examples of such risks include the maximum average or 90% quantile of the current along a transmission line in any 5-minute window. This work casts the time-windowed risk analysis problem as an infinite-dimensional linear program in occupation measures. In particular, we employ the coherent… ▽ More This paper develops a method to upper-bound extreme-values of time-windowed risks for stochastic processes. Examples of such risks include the maximum average or 90% quantile of the current along a transmission line in any 5-minute window. This work casts the time-windowed risk analysis problem as an infinite-dimensional linear program in occupation measures. In particular, we employ the coherent risk measures of the mean and the expected shortfall (conditional value at risk) to define the maximal time-windowed risk along trajectories. The infinite-dimensional linear program must then be truncated into finite-dimensional optimization problems, such as by using the moment-sum of squares hierarchy of semidefinite programs. The infinite-dimensional linear program will have the same optimal value as the original nonconvex risk estimation task under compactness and regularity assumptions, and the sequence of semidefinite programs will converge to the true value under additional properties of algebraic characterization. The scheme is demonstrated for risk analysis of example stochastic processes. △ Less

Submitted 11 April, 2024; v1 submitted 10 April, 2024; originally announced April 2024.

Comments: 26 pages, 11 figures

Emergent Braitenberg-style Behaviours for Navigating the ViZDoom `My Way Home' Labyrinth

Authors: Caleidgh Bayer, Robert J. Smith, Malcolm I. Heywood

Abstract: The navigation of complex labyrinths with tens of rooms under visual partially observable state is typically addressed using recurrent deep reinforcement learning architectures. In this work, we show that navigation can be achieved through the emergent evolution of a simple Braitentberg-style heuristic that structures the interaction between agent and labyrinth, i.e. complex behaviour from simple… ▽ More The navigation of complex labyrinths with tens of rooms under visual partially observable state is typically addressed using recurrent deep reinforcement learning architectures. In this work, we show that navigation can be achieved through the emergent evolution of a simple Braitentberg-style heuristic that structures the interaction between agent and labyrinth, i.e. complex behaviour from simple heuristics. To do so, the approach of tangled program graphs is assumed in which programs cooperatively coevolve to develop a modular indexing scheme that only employs 0.8\% of the state space. We attribute this simplicity to several biases implicit in the representation, such as the use of pixel indexing as opposed to deploying a convolutional kernel or image processing operators. △ Less

Submitted 9 April, 2024; originally announced April 2024.

Dark matter free dwarf galaxy formation at the the tips of the tentacles of jellyfish galaxies

Authors: V. Lora, R. Smith, J. Fritz, A. Pasquali, A. C. Raga

Abstract: When falling into a galaxy cluster, galaxies experience a loss of gas due to ram pressure stripping. In particular, disk galaxies lose gas from their disks and very large tentacles of gas can be formed. Because of the morphology of these stripped galaxies they have been referred to as Jellyfish galaxies. It has been found that star formation is triggered not only in the disk, but also in the tenta… ▽ More When falling into a galaxy cluster, galaxies experience a loss of gas due to ram pressure stripping. In particular, disk galaxies lose gas from their disks and very large tentacles of gas can be formed. Because of the morphology of these stripped galaxies they have been referred to as Jellyfish galaxies. It has been found that star formation is triggered not only in the disk, but also in the tentacles of such Jellyfish galaxies. The observed star forming regions located in the tentacles of those galaxies have been found to be as massive as $3\times10^7$ M$_{\odot}$ and with sizes $> 100$ pc. Interestingly, these parameters in mass and size agree with those of dwarf galaxies. In this work we make use of the state of the art magneto-hydrodynamical cosmological simulation Illustris TNG-50, to study massive jellyfish galaxies with long tentacles. We find that, in the tentacles of TNG-50 Jellyfish galaxies, the star formation regions (gas+stars) formed could be as massive as $\sim2\times10^8$ M$_{\odot}$. A particular star forming region was analyzed. This region has a star formation rate of $0.04$ M$_{\odot}$/yr, it is metal rich, has an average age of $0.46$ Gyr, and has a half mass radius of $\sim1$ kpc, typical of standard dwarf galaxies. Most importantly, this region is gravitationally self-bound. All and all, we identify a new type of dwarf galaxy being born from the gas tentacles of jellyfish galaxies, that by construction lacks a dark matter (hereafter DM) halo. △ Less

Submitted 8 April, 2024; originally announced April 2024.

Comments: 11 pages, 6 figures, ApJ Accepted

The persistence of high altitude non-equilibrium diffuse ionized gas in simulations of star forming galaxies

Authors: Lewis McCallum, Kenneth Wood, Robert Benjamin, Camilo Peñaloza, Dhanesh Krishnarao, Rowan Smith, Bert Vandenbroucke

Abstract: Widespread, high altitude, diffuse ionized gas with scale heights of around a kiloparsec is observed in the Milky Way and other star forming galaxies. Numerical radiation-magnetohydrodynamic simulations of a supernova-driven turbulent interstellar medium show that gas can be driven to high altitudes above the galactic midplane, but the degree of ionization is often less than inferred from observat… ▽ More Widespread, high altitude, diffuse ionized gas with scale heights of around a kiloparsec is observed in the Milky Way and other star forming galaxies. Numerical radiation-magnetohydrodynamic simulations of a supernova-driven turbulent interstellar medium show that gas can be driven to high altitudes above the galactic midplane, but the degree of ionization is often less than inferred from observations. For computational expediency, ionizing radiation from massive stars is often included as a post-processing step assuming ionization equilibrium. We extend our simulations of a Milky Way-like interstellar medium to include the combined effect of supernovae and photoionization feedback from midplane OB stars and a population of hot evolved low mass stars. The diffuse ionized gas has densities below 0.1 ${\rm cm^{-3}}$, so recombination timescales can exceed millions of years. Our simulations now follow the time-dependent ionization and recombination of low density gas. The long recombination timescales result in diffuse ionized gas that persists at large altitudes long after the deaths of massive stars that produce the vast majority of the ionized gas. The diffuse ionized gas does not exhibit the large variability inherent in simulations that adopt ionization equilibrium. The vertical distribution of neutral and ionized gas is close to what is observed in the Milky Way. The volume filling factor of ionized gas increases with altitude resulting in the scale height of free electrons being larger than that inferred from H$α$ emission, thus reconciling the observations of ionized gas made in H$α$ and from pulsar dispersion measurements. △ Less

Submitted 9 April, 2024; v1 submitted 8 April, 2024; originally announced April 2024.

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

Finite Sample Frequency Domain Identification

Authors: Anastasios Tsiamis, Mohamed Abdalmoaty, Roy S. Smith, John Lygeros

Abstract: We study non-parametric frequency-domain system identification from a finite-sample perspective. We assume an open loop scenario where the excitation input is periodic and consider the Empirical Transfer Function Estimate (ETFE), where the goal is to estimate the frequency response at certain desired (evenly-spaced) frequencies, given input-output samples. We show that under sub-Gaussian colored n… ▽ More We study non-parametric frequency-domain system identification from a finite-sample perspective. We assume an open loop scenario where the excitation input is periodic and consider the Empirical Transfer Function Estimate (ETFE), where the goal is to estimate the frequency response at certain desired (evenly-spaced) frequencies, given input-output samples. We show that under sub-Gaussian colored noise (in time-domain) and stability assumptions, the ETFE estimates are concentrated around the true values. The error rate is of the order of $\mathcal{O}((d_{\mathrm{u}}+\sqrt{d_{\mathrm{u}}d_{\mathrm{y}}})\sqrt{M/N_{\mathrm{tot}}})$, where $N_{\mathrm{tot}}$ is the total number of samples, $M$ is the number of desired frequencies, and $d_{\mathrm{u}},\,d_{\mathrm{y}}$ are the dimensions of the input and output signals respectively. This rate remains valid for general irrational transfer functions and does not require a finite order state-space representation. By tuning $M$, we obtain a $N_{\mathrm{tot}}^{-1/3}$ finite-sample rate for learning the frequency response over all frequencies in the $ \mathcal{H}_{\infty}$ norm. Our result draws upon an extension of the Hanson-Wright inequality to semi-infinite matrices. We study the finite-sample behavior of ETFE in simulations. △ Less

Submitted 1 April, 2024; originally announced April 2024.

Testing kinematic distances under a realistic Galactic potential

Authors: Glen H. Hunter, Mattia C. Sormani, Jan P. Beckmann, Eugene Vasiliev, Simon C. O. Glover, Ralf S. Klessen, Juan D. Soler, Noé Brucy, Philipp Girichidis, Junia Göller, Loke Ohlin, Robin Tress, Sergio Molinari, Ortwin Gerhard, Milena Benedettini, Rowan Smith, Patrick Hennebelle, Leonardo Testi

Abstract: Obtaining reliable distance estimates to gas clouds within the Milky Way is challenging in the absence of certain tracers. The kinematic distance approach has been used as an alternative, derived from the assumption of circular trajectories around the Galactic centre. Consequently, significant errors are expected in regions where gas flow deviates from purely circular motions. We aim to quantify t… ▽ More Obtaining reliable distance estimates to gas clouds within the Milky Way is challenging in the absence of certain tracers. The kinematic distance approach has been used as an alternative, derived from the assumption of circular trajectories around the Galactic centre. Consequently, significant errors are expected in regions where gas flow deviates from purely circular motions. We aim to quantify the systematic errors that arise from the kinematic distance method in the presence of a Galactic potential that is non-axisymmetric. We investigate how these errors differ in certain regions of the Galaxy and how they relate to the underlying dynamics. We perform 2D isothermal hydrodynamical simulation of the gas disk with the moving-mesh code Arepo, adding the capability of using an external potential provided by the Agama library for galactic dynamics. We introduce a new analytic potential of the Milky Way, taking elements from existing models and adjusting parameters to match recent observational constraints. We find significant errors in the kinematic distance estimate for gas close to the Sun, along sight lines towards the Galactic centre and anti-centre, and significant deviations associated with the Galactic bar. Kinematic distance errors are low within the spiral arms as gas resides close to local potential minima and the resulting line-of-sight velocity is close to what is expected for an axisymmetric potential. Interarm regions exhibit large deviations at any given Galactic radius. This is caused by the gas being sped up or slowed down as it travels into or out of the spiral arm. We are able to define 'zones of avoidance' in the lv-diagram, where the kinematic distance method is particularly unreliable and should only be used with caution. We report a power law relation between the kinematic distance error and the deviation of the project line-of-sight velocity from circular motion. △ Less

Submitted 26 March, 2024; originally announced March 2024.

Comments: 23 pages, 22 figures, 6 tables, submitted to A&A

Small Noise Analysis of Non-Parametric Closed-Loop Identification

Authors: Mohamed Abdalmoaty, Roy S. Smith

Abstract: We revisit the problem of non-parametric closed-loop identification in frequency domain; we give a brief survey of the literature and provide a small noise analysis of the direct, indirect, and joint input-output methods when two independent experiments with identical excitation are used. The analysis is asymptotic in the noise variance (i.e., as the standard deviation of the innovations $σ\to 0$)… ▽ More We revisit the problem of non-parametric closed-loop identification in frequency domain; we give a brief survey of the literature and provide a small noise analysis of the direct, indirect, and joint input-output methods when two independent experiments with identical excitation are used. The analysis is asymptotic in the noise variance (i.e., as the standard deviation of the innovations $σ\to 0$), for a finite data record of length $N$. We highlight the relationship between the estimators accuracy and the loop shape via asymptotic variance expressions given in terms of the sensitivity function. The results are illustrated using a numerical simulation example. △ Less

Submitted 23 March, 2024; originally announced March 2024.

Optimal Data-Driven Prediction and Predictive Control using Signal Matrix Models

Authors: Roy S. Smith, Mohamed Abdalmoaty, Mingzhou Yin

Abstract: Data-driven control uses a past signal trajectory to characterise the input-output behaviour of a system. Willems' lemma provides a data-based prediction model allowing a control designer to bypass the step of identifying a state-space or transfer function model. This paper provides a more parsimonious formulation of Willems' lemma that separates the model into initial condition matching and predi… ▽ More Data-driven control uses a past signal trajectory to characterise the input-output behaviour of a system. Willems' lemma provides a data-based prediction model allowing a control designer to bypass the step of identifying a state-space or transfer function model. This paper provides a more parsimonious formulation of Willems' lemma that separates the model into initial condition matching and predictive control design parts. This avoids the need for regularisers in the predictive control problem that are found in other data-driven predictive control methods. It also gives a closed form expression for the optimal (minimum variance) unbiased predictor of the future output trajectory and applies it for predictive control. Simulation comparisons illustrate very good control performance. △ Less

Submitted 22 March, 2024; originally announced March 2024.

Comments: 7 pages, 3 figures. Submitted to IEEE Control Systems Society Letters and the 2024 Conference on Decision and Control

Magnetic field morphology and evolution in the Central Molecular Zone and its effect on gas dynamics

Authors: R. G. Tress, M. C. Sormani, P. Girichidis, S. C. O. Glover, R. S. Klessen, R. J. Smith, E. Sobacchi, L. Armillotta, A. T. Barnes, C. Battersby, K. R. J. Bogue, N. Brucy, L. Colzi, C. Federrath, P. García, A. Ginsburg, J. Göller, H P. Hatchfield, C. Henkel, P. Hennebelle, J. D. Henshaw, M. Hirschmann, Y. Hu, J. Kauffmann, J. M. D. Kruijssen , et al. (12 additional authors not shown)

Abstract: The interstellar medium in the Milky Way's Central Molecular Zone (CMZ) is known to be strongly magnetised, but its large-scale morphology and impact on the gas dynamics are not well understood. We explore the impact and properties of magnetic fields in the CMZ using three-dimensional non-self gravitating magnetohydrodynamical simulations of gas flow in an external Milky Way barred potential. We f… ▽ More The interstellar medium in the Milky Way's Central Molecular Zone (CMZ) is known to be strongly magnetised, but its large-scale morphology and impact on the gas dynamics are not well understood. We explore the impact and properties of magnetic fields in the CMZ using three-dimensional non-self gravitating magnetohydrodynamical simulations of gas flow in an external Milky Way barred potential. We find that: (1) The magnetic field is conveniently decomposed into a regular time-averaged component and an irregular turbulent component. The regular component aligns well with the velocity vectors of the gas everywhere, including within the bar lanes. (2) The field geometry transitions from parallel to the Galactic plane near $z=0$ to poloidal away from the plane. (3) The magneto-rotational instability (MRI) causes an in-plane inflow of matter from the CMZ gas ring towards the central few parsecs of $0.01-0.1$ M$_\odot$ yr$^{-1}$ that is absent in the unmagnetised simulations. However, the magnetic fields have no significant effect on the larger-scale bar-driven inflow that brings the gas from the Galactic disc into the CMZ. (4) A combination of bar inflow and MRI-driven turbulence can sustain a turbulent vertical velocity dispersion of $σ_z \simeq 5$ km s$^{-1}$ on scales of $20$ pc in the CMZ ring. The MRI alone sustains a velocity dispersion of $σ_z \simeq 3$ km s$^{-1}$. Both these numbers are lower than the observed velocity dispersion of gas in the CMZ, suggesting that other processes such as stellar feedback are necessary to explain the observations. (5) Dynamo action driven by differential rotation and the MRI amplifies the magnetic fields in the CMZ ring until they saturate at a value that scales with the average local density as $B \simeq 102 (n/10^3 cm^{-3})^{0.33}$ $μ$G. Finally, we discuss the implications of our results within the observational context in the CMZ. △ Less

Submitted 19 March, 2024; originally announced March 2024.

Comments: submitted to A&A, comments welcome

de Leeuw representations of functionals on Lipschitz spaces

Authors: Ramón J. Aliaga, E. Pernecká, Richard J. Smith

Abstract: Let $\mathrm{Lip}_0(M)$ be the space of Lipschitz functions on a complete metric space $(M,d)$ that vanish at a point $0\in M$. We investigate its dual $\mathrm{Lip}_0(M)^*$ using the de Leeuw transform, which allows representing each functional on $\mathrm{Lip}_0(M)$ as a (non-unique) measure on $β\widetilde{M}$, where $\widetilde{M}$ is the space of pairs $(x,y)\in M\times M$, $x\neq y$. We dist… ▽ More Let $\mathrm{Lip}_0(M)$ be the space of Lipschitz functions on a complete metric space $(M,d)$ that vanish at a point $0\in M$. We investigate its dual $\mathrm{Lip}_0(M)^*$ using the de Leeuw transform, which allows representing each functional on $\mathrm{Lip}_0(M)$ as a (non-unique) measure on $β\widetilde{M}$, where $\widetilde{M}$ is the space of pairs $(x,y)\in M\times M$, $x\neq y$. We distinguish a set of points of $β\widetilde{M}$ that are "away from infinity", which can be assigned coordinates belonging to the Lipschitz realcompactification $M^{\mathcal{R}}$ of $M$. We define a natural metric $\bar{d}$ on $M^{\mathcal{R}}$ extending $d$ and we show that optimal (i.e. positive and norm-minimal) de Leeuw representations of well-behaved functionals are characterised by $\bar{d}$-cyclical monotonicity of their support, extending known results for functionals in $\mathcal{F}(M)$, the predual of $\mathrm{Lip}_0(M)$. We also extend the Kantorovich-Rubinstein theorem to normal Hausdorff spaces, in particular to $M^{\mathcal{R}}$, and use this to characterise measure-induced and majorisable functionals in $\mathrm{Lip}_0(M)^*$ as those admitting optimal representations with additional finiteness properties. Finally, we use de Leeuw representations to define a natural L-projection of $\mathrm{Lip}_0(M)^*$ onto $\mathcal{F}(M)$ under some conditions on $M$. △ Less

Submitted 9 May, 2024; v1 submitted 14 March, 2024; originally announced March 2024.

MSC Class: Primary 46B20; 49Q22; Secondary 46E15; 54D35; 54E50

Online Identification of Stochastic Continuous-Time Wiener Models Using Sampled Data

Authors: Mohamed Abdalmoaty, Efe C. Balta, John Lygeros, Roy S. Smith

Abstract: It is well known that ignoring the presence of stochastic disturbances in the identification of stochastic Wiener models leads to asymptotically biased estimators. On the other hand, optimal statistical identification, via likelihood-based methods, is sensitive to the assumptions on the data distribution and is usually based on relatively complex sequential Monte Carlo algorithms. We develop a sim… ▽ More It is well known that ignoring the presence of stochastic disturbances in the identification of stochastic Wiener models leads to asymptotically biased estimators. On the other hand, optimal statistical identification, via likelihood-based methods, is sensitive to the assumptions on the data distribution and is usually based on relatively complex sequential Monte Carlo algorithms. We develop a simple recursive online estimation algorithm based on an output-error predictor, for the identification of continuous-time stochastic parametric Wiener models through stochastic approximation. The method is applicable to generic model parameterizations and, as demonstrated in the numerical simulation examples, it is robust with respect to the assumptions on the spectrum of the disturbance process. △ Less

Submitted 9 March, 2024; originally announced March 2024.

Characterizing the 3D Structure of Molecular Cloud Envelopes in the "Cloud Factory" Simulations

Authors: Elijah Mullens, Catherine Zucker, Claire E. Murray, Rowan Smith

Abstract: We leverage recent numerical simulations of highly resolved star-forming regions in a Milky Way-like Galaxy to explore the nature of extended gaseous envelopes around molecular clouds. We extract a sample of two dozen star-forming clouds from the feedback-dominated suite of the "Cloud Factory'' simulations. With the goal of exploring the 3D thermal and chemical structure of the gas, we measure and… ▽ More We leverage recent numerical simulations of highly resolved star-forming regions in a Milky Way-like Galaxy to explore the nature of extended gaseous envelopes around molecular clouds. We extract a sample of two dozen star-forming clouds from the feedback-dominated suite of the "Cloud Factory'' simulations. With the goal of exploring the 3D thermal and chemical structure of the gas, we measure and fit the clouds' radial profiles in multiple tracers, including $n_{H_I}$, $n_{H_2}$, $n_{H_{tot}}$, $n_{CO}$, and gas temperature. We find that while solar neighborhood clouds recently detected via 3D dust mapping have radially symmetric, low-density envelopes that extend $\sim$ 10-15 pc, the simulated cloud envelopes are primarily radially asymmetric with low-density envelopes that extend only $\sim$ 2-3 pc. One potential explanation for the absence of extended envelopes in the simulated clouds may be the lack of magnetic fields, while a stronger local feedback prescription compared to solar neighborhood conditions may drive the radially asymmetric cloud morphologies. We make the pipeline to extract and characterize the radial profile of clouds publicly available, which can be used in complementary and future simulations to shed additional light on the key physics shaping the formation and evolution of star-forming structures in the Milky Way. △ Less

Submitted 5 March, 2024; originally announced March 2024.

Comments: 22 pages, 13 figures, 2 tables. Accepted for publication in ApJ (March 2024)

Ultralight vector dark matter search using data from the KAGRA O3GK run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, R. Abbott, H. Abe, I. Abouelfettouh, F. Acernese, K. Ackley, C. Adamcewicz, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, I. Aguilar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi , et al. (1778 additional authors not shown)

Abstract: Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we prese… ▽ More Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we present the result of a search for $U(1)_{B-L}$ gauge boson DM using the KAGRA data from auxiliary length channels during the first joint observation run together with GEO600. By applying our search pipeline, which takes into account the stochastic nature of ultralight DM, upper bounds on the coupling strength between the $U(1)_{B-L}$ gauge boson and ordinary matter are obtained for a range of DM masses. While our constraints are less stringent than those derived from previous experiments, this study demonstrates the applicability of our method to the lower-mass vector DM search, which is made difficult in this measurement by the short observation time compared to the auto-correlation time scale of DM. △ Less

Submitted 5 March, 2024; originally announced March 2024.

Comments: 20 pages, 5 figures

Report number: LIGO-P2300250

Fires in the deep: The luminosity distribution of early-time gamma-ray-burst afterglows in light of the Gamow Explorer sensitivity requirements

Authors: D. A. Kann, N. E. White, G. Ghirlanda, S. R. Oates, A. Melandri, M. Jelinek, A. de Ugarte Postigo, A. J. Levan, A. Martin-Carrillo, G. S. -H. Paek, L. Izzo, M. Blazek, C. Thone, J. F. Agui Fernandez, R. Salvaterra, N. R. Tanvir, T. -C. Chang, P. O'Brien, A. Rossi, D. A. Perley, M. Im, D. B. Malesani, A. Antonelli, S. Covino, C. Choi , et al. (36 additional authors not shown)

Abstract: Gamma-ray bursts (GRBs) are ideal probes of the Universe at high redshift (z > 5), pinpointing the locations of the earliest star-forming galaxies and providing bright backlights that can be used to spectrally fingerprint the intergalactic medium and host galaxy during the period of reionization. Future missions such as Gamow Explorer are being proposed to unlock this potential by increasing the r… ▽ More Gamma-ray bursts (GRBs) are ideal probes of the Universe at high redshift (z > 5), pinpointing the locations of the earliest star-forming galaxies and providing bright backlights that can be used to spectrally fingerprint the intergalactic medium and host galaxy during the period of reionization. Future missions such as Gamow Explorer are being proposed to unlock this potential by increasing the rate of identification of high-z GRBs to rapidly trigger observations from 6-10 m ground telescopes, JWST, and the Extremely Large Telescopes. Gamow was proposed to the NASA 2021 Medium-Class Explorer (MIDEX) program as a fast-slewing satellite featuring a wide-field lobster-eye X-ray telescope (LEXT) to detect and localize GRBs, and a 30 cm narrow-field multi-channel photo-z infrared telescope (PIRT) to measure their photometric redshifts using the Lyman-alpha dropout technique. To derive the PIRT sensitivity requirement we compiled a complete sample of GRB optical-near-infrared afterglows from 2008 to 2021, adding a total of 66 new afterglows to our earlier sample, including all known high-z GRB afterglows. We performed full light-curve and spectral-energy-distribution analyses of these afterglows to derive their true luminosity at very early times. For all the light curves, where possible, we determined the brightness at the time of the initial finding chart of Gamow, at different high redshifts and in different NIR bands. We then followed the evolution of the luminosity to predict requirements for ground and space-based follow-up. We find that a PIRT sensitivity of 15 micro-Jy (21 mag AB) in a 500 s exposure simultaneously in five NIR bands within 1000s of the GRB trigger will meet the Gamow mission requirement to recover > 80% of all redshifts at z > 5. △ Less

Submitted 29 February, 2024; originally announced March 2024.

Comments: 44 pages, 10 figures, 1 table. Accepted for publication in Astronomy and Astrophysics 15 Feb 2024. Abstract abridged for arXiv

Journal ref: A&A 686, A56 (2024)

Multi-messenger Astrophysics of Black Holes and Neutron Stars as Probed by Ground-based Gravitational Wave Detectors: From Present to Future

Authors: Alessandra Corsi, Lisa Barsotti, Emanuele Berti, Matthew Evans, Ish Gupta, Konstantinos Kritos, Kevin Kuns, Alexander H. Nitz, Benjamin J. Owen, Binod Rajbhandari, Jocelyn Read, Bangalore S. Sathyaprakash, David H. Shoemaker, Joshua R. Smith, Salvatore Vitale

Abstract: The ground-based gravitational wave (GW) detectors LIGO and Virgo have enabled the birth of multi-messenger GW astronomy via the detection of GWs from merging stellar-mass black holes (BHs) and neutron stars (NSs). GW170817, the first binary NS merger detected in GWs and all bands of the electromagnetic spectrum, is an outstanding example of the impact that GW discoveries can have on multi-messeng… ▽ More The ground-based gravitational wave (GW) detectors LIGO and Virgo have enabled the birth of multi-messenger GW astronomy via the detection of GWs from merging stellar-mass black holes (BHs) and neutron stars (NSs). GW170817, the first binary NS merger detected in GWs and all bands of the electromagnetic spectrum, is an outstanding example of the impact that GW discoveries can have on multi-messenger astronomy. Yet, GW170817 is only one of the many and varied multi-messenger sources that can be unveiled using ground-based GW detectors. In this contribution, we summarize key open questions in the astrophysics of stellar-mass BHs and NSs that can be answered using current and future-generation ground-based GW detectors, and highlight the potential for new multi-messenger discoveries ahead. △ Less

Submitted 20 February, 2024; originally announced February 2024.

Comments: Review submitted to the Frontiers Research Topic "The Dynamic Universe: Realizing the Potential of Time Domain and Multimessenger Astrophysics"

New constraints on ultraheavy dark matter from the LZ experiment

Authors: J. Aalbers, D. S. Akerib, A. K. Al Musalhi, C. S. Amarasinghe, A. Ames, T. J. Anderson, N. Angelides, H. M. Araújo, J. E. Armstrong, M. Arthurs, A. Baker, S. Balashov, J. Bang, J. W. Bargemann, A. Baxter, K. Beattie, T. Benson, A. Bhatti, A. Biekert, T. P. Biesiadzinski, H. J. Birch, E. Bishop, G. M. Blockinger, B. Boxer, C. A. J. Brew , et al. (174 additional authors not shown)

Abstract: Searches for dark matter with liquid xenon time projection chamber experiments have traditionally focused on the region of the parameter space that is characteristic of weakly interacting massive particles, ranging from a few GeV/$c^2$ to a few TeV/$c^2$. Models of dark matter with a mass much heavier than this are well motivated by early production mechanisms different from the standard thermal f… ▽ More Searches for dark matter with liquid xenon time projection chamber experiments have traditionally focused on the region of the parameter space that is characteristic of weakly interacting massive particles, ranging from a few GeV/$c^2$ to a few TeV/$c^2$. Models of dark matter with a mass much heavier than this are well motivated by early production mechanisms different from the standard thermal freeze-out, but they have generally been less explored experimentally. In this work, we present a re-analysis of the first science run (SR1) of the LZ experiment, with an exposure of $0.9$ tonne$\times$year, to search for ultraheavy particle dark matter. The signal topology consists of multiple energy deposits in the active region of the detector forming a straight line, from which the velocity of the incoming particle can be reconstructed on an event-by-event basis. Zero events with this topology were observed after applying the data selection calibrated on a simulated sample of signal-like events. New experimental constraints are derived, which rule out previously unexplored regions of the dark matter parameter space of spin-independent interactions beyond a mass of 10$^{17}$ GeV/$c^2$. △ Less

Submitted 13 February, 2024; originally announced February 2024.

Comments: 9 pages, 7 figures

Journal ref: Phys. Rev. D 109, 112010 (2024)

Low frequency coherent Raman imaging robust to optical scattering

Authors: David R. Smith, Jesse W. Wilson, Siddarth Shivkumar, Herve Rigneault, Randy A. Bartels

Abstract: We demonstrate low-frequency interferometric impulsive stimulated Raman scattering (ISRS) imaging with high robustness to distortions by optical scattering. ISRS is a pump-probe coherent Raman spectroscopy that can capture Raman vibrational spectra. Recording of ISRS spectra requires isolation of a probe pulse from the pump pulse. While this separation is simple in non-scattering specimens, such a… ▽ More We demonstrate low-frequency interferometric impulsive stimulated Raman scattering (ISRS) imaging with high robustness to distortions by optical scattering. ISRS is a pump-probe coherent Raman spectroscopy that can capture Raman vibrational spectra. Recording of ISRS spectra requires isolation of a probe pulse from the pump pulse. While this separation is simple in non-scattering specimens, such as liquids, scattering leads to significant pump pulse contamination and prevent the extraction of a Raman spectrum. We introduce a robust method for ISRS microscopy that works in complex scattering samples. High signal-to-noise ISRS spectra are obtained even when the pump and probe pulses pass through many scattering layers. △ Less

Submitted 10 February, 2024; originally announced February 2024.

Comments: 13 pages; 5 figures

Measuring line tension: thermodynamic integration during detachment of a molecular dynamics droplet

Authors: Minori Shintaku, Haruki Oga, Hiroki Kusudo, Edward R. Smith, Takeshi Omori, Yasutaka Yamaguchi

Abstract: The contact line (CL) is where solid, liquid and vapor phases meet, and Young's equation describes the macroscopic force balance of the interfacial tensions between these three phases. These interfacial tensions are related to the nanoscale stress inhomogeneity appearing around the interface, and for curved CLs, eg a three-dimensional droplet, another force known as the line tension must be includ… ▽ More The contact line (CL) is where solid, liquid and vapor phases meet, and Young's equation describes the macroscopic force balance of the interfacial tensions between these three phases. These interfacial tensions are related to the nanoscale stress inhomogeneity appearing around the interface, and for curved CLs, eg a three-dimensional droplet, another force known as the line tension must be included in Young's equation. The line tension has units of force, acting parallel to the CL, and is required to incorporate the extra stress inhomogeneity around the CL into the force balance. Considering this feature, Bey et al. [J. Chem. Phys. \textbf{152}, 094707 (2020)] reported a mechanical approach to extract the value of line tension $\taul$ from molecular dynamics (MD) simulations. In this study, we show a novel thermodynamics interpretation of the line tension as the free energy per CL length, and based on this interpretation, through MD simulations of a quasi-static detachment process of a quasi-two-dimensional droplet from a solid surface, we obtained the value $\taul$ as a function of the contact angle. The simulation scheme is considered to be an extension of a thermodynamic integration method, previously used to calculate the solid-liquid and solid-vapor interfacial tensions through a detachment process, extended here to the three phase system. The obtained value agreed well with the result by Bey et al. and show the validity of thermodynamic integration at the three-phase interface. △ Less

Submitted 9 February, 2024; originally announced February 2024.

Comments: 9 figures

On the evolution of the observed Mass-to-Length relationship for star-forming filaments

Authors: Jiancheng Feng, Rowan J. Smith, Alvaro Hacar, Susan E. Clark, Daniel Seifried

Abstract: The interstellar medium is threaded by a hierarchy of filaments from large scales (~ 100 pc) to small scales (~ 0.1pc). The masses and lengths of these nested structures may reveal important constraints for cloud formation and evolution, but it is difficult to investigate from an evolutionary perspective using single observations. In this work, we extract simulated molecular clouds from the Cloud… ▽ More The interstellar medium is threaded by a hierarchy of filaments from large scales (~ 100 pc) to small scales (~ 0.1pc). The masses and lengths of these nested structures may reveal important constraints for cloud formation and evolution, but it is difficult to investigate from an evolutionary perspective using single observations. In this work, we extract simulated molecular clouds from the Cloud Factory galactic-scale ISM suite in combination with 3D Monte Carlo radiative transfer code POLARIS to investigate how filamentary structure evolves over time. We produce synthetic dust continuum observations in three regions with a series of snapshots and use the Filfinder algorithm to identify filaments in the dust derived column density maps. When the synthetic filaments mass and length are plotted on an M-L plot, we see a scaling relation of $L\propto M^{0.45}$ similar to that seen in observations, and find that the filaments are thermally supercritical. Projection effects systematically affect the masses and lengths measured for the filaments, and are particularly severe in crowded regions. In the filament Mass-Length (M-L) diagram we identify three main evolutionary mechanisms: accretion, segmentation, and dispersal. In particular we find that the filaments typically evolve from smaller to larger masses in the observational M-L plane, indicating the dominant role of accretion in filament evolution. Moreover, we find a potential correlation between line mass and filament growth rate. Once filaments are actively star forming they then segment into smaller sections, or are dispersed by internal or external forces. △ Less

Submitted 7 February, 2024; originally announced February 2024.

Comments: 19 pages, 18 figures

Context-Aware Spectrum Coexistence of Terrestrial Beyond 5G Networks in Satellite Bands

Authors: Ta Seen Reaz Niloy, Zoheb Hasan, Rob Smith, Vikram R. Anapana, Vijay K. Shah

Abstract: Spectrum sharing between terrestrial 5G and incumbent networks in the satellite bands presents a promising avenue to satisfy the ever-increasing bandwidth demand of the next-generation wireless networks. However, protecting incumbent operations from harmful interference poses a fundamental challenge in accommodating terrestrial broadband cellular networks in the satellite bands. State-of-the-art s… ▽ More Spectrum sharing between terrestrial 5G and incumbent networks in the satellite bands presents a promising avenue to satisfy the ever-increasing bandwidth demand of the next-generation wireless networks. However, protecting incumbent operations from harmful interference poses a fundamental challenge in accommodating terrestrial broadband cellular networks in the satellite bands. State-of-the-art spectrum-sharing policies usually consider several worst-case assumptions and ignore site-specific contextual factors in making spectrum-sharing decisions, and thus, often results in under-utilization of the shared band for the secondary licensees. To address such limitations, this paper introduces CAT3S (Context-Aware Terrestrial-Satellite Spectrum Sharing) framework that empowers the coexisting terrestrial 5G network to maximize utilization of the shared satellite band without creating harmful interference to the incumbent links by exploiting the contextual factors. CAT3S consists of the following two components: (i) context-acquisition unit to collect and process essential contextual information for spectrum sharing and (ii) context-aware base station (BS) control unit to optimize the set of operational BSs and their operation parameters (i.e., transmit power and active beams per sector). To evaluate the performance of the CAT3S, a realistic spectrum coexistence case study over the 12 GHz band is considered. Experiment results demonstrate that the proposed CAT3S achieves notably higher spectrum utilization than state-of-the-art spectrum-sharing policies in different weather contexts. △ Less

Submitted 14 February, 2024; v1 submitted 6 February, 2024; originally announced February 2024.