Search

This Viewpoint discusses generative artificial intelligence and safeguarding privacy by using synthetic data as a substitute for private health data., Comment: JAMA Cardiology, 2024

Background: The MAJORANA DEMONSTRATOR , a modular array of isotopically enriched high-purity germanium (HPGe) detectors, was constructed to demonstrate backgrounds low enough to justify building a tonne-scale experiment to search for the neutrinoless double-beta decay ($\beta\beta(0\nu)$) of $^{76}\mathrm{Ge}$. Purpose: This paper presents a description of the instrument, its commissioning, and operations. It covers the electroforming, underground infrastructure, enrichment, detector fabrication, low-background and construction techniques, electronics, data acquisition, databases, and data processing of the MAJORANA DEMONSTRATOR. Method: The MAJORANA DEMONSTRATOR operated inside an ultra-low radioactivity passive shield at the 4850-foot~level of the Sanford Underground Research Facility (SURF) from 2015-2021. Results and Conclusions: The MAJORANA DEMONSTRATOR achieved the best energy resolution and second-best background level of any $\beta\beta(0\nu)$ search. This enabled it to achieve an ultimate half-life limit on $\beta\beta(0\nu)$ in $^{76}\mathrm{Ge}$ of $8.3\times 10^{25}$~yr (90\% C.L.) and perform a rich set of searches for other physics beyond the Standard Model., Comment: 72 pages

Sparse static detector networks in urban environments can be used in efforts to detect illicit radioactive sources, such as stolen nuclear material or radioactive "dirty bombs". We use detailed simulations to evaluate multiple configurations of detector networks and their ability to detect sources moving through a $6\times6$ km$^2$ area of downtown Chicago. A detector network's probability of detecting a source increases with detector density but can also be increased with strategic node placement. We show that the ability to fuse correlated data from a source-carrying vehicle passing by multiple detectors can significantly contribute to the overall detection probability. In this paper we distinguish static sensor deployments operated as networks able to correlate signals between sensors, from deployments operated as arrays where each sensor is operated individually. In particular, we show that additional visual attributes of source-carrying vehicles, such as vehicle color and make, can greatly improve the ability of a detector network to detect illicit sources.

The PAMPA (Point-Average-Moment PolynomiAl-interpreted) method, proposed in [R. Abgrall, Commun. Appl. Math. Comput., 5: 370-402, 2023], combines conservative and non-conservative formulations of hyperbolic conservation laws to evolve cell averages and point values. Solutions to hyperbolic conservation laws typically have an invariant domain, and ensuring numerical solutions stay within this domain is essential yet nontrivial. This paper presents a novel framework for designing efficient Invariant-Domain-Preserving (IDP) PAMPA schemes. We first analyze the IDP property for updated cell averages in the original PAMPA scheme, revealing the role of cell average decomposition and midpoint values in maintaining the invariant domain. This analysis highlights the difficulty of relying on continuous fluxes alone to preserve the invariant domain. Building on these insights, we introduce a simple IDP limiter for cell midpoint values, and propose a provably IDP PAMPA scheme that guarantees the preservation of the invariant domain for updated cell averages without requiring post-processing limiters. This approach contrasts with existing bound-preserving PAMPA schemes, which often require additional convex limiting to blend high-order and low-order solutions. Most notably, inspired by the Softplus and Clipped ReLU functions from machine learning, we propose an automatic IDP reformulation of the governing equations, resulting in an unconditionally limiter-free IDP scheme for evolving point values. We also introduce techniques to suppress spurious oscillations, enabling the scheme to capture strong shocks effectively. Numerical experiments on 1D problems, including the linear convection equation, Burgers equation, the compressible Euler equations, and MHD equations, demonstrate the accuracy and robustness of the proposed IDP PAMPA scheme., Comment: 28 pages

In this paper, we explore the use of the Virtual Element Method concepts to solve scalar and system hyperbolic problems on general polygonal grids. The new schemes stem from the active flux approach \cite{AF1}, which combines the usage of point values at the element boundaries with an additional degree of freedom representing the average of the solution within each control volume. Along the lines of the family of residual distribution schemes introduced in \cite{Abgrall_AF,abgrall2023activefluxtriangularmeshes} to bridge the active flux technique, we devise novel third order accurate methods that rely on the VEM technology to discretize gradients of the numerical solution by means of a polynomial-free approximation, hence adopting a virtual basis that is locally defined for each element. The obtained discretization is globally continuous, and for nonlinear problems it needs a stabilization which is provided by the \textit{a posteriori} MOOD paradigm \cite{Mood1}. This is applied to both point and average values of the discrete solution. We show applications to scalar problems, as well as to the acoustics and Euler equations in 2D. The accuracy and the robustness of the proposed schemes are assessed against a suite of benchmarks involving smooth solutions, shock waves and other discontinuities.

