Search

Benchmarks are important tools for tracking the rapid advancements in large language model (LLM) capabilities. However, benchmarks are not keeping pace in difficulty: LLMs now achieve over 90\% accuracy on popular benchmarks like MMLU, limiting informed measurement of state-of-the-art LLM capabilities. In response, we introduce Humanity's Last Exam (HLE), a multi-modal benchmark at the frontier of human knowledge, designed to be the final closed-ended academic benchmark of its kind with broad subject coverage. HLE consists of 3,000 questions across dozens of subjects, including mathematics, humanities, and the natural sciences. HLE is developed globally by subject-matter experts and consists of multiple-choice and short-answer questions suitable for automated grading. Each question has a known solution that is unambiguous and easily verifiable, but cannot be quickly answered via internet retrieval. State-of-the-art LLMs demonstrate low accuracy and calibration on HLE, highlighting a significant gap between current LLM capabilities and the expert human frontier on closed-ended academic questions. To inform research and policymaking upon a clear understanding of model capabilities, we publicly release HLE at https://lastexam.ai., Comment: 25 pages, 6 figures

We present progress on the calculation of scalar, vector, and tensor form factors for the following meson decays: $B\rightarrow\pi$, $B_s\rightarrow K$, $D\rightarrow\pi$ and $D_s\rightarrow K$. This calculation uses the MILC HISQ gluon field ensembles with HISQ valence quarks. We generate ensembles of correlator data with varying lattice spacings, some as small as 0.044 fm. Some ensembles have a strange-to-light quark mass ratio of 5:1 and others use the physical light quark mass. The fully-relativistic, heavy-HISQ approach is used for the heavy quark, with simulation masses ranging from the charm to near the bottom. This heavy-HISQ approach provides nearly full coverage of the kinematic range., Comment: 10 pages, 5 figures. Presentation on this work was given at the 41st International Symposium on Lattice Field Theory 28 July - 3 August 2024, Liverpool, UK. Submission ID: PoS(LATTICE2024)262

Out-of-equilibrium phases in many-body systems constitute a new paradigm in quantum matter - they exhibit dynamical properties that may otherwise be forbidden by equilibrium thermodynamics. Among these non-equilibrium phases are periodically driven (Floquet) systems [1-5], that are generically difficult to simulate classically due to their high entanglement. Using an array of superconducting qubits, we realize a Floquet topologically ordered state theoretically proposed in Ref. [6], image the characteristic dynamics of its chiral edge modes, and characterize its emergent anyonic excitations. Devising an interferometric algorithm allows us to introduce and measure a bulk topological invariant to probe the dynamical transmutation of anyons for system sizes up to 58 qubits. Our work demonstrates that quantum processors can provide key insights into the thus-far largely unexplored landscape of highly entangled non-equilibrium phases of matter., Comment: 20 pages, 22 figures

We analyze the asymptotic number and typical structure of claw-free graphs at constant edge densities. The first of our main results is a formula for the asymptotics of the logarithm of the number of claw-free graphs of edge density $\gamma \in (0,1)$. We show that the problem exhibits a second-order phase transition at edge density $\gamma^\ast=\frac{5-\sqrt{5}}{4}$. The asymptotic formula arises by solving a variational problem over graphons. For $\gamma\geq\gamma^\ast$ there is a unique optimal graphon, while for $\gamma<\gamma^\ast$ there is an infinite set of optimal graphons. By analyzing more detailed structure, we prove that for $\gamma<\gamma^\ast$, there is in fact a unique graphon $W$ such that almost all claw-free graphs at edge density $\gamma$ are close in cut metric to $W$. We also analyze the probability of claw-freeness in the Erd\H{o}s-R\'enyi random graph $G(n,p)$ for constant $p$, obtaining a formula for the large-deviation rate function for claw-freeness. In this case, the problem exhibits a first-order phase transition at $p^\ast=\frac{3-\sqrt{5}}{2}$, separating distinct structural regimes. At the critical point $p^\ast$, the corresponding graphon variational problem has infinitely many solutions, and we again pinpoint a unique optimal graphon that describes the typical structure of $G(n,p^\ast)$ conditioned on being claw-free.

We consider multi-agent systems with general information networks where an agent may only observe a subset of other agents. A system designer assigns local utility functions to the agents guiding their actions towards an outcome which determines the value of a given system objective. The aim is to design these local utility functions such that the Price of Anarchy (PoA), which equals the ratio of system objective at worst possible outcome to that at the optimal, is maximized. Towards this, we first develop a linear program (LP) that characterizes the PoA for any utility design and any information network. This leads to another LP that optimizes the PoA and derives the optimal utility design. Our work substantially generalizes existing approaches to the utility design problem. We also numerically show the robustness of proposed framework against unanticipated communication failures.

Mining and conveying actionable insights from complex data is a key challenge of exploratory data analysis (EDA) and storytelling. To address this challenge, we present a design space for actionable EDA and storytelling. Synthesizing theory and expert interviews, we highlight how semantic precision, rhetorical persuasion, and pragmatic relevance underpin effective EDA and storytelling. We also show how this design space subsumes common challenges in actionable EDA and storytelling, such as identifying appropriate analytical strategies and leveraging relevant domain knowledge. Building on the potential of LLMs to generate coherent narratives with commonsense reasoning, we contribute Jupybara, an AI-enabled assistant for actionable EDA and storytelling implemented as a Jupyter Notebook extension. Jupybara employs two strategies -- design-space-aware prompting and multi-agent architectures -- to operationalize our design space. An expert evaluation confirms Jupybara's usability, steerability, explainability, and reparability, as well as the effectiveness of our strategies in operationalizing the design space framework with LLMs., Comment: Conditionally Accepted at CHI '25

