Search

Bayesian phylogenetics typically estimates a posterior distribution, or aspects thereof, using Markov chain Monte Carlo methods. These methods integrate over tree space by applying local rearrangements to move a tree through its space as a random walk. Previous work explored the possibility of replacing this random walk with a systematic search, but was quickly overwhelmed by the large number of probable trees in the posterior distribution. In this paper we develop methods to sidestep this problem using a recently introduced structure called the subsplit directed acyclic graph (sDAG). This structure can represent many trees at once, and local rearrangements of trees translate to methods of enlarging the sDAG. Here we propose two methods of introducing, ranking, and selecting local rearrangements on sDAGs to produce a collection of trees with high posterior density. One of these methods successfully recovers the set of high posterior density trees across a range of data sets. However, we find that a simpler strategy of aggregating trees into an sDAG in fact is computationally faster and returns a higher fraction of probable trees.

State-of-the-art noisy-intermediate-scale quantum (NISQ) processors are currently implemented across a variety of hardware platforms, each with their own distinct gatesets. As such, circuit compilation should not only be aware of, but also deeply connect to, the native gateset and noise properties of each. Trapped-ion processors are one such platform that provides a gateset that can be continuously parameterized across both one- and two-qubit gates. Here we use the Quantum Scientific Computing Open User Testbed (QSCOUT) to study noise-aware compilations focused on continuously parameterized two-qubit $\mathcal{ZZ}$ gates (based on the M{\o}lmer-S{\o}rensen interaction) using $\textbf{Superstaq}$, a quantum software platform for hardware-aware circuit compiler optimizations. We discuss the realization of $\mathcal{ZZ}$ gates with arbitrary angle on the all-to-all connected trapped-ion system. Then we discuss a variety of different compiler optimizations that innately target these $\mathcal{ZZ}$ gates and their noise properties. These optimizations include moving from a restricted maximally entangling gateset to a continuously parameterized one, swap mirroring to further reduce total entangling angle of the operations, focusing the heaviest $\mathcal{ZZ}$ angle participation on the best performing gate pairs, and circuit approximation to remove the least impactful $\mathcal{ZZ}$ gates. We demonstrate these compilation approaches on the hardware with randomized quantum volume circuits, observing the potential to realize a larger quantum volume as a result of these optimizations. Using differing yet complementary analysis techniques, we observe the distinct improvements in system performance provided by these noise-aware compilations and study the role of stochastic and coherent error channels for each compilation choice., Comment: 17 pages, 13 figures

Selecting materials for hybrid cathodes for batteries, which combine intercalation and conversion materials, has gained interest due to their unique synergistic properties, which are not achievable by homogeneous materials. Here, we present a data-driven, chemistry-agnostic, inverse material design framework to discover hybrid cathode materials (HCMs) for metal-based batteries. This framework systematically explores the material space for any working ion, evaluates the stability of candidates, and identifies the growth modes and adsorption of components for a stable hybrid cathode. To demonstrate the framework's application, we conducted a case study aimed at discovering HCMs with an average gravimetric energy density that exceeds the widely used high-energy NMC333 cathode material. The framework identified LiCr$_4$GaS$_8$-Li$_2$S as a promising HCM, achieving an average energy density of 1,424 Wh/kg (on a lithiated cathode basis), surpassing NMC333's theoretical maximum of 1,028 Wh/kg. Additional features of the identified material include: (1) thermodynamically stable lithiated and delithiated intercalation and conversion phases; (2) minimal volume change upon (de)lithiation, which mitigates the conversion material's high volume change; (3) high energy density, compensating for the lower energy density of the intercalation material; (4) the intercalation component acting as both a conductive additive and a sulfur immobilizer while contributing to the total energy density; (5) serving as an ideal support for soft sulfur species; and (6) potential improvements in lifespan, self-discharge, mechanical integrity, and capacity fading compared to conventional Li-S batteries. The developed framework was crucial in exploring materials within the vast potential HCM space with predefined battery material design objectives.

