Search

We present a fully implicit formulation of coupled fluid flow and geomechanics for fluid injection/withdrawal in fractured reservoirs in the context of CO2storage. Utilizing a Galerkin finite-element approach, both flow and poroelasticity equations are discretized on a shared three-dimensional mesh. The fluid flow is assumed to be single-phase. The hydraulic behaviour of fractures is represented through a double-nodes flow element, which allows to efficiently model longitudinal and transversal fracture permeabilities. In addressing the mechanical subproblem, fractures are explicitly modelled using cohesive elements to account for contact, friction and opening phenomena. The nonlinear set of equations is solved implicitly through an iterative partitioned conjugate gradient procedure, extending its traditional application to continuous problems to those involving explicit discontinuities such as faults and fractures. The model's accuracy is verified against analytical solutions for different geomechanical problems, notably for the growth of a frictional slip rupture along a fault due to fluid injection. Such a particularly challenging benchmark for a critically stressed fault is here reproduced for the first time by a finite-element based scheme. The capabilities of the developed parallel solver are then illustrated through a scenario involving CO2 injection into a faulted aquifer. The original solver code, tutorials, and data visualization routines are publicly accessible.

The synthesis of binary nanocrystal superlattices (BNSLs) enables the targeted integration of orthogonal physical properties, like photoluminescence and magnetism, into a single superstructure, unlocking a vast design space for multifunctional materials. Yet, the formation mechanism of BNSLs remains poorly understood, restricting the use of simulation to predict the structure and properties of the superlattices. Here, we use a combination of in situ scattering and molecular simulation to elucidate the self-assembly of two common BNSLs through emulsion templating. Our self-assembly experiments reveal that no intermediate structures precede the formation of the final binary phases, indicating that their formation proceeds through classical nucleation. Using simulations, we find that, despite the formation of AlB2 and NaZn13 typically being attributed to entropy, their self-assembly is most consistent with the nanocrystals possessing short-range interparticle attraction, which we find can dramatically accelerate nucleation kinetics in BNSLs. We also find homogenous, classical nucleation in simulation, corroborating our experiments. These results establish a robust correspondence between experiment and theory, paving the way towards a priori prediction of BNSLs.

Plug&Play (PnP) diffusion models are state-of-the-art methods in computed tomography (CT) reconstruction. Such methods usually consider applications where the sinogram contains a sufficient amount of information for the posterior distribution to be peaked, and consequently are evaluated using image-to-image metrics such as PSNR/SSIM. Instead, we are interested in reconstructing compressible flow images from sinograms having a small number of projections, which results in a posterior distribution no longer peaked or even multimodal. Thus, in this paper, we aim at evaluating the approximate posterior of PnP diffusion models and introduce two posterior evaluation criteria. We quantitatively evaluate three PnP diffusion methods on three different datasets for several numbers of projections. We surprisingly find that, for each method, the approximate posterior deviates from the true posterior when the number of projections decreases.

Typical legged locomotion controllers are designed or trained offline. This is in contrast to many animals, which are able to locomote at birth, and rapidly improve their locomotion skills with few real-world interactions. Such motor control is possible through oscillatory neural networks located in the spinal cord of vertebrates, known as Central Pattern Generators (CPGs). Models of the CPG have been widely used to generate locomotion skills in robotics, but can require extensive hand-tuning or offline optimization of inter-connected parameters with genetic algorithms. In this paper, we present a framework for the \textit{online} optimization of the CPG parameters through Bayesian Optimization. We show that our framework can rapidly optimize and adapt to varying velocity commands and changes in the terrain, for example to varying coefficients of friction, terrain slope angles, and added mass payloads placed on the robot. We study the effects of sensory feedback on the CPG, and find that both force feedback in the phase equations, as well as posture control (Virtual Model Control) are both beneficial for robot stability and energy efficiency. In hardware experiments on the Unitree Go1, we show rapid optimization (in under 3 minutes) and adaptation of energy-efficient gaits to varying target velocities in a variety of scenarios: varying coefficients of friction, added payloads up to 15 kg, and variable slopes up to 10 degrees. See demo at: https://youtu.be/4qq5leCI2AI, Comment: Accepted by IROS2024

Quantum computing presents a promising alternative for the direct simulation of quantum systems with the potential to explore chemical problems beyond the capabilities of classical methods. However, current quantum algorithms are constrained by hardware limitations and the increased number of measurements required to achieve chemical accuracy. To address the measurement challenge, techniques for grouping commuting and anti-commuting terms, driven by heuristics, have been developed to reduce the number of measurements needed in quantum algorithms on near-term quantum devices. In this work, we propose a probabilistic framework using GFlowNets to group fully (FC) or qubit-wise commuting (QWC) terms within a given Hamiltonian. The significance of this approach is demonstrated by the reduced number of measurements for the found groupings; 51% and 67% reduction factors respectively for FC and QWC partitionings with respect to greedy coloring algorithms, highlighting the potential of GFlowNets for future applications in the measurement problem. Furthermore, the flexibility of our algorithm extends its applicability to other resource optimization problems in Hamiltonian simulation, such as circuit design., Comment: 8 pages, 2 figures. Accepted for Machine Learning and the Physical Sciences Workshop, NeurIPS 2024. Submission Number: 167

Compositionality is believed to be fundamental to intelligence. In humans, it underlies the structure of thought, language, and higher-level reasoning. In AI, compositional representations can enable a powerful form of out-of-distribution generalization, in which a model systematically adapts to novel combinations of known concepts. However, while we have strong intuitions about what compositionality is, there currently exists no formal definition for it that is measurable and mathematical. Here, we propose such a definition, which we call representational compositionality, that accounts for and extends our intuitions about compositionality. The definition is conceptually simple, quantitative, grounded in algorithmic information theory, and applicable to any representation. Intuitively, representational compositionality states that a compositional representation satisfies three properties. First, it must be expressive. Second, it must be possible to re-describe the representation as a function of discrete symbolic sequences with re-combinable parts, analogous to sentences in natural language. Third, the function that relates these symbolic sequences to the representation, analogous to semantics in natural language, must be simple. Through experiments on both synthetic and real world data, we validate our definition of compositionality and show how it unifies disparate intuitions from across the literature in both AI and cognitive science. We also show that representational compositionality, while theoretically intractable, can be readily estimated using standard deep learning tools. Our definition has the potential to inspire the design of novel, theoretically-driven models that better capture the mechanisms of compositional thought.

