Showing 1–50 of 285 results for author: Werner, F

H.E.S.S. observations of the 2021 periastron passage of PSR B1259-63/LS 2883

Authors: H. E. S. S. Collaboration, F. Aharonian, F. Ait Benkhali, J. Aschersleben, H. Ashkar, M. Backes, V. Barbosa Martins, R. Batzofin, Y. Becherini, D. Berge, K. Bernlöhr, M. Böttcher, C. Boisson, J. Bolmont, M. de Bony de Lavergne, J. Borowska, M. Bouyahiaoui, R. Brose, A. Brown, F. Brun, B. Bruno, T. Bulik, C. Burger-Scheidlin, S. Caroff, S. Casanova , et al. (119 additional authors not shown)

Abstract: PSR B1259-63 is a gamma-ray binary system that hosts a pulsar in an eccentric orbit, with a 3.4 year period, around an O9.5Ve star. At orbital phases close to periastron passages, the system radiates bright and variable non-thermal emission. We report on an extensive VHE observation campaign conducted with the High Energy Stereoscopic System, comprised of ~100 hours of data taken from $t_p-24$ day… ▽ More PSR B1259-63 is a gamma-ray binary system that hosts a pulsar in an eccentric orbit, with a 3.4 year period, around an O9.5Ve star. At orbital phases close to periastron passages, the system radiates bright and variable non-thermal emission. We report on an extensive VHE observation campaign conducted with the High Energy Stereoscopic System, comprised of ~100 hours of data taken from $t_p-24$ days to $t_p+127$ days around the system's 2021 periastron passage. We also present the timing and spectral analyses of the source. The VHE light curve in 2021 is consistent with the stacked light curve of all previous observations. Within the light curve, we report a VHE maximum at times coincident with the third X-ray peak first detected in the 2021 X-ray light curve. In the light curve -- although sparsely sampled in this time period -- we see no VHE enhancement during the second disc crossing. In addition, we see no correspondence to the 2021 GeV flare in the VHE light curve. The VHE spectrum obtained from the analysis of the 2021 dataset is best described by a power law of spectral index $Γ= 2.65 \pm 0.04_{\text{stat}}$ $\pm 0.04_{\text{sys}}$, a value consistent with the previous H.E.S.S. observations of the source. We report spectral variability with a difference of $ΔΓ= 0.56 ~\pm~ 0.18_{\text{stat}}$ $~\pm~0.10_{\text{sys}}$ at 95% c.l., between sub-periods of the 2021 dataset. We also find a linear correlation between contemporaneous flux values of X-ray and TeV datasets, detected mainly after $t_p+25$ days, suggesting a change in the available energy for non-thermal radiation processes. We detect no significant correlation between GeV and TeV flux points, within the uncertainties of the measurements, from $\sim t_p-23$ days to $\sim t_p+126$ days. This suggests that the GeV and TeV emission originate from different electron populations. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Comments: accepted to A&A

An overview of some single machine scheduling problems: polynomial algorithms, complexity and approximability

Authors: Nodari Vakhania, Frank Werner, Kevin Johedan Ramírez-Fuentes, Víctor Pacheco-Valencia

Abstract: Since the publication of the first scheduling paper in 1954, a huge number of works dealing with different types of single machine problems appeared. They addressed many heuristics and enumerative procedures, complexity results or structural properties of certain problems. Regarding surveys, often particular subjects like special objective functions are discussed, or more general scheduling proble… ▽ More Since the publication of the first scheduling paper in 1954, a huge number of works dealing with different types of single machine problems appeared. They addressed many heuristics and enumerative procedures, complexity results or structural properties of certain problems. Regarding surveys, often particular subjects like special objective functions are discussed, or more general scheduling problems were surveyed, where a substantial part is devoted to single machine problems. In this paper we present some results on polynomial algorithms, complexity and approximation issues, where the main focus is on results, which have been published during the last decades in papers, where at least one of the first two authors of this paper was involved. We hope that the reviewed results will stimulate further investigation in related research fields. △ Less

Submitted 17 June, 2024; v1 submitted 29 May, 2024; originally announced May 2024.

MSC Class: 90B35 (Primary); 90C59 (Secondary)

Wiener's Tauberian theorem in classical and quantum harmonic analysis

Authors: Robert Fulsche, Franz Luef, Reinhard F. Werner

Abstract: We investigate Wiener's Tauberian theorem from the perspective of limit functions, which results in several new versions of the Tauberian theorem. Based on this, we formulate and prove analogous Tauberian theorems for operators in the sense of quantum harmonic analysis. Using these results, we characterize the class of slowly oscillating operators and show that this class is strictly larger than t… ▽ More We investigate Wiener's Tauberian theorem from the perspective of limit functions, which results in several new versions of the Tauberian theorem. Based on this, we formulate and prove analogous Tauberian theorems for operators in the sense of quantum harmonic analysis. Using these results, we characterize the class of slowly oscillating operators and show that this class is strictly larger than the class of compact operators. Finally, we discuss uniform versions of Wiener's Tauberian theorem and its operator analogue and provide an application of this in operator theory. △ Less

Submitted 14 May, 2024; originally announced May 2024.

Comments: 36 pages, comments welcome

Accelerating Autonomy: Insights from Pro Racers in the Era of Autonomous Racing - An Expert Interview Study

Authors: Frederik Werner, René Oberhuber, Johannes Betz

Abstract: This research aims to investigate professional racing drivers' expertise to develop an understanding of their cognitive and adaptive skills to create new autonomy algorithms. An expert interview study was conducted with 11 professional race drivers, data analysts, and racing instructors from across prominent racing leagues. The interviews were conducted using an exploratory, non-standardized exper… ▽ More This research aims to investigate professional racing drivers' expertise to develop an understanding of their cognitive and adaptive skills to create new autonomy algorithms. An expert interview study was conducted with 11 professional race drivers, data analysts, and racing instructors from across prominent racing leagues. The interviews were conducted using an exploratory, non-standardized expert interview format guided by a set of prepared questions. The study investigates drivers' exploration strategies to reach their vehicle limits and contrasts them with the capabilities of state-of-the-art autonomous racing software stacks. Participants were questioned about the techniques and skills they have developed to quickly approach and maneuver at the vehicle limit, ultimately minimizing lap times. The analysis of the interviews was grounded in Mayring's qualitative content analysis framework, which facilitated the organization of the data into multiple categories and subcategories. Our findings create insights into human behavior regarding reaching a vehicle's limit and minimizing lap times. We conclude from the findings the development of new autonomy software modules that allow for more adaptive vehicle behavior. By emphasizing the distinct nuances between manual and autonomous driving techniques, the paper encourages further investigation into human drivers' strategies to maximize their vehicles' capabilities. △ Less

Submitted 4 May, 2024; originally announced May 2024.

Comments: 8 pages, 6 figures

Induced Subforests and Superforests

Authors: Dieter Rautenbach, Florian Werner

Abstract: Graph isomorphism, subgraph isomorphism, and maximum common subgraphs are classical well-investigated objects. Their (parameterized) complexity and efficiently tractable cases have been studied. In the present paper, for a given set of forests, we study maximum common induced subforests and minimum common induced superforests. We show that finding a maximum subforest is NP-hard already for two sub… ▽ More Graph isomorphism, subgraph isomorphism, and maximum common subgraphs are classical well-investigated objects. Their (parameterized) complexity and efficiently tractable cases have been studied. In the present paper, for a given set of forests, we study maximum common induced subforests and minimum common induced superforests. We show that finding a maximum subforest is NP-hard already for two subdivided stars while finding a minimum superforest is tractable for two trees but NP-hard for three trees. For a given set of $k$ trees, we present an efficient greedy $\left(\frac{k}{2}-\frac{1}{2}+\frac{1}{k}\right)$-approximation algorithm for the minimum superforest problem. Finally, we present a polynomial time approximation scheme for the maximum subforest problem for any given set of forests. △ Less

Submitted 21 March, 2024; originally announced March 2024.

Largest common subgraph of two forests

Authors: Dieter Rautenbach, Florian Werner

Abstract: A common subgraph of two graphs $G_1$ and $G_2$ is a graph that is isomorphic to subgraphs of $G_1$ and $G_2$. In the largest common subgraph problem the task is to determine a common subgraph for two given graphs $G_1$ and $G_2$ that is of maximum possible size ${\rm lcs}(G_1,G_2)$. This natural problem generalizes the well-studied graph isomorphism problem, has many applications, and remains NP-… ▽ More A common subgraph of two graphs $G_1$ and $G_2$ is a graph that is isomorphic to subgraphs of $G_1$ and $G_2$. In the largest common subgraph problem the task is to determine a common subgraph for two given graphs $G_1$ and $G_2$ that is of maximum possible size ${\rm lcs}(G_1,G_2)$. This natural problem generalizes the well-studied graph isomorphism problem, has many applications, and remains NP-hard even restricted to unions of paths. We present a simple $4$-approximation algorithm for forests, and, for every fixed $ε\in (0,1)$, we show that, for two given forests $F_1$ and $F_2$ of order at most $n$, one can determine in polynomial time a common subgraph $F$ of $F_1$ and $F_2$ with at least ${\rm lcs}(F_1,F_2)-εn$ edges. Restricted to instances with ${\rm lcs}(F_1,F_2)\geq cn$ for some fixed positive $c$, this yields a polynomial time approximation scheme. Our approach relies on the approximation of the given forests by structurally simpler forests that are composed of copies of only $O(\log (n))$ different starlike rooted trees and iterative quantizations of the options for the solutions. △ Less

Submitted 6 March, 2024; originally announced March 2024.

On k-ampleness equivalence

Authors: Laytimi Fatima Nahm Werner

Abstract: For a partition $a$ and a vector bundle $E$ on a projective variety $X$ let $\mathcal{F}l_s(E)$ be the corresponding flag manifold. There is a line bundle $\it Q_a^s$ on $\mathcal{F}l_s(E)$ with $p:\mathcal{F}l_s(E)\to X $ and $\it p_*Q_a^s = \mathcal{S}_aE$. We prove, if $\mathcal{S}_aE $ is $k$-ample (in the sense of Sommese) then $\it Q_a^s$ is $k$-ample. For the inverse if $\it Q_a^s$ is $k$-a… ▽ More For a partition $a$ and a vector bundle $E$ on a projective variety $X$ let $\mathcal{F}l_s(E)$ be the corresponding flag manifold. There is a line bundle $\it Q_a^s$ on $\mathcal{F}l_s(E)$ with $p:\mathcal{F}l_s(E)\to X $ and $\it p_*Q_a^s = \mathcal{S}_aE$. We prove, if $\mathcal{S}_aE $ is $k$-ample (in the sense of Sommese) then $\it Q_a^s$ is $k$-ample. For the inverse if $\it Q_a^s$ is $k$-ample, we prove that one of two the conditions of k-ampleness namely the cohomological vanishing is proved here but not yet the condition of semiamplenes of $\mathcal{S}_aE $ . △ Less