Early turbulence in periodic cylinder arrays is of particular interest in many practical applications to enhance mixing and material/heat exchange. In this study, we reveal a new early transition pathway to a chaotic wavy state and drag enhancement with the addition of polymers. Using 2D direct numerical simulations with sufficiently small polymer diffusion ($\epsilon=10^{-5}$), we show that viscoelastic flow past periodic cylinder arrays become unstable at a Reynolds number $Re\approx 10$, significantly lower than the Newtonian counterpart of $Re\approx 150-200$. The chaotic wavy state which ensues exhibits sheets and `arrowhead' polymer conformation structures, consistent with the saturated centre-mode instability observed in wall-bounded parallel flows (Page et al., Phys. Rev. Lett., 125, 154501, 2020). Analysis of the kinetic energy budget reveals the purely elastic origin of the chaos. However, inertial forces, in conjunction with elastic forces, can reshape the base state, affecting the formation of an invariant polymer sheet. This sheet facilitates the stretching and recoiling of polymers, which in turn induces flow fluctuations and maintains the chaos. Exploring various maximum polymer extensions $b$, and polymer concentrations $\beta$ highlights the role of elastic forces in stretching the upstream separation zone while suppressing the downstream separation zone, resulting in drag enhancement at finite $Re$. Surprisingly, these modifications to the base state can suppress the invariant polymer sheet under large elastic forces (large $b$ or small $\beta$), thereby achieving a stable polymer-modified laminar state.

We study the stellar mass function (SMF) and the co-evolution with dark matter halos via abundance matching in the largest redshift range to date $0.25$, we find increased abundances of massive (log$\, M_{\star}/M_{\odot}>10.5$) implying integrated star formation efficiencies (SFE) $\epsilon_{\star}\equiv M_{\star}\, f_{\rm b}^{-1} M_{\rm halo}^{-1} \gtrsim 0.5$. We find a flattening of the SMF at the high-mass end that is better described by a double power law at $z>5.5$. At $z \lesssim 5.5$ it transitions to a Schechter law which coincides with the emergence of the first massive quiescent galaxies in the Universe. We trace the cosmic stellar mass density (SMD) and infer the star formation rate density (SFRD), which at $z>7.5$ agrees remarkably with recent \JWST{} UV luminosity function-derived estimates. However, at $z \lesssim 3.5$, we find significant tension ($\sim 0.3$ dex) with the cosmic star formation (SF) history from instantaneous SF measures, the causes of which remain poorly understood. We infer the stellar-to-halo mass relation (SHMR) and the SFE from abundance matching out to $z=12$, finding a non-monotonic evolution. The SFE has the characteristic strong dependence with mass in the range of $0.02 - 0.2$, and mildly decreases at the low mass end out to $z\sim3.5$. At $z\sim3.5$ the SFE increases sharply from $\sim 0.1$ to approach high SFE of $0.8-1$ by $z\sim 10$ for log$(M_{\rm h}/M_{\odot})\approx11.5$, albeit with large uncertainties. Finally, we use the SHMR to track the SFE and stellar mass growth throughout the halo history and find that they do not grow at the same rate -- from the earliest times up until $z\sim3.5$ the halo growth rate outpaces galaxy assembly, but at $z>3.5$ halo growth stagnates and accumulated gas reservoirs keep the SF going and galaxies outpace halos.

We revisit viscoelastic Kolmogorov flow to show that the elastic linear instability of an Oldroyd-B fluid at vanishing Reynolds numbers ($Re$) found by Boffetta et al. (J. Fluid Mech. 523, 161-170, 2005) is the same `centre-mode' instability found at much higher $Re$ by Garg et al. (Phys. Rev. Lett. 121, 024502, 2018) in a pipe and Khalid et al. (J. Fluid Mech. 915, A43, 2021) in a channel. In contrast to these wall-bounded flows, the centre-mode instability exists even when the solvent viscosity vanishes (e.g. it exists in the upper-convective Maxwell limit with $Re=0$). Floquet analysis reveals that the preferred centre-mode instability almost always has a wavelength twice that of the forcing. All elastic instabilities give rise to familiar `arrowheads' (Page et al. Phys. Rev. Lett. 125, 154501, 2020) which in sufficiently large domains and at sufficient Weissenberg number ($W$) interact chaotically in 2D to give elastic turbulence via a bursting scenario. Finally, it is found that the $k^{-4}$ scaling of the kinetic energy spectrum seen in this 2D elastic turbulence is already contained within the component arrowhead structures., Comment: 24 pages, 17 figures

The stellar mass-star formation rate (M$_\star$ - SFR) plane is essential for distinguishing galaxy populations, but how galaxies move within this plane over cosmic time remains unclear. This study aims to describe galaxy migrations in the M$_\star$ - SFR plane by reconstructing star formation histories (SFHs) for a sample of galaxies out to redshift $z < 4$. This provides insights into the physical processes driving star formation. We use data from the COSMOS field, selecting 299131 galaxies at $z < 4$ with COSMOS-Web NIRCam data (m$_\mathrm{F444W} < 27$) over 0.54 deg$^2$. Using the SED modeling code CIGALE with non-parametric SFHs, we derive physical properties and migration vectors for these galaxies. These vectors describe the direction and velocity of evolutionary paths across the M$_\star$ - SFR plane. To assess the accuracy of these vectors, we compare them to results from the Horizon-AGN simulation. Galaxies within the main sequence show low migration amplitudes and dispersed movement directions, indicating oscillation within the main sequence. Most progenitors were already on the main sequence a billion years earlier. Starburst galaxies assembled half their mass in the last 350 Myr and originated from the main sequence. Passive galaxies show uniformly declining SFHs and include massive galaxies already in the passive region at $z = 3.5-4$, with increasing density over time. Using reconstructed SFHs up to $z < 4$, we propose a coherent picture of how galaxies migrate over cosmic time in the M$_\star$ - SFR plane, highlighting the connection between major phases in the star-formation history.

Generalized Additive Models (GAMs) are widely recognized for their ability to create fully interpretable machine learning models for tabular data. Traditionally, training GAMs involves iterative learning algorithms, such as splines, boosted trees, or neural networks, which refine the additive components through repeated error reduction. In this paper, we introduce GAMformer, the first method to leverage in-context learning to estimate shape functions of a GAM in a single forward pass, representing a significant departure from the conventional iterative approaches to GAM fitting. Building on previous research applying in-context learning to tabular data, we exclusively use complex, synthetic data to train GAMformer, yet find it extrapolates well to real-world data. Our experiments show that GAMformer performs on par with other leading GAMs across various classification benchmarks while generating highly interpretable shape functions., Comment: 20 pages, 12 figures