We present a novel positive kinetic scheme built on the efficient collide-and-stream algorithm of the lattice Boltzmann method (LBM) to address hyperbolic conservation laws. We focus on the compressible Euler equations with strong discontinuities. Starting from the work of Jin and Xin [20] and then [4,8], we show how the LBM discretization procedure can yield both first- and second-order schemes, referred to as vectorial LBM. Noticing that the first-order scheme is convex preserving under a specific CFL constraint, we develop a blending strategy that preserves both the conservation and simplicity of the algorithm. This approach employs convex limiters, carefully designed to ensure either positivity (of the density and the internal energy) preservation (PP) or well-defined local maximum principles (LMP), while minimizing numerical dissipation. On challenging test cases involving strong discontinuities and near-vacuum regions, we demonstrate the scheme accuracy, robustness, and ability to capture sharp discontinuities with minimal numerical oscillations., Comment: 23 pages, 10 figures

We provide a new geometric characterization for the spine of untwisted outer space of a right-angled Artin group, constructed by Charney, Stambaugh, and Vogtmann. We realize the spine as the natural simplicial complex associated with the category of spatial cube complexes, a new class of locally CAT(0) cube complexes with a simple definition in terms of hyperplane collapses., Comment: Corrected main definition. 26 pages

We propose a bound-preserving (BP) Point-Average-Moment PolynomiAl-interpreted (PAMPA) scheme by blending third-order and first-order constructions. The originality of the present construction is that it does not need any explicit reconstruction within each element, and therefore the construction is very flexible. The scheme employs a classical blending approach between a first-order BP scheme and a high-order scheme that does not inherently preserve bounds. The proposed BP PAMPA scheme demonstrates effectiveness across a range of problems, from scalar cases to systems such as the Euler equations of gas dynamics. We derive optimal blending parameters for both scalar and system cases, with the latter based on the recent geometric quasi-linearization (GQL) framework of [Wu \& Shu, {\em SIAM Review}, 65 (2023), pp. 1031--1073]. This yields explicit, optimal blending coefficients that ensure positivity and control spurious oscillations in both point values and cell averages. This framework incorporates a convex blending of fluxes and residuals from both high-order and first-order updates, facilitating a rigorous BP property analysis. Sufficient conditions for the BP property are established, ensuring robustness while preserving high-order accuracy. Numerical tests confirm the effectiveness of the BP PAMPA scheme on several challenging problems.

The ability to coherently control and read out qubits with long coherence times in a scalable system is a crucial requirement for any quantum processor. Nuclear spins in the solid state have shown great promise as long-lived qubits. Control and readout of individual nuclear spin qubit registers has made major progress in the recent years using individual electron spin ancilla addressed either electrically or optically. Here, we present a new platform for quantum information processing, consisting of $^{183}$W nuclear spin qubits adjacent to an Er$^{3+}$ impurity in a CaWO$_4$ crystal, interfaced via a superconducting resonator and detected using a microwave photon counter at 10mK. We study two nuclear spin qubits with $T_2^*$ of $0.8(2)~$s and $1.2(3)~$s, $T_2$ of $3.4(4)~$s and $4.4(6)~$ s, respectively. We demonstrate single-shot quantum non-demolition readout of each nuclear spin qubit using the Er$^{3+}$ spin as an ancilla. We introduce a new scheme for all-microwave single- and two-qubit gates, based on stimulated Raman driving of the coupled electron-nuclear spin system. We realize single- and two-qubit gates on a timescale of a few milliseconds, and prepare a decoherence-protected Bell state with 88% fidelity and $T_2^*$ of $1.7(2)~$s. Our results are a proof-of-principle demonstrating the potential of solid-state nuclear spin qubits as a promising platform for quantum information processing. With the potential to scale to tens or hundreds of qubits, this platform has prospects for the development of scalable quantum processors with long-lived qubits., Comment: 14 pages, 4 main figures, 7 supplementary figures

We characterize which groups splitting as finite graphs of free groups with cyclic edge groups are residually finite. Such a group $G$ is residually finite if and only if all its Baumslag-Solitar subgroups are residually finite. From a presentation of $G$, we construct a finite labeled graph $\Gamma$, and show that residual finiteness of $G$ is equivalent to an easily-detectable property of this graph. This characterization proves a conjecture of Wise., Comment: 36 pages

We report magnetic resonance spectroscopy measurements of individual nuclear spins in a crystal coupled to a neighbouring paramagnetic center, detected using microwave fluorescence at millikelvin temperatures. We observe real-time quantum jumps of the nuclear spin state, a proof of their individual nature. By driving the forbidden transitions of the coupled electron-nuclear spin system, we also achieve single-spin solid-effect dynamical nuclear polarization. Relying exclusively on microwave driving and microwave photon counting, the methods reported here are in principle applicable to a large number of electron-nuclear spin systems, in a wide variety of samples.