Contextual pricing strategies are prevalent in online retailing, where the seller adjusts prices based on products' attributes and buyers' characteristics. Although such strategies can enhance seller's profits, they raise concerns about fairness when significant price disparities emerge among specific groups, such as gender or race. These disparities can lead to adverse perceptions of fairness among buyers and may even violate the law and regulation. In contrast, price differences can incentivize disadvantaged buyers to strategically manipulate their group identity to obtain a lower price. In this paper, we investigate contextual dynamic pricing with fairness constraints, taking into account buyers' strategic behaviors when their group status is private and unobservable from the seller. We propose a dynamic pricing policy that simultaneously achieves price fairness and discourages strategic behaviors. Our policy achieves an upper bound of $O(\sqrt{T}+H(T))$ regret over $T$ time horizons, where the term $H(T)$ arises from buyers' assessment of the fairness of the pricing policy based on their learned price difference. When buyers are able to learn the fairness of the price policy, this upper bound reduces to $O(\sqrt{T})$. We also prove an $\Omega(\sqrt{T})$ regret lower bound of any pricing policy under our problem setting. We support our findings with extensive experimental evidence, showcasing our policy's effectiveness. In our real data analysis, we observe the existence of price discrimination against race in the loan application even after accounting for other contextual information. Our proposed pricing policy demonstrates a significant improvement, achieving 35.06% reduction in regret compared to the benchmark policy.

Distributed dataflow systems like Spark and Flink enable data-parallel processing of large datasets on clusters. Yet, selecting appropriate computational resources for dataflow jobs is often challenging. For efficient execution, individual resource allocations, such as memory and CPU cores, must meet the specific resource requirements of the job. An alternative to selecting a static resource allocation for a job execution is autoscaling as implemented for example by Spark. In this paper, we evaluate the resource efficiency of autoscaling batch data processing jobs based on resource demand both conceptually and experimentally by analyzing a new dataset of Spark job executions on Google Dataproc Serverless. In our experimental evaluation, we show that there is no significant resource efficiency gain over static resource allocations. We found that the inherent conceptual limitations of such autoscaling approaches are the inelasticity of node size as well as the inelasticity of the ratio of memory to CPU cores., Comment: 6 pages, 3 figures, 6 tables. IEEE Big Data 2024

Statistical experiments often seek to identify random variables with the largest population means. This inferential task, known as rank verification, has been well-studied on Gaussian data with equal variances. This work provides the first treatment of the unequal variances case, utilizing ideas from the selective inference literature. We design a hypothesis test that verifies the rank of the largest observed value without losing power due to multiple testing corrections. This test is subsequently extended for two procedures: Identifying some number of correctly-ordered Gaussian means, and validating the top-K set. The testing procedures are validated on NHANES survey data.

We introduce distributional dynamic programming (DP) methods for optimizing statistical functionals of the return distribution, with standard reinforcement learning as a special case. Previous distributional DP methods could optimize the same class of expected utilities as classic DP. To go beyond expected utilities, we combine distributional DP with stock augmentation, a technique previously introduced for classic DP in the context of risk-sensitive RL, where the MDP state is augmented with a statistic of the rewards obtained so far (since the first time step). We find that a number of recently studied problems can be formulated as stock-augmented return distribution optimization, and we show that we can use distributional DP to solve them. We analyze distributional value and policy iteration, with bounds and a study of what objectives these distributional DP methods can or cannot optimize. We describe a number of applications outlining how to use distributional DP to solve different stock-augmented return distribution optimization problems, for example maximizing conditional value-at-risk, and homeostatic regulation. To highlight the practical potential of stock-augmented return distribution optimization and distributional DP, we combine the core ideas of distributional value iteration with the deep RL agent DQN, and empirically evaluate it for solving instances of the applications discussed.

The increasing demand for larger and higher fidelity simulations has made Adaptive Mesh Refinement (AMR) and unstructured mesh techniques essential to focus compute effort and memory cost on just the areas of interest in the simulation domain. The distribution of these meshes over the compute nodes is often determined by balancing compute, memory, and network costs, leading to distributions with jagged nonconvex boundaries that fit together much like puzzle pieces. It is expensive, and sometimes impossible, to re-partition the data posing a challenge for in situ and post hoc visualization as the data cannot be rendered using standard sort-last compositing techniques that require a convex and disjoint data partitioning. We present a new distributed volume rendering and compositing algorithm, Approximate Puzzlepiece Compositing, that enables fast and high-accuracy in-place rendering of AMR and unstructured meshes. Our approach builds on Moment-Based Ordered-Independent Transparency to achieve a scalable, order-independent compositing algorithm that requires little communication and does not impose requirements on the data partitioning. We evaluate the image quality and scalability of our approach on synthetic data and two large-scale unstructured meshes on HPC systems by comparing to state-of-the-art sort-last compositing techniques, highlighting our approach's minimal overhead at higher core counts. We demonstrate that Approximate Puzzlepiece Compositing provides a scalable, high-performance, and high-quality distributed rendering approach applicable to the complex data distributions encountered in large-scale CFD simulations.

We study the long-term evolution of selected hierarchical triple systems in Newtonian gravity. We employ analytic equations derived in Paper II for the evolution of orbit-averaged orbital elements for both inner and outer orbits, which include two classes of contributions. One class consists of linear-order contributions, including quadrupole, octupole, hexadecapole and dotriacontapole orders, the latter scaling as $\epsilon^6$, where $\epsilon = a/A$, the ratio of the semimajor axes of the inner and outer orbits. The second class consists of contributions at {\em second} order in the fundamental perturbation parameter; they contribute at orders $\epsilon^{9/2}$, $\epsilon^{5}$, $\epsilon^{11/2}$ and $\epsilon^{6}$. For well studied triples such as star-planet systems perturbed by a low-mass third body (``hot Jupiters''), second-order and dotriacontapole (SOD) effects induce only small corrections. For stellar-mass binaries orbiting supermassive black holes, SOD corrections can suppress orbital flips that are generated by purely first-order effects. Planets orbiting binary star systems are susceptible to significant variations in the planetary semimajor axis, an effect that does not occur at first perturbative order. SOD effects in triple black hole systems can induce migrations of the eccentricity to significantly larger values than predicted by first-order perturbations, with implications for the gravitational-wave induced inspiral of the inner binary. We also show that in most cases, evolutions using our SOD equations are in better agreement with those from direct integration of the N-body equations of motion than those from first-order perturbations through hexadecapole order., Comment: 12 pages, 11 figures, to be submitted to the Physical Review D15

