title: Discovery of Nb hydride precipitates in superconducting qubits

abstract: We report the first evidence of the formation of niobium hydrides within niobium films on silicon substrates in superconducting qubits fabricated at Rigetti Computing. We combine complementary techniques including room and cryogenic temperature atomic scale high-resolution and scanning transmission electron microscopy (HR-TEM and STEM), atomic force microscopy (AFM), and the time-of-flight secondary ion mass spectroscopy (TOF-SIMS) to reveal the existence of the niobium hydride precipitates directly in the Rigetti chip areas. Electron diffraction and high-resolution transmission electron microscopy (HR-TEM) analyses are performed at room and cryogenic temperatures (~106 K) on superconducting qubit niobium film areas, and reveal the formation of three types of Nb hydride domains with different crystalline orientations and atomic structures. There is also variation in their size and morphology from small (~5 nm) irregular shape domains within the Nb grains to large (~10-100 nm) Nb grains fully converted to niobium hydride. As niobium hydrides are non-superconducting and can easily change in size and location upon different cooldowns to cryogenic temperatures, our findings highlight a new previously unknown source of decoherence in superconducting qubits, contributing to both quasiparticle and two-level system (TLS) losses, and offering a potential explanation for qubit performance changes upon cooldowns. A pathway to mitigate the formation of the Nb hydrides for superconducting qubit applications is also discussed.


title: Critical points of Laplace eigenfunctions on polygons

abstract: We study the critical points of Laplace eigenfunctions on polygonal domains with a focus on the second Neumann eigenfunction. We show that if each convex quadrilaterals has no second Neumann eigenfunction with an interior critical point, then there exists a convex quadrilateral with an unstable critical point. We also show that each critical point of a second-Neumann eigenfunction on a Lip-1 polygon with no orthogonal sides is an acute vertex.


title: A Sensitivity Matrix Approach Using Two-Stage Optimization for Voltage Regulation of LV Networks with High PV Penetration

abstract: The occurrence of voltage violations are a major deterrent for absorbing more roof-top solar power to smart Low Voltage Distribution Grids (LVDG). Recent studies have focused on decentralized control methods to solve this problem due to the high computational time in performing load flows in centralized control techniques. To address this issue a novel sensitivity matrix is developed to estimate voltages of the network by replacing load flow simulations. In this paper, a Centralized Active, Reactive Power Management System (CARPMS) is proposed to optimally utilize the reactive power capability of smart photo-voltaic inverters with minimal active power curtailment to mitigate the voltage violation problem. The developed sensitivity matrix is able to reduce the time consumed by 48% compared to load flow simulations, enabling near real-time control optimization. Given the large solution space of power systems, a novel two-stage optimization is proposed, where the solution space is narrowed down by a Feasible Region Search (FRS) step, followed by Particle Swarm Optimization (PSO). The performance of the proposed methodology is analyzed in comparison to the load flow method to demonstrate the accuracy and the capability of the optimization algorithm to mitigate voltage violations in near real-time. The deviation of mean voltages of the proposed methodology from load flow method was; 6.5*10^-3 p.u for reactive power control using Q-injection, 1.02*10^-2 p.u for reactive power control using Q-absorption, and 0 p.u for active power curtailment case.


title: Implementation and Optimization of the PTOLEMY Transverse Drift Electromagnetic Filter

abstract: The PTOLEMY transverse drift filter is a new concept to enable precision analysis of the energy spectrum of electrons near the tritium beta-decay endpoint. This paper details the implementation and optimization methods for successful operation of the filter. We present the first demonstrator that produces the required magnetic field properties with an iron return-flux magnet. Two methods for the setting of filter electrode voltages are detailed. The challenges of low-energy electron transport in cases of low field are discussed, such as the growth of the cyclotron radius with decreasing magnetic field, which puts a ceiling on filter performance relative to fixed filter dimensions. Additionally, low pitch angle trajectories are dominated by motion parallel to the magnetic field lines and introduce non-adiabatic conditions and curvature drift. To minimize these effects and maximize electron acceptance into the filter, we present a three-potential-well design to simultaneously drain the parallel and transverse kinetic energies throughout the length of the filter. These optimizations are shown, in simulation, to achieve low-energy electron transport from a 1 T iron core (or 3 T superconducting) starting field with initial kinetic energy of 18.6 keV drained to <10 eV (<1 eV) in about 80 cm. This result for low field operation paves the way for the first demonstrator of the PTOLEMY spectrometer for measurement of electrons near the tritium endpoint to be constructed at the Gran Sasso National Laboratary (LNGS) in Italy.


title: Entanglement and fermionization of two distinguishable fermions in a strict and non strict one-dimensional space

abstract: The fermionization regime and entanglement correlations of two distinguishable harmonically confined fermions interacting via a zero-range potential is addressed. We present two alternative representations of the ground state that we associate with two different types of one-dimensional spaces. These spaces, in turn, induce different correlations between particles and thus require a suitable definition of entanglement. We find that the entanglement of the ground state is strongly conditioned by those one-dimensional space features. We also find that in the strongly attractive regime the relative ground state is a highly localized state leading to maximum entanglement. Our analysis shows that in the strongly repulsive regime the ground state changes smoothly from a superposition of Slater-like states to a finite superposition of Slaters, this lack of accessible states yields to Pauli blocking as a strong signature of fermionization. Our results indicate that entangled states could be obtained in current experiments by reaching the non-interacting regime from the interacting regime. Entangled states could also be obtained when a state is brought from the interacting regime into the strongly repulsive regime by changing the scattering length near the confinement-induced resonance. Finally, we show that the first excited state obtained in the absence of interactions and the third excited fermionized state are maximally entangled.


title: Reachability of weakly nonlinear systems using Carleman linearization

abstract: In this article we introduce a solution method for a special class of nonlinear initial-value problems using set-based propagation techniques. The novelty of the approach is that we employ a particular embedding (Carleman linearization) to leverage recent advances of high-dimensional reachability solvers for linear ordinary differential equations based on the support function. Using a global error bound for the Carleman linearization abstraction, we are able to describe the full set of behaviors of the system for sets of initial conditions and in dense time.


title: Scaling Laws and Critical Properties for $FCC$ and $HCP$ Metals

abstract: The determination of the critical parameters of metals has remained particularly challenging both experimentally, because of the very large temperatures involved, and theoretically, because of the many-body interactions that take place in metals. Moreover, experiments have shown that these systems exhibit an unusually strong asymmetry of their binodal. Recent theoretical work has led to new similarity laws, based on the calculation of the Zeno line and of the underlying Boyle parameters, which provided results for the critical properties of atomic and molecular systems in excellent agreement with experiments. Using the recently developed Expanded Wang-Landau (EWL) simulation method, we evaluate the grand-canonical partition function, over a wide range of conditions, for $11$ $FCC$ and $HCP$ metals ($Ag$, $Al$, $Au$, $Be$, $Cu$, $Ir$, $Ni$, $Pb$, $Pd$, $Pt$ and $Rh$), modeled with a many-body interaction potential. This allows us to calculate the binodal, Zeno line, Boyle parameters and, in turn, obtain the critical properties for these systems. We also propose two scaling laws for the enthalpy and entropy of vaporization, and identify critical exponents of $0.4$ and $1.22$ for these two laws, respectively.


title: A generalized stacked reinforcement learning method for sampled systems

abstract: A common setting of reinforcement learning (RL) is a Markov decision process (MDP) in which the environment is a stochastic discrete-time dynamical system. Whereas MDPs are suitable in such applications as video-games or puzzles, physical systems are time-continuous. A general variant of RL is of digital format, where updates of the value (or cost) and policy are performed at discrete moments in time. The agent-environment loop then amounts to a sampled system, whereby sample-and-hold is a specific case. In this paper, we propose and benchmark two RL methods suitable for sampled systems. Specifically, we hybridize model-predictive control (MPC) with critics learning the optimal Q- and value (or cost-to-go) function. Optimality is analyzed and performance comparison is done in an experimental case study with a mobile robot.


title: Anisotropic nematic fluctuations above the ferroquadrupolar transition in TmVO$_4$

abstract: TmVO$_4$ exhibits ferroquadrupolar order below 2.15 K with a well-isolated non-Kramers ground state doublet, and is a model system to understand Ising nematic order. We present $^{51}$V nuclear magnetic resonance data as a function of field orientation in a single crystal. Although the spectra are well understood in terms of direct dipolar hyperfine couplings, the spin lattice relaxation rate exhibits strong anisotropy that cannot be understood in terms of magnetic fluctuations. We find that the spin lattice relaxation rate scales with the shear elastic constant associated with the ferroquadrupole phase transition, suggesting that quadrupole (nematic) fluctuations dominate the spin lattice relaxation for in-plane fields.


title: Dynamic Network Quantization for Efficient Video Inference

abstract: Deep convolutional networks have recently achieved great success in video recognition, yet their practical realization remains a challenge due to the large amount of computational resources required to achieve robust recognition. Motivated by the effectiveness of quantization for boosting efficiency, in this paper, we propose a dynamic network quantization framework, that selects optimal precision for each frame conditioned on the input for efficient video recognition. Specifically, given a video clip, we train a very lightweight network in parallel with the recognition network, to produce a dynamic policy indicating which numerical precision to be used per frame in recognizing videos. We train both networks effectively using standard backpropagation with a loss to achieve both competitive performance and resource efficiency required for video recognition. Extensive experiments on four challenging diverse benchmark datasets demonstrate that our proposed approach provides significant savings in computation and memory usage while outperforming the existing state-of-the-art methods.


title: Using Neighborhood Context to Improve Information Extraction from Visual Documents Captured on Mobile Phones

abstract: Information Extraction from visual documents enables convenient and intelligent assistance to end users. We present a Neighborhood-based Information Extraction (NIE) approach that uses contextual language models and pays attention to the local neighborhood context in the visual documents to improve information extraction accuracy. We collect two different visual document datasets and show that our approach outperforms the state-of-the-art global context-based IE technique. In fact, NIE outperforms existing approaches in both small and large model sizes. Our on-device implementation of NIE on a mobile platform that generally requires small models showcases NIE's usefulness in practical real-world applications.