We present results of our study devoted to the development of a time correction algorithm needed to precisely synchronize a free-running Rubidium atomic clock with the Coordinated Universal Time (UTC). This R&D is performed in view of the Hyper-Kamiokande (HK) experiment currently under construction in Japan, which requires a synchronization with UTC and between its different experimental sites with a precision better than $100$ ns. We use a Global Navigation Satellite System (GNSS) receiver to compare a PPS and a $10$ MHz signal, generated by a free-running Rubidium clock, to the Global Positioning System (GPS) Time signal. We use these comparisons to correct the time series (time stamps) provided by the Rubidium clock signal. We fit the difference between Rubidium and GPS Time with polynomial functions of time over a certain integration time window to extract a correction of the Rubidium time stamps in offline or online mode. In online mode, the latest fit results are used for the correction until a new comparison to GPS Time becomes available. We show that with an integration time window of around $10^4$ seconds, we can correct the time stamps drift, caused by the frequency random walk noise and the deterministic frequency drift of the free running Rubidium clock, so that the time difference with respect to GPS Time stays within a $\pm5$ ns range in both offline or online correction mode. Presented results could be of interest for other experiments in the field of neutrino physics and multi-messenger astrophysics., Comment: 32 pages, 15 figures, manuscript prepared for submission to NIM A

Nowadays we observe an evolving landscape of data management and analytics, emphasising the significance of meticulous data management practices, semantic modelling, and bridging business-technical divides, to optimise data utilisation and enhance value from datasets in modern data environments. In this paper we introduce and explain the basic formalisation of the Semantic SQL Transducer, a well-founded but practical tool providing the materialised lossless conceptual view of an arbitrary relational source data, contributing to a knowledge-centric data stack.

This work presents a procedure to solve the Euler equations by explicitly updating, in a conservative manner, a generic thermodynamic variable such as temperature, pressure or entropy instead of the total energy. The presented procedure is valid for any equation of state and spatial discretization. When using complex equations of state such as Span-Wagner, choosing the temperature as the generic thermodynamic variable yields great reductions in the computational costs associated to thermodynamic evaluations. Results computed with a state of the art thermodynamic model are presented, and computational times are analyzed. Particular attention is dedicated to the conservation of total energy, the propagation speed of shock waves and jump conditions. The procedure is thoroughly tested using the Span-Wagner equation of state through the CoolProp thermodynamic library and the Van der Waals equation of state, both in the ideal and non-ideal compressible fluid-dynamics regimes, by comparing it to the standard total energy update and analytical solutions where available., Comment: Final version. Accepted version of: Computers & Fluids (2024) 279:106340 Published version available at 10.1016/j.compfluid.2024.106340

We propose a novel method for geolocalizing Unmanned Aerial Vehicles (UAVs) in environments lacking Global Navigation Satellite Systems (GNSS). Current state-of-the-art techniques employ an offline-trained encoder to generate a vector representation (embedding) of the UAV's current view, which is then compared with pre-computed embeddings of geo-referenced images to determine the UAV's position. Here, we demonstrate that the performance of these methods can be significantly enhanced by preprocessing the images to extract their edges, which exhibit robustness to seasonal and illumination variations. Furthermore, we establish that utilizing edges enhances resilience to orientation and altitude inaccuracies. Additionally, we introduce a confidence criterion for localization. Our findings are substantiated through synthetic experiments., Comment: Accepted for publication in IGARSS2024. 4 pages, 3 figures, 3 tables

Measurements of $K_S^0$ meson production via its $\pi^{+} \pi^{-}$ decay mode in inelastic $\textit{p+p}$ interactions at incident projectile momenta of 31, 40 and 80 GeV/$c$ ($\sqrt{s_{NN}}=7.7, 8.8$ and $12.3$ GeV, respectively) are presented. The data were recorded by the NA61/SHINE spectrometer at the CERN Super Proton Synchrotron. Double-differential distributions were obtained in transverse momentum and rapidity. The mean multiplicities of $K_S^0$ mesons were determined to be $(5.95 \pm 0.19 (stat) \pm 0.22 (sys)) \times 10^{-2}$ at 31 GeV/$c$, $(7.61 \pm 0.13 (stat) \pm 0.31 (sys)) \times 10^{-2}$ at 40 GeV/$c$ and $(11.58 \pm 0.12 (stat) \pm 0.37 (sys)) \times 10^{-2}$ at 80 GeV/$c$. The results on $K^{0}_{S}$ production are compared with model calculations (Epos1.99, SMASH 2.0 and PHSD) as well as with published data from other experiments., Comment: arXiv admin note: substantial text overlap with arXiv:2106.07535