Understanding the variation of the optimal value with respect to change in the data is an old problem of mathematical optimisation. This paper focuses on the linear problem $f(\lambda) = \min c^t x$ such that $(A+\lambda D)x \leq b$, where $\lambda$ is an unknown parameter that varies within an interval and $D$ is a matrix modifying the coefficients of the constraint matrix $A$. This problem is used to analyse the impact of multiple affine changes in the constraint matrix on the objective function. The function $f(\lambda)$ can have an erratic behaviour and computing it for a large number of points is computationally heavy while not providing any guarantees in between the computed points. As a new approach to the problem, we derive several bounding methods that provide guarantees on the objective function's behaviour. Those guarantees can be exploited to avoid recomputing the problem for numerous $\lambda$. The bounding methods are based on approaches from robust optimisation or Lagrangian relaxations. For each approach, we derive three methods of increasing complexity and precision, one that provides constant bounds, one that provides $\lambda$-dependant bounds and envelope bounds. We assess each bounding method in terms of precision, availability and timing. We show that for a large number of problems, the bound approach outperforms the naive sampling approach considered with 100 points while still providing a good precision and stronger guarantees on a large dataset of problems. We also introduce an iterative algorithm that uses these bounds to compute an approximation of the original function within a given error threshold and discuss its performances.

Topological physics has broadened its scope from the study of topological insulating phases to include nodal phases containing band structure singularities. The geometry of the corresponding quantum states is described by the quantum metric which provides a theoretical framework for explaining phenomena that conventional approaches fail to address. The field has become even broader by encompassing non-Hermitian singularities: in addition to Dirac, Weyl nodes, or nodal lines, it is now common to encounter exceptional points, exceptional or Weyl rings, and even Weyl spheres. They give access to fascinating effects that cannot be reached within the Hermitian picture. However, the quantum geometry of non-Hermitian singularities is not a straightforward extension of the Hermitian one, remaining far less understood. Here, we study experimentally and theoretically the dynamics of wave packets at exceptional rings stemming from Dirac points in a photonic honeycomb lattice. First, we demonstrate a transition between conical diffraction and non-Hermitian broadening in real space. Next, we predict and demonstrate a new non-Hermitian effect in the reciprocal space, induced by the non-orthogonality of the eigenstates. We call it transverse non-Hermitian drift, and its description requires biorthogonal quantum metric. The non-Hermitian drift can be used for applications in beam steering.

The goal of machine learning is generalization. While the No Free Lunch Theorem states that we cannot obtain theoretical guarantees for generalization without further assumptions, in practice we observe that simple models which explain the training data generalize best: a principle called Occam's razor. Despite the need for simple models, most current approaches in machine learning only minimize the training error, and at best indirectly promote simplicity through regularization or architecture design. Here, we draw a connection between Occam's razor and in-context learning: an emergent ability of certain sequence models like Transformers to learn at inference time from past observations in a sequence. In particular, we show that the next-token prediction loss used to train in-context learners is directly equivalent to a data compression technique called prequential coding, and that minimizing this loss amounts to jointly minimizing both the training error and the complexity of the model that was implicitly learned from context. Our theory and the empirical experiments we use to support it not only provide a normative account of in-context learning, but also elucidate the shortcomings of current in-context learning methods, suggesting ways in which they can be improved. We make our code available at https://github.com/3rdCore/PrequentialCode.

Photonics provides a viable path to a scalable fault-tolerant quantum computer. The natural framework for this platform is measurement-based quantum computation, where fault-tolerant graph states supersede traditional quantum error-correcting codes. However, the existing formalism for foliation - the construction of fault-tolerant graph states - does not reveal how certain properties, such as single-shot error correction, manifest in the measurement-based setting. We introduce the fault complex, a representation of dynamic quantum error correction protocols particularly well-suited to describe foliation. Our approach enables precise computation of fault tolerance properties of foliated codes and provides insights into circuit-based quantum computation. Analyzing the fault complex leads to improved thresholds for three- and four-dimensional toric codes, a generalization of stability experiments, and the existence of single-shot lattice surgery with higher-dimensional topological codes., Comment: 13 pages, 5 figures

Rapid growth of high-dimensional datasets in fields such as single-cell RNA sequencing and spatial genomics has led to unprecedented opportunities for scientific discovery, but it also presents unique computational and statistical challenges. Traditional methods struggle with geometry-aware data generation, interpolation along meaningful trajectories, and transporting populations via feasible paths. To address these issues, we introduce Geometry-Aware Generative Autoencoder (GAGA), a novel framework that combines extensible manifold learning with generative modeling. GAGA constructs a neural network embedding space that respects the intrinsic geometries discovered by manifold learning and learns a novel warped Riemannian metric on the data space. This warped metric is derived from both the points on the data manifold and negative samples off the manifold, allowing it to characterize a meaningful geometry across the entire latent space. Using this metric, GAGA can uniformly sample points on the manifold, generate points along geodesics, and interpolate between populations across the learned manifold using geodesic-guided flows. GAGA shows competitive performance in simulated and real-world datasets, including a 30% improvement over the state-of-the-art methods in single-cell population-level trajectory inference.

Semiconductor membranes find their widespread use in various research fields targeting medical, biological, environmental, and optical applications. Often such membranes derive their functionality from an inherent nanopatterning, which renders the determination of their, e.g., optical, electronic, mechanical, and thermal properties a challenging task. In this work we demonstrate the non-invasive, all-optical thermal characterization of around 800-nm-thick and 150-$\mu$m-wide membranes that consist of wurtzite GaN and a stack of In$_{0.15}$Ga$_{0.85}$N quantum wells as a built-in light source. Due to their application in photonics such membranes are bright light emitters, which challenges their non-invasive thermal characterization by only optical means. As a solution, we combine two-laser Raman thermometry with (time-resolved) photoluminescence measurements to extract the in-plane (i.e., $c$-plane) thermal conductivity $\kappa_{\text{in-plane}}$ of our membranes. Based on this approach, we can disentangle the entire laser-induced power balance during our thermal analysis, meaning that all fractions of reflected, scattered, transmitted, and reemitted light are considered. As a result of our thermal imaging via Raman spectroscopy, we obtain $\kappa_{\text{in-plane}}\,=\,165^{+16}_{-14}\,$Wm$^{-1}$K$^{-1}$ for our best membrane, which compares well to our simulations yielding $\kappa_{\text{in-plane}}\,=\,177\,$Wm$^{-1}$K$^{-1}$ based on an ab initio solution of the linearized phonon Boltzmann transport equation. Our work presents a promising pathway towards thermal imaging at cryogenic temperatures, e.g., when aiming to elucidate experimentally different phonon transport regimes via the recording of non-Fourier temperature distributions., Comment: Main text (4 figures and 15 pages) and Supplemental Material (3 supplemental figures and 4 pages)