title: Multi-Dimensional Multiple Access with Resource Utilization Cost Awareness for Individualized Service Provisioning in 6G

abstract: The increasingly diversified Quality-of-Service (QoS) requirements envisioned for future wireless networks call for more flexible and inclusive multiple access techniques for supporting emerging applications and communication scenarios. To achieve this, we propose a multi-dimensional multiple access (MDMA) protocol to meet individual User Equipment's (UE's) unique QoS demands while utilizing multi-dimensional radio resources cost-effectively. In detail, the proposed scheme consists of two new aspects, i.e., selection of a tailored multiple access mode for each UE while considering the UE-specific radio resource utilization cost; and multi-dimensional radio resource allocation among coexisting UEs under dynamic network conditions. To reduce the UE-specific resource utilization cost, the base station (BS) organizes UEs with disparate multi-domain resource constraints as UE coalition by considering each UE's specific resource availability, perceived quality, and utilization capability. Each UE within a coalition could utilize its preferred radio resources, which leads to low utilization cost while avoiding resource-sharing conflicts with remaining UEs. Furthermore, to meet UE-specific QoS requirements and varying resource conditions at the UE side, the multi-dimensional radio resource allocation among coexisting UEs is formulated as an optimization problem to maximize the summation of cost-aware utility functions of all UEs. A solution to solve this NP-hard problem with low complexity is developed using the successive convex approximation and the Lagrange dual decomposition methods. The effectiveness of our proposed scheme is validated by numerical simulation and performance comparison with state-of-the-art schemes. In particular, the simulation results demonstrate that our proposed scheme outperforms these benchmark schemes by large margins.


title: Predicting Vehicles' Longitudinal Trajectories and Lane Changes on Highway On-Ramps

abstract: Vehicles on highway on-ramps are one of the leading contributors to congestion. In this paper, we propose a prediction framework that predicts the longitudinal trajectories and lane changes (LCs) of vehicles on highway on-ramps and tapers. Specifically, our framework adopts a combination of prediction models that inputs a 4 seconds duration of a trajectory to output a forecast of the longitudinal trajectories and LCs up to 15 seconds ahead. Training and Validation based on next generation simulation (NGSIM) data show that the prediction power of the developed model and its accuracy outperforms a traditional long-short term memory (LSTM) model. Ultimately, the work presented here can alleviate the congestion experienced on on-ramps, improve safety, and guide effective traffic control strategies.


title: Approximation and parameterized algorithms to find balanced connected partitions of graphs

abstract: Partitioning a connected graph into $k$~vertex-disjoint connected subgraphs of similar (or given) orders is a classical problem that has been intensively investigated since late seventies. Given a connected graph $G=(V,E)$ and a weight function $w : V \to \mathbb{Q}_\geq$, a connected $k$-partition of $G$ is a partition of $V$ such that each class induces a connected subgraph. The balanced connected $k$-partition problem consists in finding a connected $k$-partition in which every class has roughly the same weight. To model this concept of balance, one may seek connected $k$-partitions that either maximize the weight of a lightest class $(\text{max-min BCP}_k)$ or minimize the weight of a heaviest class $(\text{min-max BCP}_k)$. Such problems are equivalent when $k=2$, but they are different when $k\geq 3$. In this work, we propose a simple pseudo-polynomial $\frac{k}{2}$-approximation algorithm for $\text{min-max BCP}_k$ which runs in time $\mathcal{O}(W|V||E|)$, where $W = \sum_{v \in V} w(v)$. Based on this algorithm and using a scaling technique, we design a (polynomial) $(\frac{k}{2} +\varepsilon)$-approximation for the same problem with running-time $\mathcal{O}(|V|^3|E|/\varepsilon)$, for any fixed $\varepsilon>0$. Additionally, we propose a fixed-parameter tractable algorithm based on integer linear programming for the unweighted $\text{max-min BCP}_k$ parameterized by the size of a vertex cover.


title: Learning Motion Priors for 4D Human Body Capture in 3D Scenes

abstract: Recovering high-quality 3D human motion in complex scenes from monocular videos is important for many applications, ranging from AR/VR to robotics. However, capturing realistic human-scene interactions, while dealing with occlusions and partial views, is challenging; current approaches are still far from achieving compelling results. We address this problem by proposing LEMO: LEarning human MOtion priors for 4D human body capture. By leveraging the large-scale motion capture dataset AMASS, we introduce a novel motion smoothness prior, which strongly reduces the jitters exhibited by poses recovered over a sequence. Furthermore, to handle contacts and occlusions occurring frequently in body-scene interactions, we design a contact friction term and a contact-aware motion infiller obtained via per-instance self-supervised training. To prove the effectiveness of the proposed motion priors, we combine them into a novel pipeline for 4D human body capture in 3D scenes. With our pipeline, we demonstrate high-quality 4D human body capture, reconstructing smooth motions and physically plausible body-scene interactions. The code and data are available at https://sanweiliti.github.io/LEMO/LEMO.html.


title: Homology cobordism, knot concordance, and Heegaard Floer homology

abstract: We review some recent results in knot concordance and homology cobordism. The proofs rely on various forms of Heegaard Floer homology. We also discuss related open problems.


title: Quadratic forms in 8 prime variables

abstract: We give an asymptotic for the number of prime solutions to $Q(x_1,\dots, x_8) = N$, subject to a mild non-degeneracy condition on the homogeneous quadratic form $Q$. The argument initially proceeds via the circle method, but this does not suffice by itself. To obtain a nontrivial bound on certain averages of exponential sums, we interpret these sums as matrix coefficients for the Weil representation of the symplectic group $\operatorname{Sp}_8(\mathbf{Z}/q\mathbf{Z})$. Averages of such matrix coefficients are then bounded using an amplification argument and a convergence result for convolutions of measures, which reduces matters to understanding the action of certain 12-dimensional subgroups in the Weil representation. Sufficient understanding can be gained by using the basic represention theory of $\operatorname{SL}_2(k)$, $k$ a finite field.


title: Matching Bayesian and frequentist coverage probabilities when using an approximate data covariance matrix

abstract: Observational astrophysics consists of making inferences about the Universe by comparing data and models. The credible intervals placed on model parameters are often as important as the maximum a posteriori probability values, as the intervals indicate concordance or discordance between models and with measurements from other data. Intermediate statistics (e.g. the power spectrum) are usually measured and inferences made by fitting models to these rather than the raw data, assuming that the likelihood for these statistics has multivariate Gaussian form. The covariance matrix used to calculate the likelihood is often estimated from simulations, such that it is itself a random variable. This is a standard problem in Bayesian statistics, which requires a prior to be placed on the true model parameters and covariance matrix, influencing the joint posterior distribution. As an alternative to the commonly-used Independence-Jeffreys prior, we introduce a prior that leads to a posterior that has approximately frequentist matching coverage. This is achieved by matching the covariance of the posterior to that of the distribution of true values of the parameters around the maximum likelihood values in repeated trials, under certain assumptions. Using this prior, credible intervals derived from a Bayesian analysis can be interpreted approximately as confidence intervals, containing the truth a certain proportion of the time for repeated trials. Linking frequentist and Bayesian approaches that have previously appeared in the astronomical literature, this offers a consistent and conservative approach for credible intervals quoted on model parameters for problems where the covariance matrix is itself an estimate.


title: Robust Risk-Aware Reinforcement Learning

abstract: We present a reinforcement learning (RL) approach for robust optimisation of risk-aware performance criteria. To allow agents to express a wide variety of risk-reward profiles, we assess the value of a policy using rank dependent expected utility (RDEU). RDEU allows the agent to seek gains, while simultaneously protecting themselves against downside risk. To robustify optimal policies against model uncertainty, we assess a policy not by its distribution, but rather, by the worst possible distribution that lies within a Wasserstein ball around it. Thus, our problem formulation may be viewed as an actor/agent choosing a policy (the outer problem), and the adversary then acting to worsen the performance of that strategy (the inner problem). We develop explicit policy gradient formulae for the inner and outer problems, and show its efficacy on three prototypical financial problems: robust portfolio allocation, optimising a benchmark, and statistical arbitrage.

categories: cs.LG q-fin.CP q-fin.PM q-fin.RM q-fin.ST


title: Methods to stabilize aqueous supercooling identified by use of an isochoric nucleation detection (INDe) device

abstract: Stable aqueous supercooling has shown significant potential as a technique for human tissue preservation, food cold storage, conservation biology, and beyond, but its stochastic nature has made its translation outside the laboratory difficult. In this work, we present an isochoric nucleation detection (INDe) platform for automated, high-throughput characterization of aqueous supercooling at >1 mL volumes, which enables statistically-powerful determination of the temperatures and time periods for which supercooling in a given aqueous system will remain stable. We employ the INDe to investigate the effects of thermodynamic, surface, and chemical parameters on aqueous supercooling, and demonstrate that various simple system modifications can significantly enhance supercooling stability, including isochoric (constant-volume) confinement, hydrophobic container walls, and the addition of even mild concentrations of solute. Finally, in order to enable informed design of stable supercooled biopreservation protocols, we apply a statistical model to estimate stable supercooling durations as a function of temperature and solution chemistry, producing first-in-field supercooling stability maps for four common cryoprotective solutes.


title: Spectral Properties of Symmetric Quantum States and Symmetric Entanglement Witnesses

abstract: We introduce and explore two questions concerning spectra of operators that are of interest in the theory of entanglement in symmetric (i.e., bosonic) quantum systems. First, we investigate the inverse eigenvalue problem for symmetric entanglement witnesses -- that is, we investigate what their possible spectra are. Second, we investigate the problem of characterizing which separable symmetric quantum states remain separable after conjugation by an arbitrary unitary acting on symmetric space -- that is, which states are separable in every orthonormal symmetric basis. Both of these questions have been investigated thoroughly in the non-symmetric setting, and we contrast the answers that we find with their non-symmetric counterparts.


title: Hypergraph Tur\'an Problems in $\ell_2$-Norm