We extend to multi-dimensions the work of [1], where new fully explicit kinetic methods were built for the approximation of linear and non-linear convection-diffusion problems. The fundamental principles from the earlier work are retained: (1) rather than aiming for the desired equations in the strict limit of a vanishing relaxation parameter, as is commonly done in the diffusion limit of kinetic methods, diffusion terms are sought as a first-order correction of this limit in a Chapman-Enskog expansion, (2) introducing a coupling between the conserved variables within the relaxation process by a specifically designed collision matrix makes it possible to systematically match a desired diffusion. Extending this strategy to multi-dimensions cannot, however, be achieved through simple directional splitting, as diffusion is likely to couple space directions with each other, such as with shear viscosity in the Navier-Stokes equations. In this work, we show how rewriting the collision matrix in terms of moments can address this issue, regardless of the number of kinetic waves, while ensuring conservation systematically. This rewriting allows for introducing a new class of kinetic models called \emph{regularized} models, simplifying the numerical methods and establishing connections with Jin-Xin models. Subsequently, new explicit arbitrary high-order kinetic schemes are formulated and validated on standard two-dimensional cases from the literature. Excellent results are obtained in the simulation of a shock-boundary layer interaction, validating their ability to approximate the Navier-Stokes equations with kinetic speeds obeying nothing but a subcharacteristic condition along with a hyperbolic constraint on the time step., Comment: 32 pages, 8 figures

In this article, we show how to construct a numerical method for solving hyperbolic problems, whether linear or nonlinear, using a continuous representation of the variables and their mean value in each triangular element. This type of approach has already been introduced by Roe, and others, in the multidimensional framework under the name of Active flux, see \cite{AF1,AF2,AF3,AF4,AF5}. Here, the presentation is more general and follows \cite{Abgrall_AF,BarzukowAbgrall}. { Various} examples show the good behavior of the method in both linear and nonlinear cases, including non-convex problems. The expected order of precision is obtained in both the linear and nonlinear cases. This work represents a step towards the development of methods in the spirit of virtual finite elements for linear or nonlinear hyperbolic problems, including the case where the solution is not regular.

We report the upgrade and operation of a frequency-comb-assisted high-resolution mid-infrared molecular spectrometer allowing us to combine high spectral purity, SI-traceability, wide tunability and high sensitivity. An optical frequency comb is used to transfer the spectral purity of a SI-traceable 1.54 $\mu$m metrology-grade frequency reference to a 10.3 $\mu$m quantum cascade laser (QCL). The near-infrared reference is operated at the French time/frequency metrology institute, calibrated there to primary frequency standards, and transferred to Laboratoire de Physique des Lasers via the REFIMEVE fiber network. The QCL exhibits a line width of $\delta\nu\sim 0.1$ Hz, a sub-$10^{-15}$ relative frequency stability from 0.1 to 10 s and its frequency is traceable to the SI with a total relative uncertainty better than $4\times10^{-14}$ after 1-s averaging time. We have developed the instrumentation allowing comb modes to be continuously tuned over 9 GHz resulting in a QCL of record spectral purity uninterruptedly tunable at the precision of the reference over an unprecedented span of $\Delta\nu$ = 1.4 GHz. We have used our apparatus to conduct sub-Doppler spectroscopy of methanol in a multi-pass cell, demonstrating state-of-art frequency uncertainties down to the few kilohertz level ($\sim10^{-10}$ in relative value). We have observed weak intensity resonances unreported so far, resolved subtle doublets never seen before and brought to light discrepancies with the HITRAN database. This demonstrates the potential of our apparatus for probing subtle internal molecular processes, building accurate spectroscopic models of polyatomic molecules of atmospheric or astrophysical interest, and carrying out precise spectroscopic tests of fundamental physics., Comment: 14 pages, 5 figures

We propose a numerical approach, of the BGK kinetic type, that is able to approximate with a given, but arbitrary, order of accuracy the solution of linear and non-linear convection-diffusion type problems: scalar advection-diffusion, non-linear scalar problems of this type and the compressible Navier-Stokes equations. Our kinetic model can use \emph{finite} advection speeds that are independent of the relaxation parameter, and the time step does not suffer from a parabolic constraint. Having finite speeds is in contrast with many of the previous works about this kind of approach, and we explain why this is possible: paraphrasing more or less \cite{golse:hal-00859451}, the convection-diffusion like PDE is not a limit of the BGK equation, but a correction of the same PDE without the parabolic term at the second order in the relaxation parameter that is interpreted as Knudsen number. We then show that introducing a matrix collision instead of the well-known BGK relaxation makes it possible to target a desired convection-diffusion system. Several numerical examples, ranging from a simple pure diffusion model to the compressible Navier-Stokes equations illustrate our approach

Active Flux is an extension of the Finite Volume method and additionally incorporates point values located at cell boundaries. This gives rise to a globally continuous approximation of the solution. Originally, the Active Flux method emerged as a fully discrete method, and required an exact or approximate evolution operator for the point value update. For nonlinear problems such an operator is often difficult to obtain, in particular for multiple spatial dimensions. We demonstrate that a new semi-discrete Active Flux method (first described in Abgrall&Barsukow, 2023 for one space dimension) can be used to solve nonlinear hyperbolic systems in multiple dimensions without requiring evolution operators. We focus here on the compressible Euler equations of inviscid hydrodynamics and third-order accuracy. We introduce a multi-dimensional limiting strategy and demonstrate the performance of the new method on both Riemann problems and subsonic flows.