Submitted 12 February, 2024; originally announced February 2024.

MSC Class: 14F17

Maximum a posteriori testing in statistical inverse problems

Authors: Remo Kretschmann, Frank Werner

Abstract: This paper is concerned with a Bayesian approach to testing hypotheses in statistical inverse problems. Based on the posterior distribution $Π\left(\cdot |Y = y\right)$, we want to infer whether a feature $\langle\varphi, u^\dagger\rangle$ of the unknown quantity of interest $u^\dagger$ is positive. This can be done by the so-called maximum a posteriori test. We provide a frequentistic analysis of… ▽ More This paper is concerned with a Bayesian approach to testing hypotheses in statistical inverse problems. Based on the posterior distribution $Π\left(\cdot |Y = y\right)$, we want to infer whether a feature $\langle\varphi, u^\dagger\rangle$ of the unknown quantity of interest $u^\dagger$ is positive. This can be done by the so-called maximum a posteriori test. We provide a frequentistic analysis of this test's properties such as level and power, and prove that it is a regularized test in the sense of Kretschmann et al. (2024). Furthermore we provide lower bounds for its power under classical spectral source conditions in case of Gaussian priors. Numerical simulations illustrate its superior performance both in moderately and severely ill-posed situations. △ Less

Submitted 8 April, 2024; v1 submitted 1 February, 2024; originally announced February 2024.

MSC Class: 47A52; 62F15; 62G10; 62G20; 65J20; 65J22

Acceleration and transport of relativistic electrons in the jets of the microquasar SS 433

Authors: F. Aharonian, F. Ait Benkhali, J. Aschersleben, H. Ashkar, M. Backes, V. Barbosa Martins, R. Batzofin, Y. Becherini, D. Berge, K. Bernlöhr, B. Bi, M. Böttcher, C. Boisson, J. Bolmont, M. de Bony de Lavergne, J. Borowska, M. Bouyahiaou, M. Breuhau, R. Brose, A. M. Brown, F. Brun, B. Bruno, T. Bulik, C. Burger-Scheidlin, S. Caroff , et al. (140 additional authors not shown)

Abstract: SS 433 is a microquasar, a stellar binary system with collimated relativistic jets. We observed SS 433 in gamma rays using the High Energy Stereoscopic System (H.E.S.S.), finding an energy-dependent shift in the apparent position of the gamma-ray emission of the parsec-scale jets. These observations trace the energetic electron population and indicate the gamma rays are produced by inverse-Compton… ▽ More SS 433 is a microquasar, a stellar binary system with collimated relativistic jets. We observed SS 433 in gamma rays using the High Energy Stereoscopic System (H.E.S.S.), finding an energy-dependent shift in the apparent position of the gamma-ray emission of the parsec-scale jets. These observations trace the energetic electron population and indicate the gamma rays are produced by inverse-Compton scattering. Modelling of the energy-dependent gamma-ray morphology constrains the location of particle acceleration and requires an abrupt deceleration of the jet flow. We infer the presence of shocks on either side of the binary system at distances of 25 to 30 parsecs and conclude that self-collimation of the precessing jets forms the shocks, which then efficiently accelerate electrons. △ Less

Submitted 29 January, 2024; originally announced January 2024.

Comments: Submitted 20th Apr. 2023, published 25th January 2024 (accepted version)

Journal ref: Science383,402-406(2024)

Embezzlement of entanglement, quantum fields, and the classification of von Neumann algebras

Authors: Lauritz van Luijk, Alexander Stottmeister, Reinhard F. Werner, Henrik Wilming

Abstract: We study the quantum information theoretic task of embezzlement of entanglement in the setting of von Neumann algebras. Given a shared entangled resource state, this task asks to produce arbitrary entangled states using local operations without communication while perturbing the resource arbitrarily little. We quantify the performance of a given resource state by the worst-case error. States for w… ▽ More We study the quantum information theoretic task of embezzlement of entanglement in the setting of von Neumann algebras. Given a shared entangled resource state, this task asks to produce arbitrary entangled states using local operations without communication while perturbing the resource arbitrarily little. We quantify the performance of a given resource state by the worst-case error. States for which the latter vanishes are 'embezzling states' as they allow to embezzle arbitrary entangled states with arbitrarily small error. The best and worst performance among all states defines two algebraic invariants for von Neumann algebras. The first invariant takes only two values. Either it vanishes and embezzling states exist, which can only happen in type III, or no state allows for nontrivial embezzlement. In the case of factors not of finite type I, the second invariant equals the diameter of the state space. This provides a quantitative operational interpretation of Connes' classification of type III factors within quantum information theory. Type III$_1$ factors are 'universal embezzlers' where every state is embezzling. Our findings have implications for relativistic quantum field theory, where type III algebras naturally appear. For instance, they explain the maximal violation of Bell inequalities in the vacuum. Our results follow from a one-to-one correspondence between embezzling states and invariant probability measures on the flow of weights. We also establish that universally embezzling ITPFI factors are of type III$_1$ by elementary arguments. △ Less

Submitted 4 June, 2024; v1 submitted 14 January, 2024; originally announced January 2024.

Comments: See arXiv:2401.07292 for an overview article; 68 pages + 1 table + 1 figure; comments welcome; v3: resolved open problems

Embezzling entanglement from quantum fields

Authors: Lauritz van Luijk, Alexander Stottmeister, Reinhard F. Werner, Henrik Wilming

Abstract: Embezzlement of entanglement refers to the counterintuitive possibility of extracting entangled quantum states from a reference state of an auxiliary system (the "embezzler") via local quantum operations while hardly perturbing the latter. We uncover a deep connection between the operational task of embezzling entanglement and the mathematical classification of von Neumann algebras. Our result imp… ▽ More Embezzlement of entanglement refers to the counterintuitive possibility of extracting entangled quantum states from a reference state of an auxiliary system (the "embezzler") via local quantum operations while hardly perturbing the latter. We uncover a deep connection between the operational task of embezzling entanglement and the mathematical classification of von Neumann algebras. Our result implies that relativistic quantum fields are universal embezzlers: Any entangled state of any dimension can be embezzled from them with arbitrary precision. This provides an operational characterization of the infinite amount of entanglement present in the vacuum state of relativistic quantum field theories. △ Less

Submitted 4 June, 2024; v1 submitted 14 January, 2024; originally announced January 2024.

Comments: Overview paper accompanying the technical manuscript arXiv:2401.07299

Discovery of a Radiation Component from the Vela Pulsar Reaching 20 Teraelectronvolts

Authors: The H. E. S. S. Collaboration, :, F. Aharonian, F. Ait Benkhali, J. Aschersleben, H. Ashkar, M. Backes, V. Barbosa Martins, R. Batzofin, Y. Becherini, D. Berge, K. Bernlöhr, B. Bi, M. Böttcher, C. Boisson, J. Bolmont, M. de Bony de Lavergne, J. Borowska, F. Bradascio, M. Breuhaus, R. Brose, F. Brun, B. Bruno, T. Bulik, C. Burger-Scheidlin , et al. (157 additional authors not shown)

Abstract: Gamma-ray observations have established energetic isolated pulsars as outstanding particle accelerators and antimatter factories in the Galaxy. There is, however, no consensus regarding the acceleration mechanisms and the radiative processes at play, nor the locations where these take place. The spectra of all observed gamma-ray pulsars to date show strong cutoffs or a break above energies of a fe… ▽ More Gamma-ray observations have established energetic isolated pulsars as outstanding particle accelerators and antimatter factories in the Galaxy. There is, however, no consensus regarding the acceleration mechanisms and the radiative processes at play, nor the locations where these take place. The spectra of all observed gamma-ray pulsars to date show strong cutoffs or a break above energies of a few gigaelectronvolt (GeV). Using the H.E.S.S. array of Cherenkov telescopes, we discovered a novel radiation component emerging beyond this generic GeV cutoff in the Vela pulsar's broadband spectrum. The extension of gamma-ray pulsation energies up to at least 20 teraelectronvolts (TeV) shows that Vela pulsar can accelerate particles to Lorentz factors higher than $4\times10^7$. This is an order of magnitude larger than in the case of the Crab pulsar, the only other pulsar detected in the TeV energy range. Our results challenge the state-of-the-art models for high-energy emission of pulsars while providing a new probe, i.e. the energetic multi-TeV component, for constraining the acceleration and emission processes in their extreme energy limit. △ Less

Submitted 9 October, 2023; originally announced October 2023.

Comments: 38 pages, 6 figures. This preprint has not undergone peer review or any post-submission improvements or corrections. The Version of Record of this article is published in Nature Astronomy, Nat Astron (2023), and is available online at https://doi.org/10.1038/s41550-023-02052-3

Ensemble Laplacian Biogeography-Based Sine Cosine Algorithm for Structural Engineering Design Optimization Problems

Authors: Vanita Garg, Kusum Deep, Khalid Abdulaziz Alnowibet, Ali Wagdy Mohamed, Mohammad Shokouhifar, Frank Werner

Abstract: In this paper, an ensemble metaheuristic algorithm (denoted as LX-BBSCA) is introduced. It combines the strengths of Laplacian Biogeography-Based Optimization (LX-BBO) and the Sine Cosine Algorithm (SCA) to address structural engineering design optimization problems. Our primary objective is to mitigate the risk of getting stuck in local minima and accelerate the algorithm's convergence rate. We e… ▽ More In this paper, an ensemble metaheuristic algorithm (denoted as LX-BBSCA) is introduced. It combines the strengths of Laplacian Biogeography-Based Optimization (LX-BBO) and the Sine Cosine Algorithm (SCA) to address structural engineering design optimization problems. Our primary objective is to mitigate the risk of getting stuck in local minima and accelerate the algorithm's convergence rate. We evaluate the proposed LX-BBSCA algorithm on a set of 23 benchmark functions, including both unimodal and multimodal problems of varying complexity and dimensions. Additionally, we apply LX-BBSCA to tackle five real-world structural engineering design problems, comparing the results with those obtained using other metaheuristics in terms of objective function values and convergence behavior. To ensure the statistical validity of our findings, we employ rigorous tests such as the t-test and the Wilcoxon rank test. The experimental outcomes consistently demonstrate that the ensemble LX-BBSCA algorithm outperforms not only the basic versions of BBO, SCA, and LX-BBO but also other state-of-the-art metaheuristic algorithms. △ Less