Homomorphic encryption is one of the main solutions for building secure and privacy-preserving solutions for Machine Learning as a Service. This motivates the development of homomorphic algorithms for the main building blocks of AI, typically for the components of the various types of neural networks architectures. Among those components, we focus on the Softmax function, defined by $\mathrm{SM}(\mathbf{x}) = \left(\exp(x_i) / \sum_{j=1}^n \exp(x_j) \right)_{1\le i\le n}$. This function is deemed to be one of the most difficult to evaluate homomorphically, because of its multivariate nature and of the very large range of values for $\exp(x_i)$. The available homomorphic algorithms remain restricted, especially in large dimensions, while important applications such as Large Language Models (LLM) require computing Softmax over large dimensional vectors. In terms of multiplicative depth of the computation (a suitable measure of cost for homomorphic algorithms), our algorithm achieves $O(\log n)$ complexity for a fixed range of inputs, where $n$ is the Softmax dimension. Our algorithm is especially adapted to the situation where we must compute many Softmax at the same time, for instance, in the LLM situation. In that case, assuming that all Softmax calls are packed into $m$ ciphtertexts, the asymptotic amortized multiplicative depth cost per ciphertext is, again over a fixed range, $O(1 + m/N)$ for $N$ the homomorphic ring degree. The main ingredient of our algorithms is a normalize-and-square strategy, which interlaces the exponential computation over a large range and normalization, decomposing both in stabler and cheaper smaller steps. Comparing ourselves to the state of the art, our experiments show, in practice, a good accuracy and a gain of a factor 2.5 to 8 compared to state of the art solutions., Comment: ACM Conference on Computer and Communications Security (CCS) 2024

LiDAR odometry is essential for many robotics applications, including 3D mapping, navigation, and simultaneous localization and mapping. LiDAR odometry systems are usually based on some form of point cloud registration to compute the ego-motion of a mobile robot. Yet, few of today's LiDAR odometry systems consider the domain-specific knowledge and the kinematic model of the mobile platform during the point cloud alignment. In this paper, we present Kinematic-ICP, a LiDAR odometry system that focuses on wheeled mobile robots equipped with a 3D LiDAR and moving on a planar surface, which is a common assumption for warehouses, offices, hospitals, etc. Our approach introduces kinematic constraints within the optimization of a traditional point-to-point iterative closest point scheme. In this way, the resulting motion follows the kinematic constraints of the platform, effectively exploiting the robot's wheel odometry and the 3D LiDAR observations. We dynamically adjust the influence of LiDAR measurements and wheel odometry in our optimization scheme, allowing the system to handle degenerate scenarios such as feature-poor corridors. We evaluate our approach on robots operating in large-scale warehouse environments, but also outdoors. The experiments show that our approach achieves top performances and is more accurate than wheel odometry and common LiDAR odometry systems. Kinematic-ICP has been recently deployed in the Dexory fleet of robots operating in warehouses worldwide at their customers' sites, showing that our method can run in the real world alongside a complete navigation stack., Comment: Submitted to ICRA 2025

Computing the diameter, and more generally, all eccentricities of an undirected graph is an important problem in algorithmic graph theory and the challenge is to identify graph classes for which their computation can be achieved in subquadratic time. Using a new recursive scheme based on the structural properties of median graphs, we provide a quasilinear-time algorithm to determine all eccentricities for this well-known family of graphs. Our recursive technique manages specifically balanced and unbalanced parts of the $\Theta$-class decomposition of median graphs. The exact running time of our algorithm is O(n log^4 n). This outcome not only answers a question asked by B{\'e}n{\'e}teau et al. (2020) but also greatly improves a recent result which presents a combinatorial algorithm running in time O(n^1.6408 log^{O(1)} n) for the same problem.Furthermore we also propose a distance oracle for median graphs with both polylogarithmic size and query time. Speaking formally, we provide a combinatorial algorithm which computes for any median graph G, in quasilinear time O(n log^4 n), vertex-labels of size O(log^3 n) such that any distance of G can be retrieved in time O(log^4 n) thanks to these labels.

The Large Magellanic Cloud (LMC) is a Milky Way (MW) satellite that is massive enough to gravitationally attract the MW disc and inner halo, causing significant motion of the inner MW with respect to the outer halo. In this work, we probe this interaction by constructing a sample of 9,866 blue horizontal branch (BHB) stars with radial velocities from the DESI spectroscopic survey out to 120 kpc from the Galactic centre. This is the largest spectroscopic set of BHB stars in the literature to date, and it contains four times more stars with Galactocentric distances beyond 50 kpc than previous BHB catalogues. Using the DESI BHB sample combined with SDSS BHBs, we measure the bulk radial velocity of stars in the outer halo and observe that the velocity in the Southern Galactic hemisphere is different by 3.7$\sigma$ from the North. Modelling the projected velocity field shows that its dipole component is directed at a point 22 degrees away from the LMC along its orbit, which we interpret as the travel direction of the inner MW. The velocity field includes a monopole term that is -24 km/s, which we refer to as compression velocity. This velocity is significantly larger than predicted by the current models of the MW and LMC interaction. This work uses DESI data from its first two years of observations, but we expect that with upcoming DESI data releases, the sample of BHB stars will increase and our ability to measure the MW-LMC interaction will improve significantly., Comment: 22 pages, 19 figures. Submitted to MNRAS

Many quantum chemical similarity measures have been derived and substantiated by applying concepts and quantities from information theory to the electron density. To justify the use of information theory, the electron density is usually equated to a probability distribution, despite the fact that such an assumption can give rise to inconsistencies such as negative Kullback-Leibler divergences. In this work we show, using quantum information theory, that both interpreting electron densities as probability distributions as well as any pragmatic normalization thereof gives rise to far-reaching information theoretical inconsistencies. By applying the theory of open quantum subsystems to quantum chemical partitions, we show using thought experiments that many statements concerning the information content of atoms in molecules need to be reconsidered. This study represents a step towards generalizing and reformulating the information theoretical approach to conceptual quantum chemistry beyond its classical perspective., Comment: 37 pages

Initially designed to detect and characterize exoplanets, extreme adaptive optics systems (AO) open a new window on the solar system by resolving its small bodies. Nonetheless, despite the always increasing performances of AO systems, the correction is not perfect, degrading their image and producing a bright halo that can hide faint and close moons. Using a reference point spread function (PSF) is not always sufficient due to the random nature of the turbulence. In this work, we present our method to overcome this limitation. It blindly reconstructs the AO-PSF directly in the data of interest, without any prior on the instrument nor the asteroid's shape. This is done by first estimating the PSF core parameters under the assumption of a sharp-edge and flat object, allowing the image of the main body to be deconvolved. Then, the PSF faint extensions are reconstructed with a robust penalization optimization, discarding outliers on-the-fly such as cosmic rays, defective pixels and moons. This allows to properly model and remove the asteroid's halo. Finally, moons can be detected in the residuals, using the reconstructed PSF and the knowledge of the outliers learned with the robust method. We show that our method can be easily applied to different instruments (VLT/SPHERE, Keck/NIRC2), efficiently retrieving the features of AO-PSFs. Compared with state-of-the-art moon enhancement algorithms, moon signal is greatly improved and our robust detection method manages to discriminate faint moons from outliers., Comment: arXiv admin note: text overlap with arXiv:2407.21548

This paper presents a framework for dynamic object catching using a quadruped robot's front legs while it stands on its rear legs. The system integrates computer vision, trajectory prediction, and leg control to enable the quadruped to visually detect, track, and successfully catch a thrown object using an onboard camera. Leveraging a fine-tuned YOLOv8 model for object detection and a regression-based trajectory prediction module, the quadruped adapts its front leg positions iteratively to anticipate and intercept the object. The catching maneuver involves identifying the optimal catching position, controlling the front legs with Cartesian PD control, and closing the legs together at the right moment. We propose and validate three different methods for selecting the optimal catching position: 1) intersecting the predicted trajectory with a vertical plane, 2) selecting the point on the predicted trajectory with the minimal distance to the center of the robot's legs in their nominal position, and 3) selecting the point on the predicted trajectory with the highest likelihood on a Gaussian Mixture Model (GMM) modelling the robot's reachable space. Experimental results demonstrate robust catching capabilities across various scenarios, with the GMM method achieving the best performance, leading to an 80% catching success rate. A video demonstration of the system in action can be found at https://youtu.be/sm7RdxRfIYg ., Comment: Accepted to IROS 2024