The detection, identification, and localization of illicit radiological and nuclear material continue to be key components of nuclear non-proliferation and nuclear security efforts around the world. Networks of radiation detectors deployed at strategic locations in urban environments have the potential to provide continuous radiological/nuclear (R/N) surveillance and provide high probabilities of intercepting threat sources. The integration of contextual information from sensors such as video, Lidar, and meteorological sensors can provide significantly enhanced situational awareness, and improved detection and localization performance through the fusion of the radiological and contextual data. In this work, we present details of our work to establish a city-scale multi-sensor network testbed for intelligent, adaptive R/N detection in urban environments, and develop new techniques that enable city-scale source detection, localization, and tracking.

In this work, the high order accuracy and the well-balanced (WB) properties of some novel continuous interior penalty (CIP) stabilizations for the Shallow Water (SW) equations are investigated. The underlying arbitrary high order numerical framework is given by a Residual Distribution (RD)/continuous Galerkin (CG) finite element method (FEM) setting for the space discretization coupled with a Deferred Correction (DeC) time integration, to have a fully-explicit scheme. If, on the one hand, the introduced CIP stabilizations are all specifically designed to guarantee the exact preservation of the lake at rest steady state, on the other hand, some of them make use of general structures to tackle the preservation of general steady states, whose explicit analytical expression is not known. Several basis functions have been considered in the numerical experiments and, in all cases, the numerical results confirm the high order accuracy and the ability of the novel stabilizations to exactly preserve the lake at rest steady state and to capture small perturbations of such equilibrium. Moreover, some of them, based on the notions of space residual and global flux, have shown very good performances and superconvergences in the context of general steady solutions not known in closed-form. Many elements introduced here can be extended to other hyperbolic systems, e.g., to the Euler equations with gravity.

Since the celebrated theorem of Lax and Wendroff, we know a necessary condition that any numerical scheme for hyperbolic problem should satisfy: it should be written in flux form. A variant can also be formulated for the entropy. Even though some schemes, as for example those using continuous finite element, do not formally cast into this framework, it is a very convenient one. In this paper, we revisit this, introduce a different notion of local conservation which contains the previous one in one space dimension, and explore its consequences. This gives a more flexible framework that allows to get, systematically, entropy stable schemes, entropy dissipative ones, or accomodate more constraints. In particular, we can show that continuous finite element method can be rewritten in the finite volume framework, and all the quantities involved are explicitly computable. We end by presenting the only counter example we are aware of, i.e a scheme that seems not to be rewritten as a finite volume scheme., Comment: Proceedings of FVCA10, https://indico.math.cnrs.fr/event/8972/

In this paper, we introduce a new approach for constructing robust well-balanced numerical methods for the one-dimensional Saint-Venant system with and without the Manning friction term. Following the idea presented in [R. Abgrall, Commun. Appl. Math. Comput. 5(2023), pp. 370-402], we first combine the conservative and non-conservative (primitive) formulations of the studied conservative hyperbolic system in a natural way. The solution is globally continuous and described by a combination of point values and average values. The point values and average values will then be evolved by two different forms of PDEs: a conservative version of the cell averages and a possibly non-conservative one for the points. We show how to deal with both the conservative and non-conservative forms of PDEs in a well-balanced manner. The developed schemes are capable of exactly preserving both the still-water and moving-water equilibria. Compared with existing well-balanced methods, this new class of scheme is nonlinear-equations-solver-free. This makes the developed schemes less computationally costly and easier to extend to other models. We demonstrate the behavior of the proposed new scheme on several challenging examples.

We present a class of arbitrarily high order fully explicit kinetic numerical methods in compressible fluid dynamics, both in time and space, which include the relaxation schemes by S. Jin and Z. Xin. These methods can use CFL number larger or equal to unity on regular Cartesian meshes for multi-dimensional case. These kinetic models depend on a small parameter that can be seen as a "Knudsen" number. The method is asymptotic preserving in this Knudsen number. Also, the computational costs of the method are of the same order of a fully explicit scheme. This work is the extension of Abgrall et al. (2022) \cite{Abgrall} to multi-dimensional systems. We have assessed our method on several problems for two dimensional scalar problems and Euler equations and the scheme has proven to be robust and to achieve the theoretically predicted high order of accuracy on smooth solutions.