We present the spatially resolved relationship between the dust-to-gas mass ratio (DGR) and gas-phase metallicity (Zgas or 12+log(O/H)) (i.e., DGR-Zgas relation) of 11 nearby galaxies with a large metallicity range (1.5 dex of 12+log(O/H)) at (sub-)kpc scales. We used the large field-of-view (> 3') optical pseudo-Integral Field Spectroscopy data taken by the TYPHOON/PrISM survey, covering the optical size of galaxies, combining them with multi-wavelength data (far-UV to far-IR, CO, and HI 21 cm radio). A large scatter of DGR in the intermediate metallicity galaxies (8.0 < 12+log(O/H) < 8.3) is found, which is in line with dust evolution models, where grain growth begins to dominate the mechanism of dust mass accumulation. In the lowest metallicity galaxy of our sample, Sextans A (12+log(O/H) < 7.6), the star-forming regions have significantly higher DGR values (by 0.5-2 dex) than the global estimates from literature at the same metallicity but aligns with the DGR values from metal depletion method from Damped Lyman Alpha systems and high hydrogen gas density regions of Sextans A. Using dust evolution models with a Bayesian MCMC approach suggests: 1) a high SN dust yield and 2) a negligible amount of photofragmentation by UV radiation, although we note that our sample in the low-metallicity regime is limited to Sextans A. On the other hand, it is also possible that while metallicity influences DGR, gas density also plays a role, indicating an early onset of dust grain growth in the dust mass build-up process despite its low metallicity., Comment: 20 pages, 10 figures, 7 tables, accepted for publication in MNRAS

We present the results of a James Webb Space Telescope (JWST) NIRCam and NIRSpec investigation into the young massive star cluster (YMC) population of NGC 3256, the most cluster-rich luminous infrared galaxy (LIRG) in the Great Observatories All Sky LIRG Survey. We detect 3061 compact YMC candidates with a $S/N \geq 3$ at F150W, F200W, and F335M. Based on yggdrasil stellar population models, we identify 116/3061 sources with F150W - F200W $> 0.47$ and F200W - F355M $> -1.37$ colors suggesting they are young (t $\leq 5$ Myr), dusty ($A_{V} = 5 - 15$), and massive ($M_{\odot} > 10^{5}$). This increases the sample of dust-enshrouded YMCs detected in this system by an order of magnitude relative to previous HST studies. With NIRSpec IFU pointings centered on the northern and southern nucleus, we extract the Pa$\alpha$ and 3.3$\mu$m PAH equivalent widths for 8 bright and isolated YMCs. Variations in both the F200W - F335M color and 3.3$\mu$m PAH emission with the Pa$\alpha$ line strength suggest a rapid dust clearing ($< 3 - 4$ Myr) for the emerging YMCs in the nuclei of NGC 3256. Finally, with both the age and dust emission accurately measured we use yggdrasil to derive the color excess (E(B - V)) for all 8 YMCs. We demonstrate that YMCs with strong 3.3$\mu$m PAH emission (F200W - F335M $> 0$) correspond to sources with E(B - V) $> 3$, which are typically missed in UV-optical studies. This underscores the importance of deep near-infrared imaging for finding and characterizing these very young and dust-embedded sources., Comment: 18 pages, 7 figures, 1 table; Accepted for publication in ApJL

The nuclear stellar disc (NSD) is a flat and dense stellar structure at the centre of the Milky Way. Previous work has identified the presence of metal-rich and metal-poor stars in the NSD, suggesting that they have different origins. The recent publication of photometric, metallicity, proper motion, and orbital catalogues allows the NSD stellar population to be characterised with unprecedented detail. We aim to explore the proper motions and orbits of NSD stars with different metallicities to assess whether they have different origins and to better understand the metallicity distribution in the NSD. We distinguished between metal-rich and metal-poor stars by applying a Gaussian mixture model, as done in previous work, and analysed the proper motions, orbits, and spatial distribution of stars with different metallicities. We find that metal-rich stars exhibit a lower velocity dispersion, suggesting that they trace a kinematically cooler component compared to metal-poor ones. Furthermore, z-tube orbits are predominant among metal-rich stars, while chaotic/box orbits are more common among metal-poor ones. We also find that metal-rich and metal-poor stars show a similar extinction and are present throughout the analysed regions. As a secondary result, we detected a metallicity gradient in the metal-rich population with higher metallicity towards the centre of the NSD and a tentative gradient for the metal-poor stars, which is consistent with previous studies that did not distinguish between the two populations. Our results suggest that metal-rich stars trace the NSD, whereas metal-poor ones are related to the Galactic bar and probably constitute Galactic bar interlopers and/or are NSD stars that originated from accreted clusters. The detected metallicity gradients aligns with the currently accepted inside-out formation of the NSD., Comment: Accepted for publication in A&A. 11 pages, 12 figures

Building generalist robotic systems involves effectively endowing robots with the capabilities to handle novel objects in an open-world setting. Inspired by the advances of large pre-trained models, we propose Keypoint Affordance Learning from Imagined Environments (KALIE), which adapts pre-trained Vision Language Models (VLMs) for robotic control in a scalable manner. Instead of directly producing motor commands, KALIE controls the robot by predicting point-based affordance representations based on natural language instructions and visual observations of the scene. The VLM is trained on 2D images with affordances labeled by humans, bypassing the need for training data collected on robotic systems. Through an affordance-aware data synthesis pipeline, KALIE automatically creates massive high-quality training data based on limited example data manually collected by humans. We demonstrate that KALIE can learn to robustly solve new manipulation tasks with unseen objects given only 50 example data points. Compared to baselines using pre-trained VLMs, our approach consistently achieves superior performance., Comment: 8 pages, 7 figures

We present results of our dynamical stream modelling for the North West Stream in the outer halo of the Andromeda galaxy (M31). Comprising two main segments, the North West Stream was thought to be a single structured arching around M31. However, recent evidence suggests that it is two separate, unrelated, streams. To test this hypothesis we use observational data from 6 fields associated with the upper segment of the North West Stream together with 8 fields and 5 globular clusters associated with the lower segment to constrain model orbits. We fit both segments of the stream using a fixed potential model for M31 and an orbit integrator to compare orbits with the observed streams. We measure the central tracks and predict proper motions for for the upper segment (lower segment) finding ${\mu^*_{\alpha}}$ = 0.078$^{+0.015}_{-0.012}$ (0.085$^{+0.001}_{-0.002}$) mas/yr and ${\mu_{\delta}}$ = $-$0.05$^{+0.008}_{-0.009}$ ($-$0.095$^{+0.003}_{-0.005}$) mas/yr. Our results support the hypothesis that the dwarf spheroidal galaxy Andromeda XXVII is the progenitor of the upper segment of the North West Stream and that the upper and lower segments do not comprise a single structure. We propose that the upper segment, which appears to be on an infall trajectory with M31, be renamed the "Andromeda XXVII Stream" and the lower segment, also apparently infalling towards M31, retain the name "North West Stream"., Comment: 20 pages, 14 figures

We present a kinematic and spectroscopic analysis of 40 red giant branch stars, in 9 fields, exquisitely delineating the lower segment of the North West Stream (NW-K2), which extends for $\sim$80 kpc from the centre of the Andromeda galaxy. We measure the stream's systemic velocity as -439.3$^{+4.1}_{-3.8}$ km/s with a velocity dispersion = 16.4$^{+5.6}_{-3.8}$ km/s that is in keeping with its progenitor being a dwarf galaxy. We find no detectable velocity gradient along the stream. We determine $-$1.3$\pm$0.1 $\le$ <[Fe/H]$_{\rm spec}$> $\le$ $-$1.2$\pm$0.8 but find no metallicity gradient along the stream. We are able to plausibly associate NW-K2 with the globular clusters PandAS-04, PandAS-09, PAndAS-10, PAndAS-11, PandAS-12 but not with PandAS-13 or PandAS-15 which we find to be superimposed on the stream but not kinematically associated with it., Comment: 14 pages, 8 figures

Most star formation in the local Universe occurs in spiral galaxies, but their origin remains an unanswered question. Various theories have been proposed to explain the development of spiral arms, each predicting different spatial distributions of the interstellar medium. This study maps the star formation rate (SFR) and gas-phase metallicity of nine spiral galaxies with the TYPHOON survey to test two dominating theories: density wave theory and dynamic spiral theory. We discuss the environmental effects on our galaxies, considering reported environments and merging events. Taking advantage of the large field of view covering the entire optical disk, we quantify the fluctuation of SFR and metallicity relative to the azimuthal distance from the spiral arms. We find higher SFR and metallicity in the trailing edge of NGC~1365 (by 0.117~dex and 0.068~dex, respectively) and NGC~1566 (by 0.119~dex and 0.037~dex, respectively), which is in line with density wave theory. NGC~2442 shows a different result with higher metallicity (0.093~dex) in the leading edge, possibly attributed to an ongoing merging. The other six spiral galaxies show no statistically significant offset in SFR or metallicity, consistent with dynamic spiral theory. We also compare the behaviour of metallicity inside and outside the co-rotation radius (CR) of NGC~1365 and NGC~1566. We find comparable metallicity fluctuations near and beyond the CR of NGC~1365, indicating gravitational perturbation. NGC~1566 shows the greatest fluctuation near the CR, in line with the analytic spiral arms. Our work highlights that a combination of mechanisms explains the origin of spiral features in the local Universe., Comment: 22 pages, 10 figures, 4 tables, accepted for publication in MNRAS

Parametric strong lensing studies of galaxy clusters often display "misleading features". This is the case in the galaxy cluster Abell 370. Using strong lensing techniques, it has been described parametrically by a four dark matter clumps model and galaxy scale perturbers, as well as a significant external shear component, which physical origin remains a challenge. The dark matter distribution features a mass clump with no stellar counterpart and a significant offset between one of the dark matter clumps and its associated stellar counterpart. In this paper, based on BUFFALO data, we begin by revisiting this mass model. We find a four dark matter clumps solution which does not require any external shear and provides a slightly better RMS compared to previous models. Investigating further this new solution, we present a class of models which can accurately reproduce the strong lensing data, but whose parameters for the dark matter component are poorly constrained. We then develop a model where each large scale dark matter component must be associated with a stellar counterpart. This model is unable to reproduce the observational constraints with an RMS smaller than 2.3", and the parameters describing this dark matter component remain poorly constrained. Examining the total projected mass maps, we find a good agreement between the total mass and the stellar distribution, both being bimodal. We interpret the "misleading features" of the four dark matter clumps mass model and the failure of the three dark matter clumps mass model as being symptomatic of the lack of realism of a parametric description of the dark matter distribution, and encourage caution and criticism on the outputs of parametric strong lensing modelling. We briefly discuss the implications of our results for using Abell 370 as a gravitational telescope., Comment: Abstract largely abridged - A&A accepted on 25/10/2024

The price of oil can rise because of a disruption to supply or an increase in demand. The nature of the price change determines the dynamic effects. As Kilian (2009) put it: "not all oil price shocks are alike." Using the latest available data, we extend Kilian's (2009) analysis using the R ecosystem and provide more evidence for Kilian's (2009) conclusions. Inference based on unknown conditional heteroskedasticity strengthens the conclusions. With the updated shocks, we assess how a local economy responds to the global oil market, an application that is relevant to policymakers concerned with the transition away from fossil fuels., Comment: 20 pages, 8 figures

We investigated the effectiveness of a sustained and spiraled content literacy intervention that emphasizes building domain and topic knowledge schemas and vocabulary for elementary-grade students. The Model of Reading Engagement (MORE) intervention underscores thematic lessons that provide an intellectual structure for helping students connect new learning to a general schema in Grade 1 (animal survival), Grade 2 (scientific investigation of past events like dinosaur mass extinctions), and Grade 3 (scientific investigation of living systems). A total of 30 elementary schools (N = 2,870 students) were randomized to a treatment or control condition. In the treatment condition (i.e., full spiral curriculum), students participated in content literacy lessons from Grades 1 to 3 during the school year and wide reading of thematically related informational texts in the summer following Grades 1 and 2. In the control condition (i.e., partial spiral curriculum), students participated in lessons in only Grade 3. The Grade 3 lessons for both conditions were implemented online during the COVID-19 pandemic school year. Results reveal that treatment students outperformed control students on science vocabulary knowledge across all three grades. Furthermore, intent-to-treat analyses revealed positive transfer effects on Grade 3 science reading (ES = 0.14), domain-general reading comprehension (ES = 0.11), and mathematics achievement (ES = 0.12). Treatment impacts were sustained at 14-month follow-up on Grade 4 reading comprehension (ES = 0.12) and mathematics achievement (ES = 0.16). Findings indicate that a content literacy intervention that spirals topics and vocabulary across grades can improve students' long-term academic achievement outcomes.

We present volume-limited samples of cataclysmic variables (CVs) and AM CVn binaries jointly selected from SRG/eROSITA eRASS1 and \textit{Gaia} DR3 using an X-ray + optical color-color diagram (the ``X-ray Main Sequence"). This tool identifies all CV subtypes, including magnetic and low-accretion rate systems, in contrast to most previous surveys. We find 23 CVs, 3 of which are AM CVns, out to 150 pc in the Western Galactic Hemisphere. Our 150 pc sample is spectroscopically verified and complete down to $L_X = 1.3\times 10^{29} \;\textrm{erg s}^{-1}$ in the 0.2--2.3 keV band, and we also present CV candidates out to 300 pc and 1000 pc. We discovered two previously unknown systems in our 150 pc sample: the third nearest AM CVn and a magnetic period bouncer. We find the mean $L_X$ of CVs to be $\langle L_X \rangle \approx 4.6\times 10^{30} \;\textrm{erg s}^{-1}$, in contrast to previous surveys which yielded $\langle L_X \rangle \sim 10^{31}-10^{32} \;\textrm{erg s}^{-1}$. We construct X-ray luminosity functions that, for the first time, flatten out at $L_X\sim 10^{30} \; \textrm{erg s}^{-1}$. We find average number, mass, and luminosity densities of $\rho_\textrm{N, CV} = (3.7 \pm 0.7) \times 10^{-6} \textrm{pc}^{-3}$, $\rho_M = (5.0 \pm 1.0) \times 10^{-5} M_\odot^{-1}$, and $\rho_{L_X} = (2.3 \pm 0.4) \times 10^{26} \textrm{erg s}^{-1}M_\odot^{-1}$, respectively, in the solar neighborhood. Our uniform selection method also allows us to place meaningful estimates on the space density of AM CVns, $\rho_\textrm{N, AM CVn} = (5.5 \pm 3.7) \times 10^{-7} \textrm{pc}^{-3}$. Magnetic CVs and period bouncers make up $35\%$ and $25\%$ of our sample, respectively. This work, through a novel discovery technique, shows that the observed number densities of CVs and AM CVns, as well as the fraction of period bouncers, are still in tension with population synthesis estimates., Comment: Submitted to PASP, comments welcome

We present a high resolution chemical study of a representative sample of 21 luminous giant stars of Liller~1, a complex stellar system in the Galactic bulge, based on H band spectra acquired with the Near InfraRed Spectrograph at KeckII. 15 stars turn out to have a subsolar iron abundance and enhanced [$\alpha$/Fe] and [Al/Fe], likely old that formed early and quickly from gas mainly enriched by type~II supernovae, and 6 stars with supersolar iron and roughly solar-scaled [$\alpha$/Fe] and [Al/Fe], likely younger, thus formed at later epochs from gas also enriched by type~Ia supernovae. Moreover, both subpopulations show enhanced [N/Fe], as in the bulge field, about solar-scaled [V/Fe], and depletion of [C/Fe] and $^{12}$C/$^{13}$C with respect to the solar values, indicating the occurrence of significant mixing in the stellar interiors of these evolved stars. The current study has also made evident that the sub-solar subpopulation shows some structuring, and the presence of a third subcomponent with iron content and [$\alpha$/Fe] enhancement somewhat intermediate between the metal-poor and metal-rich main subpopulations, has been statistically assessed, providing the chemical signature of an extended star formation with multiple bursts and of some self-enrichment., Comment: Accepted for pubblication in A&A

Viscoelastic shear flows support additional chaotic states beyond simple Newtonian turbulence. In vanishing Reynolds number flows, the nonlinearity in the polymer evolution equation alone can sustain inertialess 'elastic' turbulence (ET) while 'elasto-inertial' turbulence (EIT) appears to rely on an interplay between elasticity and finite-$Re$ effects. Despite their distinct phenomenology and industrial significance, transition routes and possible connections between these states are unknown. We identify here a common Ruelle-Takens transition scenario for both of these chaotic regimes in two-dimensional direct numerical simulations of FENE-P fluids in a straight channel. The primary bifurcation is caused by a recently-discovered 'polymer diffusive instability' associated with small but non-vanishing polymer stress diffusion which generates a finite-amplitude, small-scale travelling wave localised at the wall. This is found to be unstable to a large-scale secondary instability which grows to modify the whole flow before itself breaking down in a third bifurcation to either ET or EIT. The secondary large-scale instability waves resemble 'centre' and 'wall' modes respectively - instabilities which have been conjectured to play a role in viscoelastic chaotic dynamics but were previously only thought to exist far from relevant areas of the parameter space., Comment: 6 Figs, main paper and supplementary material

We present intial results associating galaxies in the MUSE Ultra Deep Field (MUDF) with gas seen in absorption along the line-of-sight to two bright quasars in this field, to explore the dependence of metals in the circumgalactic medium (CGM) on galaxy properties. The MUDF includes $\sim$140h of VLT/MUSE data and 90 orbits of HST/G141M grism observations alongside VLT/UVES spectroscopy of the two quasars and several bands of HST imaging. We compare the metal absorption around galaxies in this field as a function of impact parameter, azimuthal angle and galaxy metallicity across redshifts 0.5 $<$ z $<$ 3.2. Due to the depth of our data and a large field-of-view, our analysis extends to low stellar masses ($<$ $10^{7}$ M$_{\odot}$) and high impact parameters ($>$ 600 kpc). We find a correlation between absorber equivalent width and number of nearby galaxies, but do not detect a significant anti-correlation with impact parameter. Our full sample does not show any significant change in absorber incidence as a function of azimuthal angle. However, we do find a bimodality in the azimuthal angle distribution of absorption at small impact parameters ($<$2 r$_{vir}$) and around highly-star-forming galaxies, possibly indicating disk-like accretion and biconical outflows. Finally, we do not detect any systematic deviation from the fundamental metallicity relation (FMR) among galaxies with detected absorption. This work is limited by gaps in the wavelength coverage of our current data; broader-wavelength observations with JWST will allow us to unlock the full potential of the MUDF for studying the CGM., Comment: Accepted for publication in The Astrophysical Journal. 31 pages, 17 figures

We present the properties of two candidate massive ($M_\star\sim10^{11}M_\odot$) and dusty ($A_{\rm v}>2.5$ mag) galaxies at $z=5-7$ in the first 0.28 deg$^2$ of the COSMOS-Web survey. One object is spectroscopically confirmed at $z_{\rm spec}=5.051$, while the other has a robust $z_{\rm phot}=6.7\pm0.3$. Thanks to their extremely red colors ($F277W-F444W\sim1.7$ mag), these galaxies satisfy the nominal color-selection for the widely-studied ``little red dot" (LRD) population with the exception of their spatially-resolved morphologies. The morphology of our targets allows us to conclude that their red continuum is dominated by highly obscured stellar emission and not by reddened nuclear activity. Using a variety of SED-fitting tools and star formation histories, we estimate the stellar masses to be $\log(M_\star)=11.32^{+0.07}_{-0.15}$ $M_\odot$ and $\log(M_\star)=11.2^{+0.1}_{-0.2}$ $M_\odot$, respectively, with a red continuum emission dominated by a recent episode of star formation. We then compare their number density to the halo mass function to infer stellar baryon fractions of $\epsilon_\star\sim0.25$ and $\epsilon_\star\sim0.5$. Both are significantly higher than what is commonly observed in lower-z galaxies or more dust-obscured galaxies at similar redshifts. With very bright ultra-high-z Lyman-Break Galaxies and some non-AGN dominated LRDs, such ``extended" LRDs represent another population that may require very efficient star formation at early times., Comment: 11 pages, 4 figures, 1 table. Accepted for publication in ApJL

We present new free-form and hybrid mass reconstructions of the galaxy cluster lens MACS J0416.1-2403 at $z=0.396$ using the lens inversion method GRALE. The reconstructions use 237 spectroscopically confirmed multiple images from Bergamini et. al. 2023 as the main input. Our primary model reconstructs images to a positional accuracy of 0.191", thus representing one of the most precise reconstructions of this lens to date. Our models find broad agreement with previous reconstructions, and identify two $\sim 10^{11} M_{\odot}$ light-unaffiliated substructures. We focus on two highly magnified arcs: Spock and Mothra. Our model features a unique critical curve structure around the Spock arc with 2 crossings. This structure enables sufficient magnification across this arc to potentially explain the large number of transients as microlensing events of supergiant stars. Additionally, we develop a model of the millilens substructure expected to be magnifying Mothra, which may be a binary pair of supergiants with $\mu \sim 6000$. This model accounts for flexibility in the millilens position while preserving the observed flux and minimizing image position displacements along the Mothra arc. We constrain the millilens mass and core radius to $\lesssim 10^6 M_{\odot}$ and $\lesssim 12$ pc, respectively, which would render it one of the smallest and most compact substructures constrained by lensing. If the millilens is dominated by wave dark matter, the axion mass is constrained to be $< 3.0 \times 10^{-21}$ eV. Further monitoring of this lens with JWST will uncover more transients, permitting tighter constraints on the structure surrounding these two arcs., Comment: Submitted. Comments Welcome

Astronomy is at a turning point in its history and in its relations with the Indigenous peoples who are the generational stewards of land where several of our main observatories are located. The controversy regarding the further development of astronomy facilities on Maunakea is probably the most significant and publicized conflict about the use of such land in the name of science. Thousands have stood in resistance, elders were arrested, and the community is divided. Astronomy's access to one of its most emblematic sites is at risk. This situation challenges our professional practice, the projects we build on Indigenous lands, and our relationships with the people who live within these lands and with society in general. This paper attempts to share the perspective of the authors on the historical events, including the very recent past, through the lens of our understanding and opinions; to provide transparency, with humility, into our process of introspection and transformation; and to share our hopes and ambitions as leaders from Maunakea Observatories for the future of astronomy in Hawai'i, as advocated by the Astro2020 report from the U.S. National Academies of Sciences, Engineering, and Medicine; and to suggest ways for the profession to commit to this long-term vision.

With a luminosity similar to that of Milky Way dwarf spheroidal (dSph) systems like Sextans, but a spatial extent similar to that of ultradiffuse galaxies (UDGs), Andromeda (And) XIX is an unusual satellite of M31. To investigate the origin of this galaxy, we measure chemical abundances for AndXIX derived from medium-resolution (R$\sim$6000) spectra from Keck II/DEIMOS. We coadd 79 red giant branch stars, grouped by photometric metallicity, in order to obtain a sufficiently high signal-to-noise ratio (S/N) to measure 20 [Fe/H] and [$\alpha$/Fe] abundances via spectral synthesis. The latter are the first such measurements for AndXIX. The mean metallicity we derive for AndXIX places it $\sim2\sigma$ higher than the present-day stellar mass-metallicity relation for Local Group dwarf galaxies, potentially indicating it has experienced tidal stripping. A loss of gas and associated quenching during such a process, which prevents the extended star formation necessary to produce shallow [$\alpha$/Fe]--[Fe/H] gradients in massive systems, is also consistent with the steeply decreasing [$\alpha$/Fe]--[Fe/H] trend we observe. In combination with the diffuse structure and disturbed kinematic properties of AndXIX, this suggests tidal interactions, rather than galaxy mergers, are strong contenders for its formation., Comment: 26 pages, 10 figures. Accepted by ApJ

The NOEMA formIng Cluster survEy (NICE) is a large program targeting 69 massive galaxy group candidates at $z>2$ in six deep fields. We report spectroscopic confirmation of eight groups at $1.65\leq z\leq3.61$ in COSMOS. Homogeneously selected as significant overdensities of red IRAC sources with red Herschel colors, four groups are confirmed by CO and [CI] with NOEMA 3mm observations, three are confirmed with ALMA, and one is confirmed by H$\alpha$ from Subaru/FMOS. We constructed the integrated FIR SEDs for the eight groups, obtaining total IR SFR $=260-1300~{\rm M_\odot}$~yr$^{-1}$. We adopted six methods to estimate the dark matter masses, including stellar mass to halo mass relations, overdensity with galaxy bias, and NFW profile fitting to radial stellar mass density. We found the radial stellar mass density are consistent with a NFW profile, supporting that they are collapsed structures hosted by a single dark matter halo. The best halo mass estimates are $\log(M_{\rm h}/{\rm M_\odot})=12.8-13.7$ with uncertainty of 0.3 dex. From halo mass estimates, we derive baryonic accretion rate ${\rm BAR}=(1-8)\times10^{3}\,{\rm M_{\odot}/yr}$ for this sample. We find a quasi-linear correlation between the integrated SFR/BAR and the theoretical halo mass limit for cold streams, $M_{\rm stream}/M_{\rm h}$, with ${\rm SFR/BAR}=10^{-0.46\pm0.22}\left({M_{\rm stream}/M_{\rm h}}\right)^{0.71\pm0.16}$ with a scatter of $0.40\,{\rm dex}$. Further, we compare halo masses and stellar masses with simulations, and find all structures are consistent with being progenitors of $M_{\rm h}(z=0)>10^{14}\,{\rm M_{\odot}}$ galaxy clusters, and the most massive central galaxies have stellar masses consistent with brightest cluster galaxies (BCGs) progenitors in the TNG300 simulation. The results strongly suggest these structures are forming massive galaxy clusters via baryonic and dark matter accretion., Comment: 44 pages (27pp appendix), 32 figures, 18 tables, accepted for publication in A&A

Traditional Monte Carlo methods for particle transport utilize source iteration to express the solution, the flux density, of the transport equation as a Neumann series. Our contribution is to show that the particle paths simulated within source iteration are associated with the adjoint flux density and the adjoint particle paths are associated with the flux density. We make our assertion rigorous through the use of stochastic calculus by representing the particle path used in source iteration as a solution to a stochastic differential equation (SDE). The solution to the adjoint Boltzmann equation is then expressed in terms of the same SDE and the solution to the Boltzmann equation is expressed in terms of the SDE associated with the adjoint particle process. An important consequence is that the particle paths used within source iteration simultaneously provide Monte Carlo approximations of the flux density and adjoint flux density. Monte Carlo simulations are presented to support the simultaneous use of the particle paths., Comment: Original submitted manuscript. Published manuscript differs in: new expanded related work section; generalized formulas for non-unit velocity; and text corrections for typos and clarity. Accepted manuscript will be posted on OSTI and this comment will be updated with that information once available. This submitted version contains 24 pages, 1 table, and 2 figures. SAND2024-08017J

We report results of intensive time-resolved imaging photometry and synoptic deep imaging of the comet C/2012 S1 (ISON) performed in February 2013. The data were obtained at the Wise Observatory in Israel (WO), at the Himalayan Chandra Telescope (HCT) in India, and at the Polaris Observatory Association in California, USA. During this period, the comet's heliocentric distance changed from 4.9 to 4.6 AU, just within the orbit of Jupiter. We analyze these early images in an attempt to determine the nuclear rotation period, assuming that at these relatively large heliocentric distances it would be possible to detect the photometric modulation of a rotating nucleus against an underdeveloped coma. Since this is not evident in our February 2013 data, with more than 400 independent photometric measurements analyzed, we can only set upper limits of 0.05 mag for periodic brightness modulations. We discuss (and discount) a possible brightening event (minor outburst) that occurred on $15-16$ February 2013. We also present deep synoptic images of the comet, obtained by combining our exposures for each night, and analyze them. We find that during the period of our observations the comet exhibited a $\sim$$30^{\prime\prime}\simeq 60000$-km tail with no substructures visible and that this appearance did not change throughout our campaign. The comet, as indicated by a single spectroscopic measurement obtained during this observation period, showed a dust coma reflecting the solar light. Our observations indicate that during February 2013, comet ISON was relatively quiet, with the dust coma presumably hiding any light modulation by a spinning nucleus., Comment: 15 pages, 10 figures

We report the first detection of coherent elastic neutrino-nucleus scattering (CEvNS) on germanium, measured at the Spallation Neutron Source at Oak Ridge National Laboratory. The Ge-Mini detector of the COHERENT collaboration employs large-mass, low-noise, high-purity germanium spectrometers, enabling excellent energy resolution, and an analysis threshold of 1.5 keV electron-equivalent ionization energy. We observe a on-beam excess of 20.6$_{+7.1}^{-6.3}$ counts with a total exposure of 10.22 GWhkg and we reject the no-CEvNS hypothesis with 3.9 sigma significance. The result agrees with the predicted standard model of particle physics signal rate within 2 sigma., Comment: 7 pages, 5 figures

We used near-infrared spectroscopy at medium-high resolution (R=8,000$-$25,000) to perform the first comprehensive chemical study of the intermediate luminosity bulge globular cluster Terzan~6. We derived detailed abundances and abundance patterns of 27 giant stars, likely members of Terzan~6, based on their accurate Hubble Space Telescope proper motions and line-of-sight radial velocities. From the spectral analysis of these stars, we determined an average heliocentric radial velocity of 143.3$\pm$1.0 km s$^{-1}$ with a velocity dispersion of 5.1$\pm$0.7 km s$^{-1}$ and an average [Fe/H]=$-0.65\pm0.01$ and a low 1$\sigma$ dispersion of 0.03 dex. We also measured some depletion of [Mn/Fe] with respect to the solar-scaled values and enhancement of for [Ca/Fe], [Si/Fe], [Mg/Fe], [Ti/Fe], [O/Fe], [Al/Fe], [Na/Fe], and, to a lower extent, for [K/Fe], consistent with previous measurements of other bulge globular clusters and favoring the scenario of a rapid bulge formation and chemical enrichment. Some spread in the light element abundances suggest the presence of first- and second-generation stars, typical of genuine globulars. Finally, we measured some depletion of carbon and low $\rm ^{12}C/^{13}C$ isotopic ratios, as in previous studies of field and cluster bulge giants, indicating that extra-mixing mechanisms should be at work during the post main sequence evolution in the high metallicity regime as well., Comment: 8 pages, 5 figures, accepted for publication in Astronomy and Astrophysics

This study examined which variables predicted the sense of presence (being there) in an immersive audio experience, with a focus on the impacts of immersion technology (headphones with spatialised sound versus speaker with 2D stereo sound), the nature of the audio experience, the narrative content, participants emotional and cognitive engagement and personal characteristics such as age and gender. Museum visitors listened to a story on the One Story, Many Voices immersive audio installation. A convergent mixed-methods design was used, including multiple regression analysis of survey data (n 185) and relational analysis of interview data (n 9). This study found mixed methods support from both surveys and interviews to suggest that presence was predicted by the audio quality (especially the audio being perceived as different from other audio stories they had listened to in the past), cognitive engagement (especially being totally absorbed), the narrative (the story and how it was told) and some support was found for emotional engagement (especially feeling a connection with the storyteller). Effects of the environment and immersion technology on presence received only either quantitative or qualitative support. No influences of personal characteristics were found. Findings are relevant for academics, immersive experience (XR) commissioners and developers.

We present BEAST DB, an open-source database comprised of ab initio electrochemical data computed using grand-canonical density functional theory in implicit solvent at consistent calculation parameters. The database contains over 20,000 surface calculations and covers a broad set of heterogeneous catalyst materials and electrochemical reactions. Calculations were performed at self-consistent fixed potential as well as constant charge to facilitate comparisons to the computational hydrogen electrode. This article presents common use cases of the database to rationalize trends in catalyst activity, screen catalyst material spaces, understand elementary mechanistic steps, analyze electronic structure, and train machine learning models to predict higher fidelity properties. Users can interact graphically with the database by querying for individual calculations to gain granular understanding of reaction steps or by querying for an entire reaction pathway on a given material using an interactive reaction pathway tool. BEAST DB will be periodically updated, with planned future updates to include advanced electronic structure data, surface speciation studies, and greater reaction coverage., Comment: 24 pages, 8 figures

The metallicity distribution function (MDF) of the Galactic bulge features a multi-peak shape, with a metal-poor peak at [Fe/H]=-0.3 dex and a metal-rich peak at [Fe/H]=+0.3 dex. This bimodality is also seen in [alpha/Fe] versus [Fe/H] ratios, indicating different stellar populations in the bulge. We aim to replicate the observed MDF by proposing a scenario where the metal-poor bulge stars formed in situ during an intense star formation burst, while the metal-rich stars formed during a second burst and/or were accreted from the inner Galactic disk due to a growing bar. We used a chemical evolution model that tracks various chemical species with detailed nucleosynthesis, focusing on Fe production from both Type Ia supernovae and massive stars, including rotating massive stars with varying velocities. Our model also accounts for gas infall, outflow, and the effect of stellar migration. Results are compared to 13,000 stars from the SDSS/APOGEE survey within 3.5 kpc of the Galactic center. Our model successfully reproduces the double-peak shape of the bulge MDF and the alpha-element abundance trends relative to Fe by assuming (i) a multi-burst star formation history with a 250 Myr quenching of the first burst and (ii) stellar migration from the inner disk due to a growing bar. We estimate that about 40% of the bulge-bar's stellar mass originates from the inner disk. Nucleosynthesis models that assume either no rotation for massive stars or a rotational velocity distribution favoring slow rotation at high metallicities best match the observed MDF and [alpha/Fe] and [Ce/Fe] versus [Fe/H] abundance patterns., Comment: Accepted for publication in A&A, 12 pages, 12 figures

We study of the role of galaxy-galaxy interactions and disk instabilities in producing starburst activity in galaxies out to z=4. For this, we use a sample of 387 galaxies with robust total star formation rate measurements from Herschel, gas masses from ALMA, stellar masses and redshifts from multi-band photometry, and JWST/NIRCam rest-frame optical imaging. Using mass-controlled samples, we find an increased fraction of interacting galaxies in the starburst regime at all redshifts out to z=4. This increase correlates with star formation efficiency (SFE), but not with gas fraction. However, the correlation is weak (and only significant out to z=2), which could be explained by the short duration of SFE increase during interaction. In addition, we find that isolated disk galaxies make up a significant fraction of the starburst population. The fraction of such galaxies with star-forming clumps ("clumpy disks") is significantly increased compared to the main-sequence disk population. Furthermore, this fraction directly correlates with SFE. This is direct observational evidence for a long-term increase of SFE maintained due to disk instabilities, contributing to the majority of starburst galaxies in our sample and hence to substantial mass growth in these systems. This result could also be of importance for explaining the growth of the most massive galaxies at z>6., Comment: 14 pages, 5 figures. Submitted to ApJ. Comments welcome

The field of robotics has made significant advances towards generalist robot manipulation policies. However, real-world evaluation of such policies is not scalable and faces reproducibility challenges, which are likely to worsen as policies broaden the spectrum of tasks they can perform. We identify control and visual disparities between real and simulated environments as key challenges for reliable simulated evaluation and propose approaches for mitigating these gaps without needing to craft full-fidelity digital twins of real-world environments. We then employ these approaches to create SIMPLER, a collection of simulated environments for manipulation policy evaluation on common real robot setups. Through paired sim-and-real evaluations of manipulation policies, we demonstrate strong correlation between policy performance in SIMPLER environments and in the real world. Additionally, we find that SIMPLER evaluations accurately reflect real-world policy behavior modes such as sensitivity to various distribution shifts. We open-source all SIMPLER environments along with our workflow for creating new environments at https://simpler-env.github.io to facilitate research on general-purpose manipulation policies and simulated evaluation frameworks.