Submitted 8 October, 2023; originally announced October 2023.

Comments: 25 pages, 9 tables, 5 figures

MSC Class: 90 C 59

Multiscale scanning with nuisance parameters

Authors: Claudia König, Axel Munk, Frank Werner

Abstract: We develop a multiscale scanning method to find anomalies in a $d$-dimensional random field in the presence of nuisance parameters. This covers the common situation that either the baseline-level or additional parameters such as the variance are unknown and have to be estimated from the data. We argue that state of the art approaches to determine asymptotically correct critical values for multisca… ▽ More We develop a multiscale scanning method to find anomalies in a $d$-dimensional random field in the presence of nuisance parameters. This covers the common situation that either the baseline-level or additional parameters such as the variance are unknown and have to be estimated from the data. We argue that state of the art approaches to determine asymptotically correct critical values for multiscale scanning statistics will in general fail when such parameters are naively replaced by plug-in estimators. Opposed to this, we suggest to estimate the nuisance parameters on the largest scale and to use (only) smaller scales for multiscale scanning. We prove a uniform invariance principle for the resulting adjusted multiscale statistic (AMS), which is widely applicable and provides a computationally feasible way to simulate asymptotically correct critical values. We illustrate the implications of our theoretical results in a simulation study and in a real data example from super-resolution STED microscopy. This allows us to identify interesting regions inside a specimen in a pre-scan with controlled family-wise error rate. △ Less

Submitted 30 April, 2024; v1 submitted 25 July, 2023; originally announced July 2023.

MSC Class: 60F17; 62H10; 60G50; 62F03

The Schmidt rank for the commuting operator framework

Authors: Lauritz van Luijk, René Schwonnek, Alexander Stottmeister, Reinhard F. Werner

Abstract: In quantum information theory, the Schmidt rank is a fundamental measure for the entanglement dimension of a pure bipartite state. Its natural definition uses the Schmidt decomposition of vectors on bipartite Hilbert spaces, which does not exist (or at least is not canonically given) if the observable algebras of the local systems are allowed to be general C*-algebras. In this work, we generalize… ▽ More In quantum information theory, the Schmidt rank is a fundamental measure for the entanglement dimension of a pure bipartite state. Its natural definition uses the Schmidt decomposition of vectors on bipartite Hilbert spaces, which does not exist (or at least is not canonically given) if the observable algebras of the local systems are allowed to be general C*-algebras. In this work, we generalize the Schmidt rank to the commuting operator framework where the joint system is not necessarily described by the minimal tensor product but by a general bipartite algebra. We give algebraic and operational definitions for the Schmidt rank and show their equivalence. We analyze bipartite states and compute the Schmidt rank in several examples: The vacuum in quantum field theory, Araki-Woods-Powers states, as well as ground states and translation invariant states on spin chains which are viewed as bipartite systems for the left and right half chains. We conclude with a list of open problems for the commuting operator framework. △ Less

Submitted 21 July, 2023; originally announced July 2023.

Comments: 44 pages, 3 figures

Journal ref: Commun. Math. Phys. 405, 152 (2024)

Scheduling on parallel machines with a common server in charge of loading and unloading operations

Authors: Abdelhak Elidrissi, Rachid Banmansour, Keramat Hasani, Frank Werner

Abstract: This paper addresses the scheduling problem on two identical parallel machines with a single server in charge of loading and unloading operations of jobs. Each job has to be loaded by the server before being processed on one of the two machines and unloaded by the same server after its processing. No delay is allowed between loading and processing, and between processing and unloading. The objecti… ▽ More This paper addresses the scheduling problem on two identical parallel machines with a single server in charge of loading and unloading operations of jobs. Each job has to be loaded by the server before being processed on one of the two machines and unloaded by the same server after its processing. No delay is allowed between loading and processing, and between processing and unloading. The objective function involves the minimization of the makespan. This problem referred to as P2, S1|sj , tj |Cmax generalizes the classical parallel machine scheduling problem with a single server which performs only the loading (i.e., setup) operation of each job. For this NP-hard problem, no solution algorithm was proposed in the literature. Therefore, we present two mixedinteger linear programming (MILP) formulations, one with completion-time variables along with two valid inequalities and one with time-indexed variables. In addition, we propose some polynomial-time solvable cases and a tight theoretical lower bound. In addition, we show that the minimization of the makespan is equivalent to the minimization of the total idle times on the machines. To solve large-sized instances of the problem, an efficient General Variable Neighborhood Search (GVNS) metaheuristic with two mechanisms for finding an initial solution is designed. The GVNS is evaluated by comparing its performance with the results provided by the MILPs and another metaheuristic. The results show that the average percentage deviation from the theoretical lower-bound of GVNS is within 0.642%. Some managerial insights are presented and our results are compared with the related literature. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 40 pages, 4 figures, 16 tables

MSC Class: 90B35

Convergence of Dynamics on Inductive Systems of Banach Spaces

Authors: Lauritz van Luijk, Alexander Stottmeister, Reinhard F. Werner

Abstract: Many features of physical systems, both qualitative and quantitative, become sharply defined or tractable only in some limiting situation. Examples are phase transitions in the thermodynamic limit, the emergence of classical mechanics from quantum theory at large action, and continuum quantum field theory arising from renormalization group fixed points. It would seem that few methods can be useful… ▽ More Many features of physical systems, both qualitative and quantitative, become sharply defined or tractable only in some limiting situation. Examples are phase transitions in the thermodynamic limit, the emergence of classical mechanics from quantum theory at large action, and continuum quantum field theory arising from renormalization group fixed points. It would seem that few methods can be useful in such diverse applications. However, we here present a flexible modeling tool for the limit of theories: soft inductive limits constituting a generalization of inductive limits of Banach spaces. In this context, general criteria for the convergence of dynamics will be formulated, and these criteria will be shown to apply in the situations mentioned and more. △ Less

Submitted 6 July, 2023; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: Comments welcome

Journal ref: Annales Henri Poincaré (2024)

Mostar index and bounded maximum degree

Authors: Michael A. Henning, Johannes Pardey, Dieter Rautenbach, Florian Werner

Abstract: Došlić et al. defined the Mostar index of a graph $G$ as $Mo(G)=\sum\limits_{uv\in E(G)}|n_G(u,v)-n_G(v,u)|$, where, for an edge $uv$ of $G$, the term $n_G(u,v)$ denotes the number of vertices of $G$ that have a smaller distance in $G$ to $u$ than to $v$. For a graph $G$ of order $n$ and maximum degree at most $Δ$, we show $Mo(G)\leq \fracΔ{2}n^2-(1-o(1))c_Δn\log(\log(n)),$ where $c_Δ>0$ only depe… ▽ More Došlić et al. defined the Mostar index of a graph $G$ as $Mo(G)=\sum\limits_{uv\in E(G)}|n_G(u,v)-n_G(v,u)|$, where, for an edge $uv$ of $G$, the term $n_G(u,v)$ denotes the number of vertices of $G$ that have a smaller distance in $G$ to $u$ than to $v$. For a graph $G$ of order $n$ and maximum degree at most $Δ$, we show $Mo(G)\leq \fracΔ{2}n^2-(1-o(1))c_Δn\log(\log(n)),$ where $c_Δ>0$ only depends on $Δ$ and the $o(1)$ term only depends on $n$. Furthermore, for integers $n_0$ and $Δ$ at least $3$, we show the existence of a $Δ$-regular graph of order $n$ at least $n_0$ with $Mo(G)\geq \fracΔ{2}n^2-c'_Δn\log(n),$ where $c'_Δ>0$ only depends on $Δ$. △ Less

Submitted 15 June, 2023; originally announced June 2023.

Adaptive minimax optimality in statistical inverse problems via SOLIT -- Sharp Optimal Lepskii-Inspired Tuning

Authors: Housen Li, Frank Werner

Abstract: We consider statistical linear inverse problems in separable Hilbert spaces and filter-based reconstruction methods of the form $\hat f_α= q_α\left(T^*T\right)T^*Y$, where $Y$ is the available data, $T$ the forward operator, $\left(q_α\right)_{α\in \mathcal A}$ an ordered filter, and $α> 0$ a regularization parameter. Whenever such a method is used in practice, $α$ has to be appropriately chosen.… ▽ More We consider statistical linear inverse problems in separable Hilbert spaces and filter-based reconstruction methods of the form $\hat f_α= q_α\left(T^*T\right)T^*Y$, where $Y$ is the available data, $T$ the forward operator, $\left(q_α\right)_{α\in \mathcal A}$ an ordered filter, and $α> 0$ a regularization parameter. Whenever such a method is used in practice, $α$ has to be appropriately chosen. Typically, the aim is to find or at least approximate the best possible $α$ in the sense that mean squared error (MSE) $\mathbb E [\Vert \hat f_α- f^\dagger\Vert^2]$ w.r.t.~the true solution $f^\dagger$ is minimized. In this paper, we introduce the Sharp Optimal Lepskiĭ-Inspired Tuning (SOLIT) method, which yields an a posteriori parameter choice rule ensuring adaptive minimax rates of convergence. It depends only on $Y$ and the noise level $σ$ as well as the operator $T$ and the filter $\left(q_α\right)_{α\in \mathcal A}$ and does not require any problem-dependent tuning of further parameters. We prove an oracle inequality for the corresponding MSE in a general setting and derive the rates of convergence in different scenarios. By a careful analysis we show that no other a posteriori parameter choice rule can yield a better performance in terms of the order of the convergence rate of the MSE. In particular, our results reveal that the typical understanding of Lepski\uı-type methods in inverse problems leading to a loss of a log factor is wrong. In addition, the empirical performance of SOLIT is examined in simulations. △ Less

Submitted 11 December, 2023; v1 submitted 20 April, 2023; originally announced April 2023.

Comments: Some technical parts are polished, and a comparison with classical Lepskii is included in the simulation section

Journal ref: Inverse Problems, Volume 40, Number 2, 2024

Variational Poisson Denoising via Augmented Lagrangian Methods

Authors: Christian Kanzow, Fabius Krämer, Patrick Mehlitz, Gerd Wachsmuth, Frank Werner