Static gamma-ray detector systems that are deployed outdoors for radiological monitoring purposes experience time- and spatially-varying natural backgrounds and encounters with man-made nuisance sources. In order to be sensitive to illicit sources, such systems must be able to distinguish those sources from benign variations due to, e.g., weather and human activity. In addition to fluctuations due to non-threats, each detector has its own response and energy resolution, so providing a large network of detectors with predetermined background and source templates can be an onerous task. Instead, we propose that static detectors use simple physics-informed algorithms to automatically learn the background and nuisance source signatures, which can them be used to bootstrap and feed into more complex algorithms. Specifically, we show that non-negative matrix factorization (NMF) can be used to distinguish static background from the effects of increased concentrations of radon progeny due to rainfall. We also show that a simple process of using multiple gross count rate filters can be used in real time to classify or ``triage'' spectra according to whether they belong to static, rain, or anomalous categories for processing with other algorithms. If a rain sensor is available, we propose a method to incorporate that signal as well. Two clustering methods for anomalous spectra are proposed, one using Kullback-Leibler divergence and the other using regularized NMF, with the goal of finding clusters of similar spectral anomalies that can be used to build anomaly templates. Finally we describe the issues involved in the implementation of some of these algorithms on deployed sensor nodes, including the need to monitor the background models for long-term drifting due to physical changes in the environment or changes in detector performance., Comment: 12 pages, 6 figures, accepted for publication in IEEE Transactions on Nuclear Science

In this paper we propose a novel way of deploying vulnerable architectures for defense and research purposes, which aims to generate deception platforms based on the formal description of a scenario. An attack scenario is described by an attack graph in which transitions are labeled by ATT&CK techniques or procedures. The state of the attacker is modeled as a set of secrets he acquires and a set of nodes he controls. Descriptions of a single scenario on a technical level can then be declined into several different scenarios on a procedural level, and each of these scenarios can be deployed into its own vulnerable architecture. To achieve this goal we introduce the notion of architecture constraints, as some procedures may only be exploited on system presenting special properties, such as having a specific operating system version. Finally, we present our deployment process for converting one of these scenarios into a vulnerable infrastructure, and offer an online proof of concept demonstration of our tool, where readers may deploy locally deploy a complete scenario inspired by the threat actor APT-29., Comment: 13 pages, 9 figures

We propose a novel Monte-Carlo based ab-initio algorithm for directly computing the statistics for quantities of interest in an immiscible two-phase compressible flow. Our algorithm samples the underlying probability space and evolves these samples with a sharp interface front-tracking scheme. Consequently, statistical information is generated without resorting to any closure assumptions and information about the underlying microstructure is implicitly included. The proposed algorithm is tested on a suite of numerical experiments and we observe that the ab-initio procedure can simulate a variety of flow regimes robustly and converges with respect of refinement of number of samples as well as number of bubbles per volume. The results are also compared with a state-of-the-art discrete equation method to reveal the inherent limitations of existing macroscopic models.

Deep neural networks do not discriminate between spurious and causal patterns, and will only learn the most predictive ones while ignoring the others. This shortcut learning behaviour is detrimental to a network's ability to generalize to an unknown test-time distribution in which the spurious correlations do not hold anymore. Debiasing methods were developed to make networks robust to such spurious biases but require to know in advance if a dataset is biased and make heavy use of minority counterexamples that do not display the majority bias of their class. In this paper, we argue that such samples should not be necessarily needed because the ''hidden'' causal information is often also contained in biased images. To study this idea, we propose 3 publicly released synthetic classification benchmarks, exhibiting predictive classification shortcuts, each of a different and challenging nature, without any minority samples acting as counterexamples. First, we investigate the effectiveness of several state-of-the-art strategies on our benchmarks and show that they do not yield satisfying results on them. Then, we propose an architecture able to succeed on our benchmarks, despite their unusual properties, using an entropic adversarial data augmentation training scheme. An encoder-decoder architecture is tasked to produce images that are not recognized by a classifier, by maximizing the conditional entropy of its outputs, and keep as much as possible of the initial content. A precise control of the information destroyed, via a disentangling process, enables us to remove the shortcut and leave everything else intact. Furthermore, results competitive with the state-of-the-art on the BAR dataset ensure the applicability of our method in real-life situations.

Context. Highly excited molecular hydrogen H2 has been observed in many regions of shocked molecular gas. A recently published $K$-band spectrum of Herbig-Haro 7 (HH7) contains several vibration-rotation lines of H2 from highly excited energy levels that have not been detected elsewhere, including a line at 2.179 $\mu$m identified as arising from the $v$=2 $J$=29 level, which lies above the dissociation limit of H2. One emission line at 2.104 $\mu$m in this spectrum was unidentified. Aims. We aim to complete the analysis of the spectrum of HH7 by including previously missing molecular data that have been recently computed. Methods. We re-analysed the $K$-band spectrum, emphasising the physics of quasi-bound upper levels that can produce infrared emission lines in the $K$ band. Results. We confirm the identification of the $2-1$ $S$(27) line at 2.1785 $\mu$m and identify the line at 2.1042 $\mu$m as due to the 1-0 $S$(29) transition of H2, whose upper level energy is also higher than the dissociation limit. This latter identification, its column density, and the energy of its upper level further substantiate the existence of a hot thermal component at 5000 K in the HH7 environment.} Conclusion. The presence of the newly identified $1-0$ $S$(29) line, whose quasi-bound upper level ($v$=1, $J$=31) has a significant spontaneous dissociation probability, shows that dissociation of H2 is occurring. The mechanism by which virtually all of the H2 levels with energies from 20,000 K to 53,000 K is maintained in local thermodynamic equilibrium at a single temperature of $\sim$5,000 K remains to be understood., Comment: 5 pages, 3 figures, accepted in Astronomy and Astrophysics Letters 22/12/2022

