Search

Recent detections of gravitational waves have made black hole quasinormal modes a powerful tool in testing predictions of general relativity. Understanding the spectrum of these quasinormal modes in a broad class of theories beyond general relativity and a variety of astrophysical environments around black holes remains vital. In this work, we study the quasinormal mode spectrum of parametrized deformations of a non-rotating black hole in the vacuum. Following Vigeland and Hughes, we model these parametrized deformations as axisymmetric multipole moments in the Weyl coordinates with amplitudes much less than the amplitude of the Schwarzschild potential. These tiny bumps in the black hole geometry satisfy the linearized vacuum Einstein equations and are asymptotically flat. We use the recently developed modified Teukolsky formalism to derive one decoupled differential equation for the radiative Weyl scalar $\Psi_0$. We then use the eigenvalue perturbation method to compute the quasinormal mode frequency shifts of both even- and odd-parity modes with $\ell=2,3$ and up to the overtone number $n=2$ for the Weyl multipoles with $\ell_W=2,3$. Our calculation provides an avenue to directly connect the multipole moments of a modified black hole spacetime to the QNM frequency shifts in a parametric way., Comment: 18 pages, 5 figures, 4 tables

The Primordial Inflation Explorer (PIXIE) is an Explorer-class mission concept to measure the energy spectrum and linear polarization of the cosmic microwave background (CMB). A single cryogenic Fourier transform spectrometer compares the sky to an external blackbody calibration target, measuring the Stokes I, Q, U parameters to levels ~200 Jy/sr in each 2.65 degree diameter beam over the full sky, in each of 300 frequency channels from 28 GHz to 6 THz. With sensitivity over 1000 times greater than COBE/FIRAS, PIXIE opens a broad discovery space for the origin, contents, and evolution of the universe. Measurements of small distortions from a CMB blackbody spectrum provide a robust determination of the mean electron pressure and temperature in the universe while constraining processes including dissipation of primordial density perturbations, black holes, and the decay or annihilation of dark matter. Full-sky maps of linear polarization measure the optical depth to reionization at nearly the cosmic variance limit and constrain models of primordial inflation. Spectra with sub-percent absolute calibration spanning microwave to far-IR wavelengths provide a legacy data set for analyses including line intensity mapping of extragalactic emission and the cosmic infrared background amplitude and anisotropy. We describe the PIXIE instrument sensitivity, foreground subtraction, and anticipated science return from both the baseline 2-year mission and a potential extended mission., Comment: 44 pages with 18 figures and 4 tables, including full sensitivity curves. Sensitivity curves are available for download at https://lambda.gsfc.nasa.gov/product/pixie/pixie_baseline_noise_get.html

Continuous-time quantum walks (CTQWs) on dynamic graphs, referred to as dynamic CTQWs, are a recently introduced universal model of computation that offers a new paradigm in which to envision quantum algorithms. In this work we develop a mapping from dynamic CTQWs to the gate model of computation in the form of an algorithm to convert arbitrary single edge walks and single self loop walks, which are the fundamental building blocks of dynamic CTQWs, to their circuit model counterparts. We use this mapping to introduce an arbitrary quantum state preparation framework based on dynamic CTQWs. Our approach utilizes global information about the target state, relates state preparation to finding the optimal path in a graph, and leads to optimizations in the reduction of controls that are not as obvious in other approaches. Interestingly, classical optimization problems such as the minimal hitting set, minimum spanning tree, and shortest Hamiltonian path problems arise in our framework. We test our methods against uniformly controlled rotations methods, used by Qiskit, and find ours requires fewer CX gates when the target state has a polynomial number of non-zero amplitudes., Comment: Comments are welcome!

We present high-cadence photometric and spectroscopic observations of supernova (SN) 2024ggi, a Type II SN with flash spectroscopy features which exploded in the nearby galaxy NGC 3621 at $\sim$7 Mpc. The light-curve evolution over the first 30 hours can be fit by two power law indices with a break after 22 hours, rising from $M_V \approx -12.95$ mag at +0.66 days to $M_V \approx -17.91$ mag after 7 days. In addition, the densely sampled color curve shows a strong blueward evolution over the first few days and then behaves as a normal SN II with a redward evolution as the ejecta cool. Such deviations could be due to interaction with circumstellar material (CSM). Early high- and low-resolution spectra clearly show high-ionization flash features from the first spectrum to +3.42 days after the explosion. From the high-resolution spectra, we calculate the CSM velocity to be 37 $\pm~4~\mathrm{km\,s^{-1}} $. We also see the line strength evolve rapidly from 1.22 to 1.49 days in the earliest high-resolution spectra. Comparison of the low-resolution spectra with CMFGEN models suggests that the pre-explosion mass-loss rate of SN 2024ggi falls in a range of $10^{-3}$ to $10^{-2}$ M$_{\odot}$ yr$^{-1}$, which is similar to that derived for SN 2023ixf. However, the rapid temporal evolution of the narrow lines in the spectra of SN 2024ggi ($R_\mathrm{CSM} \sim 2.7 \times 10^{14} \mathrm{cm}$) could indicate a smaller spatial extent of the CSM than in SN 2023ixf ($R_\mathrm{CSM} \sim 5.4 \times 10^{14} \mathrm{cm}$) which in turn implies lower total CSM mass for SN 2024ggi., Comment: 23 pages, 15 figures, 4 tables, accepted for publication in ApJL

In a social dilemma situation, where individual and collective interests are in conflict, it sounds a reasonable assumption that the presence of super or smart players, who simultaneously punish defection and reward cooperation without allowing exploitation, could solve the basic problem. The behavior of such a multi-strategy system, however, is more subtle than it is firstly anticipated. When exploring the complete parameter space, we find that the emergence of cyclic dominance among strategies is rather common, which results in several counter-intuitive phenomena. For example, the defection level can be lowered at higher temptation, or weaker punishment provides better conditions for smart players. Our study indicates that smart cooperators can unexpectedly thrive under high temptation, emphasizing the complexity of strategic interactions. This study suggests that the principles governing these interactions can be applied to other moral behaviors, such as truth-telling and honesty, providing valuable insights for future research in multi-agent systems., Comment: 10 pages, 5 figures, accepted by Applied Mathematics and Computation