The closure principle is a standard tool for achieving family-wise error rate (FWER) control in multiple testing problems. In general, the computational cost for closed testing can be exponential in the number of hypotheses. The celebrated graphical approach of FWER control overcomes the computational hurdle by using weighted Bonferroni local tests on p-values with appropriately chosen weights. In this study, we extend the graphical approach to e-values. With valid e-values -- common in settings of sequential hypothesis testing or universal inference for irregular parametric models -- we can derive strictly more powerful local tests based on weighted averages of e-values. Consequently, this e-value-based closed test is more powerful than the corresponding graphical approach with inverse e-values as p-values. Although the computational shortcuts for the p-value-based graphical approach are not applicable, we develop efficient polynomial-time algorithms using dynamic programming for e-value-based graphical approaches with any directed acyclic graph. For special graphs, such as those used in the Holm's procedure and fallback procedure, we develop tailored algorithms with computation cost linear in the number of hypotheses, up to logarithmic factors., Comment: 16 pages, 5 figures

We develop double microwave shielding, which has recently enabled evaporative cooling to the first Bose-Einstein condensate of polar molecules [Bigagli et al., Nature 631, 289 (2024)]. Two microwave fields of different frequency and polarization are employed to effectively shield polar molecules from inelastic collisions and three-body recombination. Here, we describe in detail the theory of double microwave shielding. We demonstrate that double microwave shielding effectively suppresses two- and three-body losses. Simultaneously, dipolar interactions and the scattering length can be flexibly tuned, enabling comprehensive control over interactions in ultracold gases of polar molecules. We show that this approach works for a wide range of molecules. This opens the door to studying many-body physics with strongly interacting dipolar quantum matter.

Superconducting nanowire single-photon detectors (SNSPDs) are the current leading technology for the detection of single-photons in the near-infrared (NIR) and short-wave infrared (SWIR) spectral regions, due to record performance in terms of detection efficiency, low dark count rate, minimal timing jitter, and high maximum count rates. The various geometry and design parameters of SNSPDs are often carefully tailored to specific applications, resulting in challenges in optimising each performance characteristic without adversely impacting others. In particular, when scaling to larger array formats, the key challenge is to manage the heat load generated by the many readout cables in the cryogenic cooling system. Here we demonstrate a practical, self-contained 64-pixel SNSPD array system which exhibits high performance of all operational parameters, for use in the strategically important SWIR spectral region. The detector is an 8x8 array of 27.5 x 27.8 {\mu}m pixels on a 30 {\mu}m pitch, which leads to an 80 -- 85% fill factor. At a wavelength of 1550nm, a uniform average per-pixel photon detection efficiency of 77.7% was measured and the observed system detection efficiency (SDE) across the entire array was 65%. A full performance characterisation is presented, including a dark count rate of 20 cps per pixel, full-width-half-maximum (FWHM) jitter of 100 ps per pixel, a 3-dB maximum count rate of 645 Mcps and no evidence of crosstalk at the 0.1% level. This camera system therefore facilitates a variety of picosecond time-resolved measurement-based applications that include biomedical imaging, quantum communications, and long-range single-photon light detection and ranging (LiDAR) and 3D imaging., Comment: 12+5 pages, 4+5 figures

Entanglement distribution in quantum networks will enable next-generation technologies for quantum-secured communications, distributed quantum computing and sensing. Future quantum networks will require dense connectivity, allowing multiple users to share entanglement in a reconfigurable and multiplexed manner, while long-distance connections are established through the teleportation of entanglement, or entanglement swapping. While several recent works have demonstrated fully connected, local multi-user networks based on multiplexing, extending this to a global network architecture of interconnected local networks remains an outstanding challenge. Here we demonstrate the next stage in the evolution of multiplexed quantum networks: a prototype global reconfigurable network where entanglement is routed and teleported in a flexible and multiplexed manner between two local multi-user networks composed of four users each. At the heart of our network is a programmable 8x8-dimensional multi-port circuit that harnesses the natural mode-mixing process inside a multi-mode fibre to implement on-demand high-dimensional operations on two independent photons carrying eight transverse-spatial modes. Our circuit design allows us to break away from the limited planar geometry and bypass the control and fabrication challenges of conventional integrated photonic platforms. Our demonstration showcases the potential of this architecture for enabling large-scale, global quantum networks that offer versatile connectivity while being fully compatible with an existing communications infrastructure.

The ability to organically reason over and with both text and images is a pillar of human intelligence, yet the ability of Multimodal Large Language Models (MLLMs) to perform such multimodal reasoning remains under-explored. Existing benchmarks often emphasize text-dominant reasoning or rely on shallow visual cues, failing to adequately assess integrated visual and textual reasoning. We introduce EMMA (Enhanced MultiModal reAsoning), a benchmark targeting organic multimodal reasoning across mathematics, physics, chemistry, and coding. EMMA tasks demand advanced cross-modal reasoning that cannot be addressed by reasoning independently in each modality, offering an enhanced test suite for MLLMs' reasoning capabilities. Our evaluation of state-of-the-art MLLMs on EMMA reveals significant limitations in handling complex multimodal and multi-step reasoning tasks, even with advanced techniques like Chain-of-Thought prompting and test-time compute scaling underperforming. These findings underscore the need for improved multimodal architectures and training paradigms to close the gap between human and model reasoning in multimodality.

We develop the framework of Linear Simulation-based Inference (LSBI), an application of simulation-based inference where the likelihood is approximated by a Gaussian linear function of its parameters. We obtain analytical expressions for the posterior distributions of hyper-parameters of the linear likelihood in terms of samples drawn from a simulator, for both uniform and conjugate priors. This method is applied sequentially to several toy-models and tested on emulated datasets for the Cosmic Microwave Background temperature power spectrum. We find that convergence is achieved after four or five rounds of $\mathcal{O}(10^4)$ simulations, which is competitive with state-of-the-art neural density estimation methods. Therefore, we demonstrate that it is possible to obtain significant information gain and generate posteriors that agree with the underlying parameters while maintaining explainability and intellectual oversight., Comment: 12 pages, 7 figures