abstract: There are various different notions measuring extremality of hypergraphs. In this survey we compare the recently introduced notion of the codegree squared extremal function with the Tur\'an function, the minimum codegree threshold and the uniform Tur\'an density. The codegree squared sum $\textrm{co}_2(G)$ of a $3$-uniform hypergraph $G$ is defined to be the sum of codegrees squared $d(x,y)^2$ over all pairs of vertices $x,y$. In other words, this is the square of the $\ell_2$-norm of the codegree vector. We are interested in how large $\textrm{co}_2(G)$ can be if we require $G$ to be $H$-free for some $3$-uniform hypergraph $H$. This maximum value of $\textrm{co}_2(G)$ over all $H$-free $n$-vertex $3$-uniform hypergraphs $G$ is called the codegree squared extremal function, which we denote by $\textrm{exco}_2(n,H)$. We systemically study the extremal codegree squared sum of various $3$-uniform hypergraphs using various proof techniques. Some of our proofs rely on the flag algebra method while others use more classical tools such as the stability method. In particular, we (asymptotically) determine the codegree squared extremal numbers of matchings, stars, paths, cycles, and $F_5$, the $5$-vertex hypergraph with edge set $\{123,124,345\}$. Additionally, our paper has a survey format, as we state several conjectures and give an overview of Tur\'an densities, minimum codegree thresholds and codegree squared extremal numbers of popular hypergraphs. We intend to update the arXiv version of this paper regularly.


title: Measurement of the top quark mass using events with a single reconstructed top quark in pp collisions at $\sqrt{s}$ = 13 TeV

abstract: A measurement of the top quark mass is performed using a data sample enriched with single top quark events produced in the $t$ channel. The study is based on proton-proton collision data, corresponding to an integrated luminosity of 35.9 fb$^{-1}$, recorded at $\sqrt{s}$ = 13 TeV by the CMS experiment at the LHC in 2016. Candidate events are selected by requiring an isolated high-momentum lepton (muon or electron) and exactly two jets, of which one is identified as originating from a bottom quark. Multivariate discriminants are designed to separate the signal from the background. Optimized thresholds are placed on the discriminant outputs to obtain an event sample with high signal purity. The top quark mass is found to be 172.13 $^{+0.76}_{-0.77}$ GeV, where the uncertainty includes both the statistical and systematic components, reaching sub-GeV precision for the first time in this event topology. The masses of the top quark and antiquark are also determined separately using the lepton charge in the final state, from which the mass ratio and difference are determined to be 0.9952 $^{+0.0079}_{-0.0104}$ and 0.83 $^{+1.79}_{-1.35}$ GeV, respectively. The results are consistent with $CPT$ invariance.


title: Solving Tur\'an's Tetrahedron Problem for the $\ell_2$-Norm

abstract: Tur\'an's famous tetrahedron problem is to compute the Tur\'an density of the tetrahedron $K_4^3$. This is equivalent to determining the maximum $\ell_1$-norm of the codegree vector of a $K_4^3$-free $n$-vertex $3$-uniform hypergraph. We introduce a new way for measuring extremality of hypergraphs and determine asymptotically the extremal function of the tetrahedron in our notion. The codegree squared sum, $\text{co}_2(G)$, of a $3$-uniform hypergraph $G$ is the sum of codegrees squared $d(x,y)^2$ over all pairs of vertices $xy$, or in other words, the square of the $\ell_2$-norm of the codegree vector of the pairs of vertices. We define $\text{exco}_2(n,H)$ to be the maximum $\text{co}_2(G)$ over all $H$-free $n$-vertex $3$-uniform hypergraphs $G$. We use flag algebra computations to determine asymptotically the codegree squared extremal number for $K_4^3$ and $K_5^3$ and additionally prove stability results. In particular, we prove that the extremal $K_4^3$-free hypergraphs in $\ell_2$-norm have approximately the same structure as one of the conjectured extremal hypergraphs for Tur\'an's conjecture. Further, we prove several general properties about $\text{exco}_2(n,H)$ including the existence of a scaled limit, blow-up invariance and a supersaturation result.


title: Large spin-charge interconversion induced by interfacial spin-orbit coupling at a highly conducting all-metallic system

abstract: Spin-charge interconversion in systems with spin-orbit coupling has provided a new route for the generation and detection of spin currents in functional devices for memory and logic such as spin-orbit torque switching in magnetic memories or magnetic state reading in spin-based logic. Disentangling the bulk (spin Hall effect) from the interfacial (inverse spin galvanic effect) contribution has been a common issue to properly quantify the spin-charge interconversion in these systems, being the case of Au paradigmatic. Here, we obtain a large spin-charge interconversion at a highly conducting Au/Cu interface which is experimentally shown to arise from the inverse spin galvanic effect of the interface and not from the spin Hall effect of bulk Au. We use two parameters independent of the microscopic details to properly quantify the spin-charge interconversion and the spin losses due to the interfacial spin-orbit coupling, providing an adequate benchmarking to compare with any spin-charge interconversion system. The good performance of this metallic interface, not based in Bi, opens the path to the use of much simpler light/heavy metal systems.


title: Deformable hard particles particles confined in a disordered porous matrix

abstract: With suitably designed Monte Carlo simulations we have investigated the properties of mobile, impenetrable, yet deformable particles that are immersed into a porous matrix, the latter one realized via a frozen configuration of spherical particles. By virtue of a model put forward by Batista and Miller [Phys. Rev. Lett. {\bf 105}, 088305 (2010)] the fluid particles can change under the impact of their surrounding (i.e., either other fluid particles or the matrix) their shape within the class of ellipsoids of revolution; such a change in shape is related to an energy change which is fed into suitably defined selection rules in the deformation "moves" of the Monte Carlo simulations. This concept represents a simple, yet powerful model of realistic, deformable molecules with complex internal structures (such as dendrimers or polymers). For the evaluation of the properties of the system we have used the well-known quenched-annealed protocol (with its characteristic double average prescription) and have analysed the simulation data in terms of static properties (radial distribution function and aspect ratio distribution of the ellipsoids) and dynamic features (notably the mean squared displacement). Our data provide evidence that the degree of deformability of the fluid particles has a distinct impact on the aforementioned properties of the system.


title: StreaMRAK a Streaming Multi-Resolution Adaptive Kernel Algorithm

abstract: Kernel ridge regression (KRR) is a popular scheme for non-linear non-parametric learning. However, existing implementations of KRR require that all the data is stored in the main memory, which severely limits the use of KRR in contexts where data size far exceeds the memory size. Such applications are increasingly common in data mining, bioinformatics, and control. A powerful paradigm for computing on data sets that are too large for memory is the streaming model of computation, where we process one data sample at a time, discarding each sample before moving on to the next one. In this paper, we propose StreaMRAK - a streaming version of KRR. StreaMRAK improves on existing KRR schemes by dividing the problem into several levels of resolution, which allows continual refinement to the predictions. The algorithm reduces the memory requirement by continuously and efficiently integrating new samples into the training model. With a novel sub-sampling scheme, StreaMRAK reduces memory and computational complexities by creating a sketch of the original data, where the sub-sampling density is adapted to the bandwidth of the kernel and the local dimensionality of the data. We present a showcase study on two synthetic problems and the prediction of the trajectory of a double pendulum. The results show that the proposed algorithm is fast and accurate.


title: The Kato-Ponce Inequality with Polynomial Weights

abstract: We consider various versions of fractional Leibniz rules (also known as Kato-Ponce inequalities) with polynomial weights $\langle x\rangle^a = (1+|x|^2)^{a/2}$ for $a\ge 0$. We show that the weighted Kato-Ponce estimate with the inhomogeneous Bessel potential $J^s = (1- \De)^{{s}/{2}}$ holds for the full range of bilinear Lebesgue exponents, for all polynomial weights, and for the sharp range of the degree $s$. This result, in particular, demonstrates that neither the classical Muckenhoupt weight condition nor the more general multilinear weight condition is required for the weighted Kato-Ponce inequality. We also consider a few other variants such as commutator and mixed norm estimates, and analogous conclusions are derived. Our results contain strong-type inequalities for both $L^1$ and $L^\infty$ endpoints, which extend several existing results.


title: The structure of the Milky Way based on unWISE 3.4$\mu$m integrated photometry

abstract: We present a detailed analysis of the Galaxy structure using an unWISE wide-field image at $3.4\mu$m. We perform a 3D photometric decomposition of the Milky Way taking into account i) the projection of the Galaxy on the celestial sphere and ii) that the observer is located within the Galaxy at the solar radius. We consider a large set of photometric models starting with a pure disc model and ending with a complex model which consists of thin and thick discs plus a boxy-peanut-shaped bulge. In our final model, we incorporate many observed features of the Milky Way, such as the disc flaring and warping, several over-densities in the plane, and the dust extinction. The model of the bulge with the corresponding X-shape structure is obtained from N-body simulations of a Milky Way-like galaxy. This allows us to retrieve the parameters of the aforementioned stellar components, estimate their contribution to the total Galaxy luminosity, and constrain the position angle of the bar. The mass of the thick disc in our models is estimated to be 0.4-1.3 of that for the thin disc. The results of our decomposition can be directly compared to those obtained for external galaxies via multicomponent photometric decomposition.


title: Multi-Criteria Radio Spectrum Sharing With Subspace-Based Pareto Tracing

abstract: Radio spectrum is a high-demand finite resource. To meet growing demands of data throughput for forthcoming and deployed wireless networks, new wireless technologies must operate in shared spectrum over unlicensed bands (coexist). As an example application, we consider a model of Long-Term Evolution (LTE) License-Assisted Access (LAA) with incumbent IEEE 802.11 wireless-fidelity (Wi-Fi) systems in a coexistence scenario. This scenario considers multiple LAA and Wi-Fi links sharing an unlicensed band; however, a multitude of other scenarios could be applicable to our general approach. We aim to improve coexistence by maximizing the key performance indicators (KPIs) of two networks simultaneously via dimension reduction and multi-criteria optimization. These KPIs are network throughputs as a function of medium access control (MAC) protocols and physical layer parameters. We perform an exploratory analysis of coexistence behavior by approximating active subspaces to identify low-dimensional structure in the optimization criteria, i.e., few linear combinations of parameters for simultaneously maximizing LAA and Wi-Fi throughputs. We take advantage of an aggregate low-dimensional subspace parametrized by approximate active subspaces of both throughputs to regularize a multi-criteria optimization. Additionally, a choice of two-dimensional subspace enables visualizations augmenting interpretability and explainability of the results. The visualizations and approximations suggest a predominantly convex set of KPIs over active coordinates leading to a regularized manifold of approximately Pareto optimal solutions. Subsequent geodesics lead to continuous traces through parameter space constituting non-dominated solutions in contrast to random grid search.