We consider random walks on edge coloured random graphs, where the colour of an edge reflects the cost of using it. In the simplest instance, the edges are coloured red or blue. Blue edges are free to use, whereas red edges incur a unit cost every time they are traversed.

We present a hardware-reconfigurable ansatz on $N_q$-qubits for the variational preparation of many-body states of the Anderson impurity model (AIM) with $N_{\text{imp}}+N_{\text{bath}}=N_q/2$ sites, which conserves total charge and spin z-component within each variational search subspace. The many-body ground state of the AIM is determined as the minimum over all minima of $O(N_q^2)$ distinct charge-spin sectors. Hamiltonian expectation values are shown to require $\omega(N_q) < N_{\text{meas.}} \leq O(N_{\text{imp}}N_{\text{bath}})$ symmetry-preserving, parallelizable measurement circuits, each amenable to post-selection. To obtain the one-particle impurity Green's function we show how initial Krylov vectors can be computed via mid-circuit measurement and how Lanczos iterations can be computed using the symmetry-preserving ansatz. For a single-impurity Anderson model with a number of bath sites increasing from one to six, we show using numerical emulation that the ease of variational ground-state preparation is suggestive of linear scaling in circuit depth and sub-quartic scaling in optimizer complexity. We therefore expect that, combined with time-dependent methods for Green's function computation, our ansatz provides a useful tool to account for electronic correlations on early fault-tolerant processors. Finally, with a view towards computing real materials properties of interest like magnetic susceptibilities and electron-hole propagators, we provide a straightforward method to compute many-body, time-dependent correlation functions using a combination of time evolution, mid-circuit measurement-conditioned operations, and the Hadamard test.

A rainbow stacking of $r$-edge-colorings $\chi_1, \ldots, \chi_m$ of the complete graph on $n$ vertices is a way of superimposing $\chi_1, \ldots, \chi_m$ so that no edges of the same color are superimposed on each other. We determine a sharp threshold for $r$ (as a function of $m$ and $n$) governing the existence and nonexistence of rainbow stackings of random $r$-edge-colorings $\chi_1,\ldots,\chi_m$., Comment: 12 pages, 2 figures

Epidemics are inherently stochastic, and stochastic models provide an appropriate way to describe and analyse such phenomena. Given temporal incidence data consisting of, for example, the number of new infections or removals in a given time window, a continuous-time discrete-valued Markov process provides a natural description of the dynamics of each model component, typically taken to be the number of susceptible, exposed, infected or removed individuals. Fitting the SEIR model to time-course data is a challenging problem due incomplete observations and, consequently, the intractability of the observed data likelihood. Whilst sampling based inference schemes such as Markov chain Monte Carlo are routinely applied, their computational cost typically restricts analysis to data sets of no more than a few thousand infective cases. Instead, we develop a sequential inference scheme that makes use of a computationally cheap approximation of the most natural Markov process model. Crucially, the resulting model allows a tractable conditional parameter posterior which can be summarised in terms of a set of low dimensional statistics. This is used to rejuvenate parameter samples in conjunction with a novel bridge construct for propagating state trajectories conditional on the next observation of cumulative incidence. The resulting inference framework also allows for stochastic infection and reporting rates. We illustrate our approach using synthetic and real data applications., Comment: 27 pages

Femtosecond-fast and nanometre-size pulses of electrons are emerging as unique probes for ultrafast dynamics at the nanoscale. Presently, such pulses are achievable only in highly sophisticated ultrafast electron microscopes or equally complex setups involving few-cycle-pulsed lasers with stable carrier-envelope phase (CEP) and nanotip probes. Here, we show that the generation of femtosecond pulses of nanoscale tunnelling electrons can be achieved in any ultrafast optical laboratory, using any (deep-UV to mid-IR) femtosecond laser in combination with photosensitive asymmetric nanogap (PAN) diodes fabricated via easy-to-scale adhesion lithography. The dominant mechanism producing tunnelling electrons in PANs is strong-field emission, which is easily achievable without CEP locking or external bias voltage. We employ PANs to demonstrate ultrafast nanoscopy of metal-halide perovskite quantum dots immobilised inside a 10-nm Al/Au nanogap and to characterise laser pulses across the entire optical region (266-6700 nm). Short electron pulses in PANs open the way towards scalable on-chip femtosecond electron measurements and novel design approaches for integrated ultrafast sensing nanodevices.

We study free graviton entanglement between Rindler wedges in the Minkowski vacuum state via the Euclidean path integral. We follow Kabat's method for computing the conical entropy, using the heat kernel on the cone with the tip removed, whose resulting von Neumann entropy for photons correctly predicted electromagnetic edge modes. We find that, in addition to the bulk graviton contributions, the conical entropy has a contact term that can be attributed to a vector field anchored to the (d-2)-dimensional (Euclidean) Rindler horizon whose contribution equals d-2 times Kabat's contact term for photons. We suggest that graviton edge modes are hence the d-2 large diffeomorphisms which act internally within the Rindler horizon. Along the way, we address several known issues regarding graviton entanglement. We furthermore sketch how our results may be used to study edge modes in closed bosonic string theory., Comment: 31 pages + appendices