Performances of an adaptive optics (AO) system are directly linked with the quality of its alignment. During the instrument calibration, having open loop fast tools with a large capture range are necessary to quickly assess the system misalignment and to drive it towards a state allowing to close the AO loop. During operation, complex systems are prone to misalignments (mechanical flexions, rotation of optical elements, etc.) that potentially degrade the AO performances, creating a need for a monitoring tool to tackle their driftage. In this work, we first present an improved perturbative method to quickly assess large lateral errors in open loop. It uses the spatial correlation of the measured interaction matrix of a limited number of 2D spatial modes with a synthetic model. Then, we introduce a novel solution to finely measure and correct these lateral errors via the closed loop telemetry. Non-perturbative, this method consequently does not impact the science output of the instrument. It is based on the temporal correlation of 2D spatial frequencies in the deformable mirror commands. It is model-free (no need of an interaction matrix model) and sparse in the Fourier space, making it fast and easily scalable to complex systems such as future extremely large telescopes. Finally, we present some results obtained on the development bench of the GRAVITY+ extreme AO system (Cartesian grid, 1432 actuators). In addition, we show with on-sky results gathered with CHARA and GRAVITY/CIAO that the method is adaptable to non-conventional AO geometries (hexagonal grids, 60 actuators).

We introduce Pixtral-12B, a 12--billion-parameter multimodal language model. Pixtral-12B is trained to understand both natural images and documents, achieving leading performance on various multimodal benchmarks, surpassing a number of larger models. Unlike many open-source models, Pixtral is also a cutting-edge text model for its size, and does not compromise on natural language performance to excel in multimodal tasks. Pixtral uses a new vision encoder trained from scratch, which allows it to ingest images at their natural resolution and aspect ratio. This gives users flexibility on the number of tokens used to process an image. Pixtral is also able to process any number of images in its long context window of 128K tokens. Pixtral 12B substanially outperforms other open models of similar sizes (Llama-3.2 11B \& Qwen-2-VL 7B). It also outperforms much larger open models like Llama-3.2 90B while being 7x smaller. We further contribute an open-source benchmark, MM-MT-Bench, for evaluating vision-language models in practical scenarios, and provide detailed analysis and code for standardized evaluation protocols for multimodal LLMs. Pixtral-12B is released under Apache 2.0 license.

We propose a variant of the Shapley value, the group Shapley value, to interpret counterfactual simulations in structural economic models by quantifying the importance of different components. Our framework compares two sets of parameters, partitioned into multiple groups, and applying group Shapley value decomposition yields unique additive contributions to the changes between these sets. The relative contributions sum to one, enabling us to generate an importance table that is as easily interpretable as a regression table. The group Shapley value can be characterized as the solution to a constrained weighted least squares problem. Using this property, we develop robust decomposition methods to address scenarios where inputs for the group Shapley value are missing. We first apply our methodology to a simple Roy model and then illustrate its usefulness by revisiting two published papers.

Standing acoustic waves have been known to generate Eulerian time-mean 'streaming' flows at least since the seminal investigation of Lord Rayleigh in the 1880s. Nevertheless, a recent body of numerical and experimental evidence has shown that inhomogeneities in the ambient density distribution lead to much faster flows than arise in classical Rayleigh streaming. The emergence of these unusually strong flows creates new opportunities to enhance heat transfer in systems in which convective cooling cannot otherwise be easily achieved. To assess this possibility, a theoretical study of acoustic streaming in an ideal gas confined in a rectangular channel with top and bottom walls maintained at fixed but differing temperatures is performed. A two time scale system of equations is utilized to efficiently capture the coupling between the fast acoustic waves and the slowly evolving streaming flow, enabling strongly nonlinear regimes to be accessed. A large suite of numerical simulations is carried out to probe the streaming dynamics, to highlight the critical role played by baroclinically generated wave vorticity and to quantify the additional heat flux induced by the standing acoustic wave. Proper treatment of the two-way coupling between the waves and mean flow is found to be essential for convergence to a self-consistent steady state, and the variation of the resulting acoustically enhanced steady-state heat flux with both the amplitude of the acoustic wave and the O(1) aspect ratio of the channel is documented. For certain parameters, heat fluxes almost two orders of magnitude larger than those realizable by conduction alone can be attained.