Context. The supernova remnant (SNR) W44 and its surroundings are a prime target for studying the acceleration of cosmic rays (CRs). Several previous studies established an extended gamma-ray emission that is set apart from the radio shell of W44. This emission is thought to originate from escaped high-energy CRs that interact with a surrounding dense molecular cloud complex. Aims. We present a detailed analysis of Fermi-LAT data with an emphasis on the spatial and spectral properties of W44 and its surroundings. We also report the results of the observations performed with the MAGIC telescopes of the northwestern region of W44. Finally, we present an interpretation model to explain the gamma-ray emission of the SNR and its surroundings. Methods. We first performed a detailed spatial analysis of 12 years of Fermi-LAT data at energies above 1 GeV, in order to exploit the better angular resolution, while we set a threshold of 100MeV for the spectral analysis. We performed a likelihood analysis of 174 hours of MAGIC data above 130 GeV using the spatial information obtained with Fermi-LAT. Results. The combined spectra of Fermi-LAT and MAGIC, extending from 100MeV to several TeV, were used to derive constraints on the escape of CRs. Using a time-dependent model to describe the particle acceleration and escape from the SNR, we show that the maximum energy of the accelerated particles has to be ' 40 GeV. However, our gamma-ray data suggest that a small number of lower-energy particles also needs to escape. We propose a novel model, the broken-shock scenario, to account for this effect and explain the gamma-ray emission.

Mrk 421 was in its most active state around early 2010, which led to the highest TeV gamma-ray flux ever recorded from any active galactic nuclei. We aim to characterize the multiwavelength behavior during this exceptional year for Mrk 421, and evaluate whether it is consistent with the picture derived with data from other less exceptional years. We investigated the period from November 5, 2009, (MJD 55140) until July 3, 2010, (MJD 55380) with extensive coverage from very-high-energy (VHE; E$\,>\,$100$\,$GeV) gamma rays to radio with MAGIC, VERITAS, Fermi-LAT, RXTE, Swift, GASP-WEBT, VLBA, and a variety of additional optical and radio telescopes. We investigated the variability and correlation behavior among different energy bands in great detail. We find the strongest variability in X-rays and VHE gamma rays, and PSDs compatible with power-law functions. We observe strong correlations between X-rays and VHE gamma rays. We also report a marginally significant positive correlation between high-energy (HE; E$\,>\,$100$\,$MeV) gamma rays and the ultraviolet band. We detected marginally significant correlations between the HE and VHE gamma rays, and between HE gamma rays and the X-ray, that disappear when the large flare in February 2010 is excluded from the correlation study. The activity of Mrk 421 also yielded the first ejection of features in the VLBA images of the jet of Mrk 421. Yet the large uncertainties in the ejection times of these radio features prevent us from firmly associating them to the specific flares recorded during the campaign. We also show that the collected multi-instrument data are consistent with a scenario where the emission is dominated by two regions, a compact and extended zone, which could be considered as a simplified implementation of an energy-stratified jet as suggested by recent IXPE observations., Comment: Accepted for publication in Astronomy & Astrophysics. Corresponding authors: Felix Schmuckermaier, David Paneque, Axel Arbet Engels

We introduce the Bayesian Global Sky Model (B-GSM), a novel data-driven Bayesian approach to modelling radio foregrounds at frequencies <400~MHz. B-GSM aims to address the limitations of previous models by incorporating robust error quantification and calibration. Using nested sampling, we compute Bayesian evidence and posterior distributions for the spectral behaviour and spatial amplitudes of diffuse emission components. Bayesian model comparison is used to determine the optimal number of emission components and their spectral parametrisation. Posterior sky predictions are conditioned on both diffuse emission and absolute temperature datasets, enabling simultaneous component separation and calibration. B-GSM is validated against a synthetic dataset designed to mimic the partial sky coverage, thermal noise, and calibration uncertainties present in real observations of the diffuse sky at low frequencies. B-GSM correctly identifies a model parametrisation with two emission components featuring curved power-law spectra. The posterior sky predictions agree with the true synthetic sky within statistical uncertainty. We find that the root-mean-square (RMS) residuals between the true and posterior predictions for the sky temperature as a function of LST are significantly reduced, when compared to the uncalibrated dataset. This indicates that B-GSM is able to correctly calibrate its posterior sky prediction to the independent absolute temperature dataset. We find that while the spectral parameters and component amplitudes exhibit some sensitivity to prior assumptions, the posterior sky predictions remain robust across a selection of different priors. This is the first of two papers, and is focused on validation of B-GSMs Bayesian framework, the second paper will present results of deployment on real data and introduce the low-frequency sky model which will be available for public download.

Quinn et al propose challenge datasets in their work called ``Kryptonite-N". These datasets aim to counter the universal function approximation argument of machine learning, breaking the notation that machine learning can ``approximate any continuous function" \cite{original_paper}. Our work refutes this claim and shows that universal function approximations can be applied successfully; the Kryptonite datasets are constructed predictably, allowing logistic regression with sufficient polynomial expansion and L1 regularization to solve for any dimension N.

Reinforcement learning from human feedback (RLHF) has become a cornerstone for aligning large language models with human preferences. However, the heterogeneity of human feedback, driven by diverse individual contexts and preferences, poses significant challenges for reward learning. To address this, we propose a Low-rank Contextual RLHF (LoCo-RLHF) framework that integrates contextual information to better model heterogeneous feedback while maintaining computational efficiency. Our approach builds on a contextual preference model, leveraging the intrinsic low-rank structure of the interaction between user contexts and query-answer pairs to mitigate the high dimensionality of feature representations. Furthermore, we address the challenge of distributional shifts in feedback through our Pessimism in Reduced Subspace (PRS) policy, inspired by pessimistic offline reinforcement learning techniques. We theoretically demonstrate that our policy achieves a tighter sub-optimality gap compared to existing methods. Extensive experiments validate the effectiveness of LoCo-RLHF, showcasing its superior performance in personalized RLHF settings and its robustness to distribution shifts.

When you have a question, the most effective way to have the question answered is to directly connect with experts on the topic and have a conversation with them. Prior to the invention of writing, this was the only way. Although effective, this solution exhibits scalability challenges. Writing allowed knowledge to be materialized, preserved, and replicated, enabling the development of different technologies over the centuries to connect information seekers with relevant information. This progression ultimately culminated in the ten-blue-links web search paradigm we're familiar with, just before the recent emergence of generative AI. However, we often forget that consuming static content is an imperfect solution. With the advent of large language models, it has become possible to develop a superior experience by allowing users to directly engage with experts. These interactions can of course satisfy information needs, but expert models can do so much more. This coming future requires reimagining search.