title: Honeycomb Hubbard Model at van Hove Filling II: Lower Bounds of the Self-Energy

abstract: In this paper we complete the proof that the honeycomb Hubbard model with renormalized chemical potential $\mu=1$ is {\it not} a Fermi liquid in the mathematically precise sense of Salmhofer, by establishing the necessary lower bounds for second derivatives of the self-energy w.r.t. the external momentum.


title: Internal rubble properties of asteroid (101955) Bennu

abstract: Exploration of asteroid (101955) Bennu by the OSIRIS-REx mission has provided an in-depth look at this rubble-pile near-Earth asteroid. In particular, the measured gravity field and the detailed shape model of Bennu indicate significant heterogeneities in its interior structure, compatible with a lower density at its center. Here we combine gravity inversion methods with a statistical rubble-pile model to determine the density and size-frequency distribution (SFD) index of the rubble that constitutes Bennu. The best-fitting models indicate that the SFD of the interior is consistent with that observed on the surface, with a cumulative SFD index of approximately $-2.9$. The rubble bulk density is approximately $1.35$ g/cm$^3$, corresponding to a $12$% macro-porosity. We find the largest rubble particle to be approximately $145$ m, whereas the largest void is approximately $10$ m.


title: Recurrent multiple shared layers in Depth for Neural Machine Translation

abstract: Learning deeper models is usually a simple and effective approach to improve model performance, but deeper models have larger model parameters and are more difficult to train. To get a deeper model, simply stacking more layers of the model seems to work well, but previous works have claimed that it cannot benefit the model. We propose to train a deeper model with recurrent mechanism, which loops the encoder and decoder blocks of Transformer in the depth direction. To address the increasing of model parameters, we choose to share parameters in different recursive moments. We conduct our experiments on WMT16 English-to-German and WMT14 English-to-France translation tasks, our model outperforms the shallow Transformer-Base/Big baseline by 0.35, 1.45 BLEU points, which is 27.23% of Transformer-Big model parameters. Compared to the deep Transformer(20-layer encoder, 6-layer decoder), our model has similar model performance and infer speed, but our model parameters are 54.72% of the former.


title: On the Efficiency of 5(4) RK-Embedded Pairs with High Order Compact Scheme and Robin Boundary Condition for Options Valuation

abstract: When solving the American options with or without dividends, numerical methods often obtain lower convergence rates if further treatment is not implemented even using high-order schemes. In this article, we present a fast and explicit fourth-order compact scheme for solving the free boundary options. In particular, the early exercise features with the asset option and option sensitivity are computed based on a coupled of nonlinear PDEs with fixed boundaries for which a high order analytical approximation is obtained. Furthermore, we implement a new treatment at the left boundary by introducing a third-order Robin boundary condition. Rather than computing the optimal exercise boundary from the analytical approximation, we simply obtain it from the asset option based on the linear relationship at the left boundary. As such, a high order convergence rate can be achieved. We validate by examples that the improvement at the left boundary yields a fourth-order convergence rate without further implementation of mesh refinement, Rannacher time-stepping, and/or smoothing of the initial condition. Furthermore, we extensively compare, the performance of our present method with several 5(4) Runge-Kutta pairs and observe that Dormand and Prince and Bogacki and Shampine 5(4) pairs are faster and provide more accurate numerical solutions. Based on numerical results and comparison with other existing methods, we can validate that the present method is very fast and provides more accurate solutions with very coarse grids.

categories: q-fin.CP cs.NA math.NA q-fin.MF q-fin.PR


title: A Seascape Origin of Richards Growth

abstract: First proposed as an empirical rule over half a century ago, the Richards growth equation has been frequently invoked in population modeling and pandemic forecasting. Central to this model is the advent of a fractional exponent $\gamma$, typically fitted to the data. While various motivations for this non-analytical form have been proposed, it is still considered foremost an empirical fitting procedure. Here, we find that Richards-like growth laws emerge naturally from generic analytical growth rules in a distributed population, upon inclusion of {\bf (i)} migration (spatial diffusion) amongst different locales, and {\bf (ii)} stochasticity in the growth rate, also known as "seascape noise." The latter leads to a wide (power-law) distribution in local population number that, while smoothened through the former, can still result in a fractional growth law for the overall population. This justification of the Richards growth law thus provides a testable connection to the distribution of constituents of the population.


title: Jointly Learnable Data Augmentations for Self-Supervised GNNs

abstract: Self-supervised Learning (SSL) aims at learning representations of objects without relying on manual labeling. Recently, a number of SSL methods for graph representation learning have achieved performance comparable to SOTA semi-supervised GNNs. A Siamese network, which relies on data augmentation, is the popular architecture used in these methods. However, these methods rely on heuristically crafted data augmentation techniques. Furthermore, they use either contrastive terms or other tricks (e.g., asymmetry) to avoid trivial solutions that can occur in Siamese networks. In this study, we propose, GraphSurgeon, a novel SSL method for GNNs with the following features. First, instead of heuristics we propose a learnable data augmentation method that is jointly learned with the embeddings by leveraging the inherent signal encoded in the graph. In addition, we take advantage of the flexibility of the learnable data augmentation and introduce a new strategy that augments in the embedding space, called post augmentation. This strategy has a significantly lower memory overhead and run-time cost. Second, as it is difficult to sample truly contrastive terms, we avoid explicit negative sampling. Third, instead of relying on engineering tricks, we use a scalable constrained optimization objective motivated by Laplacian Eigenmaps to avoid trivial solutions. To validate the practical use of GraphSurgeon, we perform empirical evaluation using 14 public datasets across a number of domains and ranging from small to large scale graphs with hundreds of millions of edges. Our finding shows that GraphSurgeon is comparable to six SOTA semi-supervised and on par with five SOTA self-supervised baselines in node classification tasks. The source code is available at https://github.com/zekarias-tilahun/graph-surgeon.


title: Mid-Cycle Observations of CR Boo and Estimation of the System's Parameters

abstract: We present observations (with NAO Rozhen and AS Vidojevica telescopes) of the AM Canum Venaticorum (AM CVn) type binary star CR Bootis (CR Boo) in the UBV bands. The data were obtained in two nights in July 2019, when the V band brightness was in the range of 16.1-17.0 mag. In both nights, a variability for a period of $25 (\pm 1)$ min and amplitude of about 0.2 magnitudes was visible. These brightness variations are most likely indications of "humps". During our observational time, they appear for a period similar to the CR Boo orbital period. A possible reason of their origin is the phase rotation of the bright spot, placed in the contact point of the infalling matter and the outer disc edge. We estimated some of the parameters of the binary system, on the base of the observational data.


title: An Empirical Assessment of Endpoint Security Systems Against Advanced Persistent Threats Attack Vectors

abstract: Advanced persistent threats pose a significant challenge for blue teams as they apply various attacks over prolonged periods, impeding event correlation and their detection. In this work, we leverage various diverse attack scenarios to assess the efficacy of EDRs and other endpoint security solutions against detecting and preventing APTs. Our results indicate that there is still a lot of room for improvement as state of the art endpoint security systems fail to prevent and log the bulk of the attacks that are reported in this work. Additionally, we discuss methods to tamper with the telemetry providers of EDRs, allowing an adversary to perform a more stealth attack.


title: On the Foundation of Sparse Sensing (Part I): Necessary and Sufficient Sampling Theory and Robust Remaindering Problem

abstract: In the first part of the series papers, we set out to answer the following question: given specific restrictions on a set of samplers, what kind of signal can be uniquely represented by the corresponding samples attained, as the foundation of sparse sensing. It is different from compressed sensing, which exploits the sparse representation of a signal to reduce sample complexity (compressed sampling or acquisition). We use sparse sensing to denote a board concept of methods whose main focus is to improve the efficiency and cost of sampling implementation itself. The "sparse" here is referred to sampling at a low temporal or spatial rate (sparsity constrained sampling or acquisition), which in practice models cheaper hardware such as lower power, less memory and throughput. We take frequency and direction of arrival (DoA) estimation as concrete examples and give the necessary and sufficient requirements of the sampling strategy. Interestingly, we prove that these problems can be reduced to some (multiple) remainder model. As a straightforward corollary, we supplement and complete the theory of co-prime sampling, which receives considerable attention over last decade. On the other hand, we advance the understanding of the robust multiple remainder problem, which models the case when sampling with noise. A sharpened tradeoff between the parameter dynamic range and the error bound is derived. We prove that, for N-frequency estimation in either complex or real waveforms, once the least common multiple (lcm) of the sampling rates selected is sufficiently large, one may approach an error tolerance bound independent of N.


title: Power Grid Cascading Failure Mitigation by Reinforcement Learning

abstract: This paper proposes a cascading failure mitigation strategy based on Reinforcement Learning (RL). The motivation of the Multi-Stage Cascading Failure (MSCF) problem and its connection with the challenge of climate change are introduced. The bottom-level corrective control of the MCSF problem is formulated based on DCOPF (Direct Current Optimal Power Flow). Then, to mitigate the MSCF issue by a high-level RL-based strategy, physics-informed reward, action, and state are devised. Besides, both shallow and deep neural network architectures are tested. Experiments on the IEEE 118-bus system by the proposed mitigation strategy demonstrate a promising performance in reducing system collapses.


title: On the Foundation of Sparse Sensing (Part II): Diophantine Sampling and Array Configuration