While research of reinforcement learning applied to financial markets predominantly concentrates on finding optimal behaviours, it is worth to realize that the reinforcement learning returns $G_t$ and state value functions themselves are of interest and play a pivotal role in the evaluation of assets. Instead of focussing on the more complex task of finding optimal decision rules, this paper studies and applies the power of distributional state value functions in the context of financial market valuation and machine learning based trading algorithms. Accurate and trustworthy estimates of the distributions of $G_t$ provide a competitive edge leading to better informed decisions and more optimal behaviour. Herein, ideas from predictive knowledge and deep reinforcement learning are combined to introduce a novel family of models called CDG-Model, resulting in a highly flexible framework and intuitive approach with minimal assumptions regarding underlying distributions. The models allow seamless integration of typical financial modelling pitfalls like transaction costs, slippage and other possible costs or benefits into the model calculation. They can be applied to any kind of trading strategy or asset class. The frameworks introduced provide concrete business value through their potential in market valuation of single assets and portfolios, in the comparison of strategies as well as in the improvement of market timing. They can positively impact the performance and enhance the learning process of existing or new trading algorithms. They are of interest from a scientific point-of-view and open up multiple areas of future research. Initial implementations and tests were performed on real market data. While the results are promising, applying a robust statistical framework to evaluate the models in general remains a challenge and further investigations are needed., Comment: 22 pages, 13 figures

Recently, there has been a surge of international interest in extraterrestrial exploration targeting the Moon, Mars, the moons of Mars, and various asteroids. This contribution discusses how current state-of-the-art Earth-based testing for designing rovers and landers for these missions currently leads to overly optimistic conclusions about the behavior of these devices upon deployment on the targeted celestial bodies. The key misconception is that gravitational offset is necessary during the \textit{terramechanics} testing of rover and lander prototypes on Earth. The body of evidence supporting our argument is tied to a small number of studies conducted during parabolic flights and insights derived from newly revised scaling laws. We argue that what has prevented the community from fully diagnosing the problem at hand is the absence of effective physics-based models capable of simulating terramechanics under low gravity conditions. We developed such a physics-based simulator and utilized it to gauge the mobility of early prototypes of the Volatiles Investigating Polar Exploration Rover (VIPER), which is slated to depart for the Moon in November 2024. This contribution discusses the results generated by this simulator, how they correlate with physical test results from the NASA-Glenn SLOPE lab, and the fallacy of the gravitational offset in rover and lander testing. The simulator developed is open sourced and made publicly available for unfettered use; it can support principled studies that extend beyond trafficability analysis to provide insights into in-situ resource utilization activities, e.g., digging, bulldozing, and berming in low gravity.