The BL Lacertae object VER J0521+211 underwent a notable flaring episode in February 2020. A short-term monitoring campaign, led by the MAGIC (Major Atmospheric Gamma Imaging Cherenkov) collaboration, covering a wide energy range from radio to very-high-energy (VHE, 100 GeV < E < 100 TeV) gamma rays was organised to study its evolution. These observations resulted in a consistent detection of the source over six consecutive nights in the VHE gamma-ray domain. Combining these nightly observations with an extensive set of multiwavelength data made modelling of the blazar's spectral energy distribution (SED) possible during the flare. This modelling was performed with a focus on two plausible emission mechanisms: i) a leptonic two-zone synchrotron-self-Compton scenario, and ii) a lepto-hadronic one-zone scenario. Both models effectively replicated the observed SED from radio to the VHE gamma-ray band. Furthermore, by introducing a set of evolving parameters, both models were successful in reproducing the evolution of the fluxes measured in different bands throughout the observing campaign. Notably, the lepto-hadronic model predicts enhanced photon and neutrino fluxes at ultra-high energies (E > 100 TeV). While the photon component, generated via decay of neutral pions, is not directly observable as it is subject to intense pair production (and therefore extinction) through interactions with the cosmic microwave background photons, neutrino detectors (e.g. IceCube) can probe the predicted neutrino component. Finally, the analysis of the gamma-ray spectra, as observed by MAGIC and the Fermi-LAT telescopes, yielded a conservative 95\% confidence upper limit of z \leq 0.244 for the redshift of this blazar., Comment: Accepted for publication in A&A

The cosmology of metric-affine gravity is studied for the general, parity preserving action quadratic in curvature, torsion and non-metricity. The model contains 27 a priori independent couplings in addition to the Einstein constant. Linear and higher order relations between the quadratic operators in a Friedmann--Lemaitre--Robertson--Walker spacetime are obtained, along with the modified Friedmann, torsion and non-metricity equations. Extra parameter constraints lead to two special branches of the model. Firstly, a branch is found in which the Riemannian spatial curvature (thought to be slightly closed or flat in the Lambda-CDM model of our Universe) is entirely screened from all the field equations, regardless of its true value. Secondly, an integrable branch is found which yields (anti) de Sitter expansion at late times. The particle spectra of these two branches are studied, and the need to eliminate higher-spin particles as well as ghosts and tachyons motivates further parameter constraints in each case. The most general model is also found which reproduces the exact Friedmann equations of general relativity. The full set of equations describing closed, open or flat cosmologies, for general parity-even quadratic metric-affine gravity, is made available for SymPy, Mathematica and Maple platforms., Comment: 23 pages, 8 figures

We prove that the set of limit points of the set of all spectral gaps of closed arithmetic hyperbolic surfaces equals $[0,\frac{1}{4}]$., Comment: 21 pages, 3 figures, 1 ancillary file

We give results on the asymptotic in Waring's problem over function fields that are stronger than the results obtained over the integers using the main conjecture in Vinogradov's mean value theorem. Similar estimates apply to Manin's conjecture for Fermat hypersurfaces over function fields. Following an idea of Pugin, rather than applying analytic methods to estimate the minor arcs, we treat them as complete exponential sums over finite fields and apply results of Katz, which bound the sum in terms of the dimension of a certain singular locus, which we estimate by tangent space calculations., Comment: 48 pages

Diffusion probabilistic models have achieved mainstream success in many generative modeling tasks, from image generation to inverse problem solving. A distinct feature of these models is that they correspond to deep hierarchical latent variable models optimizing a variational evidence lower bound (ELBO) on the data likelihood. Drawing on a basic connection between likelihood modeling and compression, we explore the potential of diffusion models for progressive coding, resulting in a sequence of bits that can be incrementally transmitted and decoded with progressively improving reconstruction quality. Unlike prior work based on Gaussian diffusion or conditional diffusion models, we propose a new form of diffusion model with uniform noise in the forward process, whose negative ELBO corresponds to the end-to-end compression cost using universal quantization. We obtain promising first results on image compression, achieving competitive rate-distortion and rate-realism results on a wide range of bit-rates with a single model, bringing neural codecs a step closer to practical deployment., Comment: 20 pages, 10 figures

Direct observation of proton emission with very small emission energy is often unfeasible due to the long partial half-lives associated with tunneling through the Coulomb barrier. Therefore proton emitters with very small Q-values may require masses of both parent and daughter nuclei to establish them as proton unbound. Nuclear mass models have been used to predict the proton drip-line of the thulium (Tm) isotopic chain ($Z=69$), but up until now the proton separation energy has not been experimentally tested. Mass measurements were therefore performed using a Multiple Reflection Time-Of-Flight Mass Spectrometer (MR-TOF-MS) at TRIUMF's TITAN facility to definitively map the limit of proton-bound Tm. The masses of neutron-deficient, $^{149}$Tm and $^{150}$Tm, combined with measurements of $^{149m,g}$Er (which were found to deviate from literature by $\sim$150 keV), provide the first experimental confirmation that $^{149}$Tm is the first proton-unbound nuclide in the Tm chain. Our measurements also enable the strength of the $N=82$ neutron shell gap to be determined at the Tm proton drip-line, providing evidence supporting its continued existence.

We present a novel approach to mesh shape editing, building on recent progress in 3D reconstruction from multi-view images. We formulate shape editing as a conditional reconstruction problem, where the model must reconstruct the input shape with the exception of a specified 3D region, in which the geometry should be generated from the conditional signal. To this end, we train a conditional Large Reconstruction Model (LRM) for masked reconstruction, using multi-view consistent masks rendered from a randomly generated 3D occlusion, and using one clean viewpoint as the conditional signal. During inference, we manually define a 3D region to edit and provide an edited image from a canonical viewpoint to fill in that region. We demonstrate that, in just a single forward pass, our method not only preserves the input geometry in the unmasked region through reconstruction capabilities on par with SoTA, but is also expressive enough to perform a variety of mesh edits from a single image guidance that past works struggle with, while being 10x faster than the top-performing competing prior work., Comment: Project Page: https://chocolatebiscuit.github.io/MaskedLRM/

We develop a ray-tracing model for laser-plasma interaction suitable for coupling in-line into kinetic particle-in-cell plasma simulation. The model is based on inverse Bremsstrahlung absorption and includes oblique incidence effects and reflection at the critical surface. The energy deposition is given to electrons by randomized kicks to momentum. The model is verified against analytic solutions and a 2-D laser ray-tracing code., Comment: 8 pages, 5 figures