abstract: In the second part of the series papers, we set out to study the algorithmic efficiency of sparse sensing. Stemmed from co-prime sensing, we propose a generalized framework, termed Diophantine sensing, which utilizes generic Diophantine equation theory and higher-order sparse ruler to strengthen the sampling time, the degree of freedom (DoF), and the sampling sparsity, simultaneously. Resorting to higher-moment statistics, the proposed Diophantine framework presents two fundamental improvements. First, on frequency estimation, we prove that given arbitrarily large down-sampling rates, there exist sampling schemes where the number of samples needed is only proportional to the sum of DoF and the number of snapshots required, which implies a linear sampling time. Second, on Direction-of-arrival (DoA) estimation, we propose two generic array constructions such that given N sensors, the minimal distance between sensors can be as large as a polynomial of N, O(N^q), which indicates that an arbitrarily sparse array (with arbitrarily small mutual coupling) exists given sufficiently many sensors. In addition, asymptotically, the proposed array configurations produce the best known DoF bound compared to existing sparse array designs.


title: Ultrafast preparation and strong-field ionization of an atomic Bell-like state

abstract: Molecules are many body systems with a substantial amount of entanglement between their electrons. Is there a way to break the molecular bond of a diatomic molecule and obtain two atoms in their ground state which are still entangled and form a Bell-like state? We present a scheme that allows for the preparation of such entangled atomic states from single oxygen molecules on femtosecond time scales. The two neutral oxygen atoms are entangled in the magnetic quantum number of their valence electrons. In a time-delayed probe step, we employ non-adiabatic tunnel ionization, which is a magnetic quantum number-sensitive mechanism. We then investigate correlations by comparing single and double ionization probabilities of the Bell-like state. The experimental results agree with the predictions for an entangled state.


title: Graph-LDA: Graph Structure Priors to Improve the Accuracy in Few-Shot Classification

abstract: It is very common to face classification problems where the number of available labeled samples is small compared to their dimension. These conditions are likely to cause underdetermined settings, with high risk of overfitting. To improve the generalization ability of trained classifiers, common solutions include using priors about the data distribution. Among many options, data structure priors, often represented through graphs, are increasingly popular in the field. In this paper, we introduce a generic model where observed class signals are supposed to be deteriorated with two sources of noise, one independent of the underlying graph structure and isotropic, and the other colored by a known graph operator. Under this model, we derive an optimal methodology to classify such signals. Interestingly, this methodology includes a single parameter, making it particularly suitable for cases where available data is scarce. Using various real datasets, we showcase the ability of the proposed model to be implemented in real world scenarios, resulting in increased generalization accuracy compared to popular alternatives.


title: Coherent energy loss effects in dihadron azimuthal angular correlations in Deep Inelastic Scattering at small $x$

abstract: We perform an exploratory study of the role of coherent, medium-induced energy loss in azimuthal angular correlations in dihadron production in Deep Inelastic Scattering (DIS) at small $x$ where the target proton/nucleus is modeled as a Color Glass Condensate. In this approach coherent radiative energy loss is part of the higher order corrections to the leading order dihadron production cross section. We include the effects of both gluon saturation and coherent radiative energy loss and show that radiative cold-matter energy loss has a significant effect on the so-called coincidence probability for the back to back production of dihadrons in DIS. We also define a double ratio of coincidence probabilities for a nucleus and proton targets and show that it is very robust against higher order radiative corrections.


title: High-accuracy simulations of highly spinning binary neutron star systems

abstract: With an increasing number of expected gravitational-wave detections of binary neutron star mergers, it is essential that gravitational-wave models employed for the analysis of observational data are able to describe generic compact binary systems. This includes systems in which the individual neutron stars are millisecond pulsars for which spin effects become essential. In this work, we perform numerical-relativity simulations of binary neutron stars with aligned and anti-aligned spins within a range of dimensionless spins of $\chi \sim [-0.28,0.58]$. The simulations are performed with multiple resolutions, show a clear convergence order and, consequently, can be used to test existing waveform approximants. We find that for very high spins gravitational-wave models that have been employed for the interpretation of GW170817 and GW190425 are not capable of describing our numerical-relativity dataset. We verify through a full parameter estimation study in which clear biases in the estimate of the tidal deformability and effective spin are present. We hope that in preparation of the next gravitational-wave observing run of the Advanced LIGO and Advanced Virgo detectors our new set of numerical-relativity data can be used to support future developments of new gravitational-wave models.


title: CoverTheFace: face covering monitoring and demonstrating using deep learning and statistical shape analysis

abstract: Wearing a mask is a strong protection against the COVID-19 pandemic, even though the vaccine has been successfully developed and is widely available. However, many people wear them incorrectly. This observation prompts us to devise an automated approach to monitor the condition of people wearing masks. Unlike previous studies, our work goes beyond mask detection; it focuses on generating a personalized demonstration on proper mask-wearing, which helps people use masks better through visual demonstration rather than text explanation. The pipeline starts from the detection of face covering. For images where faces are improperly covered, our mask overlay module incorporates statistical shape analysis (SSA) and dense landmark alignment to approximate the geometry of a face and generates corresponding face-covering examples. Our results show that the proposed system successfully identifies images with faces covered properly. Our ablation study on mask overlay suggests that the SSA model helps to address variations in face shapes, orientations, and scales. The final face-covering examples, especially half profile face images, surpass previous arts by a noticeable margin.


title: Theory of mirror benchmarking and demonstration on a quantum computer

abstract: A new class of protocols called mirror benchmarking was recently proposed to measure the system-level performance of quantum computers. These protocols involve circuits with random sequences of gates followed by mirroring, that is, inverting each gate in the sequence. We give a simple proof that mirror benchmarking leads to an exponential decay of the survival probability with sequence length, under the uniform noise assumption, provided the twirling group forms a 2-design. The decay rate is determined by a quantity that is a quadratic function of the error channel, and for certain types of errors is equal to the unitarity. This result yields a new method for estimating the coherence of noise. We present data from mirror benchmarking experiments run on the Honeywell System Model H1. This data constitutes a set of performance curves, indicating the success probability for random circuits as a function of qubit number and circuit depth.


title: Resource Allocation in Heterogeneously-Distributed Joint Radar-Communications under Asynchronous Bayesian Tracking Framework

abstract: Optimal allocation of shared resources is key to deliver the promise of jointly operating radar and communications systems. In this paper, unlike prior works which examine synergistic access to resources in colocated joint radar-communications or among identical systems, we investigate this problem for a distributed system comprising heterogeneous radars and multi-tier communications. In particular, we focus on resource allocation in the context of multi-target tracking (MTT) while maintaining stable communication connections. By simultaneously allocating the available power, dwell time and shared bandwidth, we improve the MTT performance under a Bayesian tracking framework and guarantee the communications throughput. Our alternating allocation of heterogeneous resources (ANCHOR) approach solves the resulting nonconvex problem based on the alternating optimization method that monotonically improves the Bayesian Cram\'er-Rao bound. Numerical experiments demonstrate that ANCHOR significant improves the tracking error over two baseline allocations and stability under different target scenarios and radar-communications network distributions.


title: Computational multiphase periporomechanics for unguided cracking in unsaturated porous media

abstract: In this article we formulate and implement a computational multiphase periporomechanics model for unguided fracturing in unsaturated porous media. The same governing equation for the solid phase applies on and off cracks. Crack formation in this framework is autonomous, requiring no prior estimates of crack topology. As a new contribution, an energy-based criterion for arbitrary crack formation is formulated using the peridynamic effective force state for unsaturated porous media. Unsaturated fluid flow in the fracture space is modeled in a simplified way in line with the nonlocal formulation of unsaturated fluid flow in bulk. The formulated unsaturated fracturing periporomechanics is numerically implemented through a fractional step algorithm in time and a two-phase mixed meshless method in space. The two-stage operator split converts the coupled periporomechanics problem into an undrained deformation and fracture problem and an unsaturated fluid flow in the deformed skeleton configuration. Numerical simulations of in-plane open and shear cracking are conducted to validate the accuracy and robustness of the fracturing unsaturated periporomechanics model. Then numerical examples of wing cracking and non-planar cracking in unsaturated soil specimens are presented to demonstrate the efficacy of the proposed multiphase periporomechanics model for unguided cracking in unsaturated porous media.


title: Adaptive shot allocation for fast convergence in variational quantum algorithms

abstract: Variational Quantum Algorithms (VQAs) are a promising approach for practical applications like chemistry and materials science on near-term quantum computers as they typically reduce quantum resource requirements. However, in order to implement VQAs, an efficient classical optimization strategy is required. Here we present a new stochastic gradient descent method using an adaptive number of shots at each step, called the global Coupled Adaptive Number of Shots (gCANS) method, which improves on prior art in both the number of iterations as well as the number of shots required. These improvements reduce both the time and money required to run VQAs on current cloud platforms. We analytically prove that in a convex setting gCANS achieves geometric convergence to the optimum. Further, we numerically investigate the performance of gCANS on some chemical configuration problems. We also consider finding the ground state for an Ising model with different numbers of spins to examine the scaling of the method. We find that for these problems, gCANS compares favorably to all of the other optimizers we consider.


title: Electric circuit emulation of topological transitions driven by quantum statistics

abstract: Topological phases exhibit a plethora of striking phenomena including disorder-robust localization and propagation of waves of various nature. Of special interest are the transitions between the different topological phases which are typically controlled by the external parameters. In contrast, in this Letter, we predict the topological transition in the two-particle interacting system driven by the particles' quantum statistics. As a toy model, we investigate an extended one-dimensional Hubbard model with two anyonic excitations obeying fractional quantum statistics in-between bosons and fermions. As we demonstrate, the interplay of two-particle interactions and tunneling processes enables topological edge states of anyon pairs whose existence and localization at one or another edge of the one-dimensional system is governed by the quantum statistics of particles. Since a direct realization of the proposed system is challenging, we develop a rigorous method to emulate the eigenmodes and eigenenergies of anyon pairs with resonant electric circuits.


title: Rewrite Rule Inference Using Equality Saturation