In this technical report, we present Falcon Mamba 7B, a new base large language model based on the novel Mamba architecture. Falcon Mamba 7B is trained on 5.8 trillion tokens with carefully selected data mixtures. As a pure Mamba-based model, Falcon Mamba 7B surpasses leading open-weight models based on Transformers, such as Mistral 7B, Llama3.1 8B, and Falcon2 11B. It is on par with Gemma 7B and outperforms models with different architecture designs, such as RecurrentGemma 9B and RWKV-v6 Finch 7B/14B. Currently, Falcon Mamba 7B is the best-performing Mamba model in the literature at this scale, surpassing both existing Mamba and hybrid Mamba-Transformer models, according to the Open LLM Leaderboard. Due to its architecture, Falcon Mamba 7B is significantly faster at inference and requires substantially less memory for long sequence generation. Despite recent studies suggesting that hybrid Mamba-Transformer models outperform pure architecture designs, we demonstrate that even the pure Mamba design can achieve similar, or even superior results compared to the Transformer and hybrid designs. We make the weights of our implementation of Falcon Mamba 7B publicly available on https://huggingface.co/tiiuae/falcon-mamba-7b, under a permissive license.

The new NHSTA regulations require Automatic Emergency Braking (AEB) systems to operate at night, to protect pedestrians in the deadliest conditions. We propose thermal imaging as a new sensor to complement the AEB sensor suite, alongside the visible front camera and RADAR. This paper explores the benefits of visible-thermal fusion, proposes a windshield integration for such a system, and evaluates the minimum performance requirements for a thermal camera compliant with the NHSTA standards, based on a field study of pedestrian detection range., Comment: SIA VISION 2024, Oct 2024, Paris, France

Solution-based soft matter self-assembly (SA) promises unique materials properties from approaches including additive manufacturing/three-dimensional (3D) printing. We report direct ink writing derived, hierarchically porous transition metal nitride superconductors (SCs) and precursor oxides, structure-directed by Pluronics-family block copolymer (BCP) SA and heat treated in various environments. SCs with periodic lattices on three length scales show record nanoconfinement-induced upper critical field enhancements correlated with BCP molar mass, attaining values of 50 T for NbN and 8.1 T for non-optimized TiN samples, the first mapping of a tailorable SC property onto a macromolecular parameter. They reach surface areas above 120 m$^2$/g, the highest reported for compound SCs to date. Embedded printing enables the first BCP directed mesoporous non-self-supporting helical SCs. Results suggest that additive manufacturing may open pathways to mesoporous SCs with not only a variety of macroscopic form factors but enhanced properties from intrinsic, SA-derived mesostructures with substantial academic and technological promise.

Groundwater management, especially in regions like Tunisia, is challenging due to diverse stakeholder interests and the dry structure of climate, which is extremely challenging for the sustainability of water resources. This paper proposes an innovative approach to policy design by merging Problem Structuring Methods (PSMs) and the Policy-Knowledge, Concepts, Proposals (P-KCP) methodology. Utilizing cognitive maps and value trees, the study aims to generate new collective groundwater management practices. Bridging decision theory and design theory, the study addresses the gap in new alternative generation and highlights the P-KCP's role in innovative policy design. Integrating PSMs and C-K theory, the framework expands policy alternatives and advocates for participatory approaches. It emphasizes adaptability across contexts, provides replicable process descriptions, and encourages the creation of unconventional policy solutions. Ultimately, this comprehensive framework offers a practical guide for policy innovation and collaboration.

We investigate the Lyman-$\alpha$ (Ly$\alpha$) and Lyman continuum (LyC) properties of the Sunburst Arc, a $z=2.37$ gravitationally lensed galaxy with a multiply-imaged, compact region leaking LyC and a triple-peaked Ly$\alpha$ profile indicating direct Ly$\alpha$ escape. Non-LyC-leaking regions show a redshifted Ly$\alpha$ peak, a redshifted and central Ly$\alpha$ peak, or a triple-peaked Ly$\alpha$ profile. We measure the properties of the Ly$\alpha$ profile from different regions of the galaxy using $R\sim5000$ Magellan/MagE spectra. We compare the Ly$\alpha$ spectral properties to LyC and narrowband Ly$\alpha$ maps from Hubble Space Telescope (HST) imaging to explore the subgalactic Ly$\alpha-$LyC connection. We find strong correlations (Pearson correlation coefficient $r>0.6$) between the LyC escape fraction ($f_{\rm esc}^{\rm LyC}$) and Ly$\alpha$ (1) peak separation $v_{\rm{sep}}$, (2) ratio of the minimum flux density between the redshifted and blueshifted Ly$\alpha$ peaks to continuum flux density $f_{\rm{min}}/f_{\rm{cont}}$, and (3) equivalent width. We favor a complex \ion{H}{1} geometry to explain the Ly$\alpha$ profiles from non-LyC-leaking regions and suggest two \ion{H}{1} geometries that could diffuse and/or rescatter the central Ly$\alpha$ peak from the LyC-leaking region into our sightline across transverse distances of several hundred parsecs. Our results emphasize the complexity of Ly$\alpha$ radiative transfer and its sensitivity to the anisotropies of \ion{H}{1} gas on subgalactic scales. Large differences in the physical scales on which we observe spatially variable direct escape Ly$\alpha$, blueshifted Ly$\alpha$, and escaping LyC photons in the Sunburst Arc underscore the importance of resolving the physical scales that govern Ly$\alpha$ and LyC escape., Comment: Submitted to The Astrophysical Journal with revisions from the first referee report. Comments welcome