A solution is provided to the Bruxelles Problem, a geometric decision problem originally posed in 1825, that asks for a synthetic construction to determine when ten points in 3-space lie on a quadric surface, a surface given by the vanishing of a degree-2 polynomial. The solution constructs four new points that are coplanar precisely when the ten original points lie on a quadric surface. The solution uses only lines constructed through two known points, planes constructed through three known points, and intersections of these objects. The tools involved include an extension of the Area Principle to three-dimensional space, bracket polynomials and the Grassmann-Cayley algebra, and von Staudt's results on geometric arithmetic. Many special cases are treated directly, leading to the generic case, where three pairs of the points generate skew lines and the remaining four points are in general position. A key step in the generic case involves finding a nice basis for the quadrics that pass through six of the ten points, which uses insights derived from Macaulay2, a computational algebra package not available in the nineteenth century., Comment: 14 pages, 4 figures

We introduce JPC, a JAX library for training neural networks with Predictive Coding. JPC provides a simple, fast and flexible interface to train a variety of PC networks (PCNs) including discriminative, generative and hybrid models. Unlike existing libraries, JPC leverages ordinary differential equation solvers to integrate the gradient flow inference dynamics of PCNs. We find that a second-order solver achieves significantly faster runtimes compared to standard Euler integration, with comparable performance on a range of tasks and network depths. JPC also provides some theoretical tools that can be used to study PCNs. We hope that JPC will facilitate future research of PC. The code is available at https://github.com/thebuckleylab/jpc., Comment: 9 pages, 7 figures

Atmospheric mass loss is thought to have strongly shaped the sample of close-in exoplanets. These atmospheres should be lost isotropically, leading to no net migration on the planetary orbit. However, strong stellar winds can funnel the escaping atmosphere into a tail trailing the planet. We derive a simple kinematic model of the gravitational interaction between the planet and this anisotropic wind, and derive expressions for the expected migration of the planet. Over the expected range of parameters, we find typical migrations of a few tenths to a few percent inward. We argue that this modest migration may be observable for planet pairs near mean motion resonances, which would provide an independent observational constraint on atmospheric mass loss models., Comment: Accepted to the Astronomical Journal on 18 Nov. 2024

We study the probability that the random graph $G(n,p)$ is triangle-free. When $p =o(n^{-1/2})$ or $p = \omega(n^{-1/2})$ the asymptotics of the logarithm of this probability are known via Janson's inequality in the former case and via regularity or hypergraph container methods in the latter case. We prove for the first time an asymptotic formula for the logarithm of this probability when $p = c n^{-1/2}$ for $c$ a sufficiently small constant. More generally, we study lower-tail large deviations for triangles in random graphs: the probability that $G(n,p)$ has at most $\eta$ times its expected number of triangles, when $p = c n^{-1/2}$ for $c$ and $\eta \in [0,1)$ constant. Our results apply for all $c$ if $\eta \ge .4993$ and for $c$ small enough otherwise. For $\eta$ small (including the case of triangle-freeness), we prove that a phase transition occurs as $c$ varies, in the sense of a non-analyticity of the rate function, while for $\eta \ge .4993$ we prove that no phase transition occurs. On the other hand for the random graph $G(n,m)$, with $m = b n^{3/2}$, we show that a phase transition occurs in the lower-tail problem for triangles as $b$ varies for \emph{every} $\eta \in [0,1)$. Our method involves ingredients from algorithms and statistical physics including the cluster expansion and concentration inequalities for contractive Markov chains.

The Airy line ensemble is a random collection of continuous ordered paths that plays an important role within random matrix theory and the Kardar-Parisi-Zhang universality class. The aim of this paper is to prove a universality property of the Airy line ensemble. We study growing numbers of i.i.d. continuous-time random walks which are then conditioned to stay in the same order for all time using a Doob h-transform. We consider a general class of increment distributions; a sufficient condition is the existence of an exponential moment and a log-concave density. We prove that the top particles in this system converge in an edge scaling limit to the Airy line ensemble in a regime where the number of random walks is required to grow slower than a certain power (with a non-optimal exponent 3/50) of the expected number of random walk steps. Furthermore, in a similar regime we prove that the law of large numbers and fluctuations of linear statistics agree with non-intersecting Brownian motions., Comment: 34 pages

There is an ever-growing need in the gravitational wave community for fast and reliable inference methods, accompanied by an informative error bar. Nested sampling satisfies the last two requirements, but its computational cost can become prohibitive when using the most accurate waveform models. In this paper, we demonstrate the acceleration of nested sampling using a technique called posterior repartitioning. This method leverages nested sampling's unique ability to separate prior and likelihood contributions at the algorithmic level. Specifically, we define a `repartitioned prior' informed by the posterior from a low-resolution run. To construct this repartitioned prior, we use a $\beta$-flow, a novel type of conditional normalizing flow designed to better learn deep tail probabilities. $\beta$-flows are trained on the entire nested sampling run and conditioned on an inverse temperature $\beta$. Applying our methods to simulated and real binary black hole mergers, we demonstrate how they can reduce the number of likelihood evaluations required for convergence by up to an order of magnitude, enabling faster model comparison and parameter estimation. Furthermore, we highlight the robustness of using $\beta$-flows over standard normalizing flows to accelerate nested sampling. Notably, $\beta$-flows successfully recover the same posteriors and evidences as traditional nested sampling, even in cases where standard normalizing flows fail., Comment: 12 pages, 13 figures

Super sample covariance (SSC) is important when estimating covariance matrices using a set of mock catalogues for galaxy surveys. If the underlying cosmological simulations do not include the variation in background parameters appropriate for the simulation sizes, then the scatter between mocks will be missing the SSC component. The coupling between large and small modes due to non-linear structure growth makes this pernicious on small scales. We compare different methods for generating ensembles of mocks with SSC built in to the covariance, and contrast against methods where the SSC component is computed and added to the covariance separately. We find that several perturbative expansions, developed to derive background fluctuations, give similar results. We then consider scaling covariance matrices calculated for simulations of different volumes to improve the accuracy of covariance matrix estimation for a given amount of computational time. On large scales, we find that the primary limitation is from the discrete number of modes contributing to the measured power spectrum, and we propose a new method for correcting this effect. Correct implementation of SSC and the effect of discrete mode numbers allows covariance matrices created from mocks to be scaled between volumes, potentially leading to a significant saving on computational resources when producing covariance matrices. We argue that a sub-percent match is difficult to achieve because of the effects of modes on scales between the box sizes, which cannot be easily included. Even so, a 3% match is achievable on scales of interest for current surveys scaling the simulation volume by 512x, costing a small fraction of the computational time of running full-sized simulations. This is comparable to the agreement between analytic and mock-based covariance estimates to be used with DESI Y1 results., Comment: 26 pages, 9 figures. For submission to JCAP