abstract: Many compilers, synthesizers, and theorem provers rely on rewrite rules to simplify expressions or prove equivalences. Developing rewrite rules can be difficult: rules may be subtly incorrect, profitable rules are easy to miss, and rulesets must be rechecked or extended whenever semantics are tweaked. Large rulesets can also be challenging to apply: redundant rules slow down rule-based search and frustrate debugging. This paper explores how equality saturation, a promising technique that uses e-graphs to apply rewrite rules, can also be used to infer rewrite rules. E-graphs can compactly represent the exponentially large sets of enumerated terms and potential rewrite rules. We show that equality saturation efficiently shrinks both sets, leading to faster synthesis of smaller, more general rulesets. We prototyped these strategies in a tool dubbed ruler. Compared to a similar tool built on CVC4, ruler synthesizes 5.8X smaller rulesets 25X faster without compromising on proving power. In an end-to-end case study, we show ruler-synthesized rules which perform as well as those crafted by domain experts, and addressed a longstanding issue in a popular open source tool.


title: Longitudinal Distance: Towards Accountable Instance Attribution

abstract: Previous research in interpretable machine learning (IML) and explainable artificial intelligence (XAI) can be broadly categorized as either focusing on seeking interpretability in the agent's model (i.e., IML) or focusing on the context of the user in addition to the model (i.e., XAI). The former can be categorized as feature or instance attribution. Example- or sample-based methods such as those using or inspired by case-based reasoning (CBR) rely on various approaches to select instances that are not necessarily attributing instances responsible for an agent's decision. Furthermore, existing approaches have focused on interpretability and explainability but fall short when it comes to accountability. Inspired in case-based reasoning principles, this paper introduces a pseudo-metric we call Longitudinal distance and its use to attribute instances to a neural network agent's decision that can be potentially used to build accountable CBR agents.


title: Unidirectional valley-contrasting photo-current in strained transition metal dichalcogenide monolayers

abstract: We examine the full static non-linear optical response of uniaxially strained transition metal dichalcogenide monolayers doped with a finite carrier density in the conduction band, in the presence of disorder. We find that the customary shift current is suppressed, yet we identify a strong, valley-dependent non-reciprocal response, which we term a \textit{unidirectional valley-contrasting photo-current} (UVCP). This DC current originates from the combined effect of strain and Kramers symmetry breaking by trigonal warping, while the contributions due to individual valleys can be separated by introducing an energy offset between them by means of a magnetization. This latter fact enables one to monitor inter-valley transitions. The UVCP is proportional to the mobility and is enhanced by the excitonic Coulomb interaction and inter-valley scattering, as well as by a top gate bias. We discuss detection strategies in state-of-the-art experiments.


title: Large Weyl sums and Hausdorff dimension

abstract: We obtain the exact value of the Hausdorff dimension of the set of coefficients of Gauss sums which for a given $\alpha \in (1/2,1)$ achieve the order at least $N^{\alpha}$ for infinitely many sum lengths $N$. For Weyl sums with polynomials of degree $d\ge 3$ we obtain a new upper bound on the Hausdorff dimension of the set of polynomial coefficients corresponding to large values of Weyl sums. Our methods also work for monomial sums, match the previously known lower bounds, just giving exact value for the corresponding Hausdorff dimension when $\alpha$ is close to $1$. We also obtain a nearly tight bound in a similar question with arbitrary integer sequences of polynomial growth.


title: Foaming and antifoaming in non-aqueous liquids

abstract: Investigation into the physics of foaming has traditionally been focused on aqueous systems. Non-aqueous foams, by contrast, are not well understood, but have been the subject of a recent surge in interest motivated by the need to manage foaming across industrial applications. In this review, we provide a comprehensive discussion of the current state-of-the-art methods for characterizing non-aqueous foams, with a critical evaluation of the advantages and limitations of each. Subsequently we present a concise overview of the current understanding of the mechanisms and methods used for stabilizing and destabilizing non-aqueous foams. We conclude the review by discussing open questions to guide future investigations.


title: Interseries dipole transitions from yellow to green excitons in cuprous oxide

abstract: We study dipole interseries transitions between the yellow and green exciton series in cuprous oxide including the complex valence band structure. To this end, we extend previous studies of the spectrum of complex green exciton resonances [Phys. Rev. B 101, 075208 (2020)] to optical transitions between different exciton states in addition to transitions from the crystal ground state. This allows us to augment the calculations on interseries transitions using a hydrogen-like model [Phys. Rev. B 100, 085201 (2019)] by a more comprehensive treatment of the valence band structure.


title: Evaluating Fairness in Argument Retrieval

abstract: Existing commercial search engines often struggle to represent different perspectives of a search query. Argument retrieval systems address this limitation of search engines and provide both positive (PRO) and negative (CON) perspectives about a user's information need on a controversial topic (e.g., climate change). The effectiveness of such argument retrieval systems is typically evaluated based on topical relevance and argument quality, without taking into account the often differing number of documents shown for the argument stances (PRO or CON). Therefore, systems may retrieve relevant passages, but with a biased exposure of arguments. In this work, we analyze a range of non-stochastic fairness-aware ranking and diversity metrics to evaluate the extent to which argument stances are fairly exposed in argument retrieval systems. Using the official runs of the argument retrieval task Touch\'e at CLEF 2020, as well as synthetic data to control the amount and order of argument stances in the rankings, we show that systems with the best effectiveness in terms of topical relevance are not necessarily the most fair or the most diverse in terms of argument stance. The relationships we found between (un)fairness and diversity metrics shed light on how to evaluate group fairness -- in addition to topical relevance -- in argument retrieval settings.


title: Gradient estimates for parabolic problems with Orlicz growth and discontinuous coefficients

abstract: We obtain Calder\'on-Zygmund type estimates for parabolic equations with Orlicz growth, where nonlinearities involved in the equations may be discontinuous for the space and time variables. In addition, we consider parabolic systems with the Uhlenbeck structure.


title: Shape multistability in flexible tubular crystals through interactions of mobile dislocations

abstract: We study avenues to shape multistability and shape-morphing in flexible crystalline membranes of cylindrical topology, enabled by glide mobility of dislocations. Using computational modeling, we obtain states of mechanical equilibrium presenting a wide variety of tubular crystal deformation geometries, due to an interplay of effective defect interactions with out-of-tangent-plane deformations that reorient the tube axis. Importantly, this interplay often stabilizes defect configurations quite distinct from those predicted for a two-dimensional crystal confined to the surface of a rigid cylinder. We find that relative and absolute stability of competing states depend strongly on control parameters such as bending rigidity, applied stress, and spontaneous curvature. Using stable dislocation pair arrangements as building blocks, we demonstrate that targeted macroscopic three-dimensional conformations of thin crystalline tubes can be programmed by imposing certain sparse patterns of defects. Our findings reveal a broad design space for controllable and reconfigurable colloidal tube geometries, with potential relevance also to carbon nanotubes and microtubules.


title: A VLA Survey of Late-time Radio Emission from Superluminous Supernovae and the Host Galaxies

abstract: We present the results of 3 GHz radio continuum observations of 23 superluminous supernovae (SLSNe) and their host galaxies by using the Karl G. Jansky Very Large Array conducted 5-21 years after the explosions. The sample consists of 15 Type I and 8 Type II SLSNe at z < 0.3, providing one of the largest sample of SLSNe with late-time radio data. We detected radio emission from one SLSN (PTF10hgi) and 5 hosts with a significance of >5$\sigma$. No time variability is found in late-time radio light curves of the radio-detected sources in a timescale of years except for PTF10hgi, whose variability is reported in a separate study. Comparison of star-formation rates (SFRs) derived from the 3 GHz flux densities with those derived from SED modeling based on UV-NIR data shows that four hosts have an excess of radio SFRs, suggesting obscured star formation. Upper limits for undetected hosts and stacked results show that the majority of the SLSN hosts do not have a significant obscured star formation. By using the 3 GHz upper limits, we constrain the parameters for afterglows arising from interaction between initially off-axis jets and circumstellar medium (CSM). We found that the models with higher energies ($E_{\rm iso} \gtrsim$ several $\times 10^{53}$ erg) and CSM densities ($n \gtrsim 0.01$ cm$^{-3}$) are excluded, but lower energies or CSM densities are not excluded with the current data. We also constrained the models of pulsar wind nebulae powered by a newly born magnetar for a subsample of SLSNe with model predictions in the literature.


title: All You Need is Color: Image based Spatial Gene Expression Prediction using Neural Stain Learning

abstract: "Is it possible to predict expression levels of different genes at a given spatial location in the routine histology image of a tumor section by modeling its stain absorption characteristics?" In this work, we propose a "stain-aware" machine learning approach for prediction of spatial transcriptomic gene expression profiles using digital pathology image of a routine Hematoxylin & Eosin (H&E) histology section. Unlike recent deep learning methods which are used for gene expression prediction, our proposed approach termed Neural Stain Learning (NSL) explicitly models the association of stain absorption characteristics of the tissue with gene expression patterns in spatial transcriptomics by learning a problem-specific stain deconvolution matrix in an end-to-end manner. The proposed method with only 11 trainable weight parameters outperforms both classical regression models with cellular composition and morphological features as well as deep learning methods. We have found that the gene expression predictions from the proposed approach show higher correlations with true expression values obtained through sequencing for a larger set of genes in comparison to other approaches.


title: One TTS Alignment To Rule Them All

abstract: Speech-to-text alignment is a critical component of neural textto-speech (TTS) models. Autoregressive TTS models typically use an attention mechanism to learn these alignments on-line. However, these alignments tend to be brittle and often fail to generalize to long utterances and out-of-domain text, leading to missing or repeating words. Most non-autoregressive endto-end TTS models rely on durations extracted from external sources. In this paper we leverage the alignment mechanism proposed in RAD-TTS as a generic alignment learning framework, easily applicable to a variety of neural TTS models. The framework combines forward-sum algorithm, the Viterbi algorithm, and a simple and efficient static prior. In our experiments, the alignment learning framework improves all tested TTS architectures, both autoregressive (Flowtron, Tacotron 2) and non-autoregressive (FastPitch, FastSpeech 2, RAD-TTS). Specifically, it improves alignment convergence speed of existing attention-based mechanisms, simplifies the training pipeline, and makes the models more robust to errors on long utterances. Most importantly, the framework improves the perceived speech synthesis quality, as judged by human evaluators.