We present the two-dimensional unstructured grids extension of the a posteriori local subcell correction of discontinuous Galerkin (DG) schemes introduced in [52]. The technique is based on the reformulation of DG scheme as a finite volume (FV) like method through the definition of some specific numerical fluxes referred to as reconstructed fluxes. High-order DG numerical scheme combined with this new local subcell correction technique ensures positivity preservation of the solution, along with a low oscillatory and sharp shocks representation. The main idea of this correction procedure is to retain as much as possible the high accuracy and the very precise subcell resolution of DG schemes, while ensuring the robustness and stability of the numerical method, as well as preserving physical admissibility of the solution. Consequently, an a posteriori correction will only be applied locally at the subcell scale where it is needed, but still ensuring the local scheme conservation. Practically, at each time step, we compute a DG candidate solution and check if this solution is admissible (for instance positive, non-oscillating, etc). Numerical results on various type problems and test cases will be presented to assess the very good performance of the design correction algorithm.

Computing the distance function to some surface or line is a problem that occurs very frequently. There are several ways of computing a relevant approximation of this function, using for example technique originating from the approximation of Hamilton Jacobi problems, or the fast sweeping method. Here we make a link with some elliptic problem and propose a very fast way to approximate the distance function.

Deep neural networks often fail to generalize outside of their training distribution, in particular when only a single data domain is available during training. While test-time adaptation has yielded encouraging results in this setting, we argue that, to reach further improvements, these approaches should be combined with training procedure modifications aiming to learn a more diverse set of patterns. Indeed, test-time adaptation methods usually have to rely on a limited representation because of the shortcut learning phenomenon: only a subset of the available predictive patterns is learned with standard training. In this paper, we first show that the combined use of existing training-time strategies, and test-time batch normalization, a simple adaptation method, does not always improve upon the test-time adaptation alone on the PACS benchmark. Furthermore, experiments on Office-Home show that very few training-time methods improve upon standard training, with or without test-time batch normalization. We therefore propose a novel approach using a pair of classifiers and a shortcut patterns avoidance loss that mitigates the shortcut learning behavior by reducing the generalization ability of the secondary classifier, using the additional shortcut patterns avoidance loss that encourages the learning of samples specific patterns. The primary classifier is trained normally, resulting in the learning of both the natural and the more complex, less generalizable, features. Our experiments show that our method improves upon the state-of-the-art results on both benchmarks and benefits the most to test-time batch normalization.

We compute the radiative ro-vibrational emission spectrum of H2 involving quasibound states via a simple numerical method of resolution of the Schr\"odinger equation by introducing a modifed effective molecular potential. The comparison of the eigenvalues obtained with our approximation and other theoretical methods based on scattering resonance properties is excellent. Electric quadrupole and magnetic dipole contributions are calculated and we confirm the previous computations of Forrey of the electric quadrupole transition Einstein coeffcients. The astrophysical relevance of such quasibound levels is emphasized, Comment: 10 pages,1 figure, 4 tables. Molecular Physics

We propose an arbitrarily high-order accurate numerical method for conservation laws that is based on a continuous approximation of the solution. The degrees of freedom are point values at cell interfaces and moments of the solution inside the cell. To lowest ($3^\text{rd}$) order this method reduces to the Active Flux method. The update of the moments is achieved immediately by integrating the conservation law over the cell, integrating by parts and employing the continuity across cell interfaces. We propose two ways how the point values can be updated in time: either by first deriving a semi-discrete method that uses a finite-difference-type formula to approximate the spatial derivative, and integrating this method e.g. with a Runge-Kutta scheme, or by using a characteristics-based update, which is inspired by the original (fully discrete) Active Flux method. We analyze stability and accuracy of the resulting methods.

Active Flux is a recently developed numerical method for hyperbolic conservation laws. Its classical degrees of freedom are cell averages and point values at cell interfaces. These latter are shared between adjacent cells, leading to a globally continuous reconstruction. The update of the point values includes upwinding, but without solving a Riemann Problem. The update of the cell average requires a flux quadrature at the cell interface, which can be immediately performed using the point values. This paper explores different extensions of Active Flux to arbitrarily high order of accuracy, while maintaining the idea of global continuity. We propose to either increase the stencil while keeping the same degrees of freedom, or to increase the number of point values, or to include higher moments as new degrees of freedom. These extensions have different properties, and reflect different views upon the relation of Active Flux to the families of Finite Volume, Finite Difference and Finite Element methods.