In this paper, we construct and analyze preconditioners for the interior penalty discontinuous Galerkin discretization posed in the space $H(\mathrm{div})$. These discretizations are used as one component in exactly divergence-free pressure-robust discretizations for the Stokes problem. Three preconditioners are presently considered: a subspace correction preconditioner using vertex patches and the lowest-order $H^1$-conforming space as a coarse space, a fictitious space preconditioner using the degree-$p$ discontinuous Galerkin space, and an auxiliary space preconditioner using the degree-$(p-1)$ discontinuous Galerkin space and a block Jacobi smoother. On certain classes of meshes, the subspace and fictitious space preconditioners result in provably well-conditioned systems, independent of the mesh size $h$, polynomial degree $p$, and penalty parameter $\eta$. All three preconditioners are shown to be robust with respect to $h$ on general meshes, and numerical results indicate that the iteration counts grow only mildly with respect to $p$ in the general case. Numerical examples illustrate the convergence properties of the preconditioners applied to structured and unstructured meshes. These solvers are used to construct block-diagonal preconditioners for the Stokes problem, which result in uniform convergence when used with MINRES., Comment: 24 pages, 1 figure

In this paper, we consider the problem of finding weak independent sets in a distributed network represented by a hypergraph. In this setting, each edge contains a set of r vertices rather than simply a pair, as in a standard graph. A k-weak independent set in a hypergraph is a set where no edge contains more than k vertices in the independent set. We focus two variations of this problem. First, we study the problem of finding k-weak maximal independent sets, k-weak independent sets where each vertex belongs to at least one edge with k vertices in the independent set. Second we introduce a weaker variant that we call (\alpha, \beta)-independent sets where the independent set is \beta-weak, and each vertex belongs to at least one edge with at least \alpha vertices in the independent set. Finally, we consider the problem of finding a (2, k)-ruling set on hypergraphs, i.e. independent sets where no vertex is a distance of more than k from the nearest member of the set. Given a hypergraph H of rank r and maximum degree \Delta, we provide a LLL formulation for finding an (\alpha, \beta)-independent set when (\beta - \alpha)^2 / (\beta + \alpha) \geq 6 \log(16 r \Delta), an O(\Delta r / (\beta - \alpha + 1) + \log^* n) round deterministic algorithm finding an (\alpha, \beta)-independent set, and a O(\Delta^2(r - k) \log r + \Delta \log r \log^* r + \log^* n) round algorithm for finding a k-weak maximal independent set. Additionally, we provide zero round randomized algorithms for finding (\alpha, \beta) independent sets, when (\beta - \alpha)^2 / (\beta + \alpha) \geq 6 c \log n + 6 for some constant c, and finding an m-weak independent set for some m \geq r / 2k where k is a given parameter. Finally, we provide lower bounds of \Omega(\Delta + \log^* n) and \Omega(r + \log^* n) on the problems of finding a k-weak maximal independent sets for some values of k.

The estimation of uncertainties in cosmological parameters is an important challenge in Large-Scale-Structure (LSS) analyses. For standard analyses such as Baryon Acoustic Oscillations (BAO) and Full Shape, two approaches are usually considered. First: analytical estimates of the covariance matrix use Gaussian approximations and (nonlinear) clustering measurements to estimate the matrix, which allows a relatively fast and computationally cheap way to generate matrices that adapt to an arbitrary clustering measurement. On the other hand, sample covariances are an empirical estimate of the matrix based on en ensemble of clustering measurements from fast and approximate simulations. While more computationally expensive due to the large amount of simulations and volume required, these allow us to take into account systematics that are impossible to model analytically. In this work we compare these two approaches in order to enable DESI's key analyses. We find that the configuration space analytical estimate performs satisfactorily in BAO analyses and its flexibility in terms of input clustering makes it the fiducial choice for DESI's 2024 BAO analysis. On the contrary, the analytical computation of the covariance matrix in Fourier space does not reproduce the expected measurements in terms of Full Shape analyses, which motivates the use of a corrected mock covariance for DESI's Full Shape analysis., Comment: 23 pages, 5 figures 7 tables, submitted to JCAP

Algorithms with predictions is a recent framework for decision-making under uncertainty that leverages the power of machine-learned predictions without making any assumption about their quality. The goal in this framework is for algorithms to achieve an improved performance when the predictions are accurate while maintaining acceptable guarantees when the predictions are erroneous. A serious concern with algorithms that use predictions is that these predictions can be biased and, as a result, cause the algorithm to make decisions that are deemed unfair. We show that this concern manifests itself in the classical secretary problem in the learning-augmented setting -- the state-of-the-art algorithm can have zero probability of accepting the best candidate, which we deem unfair, despite promising to accept a candidate whose expected value is at least $\max\{\Omega (1) , 1 - O(\epsilon)\}$ times the optimal value, where $\epsilon$ is the prediction error. We show how to preserve this promise while also guaranteeing to accept the best candidate with probability $\Omega(1)$. Our algorithm and analysis are based on a new "pegging" idea that diverges from existing works and simplifies/unifies some of their results. Finally, we extend to the $k$-secretary problem and complement our theoretical analysis with experiments., Comment: NeurIPS 2024

Pirogov--Sinai theory is a well-developed method for understanding the low-temperature phase diagram of statistical mechanics models on lattices. Motivated by physical and algorithmic questions beyond the setting of lattices, we develop a combinatorially flexible version of Pirogov--Sinai theory for the hard-core model of independent sets. Our results illustrate that the main conclusions of Pirogov--Sinai theory can be obtained in significantly greater generality than that of $\mathbb Z^{d}$. The main ingredients in our generalization are combinatorial and involve developing appropriate definitions of contours based on the notion of cycle basis connectivity. This is inspired by works of Tim\'{a}r and Georgakopoulos--Panagiotis.