Resonances, also known as quasi normal modes (QNM) in the non-Hermitian case, play an ubiquitous role in all domains of physics ruled by wave phenomena, notably in continuum mechanics, acoustics, electrodynamics, and quantum theory. In this paper, we present a QNM expansion for dispersive systems, recently applied to photonics but based on sixty year old techniques in mechanics. The resulting numerical algorithm appears to be physically agnostic, that is independent of the considered physical problem and can therefore be implemented as a mere toolbox in a nonlinear eigenvalue computation library., Comment: 23 pages, 6 figures, dedicated to Natasha and Sasha Movchan on the occasion of their jubilee

With the ongoing energy transition, power grids are evolving fast. They operate more and more often close to their technical limit, under more and more volatile conditions. Fast, essentially real-time computational approaches to evaluate their operational safety, stability and reliability are therefore highly desirable. Machine Learning methods have been advocated to solve this challenge, however they are heavy consumers of training and testing data, while historical operational data for real-world power grids are hard if not impossible to access. This manuscript describes an algorithmic approach for generating large datasets of power injections in electric power grids. The method allows one to generate arbitrarily large time series from the knowledge of the grid -- the admittance of its lines as well as the location, type and capacity of its power generators -- and aggregated power consumption data, such as the national load data given by ENTSO-E. The obtained datasets are statistically validated against real-world data., Comment: 15 pages, 8 figures, 2 tables. Dataset available at https://zenodo.org/records/13378476

We introduce AirLetters, a new video dataset consisting of real-world videos of human-generated, articulated motions. Specifically, our dataset requires a vision model to predict letters that humans draw in the air. Unlike existing video datasets, accurate classification predictions for AirLetters rely critically on discerning motion patterns and on integrating long-range information in the video over time. An extensive evaluation of state-of-the-art image and video understanding models on AirLetters shows that these methods perform poorly and fall far behind a human baseline. Our work shows that, despite recent progress in end-to-end video understanding, accurate representations of complex articulated motions -- a task that is trivial for humans -- remains an open problem for end-to-end learning., Comment: ECCV'24, HANDS workshop

Many successful methods to learn dynamical systems from data have recently been introduced. Such methods often rely on the availability of the system's full state. However, this underlying hypothesis is rather restrictive as it is typically not confirmed in practice, leaving us with partially observed systems. Utilizing the Mori-Zwanzig (MZ) formalism from statistical physics, we demonstrate that Constant Lag Neural Delay Differential Equations (NDDEs) naturally serve as suitable models for partially observed states. In empirical evaluation, we show that such models outperform existing methods on both synthetic and experimental data.

We introduce homological measurement, a framework for measuring the logical Pauli operators encoded in CSS stabilizer codes. The framework is based on the algebraic description of such codes as chain complexes. Protocols such as lattice surgery some of its recent generalizations are shown to be special cases of homological measurement. Using this framework, we develop a specific protocol called edge expanded homological measurement for fault-tolerant measurement of arbitrary logical Pauli operators of general qLDPC codes, requiring a number of ancillary qubits growing only linearly with the weight of the logical operator measured, and guaranteed that the distance of the code is preserved. We further benchmark our protocol numerically in a photonic architecture based on GKP qubits, showing that the logical error rate of various codes are on par with other methods requiring more ancilla qubits., Comment: 20 pages, 5 figures, 5 tables. Version 2 fixes a mistake in section 3.4 and briefly addresses related recent work

Circumbinary planets - planets that orbit both stars in a binary system - offer the opportunity to study planet formation and orbital migration in a different environment compare to single stars. However, despite the fact that > 90% of binary systems in the solar neighbourhood are spectrally resolved double-lined binaries, there has been only one detection of a circumbinary planet orbitting a double-lined binary using the radial velocity method so far. Spectrally disentangling both components of a binary system is hard to do accurately. Weak spectral lines blend with one another in a time-varying way, and inaccuracy in spectral modelling can lead to an inaccurate estimation of the radial-velocity of each component. This inaccuracy adds scatter to the measurements that can hide the weak radial-velocity signature of circumbinary exoplanets. We have obtained new high signal-to-noise and high-resolution spectra with the SOPHIE spectrograph, mounted on the 193cm telescope at Observatoire de Haute-Provence (OHP) for six, bright, double-lined binaries for which a circumbinary exoplanet detection has been attempted in the past. To extract radial-velocities we use the DOLBY code, a recent method of spectral disentangling using Gaussian processes to model the time-varying components. We analyse the resulting radial-velocities with a diffusive nested sampler to seek planets, and compute sensitivity limits. We do not detect any new circumbinary planet. However, we show that the combination of new data, new radial-velocity extraction methods, and new statistical methods to determine a dataset's sensitivity to planets leads to an approximately one order of magnitude improvement compared to previous results. This improvement brings us into the range of known circumbinary exoplanets and paves the way for new campaigns of observations targeting double-lined binaries., Comment: 19 pages, 8 figures, 7 tables. Accepted for publication in MNRAS

Mycetoma is a chronic and neglected inflammatory disease prevalent in tropical and subtropical regions. It can lead to severe disability and social stigma. The disease is classified into two types based on the causative microorganisms: eumycetoma (fungal) and actinomycetoma (bacterial). Effective treatment strategies depend on accurately identifying the causative agents. Current identification methods include molecular, cytological, and histopathological techniques, as well as grain culturing. Among these, histopathological techniques are considered optimal for use in endemic areas, but they require expert pathologists for accurate identification, which can be challenging in rural areas lacking such expertise. The advent of digital pathology and automated image analysis algorithms offers a potential solution. This report introduces a novel dataset designed for the automated detection and classification of mycetoma using histopathological images. It includes the first database of microscopic images of mycetoma tissue, detailing the entire pipeline from species distribution and patient sampling to acquisition protocols through histological procedures. The dataset consists of images from 142 patients, totalling 864 images, each annotated with binary masks indicating the presence of grains, facilitating both detection and segmentation tasks.

The Gemini Planet Imager (GPI) is a high contrast imaging instrument that aims to detect and characterize extrasolar planets. GPI is being upgraded to GPI 2.0, with several subsystems receiving a re-design to improve its contrast. To enable observations on fainter targets and increase performance on brighter ones, one of the upgrades is to the adaptive optics system. The current Shack-Hartmann wavefront sensor (WFS) is being replaced by a pyramid WFS with an low-noise electron multiplying CCD (EMCCD). EMCCDs are detectors capable of counting single photon events at high speed and high sensitivity. In this work, we characterize the performance of the HN\"u 240 EMCCD from N\"uv\"u Cameras, which was custom-built for GPI 2.0. Through our performance evaluation we found that the operating mode of the camera had to be changed from inverted-mode (IMO) to non-inverted mode (NIMO) in order to improve charge diffusion features found in the detector's images. Here, we characterize the EMCCD's noise contributors (readout noise, clock-induced charges, dark current) and linearity tests (EM gain, exposure time) before and after the switch to NIMO., Comment: Proceeding of the SPIE Astronomical Telescopes+Instrumentation. 14 pages, 15 figures