Intraoperative ultrasound (iUS) imaging has the potential to improve surgical outcomes in brain surgery. However, its interpretation is challenging, even for expert neurosurgeons. In this work, we designed the first patient-specific framework that performs brain tumor segmentation in trackerless iUS. To disambiguate ultrasound imaging and adapt to the neurosurgeon's surgical objective, a patient-specific real-time network is trained using synthetic ultrasound data generated by simulating virtual iUS sweep acquisitions in pre-operative MR data. Extensive experiments performed in real ultrasound data demonstrate the effectiveness of the proposed approach, allowing for adapting to the surgeon's definition of surgical targets and outperforming non-patient-specific models, neurosurgeon experts, and high-end tracking systems. Our code is available at: \url{https://github.com/ReubenDo/MHVAE-Seg}., Comment: Early accept at MICCAI 2024 - code available at: https://github.com/ReubenDo/MHVAE-Seg

We present an exploration of the unpolarized isovector proton generalized parton distributions (GPDs) $H^{u-d}(x, \xi, t)$ and $E^{u-d}(x, \xi, t)$ in the pseudo-distribution formalism using distillation. Taking advantage of the large kinematic coverage made possible by this approach, we present results on the moments of GPDs up to the order $x^3$ -- including their skewness dependence -- at a pion mass $m_\pi = 358$ MeV and a lattice spacing $a = 0.094$ fm., Comment: added discussion on isovector D-term and stability of fits with respect to the range in non-local separation of the operator

We present the identification of the non-linear dynamics of a novel hovercraft design, employing end-to-end deep learning techniques. Our experimental setup consists of a hovercraft propelled by racing drone propellers mounted on a lightweight foam base, allowing it to float and be controlled freely on an air hockey table. We learn parametrized physics-inspired non-linear models directly from data trajectories, leveraging gradient-based optimization techniques prevalent in machine learning research. The chosen model structure allows us to control the position of the hovercraft precisely on the air hockey table. We then analyze the prediction performance and demonstrate the closed-loop control performance on the real system.

Comet 67P/Churyumov-Gerasimenko (67P) become observable for the first time in 2021 since the Rosetta rendezvous in 2014--16. Here, we present pre-perihelion polarimetric measurements of 67P from 2021 performed with the Very Large Telescope (VLT), as well as post-perihelion polarimetric measurements from 2015--16 obtained with the VLT and the William Herschel Telescope (WHT). This new data covers a phase angle range of ~4-50{\deg} and presents polarimetric measurements of unprecedentedly high S/N ratio. Complementing previous measurements, the polarimetric phase curve of 67P resembles that of other Jupiter family comets and high-polarisation, dusty comets. Comparing pre- and post-perihelion data sets, we find only a marginal difference between the polarimetric phase curves. In our imaging maps, we detect various linear structures produced by the dust in the inner coma of the comet. Despite this, we find a homogeneous spread of polarisation around the photocentre throughout the coma and tail, in contrast to previous studies. Finally, we explore the consequences of image misalignments on both polarimetric maps and aperture polarimetric measurements., Comment: 15 pages, 10 figures, 2 tables. Accepted for publication in MNRAS, May 13th 2024

There is a significant overlap between people who are supported by income-related social welfare benefits, often in precarious situations, and those who experience greater digital exclusion. We report on a study of claimants using the UK's Universal Credit online welfare benefit system designed as, and still, "digital by default". Through data collection involving remote interviews (n=11) and online surveys (n=66), we expose claimants' own lived experiences interacting with this system. The claimants explain how digital channels can contribute to an imbalance of power and agency, at a time when their own circumstances mean they have reduced abilities, resources and capacities, and where design choices can adversely affect people's utility to leverage help from their own wider socio-technical ecosystems. We contribute eight recommendations from these accounts to inform the future design and development of digital welfare benefit systems for this population, to reduce digital barriers and harms., Comment: 22 pages

Value-at-risk (VaR) and expected shortfall (ES) are two commonly utilized metrics for quantifying financial risk. In this study, we review the widely employed Generalized Autoregressive Conditional Heteroskedasticity (GARCH) models. These models are explored with diverse distributional assumptions on innovation, including parametric, non-parametric, and `semi-parametric' that incorporates a parametric tail distribution based on extreme value theory. Additionally, we introduce a non-parametric local linear quantile autoregression (LLQAR) with kernel weights depending on the distance between the current loss and past losses, and decreasing in the time lag. To evaluate the performance of different methods for VaR and ES estimation, we employ a multi-criteria approach. This involves mean squared error assessment using simulated data, backtesting on both simulated data and US stocks, and application of the ESBootstrap test. The LLQAR method, which does not necessarily require stationarity assumptions, seems to perform better for simulated non-stationary data as well as real-world data, for estimating VaR and ES., Comment: 54 pages, 21 figures

Let $S_n$ be the symmetric group on the set $[n]:=\{1,2,\ldots,n\}$. Given a permutation $\sigma=\sigma_1\sigma_2 \cdots \sigma_n \in S_n$, we say it has a descent at index $i$ if $\sigma_i>\sigma_{i+1}$. Let $\mathcal{D}(\sigma)$ be the set of all descents of $\sigma$ and define $\mathcal{D}(S;n)=\{\sigma\in S_n\, | \,\mathcal{D}(\sigma)=S\}$. We study the Hamming metric and $\ell_\infty$-metric on the sets $\mathcal{D}(S;n)$ for all possible nonempty $S\subset[n-1]$ to determine the maximum possible value that these metrics can achieve when restricted to these subsets., Comment: 10 pages, 2 tables

The United States national forest inventory (NFI) serves as the foundation for forest aboveground biomass (AGB) and carbon accounting across the nation. These data enable design-based estimates of forest carbon stocks and stock-changes at state and regional levels, but also serve as inputs to model-based approaches for characterizing forest carbon stocks and stock-changes at finer resolutions. Although NFI tree and plot-level data are often treated as truth in these models, they are in fact estimates based on regional species-group models known collectively as the Component Ratio Method (CRM). In late 2023 the Forest Inventory and Analysis (FIA) program introduced a new National Scale Volume and Biomass Estimators (NSVB) system to replace CRM nationwide and offer more precise and accurate representations of forest AGB and carbon. Given the prevalence of model-based AGB studies relying on FIA, there is concern about the transferability of methods from CRM to NSVB models, as well as the comparability of existing CRM AGB products (e.g. maps) to new and forthcoming NSVB AGB products. To begin addressing these concerns we compared previously published CRM AGB maps to new maps produced using identical methods with NSVB AGB reference data. Our results suggest that models relying on passive satellite imagery (e.g. Landsat) provide acceptable estimates of point-in-time NSVB AGB and carbon stocks, but fail to accurately quantify growth in mature closed-canopy forests. We highlight that existing estimates, models, and maps based on FIA reference data are no longer compatible with NSVB, and recommend new methods as well as updated models and maps for accommodating this step-change. Our collective ability to adopt NSVB in our modeling and mapping workflows will help us provide the most accurate spatial forest carbon data possible in order to better inform local management and decision making., Comment: Manuscript: 16 pages, 7 figures; Supplements: 3 pages, 2 figures; Submitted to: Remote Sensing of Environment

Local relaxation after a quench in 1-D quantum many-body systems is a well known and very active problem with rich phenomenology. Except for pathological cases, the local relaxation is accompanied by the local restoration of the symmetries broken by the initial state that are preserved by the unitary evolution. Recently, the entanglement asymmetry has been introduced as a probe to study the interplay between symmetry breaking and relaxation in an extended quantum system. In particular, using the asymmetry, it has been shown that the more a symmetry is initially broken, the faster it may be restored. This surprising effect, which has been also observed in trapped-ion experiments, can be seen as a quantum version of the Mpemba effect and is manifested by the crossing at a finite time of the entanglement asymmetry curves of two different initial symmetry breaking configurations. In this paper we show, how, by tuning the initial state, the symmetry dynamics in free fermionic systems can display much richer behaviour than seen previously. In particular, for certain classes of initial states, including ground states of free fermionic models with long-range couplings, the entanglement asymmetry can exhibit multiple crossings. This illustrates that the existence of the quantum Mpemba effect can only be inferred by examining the late time behaviour of the entanglement asymmetry., Comment: 25 pages, 5 figures

For $R>0$, we give a rigorous probabilistic construction on the cylinder $\mathbb{R} \times (\mathbb{R}/(2\pi R\mathbb{Z}))$ of the (massless) Sinh-Gordon model. In particular we define the $n$-point correlation functions of the model and show that these exhibit a scaling relation with respect to $R$. The construction, which relies on the massless Gaussian Free Field, is based on the spectral analysis of a quantum operator associated to the model. Using the theory of Gaussian multiplicative chaos, we prove that this operator has discrete spectrum and a strictly positive ground state., Comment: 36 pages, typos corrected

Non-adiabatic production of massive particles is a generic feature of many inflationary mechanisms. If sufficiently massive, these particles can leave features in the cosmic microwave background (CMB) that are not well-captured by traditional correlation function analyses. We consider a scenario in which particle production occurs only in a narrow time-interval during inflation, eventually leading to CMB hot- or coldspots with characteristic shapes and sizes. Searching for such features in CMB data is analogous to searching for late-Universe hot- or coldspots, such as those due to the thermal Sunyaev-Zel'dovich (tSZ) effect. Exploiting this data-analysis parallel, we perform a search for particle-production hotspots in the Planck PR4 temperature dataset, which we implement via a matched-filter analysis. Our pipeline is validated on synthetic observations and found to yield unbiased constraints on sufficiently large hotspots across $\approx 60\%$ of the sky. After removing point sources and tSZ clusters, we find no evidence for new physics and place novel bounds on the coupling between the inflaton and massive particles. These bounds are strongest for larger hotspots, produced early in inflation, whilst sensitivity to smaller hotspots is limited by noise and beam effects. Through such methods we can constrain particles with masses $\mathcal{O}(100)$ times larger than the inflationary Hubble scale, which represents possibly the highest energies ever directly probed with observational data., Comment: 12 pages, 7 figures, submitted to Phys. Rev. D

HCI and design researchers have designed, adopted, and customized a range of ethics-focused methods to inscribe values and support ethical decision making in a design process. In this work-in-progress, we add to this body of resources, constructing a method that surfaces the designer's intentions in an action-focused way, encouraging consideration of both manipulative and value-centered roles. Anti-Heroes is a card deck that allows a designer to playfully take on pairs of manipulative (Anti-Hero) and value-centered (Hero) roles during design ideation/conceptualization, evaluation, and ethical dialogue. The card deck includes twelve cards with Anti-Hero and Hero faces, along with three action cards that include reflective questions for different play modes. Alongside the creation of the Anti-Hero card deck, we describe the evaluation and iteration of the card deck through playtesting sessions with four groups of three design students. We propose implications of Anti-Heros for technology and design education and practice.

Thanks to its solid theoretical foundation, the SHAP framework is arguably one the most widely utilized frameworks for local explainability of ML models. Despite its popularity, its exact computation is known to be very challenging, proven to be NP-Hard in various configurations. Recent works have unveiled positive complexity results regarding the computation of the SHAP score for specific model families, encompassing decision trees, random forests, and some classes of boolean circuits. Yet, all these positive results hinge on the assumption of feature independence, often simplistic in real-world scenarios. In this article, we investigate the computational complexity of the SHAP score by relaxing this assumption and introducing a Markovian perspective. We show that, under the Markovian assumption, computing the SHAP score for the class of Weighted automata, Disjoint DNFs and Decision Trees can be performed in polynomial time, offering a first positive complexity result for the problem of SHAP score computation that transcends the limitations of the feature independence assumption., Comment: Accepted at ICML 2024

A new method for dark field imaging is introduced which uses scanned electron diffraction (or 4DSTEM - 4-dimensional scanning transmission electron microscopy) datasets as its input. Instead of working on simple summation of intensity, it works on a sparse representation of the diffraction patterns in terms of a list of their diffraction peaks. This is tested on a thin perovskite film containing structural ordering resulting in additional superlattice spots that reveal details of domain structures, and is shown to give much better selectivity and contrast than conventional virtual dark field imaging. It is also shown to work well in polycrystalline aggregates of CuO nanoparticles. In view of the higher contrast and selectivity, and the complete exclusion of diffuse scattering from the image formation, it is expected to be of significant benefit for characterisation of a wide variety of crystalline materials.

Patchy reionization generates kinematic Sunyaev-Zeldovich (kSZ) anisotropies in the cosmic microwave background (CMB). Large-scale velocity perturbations along the line of sight modulate the small-scale kSZ power spectrum, leading to a trispectrum (or four-point function) in the CMB that depends on the physics of reionization. We investigate the challenges in detecting this trispectrum and use tools developed for CMB lensing, such as realization-dependent bias subtraction and cross-correlation based estimators, to counter uncertainties in the instrumental noise and assumed CMB power spectrum. We also find that both lensing and extragalactic foregrounds can impart larger trispectrum contributions than the reionization kSZ signal. We present a range of mitigation methods for both of these sources of contamination, validated on microwave-sky simulations. We use ACT DR6 and Planck data to calculate an upper limit on the reionization kSZ trispectrum from a measurement dominated by foregrounds. The upper limit is about 50 times the signal predicted from recent simulations., Comment: Measurements and covariances will be made public upon publication

The quantum imaginary time evolution (QITE) methodology was developed to overcome a critical issue as regards non-unitarity in the implementation of imaginary time evolution on a quantum computer. QITE has since been used to approximate ground states of various physical systems. In this paper, we demonstrate a practical application of QITE as a quantum numerical solver for linear partial differential equations. Our algorithm takes inspiration from QITE in that the quantum state follows the same normalised trajectory in both algorithms. However, it is our QITE methodology's ability to track the scale of the state vector over time that allows our algorithm to solve differential equations. We demonstrate our methodology with numerical simulations and use it to solve the heat equation in one and two dimensions using six and ten qubits, respectively., Comment: 11 pages, 3 figures

There are many different forms of design knowledge that guide and shape a designer's ability to act and realize potential realities. Methods and schemas are examples of design knowledge commonly used by design researchers and designers alike. In this pictorial, we explore, engage, and describe the role of schemas as tools that can support design researchers in formulating methods to support design action, with our framing of method design specifically focused on ethical design complexity. We present four ways for method designers to engage with schema: 1) Systems to operationalize complex design constructs such as ethical design complexity through an A.E.I.O.YOU schema; 2) Classifiers to map existing methods and identify the possibility for new methods through descriptive semantic differentials; 3) Tools that enable the creation of methods that relate to one or more elements of the schema through creative departures from research to design; and 4) Interactive channels to playfully engage potential and new opportunities through schema interactivity.

A single-photon entangled state (or single-particle entangled state (SPES) in general) can offer a more secure way of encoding and processing quantum information than their multi-photon (or multi-particle) counterparts. The SPES generated via a 2D alternate quantum-walk setup from initially separable states can be either 3-way or 2-way entangled. This letter shows that the generated genuine three-way and nonlocal two-way SPES can be used as cryptographic keys to securely encode two distinct messages simultaneously. We detail the message encryption-decryption steps and show the resilience of the 3-way and 2-way SPES-based cryptographic protocols against eavesdropper attacks like intercept-and-resend and man-in-the-middle. We also detail the experimental realization of these protocols using a single photon, with the three degrees of freedom being OAM, path, and polarization. We have proved that the protocols have unconditional security for quantum communication tasks. The ability to simultaneously encode two distinct messages using the generated SPES showcases the versatility and efficiency of the proposed cryptographic protocol. This capability could significantly improve the throughput of quantum communication systems., Comment: 13 pages (including appendix), two figures and one table

A classical result of topos theory holds that the category of coalgebras for a Cartesian comonad on a topos is again a topos (Kock and Wraith, 1971). It is natural to refine this result to a topos-theoretic setting that includes universes. To this end, we introduce the notions of natural display topos and natural Cartesian display comonad, and show that the natural model of coalgebras for a natural Cartesian display comonad on a natural display topos is again a natural display topos. As an application, this result extends the approach to universes of Hofmann and Streicher (1997) from presheaf toposes to sheaf toposes with enough points. Whereas natural display toposes provide a categorical semantics for a form of extensional Martin-L\"of type theory, we also prove our main result in the more general setting of natural typoses, which encompasses models of intensional Martin-L\"of type theory. A natural Cartesian display comonad on a natural typos may also be used as a model for dependent type theory with an S4 box operator, or comonadic modality, as introduced by Nanevski et al. (2008). Modal contexts, which have been regarded as tricky to handle semantically, are interpreted as contexts of the natural typos of coalgebras. We sketch an interpretation within this approach. As part of the framework in which the above takes place, we introduce a refinement of the notion of natural model (see Awodey, 2018), which is (strictly 2-)equivalent to the notion of full, split comprehension category (see Jacobs, 1993), rather than the notion of category with attributes (Cartmell 1978)., Comment: PhD Thesis, Carnegie Mellon University

In this chapter, we address the challenge of exploring the posterior distributions of Bayesian inverse problems with computationally intensive forward models. We consider various multivariate proposal distributions, and compare them with single-site Metropolis updates. We show how fast, approximate models can be leveraged to improve the MCMC sampling efficiency., Comment: To appear in the Handbook of Markov Chain Monte Carlo (2nd edition)

Tracking objects in three-dimensional space is critical for autonomous driving. To ensure safety while driving, the tracker must be able to reliably track objects across frames and accurately estimate their states such as velocity and acceleration in the present. Existing works frequently focus on the association task while either neglecting the model performance on state estimation or deploying complex heuristics to predict the states. In this paper, we propose STT, a Stateful Tracking model built with Transformers, that can consistently track objects in the scenes while also predicting their states accurately. STT consumes rich appearance, geometry, and motion signals through long term history of detections and is jointly optimized for both data association and state estimation tasks. Since the standard tracking metrics like MOTA and MOTP do not capture the combined performance of the two tasks in the wider spectrum of object states, we extend them with new metrics called S-MOTA and MOTPS that address this limitation. STT achieves competitive real-time performance on the Waymo Open Dataset., Comment: ICRA 2024

Tensor processing infrastructures such as deep learning frameworks and specialized hardware accelerators have revolutionized how computationally intensive code from domains such as deep learning and image processing is executed and optimized. These infrastructures provide powerful and expressive abstractions while ensuring high performance. However, to utilize them, code must be written specifically using the APIs / ISAs of such software frameworks or hardware accelerators. Importantly, given the fast pace of innovation in these domains, code written today quickly becomes legacy as new frameworks and accelerators are developed, and migrating such legacy code manually is a considerable effort. To enable developers in leveraging such DSLs while preserving their current programming paradigm, we introduce Tenspiler, a verified lifting-based compiler that uses program synthesis to translate sequential programs written in general-purpose programming languages (e.g., C++ or Python code) into tensor operations. Central to Tenspiler is our carefully crafted yet simple intermediate language, named TensIR, that expresses tensor operations. TensIR enables efficient lifting, verification, and code generation. Currently, Tenspiler already supports \textbf{six} DSLs, spanning a broad spectrum of software and hardware environments. Furthermore, we show that new backends can be easily supported by Tenspiler by adding simple pattern-matching rules for TensIR. Using 10 real-world code benchmark suites, our experimental evaluation shows that by translating code to be executed on \textbf{6} different software frameworks and hardware devices, Tenspiler offers on average 105$\times$ kernel and 9.65$\times$ end-to-end execution time improvement over the fully-optimized sequential implementation of the same benchmarks.

The early dark energy (EDE) extension to $\Lambda$CDM has been proposed as a candidate scenario to resolve the "Hubble tension". We present new constraints on the EDE model by incorporating new data from the Dark Energy Spectroscopic Instrument (DESI) Baryon Acoustic Oscillation (BAO) survey and CMB lensing measurements from the Atacama Cosmology Telescope (ACT) DR6 and \textit{Planck} NPIPE data. We do not find evidence for EDE. The maximum fractional contribution of EDE to the total energy density is $f_\mathrm{EDE}< 0.091 \; (95\% \; \mathrm{CL} )$ from our baseline combination of \textit{Planck} CMB, CMB lensing, and DESI BAO. Our strongest constraints on EDE come from the combination of \textit{Planck} CMB and CMB lensing alone, yielding $f_\mathrm{EDE}< 0.070 \; (95\% \; \mathrm{CL} )$. We also explore extensions of $\Lambda$CDM beyond the EDE parameters by treating the total neutrino mass as a free parameter, finding $\sum m_\nu < 0.096 \,\, {\rm eV} \; (95\% \; \mathrm{CL} )$ and $f_\mathrm{EDE}< 0.087 \; (95\% \; \mathrm{CL} )$. For the first time in EDE analyses, we perform Bayesian parameter estimation using neural network emulators of cosmological observables, which are on the order of a hundred times faster than full Boltzmann solutions., Comment: 8+7 pages, 3+6 figures

We present JWST/NIRSpec high-resolution spectroscopy G395H/F290LP of MACS0647-JD, a gravitationally lensed galaxy merger at $z=10.167$. The new spectroscopy, which is acquired for the two lensed images (JD1 and JD2), detects and resolves emission lines in the rest-frame ultraviolet (UV) and blue optical, including the resolved [OII]3726,3729 doublet, [NeIII]3870, [HeI]3890, H$\delta$, H$\gamma$, and [OIII]4363. This is the first observation of the resolved [OII]3726,3729 doublet for a galaxy at $z>8$. We measure a line flux ratio [OII]3729/3726 $= 0.9 \pm 0.3$, which corresponds to an estimated electron density of $\log(n_{e} / \rm{cm}^{-3}) = 2.9 \pm 0.5$. This is significantly higher than the electron densities of local galaxies reported in the literature. We compile the measurements from the literature and further analyze the redshift evolution of $n_{e}$. We find that the redshift evolution follows the power-law form of $n_{e} = A\times (1+z)^{p}$ with $A=54^{+31}_{-23}$ cm$^{-3}$ and $p=1.2^{+0.4}_{-0.4}$. This power-law form may be explained by a combination of metallicity and morphological evolution of galaxies, which become, on average, more metal-poor and more compact with increasing redshift., Comment: 20 pages, 11 figures. Accepted for publication in ApJ

JWST spectroscopy has revolutionized our understanding of galaxies in the early universe. Covering wavelengths up to $5.3\,{\rm \mu m}$, NIRSpec can detect rest-frame optical emission lines H$\alpha$ out to $z = 7$ and [O III] to $z = 9.5$. Observing these lines in more distant galaxies requires longer wavelength spectroscopy with MIRI. Here we present MIRI MRS IFU observations of the lensed galaxy merger MACS0647$-$JD at $z = 10.165$. With exposure times of 4.2 hours in each of two bands, we detect H$\alpha$ at $9\sigma$, [O III]$\,\lambda5008$ at $11\sigma$, and [O III]$\,\lambda4960$ at $3\sigma$. Combined with previously reported NIRSpec spectroscopy that yields seven emission lines including the auroral line [O III]$\,\lambda4363$, we present the first direct metallicity measurement of a $z > 10$ galaxy: $12+{\rm log(O/H)}= 7.79\pm0.09$, or $0.13^{+0.02}_{-0.03}\,Z_{\odot}$. This is similar to galaxies at $z \sim 4 - 9$ with direct metallicity measurements, though higher than expected given the high specific star formation rate ${\rm log(sSFR / yr^{-1})} = -7.4 \pm 0.3$. We further constrain the ionization parameter ${\rm log}(U)$ = $-1.9 \pm 0.1$, ionizing photon production efficiency ${\rm log}(\xi_{\rm ion})$ = $25.3\pm0.1$, and star formation rate $5.0\pm0.6\,M_{\odot}/{\rm yr}$ within the past $10\,{\rm Myr}$. These observations demonstrate the combined power of JWST NIRSpec and MIRI for studying galaxies in the first $500$ million years., Comment: 16 pages, 8 figures, accepted for publication in ApJ

We present a multimode theory of squeezed state generation in resonant systems valid for arbitrary pump power and including pump depletion. The Hamiltonian is written in terms of asymptotic-in and -out fields from scattering theory, capable of describing a general interaction. As an example we consider the lossy generation of a highly squeezed state by an effective second-order interaction in a silicon nitride ring resonator point-coupled to a waveguide. We calculate the photon number, Schmidt number, and the second-order correlation function of the generated state in the waveguide. The treatment we present provides a path forward to study the deterministic generation of non-Gaussian states in resonant systems., Comment: 14 pages, 4 figures

In June 1976, a pristine meteorite stone weighing approximately 1 kg, fully covered with a fresh black fusion crust, was collected on a mountain road in the high-altitude Alpine environment. The recovery took place while clearing the remnants of a snow avalanche, 2 km northwest of Ischgl in Austria. Subsequent to its retrieval, the specimen remained in the finder's private residence without undergoing any scientific examination or identification until 2008, when it was brought to the University of Innsbruck. The sample was classified as a well-preserved LL6 chondrite, with a W0 weathering grade, implying a relatively short time between the meteorite fall and its retrieval. To investigate the potential connection between the Ischgl meteorite and a recorded fireball event, we have reviewed all documented fireballs ever photographed by German fireball camera stations. This examination led us to identify the fireball EN241170 observed in Germany by ten different European Network stations on the night of November 23/24, 1970, as the most likely candidate. We employed state-of-the-art techniques to reconstruct the fireball's trajectory, and to reproduce both its luminous and dark flight phases in detail. We find that the determined strewn field and the generated heat map closely align with the recovery location of the Ischgl meteorite. Furthermore, the measured radionuclide data reported here indicate that the pre-atmospheric size of the Ischgl meteoroid is consistent with the mass estimate inferred from our deceleration analysis along the trajectory. Our findings strongly support the conclusion that the Ischgl meteorite originated from the EN241170 fireball, effectively establishing it as a confirmed meteorite fall. This discovery enables to determine, along with the physical properties, also the heliocentric orbit and cosmic history of the Ischgl meteorite., Comment: Accepted for publication in Meteoritics & Planetary Science

Climate risk assessments must account for a wide range of possible futures, so scientists often use simulations made by numerous global climate models to explore potential changes in regional climates and their impacts. Some of the latest-generation models have high effective climate sensitivities or EffCS. It has been argued these so-called hot models are unrealistic and should therefore be excluded from analyses of climate change impacts. Whether this would improve regional impact assessments, or make them worse, is unclear. Here we show there is no universal relationship between EffCS and projected changes in a number of important climatic drivers of regional impacts. Analysing heavy rainfall events, meteorological drought, and fire weather in different regions, we find little or no significant correlation with EffCS for most regions and climatic drivers. Even when a correlation is found, internal variability and processes unrelated to EffCS have similar effects on projected changes in the climatic drivers as EffCS. Model selection based solely on EffCS appears to be unjustified and may neglect realistic impacts, leading to an underestimation of climate risks., Comment: Preprint, 30 pages, 4 figures and 2 tables

Metal halide perovskites have emerged as an exciting class of materials for applications in solar energy harvesting, optical devices, catalysis, and other photophysical applications. Many of the exciting properties of halide perovskites are tied to their soft, dynamic, and anharmonic lattice. In particular, the precise coupling between anharmonic lattice dynamics and electronic fluctuations is not completely understood. To build an understanding of this coupling, we use ab initio molecular dynamics simulations supplemented by the calculation of maximally localized Wannier functions to carry out a dynamic group theory analysis of local electron density fluctuations and how these fluctuations are coupled to lattice fluctuations in the model inorganic halide perovskite CsSnBr3. We detail symmetry-dependent couplings between vibrational modes, including octahedral tilting. Importantly, we suggest that the large anharmonicity of some of the vibrational modes in CsSnBr3 result from electron rotation--nuclear translation coupling, in analogy to rotation--translation coupling effects in molecular plastic crystals. We also identify electronic fluctuations in the Cs cation that couple to distortions in the surrounding Sn-Br cubic coordination environment. We anticipate that our approach and resulting insights into electronic fluctuations will aid in further understanding the role of the fluctuating lattice in determining important physical properties of halide perovskites and beyond., Comment: 15 pages, 15 figures

We present recent JWST NIRCam imaging observations of SPT0615-JD (also known as the Cosmic Gems Arc), lensed by the galaxy cluster SPT-CL J0615-5746. The 5-arcsec-long arc is the most highly magnified $z>10$ galaxy known, straddling the lensing critical curve and revealing five star clusters with radii $\sim 1$ pc or less. We measure the full arc to have F200W 24.5 AB mag, consisting of two mirror images, each 25.3 AB mag with a magnification $\mu \sim 60$ (delensed 29.7 AB mag, $M_{UV} = -17.8$). The galaxy has an extremely strong Lyman break F115W$-$F200W $>3.2$ mag ($2\sigma$ lower limit), is undetected in all bluer filters ($< 2\sigma$), and has a very blue continuum slope redward of the break ($\beta = -2.7 \pm 0.1$), resulting in a photometric redshift $z_{phot} = 10.2 \pm 0.2$ (95% confidence) with no significant likelihood below $z < 9.8$. Based on SED fitting to the total photometry, we estimate an intrinsic stellar mass of $M_{*} \sim 2.4 - 5.6 \times 10^{7} M_{\odot}$, young mass-weighted age of $\sim 21 - 79$ Myr, low dust content ($A_V < 0.15$), and a low metallicity of $\lesssim 1\%~Z_{\odot}$. We identify a fainter third counterimage candidate within 2.2 arcsec of the predicted position, lensed to AB mag 28.4 and magnified by $\mu \sim 2$, suggesting the fold arc may only show $\sim60$% of the galaxy. SPT0615-JD is a unique laboratory to study star clusters observed within a galaxy just 460 Myr after the Big Bang., Comment: 22 pages, 8 figures, 4 tables, submitted to ApJ

The attachment of redox active molecules to transition metal dichalcogenides (TMDs), such as MoS2, constitutes a promising approach for designing electrochemically switchable devices through the control of the material charge/spin transport properties by the redox state of the grafted molecule and thus the applied electrical potential. In this work, defective plasma treated MoS2 is functionalized by a ferrocene derivative and thoroughly investigated by various characterization techniques, such as Raman, photoluminescence, X-ray photoelectron spectroscopies, atomic force microscopy (AFM) and electrochemistry. Furthermore, in-plane and out-of-plane conductive-AFM measurements (I-V and first derivative dI/dV-V curves) are measured to investigate the effect of the chemical functionalization of MoS2 on the electron transport properties. While the conduction and valence bands are determined at +0.7 and -1.2 eV with respect of the electrode Fermi energy for pristineMoS2, additional states in an energy range of ca. 0.45 eV below the MoS2 conduction band are measured after plasma treatment, attributed to S-vacancies. For ferrocene functionalized MoS2, the S-vacancy states are no longer observed resulting from the defect healing. However, two bumps at lower voltages in the dI/dV-V indicate a contribution to the electron transport through ferrocene HOMO, which is located in the MoS2 band gap at ca. 0.4-0.6 eV below the Fermi energy. These results are in good agreement with theoretical density functional theory (DFT) calculations and UV photoelectron spectroscopy (UPS) measurements., Comment: Full manuscript with figures and supporting information

The demand-side control of district heating networks is notoriously challenging due to the large number of connected users and the high number of states to be considered. To overcome these challenges, this paper presents a hierarchical optimization scheme using the flexibility in heating demand provided by the users to improve the performance of the network. This hierarchical scheme relies on a low level controller to calculate the costs for a subsystem over a given set of potential pressure drops for that subsystem. The high level controller then uses these calculated costs to determine the optimal set of pressure drops for every subgraph of the partitioned network. The proposed hierarchical optimization scheme is demonstrated on a representative 20 user district heating network, resulting in a 67\% reduction in bypass mass flow while ensuring all network users stay within 2 \degree C of their desired nominal temperatures., Comment: Accepted for publication by European Control Conference 2024