Optical tweezer arrays have emerged as a key experimental platform for quantum computation, quantum simulation, and quantum metrology, enabling unprecedented levels of control over single atoms and molecules. Existing methods to generate tweezer arrays mostly rely on active beam-shaping devices, such as acousto-optic deflectors or liquid-crystal spatial light modulators. However, these approaches have fundamental limitations in array geometry, size, and scalability. Here we demonstrate the trapping of single atoms in optical tweezer arrays generated via holographic metasurfaces. We realize two-dimensional arrays with more than 250 tweezer traps, arranged in arbitrary geometries with trap spacings as small as 1.5 um. The arrays have a high uniformity in terms of trap depth, trap frequency, and positional accuracy, rivaling or exceeding existing approaches. Owing to sub-micrometer pixel sizes and high pixel densities, holographic metasurfaces open a path towards optical tweezer arrays with >100,000 traps, facilitating tweezer-array based quantum applications that require large system sizes., Comment: 15 pages, 5 main figures, 8 extended data figures

This study explores the electronic structure of the CH$_2$ molecule, modeled as a (6e, 23o) system using a 52-qubit quantum experiment, which is relevant for interstellar and combustion chemistry. We focused on calculating the dissociation energies for CH$_2$ in the ground state triplet and the first excited state singlet, applying the Sample-based Quantum Diagonalization (SQD) method within a quantum-centric supercomputing framework. We evaluated the ability of SQD to provide accurate results compared to Selected Configuration Interaction (SCI) calculations and experimental values for the singlet-triplet gap. To our knowledge, this is the first study of an open-shell system, such as the CH$_2$ triplet, using SQD. To obtain accurate energy values, we implemented post-SQD orbital optimization and employed a warm-start approach using previously converged states. While the results for the singlet state dissociation were only a few milli-Hartrees from the SCI reference values, the triplet state exhibited greater variability. This discrepancy likely arises from differences in bit-string handling within the SQD method for open- versus closed-shell systems, as well as the inherently complex wavefunction character of the triplet state. The SQD-calculated singlet-triplet energy gap matched well with experimental and SCI values. This study enhances our understanding of the SQD method for open-shell systems and lays the groundwork for future applications in large-scale electronic structure studies using quantum algorithms.

This paper considers cooperative control of robots involving two different testbed systems in remote locations with communication on the internet. This provides us the capability to exchange robots status like positions, velocities and directions needed for the swarming algorithm. The results show that all robots properly follow some leader defined one of the testbeds. Measurement of data exchange rates show no loss of packets, and average transfer delays stay within tolerance limits for practical applications. In our knowledge, the novelty of this paper concerns this kind of control over a large network like internet.

Spherically symmetric Einstein-{\ae}ther (E{\AE}) theory with a Maxwell-like kinetic term is revisited. We consider a general choice of the metric and the \ae{}ther field, finding that:~(i) there is a gauge freedom allowing one always to use a diagonal metric; and~(ii) the nature of the Maxwell equation forces the \ae{}ther field to be time-like in the coordinate basis. We derive the vacuum solution and confirm that the innermost stable circular orbit (ISCO) and photon ring are enlarged relative to general relativity (GR). Buchdahl's theorem in E\AE{} theory is derived. For a uniform physical density, we find that the upper bound on compactness is always lower than in GR. Additionally, we observe that the Newtonian and E\AE{} radial acceleration relations run parallel in the low pressure limit. Our analysis of E\AE{} theory may offer novel insights into its interesting phenomenological generalization: \AE{}ther--scalar--tensor theory ({\AE}ST)., Comment: 15 pages, 4 figures

Next-generation gravitational wave detectors are expected to detect millions of compact binary mergers across cosmological distances. The features of the mass distribution of these mergers, combined with gravitational wave distance measurements, will enable precise cosmological inferences, even without the need for electromagnetic counterparts. However, achieving accurate results requires modeling the mass spectrum, particularly considering possible redshift evolution. Binary neutron star (BNS) mergers are thought to be less influenced by changes in metallicity compared to binary black holes (BBH) or neutron star-black hole (NSBH) mergers. This stability in their mass spectrum over cosmic time reduces the chances of introducing biases in cosmological parameters caused by redshift evolution. In this study, we use the population synthesis code COMPAS to generate astrophysically motivated catalogs of BNS mergers and explore whether assuming a non-evolving BNS mass distribution with redshift could introduce biases in cosmological parameter inference. Our findings show that despite significant variations in the BNS mass distribution across binary physics assumptions and initial conditions in COMPAS, the joint mass-redshift population can be expressed as the product of the mass distribution marginalized over redshift and the redshift distribution marginalized over masses. This enables a 2% unbiased constraint on the Hubble constant - sufficient to address the Hubble tension. Additionally, we show that in the fiducial COMPAS setup, the bias from a non-evolving BNS mass model is less than 0.5% for the Hubble parameter measured at redshift 0.4. These results establish BNS mergers as strong candidates for spectral siren cosmology in the era of next-generation gravitational wave detectors., Comment: Submitted to ApJL 16 pages, 5 figures. Associated Zenodo link: https://zenodo.org/records/14704635 and GitHub link: https://github.com/SoumendraRoy/RedevolBNS

Distributed training of large models consumes enormous computation resources and requires substantial engineering efforts to compose various training techniques. This paper presents SimpleFSDP, a PyTorch-native compiler-based Fully Sharded Data Parallel (FSDP) framework, which has a simple implementation for maintenance and composability, allows full computation-communication graph tracing, and brings performance enhancement via compiler backend optimizations. SimpleFSDP's novelty lies in its unique $torch.compile$-friendly implementation of collective communications using existing PyTorch primitives, namely parametrizations, selective activation checkpointing, and DTensor. It also features the first-of-its-kind intermediate representation (IR) nodes bucketing and reordering in the TorchInductor backend for effective computation-communication overlapping. As a result, users can employ the aforementioned optimizations to automatically or manually wrap model components for minimal communication exposure. Extensive evaluations of SimpleFSDP on Llama 3 models (including the ultra-large 405B) using TorchTitan demonstrate up to 28.54% memory reduction and 68.67% throughput improvement compared to the most widely adopted FSDP2 eager framework, when composed with other distributed training techniques.