title: Analyzing cluster randomized trials designed to support generalizable inferences

abstract: Background: When planning a cluster randomized trial, evaluators often have access to an enumerated cohort representing the target population of clusters. Practicalities of conducting the trial, such as the need to oversample clusters with certain characteristics to improve trial economy or to support inference about subgroups of clusters, may preclude simple random sampling from the cohort into the trial, and thus interfere with the goal of producing generalizable inferences about the target population. Methods: We describe a nested trial design where the randomized clusters are embedded within a cohort of trial-eligible clusters from the target population and where clusters are selected for inclusion in the trial with known sampling probabilities that may depend on cluster characteristics (e.g., allowing clusters to be chosen to facilitate trial conduct or to examine hypotheses related to their characteristics). We develop and evaluate methods for analyzing data from this design to generalize causal inferences to the target population underlying the cohort. Results: We present identification and estimation results for the expectation of the average potential outcome and for the average treatment effect, in the entire target population of clusters and in its non-randomized subset. In simulation studies we show that different estimators have low bias but markedly different precision. Conclusions: Cluster randomized trials where clusters are selected for inclusion with known sampling probabilities that depend on cluster characteristics, combined with efficient estimation methods, can precisely quantify treatment effects in the target population, while addressing objectives of trial conduct that require oversampling clusters on the basis of their characteristics.


title: Multi-task Unscented Kalman Inversion for joint inversion of receiver function and surface wave dispersion

abstract: Based on the recently developed theory of Unscented Kalman Inversion in computational mathematics, we proposed a Bayesian joint inversion framework, i.e., Multi-task Unscented Kalman Inversion (MTUKI), and apply it to the joint inversion of receiver function (RF) and surface wave dispersion (SWD). This method can share information between different observations in a derivative-free way and provide an efficient Gaussian approximation to the posterior distribution of model parameters (thickness and S-wave velocity in each layer of media). The theory and experiments show that our proposed framework demonstrates superior performance in terms of robustness, accuracy, and high efficiency.


title: ECLIPSE: Efficient Long-range Video Retrieval using Sight and Sound

abstract: We introduce an audiovisual method for long-range text-to-video retrieval. Unlike previous approaches designed for short video retrieval (e.g., 5-15 seconds in duration), our approach aims to retrieve minute-long videos that capture complex human actions. One challenge of standard video-only approaches is the large computational cost associated with processing hundreds of densely extracted frames from such long videos. To address this issue, we propose to replace parts of the video with compact audio cues that succinctly summarize dynamic audio events and are cheap to process. Our method, named ECLIPSE (Efficient CLIP with Sound Encoding), adapts the popular CLIP model to an audiovisual video setting, by adding a unified audiovisual transformer block that captures complementary cues from the video and audio streams. In addition to being 2.92x faster and 2.34x memory-efficient than long-range video-only approaches, our method also achieves better text-to-video retrieval accuracy on several diverse long-range video datasets such as ActivityNet, QVHighlights, YouCook2, DiDeMo and Charades.


title: Phase transitions in de Sitter: The stochastic formalism

abstract: The stochastic spectral expansion method offers a simple framework for calculations in de Sitter spacetimes. We show how to extend its reach to metastable vacuum states, both in the case when the potential is bounded from below, and when it is unbounded from below and therefore no stable vacuum state exists. In both cases, the decay rate of the metastable vacuum is given by the lowest non-zero eigenvalue associated to the Fokker-Planck equation. We show how the corresponding eigenfunction determines the field probability distribution which can be used to compute correlation functions and other observables in the metastable vacuum state.


title: Designing globally optimal entangling gates using geometric space curves

abstract: High-fidelity entangling gates are essential for quantum computation. Currently, most approaches to designing such gates are based either on simple, analytical pulse waveforms or on ones obtained from numerical optimization techniques. In both cases, it is typically not possible to obtain a global understanding of the space of waveforms that generate a target gate operation, making it challenging to design globally optimal gates. Here, we show that in the case of weakly coupled qubits, it is possible to find all pulses that implement a target entangling gate. We do this by mapping quantum evolution onto geometric space curves. We derive the minimal conditions these curves must satisfy in order to guarantee a gate with a desired entangling power is implemented. Pulse waveforms are extracted from the curvatures of these curves. We illustrate our method by designing fast, CNOT-equivalent entangling gates for silicon quantum dot spin qubits with fidelities exceeding 99%. We show that fidelities can be further improved while maintaining low bandwidth requirements by using geometrically derived pulses as initial guesses in numerical optimization routines.


title: PAnDR: Fast Adaptation to New Environments from Offline Experiences via Decoupling Policy and Environment Representations

abstract: Deep Reinforcement Learning (DRL) has been a promising solution to many complex decision-making problems. Nevertheless, the notorious weakness in generalization among environments prevent widespread application of DRL agents in real-world scenarios. Although advances have been made recently, most prior works assume sufficient online interaction on training environments, which can be costly in practical cases. To this end, we focus on an offline-training-online-adaptation setting, in which the agent first learns from offline experiences collected in environments with different dynamics and then performs online policy adaptation in environments with new dynamics. In this paper, we propose Policy Adaptation with Decoupled Representations (PAnDR) for fast policy adaptation. In offline training phase, the environment representation and policy representation are learned through contrastive learning and policy recovery, respectively. The representations are further refined by mutual information optimization to make them more decoupled and complete. With learned representations, a Policy-Dynamics Value Function (PDVF) [Raileanu et al., 2020] network is trained to approximate the values for different combinations of policies and environments from offline experiences. In online adaptation phase, with the environment context inferred from few experiences collected in new environments, the policy is optimized by gradient ascent with respect to the PDVF. Our experiments show that PAnDR outperforms existing algorithms in several representative policy adaptation problems.


title: Thermal activation of low-density Ga implanted in Ge

abstract: The nuclear spins of low-density implanted Ga atoms in Ge are interesting candidates for solid state-based qubits. To date, activation studies of implanted Ga in Ge have focused on high densities. Here we extend activation studies into the low-density regime. We use spreading resistance profiling and secondary ion mass spectrometry to derive electrical activation of Ga ions implanted into Ge as a function of rapid thermal anneal temperature and implant density. We show that for our implant conditions the activation is best for anneal temperatures between 400 and 650 $^\circ$C, with a maximum activation of 64% at the highest fluence. Below 400 $^\circ$C, remaining implant damage results in defects that act as superfluous carriers, and above 650 $^\circ$C, surface roughening and loss of Ga ions are observed. The activation increased monotonically from 10% to 64% as the implant fluence increased from $6\times10^{10}$ to $6\times10^{12}$ cm$^{-2}$. The results provide thermal anneal conditions to be used for initial studies of using low-density Ga atoms in Ge as nuclear spin qubits.


title: Combinatorics of Integer Partitions With Prescribed Perimeter

abstract: We prove that the number of even parts and the number of times that parts are repeated have the same distribution over integer partitions with a fixed perimeter. This refines Straub's analog of Euler's Odd-Distinct partition theorem. We generalize the two concerned statistics to these of the part-difference less than $d$ and the parts not congruent to $1$ modulo $d+1$ and prove a distribution inequality, that has a similar flavor as Alder's ex-conjecture, over partitions with a prescribed perimeter. Both of our results are proved analytically and combinatorially.


title: Leading superconducting instabilities in three-dimensional models for Sr2RuO4