We consider dijet production in deep inelastic scattering at small $x$, on a purely gluonic unpolarized target. Starting from earlier results obtained at next-to-eikonal accuracy in the high-energy limit, we perform the expansion in the back-to-back dijet limit, at next-to-leading power accuracy. We rewrite our results in the language of transverse-momentum-dependent (TMD) factorization, in terms of twist-2 and twist-3 TMD gluon distributions (gluon TMDs). Among the next-to-eikonal corrections, we find in particular twist-2 contributions corresponding to the $x$ dependent phase of the twist-2 gluon TMDs. We also find two types of twist-3 unpolarized gluon TMDs, as well as correlators of three gluon field strength tensors., Comment: 38 pages

IJCLab has been leading the development and deployment of low-$\beta$ Superconducting Radio Frequency (SRF) cavities for proton and heavy ion accelerators. We are launching an electron accelerator project for sustainable Energy Recovery Linac (iSAS/PERLE) with state-of-the-art SRF cavities at 800~MHz. Our proposal includes advanced heat treatment of such cavities to reach an excellent quality factor of $3\times 10^{10}$ at $22$~MV/m. In this paper, we overview the status of this activity., Comment: Contribution to LCWS2024 proceedings

The computer vision community has explored dyadic interactions for atomic actions such as pushing, carrying-object, etc. However, with the advancement in deep learning models, there is a need to explore more complex dyadic situations such as loose interactions. These are interactions where two people perform certain atomic activities to complete a global action irrespective of temporal synchronisation and physical engagement, like cooking-together for example. Analysing these types of dyadic-interactions has several useful applications in the medical domain for social-skills development and mental health diagnosis. To achieve this, we propose a novel dual-path architecture to capture the loose interaction between two individuals. Our model learns global abstract features from each stream via a CNNs backbone and fuses them using a new Global-Layer-Attention module based on a cross-attention strategy. We evaluate our model on real-world autism diagnoses such as our Loose-Interaction dataset, and the publicly available Autism dataset for loose interactions. Our network achieves baseline results on the Loose-Interaction and SOTA results on the Autism datasets. Moreover, we study different social interactions by experimenting on a publicly available dataset i.e. NTU-RGB+D (interactive classes from both NTU-60 and NTU-120). We have found that different interactions require different network designs. We also compare a slightly different version of our method by incorporating time information to address tight interactions achieving SOTA results.

Quantum metrology, a cornerstone of quantum technologies, exploits entanglement and superposition to achieve higher precision than classical protocols in parameter estimation tasks. When combined with critical phenomena such as phase transitions, the divergence of quantum fluctuations is predicted to enhance the performance of quantum sensors. Here, we implement a critical quantum sensor using a superconducting parametric (i.e., two-photon driven) Kerr resonator. The sensor, a linear resonator terminated by a supercondicting quantum interference device, operates near the critical point of a finite-component second-order dissipative phase transition obtained by scaling the system parameters. We analyze the performance of a frequency-estimation protocol and show that quadratic precision scaling with respect to the system size can be achieved with finite values of the Kerr nonlinearity. Since each photon emitted from the cavity carries more information about the parameter to be estimated compared to its classical counterpart, our protocol opens perspectives for faster or more precise metrological protocols. Our results demonstrate that quantum advantage in a sensing protocol can be achieved by exploiting a finite-component phase transition.

The interfacial Dzyaloshinskii-Moriya interaction (DMI) can be exploited in magnetic thin films to realize lateral chirally coupled systems, providing a way to couple different sections of a magnetic racetrack and realize interconnected networks of magnetic logic gates. Here, we systematically investigate the interplay between spin-orbit torques, chiral coupling and the device design in domain wall racetracks. We show that the current-induced domain nucleation process can be tuned between single-domain nucleation and repeated nucleation of alternate domains by changing the orientation of an in-plane patterned magnetic region within an out-of-plane magnetic racetrack. Furthermore, by combining experiments and micromagnetic simulations, we show that the combination of damping-like and field-like spin-orbit torques with DMI results in selective domain wall injection in one of two arms of a Y-shaped devices depending on the current density. Such an element constitutes the basis of domain wall based demultiplexer, which is essential for distributing a single input to any one of the multiple outputs in logic circuits. Our results provide input for the design of reliable and multifunctional domain-wall circuits based on chirally coupled interfaces.

Bayesian calibration is an effective approach for ensuring that numerical simulations accurately reflect the behavior of physical systems. However, because numerical models are never perfect, a discrepancy known as model error exists between the model outputs and the observed data, and must be quantified. Conventional methods can not be implemented in transposition situations, such as when a model has multiple outputs but only one is experimentally observed. To account for the model error in this context, we propose augmenting the calibration process by introducing additional input numerical parameters through a hierarchical Bayesian model, which includes hyperparameters for the prior distribution of the calibration variables. Importance sampling estimators are used to avoid increasing computational costs. Performance metrics are introduced to assess the proposed probabilistic model and the accuracy of its predictions. The method is applied on a computer code with three outputs that models the Taylor cylinder impact test. The outputs are considered as the observed variables one at a time, to work with three different transposition situations. The proposed method is compared with other approaches that embed model errors to demonstrate the significance of the hierarchical formulation., Comment: Submitted to International Journal for Uncertainty Quantification

An equimolar NiCoFeCrGa high entropy alloy having dual-phase homogeneous components was studied, where the constituent phases exhibit distinct mechanical properties. Micropillars with various diameters were created from two differently heat treated samples, then they were compressed at slow strain rates, that revealed the material's limited sensitivity to size. On the other hand, increased strain rate sensitivity at high deformation speeds was observed, that differs substantially depending on the phase composition of the specimen. Dislocations within the two phases were studied by high resolution transmission electron microscopy and high angular resolution electron backscatter diffraction. The performed chemical analysis confirmed that slow cooling during casting create Cr-rich precipitates, that have significant impact on the global strength of the material.