title: Fast Robust Tensor Principal Component Analysis via Fiber CUR Decomposition

abstract: We study the problem of tensor robust principal component analysis (TRPCA), which aims to separate an underlying low-multilinear-rank tensor and a sparse outlier tensor from their sum. In this work, we propose a fast non-convex algorithm, coined Robust Tensor CUR (RTCUR), for large-scale TRPCA problems. RTCUR considers a framework of alternating projections and utilizes the recently developed tensor Fiber CUR decomposition to dramatically lower the computational complexity. The performance advantage of RTCUR is empirically verified against the state-of-the-arts on the synthetic datasets and is further demonstrated on the real-world application such as color video background subtraction.


title: Differential Music: Automated Music Generation Using LSTM Networks with Representation Based on Melodic and Harmonic Intervals

abstract: This paper presents a generative AI model for automated music composition with LSTM networks that takes a novel approach at encoding musical information which is based on movement in music rather than absolute pitch. Melodies are encoded as a series of intervals rather than a series of pitches, and chords are encoded as the set of intervals that each chord note makes with the melody at each timestep. Experimental results show promise as they sound musical and tonal. There are also weaknesses to this method, mainly excessive modulations in the compositions, but that is expected from the nature of the encoding. This issue is discussed later in the paper and is a potential topic for future work.


title: A Novel Solution for the General Diffusion

abstract: The Fisher-KPP equation is a reaction-diffusion equation originally proposed by Fisher to represent allele propagation in genetic hosts or population. It was also proposed by Kolmogorov for more general applications. A novel method for solving nonlinear partial differential equations is applied to produce a unique, approximate solution for the Fisher-KPP equation. Analysis proves the solution is counterintuitive. Although still satisfying the maximum principle, time dependence collapses for all time greater than zero, therefore, the solution is highly irregular and not smooth, invalidating the traveling wave approximation so often employed.


title: Adversarial Robustness of Deep Learning: Theory, Algorithms, and Applications

abstract: This tutorial aims to introduce the fundamentals of adversarial robustness of deep learning, presenting a well-structured review of up-to-date techniques to assess the vulnerability of various types of deep learning models to adversarial examples. This tutorial will particularly highlight state-of-the-art techniques in adversarial attacks and robustness verification of deep neural networks (DNNs). We will also introduce some effective countermeasures to improve the robustness of deep learning models, with a particular focus on adversarial training. We aim to provide a comprehensive overall picture about this emerging direction and enable the community to be aware of the urgency and importance of designing robust deep learning models in safety-critical data analytical applications, ultimately enabling the end-users to trust deep learning classifiers. We will also summarize potential research directions concerning the adversarial robustness of deep learning, and its potential benefits to enable accountable and trustworthy deep learning-based data analytical systems and applications.


title: Chiral Effective Field Theory's Impact on Advancing Quantum Monte Carlo Methods

abstract: Thirty years ago, Steven Weinberg published his seminal paper on "Nuclear Forces from chiral Lagrangians" which has revolutionized the field of theoretical nuclear physics. Nowadays, interactions derived from chiral effective field theory are routinely used to describe nuclear systems ranging from atomic nuclei to the dense matter explored in the core of neutron stars with theoretical uncertainty estimates. In our contribution to the special issue "Celebrating 30 years of Steven Weinberg's paper Nuclear Forces from Chiral Lagrangians", we focus on the impact that chiral effective field theory interactions have played in advancing microscopic studies of atomic nuclei and the nuclear-matter equation of state using quantum Monte Carlo methods.


title: Continuous Treatment Recommendation with Deep Survival Dose Response Function

abstract: We propose a general formulation for continuous treatment recommendation problems in settings with clinical survival data, which we call the Deep Survival Dose Response Function (DeepSDRF). That is, we consider the problem of learning the conditional average dose response (CADR) function solely from historical data in which observed factors (confounders) affect both observed treatment and time-to-event outcomes. The estimated treatment effect from DeepSDRF enables us to develop recommender algorithms with the correction for selection bias. We compared two recommender approaches based on random search and reinforcement learning and found similar performance in terms of patient outcome. We tested the DeepSDRF and the corresponding recommender on extensive simulation studies and the eICU Research Institute (eRI) database. To the best of our knowledge, this is the first time that causal models are used to address the continuous treatment effect with observational data in a medical context.


title: Equilibrium and nonequilibrium quantum correlations between two detectors in curved space time

abstract: We investigate the equilibrium and nonequilibrium quantum information correlations encoded in two-qubit system (near the horizon of a Kerr black hole). We study the impact of mass and the angular momentum, and further the local curvature or accelerations on the behaviors of the quantum correlations between two qubits. We show the quantum information of two qubits is encoded in the space time structure. In nonequilibrium case, the nonequilibrium can also contribute to the correlations.


title: Towards precise galaxy evolution: a comparison between spectral indices of $z\sim1$ galaxies in the IllustrisTNG simulation and the LEGA-C survey

abstract: We present the first comparison of observed stellar continuum spectra of high-redshift galaxies and mock galaxy spectra generated from hydrodynamical simulations. The mock spectra are produced from the IllustrisTNG TNG100 simulation combined with stellar population models and take into account dust attenuation and realistic observational effects (aperture effects and noise). We compare the simulated $D_n4000$ and EW(H$\delta$) of galaxies with $10.5 \leq \log(M_\ast/M_\odot) \leq 11.5$ at $0.6 \leq z \leq 1.0$ to the observed distributions from the LEGA-C survey. TNG100 globally reproduces the observed distributions of spectral indices, implying that the age distribution of galaxies in TNG100 is generally realistic. Yet there are small but significant differences. For old galaxies, TNG100 shows small $D_n4000$ when compared to LEGA-C, while LEGA-C galaxies have larger EW(H$\delta$) at fixed $D_n4000$. There are several possible explanations: 1) LEGA-C galaxies have overall older ages combined with small contributions (a few percent in mass) from younger ($<1$~Gyr) stars, while TNG100 galaxies may not have such young sub-populations; 2) the spectral mismatch could be due to systematic uncertainties in the stellar population models used to convert stellar ages and metallicities to observables. In conclusion, the latest cosmological galaxy formation simulations broadly reproduce the global age distribution of galaxies at $z\sim1$ and, at the same time, the high quality of the latest observed and simulated datasets help constrain stellar population synthesis models as well as the physical models underlying the simulations.


title: Long-Term, in-the-Wild Study of Feedback about Speech Intelligibility for K-12 Students Attending Class via a Telepresence Robot

abstract: Telepresence robots offer presence, embodiment, and mobility to remote users, making them promising options for homebound K-12 students. It is difficult, however, for robot operators to know how well they are being heard in remote and noisy classroom environments. One solution is to estimate the operator's speech intelligibility to their listeners in order to provide feedback about it to the operator. This work contributes the first evaluation of a speech intelligibility feedback system for homebound K-12 students attending class remotely. In our four long-term, in-the-wild deployments we found that students speak at different volumes instead of adjusting the robot's volume, and that detailed audio calibration and network latency feedback are needed. We also contribute the first findings about the types and frequencies of multimodal comprehension cues given to homebound students by listeners in the classroom. By annotating and categorizing over 700 cues, we found that the most common cue modalities were conversation turn timing and verbal content. Conversation turn timing cues occurred more frequently overall, whereas verbal content cues contained more information and might be the most frequent modality for negative cues. Our work provides recommendations for telepresence systems that could intervene to ensure that remote users are being heard.


title: A Fault-Tolerant Honeycomb Memory

abstract: Recently, Hastings & Haah introduced a quantum memory defined on the honeycomb lattice. Remarkably, this honeycomb code assembles weight-six parity checks using only two-local measurements. The sparse connectivity and two-local measurements are desirable features for certain hardware, while the weight-six parity checks enable robust performance in the circuit model. In this work, we quantify the robustness of logical qubits preserved by the honeycomb code using a correlated minimum-weight perfect-matching decoder. Using Monte Carlo sampling, we estimate the honeycomb code's threshold in different error models, and project how efficiently it can reach the "teraquop regime" where trillions of quantum logical operations can be executed reliably. We perform the same estimates for the rotated surface code, and find a threshold of $0.2\%-0.3\%$ for the honeycomb code compared to a threshold of $0.5\%-0.7\%$ for the surface code in a controlled-not circuit model. In a circuit model with native two-body measurements, the honeycomb code achieves a threshold of $1.5\% < p <2.0\%$, where $p$ is the collective error rate of the two-body measurement gate - including both measurement and correlated data depolarization error processes. With such gates at a physical error rate of $10^{-3}$, we project that the honeycomb code can reach the teraquop regime with only $600$ physical qubits.


title: Supernodes

abstract: In this paper, we present two new concepts related to subgraph counting where the focus is not on the number of subgraphs that are isomorphic to some fixed graph $H$, but on the frequency with which a vertex or an edge belongs to such subgraphs. In particular, we are interested in the case where $H$ is a complete graph. These new concepts are termed vertex participation and edge participation respectively. We combine these concepts with that of the rich-club to identify what we call a Super rich-club and rich edge-club. We show that the concept of vertex participation is a generalisation of the rich-club. We present experimental results on randomised Erd\"{o}s R\'{e}nyi and Watts-Strogatz small-world networks. We further demonstrate both concepts on a complex brain network and compare our results to the rich-club of the brain.


title: Towards automated sampling of polymorph nucleation and free energies with SGOOP and metadynamics

abstract: Understanding the driving forces behind the nucleation of different polymorphs is of great importance for material sciences and the pharmaceutical industry. This includes understanding the reaction coordinate that governs the nucleation process as well as correctly calculating the relative free energies of different polymorphs. Here we demonstrate, for the prototypical case of urea nucleation from melt, how one can learn such a 1-dimensional reaction coordinate as a function of pre-specified order parameters, and use it to perform efficient biased all-atom molecular dynamics simulations. The reaction coordinate is learnt as a function of generic thermodynamic and structural order parameters using the "Spectral Gap Optimization of Order Parameters (SGOOP)" approach [P. Tiwary and B. J. Berne, Proc. Natl. Acad. Sci. (2016)], and is biased using well-tempered metadynamics simulations. The reaction coordinate gives insight into the role played by different structural and thermodynamics order parameters, and the biased simulations obtain accurate relative free energies for different polymorphs. This includes accurate prediction of the approximate pressure at which urea undergoes a phase transition and one of the metastable polymorphs becomes the most stable conformation. We believe the ideas demonstrated in thus work will facilitate efficient sampling of nucleation in complex, generic systems.