Abstract: In this paper, we denoise a given noisy image by minimizing a smoothness promoting function over a set of local similarity measures which compare the mean of the given image and some candidate image on a large collection of subboxes. The associated convex optimization problem possesses a huge number of constraints which are induced by extended real-valued functions stemming from the Kullback--Leib… ▽ More In this paper, we denoise a given noisy image by minimizing a smoothness promoting function over a set of local similarity measures which compare the mean of the given image and some candidate image on a large collection of subboxes. The associated convex optimization problem possesses a huge number of constraints which are induced by extended real-valued functions stemming from the Kullback--Leibler divergence. Alternatively, these nonlinear constraints can be reformulated as affine ones, which makes the model seemingly more tractable. For the numerical treatment of both formulations of the model (i.e., the original one as well as the one with affine constraints), we propose a rather general augmented Lagrangian method which is capable of handling the huge amount of constraints. A self-contained, derivative-free, global convergence theory is provided, allowing an extension to other problem classes. For the solution of the resulting subproblems in the setting of our suggested image denoising models, we make use of a suitable stochastic gradient method. Results of several numerical experiments are presented in order to compare both formulations and the associated augmented Lagrangian methods. △ Less

Submitted 21 June, 2024; v1 submitted 13 April, 2023; originally announced April 2023.

Comments: 38 pages, 6 figures, 3 tables

MSC Class: 49M37; 90C30; 90C48; 90C90

A Contribution of the HAWC Observatory to the TeV era in the High Energy Gamma-Ray Astrophysics: The case of the TeV-Halos

Authors: Ramiro Torres-Escobedo, Hao Zhou, Eduardo de la Fuente, A. U. Abeysekara, A. Albert, R. Alfaro, C. Alvarez, J. D. Álvarez, J. R. Angeles Camacho, J. C. Arteaga-Velázquez, K. P. Arunbabu, D. Avila Rojas, H. A. Ayala Solares, R. Babu, V. Baghmanyan, A. S. Barber, J. Becerra Gonzalez, E. Belmont-Moreno, S. Y. BenZvi, D. Berley, C. Brisbois, K. S. Caballero-Mora, T. Capistrán, A. Carramiñana, S. Casanova , et al. (108 additional authors not shown)

Abstract: We present a short overview of the TeV-Halos objects as a discovery and a relevant contribution of the High Altitude Water Čerenkov (HAWC) observatory to TeV astrophysics. We discuss history, discovery, knowledge, and the next step through a new and more detailed analysis than the original study in 2017. TeV-Halos will contribute to resolving the problem of the local positron excess observed on th… ▽ More We present a short overview of the TeV-Halos objects as a discovery and a relevant contribution of the High Altitude Water Čerenkov (HAWC) observatory to TeV astrophysics. We discuss history, discovery, knowledge, and the next step through a new and more detailed analysis than the original study in 2017. TeV-Halos will contribute to resolving the problem of the local positron excess observed on the Earth. To clarify the latter, understanding the diffusion process is mandatory. △ Less

Submitted 13 April, 2023; originally announced April 2023.

Comments: Work presented in the 21st International Symposium on Very High Energy Cosmic Ray Interactions(ISVHECRI 2022) as part of the Ph. D. Thesis of Ramiro Torres-Escobedo (SJTU, Shanghai, China). Accepted for publication in SciPost Physics Proceedings (ISSN 2666-4003). 11 pages, 3 Figures. Short overview of HAWC and TeV Halos objects until 2022

Connected and Autonomous Vehicle Scheduling Problems: Some Models and Algorithms

Authors: Evgeny R. Gafarov, Frank Werner

Abstract: In this paper, we consider scheduling problems that arise in connected and autonomous vehicle systems. For four variants of such problems, mathematical models and solution algorithms are presented. In particular, three polynomial algorithms and a branch and bound algorithms are developed. In this paper, we consider scheduling problems that arise in connected and autonomous vehicle systems. For four variants of such problems, mathematical models and solution algorithms are presented. In particular, three polynomial algorithms and a branch and bound algorithms are developed. △ Less

Submitted 4 April, 2023; originally announced April 2023.

Comments: 8 pages, 4 figures

MSC Class: 90B35; 90C27; 68Q25; 68W40

The High-Altitude Water Cherenkov (HAWC) Observatory in México: The Primary Detector

Authors: A. U. Abeysekara, A. Albert, R. Alfaro, C. Álvarez, J. D. Álvarez, M. Araya, J. C. Arteaga-Velázquez, K. P. Arunbabu, D. Avila Rojas, H. A. Ayala Solares, R. Babu, A. S. Barber, A. Becerril, E. Belmont-Moreno, S. Y. BenZvi, O. Blanco, J. Braun, C. Brisbois, K. S. Caballero-Mora, J. I. Cabrera Martínez, T. Capistrán, A. Carramiñana, S. Casanova, M. Castillo, O. Chaparro-Amaro , et al. (118 additional authors not shown)

Abstract: The High-Altitude Water Cherenkov (HAWC) observatory is a second-generation continuously operated, wide field-of-view, TeV gamma-ray observatory. The HAWC observatory and its analysis techniques build on experience of the Milagro experiment in using ground-based water Cherenkov detectors for gamma-ray astronomy. HAWC is located on the Sierra Negra volcano in México at an elevation of 4100 meters a… ▽ More The High-Altitude Water Cherenkov (HAWC) observatory is a second-generation continuously operated, wide field-of-view, TeV gamma-ray observatory. The HAWC observatory and its analysis techniques build on experience of the Milagro experiment in using ground-based water Cherenkov detectors for gamma-ray astronomy. HAWC is located on the Sierra Negra volcano in México at an elevation of 4100 meters above sea level. The completed HAWC observatory principal detector (HAWC) consists of 300 closely spaced water Cherenkov detectors, each equipped with four photomultiplier tubes to provide timing and charge information to reconstruct the extensive air shower energy and arrival direction. The HAWC observatory has been optimized to observe transient and steady emission from sources of gamma rays within an energy range from several hundred GeV to several hundred TeV. However, most of the air showers detected are initiated by cosmic rays, allowing studies of cosmic rays also to be performed. This paper describes the characteristics of the HAWC main array and its hardware. △ Less

Submitted 10 April, 2023; v1 submitted 3 April, 2023; originally announced April 2023.