In this work, we prove the convergence of residual distribution schemes to dissipative weak solutions of the Euler equations. We need to guarantee that the residual distribution schemes are fulfilling the underlying structure preserving properties such as positivity of density and internal energy. Consequently, the residual distribution schemes lead to a consistent and stable approximation of the Euler equations. Our result can be seen as a generalization of the Lax-Richtmyer equivalence theorem to nonlinear problems that consistency plus stability is equivalent to convergence., Comment: 28 pages, 3 figures

In this work we study various continuous finite element discretization for two dimensional hyperbolic partial differential equations, varying the polynomial space (Lagrangian on equispaced, Lagrangian on quadrature points (Cubature) and Bernstein), the stabilization techniques (streamline-upwind Petrov-Galerkin, continuous interior penalty, orthogonal subscale stabilization) and the time discretization (Runge-Kutta (RK), strong stability preserving RK and deferred correction). This is an extension of the one dimensional study by Michel S. et al J. Sci. Comput. (2021), whose results do not hold in multi-dimensional frameworks. The study ranks these schemes based on efficiency (most of them are mass-matrix free), stability and dispersion error, providing the best CFL and stabilization coefficients. The challenges in two-dimensions are related to the Fourier analysis. Here, we perform it on two types of periodic triangular meshes varying the angle of the advection, and we combine all the results for a general stability analysis. Furthermore, we introduce additional high order viscosity to stabilize the discontinuities, in order to show how to use these methods for tests of practical interest. All the theoretical results are thoroughly validated numerically both on linear and non-linear problems, and error-CPU time curves are provided. Our final conclusions suggest that Cubature elements combined with SSPRK and OSS stabilization is the most promising combination., Comment: arXiv admin note: text overlap with arXiv:2103.16158

The modeling of multi-phase flow is very challenging, given the range of scales as well as the diversity of flow regimes that one encounters in this context. We revisit the discrete equation method (DEM) for two-phase flow in the absence of heat conduction and mass transfer. We analyze the resulting probability coefficients and prove their local convexity, rigorously establishing that our version of DEM can model different flow regimes ranging from the disperse to stratified (or separated) flow. Moreover, we reformulate the underlying mesoscopic model in terms of an one-parameter family of PDEs that interpolates between different flow regimes. We also propose two sets of procedures to enforce relaxation to equilibrium. We perform several numerical tests to show the flexibility of the proposed formulation, as well as to interpret different model components. The one-parameter family of PDEs provides an unified framework for modeling mean quantities for a multiphase flow, while at the same time identifying two key parameters that model the inherent uncertainty in terms of the underlying microstructure.

In this work, we propose an accurate, robust, and stable discretization of the gamma-based compressible multicomponent model by Shyue [J. Comput. Phys., 142 (1998), 208-242] where each component follows a stiffened gas equation of state (EOS). We here extend the framework proposed in Renac [J. Comput. Phys. 382 (2019), 1-26] and Coquel et al. [J. Comput. Phys. 431 (2021) 110135] for the discretization of hyperbolic systems, with both fluxes and nonconservative products, to unstructured meshes with curved elements in multiple space dimensions. The framework relies on the discontinuous Galerkin spectral element method (DGSEM) using collocation of quadrature and interpolation points. We modify the integrals over discretization elements where we replace the physical fluxes and nonconservative products by two-point numerical fluctuations. The contributions of this work are threefold. First, we analyze the semi-discrete DGSEM discretization and prove that the scheme is high-order accurate, free-stream preserving, and entropy stable when excluding material interfaces. Second, we design a three-point scheme with a HLLC solver that does not require a root-finding algorithm for approximating the nonconservative products. The scheme is proved to be robust and entropy stable for convex entropies, preserves uniform states across material interfaces, satisfies a discrete minimum principle on the specific entropy and maximum principles on the EOS parameters. Third, the HLLC solver is applied at interfaces in the DGSEM scheme, while we consider two kinds of fluctuations in the integrals over discretization elements: material interface preserving and entropy conservative. Time integration is performed using SSP Runge-Kutta schemes. The high-order accuracy, nonlinear stability, and robustness of the present scheme are assessed through several numerical experiments in one and two space dimensions.

We present an experimental investigation of the stability limits specific to optical frequency standards using a fiber-optic architecture and semiconductor lasers. A compact setup composed of a semiconductor laser frequency-locked onto an acetylene transition detected in saturated absorption has been implemented using only fiber-optic components. Fiber optic technology allows compact and reliable solutions for various applications. However, for high sensitivity and stability applications such as metrology, residual reflections induced by optical index inhomogeneities in connectors and fiber-coupled components leading to spurious interference significantly limit performance. We have examined the origin of the interference fringes superimposed on the detected signal and the limitations they cause to the frequency stability of the reference. The effects of temperature and beam power fluctuations are also examined. Our results show that the frequency stability is limited in the 10-13 range by the effect of interference fringes due to use of fiber components.

We are interested in building schemes for the compressible Euler equations that are also locally conserving the angular momentum. We present a general framework, describe a few examples of schemes and show results. These schemes can be of arbitrary order.