abstract: The superconductor Sr2RuO4 has been the subject of enormous interest over more than two decades, but until now the form of its order parameter has not been determined. Since groundbreaking NMR experiments revealed that the pairs are of dominant spin-singlet character, attention has focused on time-reversal symmetry breaking linear combinations of $s$-, $d$- and $g$-wave one-dimensional (1D) irreducible representations. However, a state of the form $d_{xz}+id_{yz}$ has also been proposed. We present a systematic study of the stability of various superconducting candidate states, assuming that pairing is driven by the fluctuation exchange mechanism, including a realistic three-dimensional Fermi surface, full treatment of both local and non-local spin-orbit couplings, and a wide range of interaction parameters $U,J,U',J'$. The leading superconducting instabilities are found to exhibit nodal even-parity $A_{1g} (s')$ or $B_{1g} (d_{x^2-y^2})$ symmetries, similar to the findings in two-dimensional models without longer-range Coulomb interaction which tends to favor $d_{xy}$ over $d_{x^2-y^2}$. Within the so-called Hund's coupling mean-field pairing scenario, the $E_g (d_{xz}/d_{yz})$ solution can be stabilized for large $J$ and specific forms of the spin-orbit coupling, but for all cases studied here the eigenvalues of other superconducting solutions are significantly larger when the full fluctuation exchange vertex is included in the pairing kernel. Additionally, we compute the spin susceptibility in relevant superconducting candidate phases and compare to recent neutron scattering and NMR Knight shift measurements.


title: Community Action on FAIR Data will Fuel a Revolution in Materials Research

abstract: Data - arguably the most important product of worldwide materials research investment - are rarely shared. The small and biased proportion of results published are buried in plots and text licensed by journals. This situation wastes resources, hinders innovation, and, in the current era of data-driven discovery, is no longer tenable. In this comment, we identify opportunities for synergistic, collaborative, and global actions to assemble large quantities of FAIR (Findable, Accessible, Interoperable, Reusable) (1) materials data.


title: Meager composants of tree-like continua

abstract: A subset $M$ of a continuum $X$ is called a \textit{meager composant} if $M$ is maximal with respect to the property that every two of its points are contained in a nowhere dense subcontinuum of $X$. Motivated by questions of Bellamy, Mouron and Ordo\~{n}ez, we show that no tree-like continuum has a proper open meager composant, and that every tree-like continuum has either $1$ or $2^{\aleph_0}$ meager composants. We also prove a decomposition theorem: If $X$ is tree-like and every indecomposable subcontinuum of $X$ is nowhere dense, then the partition of $X$ into meager composants is upper semi-continuous and the space of meager composants is a dendrite.


title: A Convex Optimization Approach for Control of Linear Quadratic Systems with Multiplicative Noise via System Level Synthesis

abstract: This paper presents a convex optimization-based solution to the design of state-feedback controllers for solving the linear quadratic regulator (LQR) problem of uncertain discrete-time systems with multiplicative noise. To synthesize a tractable solution, the recently developed system level synthesis (SLS) framework is leveraged. It is shown that SLS shifts the controller synthesis task from the design of a robust controller to the design of the entire set-valued closed-loop system responses. To this end, the closed-loop system response is entirely characterized by probabilistic set-valued maps from the additive noise to control actions and states. A bi-level convex optimization over the achievable set-valued closed-loop responses is then developed to optimize the expected value of the LQR cost against the worst-case closed-loop system response. The solution to this robust optimization problem may be too conservative since it aims at enforcing the design constraints for all possible system realizations. To deal with this issue, the presented optimization problem is next reformulated as a chance-constrained program (CCP) in which the guarantees are not intended in a deterministic sense of satisfaction against all possible closed-loop system responses, but are instead intended in a probabilistic sense of satisfaction against all but a small fraction of the system responses. To approximately solve the CCP without the requirement of knowing the probabilistic description of the uncertainty in the system matrices, the so-called scenario optimization approach is employed, which provides probabilistic guarantees based on a finite number of system realizations and results in a convex optimization program with moderate computational complexity. Finally, numerical simulations are presented to illustrate the theoretical findings.


title: Speed limits of the trace distance for open quantum system

abstract: We investigate the speed limit of the state transformation in open quantum systems described by the Lindblad type quantum master equation. We obtain universal bounds of the total entropy production described by the trace distance between the initial and final states in the interaction picture. Our bounds can be tighter than the bound of Vu and Hasegawa [Phys. Rev. Lett. 126, 010601 (2021)] which measures the distance by the eigenvalues of the initial and final states: This distance is less than or equal to the trace distance. For this reason, our results can significantly improve Vu-Hasegawa's bound. The trace distance in the Schr\"{o}dinger picture is bounded by a sum of the trace distance in the interaction picture and the trace distance for unitary dynamics described by only the Hamiltonian in the quantum master equation.


title: Model for two-body collisions between ultracold dipolar molecules around a F\"orster resonance in an electric field

abstract: We propose a one-channel, simple model to describe the dynamics of ultracold dipolar molecules around a F\"orster resonance. Slightly above a specific electric field, a collisional shielding can take place, suppressing the molecular losses in a gas. The overall description of the quantum physical mechanism comes back to the dynamics on a unique energy surface, which depends on the relative distance and angular approach of the molecules. This surface enables to interpret how the dipole moments of the molecules are induced and interlocked by the electric field and the dipole-dipole interaction during the process, especially when the shielding is triggered. Averaging the relative angular motion over a unique partial wave (the lowest one when the ultracold regime is reached), the model reproduces well the behaviour of the rate coefficients observed experimentally and predicted theoretically [Matsuda et al., Science 370, 1324 (2020); Li et al., Nat. Phys. 17, 1144 (2021)]. This economic model encapsulates the main physics of the quantum process. Therefore, it can be used as an alternative to a full quantum dynamical treatment and is promising for future studies of collisions involving more bodies.


title: Towards interferometry of neutrino electromagnetism

abstract: It is predicted within the Standard Model of elementary particles that asymmetric neutrino environments cause rotation of linear polarization of electromagnetic wave -- the birefringence. We demonstrate that this effect is strongly enhanced if additionally the photon is propagating through refractive medium, which effectively increases the photon exposure to the neutrino medium. Our estimate for infrared laser beam in $1\,\mathrm{m}$ long optical fiber exposed to reactor anti-neutrino flux results in linear polarization rotation by the angle $\sim4.6\times10^{-39}\,\mathrm{rad}$. We also derive the proper dependence of the effect on the angle between the directions of photon and neutrino propagation in the laboratory frame. For that purpose we derive the correct form of the basis of polarization four-vectors, which differs from the one widely used in literature. We also estimate the sub-leading optical effect of the neutrino medium due to the neutrino dipole magnetic moment, in terms of a variation of the refractive index and its angular dependence. A rough monochromatic approximation points towards the existence of a resonant enhancement of the effect.


title: Sampling-based Fast Gradient Rescaling Method for Highly Transferable Adversarial Attacks

abstract: Deep neural networks have shown to be very vulnerable to adversarial examples crafted by adding human-imperceptible perturbations to benign inputs. After achieving impressive attack success rates in the white-box setting, more focus is shifted to black-box attacks. In either case, the common gradient-based approaches generally use the $sign$ function to generate perturbations at the end of the process. However, only a few works pay attention to the limitation of the $sign$ function. Deviation between the original gradient and the generated noises may lead to inaccurate gradient update estimation and suboptimal solutions for adversarial transferability, which is crucial for black-box attacks. To address this issue, we propose a Sampling-based Fast Gradient Rescaling Method (S-FGRM) to improve the transferability of the crafted adversarial examples. Specifically, we use data rescaling to substitute the inefficient $sign$ function in gradient-based attacks without extra computational cost. We also propose a Depth First Sampling method to eliminate the fluctuation of rescaling and stabilize the gradient update. Our method can be used in any gradient-based optimizations and is extensible to be integrated with various input transformation or ensemble methods for further improving the adversarial transferability. Extensive experiments on the standard ImageNet dataset show that our S-FGRM could significantly boost the transferability of gradient-based attacks and outperform the state-of-the-art baselines.


title: Mixed superconducting state without applied magnetic field

abstract: A superconducting (SC) mixed state occurs in type-II superconductors where the upper critical field Hc2 is higher than the thermodynamic critical field Hc. When an applied field is in between these fields, the free energy depends weakly on the order parameter which therefore can be small (SC state) or zero (normal state) at different parts of the sample. In this paper we demonstrate how a normal state along a line traversing a superconductor can be turned on and off externally in zero field. The concept is based on a long, current-carrying excitation coil, piercing a ringshaped superconductor. The ring experiences zero field, but the vector potential produced by the coil generates a circular current that destroys superconductivity along a radial line starting at preexisting nucleation points in the sample. Unlike the destruction of superconductivity with magnetic field, the vector potential method is reversible and reproducible; full superconductivity is recovered upon removing the current from the coil, and different cooldowns yield the same normal lines. We suggest potential applications of this magnetic-field-free mixed state.


title: A Cognitive Framework for Delegation Between Error-Prone AI and Human Agents

abstract: With humans interacting with AI-based systems at an increasing rate, it is necessary to ensure the artificial systems are acting in a manner which reflects understanding of the human. In the case of humans and artificial AI agents operating in the same environment, we note the significance of comprehension and response to the actions or capabilities of a human from an agent's perspective, as well as the possibility to delegate decisions either to humans or to agents, depending on who is deemed more suitable at a certain point in time. Such capabilities will ensure an improved responsiveness and utility of the entire human-AI system. To that end, we investigate the use of cognitively inspired models of behavior to predict the behavior of both human and AI agents. The predicted behavior, and associated performance with respect to a certain goal, is used to delegate control between humans and AI agents through the use of an intermediary entity. As we demonstrate, this allows overcoming potential shortcomings of either humans or agents in the pursuit of a goal.


title: DBF: Dynamic Belief Fusion for Combining Multiple Object Detectors

abstract: In this paper, we propose a novel and highly practical score-level fusion approach called dynamic belief fusion (DBF) that directly integrates inference scores of individual detections from multiple object detection methods. To effectively integrate the individual outputs of multiple detectors, the level of ambiguity in each detection score is estimated using a confidence model built on a precision-recall relationship of the corresponding detector. For each detector output, DBF then calculates the probabilities of three hypotheses (target, non-target, and intermediate state (target or non-target)) based on the confidence level of the detection score conditioned on the prior confidence model of individual detectors, which is referred to as basic probability assignment. The probability distributions over three hypotheses of all the detectors are optimally fused via the Dempster's combination rule. Experiments on the ARL, PASCAL VOC 07, and 12 datasets show that the detection accuracy of the DBF is significantly higher than any of the baseline fusion approaches as well as individual detectors used for the fusion.


title: Stochastic volatility modeling of high-frequency CSI 300 index and dynamic jump prediction driven by machine learning

abstract: This paper models stochastic process of price time series of CSI 300 index in Chinese financial market, analyzes volatility characteristics of intraday high-frequency price data. In the new generalized Barndorff-Nielsen and Shephard model, the lag caused by asynchrony of market information is considered, and the problem of lack of long-term dependence is solved. To speed up the valuation process, several machine learning and deep learning algorithms are used to estimate parameter and evaluate forecast results. Tracking historical jumps of different magnitudes offers promising avenues for simulating dynamic price processes and predicting future jumps. Numerical results show that the deterministic component of stochastic volatility processes would always be captured over short and longer-term windows. Research finding could be suitable for influence investors and regulators interested in predicting market dynamics based on realized volatility.


title: Sub-Task Decomposition Enables Learning in Sequence to Sequence Tasks

abstract: The field of Natural Language Processing has experienced a dramatic leap in capabilities with the recent introduction of huge Language Models. Despite this success, natural language problems that involve several compounded steps are still practically unlearnable, even by the largest LMs. This complies with experimental failures for end-to-end learning of composite problems that were demonstrated in a variety of domains. An effective mitigation is to introduce intermediate supervision for solving sub-tasks of the compounded problem. Recently, several works have demonstrated high gains by taking a straightforward approach for incorporating intermediate supervision in compounded natural language problems: the sequence-to-sequence LM is fed with an augmented input, in which the decomposed tasks' labels are simply concatenated to the original input. In this paper, we prove a positive learning result that motivates these recent efforts. We show that when concatenating intermediate supervision to the input and training a sequence-to-sequence model on this modified input, unlearnable composite problems can become learnable. We show that this is true for any family of tasks which on the one hand, are unlearnable, and on the other hand, can be decomposed into a polynomial number of simple sub-tasks, each of which depends only on O(1) previous sub-task results. Beyond motivating contemporary empirical efforts for incorporating intermediate supervision in sequence-to-sequence language models, our positive theoretical result is the first of its kind in the landscape of results on the benefits of intermediate supervision for neural-network learning: Until now, all theoretical results on the subject are negative, i.e., show cases where learning is impossible without intermediate supervision, while our result is positive, showing that learning is facilitated in the presence of intermediate supervision.


title: A Potential Based Quantization Procedure of the Damped Oscillator

abstract: Nowadays, two of the most prospering fields of physics are quantum computing and spintronics. In both, the loss of information and dissipation plays a crucial role. In the present work we formulate the quantization of the dissipative oscillator, which aids understanding of the above mentioned, and creates a theoretical frame to overcome these issues in the future. Based on the Lagrangian framework of the damped spring system, the canonically conjugated pairs and the Hamiltonian of the system are obtained, by which the quantization procedure can be started and consistently applied. As a result, the damping quantum wave equation of the dissipative oscillator is deduced, by which an exact damping wave solution of this equation is obtained. Consequently, we arrive at such an irreversible quantum theory by which the quantum losses can be described.


title: Incompressible limit of a compressible Oldroyd-B model

abstract: In this paper, we consider the Cauchy problem for a compressible Oldroyd-B model in three dimensions. Under some smallness assumptions on the initial data, we obtain the global wellposedness of strong solution with uniform regularity. Moreover, we provide a rigorous justification for the link between the compressible model and the corresponding incompressible model via incompressible limits.


title: Eliminating Orbital Selectivity from the Metal-Insulator Transition by Strong Magnetic Fluctuations

abstract: The orbital-selective electronic behavior is one of the most remarkable manifestations of strong electronic correlations in multi-orbital systems. A prominent example is the orbital-selective Mott transition (OSMT), which is characterized by the coexistence of localized electrons in some orbitals, and itinerant electrons in other orbitals. The state-of-the-art theoretical description of the OSMT in two- and three-dimensional systems is based on local non-perturbative approximations to electronic correlations provided by dynamical mean-field theory or slave spin method. In this work we go beyond this local picture and focus on the effect of spatial collective electronic fluctuations on the OSMT. To this aim, we consider a half-filled Hubbard-Kanamori model on a cubic lattice with two orbitals that have different bandwidths. We show that strong magnetic fluctuations that are inherent in this system prevent the OSMT and favor the N\'eel transition that occurs at the same critical temperature for both orbitals.


title: Scalar Induced Gravitational Waves from Warm Inflation

abstract: Stochastic gravitational waves can be induced from the primordial curvature perturbations generated during inflation, through scalar-tensor mode coupling at the second order of cosmological perturbation theory. Here we discuss a model of warm inflation in which large curvature perturbations are generated at the small scales because of inflaton dissipation. These overdense perturbations then collapse at later epoch to form primordial black holes, as was studied in our earlier work (Ref. \cite{Arya:2019wck}), and therefore may also act as a source to the second-order tensor perturbations. In this study, we calculate the spectrum of these secondary gravitational waves from our warm inflationary model. We find that our model leads to a generation of scalar induced gravitational waves (SIGW) over a frequency range ($1-10^6$) Hz. Further, we discuss the detection possibilities of these SIGWs, taking in account the sensitivity of different ongoing and future gravitational wave experiments.


title: Generalized model for steady-state bifurcations without parameters in memristor-based oscillators with lines of equilibria

abstract: We demonstrate how the pitchfork, transcritical and saddle-node bifurcations of steady states observed in dynamical systems with a finite number of isolated equilibrium points occur in systems with lines of equilibria. The exploration is carried out by using the numerical simulation and linear stability analysis applied to a model of a memristor-based oscillator. First, all the discussed bifurcation scenarios are considered in the context of systems including Chua's memristor with a piecewise-smooth characteristic. Then the memristor characteristic is changed to a function that is smooth everywhere. Finally, the action of the memristor forgetting effect is taken into consideration. The presented results are obtained for electronic circuit models, but the considered bifurcation phenomena can be exhibited by systems with a line of equilibria of any nature.


title: End-to-End Instance Edge Detection

abstract: Edge detection has long been an important problem in the field of computer vision. Previous works have explored category-agnostic or category-aware edge detection. In this paper, we explore edge detection in the context of object instances. Although object boundaries could be easily derived from segmentation masks, in practice, instance segmentation models are trained to maximize IoU to the ground-truth mask, which means that segmentation boundaries are not enforced to precisely align with ground-truth edge boundaries. Thus, the task of instance edge detection itself is different and critical. Since precise edge detection requires high resolution feature maps, we design a novel transformer architecture that efficiently combines a FPN and a transformer decoder to enable cross attention on multi-scale high resolution feature maps within a reasonable computation budget. Further, we propose a light weight dense prediction head that is applicable to both instance edge and mask detection. Finally, we use a penalty reduced focal loss to effectively train the model with point supervision on instance edges, which can reduce annotation costs. We demonstrate highly competitive instance edge detection performance compared to state-of-the-art baselines, and also show that the proposed task and loss are complementary to instance segmentation and object detection.


title: Optimal steering of matrix product states and quantum many-body scars

abstract: Ongoing development of quantum simulators allows for a progressively finer degree of control of quantum many-body systems. This motivates the development of efficient approaches to facilitate the control of such systems and enable the preparation of non-trivial quantum states. Here we formulate an approach to control quantum systems based on matrix product states~(MPS). We compare counter-diabatic and leakage minimization approaches to the so-called local steering problem, that consists in finding the best value of the control parameters for generating a unitary evolution of the specific MPS state in a given direction. In order to benchmark the different approaches, we apply them to the generalization of the PXP model known to exhibit coherent quantum dynamics due to quantum many-body scars. We find that the leakage-based approach generally outperforms the counter-diabatic framework and use it to construct a Floquet model with quantum scars. We perform the first steps towards global trajectory optimization and demonstrate entanglement steering capabilities in the generalized PXP model. Finally we apply our leakage minimization approach to construct quantum scars in the periodically driven non-integrable Ising model.

categories: quant-ph cond-mat.quant-gas cond-mat.stat-mech cond-mat.str-el


title: Partial $*$-algebraic quantum groups and Drinfeld doubles of partial compact quantum groups

abstract: We introduce a notion of partial algebraic quantum group. This is an important special case of a weak multiplier Hopf algebra with integrals, as introduced in the work of Van Daele and Wang. At the same time, it generalizes the notion of partial compact quantum group as introduced by De Commer and Timmermann. As an application, we show that the Drinfeld double of a partial compact quantum group can be defined as a partial $*$-algebraic quantum group.


title: Visualizing Multidimensional Linear Programming Problems

abstract: The article proposes an n-dimensional mathematical model of the visual representation of a linear programming problem. This model makes it possible to use artificial neural networks to solve multidimensional linear optimization problems, the feasible region of which is a bounded non-empty set. To visualize the linear programming problem, an objective hyperplane is introduced, the orientation of which is determined by the gradient of the linear objective function: the gradient is the normal to the objective hyperplane. In the case of searching a maximum, the objective hyperplane is positioned in such a way that the value of the objective function at all its points exceeds the value of the objective function at all points of the feasible region, which is a bounded convex polytope. For an arbitrary point of the objective hyperplane, the objective projection onto the polytope is determined: the closer the objective projection point is to the objective hyperplane, the greater the value of the objective function at this point. Based on the objective hyperplane, a finite regular set of points is constructed, called the receptive field. Using objective projections, an image of a polytope is constructed. This image includes the distances from the receptive points to the corresponding points of the polytope surface. Based on the proposed model, parallel algorithms for visualizing a linear programming problem are constructed. An analytical estimation of its scalability is performed. Information about the software implementation and the results of large-scale computational experiments confirming the efficiency of the proposed approaches are presented.


title: Efficient Bayesian Network Structure Learning via Parameterized Local Search on Topological Orderings

abstract: In Bayesian Network Structure Learning (BNSL), one is given a variable set and parent scores for each variable and aims to compute a DAG, called Bayesian network, that maximizes the sum of parent scores, possibly under some structural constraints. Even very restricted special cases of BNSL are computationally hard, and, thus, in practice heuristics such as local search are used. A natural approach for a local search algorithm is a hill climbing strategy, where one replaces a given BNSL solution by a better solution within some pre-defined neighborhood as long as this is possible. We study ordering-based local search, where a solution is described via a topological ordering of the variables. We show that given such a topological ordering, one can compute an optimal DAG whose ordering is within inversion distance $r$ in subexponential FPT time; the parameter $r$ allows to balance between solution quality and running time of the local search algorithm. This running time bound can be achieved for BNSL without structural constraints and for all structural constraints that can be expressed via a sum of weights that are associated with each parent set. We also introduce a related distance called `window inversions distance' and show that the corresponding local search problem can also be solved in subexponential FPT time for the parameter $r$. For two further natural modification operations on the variable orderings, we show that algorithms with an FPT time for $r$ are unlikely. We also outline the limits of ordering-based local search by showing that it cannot be used for common structural constraints on the moralized graph of the network.


title: In Situ Plasma Processing of Superconducting Cavities at Jefferson Lab

abstract: Jefferson Lab began a plasma processing program start-ing in the spring of 2019. Plasma processing is a common technique for removing hydrocarbons from surfaces, which increases the work function and reduces the sec-ondary emission coefficient. Unlike helium processing which relies on ion bombardment of the field emitters, plasma processing uses free oxygen produced in the plasma to break down the hydrocarbons on the surface of the cavity. The residuals of the hydrocarbons in the form of water, carbon monoxide and carbon dioxide are re-moved from the cryomodule as part of the process gas flow. The initial focus of the effort is processing C100 cavities by injecting RF power into the HOM coupler ports. We will then start investigating processing of C50 cavities by introducing RF into the fundamental power coupler. The plan is to start processing cryomodules in the CEBAF tunnel in the mid-term future, with a goal of improving the operational gradients and the energy mar-gin of the linacs. This work will describe the systems and methods used at JLAB for processing cavities using an argon oxygen gas mixture. Before and after plasma pro-cessing results will also be presented.


title: Vecchia-approximated Deep Gaussian Processes for Computer Experiments

abstract: Deep Gaussian processes (DGPs) upgrade ordinary GPs through functional composition, in which intermediate GP layers warp the original inputs, providing flexibility to model non-stationary dynamics. Two DGP regimes have emerged in recent literature. A "big data" regime, prevalent in machine learning, favors approximate, optimization-based inference for fast, high-fidelity prediction. A "small data" regime, preferred for computer surrogate modeling, deploys posterior integration for enhanced uncertainty quantification (UQ). We aim to bridge this gap by expanding the capabilities of Bayesian DGP posterior inference through the incorporation of the Vecchia approximation, allowing linear computational scaling without compromising accuracy or UQ. We are motivated by surrogate modeling of simulation campaigns with upwards of 100,000 runs - a size too large for previous fully-Bayesian implementations - and demonstrate prediction and UQ superior to that of "big data" competitors. All methods are implemented in the "deepgp" package on CRAN.


title: EMMT: A simultaneous eye-tracking, 4-electrode EEG and audio corpus for multi-modal reading and translation scenarios

abstract: We present the Eyetracked Multi-Modal Translation (EMMT) corpus, a dataset containing monocular eye movement recordings, audio and 4-electrode electroencephalogram (EEG) data of 43 participants. The objective was to collect cognitive signals as responses of participants engaged in a number of language intensive tasks involving different text-image stimuli settings when translating from English to Czech. Each participant was exposed to 32 text-image stimuli pairs and asked to (1) read the English sentence, (2) translate it into Czech, (3) consult the image, (4) translate again, either updating or repeating the previous translation. The text stimuli consisted of 200 unique sentences with 616 unique words coupled with 200 unique images as the visual stimuli. The recordings were collected over a two week period and all the participants included in the study were Czech natives with strong English skills. Due to the nature of the tasks involved in the study and the relatively large number of participants involved, the corpus is well suited for research in Translation Process Studies, Cognitive Sciences among other disciplines.


title: Knowledge Base Index Compression via Dimensionality and Precision Reduction

abstract: Recently neural network based approaches to knowledge-intensive NLP tasks, such as question answering, started to rely heavily on the combination of neural retrievers and readers. Retrieval is typically performed over a large textual knowledge base (KB) which requires significant memory and compute resources, especially when scaled up. On HotpotQA we systematically investigate reducing the size of the KB index by means of dimensionality (sparse random projections, PCA, autoencoders) and numerical precision reduction. Our results show that PCA is an easy solution that requires very little data and is only slightly worse than autoencoders, which are less stable. All methods are sensitive to pre- and post-processing and data should always be centered and normalized both before and after dimension reduction. Finally, we show that it is possible to combine PCA with using 1bit per dimension. Overall we achieve (1) 100$\times$ compression with 75%, and (2) 24$\times$ compression with 92% original retrieval performance.


title: Scalar gauge-Higgs models with discrete Abelian symmetry groups

abstract: We investigate the phase diagram and the nature of the phase transitions of three-dimensional lattice gauge-Higgs models obtained by gauging the Z_N subgroup of the global Z_q invariance group of the Z_q clock model (N is a submultiple of q). The phase diagram is generally characterized by the presence of three different phases, separated by three distinct transition lines. We investigate the critical behavior along the two transition lines characterized by the ordering of the scalar field. Along the transition line separating the disordered-confined phase from the ordered-deconfined phase, standard arguments within the Landau-Ginzburg-Wilson framework predict that the behavior is the same as in a generic ferromagnetic model with Z_p global symmetry, p being the ratio q/N. Thus, continuous transitions belong to the Ising and to the O(2) universality class for p=2 and p>3, respectively, while for p=3 only first-order transitions are possible. The results of Monte Carlo simulations confirm these predictions. There is also a second transition line, which separates two phases in which gauge fields are essentially ordered. Along this line we observe the same critical behavior as in the Z_q clock model, as it occurs in the absence of gauge fields.


title: Question Generation for Reading Comprehension Assessment by Modeling How and What to Ask

abstract: Reading is integral to everyday life, and yet learning to read is a struggle for many young learners. During lessons, teachers can use comprehension questions to increase engagement, test reading skills, and improve retention. Historically such questions were written by skilled teachers, but recently language models have been used to generate comprehension questions. However, many existing Question Generation (QG) systems focus on generating literal questions from the text, and have no way to control the type of the generated question. In this paper, we study QG for reading comprehension where inferential questions are critical and extractive techniques cannot be used. We propose a two-step model (HTA-WTA) that takes advantage of previous datasets, and can generate questions for a specific targeted comprehension skill. We propose a new reading comprehension dataset that contains questions annotated with story-based reading comprehension skills (SBRCS), allowing for a more complete reader assessment. Across several experiments, our results show that HTA-WTA outperforms multiple strong baselines on this new dataset. We show that the HTA-WTA model tests for strong SCRS by asking deep inferential questions.


title: A Friendly Tutorial on Mean-Field Spin Glass Techniques for Non-Physicists

abstract: This tutorial is based on lecture notes written for a class taught in the Statistics Department at Stanford in the Winter Quarter of 2017. The objective was to provide a working knowledge of some of the techniques developed over the last 40 years by theoretical physicists and mathematicians to study mean field spin glasses and their applications to high-dimenensional statistics and statistical learning.


title: Shortened universal cycles for permutations

abstract: Kitaev, Potapov, and Vajnovszki [On shortening u-cycles and u-words for permutations, Discrete Appl. Math, 2019] described how to shorten universal words for permutations, to length $n!+n-1-i(n-1)$ for any $i \in [(n-2)!]$, by introducing incomparable elements. They conjectured that it is also possible to use incomparable elements to shorten universal cycles for permutations to length $n!-i(n-1)$ for any $i \in [(n-2)!]$. In this note we prove their conjecture. The proof is constructive, and, on the way, we also show a new method for constructing universal cycles for permutations.


title: First-principles density limit scaling in tokamaks based on edge turbulent transport and implications for ITER

abstract: A first-principles scaling law, based on turbulent transport considerations, and a multi-machine database of density limit discharges from the ASDEX Upgrade, JET and TCV tokamaks, show that the increase of the boundary turbulent transport with the plasma collisionality sets the maximum density achievable in tokamaks. This scaling law shows a strong dependence on the heating power, therefore predicting for ITER a significantly larger safety margin than the Greenwald empirical scaling (Greenwald et al, Nucl. Fusion, 28(12), 1988) in case of unintentional H-L transition.


title: Variational Quantum Evolution Equation Solver

abstract: Variational quantum algorithms offer a promising new paradigm for solving partial differential equations on near-term quantum computers. Here, we propose a variational quantum algorithm for solving a general evolution equation through implicit time-stepping of the Laplacian operator. The use of encoded source states informed by preceding solution vectors results in faster convergence compared to random re-initialization. Through statevector simulations of the heat equation, we demonstrate how the time complexity of our algorithm scales with the ansatz volume for gradient estimation and how the time-to-solution scales with the diffusion parameter. Our proposed algorithm extends economically to higher-order time-stepping schemes, such as the Crank-Nicolson method. We present a semi-implicit scheme for solving systems of evolution equations with non-linear terms, such as the reaction-diffusion and the incompressible Navier-Stokes equations, and demonstrate its validity by proof-of-concept results.


title: Domination ratio of a family of integer distance digraphs with arbitrary degree

abstract: An integer distance digraph is the Cayley graph $\Gamma(\mathbb{Z},S)$ of the additive group $\mathbb{Z}$ of all integers with respect to a finite subset $S\subseteq\mathbb{Z}$. The domination ratio of $\Gamma(\mathbb{Z},S)$, defined as the minimum density of its dominating sets, is related to some number theory problems, such as tiling the integers and finding the maximum density of a set of integers with missing differences. We precisely determine the domination ratio of the integer distance graph $\Gamma(\mathbb{Z},\{1,2,\ldots,d-2,s\})$ for any integers $d$ and $s$ satisfying $d\ge2$ and $s\notin[0,d-2]$. Our result generalizes a previous result on the domination ratio of the graph $\Gamma(\mathbb{Z},\{1,s\})$ with $s\in\mathbb{Z}\setminus\{0,1\}$ and also implies the domination number of certain circulant graphs $\Gamma(\mathbb{Z}_n,S)$, where $\mathbb{Z}_n$ is the finite cyclic group of integers modulo $n$ and $S$ is a subset of $\mathbb{Z}_n$.


title: Some perspectives on Eisenstein series

abstract: We review some topics in the analytic theory of Eisenstein series, including meromorphic continuation, $L^2$-spectral expansion and Fourier coefficients. We also discuss some open problems.


title: Compact Post-Quantum Signatures from Proofs of Knowledge leveraging Structure for the PKP, SD and RSD Problems

abstract: The MPC-in-the-head introduced in [IKOS07] has established itself as an important paradigm to design efficient digital signatures. It has been leveraged in the Picnic scheme [CDG+ 20] that reached the third round of the NIST PQC Standardization process. It has also been used in [Beu20] to introduce the Proof of Knowledge (PoK) with Helper paradigm. This construction permits to design shorter signatures but induces a non negligible performance overhead as it uses cut-and-choose. In this paper, we introduce the PoK leveraging structure paradigm along with its associated challenge space amplification technique. Our new approach to design PoK brings some improvements over the PoK with Helper one. Indeed, we show how one can substitute the Helper in these constructions by leveraging the underlying structure of the considered problem. This approach does not suffer from the performance overhead inherent to the PoK with Helper paradigm hence offers different trade-offs between security, signature sizes and performances. We also present four new post-quantum signature schemes. The first one is based on a new PoK with Helper for the Syndrome Decoding problem. It relies on ideas from [BGKM22] and [FJR21] and improve the latter using a new technique that can be seen as performing some cut-and-choose with a meet in the middle approach. The three other signatures are based on our new PoK leveraging structure approach and as such illustrate its versatility. We provide new PoK related to the Permuted Kernel Problem (PKP), Syndrome Decoding (SD) problem and Rank Syndrome Decoding (RSD) problem. In practice, these PoK lead to comparable or shorter signatures than existing ones. Indeed, considering (public key + signature), we get sizes below 9kB for our signature related to the PKP problem, below 15kB for our signature related to the SD problem and below 7kB for our signature related to the RSD problem.


title: Simulations of un-equal binary collisions embedded in a turbulent gas cloud far away of a massive gravitational center

abstract: We simulate the collapse of a turbulent gas cloud with the particle-based code Gadget2. We choose two sub-clouds, formed by those particles located around the centers $\vec{r}_L$ and $\vec{r}_R$ and within the radii $r_L$ and $r_R$, respectively. A translational velocity $\vec{v}_L$ or $\vec{v}_R$ is added, so that the sub-clouds move towards each other to collide. The radius and pre-collision velocity of the sub-clouds are chosen to be un-equal, and both head-on and oblique collisions are considered. The simulations are all calibrated to have the same total mass and the initial energy ratio $\alpha=0.16$, which is defined as the ratio of thermal energy to gravitational energy. We compare low-$\beta$ models to a high-$\beta$ models, where $\beta$ is defined as the ratio of kinetic energy to gravitational energy. Finally, we consider the turbulent cloud to be under the gravitational influence of an object located far enough, in order to approximate the tidal effects by means of an azimuthal velocity $V_{\rm cir}$ added to the cloud particles in addition to the translational and turbulent velocities mentioned above. We compare a low-$V_{\rm cir}$ model with a high-$V_{\rm cir}$ one.


title: Gabor frame bound optimizations

abstract: We study sharp frame bounds of Gabor systems over rectangular lattices for different windows and integer oversampling rate. In some cases we obtain optimality results for the square lattice, while in other cases the lattices optimizing the frame bounds and the condition number are rectangular lattices which are different for the respective quantities. Also, in some cases optimal lattices do not exist at all and a degenerated system is optimal.


title: How SVC enables Distributed Caching in MEC?

abstract: With an ever increasing demand for the delivery of internet video service, the service providers are facing a huge challenge to deliver ultra-HD (2k/4k) video at sub-second latency. The multi-access edge computing (MEC) platform actually helps in achieving this objective by caching popular contents at the edge of cellular network. This not only reduces the delivery latency, but also the load and the cost of the backhaul links. However, MEC platforms are afflicted by constrained resources in terms of storage and processing capabilities; and centralized caching of contents may nullify the advantage of reduced latency by lowering the offloading probability. Distributed caching at the edge not only improves the offloading probability, but also dynamically adjusts the load distribution among the MEC servers. In this article, we propose an architecture for deployment of MEC platforms by exploiting the characteristics of a scalable video encoding technique. The layered video coding techniques, such as the scalable video coding (SVC), is used by the content providers to adjust to the network dynamics, by dynamically dropping packets in order to reduce latency. We show how an SVC video easily lends itself to distributed caching at the edge. Then we investigate the latency-storage trade-off by storing the video layers at different parts of the access networks.


title: Branch Decomposition-Independent Edit Distances for Merge Trees

abstract: Edit distances between merge trees of scalar fields have many applications in scientific visualization, such as ensemble analysis, feature tracking or symmetry detection. In this paper, we propose branch mappings, a novel approach to the construction of edit mappings for merge trees. Classic edit mappings match nodes or edges of two trees onto each other, and therefore have to either rely on branch decompositions of both trees or have to use auxiliary node properties to determine a matching. In contrast, branch mappings employ branch properties instead of node similarity information, and are independent of predetermined branch decompositions. Especially for topological features, which are typically based on branch properties, this allows a more intuitive distance measure which is also less susceptible to instabilities from small-scale perturbations. We describe a quartic runtime algorithm for computing optimal branch mappings, which is faster than the only other branch decomposition-independent method in the literature by more than a linear factor. Furthermore, we compare the results of our method on synthetic and real-world examples to demonstrate its practicality and utility.


title: Continuous Variable Entanglement in an Optical Parametric Oscillator Based on a Nondegenerate Four Wave Mixing Process in Hot Alkali Atoms

abstract: We present the measurement of entanglement between twin beams generated with a doubly resonant optical parameter oscillator (OPO) based on four-wave mixing in hot $^{85}$Rb vapor above threshold. This is the first measurement of entanglement in an OPO with $\chi^{(3)}$ media above threshold. We reconstruct the covariance matrix for several configurations and with a full picture of the four side band mode state, we study entanglement between all possible bi-partitions. We show a robust generation of entanglement with stronger generation for a specific pair of modes. For this system, we show that atomic density is a determinant factor for generation and loss of quantum correlations. The generation of entangled fields by an atomic OPO close to atomic resonance of alkali atoms enables natural integration into quantum networks.


title: A survey on recently proposed activation functions for Deep Learning

abstract: Artificial neural networks (ANN), typically referred to as neural networks, are a class of Machine Learning algorithms and have achieved widespread success, having been inspired by the biological structure of the human brain. Neural networks are inherently powerful due to their ability to learn complex function approximations from data. This generalization ability has been able to impact multidisciplinary areas involving image recognition, speech recognition, natural language processing, and others. Activation functions are a crucial sub-component of neural networks. They define the output of a node in the network given a set of inputs. This survey discusses the main concepts of activation functions in neural networks, including; a brief introduction to deep neural networks, a summary of what are activation functions and how they are used in neural networks, their most common properties, the different types of activation functions, some of the challenges, limitations, and alternative solutions faced by activation functions, concluding with the final remarks.


title: Paying More Attention to Self-attention: Improving Pre-trained Language Models via Attention Guiding

abstract: Pre-trained language models (PLM) have demonstrated their effectiveness for a broad range of information retrieval and natural language processing tasks. As the core part of PLM, multi-head self-attention is appealing for its ability to jointly attend to information from different positions. However, researchers have found that PLM always exhibits fixed attention patterns regardless of the input (e.g., excessively paying attention to [CLS] or [SEP]), which we argue might neglect important information in the other positions. In this work, we propose a simple yet effective attention guiding mechanism to improve the performance of PLM by encouraging attention towards the established goals. Specifically, we propose two kinds of attention guiding methods, i.e., map discrimination guiding (MDG) and attention pattern decorrelation guiding (PDG). The former definitely encourages the diversity among multiple self-attention heads to jointly attend to information from different representation subspaces, while the latter encourages self-attention to attend to as many different positions of the input as possible. We conduct experiments with multiple general pre-trained models (i.e., BERT, ALBERT, and Roberta) and domain-specific pre-trained models (i.e., BioBERT, ClinicalBERT, BlueBert, and SciBERT) on three benchmark datasets (i.e., MultiNLI, MedNLI, and Cross-genre-IR). Extensive experimental results demonstrate that our proposed MDG and PDG bring stable performance improvements on all datasets with high efficiency and low cost.


title: A Variational Ansatz for the Ground State of the Quantum Sherrington-Kirkpatrick Model

abstract: We present an ansatz for the ground states of the Quantum Sherrington-Kirkpatrick model, a paradigmatic model for quantum spin glasses. Our ansatz, based on the concept of generalized coherent states, very well captures the fundamental aspects of the model, including the ground state energy and the position of the spin glass phase transition. It further enables us to study some previously unexplored features, such as the non-vanishing longitudinal field regime and the entanglement structure of the ground states. We find that the ground state entanglement can be captured by a simple ensemble of weighted graph states with normally distributed phase gates, leading to a volume law entanglement, contrasting with predictions based on entanglement monogamy.


title: CREIME: A Convolutional Recurrent model for Earthquake Identification and Magnitude Estimation

abstract: The detection and rapid characterisation of earthquake parameters such as magnitude are of prime importance in seismology, particularly in applications such as Earthquake Early Warning (EEW). Traditionally, algorithms such as STA/LTA are used for event detection, while frequency or amplitude domain parameters calculated from 1-3 seconds of first P-arrival data are sometimes used to provide a first estimate of (body wave) magnitude. Owing to extensive involvement of human experts in parameter determination, these approaches are often found to be insufficient. Moreover, these methods are sensitive to the signal to noise ratio and may often lead to false or missed alarms depending on the choice of parameters. We, therefore, propose a multitasking deep learning model the Convolutional Recurrent model for Earthquake Identification and Magnitude Estimation (CREIME) that: (i) detects the first earthquake signal, from background seismic noise, (ii) determines first P arrival time as well as (iii) estimates the magnitude using the raw 3-component waveform data from a single station as model input. Considering, speed is of essence in EEW, we use up to two seconds of P-wave information which, to the best of our knowledge, is a significantly smaller data window (5 second window with up to of P wave data) compared to the previous studies. To examine the robustness of CREIME we test it on two independent datasets and find that it achieves an average accuracy of 98 percent for event vs noise discrimination and is able to estimate first P arrival time and local magnitude with average root mean squared errors of 0.13 seconds and 0.65 units, respectively. We also compare CREIME architecture with architectures of other baseline models, by training them on the same data, and also with traditional algorithms such as STA/LTA, and show that our architecture outperforms these methods.


title: Kinematics of molecular gas in star-forming galaxies with large-scale ionised outflows

abstract: We investigate the kinematics of the molecular gas in a sample of seven edge-on (i>60{\deg}) galaxies identified as hosting large-scale outflows of ionised gas, using ALMA CO(1-0) observations at ~ 1kpc resolution. We build on Hogarth et al. (2021; H21), where we find that molecular gas is more centrally concentrated in galaxies which host winds than in control objects. We perform full 3-dimensional kinematic modelling with multiple combinations of kinematic components, allowing us to infer whether these objects share any similarities in their molecular gas structure. We use modelling to pinpoint the kinematic centre of each galaxy, in order to interpret their minor- and major-axis position velocity diagrams (PVDs). From the PVDs, we find that the bulk of the molecular gas in our galaxies is dynamically cold, tracing the rotation curves predicted by our symmetric, rotation-dominated models, but with minor flux asymmetries. Most notably, we find evidence of radial gas motion in a subset of our objects, which demonstrate a characteristic "twisting" in their minor-axis PVDs generally associated with gas flow along the plane of a galaxy. In our highest S/N object, we include bi-symmetric radial flow in our kinematic model, and find (via the Bayesian Information Criterion) that the presence of radial gas motion is strongly favoured. This may provide one mechanism by which molecular gas and star formation are centrally concentrated, enabling the launch of massive ionised gas winds. However, in the remainder of our sample, we do not observe evidence that gas is being driven radially, once again emphasising the variety of physical processes that may be powering the outflows in these objects, as originally noted in H21.


title: Density Fluctuations in Stochastic Kinematic Flows

abstract: At the macroscopic scale, many important models of collective motion fall into the class of kinematic flows for which both velocity and diffusion terms depend only on particle density. When total particle numbers are fixed and finite, simulations of corresponding microscopic dynamics exhibit stochastic effects which can induce a variety of interesting behaviours not present in the large system limit. In this article we undertake a systematic examination of finite-size fluctuations in a general class of particle models whose statistics correspond to those of stochastic kinematic flows. Doing so, we are able to characterise phenomena including: quasi-jams in models of traffic flow; stochastic pattern formation amongst spatially-coupled oscillators; anomalous bulk sub-diffusion in porous media; and travelling wave fluctuations in a model of bacterial swarming.


title: Six-Functor-Formalisms on Higher Stacks

abstract: In this article, it is shown that derivator six-functor-formalisms on any (classical) site canonically extend to higher geometric stacks as defined by To\"en-Vezzosi under some natural locality conditions. As an application, it is shown that the six-functor-formalisms of Morel-Voevodsky-Ayoub extend to higher (Nisnevich-)Artin stacks locally of finite type over some fixed base scheme.


title: Snowmass2021 - White Paper, Implications of Energy Peak for Collider Phenomenology: Top Quark Mass Determination and Beyond

abstract: We first review the decade-old, broad collider physics research program dubbed energy-peaks. We consider the energy distribution of a massless particle in the lab frame arising from the two-body decay of a heavy particle produced unpolarized, whose boost distribution is arbitrary. Remarkably, the location of the peak of this child particle's energy distribution is identical to its single-valued energy in the rest frame of the parent, which is a function of the parent's mass and that of the other decay product. We summarize generalizations to other types of decay and a variety of applications to BSM. The energy-peak idea can also furnish a measurement of the top quark via the energy of the bottom quark from its decay, which, based on the "parent-boost-invariance," is less sensitive to details of the production mechanism of the top quark (cf.~most other methods assume purely SM production of the top quarks, hence are subject to uncertainties therein, including a possible BSM contribution). The original proposal along this line was to simply use the $b$-jet energy as a very good approximation to the bottom quark energy. This method has been successfully implemented by the CMS collaboration. However, the $b$-jet energy-peak method is afflicted by the jet-energy scale (JES) uncertainty. Fortunately, this drawback can be circumvented by using the decay length of a $B$-hadron contained in the $b$-jet as a proxy for the bottom quark energy. An interesting proposal is to then appropriately dovetail the above two ideas resulting in a "best of both worlds" determination of the top quark mass, i.e., based on a measurement of the $B$-hadron decay length, but improved by the energy-peak concept: this would be free of JES uncertainty and largely independent of the top quark production model. We summarize here the results of such an analysis which is to appear in a forthcoming paper.


title: Beam Search: Faster and Monotonic

abstract: Beam search is a popular satisficing approach to heuristic search problems that allows one to trade increased computation time for lower solution cost by increasing the beam width parameter. We make two contributions to the study of beam search. First, we show how to make beam search monotonic; that is, we provide a new variant that guarantees non-increasing solution cost as the beam width is increased. This makes setting the beam parameter much easier. Second, we show how using distance-to-go estimates can allow beam search to find better solutions more quickly in domains with non-uniform costs. Together, these results improve the practical effectiveness of beam search.


title: Fundamental Physics with ESPRESSO, Constraining a simple parametrisation for varying $\alpha$

abstract: The spectrograph ESPRESSO recently obtained a limit on the variation of the fine-structure constant, $\alpha$, through measurements along the line of sight of a bright quasar with a precision of $1.36$ ppm at $1\sigma$ level. This imposes new constraints on cosmological models with a varying $\alpha$. We assume such a model where the electromagnetic sector is coupled to a scalar field dark energy responsible for the current acceleration of the Universe. We parametrise the variation of $\alpha$ with two extra parameters, one defining the cosmological evolution of the quintessence component and the other fixing the coupling with the electromagnetic field. The objective of this work is to constrain these parameters with both astrophysical and local probes. We also carried out a comparative analysis of how each data probe may constrain our parametrisation. We performed a Bayesian analysis by comparing the predictions of the model with observations. The astrophysical datasets are composed of quasar spectra measurements, including the latest ESPRESSO data point, as well as Planck observations of the cosmic microwave background. We combined these with local results from atomic clocks and the MICROSCOPE experiment. The constraints placed on the quintessence parameter are consistent with a null variation of the field, and are therefore compatible with a $\Lambda$CDM cosmology. The constraints on the coupling to the electromagnetic sector are dominated by the E\"otv\"os parameter local bound. More precise measurements with ESPRESSO will be extremely important to study the cosmological evolution of $\alpha$ as it probes an interval of redshift not accessible to other types of observations. However, for this particular model, current available data favour a null variation of $\alpha$ resulting mostly from the strong MICROSCOPE limits.


title: Observation of giant surface second harmonic generation coupled to nematic orders in the van der Waals antiferromagnet FePS$_3$

abstract: Second harmonic generation has been applied to study lattice, electronic and magnetic proprieties in atomically thin materials. However, inversion symmetry breaking is usually required for the materials to generate a large signal. In this work, we report a giant second-harmonic generation that arises below the N\'eel temperature in few-layer centrosymmetric FePS$_3$. Layer-dependent study indicates the detected signal is from the second-order nonlinearity of the surface. The magnetism-induced surface second-harmonic response is two orders of magnitude larger than those reported in other magnetic systems, with the surface nonlinear susceptibility reaching 0.08--0.13 nm$^2$/V in 2 L--5 L samples. By combing linear dichroism and second harmonic generation experiments, we further confirm the giant second-harmonic generation is coupled to nematic orders formed by the three possible Zigzag antiferromagnetic domains. Our study shows that the surface second-harmonic generation is also a sensitive tool to study antiferromagnetic states in centrosymmetric atomically thin materials.


title: An Empirical Study of End-to-End Temporal Action Detection

abstract: Temporal action detection (TAD) is an important yet challenging task in video understanding. It aims to simultaneously predict the semantic label and the temporal interval of every action instance in an untrimmed video. Rather than end-to-end learning, most existing methods adopt a head-only learning paradigm, where the video encoder is pre-trained for action classification, and only the detection head upon the encoder is optimized for TAD. The effect of end-to-end learning is not systematically evaluated. Besides, there lacks an in-depth study on the efficiency-accuracy trade-off in end-to-end TAD. In this paper, we present an empirical study of end-to-end temporal action detection. We validate the advantage of end-to-end learning over head-only learning and observe up to 11\% performance improvement. Besides, we study the effects of multiple design choices that affect the TAD performance and speed, including detection head, video encoder, and resolution of input videos. Based on the findings, we build a mid-resolution baseline detector, which achieves the state-of-the-art performance of end-to-end methods while running more than 4$\times$ faster. We hope that this paper can serve as a guide for end-to-end learning and inspire future research in this field. Code and models are available at \url{https://github.com/xlliu7/E2E-TAD}.


title: Quantitative differentiation and the medial axis

abstract: We study the medial axis of a set $K$ in Euclidean space (the set of points in space with more than one closest point in $K$) from a "coarse" and "quantitative" perspective. We show that on "most" balls $B(x,r)$ in the complement of $K$, the set of almost-closest points to $x$ in $K$ takes up a small angle as seen from $x$. In other words, most locations and scales in the complement of $K$ "appear" to fall outside the medial axis if one looks with only a certain finite resolution. The word "most" involves a Carleson packing condition, and our bounds are independent of the set $K$.


title: Parallel, Portable Algorithms for Distance-2 Maximal Independent Set and Graph Coarsening

abstract: Given a graph, finding the distance-2 maximal independent set (MIS-2) of the vertices is a problem that is useful in several contexts such as algebraic multigrid coarsening or multilevel graph partitioning. Such multilevel methods rely on finding the independent vertices so they can be used as seeds for aggregation in a multilevel scheme. We present a parallel MIS-2 algorithm to improve performance on modern accelerator hardware. This algorithm is implemented using the Kokkos programming model to enable performance portability. We demonstrate the portability of the algorithm and the performance on a variety of architectures (x86/ARM CPUs and NVIDIA/AMD GPUs). The resulting algorithm is also deterministic, producing an identical result for a given input across all of these platforms. The new MIS-2 implementation outperforms implementations in state of the art libraries like CUSP and ViennaCL by 3-8x while producing similar quality results. We further demonstrate the benefits of this approach by developing parallel graph coarsening scheme for two different use cases. First, we develop an algebraic multigrid (AMG) aggregation scheme using parallel MIS-2 and demonstrate the benefits as opposed to previous approaches used in the MueLu multigrid package in Trilinos. We also describe an approach for implementing a parallel multicolor "cluster" Gauss-Seidel preconditioner using this MIS-2 coarsening, and demonstrate better performance with an efficient, parallel, multicolor Gauss-Seidel algorithm.


title: Hydrogen abstraction reactions in formic and thioformic acid isomers by hydrogen and deuterium atoms. Insights on isomerism and deuteration

abstract: The isomerism of molecules in the interstellar medium and the mechanisms behind it are essential questions in the chemistry of organic molecules in space. In particular, for the simple formic and thioformic acids, the low temperatures found in molecular clouds indicate that cis-trans isomerization in the gas-phase must be impeded. Reactions happening on top of interstellar dust grains may explain the isomer interconversion at low temperatures. We studied the isomerization processes of formic and thioformic acid susceptible to happen on the surface of interstellar dust grains and initiated by H abstraction reactions. Similarly, deuterium enrichment of the acids can occur by the same mechanism. Our objective is to shed light on both topics to increase our understanding of key precursors of organic molecules in space.


title: Dependence of the glass transition and jamming densities on spatial dimension

abstract: We investigate the dynamics of soft sphere liquids through computer simulations for spatial dimensions from $d =3$ to $8$, over a wide range of temperatures and densities. Employing a scaling of density-temperature dependent relaxation times, we precisely identify the density $\phi_0$ which marks the ideal glass transition in the hard sphere limit, and a crossover from sub- to super-Arrhenius temperature dependence. The difference between $\phi_0$ and the athermal jamming density $\phi_J$, small in 3 and 4 dimensions, increases with dimension, with $\phi_0 > \phi_J$ for $d > 4$. We compare our results with recent theoretical calculations.


title: Last Layer Re-Training is Sufficient for Robustness to Spurious Correlations

abstract: Neural network classifiers can largely rely on simple spurious features, such as backgrounds, to make predictions. However, even in these cases, we show that they still often learn core features associated with the desired attributes of the data, contrary to recent findings. Inspired by this insight, we demonstrate that simple last layer retraining can match or outperform state-of-the-art approaches on spurious correlation benchmarks, but with profoundly lower complexity and computational expenses. Moreover, we show that last layer retraining on large ImageNet-trained models can also significantly reduce reliance on background and texture information, improving robustness to covariate shift, after only minutes of training on a single GPU.


title: Mass Spectra and Semileptonic Decays of Doubly Heavy $ \Xi $ and $ \Omega $ Baryons

abstract: In the framework of a non-relativistic quark model, the mass spectra of the ground and excited states of doubly heavy $\Xi$ and $\Omega$ baryons are calculated. We estimate the mass difference between the $ \Omega $ and corresponding $ \Xi $ baryons as $M_{\Omega}-M_{\Xi}\simeq178 $ MeV for all the states containing $ cc,~bc $, or $ bb $ quarks. A simple form of the universal Isgur-Wise function, as the transition form factor between the doubly heavy baryons, is introduced. Working in the close-to-zero recoil limit, we investigate the $ b \rightarrow c $ semileptonic decay widths and branching fractions of the doubly heavy baryons. The obtained results are compared with other theoretical predictions.


title: S-R2F2U-Net: A single-stage model for teeth segmentation

abstract: Precision tooth segmentation is crucial in the oral sector because it provides location information for orthodontic therapy, clinical diagnosis, and surgical treatments. In this paper, we investigate residual, recurrent, and attention networks to segment teeth from panoramic dental images. Based on our findings, we suggest three single-stage models: Single Recurrent R2U-Net (S-R2U-Net), Single Recurrent Filter Double R2U-Net (S-R2F2U-Net), and Single Recurrent Attention Enabled Filter Double (S-R2F2-Attn-U-Net). Particularly, S-R2F2U-Net outperforms state-of-the-art models in terms of accuracy and dice score. A hybrid loss function combining the cross-entropy loss and dice loss is used to train the model. In addition, it reduces around 45% of model parameters compared to the R2U-Net model. Models are trained and evaluated on a benchmark dataset containing 1500 dental panoramic X-ray images. S-R2F2U-Net achieves 97.31% of accuracy and 93.26% of dice score, showing superiority over the state-of-the-art methods. Codes are available at https://github.com/mrinal054/teethSeg_sr2f2u-net.git.


title: Variability of the Sun's luminosity places constraints on the thermal equilibrium of the convection zone

abstract: Luminosity, which is the total amount of radiant energy emitted by an object, is one of the most critical quantities in astrophysics for characterizing stars. Equally important is the temporal evolution of a star's luminosity because of its intimate connection with the stellar energy budget, large-scale convective motion, and heat storage in the stellar interior. Here, we model the solar luminosity by extending a semi-empirical total solar irradiance (TSI) model that uses solar-surface magnetism to reconstruct solar irradiance over the entire 4{\pi} solid angle around the Sun. This model was constrained by comparing its output to the irradiance in the Earth's direction with the measured TSI. Comparing the solar luminosity to the TSI on timescales from days to for cycles 23 and 24, we find poor agreement on short timescales (< solar rotation). On longer timescales, however, we find good agreement between the luminosity model and the TSI, which suggests that the extrapolation of luminosities to multi-cycle timescales based on TSI reconstructions may be possible. We show that the solar luminosity is not constant but varies in phase with the solar cycle. This variation has an amplitude of 0.14% from minimum to maximum for solar cycle 23. Considering the energetics in the solar convection zone, it is therefore obvious that a steady-state input from the radiative zone at the solar minimum level would lead to a gradual reduction in the energy content in the convection zone over multi-century timescales. We show that the luminosity at the base of the convection zone should be approximately 0.032% higher than that at the solar surface during solar minimum to maintain net energy equilibrium through the solar cycle. These different energy-input scenarios place constraints on the long-term evolution of the total solar irradiance and its impact on the solar forcing of climate variability.


title: A Note on Boundedness of Singular Integral Operators on Function Spaces over Local Fields

abstract: In this article we consider the classical singular integral operator over a local field with rough kernels. We study the boundedness of such an operator on different function spaces by relaxing the smoothness condition on kernels.


title: A Design Methodology for Fault-Tolerant Computing using Astrocyte Neural Networks

abstract: We propose a design methodology to facilitate fault tolerance of deep learning models. First, we implement a many-core fault-tolerant neuromorphic hardware design, where neuron and synapse circuitries in each neuromorphic core are enclosed with astrocyte circuitries, the star-shaped glial cells of the brain that facilitate self-repair by restoring the spike firing frequency of a failed neuron using a closed-loop retrograde feedback signal. Next, we introduce astrocytes in a deep learning model to achieve the required degree of tolerance to hardware faults. Finally, we use a system software to partition the astrocyte-enabled model into clusters and implement them on the proposed fault-tolerant neuromorphic design. We evaluate this design methodology using seven deep learning inference models and show that it is both area and power efficient.


title: Intervertebral Disc Labeling With Learning Shape Information, A Look Once Approach

abstract: Accurate and automatic segmentation of intervertebral discs from medical images is a critical task for the assessment of spine-related diseases such as osteoporosis, vertebral fractures, and intervertebral disc herniation. To date, various approaches have been developed in the literature which routinely relies on detecting the discs as the primary step. A disadvantage of many cohort studies is that the localization algorithm also yields false-positive detections. In this study, we aim to alleviate this problem by proposing a novel U-Net-based structure to predict a set of candidates for intervertebral disc locations. In our design, we integrate the image shape information (image gradients) to encourage the model to learn rich and generic geometrical information. This additional signal guides the model to selectively emphasize the contextual representation and suppress the less discriminative features. On the post-processing side, to further decrease the false positive rate, we propose a permutation invariant 'look once' model, which accelerates the candidate recovery procedure. In comparison with previous studies, our proposed approach does not need to perform the selection in an iterative fashion. The proposed method was evaluated on the spine generic public multi-center dataset and demonstrated superior performance compared to previous work. We have provided the implementation code in https://github.com/rezazad68/intervertebral-lookonce


title: "The Pedestrian next to the Lamppost" Adaptive Object Graphs for Better Instantaneous Mapping

abstract: Estimating a semantically segmented bird's-eye-view (BEV) map from a single image has become a popular technique for autonomous control and navigation. However, they show an increase in localization error with distance from the camera. While such an increase in error is entirely expected - localization is harder at distance - much of the drop in performance can be attributed to the cues used by current texture-based models, in particular, they make heavy use of object-ground intersections (such as shadows), which become increasingly sparse and uncertain for distant objects. In this work, we address these shortcomings in BEV-mapping by learning the spatial relationship between objects in a scene. We propose a graph neural network which predicts BEV objects from a monocular image by spatially reasoning about an object within the context of other objects. Our approach sets a new state-of-the-art in BEV estimation from monocular images across three large-scale datasets, including a 50% relative improvement for objects on nuScenes.


title: Probing massive neutrinos with the Minkowski functionals of large-scale structure

abstract: Massive neutrinos suppress the growth of structure under their free-streaming scales. The effect is most prominent on small scales where the widely-used two-point statistics can no longer capture the full information. In this work, we study the signatures massive neutrinos leave on large-scale structure (LSS) as revealed by its morphological properties, which are fully described by $4$ Minkowski functionals (MFs), and quantify the constraints on the summed neutrino mass $M_{\nu}$ from the MFs, by using publicly available N-body simulations. We find the MFs provide important complementary information, and give tighter constraints on $M_{\nu}$ than the power spectrum. Specifically, depending on whether massive neutrinos are included in the density field (the `m' field) or not (the `cb' field), we find the constraint on $M_{\nu}$ from the MFs with a smoothing scale of $R_G=5 h^{-1}$Mpc is $48$ or $4$ times better than that from the power spectrum. When the MFs are combined with the power spectrum, they can improve the constraint on $M_{\nu}$ from the latter by a factor of 63 for the `m' field and 5 for the `cb' field. Notably, when the `m' field is used, the constraint on $M_{\nu}$ from the MFs can reach $0.0177$eV with a volume of $1(h^{-1}\rm Gpc)^3$, while the combination of the MFs and power spectrum can tighten this constraint to be $0.0133$eV, a $4.5\sigma$ significance on detecting the minimum sum of the neutrino masses. For the `m' field, we also find the $\sigma_8$ and $M_{\nu}$ degeneracy is broken with the MFs, leading to stronger constraints on all 6 cosmological parameters considered in this work than the power spectrum.


title: Cancelling microwave crosstalk with fixed-frequency qubits

abstract: Scalable quantum information processing requires that modular gate operations can be executed in parallel. The presence of crosstalk decreases the individual addressability, causing erroneous results during simultaneous operations. For superconducting qubits which operate in the microwave regime, electromagnetic isolation is often limited due to design constraints, leading to signal crosstalk that can deteriorate the quality of simultaneous gate operations. Here, we propose and demonstrate a method based on AC Stark effect for calibrating the microwave signal crosstalk. The method is suitable for processors based on fixed-frequency qubits which are known for high coherence and simple control. The optimal compensation parameters can be reliably identified from a well-defined interference pattern. We implement the method on an array of 7 superconducting qubits, and show its effectiveness in removing the majority of crosstalk errors.


title: Marrying Fairness and Explainability in Supervised Learning

abstract: Machine learning algorithms that aid human decision-making may inadvertently discriminate against certain protected groups. We formalize direct discrimination as a direct causal effect of the protected attributes on the decisions, while induced discrimination as a change in the causal influence of non-protected features associated with the protected attributes. The measurements of marginal direct effect (MDE) and SHapley Additive exPlanations (SHAP) reveal that state-of-the-art fair learning methods can induce discrimination via association or reverse discrimination in synthetic and real-world datasets. To inhibit discrimination in algorithmic systems, we propose to nullify the influence of the protected attribute on the output of the system, while preserving the influence of remaining features. We introduce and study post-processing methods achieving such objectives, finding that they yield relatively high model accuracy, prevent direct discrimination, and diminishes various disparity measures, e.g., demographic disparity.


title: Guaranteed Bounds for Posterior Inference in Universal Probabilistic Programming

abstract: We propose a new method to approximate the posterior distribution of probabilistic programs by means of computing guaranteed bounds. The starting point of our work is an interval-based trace semantics for a recursive, higher-order probabilistic programming language with continuous distributions. Taking the form of (super-/subadditive) measures, these lower/upper bounds are non-stochastic and provably correct: using the semantics, we prove that the actual posterior of a given program is sandwiched between the lower and upper bounds (soundness); moreover the bounds converge to the posterior (completeness). As a practical and sound approximation, we introduce a weight-aware interval type system, which automatically infers interval bounds on not just the return value but also weight of program executions, simultaneously. We have built a tool implementation, called GuBPI, which automatically computes these posterior lower/upper bounds. Our evaluation on examples from the literature shows that the bounds are useful, and can even be used to recognise wrong outputs from stochastic posterior inference procedures.


title: Search for the $Z^\prime$ boson decaying to a right-handed neutrino pair in leptophobic $\mathrm{U(1)}$ models

abstract: The $U(1)$ extensions of the Standard Model contain a heavy neutral gauge boson $Z^\prime$. If leptophobic, the boson can evade the stringent bounds from the dilepton resonance searches. We consider two theoretically well-motivated examples of leptophobic $U(1)$ extensions in which the $Z'$ decays to right-handed neutrinos (RHNs) with substantial branchings. The coexistence of a leptophobic $Z^\prime$ and the RHNs opens up a new possibility of searching for these particles simultaneously through the production of a $Z^\prime$ at the LHC and its decay to a RHN pair. For this decay to occur, the RHNs need to be lighter than the $Z^\prime$. Hence, we study this process in an inverse seesaw setup where the RHNs can be in the TeV range. However, in this case, they have a pseudo-Dirac nature, i.e., a RHN pair would produce only opposite-sign lepton pairs, as opposed to the Majorana-type neutrinos, which can produce both same- and opposite-sign lepton pairs. Hence, the final state we study has a same-flavour opposite-sign lepton pair plus hadronically-decaying boosted $W$ bosons. Our analysis shows that the high luminosity LHC can discover a TeV-scale leptophobic $Z^\prime$ decaying via a RHN pair in a wide range of available parameters. Interestingly, large parameter regions beyond the reach of future dijet-resonance searches can be probed exclusively through our channel.


title: Cosmology from Strong Interactions

abstract: The wealth of theoretical and phenomenological information about Quantum Chromodynamics at short and long distances collected so far in major collider measurements has profound implications in cosmology. We provide a brief discussion on the major implications of the strongly coupled dynamics of quarks and gluons as well as on effects due to their collective motion on the physics of the early universe and in astrophysics.


title: Koopman-based spectral clustering of directed and time-evolving graphs

abstract: While spectral clustering algorithms for undirected graphs are well established and have been successfully applied to unsupervised machine learning problems ranging from image segmentation and genome sequencing to signal processing and social network analysis, clustering directed graphs remains notoriously difficult. Two of the main challenges are that the eigenvalues and eigenvectors of graph Laplacians associated with directed graphs are in general complex-valued and that there is no universally accepted definition of clusters in directed graphs. We first exploit relationships between the graph Laplacian and transfer operators and in particular between clusters in undirected graphs and metastable sets in stochastic dynamical systems and then use a generalization of the notion of metastability to derive clustering algorithms for directed and time-evolving graphs. The resulting clusters can be interpreted as coherent sets, which play an important role in the analysis of transport and mixing processes in fluid flows.


title: Inducing Positive Perspectives with Text Reframing

abstract: Sentiment transfer is one popular example of a text style transfer task, where the goal is to reverse the sentiment polarity of a text. With a sentiment reversal comes also a reversal in meaning. We introduce a different but related task called positive reframing in which we neutralize a negative point of view and generate a more positive perspective for the author without contradicting the original meaning. Our insistence on meaning preservation makes positive reframing a challenging and semantically rich task. To facilitate rapid progress, we introduce a large-scale benchmark, Positive Psychology Frames, with 8,349 sentence pairs and 12,755 structured annotations to explain positive reframing in terms of six theoretically-motivated reframing strategies. Then we evaluate a set of state-of-the-art text style transfer models, and conclude by discussing key challenges and directions for future work.


title: Scheduling to Minimize Age of Information with Multiple Sources

abstract: We consider a G/G/1 queueing system with a single server, where updates arrive from different sources stochastically with possibly different update inter-generation time distributions. The server can transmit/serve at most one update at any time, with potentially different transmission/service times for updates belonging to distinct sources. The age of information (AoI) of any source is a function of the time difference between the departure time of successive updates of that source. Each fully/partially transmitted update incurs a fixed (energy) cost, and the goal of the scheduler is to minimize the linear combination of the sum of the age of information across all sources and the total energy cost. We propose a simple non-preemptive randomized scheduling algorithm that randomly marks arriving updates from a source to be eligible for transmission with a fixed probability and discards them otherwise. Every time the server becomes free, it chooses a source for transmission randomly with another fixed probability and begins to transmit the most recently marked update of the chosen source. Both the respective probabilities are chosen by solving a convex program. The competitive ratio of the proposed algorithm (against a non-preemptive offline optimal algorithm) is shown to be 3 plus the maximum of the ratio of the variance and the mean of the inter-arrival time distribution of sources. For several common distributions such as exponential, uniform and Rayleigh, the competitive ratio is at most 4. For preemptive policies, a G/M/1 system is considered and a non-preemptive policy is shown to have competitive ratio (against a preemptive offline optimal algorithm) at most 5 plus the maximum of the ratio of the variance and the mean of the inter-arrival time distribution of sources.


title: Strongly convergent homogeneous approximations to inhomogeneous Markov jump processes and applications

abstract: The study of time-inhomogeneous Markov jump processes is a traditional topic within probability theory that has recently attracted substantial attention in various applications. However, their flexibility also incurs a substantial mathematical burden which is usually circumvented by using well-known generic distributional approximations or simulations. This article provides a novel approximation method that tailors the dynamics of a time-homogeneous Markov jump process to meet those of its time-inhomogeneous counterpart on an increasingly fine Poisson grid. Strong convergence of the processes in terms of the Skorokhod $J_1$ metric is established, and convergence rates are provided. Under traditional regularity assumptions, distributional convergence is established for unconditional proxies, to the same limit. Special attention is devoted to the case where the target process has one absorbing state and the remaining ones transient, for which the absorption times also converge. Some applications are outlined, such as univariate hazard-rate density estimation, ruin probabilities, and multivariate phase-type density evaluation.