Comments: Accepted for publications in Nuclear Inst. and Methods in Physics Research, A (2023) 168253 ( https://www.sciencedirect.com/science/article/abs/pii/S0168900223002437 ); 39 pages, 14 Figures

Journal ref: Nuclear Inst. and Methods in Physics Research, A (2023) 168253

Irregularity of Graphs respecting Degree Bounds

Authors: Dieter Rautenbach, Florian Werner

Abstract: Albertson defined the irregularity of a graph $G$ as $irr(G)=\sum\limits_{uv\in E(G)}|d_G(u)-d_G(v)|$. For a graph $G$ with $n$ vertices, $m$ edges, maximum degree $Δ$, and $d=\left\lfloor \frac{Δm}{Δn-m}\right\rfloor$, we show $$irr(G)\leq d(d+1)n+\frac{1}Δ\left(Δ^2-(2d+1)Δ-d^2-d\right)m.$$ Albertson defined the irregularity of a graph $G$ as $irr(G)=\sum\limits_{uv\in E(G)}|d_G(u)-d_G(v)|$. For a graph $G$ with $n$ vertices, $m$ edges, maximum degree $Δ$, and $d=\left\lfloor \frac{Δm}{Δn-m}\right\rfloor$, we show $$irr(G)\leq d(d+1)n+\frac{1}Δ\left(Δ^2-(2d+1)Δ-d^2-d\right)m.$$ △ Less

Submitted 22 March, 2023; originally announced March 2023.

Detailed Analysis of the TeV γ-Ray Sources 3HWC J1928+178, 3HWC J1930+188, and the New Source HAWC J1932+192

Authors: A. Albert, R. Alfaro, C. Alvarez, J. C. Arteaga-Velázquez, D. Avila Rojas, H. A. Ayala Solares, R. Babu, E. Belmont-Moreno, C. Brisbois, K. S. Caballero-Mora, T. Capistrń, A. Carramiñana, S. Casanova, O. Chaparro-Amaro, U. Cotti, J. Cotzomi, S. CoutiñodeLeón, E. De la Fuente, C. de León, R. Diaz Hernandez, J. C. Díaz-Vélez, B. L. Dingus, M. A. DuVernois, M. Durocher, K. Engel , et al. (69 additional authors not shown)

Abstract: The latest High Altitude Water Cherenkov (HAWC) point-like source catalog up to 56 TeV reported the detection of two sources in the region of the Galactic plane at galactic longitude 52°< l < 55°, 3HWC J1930+188 and 3HWC J1928+178. The first one is associated with a known TeV source, the supernova remnant SNR G054.1+00.3. It was discovered by one of the currently operating Imaging Atmospheric Cher… ▽ More The latest High Altitude Water Cherenkov (HAWC) point-like source catalog up to 56 TeV reported the detection of two sources in the region of the Galactic plane at galactic longitude 52°< l < 55°, 3HWC J1930+188 and 3HWC J1928+178. The first one is associated with a known TeV source, the supernova remnant SNR G054.1+00.3. It was discovered by one of the currently operating Imaging Atmospheric Cherenkov Telescope (IACT), the Very Energetic Radiation Imaging Telescope Array System (VERITAS), detected by the High Energy Stereoscopic System (H.E.S.S.), and identified as a composite SNR. However, the source 3HWC J1928+178, discovered by HAWC and coincident with the pulsar PSR J1928+1746, was not detected by any IACT despite their long exposure on the region, until a recent new analysis of H.E.S.S. data was able to confirm it. Moreover, no X-ray counterpart has been detected from this pulsar. We present a multicomponent fit of this region using the latest HAWC data. This reveals an additional new source, HAWC J1932+192, which is potentially associated with the pulsar PSR J1932+1916, whose gamma-ray emission could come from the acceleration of particles in its pulsar wind nebula. In the case of 3HWC J1928+178, several possible explanations are explored, in a attempt to unveil the origins of the very-high-energy gamma-ray emission. △ Less

Submitted 27 February, 2023; originally announced February 2023.

Covariant catalysis requires correlations and good quantum reference frames degrade little

Authors: Lauritz van Luijk, Reinhard F. Werner, Henrik Wilming

Abstract: Catalysts are quantum systems that open up dynamical pathways between quantum states which are otherwise inaccessible under a given set of operational restrictions while, at the same time, they do not change their quantum state. We here consider the restrictions imposed by symmetries and conservation laws, where any quantum channel has to be covariant with respect to the unitary representation of… ▽ More Catalysts are quantum systems that open up dynamical pathways between quantum states which are otherwise inaccessible under a given set of operational restrictions while, at the same time, they do not change their quantum state. We here consider the restrictions imposed by symmetries and conservation laws, where any quantum channel has to be covariant with respect to the unitary representation of a symmetry group, and present two results. First, for an exact catalyst to be useful, it has to build up correlations to either the system of interest or the degrees of freedom dilating the given process to covariant unitary dynamics. This explains why catalysts in pure states are useless. Second, if a quantum system ("reference frame") is used to simulate to high precision unitary dynamics (which possibly violates the conservation law) on another system via a global, covariant quantum channel, then this channel can be chosen so that the reference frame is approximately catalytic. In other words, a reference frame that simulates unitary dynamics to high precision degrades only very little. △ Less

Submitted 25 October, 2023; v1 submitted 24 January, 2023; originally announced January 2023.

Comments: 12+1 pages. Accepted for Publication at Quantum - the open journal for quantum science

Journal ref: Quantum 7, 1166 (2023)

Optimal regularized hypothesis testing in statistical inverse problems

Authors: Remo Kretschmann, Daniel Wachsmuth, Frank Werner

Abstract: Testing of hypotheses is a well studied topic in mathematical statistics. Recently, this issue has also been addressed in the context of Inverse Problems, where the quantity of interest is not directly accessible but only after the inversion of a (potentially) ill-posed operator. In this study, we propose a regularized approach to hypothesis testing in Inverse Problems in the sense that the underl… ▽ More Testing of hypotheses is a well studied topic in mathematical statistics. Recently, this issue has also been addressed in the context of Inverse Problems, where the quantity of interest is not directly accessible but only after the inversion of a (potentially) ill-posed operator. In this study, we propose a regularized approach to hypothesis testing in Inverse Problems in the sense that the underlying estimators (or test statistics) are allowed to be biased. Under mild source-condition type assumptions we derive a family of tests with prescribed level $α$ and subsequently analyze how to choose the test with maximal power out of this family. As one major result we prove that regularized testing is always at least as good as (classical) unregularized testing. Furthermore, using tools from convex optimization, we provide an adaptive test by maximizing the power functional, which then outperforms previous unregularized tests in numerical simulations by several orders of magnitude. △ Less

Submitted 17 October, 2023; v1 submitted 25 December, 2022; originally announced December 2022.

MSC Class: 62G10; 47A52; 65J20; 65F22; 65R30

Journal ref: Inverse Problems 40:015013 (2024)

On uniqueness and ill-posedness for the deautoconvolution problem in the multi-dimensional case

Authors: Bernd Hofmann, Frank Werner, Yu Deng

Abstract: This paper analyzes the inverse problem of deautoconvolution in the multi-dimensional case with respect to solution uniqueness and ill-posedness. Deautoconvolution means here the reconstruction of a real-valued $L^2$-function with support in the $n$-dimensional unit cube $[0,1]^n$ from observations of its autoconvolution either in the full data case (i.e. on $[0,2]^n$) or in the limited data case… ▽ More This paper analyzes the inverse problem of deautoconvolution in the multi-dimensional case with respect to solution uniqueness and ill-posedness. Deautoconvolution means here the reconstruction of a real-valued $L^2$-function with support in the $n$-dimensional unit cube $[0,1]^n$ from observations of its autoconvolution either in the full data case (i.e. on $[0,2]^n$) or in the limited data case (i.e. on $[0,1]^n$). Based on multi-dimensional variants of the Titchmarsh convolution theorem due to Lions and Mikusiński, we prove in the full data case a twofoldness assertion, and in the limited data case uniqueness of non-negative solutions for which the origin belongs to the support. The latter assumption is also shown to be necessary for any uniqueness statement in the limited data case. A glimpse of rate results for regularized solutions completes the paper. △ Less

Submitted 13 December, 2022; originally announced December 2022.

Three-body contact for fermions. I. General relations

Authors: Félix Werner, Xavier Leyronas

Abstract: We consider the resonant Fermi gas, that is, two-component fermions in three dimensions interacting by a short-range potential of large scattering length. We introduce a quantity, the three-body contact, that determines several observables. Within the zero-range model, the number of nearby fermion triplets, the large-momentum tail of the center-of-mass momentum distribution of nearby fermion pairs… ▽ More We consider the resonant Fermi gas, that is, two-component fermions in three dimensions interacting by a short-range potential of large scattering length. We introduce a quantity, the three-body contact, that determines several observables. Within the zero-range model, the number of nearby fermion triplets, the large-momentum tail of the center-of-mass momentum distribution of nearby fermion pairs, as well as the large-momentum tail of the two-particle momentum distribution, are expressed in terms of the three-body contact. For a small finite interaction range, the formation rate of deeply bound dimers by three-body recombination, as well as the three-body contribution to the finite-range correction to the energy, are expressed in terms of the three-body contact and of a three-body parameter. This three-body parameter, which vanishes in the zero-range limit, is defined through the asymptotic behavior of the zero-energy scattering state at distances intermediate between the range and the two-body scattering length. In general, the three-body contact has different contributions labeled by spin and angular momentum indices, and the three-body parameter can depend on those indices. We also include the generalization to unequal masses for $\uparrow$ and $\downarrow$ particles. With respect to the relation between three-body loss rate and number of nearby triplets stated in [Petrov, Salomon and Shlyapnikov, PRL 93, 090404 (2004)], the present work adds a derivation, expresses the proportionality factor in terms of the three-body parameter, and includes the general case where there are several contributions to the three-body contact and several three-body parameters. △ Less

Submitted 17 April, 2024; v1 submitted 17 November, 2022; originally announced November 2022.

Comments: Final version, typos corrected [additions in v2: relations on large-momentum tail of 2-particle momentum distribution and 3-body contribution to finite-range energy-correction; mass-imbalanced case; corrected and extended derivation of 3-body loss rate; links with works on 3-body losses (by Petrov et al.) and on scattering hypervolume (by Tan et al); value of C_3 for 3 trapped particles]

Journal ref: C. R. Phys. 25, 179 (2024)

Bounding the Mostar index

Authors: Štefko Miklavič, Johannes Pardey, Dieter Rautenbach, Florian Werner

Abstract: Došlić et al. defined the Mostar index of a graph $G$ as $Mo(G)=\sum\limits_{uv\in E(G)}|n_G(u,v)-n_G(v,u)|$, where, for an edge $uv$ of $G$, the term $n_G(u,v)$ denotes the number of vertices of $G$ that have a smaller distance in $G$ to $u$ than to $v$. They conjectured that $Mo(G)\leq 0.\overline{148}n^3$ for every graph $G$ of order $n$. As a natural upper bound on the Mostar index, Geneson an… ▽ More Došlić et al. defined the Mostar index of a graph $G$ as $Mo(G)=\sum\limits_{uv\in E(G)}|n_G(u,v)-n_G(v,u)|$, where, for an edge $uv$ of $G$, the term $n_G(u,v)$ denotes the number of vertices of $G$ that have a smaller distance in $G$ to $u$ than to $v$. They conjectured that $Mo(G)\leq 0.\overline{148}n^3$ for every graph $G$ of order $n$. As a natural upper bound on the Mostar index, Geneson and Tsai implicitly consider the parameter $Mo^\star(G)=\sum\limits_{uv\in E(G)}\big(n-\min\{ d_G(u),d_G(v)\}\big)$. For a graph $G$ of order $n$, they show that $Mo^\star(G)\leq \frac{5}{24}(1+o(1))n^3$. We improve this bound to $Mo^\star(G)\leq \left(\frac{2}{\sqrt{3}}-1\right)n^3$, which is best possible up to terms of lower order. Furthermore, we show that $Mo^\star(G)\leq \left(2\left(\fracΔ{n}\right)^2+\left(\fracΔ{n}\right)-2\left(\fracΔ{n}\right)\sqrt{\left(\fracΔ{n}\right)^2+\left(\fracΔ{n}\right)}\right)n^3$ provided that $G$ has maximum degree $Δ$. △ Less

Submitted 12 November, 2022; originally announced November 2022.

Symmetry-broken perturbation theory to large orders in antiferromagnetic phases

Authors: R. Garioud, F. Šimkovic IV, R. Rossi, G. Spada, T. Schäfer, F. Werner, M. Ferrero

Abstract: We introduce a spin-symmetry-broken extension of the connected determinant algorithm [Phys. Rev. Lett. 119, 045701 (2017)]. The resulting systematic perturbative expansions around an antiferromagnetic state allow for numerically exact calculations directly inside a magnetically ordered phase. We show new precise results for the magnetic phase diagram and thermodynamics of the three-dimensional cub… ▽ More We introduce a spin-symmetry-broken extension of the connected determinant algorithm [Phys. Rev. Lett. 119, 045701 (2017)]. The resulting systematic perturbative expansions around an antiferromagnetic state allow for numerically exact calculations directly inside a magnetically ordered phase. We show new precise results for the magnetic phase diagram and thermodynamics of the three-dimensional cubic Hubbard model at half-filling. With detailed computations of the order parameter in the low to intermediate-coupling regime, we establish the N{é}el phase boundary. The critical behavior in its vicinity is shown to be compatible with the $O(3)$ Heisenberg universality class. By determining the evolution of the entropy with decreasing temperature through the phase transition we identify the different physical regimes at $U/t\!=\!4$. We provide quantitative results for several thermodynamic quantities deep inside the antiferromagnetic dome up to large interaction strengths and investigate the crossover between the Slater and Heisenberg regimes. △ Less

Submitted 27 June, 2023; v1 submitted 31 October, 2022; originally announced October 2022.

Comments: 13 pages, 10 figures

Journal ref: Phys. Rev. Lett. 132, 246505 (2024)

Deautoconvolution in the two-dimensional case

Authors: Yu Deng, Bernd Hofmann, Frank Werner

Abstract: There is extensive mathematical literature on the inverse problem of deautoconvolution for a function with support in the unit interval $[0,1] \subset \mathbb R$, but little is known about the multidimensional situation. This article tries to fill this gap with analytical and numerical studies on the reconstruction of a real function of two real variables over the unit square from observations of… ▽ More There is extensive mathematical literature on the inverse problem of deautoconvolution for a function with support in the unit interval $[0,1] \subset \mathbb R$, but little is known about the multidimensional situation. This article tries to fill this gap with analytical and numerical studies on the reconstruction of a real function of two real variables over the unit square from observations of its autoconvolution on $[0,2]^2 \subset \mathbb R^2$ (full data case) or on $[0,1]^2$ (limited data case). In an $L^2$-setting, twofoldness and uniqueness assertions are proven for the deautoconvolution problem in 2D. Moreover, its ill-posedness is characterized and illustrated. Extensive numerical case studies give an overview of the behaviour of stable approximate solutions to the two-dimensional deautoconvolution problem obtained by Tikhonov-type regularization with different penalties and the iteratively regularized Gauss-Newton method. △ Less

Submitted 25 October, 2022; originally announced October 2022.

Automated dysgraphia detection by deep learning with SensoGrip

Authors: Mugdim Bublin, Franz Werner, Andrea Kerschbaumer, Gernot Korak, Sebastian Geyer, Lena Rettinger, Erna Schoenthaler, Matthias Schmid-Kietreiber

Abstract: Dysgraphia, a handwriting learning disability, has a serious negative impact on children's academic results, daily life and overall wellbeing. Early detection of dysgraphia allows for an early start of a targeted intervention. Several studies have investigated dysgraphia detection by machine learning algorithms using a digital tablet. However, these studies deployed classical machine learning algo… ▽ More Dysgraphia, a handwriting learning disability, has a serious negative impact on children's academic results, daily life and overall wellbeing. Early detection of dysgraphia allows for an early start of a targeted intervention. Several studies have investigated dysgraphia detection by machine learning algorithms using a digital tablet. However, these studies deployed classical machine learning algorithms with manual feature extraction and selection as well as binary classification: either dysgraphia or no dysgraphia. In this work, we investigated fine grading of handwriting capabilities by predicting SEMS score (between 0 and 12) with deep learning. Our approach provide accuracy more than 99% and root mean square error lower than one, with automatic instead of manual feature extraction and selection. Furthermore, we used smart pen called SensoGrip, a pen equipped with sensors to capture handwriting dynamics, instead of a tablet, enabling writing evaluation in more realistic scenarios. △ Less

Submitted 3 January, 2023; v1 submitted 14 October, 2022; originally announced October 2022.

Maximizing the Mostar index for bipartite graphs and split graphs

Authors: Štefko Miklavič, Johannes Pardey, Dieter Rautenbach, Florian Werner

Abstract: Došlić et al.~defined the Mostar index of a graph $G$ as $\sum\limits_{uv\in E(G)}|n_G(u,v)-n_G(v,u)|$, where, for an edge $uv$ of $G$, the term $n_G(u,v)$ denotes the number of vertices of $G$ that have a smaller distance in $G$ to $u$ than to $v$. Contributing to conjectures posed by Došlić et al., we show that the Mostar index of bipartite graphs of order $n$ is at most… ▽ More Došlić et al.~defined the Mostar index of a graph $G$ as $\sum\limits_{uv\in E(G)}|n_G(u,v)-n_G(v,u)|$, where, for an edge $uv$ of $G$, the term $n_G(u,v)$ denotes the number of vertices of $G$ that have a smaller distance in $G$ to $u$ than to $v$. Contributing to conjectures posed by Došlić et al., we show that the Mostar index of bipartite graphs of order $n$ is at most $\frac{\sqrt{3}}{18}n^3$, and that the Mostar index of split graphs of order $n$ is at most $\frac{4}{27}n^3$. △ Less

Submitted 7 October, 2022; originally announced October 2022.

A Double Layered Water Cherenkov Detector Array for Gamma-Ray Astronomy

Authors: Samridha Kunwar, Hazal Goksu, Jim Hinton, Harm Schoorlemmer, Andrew Smith, Werner Hofmann, Felix Werner

Abstract: Ground-level particle detection is now a well-established approach to TeV gamma-ray astronomy. Detection of Cherenkov light produced in water-filled detection units is a proven and cost-effective method. Here we discuss the optimization of the units towards the future Southern Wide-field Gamma-ray Observatory (SWGO). In this context, we investigate a new type of configuration in which each water C… ▽ More Ground-level particle detection is now a well-established approach to TeV gamma-ray astronomy. Detection of Cherenkov light produced in water-filled detection units is a proven and cost-effective method. Here we discuss the optimization of the units towards the future Southern Wide-field Gamma-ray Observatory (SWGO). In this context, we investigate a new type of configuration in which each water Cherenkov detector (WCD) unit in the array comprises two chambers with black or reflective walls and a single photomultiplier tube (PMT) in each chamber. We find that this is a cost-effective approach that improves the performance of the WCD array with respect to current approaches. A shallow lower chamber with a PMT facing downwards enables muon tagging and the identification of hadron-induced air showers, which are the primary source of background in gamma-ray astronomy. We investigate how gamma/hadron separation power and achievable angular resolution depend on the geometry and wall reflectivity of the detector units in this configuration. We find that excellent angular resolution, background rejection power and low-energy response are achievable in this double-layer configuration, with the aid of reflective surfaces in both chambers. △ Less

Submitted 23 January, 2023; v1 submitted 19 September, 2022; originally announced September 2022.

Comments: 17 pages, 20 figures

Efficiency Evaluation of Banks with Many Branches using a Heuristic Framework and Dynamic Data Envelopment Optimization Approach: A Real Case Study

Authors: Vahid Kayvanfar, Hamed Baziyad, Shaya Sheikh, Frank Werner

Abstract: Evaluating the efficiency of organizations and branches within an organization is a challenging issue for managers. Evaluation criteria allow organizations to rank their internal units, identify their position concerning their competitors, and implement strategies for improvement and development purposes. Among the methods that have been applied in the evaluation of bank branches, non-parametric m… ▽ More Evaluating the efficiency of organizations and branches within an organization is a challenging issue for managers. Evaluation criteria allow organizations to rank their internal units, identify their position concerning their competitors, and implement strategies for improvement and development purposes. Among the methods that have been applied in the evaluation of bank branches, non-parametric methods have captured the attention of researchers in recent years. One of the most widely used non-parametric methods is the data envelopment analysis (DEA) which leads to promising results. However, the static DEA approaches do not consider the time in the model. Therefore, this paper uses a dynamic DEA (DDEA) method to evaluate the branches of a private Iranian bank over three years (2017-2019). The results are then compared with static DEA. After ranking the branches, they are clustered using the K-means method. Finally, a comprehensive sensitivity analysis approach is introduced to help the managers to decide about changing variables to shift a branch from one cluster to a more efficient one. △ Less

Submitted 11 September, 2022; originally announced September 2022.

Comments: 22 pages, 4 figures, 11 tables

MSC Class: 68T09; 90C59; 91C20 ACM Class: E.m

Irrational quantum walks

Authors: Gabriel Coutinho, Pedro Ferreira Baptista, Chris Godsil, Thomás Jung Spier, Reinhard Werner

Abstract: The adjacency matrix of a graph G is the Hamiltonian for a continuous-time quantum walk on the vertices of G. Although the entries of the adjacency matrix are integers, its eigenvalues are generally irrational and, because of this, the behaviour of the walk is typically not periodic. In consequence we can usually only compute numerical approximations to parameters of the walk. In this paper, we de… ▽ More The adjacency matrix of a graph G is the Hamiltonian for a continuous-time quantum walk on the vertices of G. Although the entries of the adjacency matrix are integers, its eigenvalues are generally irrational and, because of this, the behaviour of the walk is typically not periodic. In consequence we can usually only compute numerical approximations to parameters of the walk. In this paper, we develop theory to exactly study any quantum walk generated by an integral Hamiltonian. As a result, we provide exact methods to compute the average of the mixing matrices, and to decide whether pretty good (or almost) perfect state transfer occurs in a given graph. We also use our methods to study geometric properties of beautiful curves arising from entries of the quantum walk matrix, and discuss possible applications of these results. △ Less

Submitted 18 August, 2022; originally announced August 2022.

Comments: 22 pages, 5 figures. Code to test pretty good state transfer is available upon request

MSC Class: 81P68; 05C50; 42A75

Quantum-Classical Hybrid Systems and their Quasifree Transformations

Authors: Lars Dammeier, Reinhard F. Werner

Abstract: We study continuous variable systems, in which quantum and classical degrees of freedom are combined and treated on the same footing. Thus all systems, including the inputs or outputs to a channel, may be quantum-classical hybrids. This allows a unified treatment of a large variety of quantum operations involving measurements or dependence on classical parameters. The basic variables are given by… ▽ More We study continuous variable systems, in which quantum and classical degrees of freedom are combined and treated on the same footing. Thus all systems, including the inputs or outputs to a channel, may be quantum-classical hybrids. This allows a unified treatment of a large variety of quantum operations involving measurements or dependence on classical parameters. The basic variables are given by canonical operators with scalar commutators. Some variables may commute with all others and hence generate a classical subsystem. We systematically study the class of "quasifree" operations, which are characterized equivalently either by an intertwining condition for phase-space translations or by the requirement that, in the Heisenberg picture, Weyl operators are mapped to multiples of Weyl operators. This includes the well-known Gaussian operations, evolutions with quadratic Hamiltonians, and "linear Bosonic channels", but allows for much more general kinds of noise. For example, all states are quasifree. We sketch the analysis of quasifree preparation, measurement, repeated observation, cloning, teleportation, dense coding, the setup for the classical limit, and some aspects of irreversible dynamics, together with the precise salient tradeoffs of uncertainty, error, and disturbance. Although the spaces of observables and states are infinite dimensional for every non-trivial system that we consider, we treat the technicalities related to this in a uniform and conclusive way, providing a calculus that is both easy to use and fully rigorous. △ Less

Submitted 19 July, 2023; v1 submitted 9 August, 2022; originally announced August 2022.

Comments: 63 pages, 6 figures

Journal ref: Quantum 7, 1068 (2023)

On a Dynamic Variant of the Iteratively Regularized Gauss-Newton Method with Sequential Data

Authors: Neil K. Chada, Marco A. Iglesias, Shuai Lu, Frank Werner

Abstract: For numerous parameter and state estimation problems, assimilating new data as they become available can help produce accurate and fast inference of unknown quantities. While most existing algorithms for solving those kind of ill-posed inverse problems can only be used with a single instance of the observed data, in this work we propose a new framework that enables existing algorithms to invert mu… ▽ More For numerous parameter and state estimation problems, assimilating new data as they become available can help produce accurate and fast inference of unknown quantities. While most existing algorithms for solving those kind of ill-posed inverse problems can only be used with a single instance of the observed data, in this work we propose a new framework that enables existing algorithms to invert multiple instances of data in a sequential fashion. Specifically we will work with the well-known iteratively regularized Gauss-Newton method (IRGNM), a variational methodology for solving nonlinear inverse problems. We develop a theory of convergence analysis for a proposed dynamic IRGNM algorithm in the presence of Gaussian white noise. We combine this algorithm with the classical IRGNM to deliver a practical (hybrid) algorithm that can invert data sequentially while producing fast estimates. Our work includes the proof of well-definedness of the proposed iterative scheme, as well as various error bounds that rely on standard assumptions for nonlinear inverse problems. We use several numerical experiments to verify our theoretical findings, and to highlight the benefits of incorporating sequential data. The context of the numerical experiments comprises various parameter identification problems including a Darcy flow elliptic PDE example, and that of electrical impedance tomography. △ Less

Submitted 27 July, 2022; originally announced July 2022.

Towards quantitative super-resolution microscopy: Molecular maps with statistical guarantees

Authors: Katharina Proksch, Frank Werner, Jan Keller-Findeisen, Haisen Ta, Axel Munk

Abstract: Quantifying the number of molecules from fluorescence microscopy measurements is an important topic in cell biology and medical research. In this work, we present a consecutive algorithm for super-resolution (STED) scanning microscopy that provides molecule counts in automatically generated image segments and offers statistical guarantees in form of asymptotic confidence intervals. To this end, we… ▽ More Quantifying the number of molecules from fluorescence microscopy measurements is an important topic in cell biology and medical research. In this work, we present a consecutive algorithm for super-resolution (STED) scanning microscopy that provides molecule counts in automatically generated image segments and offers statistical guarantees in form of asymptotic confidence intervals. To this end, we first apply a multiscale scanning procedure on STED microscopy measurements of the sample to obtain a system of significant regions, each of which contains at least one molecule with prescribed uniform probability. This system of regions will typically be highly redundant and consists of rectangular building blocks. To choose an informative but non-redundant subset of more naturally shaped regions, we hybridize our system with the result of a generic segmentation algorithm. The diameter of the segments can be of the order of the resolution of the microscope. Using multiple photon coincidence measurements of the same sample in confocal mode, we are then able to estimate the brightness and number of the molecules and give uniform confidence intervals on the molecule counts for each previously constructed segment. In other words, we establish a so-called molecular map with uniform error control. The performance of the algorithm is investigated on simulated and real data. △ Less

Submitted 2 October, 2023; v1 submitted 27 July, 2022; originally announced July 2022.

A multi-objective sustainable planning for a real hazardous waste production problem

Authors: Abed Zabihian-Bisheh, Hadi Rezaei Vandchali, Vahid Kayvanfar, Frank Werner

Abstract: A significant amount of hazardous waste generated from health sectors and industrial processes has posed a major threat to human health by causing environmental issues and contamination of air, soil, and water resources. This paper presents a multi-objective mixed-integer nonlinear programming (MINLP) formulation for a sustainable hazardous waste location-routing problem. The location of the facil… ▽ More A significant amount of hazardous waste generated from health sectors and industrial processes has posed a major threat to human health by causing environmental issues and contamination of air, soil, and water resources. This paper presents a multi-objective mixed-integer nonlinear programming (MINLP) formulation for a sustainable hazardous waste location-routing problem. The location of the facilities and routing decisions for transporting hazardous waste and the waste residue is considered to design a suitable waste collection system. The presented model simultaneously minimizes the total costs of the waste management system, total risks from transportation and facilities, along with CO2 emissions. A real-world case study is presented to illustrate the applicability of the proposed model. To illustrate the significance of sustainability, the results of the original model are compared with the results of the model without considering sustainability. It indicates that, under the condition when sustainability is not taken into account, total cost, transportation, and site risk along with CO2 emission increased, which in turn demonstrated the importance of sustainability. Furthermore, the managerial insights gained from the optimal results would enable the managers to make better decisions in the hazardous waste management system. △ Less

Submitted 3 July, 2022; originally announced July 2022.

Comments: 25 pages, 6 figurs, 12 tables

MSC Class: 90B50; 90B30; 90C11 ACM Class: G.m

TUM autonomous motorsport: An autonomous racing software for the Indy Autonomous Challenge

Authors: Johannes Betz, Tobias Betz, Felix Fent, Maximilian Geisslinger, Alexander Heilmeier, Leonhard Hermansdorfer, Thomas Herrmann, Sebastian Huch, Phillip Karle, Markus Lienkamp, Boris Lohmann, Felix Nobis, Levent Ögretmen, Matthias Rowold, Florian Sauerbeck, Tim Stahl, Rainer Trauth, Frederik Werner, Alexander Wischnewski

Abstract: For decades, motorsport has been an incubator for innovations in the automotive sector and brought forth systems like disk brakes or rearview mirrors. Autonomous racing series such as Roborace, F1Tenth, or the Indy Autonomous Challenge (IAC) are envisioned as playing a similar role within the autonomous vehicle sector, serving as a proving ground for new technology at the limits of the autonomous… ▽ More For decades, motorsport has been an incubator for innovations in the automotive sector and brought forth systems like disk brakes or rearview mirrors. Autonomous racing series such as Roborace, F1Tenth, or the Indy Autonomous Challenge (IAC) are envisioned as playing a similar role within the autonomous vehicle sector, serving as a proving ground for new technology at the limits of the autonomous systems capabilities. This paper outlines the software stack and approach of the TUM Autonomous Motorsport team for their participation in the Indy Autonomous Challenge, which holds two competitions: A single-vehicle competition on the Indianapolis Motor Speedway and a passing competition at the Las Vegas Motor Speedway. Nine university teams used an identical vehicle platform: A modified Indy Lights chassis equipped with sensors, a computing platform, and actuators. All the teams developed different algorithms for object detection, localization, planning, prediction, and control of the race cars. The team from TUM placed first in Indianapolis and secured second place in Las Vegas. During the final of the passing competition, the TUM team reached speeds and accelerations close to the limit of the vehicle, peaking at around 270 km/h and 28 ms2. This paper will present details of the vehicle hardware platform, the developed algorithms, and the workflow to test and enhance the software applied during the two-year project. We derive deep insights into the autonomous vehicle's behavior at high speed and high acceleration by providing a detailed competition analysis. Based on this, we deduce a list of lessons learned and provide insights on promising areas of future work based on the real-world evaluation of the displayed concepts. △ Less

Submitted 13 January, 2023; v1 submitted 31 May, 2022; originally announced May 2022.

Comments: 37 pages, 18 figures, 2 tables

Journal ref: Journal of Field Robotics, 2023, 1-27

Gamma/Hadron Separation with the HAWC Observatory

Authors: R. Alfaro, C. Alvarez, J. D. Álvarez, J. R. Angeles Camacho, J. C. Arteaga-Velázquez, D. Avila Rojas, H. A. Ayala Solares, R. Babu, E. Belmont-Moreno, C. Brisbois, K. S. Caballero-Mora, T. Capistrán, A. Carramiñana, S. Casanova, O. Chaparro-Amaro, U. Cotti, J. Cotzomi, S. Coutiño de León, E. De la Fuente, C. de León, R. Diaz Hernandez, B. L. Dingus, M. A. DuVernois, M. Durocher, J. C. Díaz-Vélez , et al. (68 additional authors not shown)

Abstract: The High Altitude Water Cherenkov (HAWC) gamma-ray observatory observes atmospheric showers produced by incident gamma rays and cosmic rays with energy from 300 GeV to more than 100 TeV. A crucial phase in analyzing gamma-ray sources using ground-based gamma-ray detectors like HAWC is to identify the showers produced by gamma rays or hadrons. The HAWC observatory records roughly 25,000 events per… ▽ More The High Altitude Water Cherenkov (HAWC) gamma-ray observatory observes atmospheric showers produced by incident gamma rays and cosmic rays with energy from 300 GeV to more than 100 TeV. A crucial phase in analyzing gamma-ray sources using ground-based gamma-ray detectors like HAWC is to identify the showers produced by gamma rays or hadrons. The HAWC observatory records roughly 25,000 events per second, with hadrons representing the vast majority ($>99.9\%$) of these events. The standard gamma/hadron separation technique in HAWC uses a simple rectangular cut involving only two parameters. This work describes the implementation of more sophisticated gamma/hadron separation techniques, via machine learning methods (boosted decision trees and neural networks), and summarizes the resulting improvements in gamma/hadron separation obtained in HAWC. △ Less

Submitted 24 May, 2022; originally announced May 2022.

Comments: 35 pages, 9 figures, published in Nuclear Instruments and Methods in Physics Research Section A

Simultaneous Imaging of Widely Differing Particle Concentrations in MPI: Problem Statement and Algorithmic Proposal for Improvement

Authors: Marija Boberg, Nadine Gdaniec, Patryk Szwargulski, Franziska Werner, Martin Möddel, Tobias Knopp

Abstract: Magnetic Particle Imaging (MPI) is a tomographic imaging technique for determining the spatial distribution of superparamagnetic nanoparticles. Current MPI systems are capable of imaging iron masses over a wide dynamic range of more than four orders of magnitude. In theory, this range could be further increased using adaptive amplifiers, which prevent signal clipping. While this applies to a singl… ▽ More Magnetic Particle Imaging (MPI) is a tomographic imaging technique for determining the spatial distribution of superparamagnetic nanoparticles. Current MPI systems are capable of imaging iron masses over a wide dynamic range of more than four orders of magnitude. In theory, this range could be further increased using adaptive amplifiers, which prevent signal clipping. While this applies to a single sample, the dynamic range is severely limited if several samples with different concentrations or strongly inhomogeneous particle distributions are considered. One scenario that occurs quite frequently in pre-clinical applications is that a highly concentrated tracer bolus in the vascular system "shadows" nearby organs with lower effective tracer concentrations. The root cause of the problem is the ill-posedness of the MPI imaging operator, which requires regularization for stable reconstruction. In this work, we introduce a simple two-step algorithm that increases the dynamic range by a factor of four. Furthermore, the algorithm enables spatially adaptive regularization, i.e. highly concentrated signals can be reconstructed with maximum spatial resolution, while low concentrated signals are strongly regularized to prevent noise amplification. △ Less

Submitted 3 May, 2022; originally announced May 2022.

Journal ref: Physics in Medicine & Biology 66(9), 095004 (2021)

Home health care planning with considering flexible starting/ending points and service features

Authors: Pouria Khodabandeh, Vahid Kayvanfar, Majid Rafiee, Frank Werner

Abstract: One of the recently proposed strategies in health systems is providing services to patients at home, improving the service quality, besides reducing the health system costs. In the real world, some services, such as biological tests or blood sampling, force the nurses to start or end his/her route from/at the laboratory instead of the depot, changing the whole optimal planning. The effect of these… ▽ More One of the recently proposed strategies in health systems is providing services to patients at home, improving the service quality, besides reducing the health system costs. In the real world, some services, such as biological tests or blood sampling, force the nurses to start or end his/her route from/at the laboratory instead of the depot, changing the whole optimal planning. The effect of these special service requirements and features has not been considered so far. In this study, a new mathematical model is suggested considering the flexibility of starting/ending places of each nurse's route according to the specific characteristics of each service. Then several sets of problems in various sizes are solved using the proposed model, where the results confirm the efficiency of the proposed approach. In addition, some sensitivity analyses are performed on the parameters of the required features of the services, followed by some managerial insights and directions for future studies. △ Less

Submitted 2 May, 2022; originally announced May 2022.

Comments: 17 pages, 2 figures, 11 tables

MSC Class: 90B35

Scheduling a single machine with compressible jobs to minimize maximum lateness

Authors: Nodari Vakhania, Frank Werner, Alejandro Reynoso

Abstract: The problem of scheduling non-simultaneously released jobs with due dates on a single machine with the objective to minimize the maximum job lateness is known to be strongly NP-hard. Here we consider an extended model in which the compression of the job processing times is allowed. The compression is accomplished at the cost of involving additional emerging resources, whose use, however, yields so… ▽ More The problem of scheduling non-simultaneously released jobs with due dates on a single machine with the objective to minimize the maximum job lateness is known to be strongly NP-hard. Here we consider an extended model in which the compression of the job processing times is allowed. The compression is accomplished at the cost of involving additional emerging resources, whose use, however, yields some cost. With a given upper limit $U$ on the total allowable cost, one wishes to minimize the maximum job lateness. It is clear that, by using the available resources, some jobs may complete earlier and the objective function value may respectively be decreased. As we show here, for minimizing the maximum job lateness, by shortening the processing time of some specially determined jobs, the objective value can be decreased. Although the generalized problem is harder than the generic non-compressible version, given a ``sufficient amount'' of additional resources, we can solve the problem optimally. We determine the compression rate for some specific jobs and develop an algorithm that obtains an optimal solution. Such an approach can be beneficial in practice since the manufacturer can be provided with an information about the required amount of additional resources in order to solve the problem optimally. In case the amount of the available additional resources is less than used in the above solution, i.e., it is not feasible, it is transformed to a tight minimal feasible solution. △ Less

Submitted 14 June, 2023; v1 submitted 18 March, 2022; originally announced March 2022.

Comments: 20 pages

MSC Class: 90B35 ACM Class: G.2

Time-resolved hadronic particle acceleration in the recurrent Nova RS Ophiuchi

Authors: H. E. S. S. Collaboration, F. Aharonian, F. Ait Benkhali, E. O. Angüner, H. Ashkar, M. Backes, V. Baghmanyan, V. Barbosa Martins, R. Batzofin, Y. Becherini, D. Berge, K. Bernlöhr, B. Bi, M. Böttcher, C. Boisson, J. Bolmont, M. de Bony de Lavergne, M. Breuhaus, R. Brose, F. Brun, S. Caroff, S. Casanova, M. Cerruti, T. Chand, A. Chen , et al. (150 additional authors not shown)

Abstract: Recurrent Novae are repeating thermonuclear explosions in the outer layers of white dwarfs, due to the accretion of fresh material from a binary companion. The shock generated by ejected material slamming into the companion star's wind, accelerates particles to very-high-energies. We report very-high-energy (VHE, $\gtrsim100$\,GeV) gamma rays from the recurrent nova RS\,Ophiuchi up to a month afte… ▽ More Recurrent Novae are repeating thermonuclear explosions in the outer layers of white dwarfs, due to the accretion of fresh material from a binary companion. The shock generated by ejected material slamming into the companion star's wind, accelerates particles to very-high-energies. We report very-high-energy (VHE, $\gtrsim100$\,GeV) gamma rays from the recurrent nova RS\,Ophiuchi up to a month after its 2021 outburst, using the High Energy Stereoscopic System. The VHE emission has a similar temporal profile to lower-energy GeV emission, indicating a common origin, with a two-day delay in peak flux. These observations constrain models of time-dependent particle energization, favouring a hadronic emission scenario over the leptonic alternative. This confirms that shocks in dense winds provide favourable environments for efficient cosmic-ray acceleration to very-high-energies. △ Less

Submitted 28 March, 2022; v1 submitted 16 February, 2022; originally announced February 2022.

Comments: Submitted 1st Nov. 2021, first release 10th March 2022 (accepted version)

Self-Adjointness of Toeplitz Operators on the Segal-Bargmann Space

Authors: Wolfram Bauer, Lauritz van Luijk, Alexander Stottmeister, Reinhard F. Werner

Abstract: We prove a new criterion that guarantees self-adjointness of Toeplitz operator with unbounded operator-valued symbols. Our criterion applies, in particular, to symbols with Lipschitz continuous derivatives, which is the natural class of Hamiltonian functions for classical mechanics. For this we extend the Berger-Coburn estimate to the case of vector-valued Segal-Bargmann spaces. Finally, we apply… ▽ More We prove a new criterion that guarantees self-adjointness of Toeplitz operator with unbounded operator-valued symbols. Our criterion applies, in particular, to symbols with Lipschitz continuous derivatives, which is the natural class of Hamiltonian functions for classical mechanics. For this we extend the Berger-Coburn estimate to the case of vector-valued Segal-Bargmann spaces. Finally, we apply our result to prove self-adjointness for a class of (operator-valued) quadratic forms on the space of Schwartz functions in the Schrödinger representation. △ Less

Submitted 5 October, 2022; v1 submitted 9 February, 2022; originally announced February 2022.

Comments: 17 pages

Journal ref: Journal of Functional Analysis, Volume 284, Issue 4, 2023

Evidence for gamma-ray emission from the remnant of Kepler's supernova based on deep H.E.S.S. observations

Authors: H. E. S. S. Collaboration, F. Aharonian, F. Ait Benkhali, E. O. Anguner, H. Ashkar, M. Backes, V. Barbosa Martins, R. Batzofin, Y. Becherini, D. Berge, K. Bernloehr, M. Boettcher, C. Boisson, J. Bolmont, M. de Bony de Lavergne, M. Breuhaus, R. Brose, F. Brun, T. Bulik, T. Bylund, F. Cangemi, S. Caroff, S. Casanova, M. Cerruti, T. Chand , et al. (136 additional authors not shown)

Abstract: Observations with imaging atmospheric Cherenkov telescopes (IACTs) have enhanced our knowledge of nearby supernova (SN) remnants with ages younger than 500 years by establishing Cassiopeia A and the remnant of Tycho's SN as very-high-energy (VHE) gamma-ray sources. The remnant of Kepler's SN, which is the product of the most recent naked-eye supernova in our Galaxy, is comparable in age to the oth… ▽ More Observations with imaging atmospheric Cherenkov telescopes (IACTs) have enhanced our knowledge of nearby supernova (SN) remnants with ages younger than 500 years by establishing Cassiopeia A and the remnant of Tycho's SN as very-high-energy (VHE) gamma-ray sources. The remnant of Kepler's SN, which is the product of the most recent naked-eye supernova in our Galaxy, is comparable in age to the other two, but is significantly more distant. If the gamma-ray luminosities of the remnants of Tycho's and Kepler's SNe are similar, then the latter is expected to be one of the faintest gamma-ray sources within reach of the current generation IACT arrays. Here we report evidence at a statistical level of 4.6 sigma for a VHE signal from the remnant of Kepler's SN based on deep observations by the High Energy Stereoscopic System (H.E.S.S.) with an exposure of 152 hours. The measured integral flux above an energy of 226 GeV is ~0.3% of the flux of the Crab Nebula. The spectral energy distribution (SED) reveals a gamma-ray emitting component connecting the VHE emission observed with H.E.S.S. to the emission observed at GeV energies with Fermi-LAT. The overall SED is similar to that of the remnant of Tycho's SN, possibly indicating the same non-thermal emission processes acting in both these young remnants of thermonuclear SNe. △ Less

Submitted 23 March, 2024; v1 submitted 15 January, 2022; originally announced January 2022.

Comments: Published in Astronomy and Astrophysics. A typo in the normalization constant in Sect. 4 was corrected

Journal ref: A&A 662, A65 (2022)

H.E.S.S. follow-up observations of Binary Black Hole Coalescence events during the second and third Gravitational Waves observing runs of Advanced LIGO and Advanced Virgo

Authors: H. E. S. S. collaboration, :, H. Abdalla, F. Aharonian, F. Ait Benkhali, E. O. Angüner, H. Ashkar, M. Backes, V. Baghmanyan, V. Barbosa Martins, R. Batzofin, Y. Becherini, D. Berge, K. Bernlöhr, B. Bi, M. Böttcher, C. Boisson, J. Bolmont, M. de Bony de Lavergne, R. Brose, F. Brun, T. Bulik, T. Bylund, F. Cangemi, S. Caroff , et al. (129 additional authors not shown)

Abstract: We report on the observations of four well-localized binary black hole (BBH) mergers by the High Energy Stereoscopic System (H.E.S.S.) during the second and third observing runs of Advanced LIGO and Advanced Virgo, O2 and O3. H.E.S.S. can observe $\mathrm{20\,deg^2}$ of the sky at a time and follows up gravitational-wave (GW) events by ``tiling'' localization regions to maximize the covered locali… ▽ More We report on the observations of four well-localized binary black hole (BBH) mergers by the High Energy Stereoscopic System (H.E.S.S.) during the second and third observing runs of Advanced LIGO and Advanced Virgo, O2 and O3. H.E.S.S. can observe $\mathrm{20\,deg^2}$ of the sky at a time and follows up gravitational-wave (GW) events by ``tiling'' localization regions to maximize the covered localization probability. During O2 and O3, H.E.S.S. observed large portions of the localization regions, between 35\% and 75\%, for four BBH mergers (GW170814, GW190512\_180714, GW190728\_064510, and S200224ca). For these four GW events, we find no significant signal from a pointlike source in any of the observations, and set upper limits on the very high energy ($>$100 GeV) $γ$-ray emission. The 1-10 TeV isotropic luminosity of these GW events is below $10^{45}$ erg s$^{-1}$ at the times of the H.E.S.S. observations, around the level of the low-luminosity GRB 190829A. Assuming no changes are made to how follow-up observations are conducted, H.E.S.S. can expect to observe over 60 GW events per year in the fourth GW observing run, O4, of which eight would be observable with minimal latency. △ Less

Submitted 15 December, 2021; originally announced December 2021.

Journal ref: ApJ, 923, 109 (2021)