categories: cond-mat.soft cond-mat.stat-mech physics.chem-ph physics.comp-ph


title: A Refined Einstein-Gauss-Bonnet Inflationary Theoretical Framework

abstract: We provide a refined and much more simplified Einstein-Gauss-Bonnet inflationary theoretical framework, which is compatible with the GW170817 observational constraints on the gravitational wave speed. As in previous works, the constraint that the gravitational wave speed is $c_T^2=1$ in natural units, results to a constraint differential equation that relates the coupling function of the scalar field to the Gauss-Bonnet invariant $\xi(\phi)$ and the scalar potential $V(\phi)$. Adopting the slow-roll conditions for the scalar field and the Hubble rate, and in contrast to previous works, by further assuming that $\kappa \frac{\xi '}{\xi''}\ll 1$, which is motivated by slow-roll arguments, we succeed in providing much more simpler expressions for the slow-roll indices and for the tensor and scalar spectral indices and for the tensor-to-scalar ratio. We exemplify our refined theoretical framework by using an illustrative example with a simple power-law scalar coupling function $\xi(\phi)\sim \phi^{\nu}$ and as we demonstrate the resulting inflationary phenomenology is compatible with the latest Planck data. Moreover, this particular model produces a blue-tilted tensor spectral index, so we discuss in brief the perspective of describing the NANOGrav result with this model as is indicated in the recent literature.


title: Deterministic Dynamic Matching In Worst-Case Update Time

abstract: We present deterministic algorithms for maintaining a $(3/2 + \epsilon)$ and $(2 + \epsilon)$-approximate maximum matching in a fully dynamic graph with worst-case update times $\hat{O}(\sqrt{n})$ and $\tilde{O}(1)$ respectively. The fastest known deterministic worst-case update time algorithms for achieving approximation ratio $(2 - \delta)$ (for any $\delta > 0$) and $(2 + \epsilon)$ were both shown by Roghani et al. [2021] with update times $O(n^{3/4})$ and $O_\epsilon(\sqrt{n})$ respectively. We close the gap between worst-case and amortized algorithms for the two approximation ratios as the best deterministic amortized update times for the problem are $O_\epsilon(\sqrt{n})$ and $\tilde{O}(1)$ which were shown in Bernstein and Stein [SODA'2021] and Bhattacharya and Kiss [ICALP'2021] respectively. In order to achieve both results we explicitly state a method implicitly used in Nanongkai and Saranurak [STOC'2017] and Bernstein et al. [arXiv'2020] which allows to transform dynamic algorithms capable of processing the input in batches to a dynamic algorithms with worst-case update time. \textbf{Independent Work:} Independently and concurrently to our work Grandoni et al. [arXiv'2021] has presented a fully dynamic algorithm for maintaining a $(3/2 + \epsilon)$-approximate maximum matching with deterministic worst-case update time $O_\epsilon(\sqrt{n})$.


title: Mass spectra of doubly heavy tetraquarks in an improved chromomagnetic interaction model

abstract: Doubly heavy tetraquark states are the prime candidates for tightly bound exotic states. We present a systematic study of the mass spectra of the $S$-wave doubly heavy tetraquark states $QQ\bar{q}\bar{q}$ ($q=u, d, s$ and $Q=c, b$) with different quantum numbers $J^P=0^+$, $1^+$, and $2^+$ in the framework of the improved chromomagnetic interaction (ICMI) model. The parameters in the ICMI model are obtained by fitting the conventional hadron spectra and are used directly to predict the masses of the tetraquark states. For heavy quarks, the uncertainties of the parameters are obtained by comparing the masses of doubly and triply heavy baryons with those given by lattice QCD, QCD sum rules, and potential models. Several compact and stable bound states are found in both the doubly charmed and doubly bottomed tetraquark systems. The predicted mass of the $cc\bar u\bar d$ state is consistent with the recent measurement from the LHCb collaboration.


title: Semiclassical measures for higher dimensional quantum cat maps

abstract: Consider a quantum cat map $M$ associated to a matrix $A\in\mathop{\mathrm{Sp}}(2n,\mathbb Z)$, which is a common toy model in quantum chaos. We show that the mass of eigenfunctions of $M$ on any nonempty open set in the position-frequency space satisfies a lower bound which is uniform in the semiclassical limit, under two assumptions: (1) there is a unique simple eigenvalue of $A$ of largest absolute value and (2) the characteristic polynomial of $A$ is irreducible over the rationals. This is similar to previous work [arXiv:1705.05019], [arXiv:1906.08923] on negatively curved surfaces and [arXiv:2103.06633] on quantum cat maps with $n=1$, but this paper gives the first results of this type which apply in any dimension. When condition (2) fails we provide a weaker version of the result and discuss relations to existing counterexamples. We also obtain corresponding statements regarding semiclassical measures and damped quantum cat maps.


title: The Case for Task Sampling based Learning for Cluster Job Scheduling

abstract: The ability to accurately estimate job runtime properties allows a scheduler to effectively schedule jobs. State-of-the-art online cluster job schedulers use history-based learning, which uses past job execution information to estimate the runtime properties of newly arrived jobs. However, with fast-paced development in cluster technology (in both hardware and software) and changing user inputs, job runtime properties can change over time, which lead to inaccurate predictions. In this paper, we explore the potential and limitation of real-time learning of job runtime properties, by proactively sampling and scheduling a small fraction of the tasks of each job. Such a task-sampling-based approach exploits the similarity among runtime properties of the tasks of the same job and is inherently immune to changing job behavior. Our study focuses on two key questions in comparing task-sampling-based learning (learning in space) and history-based learning (learning in time): (1) Can learning in space be more accurate than learning in time? (2) If so, can delaying scheduling the remaining tasks of a job till the completion of sampled tasks be more than compensated by the improved accuracy and result in improved job performance? Our analytical and experimental analysis of 3 production traces with different skew and job distribution shows that learning in space can be substantially more accurate. Our simulation and testbed evaluation on Azure of the two learning approaches anchored in a generic job scheduler using 3 production cluster job traces shows that despite its online overhead, learning in space reduces the average Job Completion Time (JCT) by 1.28x, 1.56x, and 1.32x compared to the prior-art history-based predictor.


abstract: We construct a Weyl transverse diffeomorphism invariant theory of symmetric teleparallel gravity by employing the Weyl compensator formalism. The low-energy dynamics has a single spin two gravition without a scalar degree of freedom. By construction, it is equivalent to the unimodular gravity (as well as the Einstein gravity) at the non-linear level.


title: Asymptotic additivity of the Turaev-Viro invariants for a family of $3$-manifolds

abstract: In this paper, we show that the Turaev-Viro invariant volume conjecture posed by Chen and Yang is preserved under gluings of toroidal boundary components for a family of $3$-manifolds. In particular, we show that the asymptotics of the Turaev-Viro invariants are additive under certain gluings of elementary pieces arising from a construction of hyperbolic cusped $3$-manifolds due to Agol. The gluings of the elementary pieces are known to be additive with respect to the simplicial volume. This allows us to construct families of manifolds with an arbitrary number of hyperbolic pieces such that the resultant manifolds satisfy an extended version of the Turaev-Viro invariant volume conjecture.


title: Simulations of hybrid charge-sensing single-electron-transistors and CMOS circuits

abstract: Single-electron transistors (SETs) have been extensively used as charge sensors in many areas such as quantum computations. In general, the signals of SETs are smaller than those of complementary metal-oxide semiconductor (CMOS) devices, and many amplifying circuits are required to enlarge these signals. Instead of amplifying a single small output, we theoretically consider the amplification of pairs of SETs, such that one of the SETs is used as a reference. We simulate the two-stage amplification process of SETs and CMOS devices using a conventional SPICE (Simulation Program with Integrated Circuit Emphasis) circuit simulator. Implementing the pairs of SETs into CMOS circuits makes the integration of SETs more feasible because of direct signal transfer from the SET to the CMOS circuits.


title: On the narrow spectral feature at $\sim$3 TeV in the MAGIC spectrum of Mrk 501

abstract: Using a time-dependent one-zone leptonic model that incorporates both shock acceleration and stochastic acceleration processes, we investigate the formation of the narrow spectral feature at $\sim3$ TeV of Mrk 501 which was observed during the X-ray and TeV flaring activity in July 2014. It is found that the broadband spectral energy distribution (SED) can be well interpreted as the synchrotron and synchrotron-self-Compton emission from the electron energy distribution (EED) that is composed by a power-law (PL) branch and a pileup branch. The PL branch produces synchrotron photons which are scattered by the electrons of the pileup branch via inverse-Compton scattering and form the narrow spectral feature observed at the TeV energies. The EED is produced by two injection episodes, and the pileup branch in EED is caused by shock acceleration rather than stochastic acceleration.


title: Thermometric machine for ultraprecise thermometry of low temperatures

abstract: Thermal equilibrium states are exponentially hard to distinguish at very low temperatures, making equilibrium quantum thermometry in this regime a formidable task. We present a thermometric scheme that circumvents this limitation, by using a two-level probe that does not thermalize with the sample whose temperature is measured. This is made possible thanks to a suitable interaction that couples the probe to the sample and to an auxiliary thermal bath known to be at a higher temperature. Provided a reasonable upper bound on the temperature of the sample, the resulting 'thermometric machine' drives the probe towards a steady state whose signal-to-noise ratio can achieve values as high as $\mathcal{O}(1/T)$. We also characterize the transient state of the probe and numerically illustrate an extreme reduction in the number of measurements to attain a given precision, as compared to optimal measurements on a thermalized probe.