Neuroscience employs diverse neuroimaging techniques, each offering distinct insights into brain activity, from electrophysiological recordings such as EEG, which have high temporal resolution, to hemodynamic modalities such as fMRI, which have increased spatial precision. However, integrating these heterogeneous data sources remains a challenge, which limits a comprehensive understanding of brain function. We present the Spatiotemporal Alignment of Multimodal Brain Activity (SAMBA) framework, which bridges the spatial and temporal resolution gaps across modalities by learning a unified latent space free of modality-specific biases. SAMBA introduces a novel attention-based wavelet decomposition for spectral filtering of electrophysiological recordings, graph attention networks to model functional connectivity between functional brain units, and recurrent layers to capture temporal autocorrelations in brain signal. We show that the training of SAMBA, aside from achieving translation, also learns a rich representation of brain information processing. We showcase this classify external stimuli driving brain activity from the representation learned in hidden layers of SAMBA, paving the way for broad downstream applications in neuroscience research and clinical contexts.

Cat qubits, a type of bosonic qubit encoded in a harmonic oscillator, can exhibit an exponential noise bias against bit-flip errors with increasing mean photon number. Here, we focus on cat qubits stabilized by two-photon dissipation, where pairs of photons are added and removed from a harmonic oscillator by an auxiliary, lossy buffer mode. This process requires a large loss rate and strong nonlinearities of the buffer mode that must not degrade the coherence and linearity of the oscillator. In this work, we show how to overcome this challenge by coloring the loss environment of the buffer mode with a multi-pole filter and optimizing the circuit to take into account additional inductances in the buffer mode. Using these techniques, we achieve near-ideal enhancement of cat-qubit bit-flip times with increasing photon number, reaching over $0.1$ seconds with a mean photon number of only $4$. Concurrently, our cat qubit remains highly phase coherent, with phase-flip times corresponding to an effective lifetime of $T_{1,\text{eff}} \simeq 70$ $\mu$s, comparable with the bare oscillator lifetime. We achieve this performance even in the presence of an ancilla transmon, used for reading out the cat qubit states, by engineering a tunable oscillator-ancilla dispersive coupling. Furthermore, the low nonlinearity of the harmonic oscillator mode allows us to perform pulsed cat-qubit stabilization, an important control primitive, where the stabilization can remain off for a significant fraction (e.g., two thirds) of a $3~\mathrm{\mu s}$ cycle without degrading bit-flip times. These advances are important for the realization of scalable error-correction with cat qubits, where large noise bias and low phase-flip error rate enable the use of hardware-efficient outer error-correcting codes., Comment: Comments welcome!

We present a novel approach for benchmarking and validating quantitative phase tomography (QPT) systems using three-dimensional microphantoms. These microphantoms, crafted from biological and imaging data, replicate the optical and structural properties of multicellular biological samples. Their fabrication featuring refractive index modulation at sub-micrometer details is enabled by two-photon polymerization. We showcase the effectiveness of our technique via a round-robin test of healthy and tumoral liver organoid phantoms across three different QPT systems. This test reveals sample- and system-dependent errors in measuring dry mass and morphology. Tailored microphantoms establish the gold standard for the optimization of hardware setups, reconstruction strategies and error assessment, paving the way for novel non-invasive, label-free measurements of 3D biological samples.

We consider random trigonometric polynomials with general dependent coefficients. We show that under mild hypotheses on the structure of dependence, the asymptotics as the degree goes to infinity of the expected number of real zeros coincides with the independent case. To the best of our knowledge, this universality result is the first obtained in a non-Gaussian dependent context. Our proof highlights the robustness of real zeros, even in the presence of dependencies. These findings bring the behavior of random polynomials closer to real-world models, where dependencies between coefficients are common., Comment: 29 pages

In quasi-fission, it is unclear what the interaction between the relative orbital angular momentum and the spin of the fragments is from a microscopic perspective. In macroscopic approaches, it is expected that the large value of the relative orbital angular momentum, of the order of 100~$\hbar$ is transferred through tangential dissipation to the fragments' intrinsic spin by sliding and rolling friction. The goal is to investigate the angular momentum transfer from the initial relative orbital angular momentum to the fragments' spin. How is the transferred spin shared between the fragments? What is the time scale associated with the different mechanisms? How does deformation play a role? A TDDFT simulation in the TDHF-Skyrme framework is used to describe several reactions at different impact parameters with increasing complexity. A method is proposed to study the evolution of the fragments' total spin as a function of time and angular velocity. The increasingly complex reaction allows for careful analysis of all the mechanisms responsible for the transfer of spin. In particular, it is shown that the transfer of nucleons, and neck formation can significantly affect the transfer of spin through tangential friction. Several mechanisms are in contradiction with previous macroscopic calculations. In particular, the spin of the fragments does not always increase during the collision which prevents it from being used to estimate the collision time of the reaction., Comment: 11 pages, 14 figures

We discuss in this article the creation and organization of a national core facility for bioimage analysis, based on a distributed team. F-BIAS federates bioimage analysts across France and relies on them to deliver services to the researchers of this territory. The main challenge in implementing this structure is to provide significant scientific value to its members, thereby encouraging their active participation and persuading their respective host teams to support their involvement. F-BIAS accomplished this by creating a professional network that mitigates the negative effects of isolation experienced by its members, who are often the sole bioimage analyst within their local teams, and fosters the development of their technical skills. In a second phase we capitalized on F-BIAS to create a virtual, remotely-operating core facility for bioimage analysis, offering consultations and collaborative project services to the scientific community of France. The services are organized so that they also contribute to the technical proficiency of the analysts. To promote the creation of similar structures, we present and discuss here the organization of this nationally distributed bioimage analysis service core, highlighting successes and challenges.

In this paper, a new classification model based on covariance matrices is built in order to classify buried objects. The inputs of the proposed models are the hyperbola thumbnails obtained with a classical Ground Penetrating Radar (GPR) system. These thumbnails are then inputs to the first layers of a classical CNN, which then produces a covariance matrix using the outputs of the convolutional filters. Next, the covariance matrix is given to a network composed of specific layers to classify Symmetric Positive Definite (SPD) matrices. We show in a large database that our approach outperform shallow networks designed for GPR data and conventional CNNs typically used in computer vision applications, particularly when the number of training data decreases and in the presence of mislabeled data. We also illustrate the interest of our models when training data and test sets are obtained from different weather modes or considerations.