title: Distribution-free calibration guarantees for histogram binning without sample splitting

abstract: We prove calibration guarantees for the popular histogram binning (also called uniform-mass binning) method of Zadrozny and Elkan [2001]. Histogram binning has displayed strong practical performance, but theoretical guarantees have only been shown for sample split versions that avoid 'double dipping' the data. We demonstrate that the statistical cost of sample splitting is practically significant on a credit default dataset. We then prove calibration guarantees for the original method that double dips the data, using a certain Markov property of order statistics. Based on our results, we make practical recommendations for choosing the number of bins in histogram binning. In our illustrative simulations, we propose a new tool for assessing calibration -- validity plots -- which provide more information than an ECE estimate. Code for this work will be made publicly available at https://github.com/aigen/df-posthoc-calibration.


title: Hydride Conductivity in Nitride Hydrides

abstract: Nitride hydrides are a largely unexplored class of materials with promising applications in solid-state hydrogen fuel cells. Here, we use first-principles calculations to characterize defects and ionic mobility in Sr$_2$LiH$_2$N (SLHN), a nitride hydride with high hydride conductivity. Calculating defect formation energies, we find that SLHN contains high concentrations of hydrogen interstitials (H$_i$). H$_i^-$ migrates with very low energetic barriers, which, together with its low formation energy, implies that SLHN will have excellent hydride kinetics, potentially surpassing those of other known hydride electrolytes. Oxygen contamination is a concern, meaning that encapsulation will be critical. By direct analogy to the La/Sr-based oxyhydrides, which have similar crystal structures, we also investigate the La-based nitride hydride La$_2$LiHN$_2$ but find that it will be significantly less conductive, and thus not as technologically useful. Our findings buttress the exploration of SLHN and similar nitride hydrides for use in solid-state hydrogen fuel cells.


title: The Schwarz-Milnor lemma for braids and area-preserving diffeomorphisms

abstract: We prove a number of new results on the large-scale geometry of the $L^p$-metrics on the group of area-preserving diffeomorphisms of each orientable surface. Our proofs use in a key way the Fulton-MacPherson type compactification of the configuration space of $n$ points on the surface due to Axelrod-Singer and Kontsevich. This allows us to apply the Schwarz-Milnor lemma to configuration spaces, a natural approach which we carry out successfully for the first time. As sample results, we prove that all right-angled Artin groups admit quasi-isometric embeddings into the group of area-preserving diffeomorphisms endowed with the $L^p$-metric, and that all Gambaudo-Ghys quasi-morphisms on this metric group coming from the braid group on $n$ strands are Lipschitz. This was conjectured to hold, yet proven only for low values of $n$ and the genus $g$ of the surface.


title: Different Environment Feedback in Fast-slow Eco-evolutionary Dynamics

abstract: The fast-slow dynamics of an eco-evolutionary system are studied, where we consider the feedback actions of environmental resources that are classified into those that are self-renewing and those externally supplied. We show although these two types of resources are drastically different, the resulting closed-loop systems bear close resemblances, which include the same equilibria and their stability conditions on the boundary of the phase space, and the similar appearances of equilibria in the interior. After closer examination of specific choices of parameter values, we disclose that the global dynamical behaviors of the two types of closed-loop systems can be fundamentally different in terms of limit cycles: the system with self-renewing resources undergoes a generalized Hopf bifurcation such that one stable limit cycle and one unstable limit cycle can coexist; the system with externally supplied resources can only have the stable limit cycle induced by a supercritical Hopf bifurcation. Finally, the explorative analysis is carried out to show the discovered dynamic behaviors are robust in even larger parameter space.


title: Deletion to Scattered Graph Classes I -- case of finite number of graph classes

abstract: Graph-modification problems, where we modify a graph by adding or deleting vertices or edges or contracting edges to obtain a graph in a {\it simpler} class, is a well-studied optimization problem in all algorithmic paradigms including classical, approximation and parameterized complexity. Specifically, graph-deletion problems, where one needs to delete a small number of vertices to make the resulting graph to belong to a given non-trivial hereditary graph class, captures several well-studied problems including {\sc Vertex Cover}, {\sc Feedback Vertex Set}, {\sc Odd Cycle Transveral}, {\sc Cluster Vertex Deletion}, and {\sc Perfect Deletion}. Investigation into these problems in parameterized complexity has given rise to powerful tools and techniques. We initiate a study of a natural variation of the problem of deletion to {\it scattered graph classes}. We want to delete at most $k$ vertices so that in the resulting graph, each connected component belongs to one of a constant number of graph classes. As our main result, we show that this problem is fixed-parameter tractable (FPT) when the deletion problem corresponding to each of the finite number of graph classes is known to be FPT and the properties that a graph belongs to any of the classes is expressible in Counting Monodic Second Order (CMSO) logic. While this is shown using some black box theorems in parameterized complexity, we give a faster FPT algorithm when each of the graph classes has a finite forbidden set.


title: No speedup for geometric theories

abstract: Geometric theories based on classical logic are conservative over their intuitionistic counterparts for geometric implications. The latter result (sometimes referred to as Barr's theorem) is squarely a consequence of Gentzen's Hauptsatz. Prima facie though, cut elimination can result in superexponentially longer proofs. In this paper it is shown that the transformation of a classical proof of a geometric implication in a geometric theory into an intuitionistic proof can be achieved in feasibly many steps.


title: Multi-Objective Controller Synthesis with Uncertain Human Preferences

abstract: Complex real-world applications of cyber-physical systems give rise to the need for multi-objective controller synthesis, which concerns the problem of computing an optimal controller subject to multiple (possibly conflicting) criteria. The relative importance of objectives is often specified by human decision-makers. However, there is inherent uncertainty in human preferences (e.g., due to artifacts resulting from different preference elicitation methods). In this paper, we formalize the notion of uncertain human preferences and present a novel approach that accounts for this uncertainty in the context of multi-objective controller synthesis for Markov decision processes (MDPs). Our approach is based on mixed-integer linear programming and synthesizes an optimally permissive multi-strategy that satisfies uncertain human preferences with respect to a multi-objective property. Experimental results on a range of large case studies show that the proposed approach is feasible and scalable across varying MDP model sizes and uncertainty levels of human preferences. Evaluation via an online user study also demonstrates the quality and benefits of the synthesized controllers.


title: GSPMD: General and Scalable Parallelization for ML Computation Graphs

abstract: We present GSPMD, an automatic, compiler-based parallelization system for common machine learning computations. It allows users to write programs in the same way as for a single device, then give hints through a few annotations on how to distribute tensors, based on which GSPMD will parallelize the computation. Its representation of partitioning is simple yet general, allowing it to express different or mixed paradigms of parallelism on a wide variety of models. GSPMD infers the partitioning for every operator based on limited user annotations, making it convenient to scale existing single-device programs. It solves several technical challenges for production usage, allowing GSPMD to achieve 50% to 62% compute utilization on up to 2048 Cloud TPUv3 cores for models with up to one trillion parameters.


title: On the Theoretical Foundation of Overset Grid Methods for Hyperbolic Problems: Well-Posedness and Conservation

abstract: We use the energy method to study the well-posedness of initial-boundary value problems approximated by overset mesh methods in one and two space dimensions for linear constant-coefficient hyperbolic systems. We show that in one space dimension, for both scalar equations and systems of equations, the problem where one domain partially oversets another is well-posed when characteristic coupling conditions are used. If a system cannot be diagonalized, as is ususally the case in multiple space dimensions, then the energy method does not give proper bounds in terms of initial and boundary data. For those problems, we propose a novel penalty approach. We show, by using a global energy that accounts for the energy in the overlap region of the domains, that under well-defined conditions on the coupling matrices the penalized overset domain problems are energy bounded, conservative, well-posed and have solutions equivalent to the original single domain problem.


title: Correction of the Lusztig-Williamson Billiards Conjecture

abstract: A new algorithm allows us to calculate many new tilting characters for $SL_3$, $SP_4$, $G_2$, $SL_4$ and potentially many other groups. These calculations show that the Lusztig-Williamson Billiards Conjecture needs to be corrected. In this paper we present the new results calculated for $SL_3$ and a correction of the conjecture.


title: The Influence of Memory in Multi-Agent Consensus

abstract: Multi-agent consensus problems can often be seen as a sequence of autonomous and independent local choices between a finite set of decision options, with each local choice undertaken simultaneously, and with a shared goal of achieving a global consensus state. Being able to estimate probabilities for the different outcomes and to predict how long it takes for a consensus to be formed, if ever, are core issues for such protocols. Little attention has been given to protocols in which agents can remember past or outdated states. In this paper, we propose a framework to study what we call \emph{memory consensus protocol}. We show that the employment of memory allows such processes to always converge, as well as, in some scenarios, such as cycles, converge faster. We provide a theoretical analysis of the probability of each option eventually winning such processes based on the initial opinions expressed by agents. Further, we perform experiments to investigate network topologies in which agents benefit from memory on the expected time needed for consensus.


title: Can an Agency Role-Reversal Lead to an Organizational Collapse?; A Study Proposal

abstract: The Principal-Agent Theory model is widely used to explain governance role where there is a separation of ownership and control, as it defines clear boundaries between governance and executives. However, examination of recent corporate failure reveals the concerning contribution of the Board of Directors to such failures and calls into question governance effectiveness in the presence of a powerful and charismatic CEO. This study proposes a framework for analyzing the relationship between the Board of Directors and the CEO, and how certain relationships affect the power structure and behavior of the Board, which leads to a role reversal in the Principal-Agent Theory, as the Board assumes the role of the CEO's agent. This study's results may help create a red flag for a board and leader's behavior that may result in governance failure.


title: HuMoR: 3D Human Motion Model for Robust Pose Estimation

abstract: We introduce HuMoR: a 3D Human Motion Model for Robust Estimation of temporal pose and shape. Though substantial progress has been made in estimating 3D human motion and shape from dynamic observations, recovering plausible pose sequences in the presence of noise and occlusions remains a challenge. For this purpose, we propose an expressive generative model in the form of a conditional variational autoencoder, which learns a distribution of the change in pose at each step of a motion sequence. Furthermore, we introduce a flexible optimization-based approach that leverages HuMoR as a motion prior to robustly estimate plausible pose and shape from ambiguous observations. Through extensive evaluations, we demonstrate that our model generalizes to diverse motions and body shapes after training on a large motion capture dataset, and enables motion reconstruction from multiple input modalities including 3D keypoints and RGB(-D) videos.


title: European Spallation Source: a future for Coherent Neutrino Nucleus Scattering

abstract: The European Spallation Source (ESS), currently finishing its construction, will soon provide the most intense neutron beams for multi-disciplinary science. At the same time, it will also produce a high-intensity neutrino flux with an energy suitable for precision measurements of Coherent Elastic Neutrino-Nucleus Scattering. We describe some physics prospects, within and beyond the Standard Model, of employing innovative detector technologies to take the most out of this large flux. We show that, compared to current measurements, the ESS will provide a much more precise understanding of neutrino and nuclear properties.


title: Holographic correlators with multi-particle states

abstract: We derive the connected tree-level part of 4-point holographic correlators in AdS$_3\times S^3\times \mathcal{M}$ (where ${\cal M}$ is $T^4$ or $K3$) involving two multi-trace and two single-trace operators. These connected correlators are obtained by studying a heavy-heavy-light-light correlation function in the formal limit where the heavy operators become light. These results provide a window into higher-point holographic correlators of single-particle operators. We find that the correlators involving multi-trace operators are compactly written in terms of Bloch-Wigner-Ramakrishnan functions -- particular linear combinations of higher-order polylogarithm functions. Several consistency checks of the derived expressions are performed in various OPE channels. We also extract the anomalous dimensions and 3-point couplings of the non-BPS double-trace operators of lowest twist at order $1/c$ and find some positive anomalous dimensions at spin zero and two in the K3 case.


title: Extending Python for Quantum-Classical Computing via Quantum Just-in-Time Compilation

abstract: Python is a popular programming language known for its flexibility, usability, readability, and focus on developer productivity. The quantum software community has adopted Python on a number of large-scale efforts due to these characteristics, as well as the remote nature of near-term quantum processors. The use of Python has enabled quick prototyping for quantum code that directly benefits pertinent research and development efforts in quantum scientific computing. However, this rapid prototyping ability comes at the cost of future performant integration for tightly-coupled CPU-QPU architectures with fast-feedback. Here we present a language extension to Python that enables heterogeneous quantum-classical computing via a robust C++ infrastructure for quantum just-in-time (QJIT) compilation. Our work builds off the QCOR C++ language extension and compiler infrastructure to enable a single-source, quantum hardware-agnostic approach to quantum-classical computing that retains the performance required for tightly coupled CPU-QPU compute models. We detail this Pythonic extension, its programming model and underlying software architecture, and provide a robust set of examples to demonstrate the utility of our approach.


title: New Integral Estimates in Substatic Riemannian Manifolds and the Alexandrov Theorem

abstract: We derive new integral estimates on substatic manifolds with boundary of horizon type, naturally arising in General Relativity. In particular, we generalize to this setting an identity due to Magnanini-Poggesi leading to the Alexandrov Theorem in the Euclidean space and improve on a Heintze-Karcher type inequality due to Li-Xia. Our method relies on the introduction of a new vector field with nonnegative divergence, generalizing to this setting the P-function technique of Weinberger.


title: Doubled patterns with reversal and square-free doubled patterns

abstract: In combinatorics on words, a word $w$ over an alphabet $\Sigma$ is said to avoid a pattern $p$ over an alphabet $\Delta$ if there is no factor $f$ of $w$ such that $f=h(p)$ where $h:\Delta^*\to\Sigma^*$ is a non-erasing morphism. A pattern $p$ is said to be $k$-avoidable if there exists an infinite word over a $k$-letter alphabet that avoids $p$. A pattern is \emph{doubled} if every variable occurs at least twice. Doubled patterns are known to be $3$-avoidable. Currie, Mol, and Rampersad have considered a generalized notion which allows variable occurrences to be reversed. That is, $h(V^R)$ is the mirror image of $h(V)$ for every $V\in\Delta$. We show that doubled patterns with reversal are $3$-avoidable. We also conjecture that (classical) doubled patterns that do not contain a square are $2$-avoidable. We confirm this conjecture for patterns with at most 4 variables. This implies that for every doubled pattern $p$, the growth rate of ternary words avoiding $p$ is at least the growth rate of ternary square-free words. A previous version of this paper containing only the first result has been presented at WORDS 2021.


title: What shall we do with an hour of data? Speech recognition for the un- and under-served languages of Common Voice

abstract: This technical report describes the methods and results of a three-week sprint to produce deployable speech recognition models for 31 under-served languages of the Common Voice project. We outline the preprocessing steps, hyperparameter selection, and resulting accuracy on official testing sets. In addition to this we evaluate the models on multiple tasks: closed-vocabulary speech recognition, pre-transcription, forced alignment, and key-word spotting. The following experiments use Coqui STT, a toolkit for training and deployment of neural Speech-to-Text models.


title: Closed geodesics and Fr{\o}yshov invariants of hyperbolic three-manifolds

abstract: Froyshov invariants are numerical invariants of rational homology three-spheres derived from gradings in monopole Floer homology. In the past few years, they have been employed to solve a wide range of problems in three and four-dimensional topology. In this paper, we look at connections with hyperbolic geometry for the class of minimal $L$-spaces. In particular, we study relations between Froyshov invariants and closed geodesics using ideas from analytic number theory. We discuss two main applications of our approach. First, we derive effective upper bounds for the Froyshov invariants of minimal hyperbolic $L$-spaces purely in terms of volume and injectivity radius. Second, we describe an algorithm to compute Froyshov invariants of minimal $L$-spaces in terms of data arising from hyperbolic geometry. As a concrete example of our method, we compute the Froyshov invariants for all spin$^c$ structures on the Seifert-Weber dodecahedral space. Along the way, we also prove several results about the eta invariants of the odd signature and Dirac operators on hyperbolic three-manifolds which might be of independent interest.


title: Some inequalities and applications of Simons' type formulas in Riemannian, affine and statistical geometry

abstract: A few formulas and theorems for statistical structures are proved. They deal with various curvatures as well as with metric properties of the cubic form or its covariant derivative. Some of them generalize formulas and theorems known in the case of Lagrangian submanifolds or affine hypersurfaces.


title: Continuous wave laser thermal restoration of oxidized lead-based pigments in mural paintings

abstract: Red lead and lead white are some of the most ancient and common pigments in mural paintings. However, they tend to blacken with time due to their oxidation to plattnerite (\b{eta}-PbO2). The possibility to induce the reconversion reactions by CW laser heating is hereby discussed. A thermodynamic study by TGA showed that direct cerussite or hydrocerussite formation from plattnerite are not suitable reconversion routes, which was confirmed by laser irradiation trials under CO2 and CO2/H2O fluxes. Minium (Pb3O4) and subsequent massicot (\b{eta}-PbO) formation from plattnerite were achieved (confirmed by SEM-EDS, XRD and micro-Raman) under Ar+, 810 nm diode and Nd:YAG lasers. The latter appears to be the most suited for restauration purposes, given the broad minium reconversion irradiance range. This is confirmed by successful trials on macroscopic areas of naturally darkened red lead containing samples.


title: Sample selection for efficient image annotation

abstract: Supervised object detection has been proven to be successful in many benchmark datasets achieving human-level performances. However, acquiring a large amount of labeled image samples for supervised detection training is tedious, time-consuming, and costly. In this paper, we propose an efficient image selection approach that samples the most informative images from the unlabeled dataset and utilizes human-machine collaboration in an iterative train-annotate loop. Image features are extracted by the CNN network followed by the similarity score calculation, Euclidean distance. Unlabeled images are then sampled into different approaches based on the similarity score. The proposed approach is straightforward, simple and sampling takes place prior to the network training. Experiments on datasets show that our method can reduce up to 80% of manual annotation workload, compared to full manual labeling setting, and performs better than random sampling.


title: Photonic modal circulator using temporal refractive-index modulation with spatial inversion symmetry

abstract: It has been demonstrated that dynamic refractive index modulation, which breaks time-reversal symmetry, can be used to create on-chip non-reciprocal photonic devices. In order to achieve amplitude non-reciprocity, all such devices moreover require modulations that break spatial symmetries, which adds complexity in implementations. Here we introduce a modal circulator, which achieves amplitude non-reciprocity through a circulation motion among three modes. We show that such a circulator can be achieved in a dynamically-modulated structure that preserves mirror symmetry, and as a result can be implemented using only a single standing-wave modulator, which significantly simplifies the implementation of dynamically-modulated non-reciprocal device. We also prove that in terms of the number of modes involved in the transport process, the modal circulator represents the minimum configuration in which complete amplitude non-reciprocity can be achieved while preserving spatial symmetry.


title: Gauge invariants of linearized gravity with a general background metric

abstract: In linearized gravity with distributed matter, the background metric has no generic symmetries, and decomposition of the metric perturbation into global normal modes is generally impractical. This complicates the identification of the gauge-invariant part of the perturbation, which is a concern, for example, in the theory of dispersive gravitational waves whose energy--momentum must be gauge-invariant. Here, we propose how to identify the gauge-invariant part of the metric perturbation and the six independent gauge invariants \textit{per~se} for an arbitrary background metric. For the Minkowski background, the operator that projects the metric perturbation on the invariant subspace is proportional to the well-known dispersion operator of linear gravitational waves in~vacuum. For a general background, this operator is expressed in terms of the Green's operator of the vacuum wave equation. If the background is smooth, it can be found asymptotically using the inverse scale of the background metric as a small parameter.


title: An ALMA Survey of Chemistry in Disks around M4-M5 Stars

abstract: M-stars are the most common hosts of planetary systems in the Galaxy. Protoplanetary disks around M-stars thus offer a prime opportunity to study the chemistry of planet-forming environments. We present an ALMA survey of molecular line emission toward a sample of five protoplanetary disks around M4-M5 stars (FP Tau, J0432+1827, J1100-7619, J1545-3417, and Sz 69). These observations can resolve chemical structures down to tens of AU. Molecular lines of $^{12}$CO, $^{13}$CO, C$^{18}$O, C$_2$H, and HCN are detected toward all five disks. Lines of H$_2$CO and DCN are detected toward 2/5 and 1/5 disks, respectively. For disks with resolved C$^{18}$O, C$_2$H, HCN, and H$_2$CO emission, we observe substructures similar to those previously found in disks around solar-type stars (e.g., rings, holes, and plateaus). C$_2$H and HCN excitation conditions estimated interior to the pebble disk edge for the bright disk J1100-7619 are consistent with previous measurements around solar-type stars. The correlation previously found between C$_2$H and HCN fluxes for solar-type disks extends to our M4-M5 disk sample, but the typical C$_2$H/HCN ratio is higher for the M4-M5 disk sample. This latter finding is reminiscent of the hydrocarbon enhancements found by previous observational infrared surveys in the innermost ($<$10AU) regions of M-star disks, which is intriguing since our disk-averaged fluxes are heavily influenced by flux levels in the outermost disk, exterior to the pebble disk edge. Overall, most of the observable chemistry at 10-100AU appears similar for solar-type and M4-M5 disks, but hydrocarbons may be more abundant around the cooler stars.


title: Value Iteration in Continuous Actions, States and Time

abstract: Classical value iteration approaches are not applicable to environments with continuous states and actions. For such environments, the states and actions are usually discretized, which leads to an exponential increase in computational complexity. In this paper, we propose continuous fitted value iteration (cFVI). This algorithm enables dynamic programming for continuous states and actions with a known dynamics model. Leveraging the continuous-time formulation, the optimal policy can be derived for non-linear control-affine dynamics. This closed-form solution enables the efficient extension of value iteration to continuous environments. We show in non-linear control experiments that the dynamic programming solution obtains the same quantitative performance as deep reinforcement learning methods in simulation but excels when transferred to the physical system. The policy obtained by cFVI is more robust to changes in the dynamics despite using only a deterministic model and without explicitly incorporating robustness in the optimization. Videos of the physical system are available at \url{https://sites.google.com/view/value-iteration}.


title: Deep Bandits Show-Off: Simple and Efficient Exploration with Deep Networks

abstract: Designing efficient exploration is central to Reinforcement Learning due to the fundamental problem posed by the exploration-exploitation dilemma. Bayesian exploration strategies like Thompson Sampling resolve this trade-off in a principled way by modeling and updating the distribution of the parameters of the action-value function, the outcome model of the environment. However, this technique becomes infeasible for complex environments due to the computational intractability of maintaining probability distributions over parameters of outcome models of corresponding complexity. Moreover, the approximation techniques introduced to mitigate this issue typically result in poor exploration-exploitation trade-offs, as observed in the case of deep neural network models with approximate posterior methods that have been shown to underperform in the deep bandit scenario. In this paper we introduce Sample Average Uncertainty (SAU), a simple and efficient uncertainty measure for contextual bandits. While Bayesian approaches like Thompson Sampling estimate outcomes uncertainty indirectly by first quantifying the variability over the parameters of the outcome model, SAU is a frequentist approach that directly estimates the uncertainty of the outcomes based on the value predictions. Importantly, we show theoretically that the uncertainty measure estimated by SAU asymptotically matches the uncertainty provided by Thompson Sampling, as well as its regret bounds. Because of its simplicity SAU can be seamlessly applied to deep contextual bandits as a very scalable drop-in replacement for epsilon-greedy exploration. We confirm empirically our theory by showing that SAU-based exploration outperforms current state-of-the-art deep Bayesian bandit methods on several real-world datasets at modest computation cost. Code is available at \url{https://github.com/ibm/sau-explore}.


title: An automatic system to detect equivalence between iterative algorithms

abstract: When are two algorithms the same? How can we be sure a recently proposed algorithm is novel, and not a minor twist on an existing method? In this paper, we present a framework for reasoning about equivalence between a broad class of iterative algorithms, with a focus on algorithms designed for convex optimization. We propose several notions of what it means for two algorithms to be equivalent, and provide computationally tractable means to detect equivalence. Our main definition, oracle equivalence, states that two algorithms are equivalent if they result in the same sequence of calls to the function oracles (for suitable initialization). Borrowing from control theory, we use state-space realizations to represent algorithms and characterize algorithm equivalence via transfer functions. Our framework can also identify and characterize some algorithm transformations including permutations of the update equations, repetition of the iteration, and conjugation of some of the function oracles in the algorithm. To support the paper, we have developed a software package named Linnaeus that implements the framework to identify other iterative algorithms that are equivalent to an input algorithm. More broadly, this framework and software advances the goal of making mathematics searchable.


title: Large deviation principles induced by the Stiefel manifold, and random multi-dimensional projections

abstract: Given an $n$-dimensional random vector $X^{(n)}$ , for $k < n$, consider its $k$-dimensional projection $\mathbf{a}_{n,k}X^{(n)}$, where $\mathbf{a}_{n,k}$ is an $n \times k$-dimensional matrix belonging to the Stiefel manifold $\mathbb{V}_{n,k}$ of orthonormal $k$-frames in $\mathbb{R}^n$. For a class of sequences $\{X^{(n)}\}$ that includes the uniform distributions on scaled $\ell_p^n$ balls, $p \in (1,\infty]$, and product measures with sufficiently light tails, it is shown that the sequence of projected vectors $\{\mathbf{a}_{n,k}^\intercal X^{(n)}\}$ satisfies a large deviation principle whenever the empirical measures of the rows of $\sqrt{n} \mathbf{a}_{n,k}$ converge, as $n \rightarrow \infty$, to a probability measure on $\mathbb{R}^k$. In particular, when $\mathbf{A}_{n,k}$ is a random matrix drawn from the Haar measure on $\mathbb{V}_{n,k}$, this is shown to imply a large deviation principle for the sequence of random projections $\{\mathbf{A}_{n,k}^\intercal X^{(n)}\}$ in the quenched sense (that is, conditioned on almost sure realizations of $\{\mathbf{A}_{n,k}\}$). Moreover, a variational formula is obtained for the rate function of the large deviation principle for the annealed projections $\{\mathbf{A}_{n,k}^\intercal X^{(n)}\}$, which is expressed in terms of a family of quenched rate functions and a modified entropy term. A key step in this analysis is a large deviation principle for the sequence of empirical measures of rows of $\sqrt{n} \mathbf{A}_{n,k}$, which may be of independent interest. The study of multi-dimensional random projections of high-dimensional measures is of interest in asymptotic functional analysis, convex geometry and statistics. Prior results on quenched large deviations for random projections of $\ell_p^n$ balls have been essentially restricted to the one-dimensional setting.


title: Stability of pulsatile quasi-two-dimensional duct flows under a transverse magnetic field

abstract: This manuscript has been accepted for publication in Physical Review Fluids, see https://journals.aps.org/prfluids/accepted/53075Se8O0b1b109b1cc0061b280aaa122f0f92dc. The stability of a pulsatile quasi-two-dimensional duct flow was numerically investigated. Flow was driven, in concert, by a constant pressure gradient and by the synchronous oscillation of the lateral walls. This prototypical setup serves to aid understanding of unsteady magnetohydrodynamic flows in liquid metal coolant ducts subjected to transverse magnetic fields, motivated by the conditions expected in magnetic confinement fusion reactors. A wide range of wall oscillation frequencies and amplitudes were simulated. Focus was placed on the driving pulsation optimized for the greatest reduction in the critical Reynolds number, for a range of friction parameters $H$ (proportional to magnetic field strength). An almost $70$% reduction in the critical Reynolds number, relative to that for the steady base flow, was obtained toward the hydrodynamic limit ($H=10^{-7}$), while just over a $90$% reduction was obtained by $H=10$. For all oscillation amplitudes, increasing $H$ consistently led to an increasing percentage reduction in the critical Reynolds number. This is a promising result, given fusion relevant conditions of $H \geq 10^4$. These reductions were obtained by selecting a frequency that both ensures prominent inflection points, and a growth in perturbation energy in phase with deceleration of the base flow. Nonlinear simulations at the optimized frequency and amplitude still satisfied the no net growth condition at the greatly reduced critical Reynolds numbers. However, although the linear mode undergoes a symmetry breaking process, turbulence was not triggered. Nonlinear base flow modulation also arrested the linear decay of the perturbation, with exponential growth not observed at supercritical Reynolds numbers.


title: Gravity with a generalized conformal scalar field: Theory and solutions

abstract: We naturally extend the theory of gravity with a conformally coupled scalar field by only requiring conformal invariance of the scalar field equation of motion and not of the action. The classically extended theory incorporates a scalar-Gauss-Bonnet sector and has second-order equations of motion, belonging to the Horndeski class. Remarkably, the theory features a purely geometrical field equation that allows for closed-form black hole solutions and cosmologies to be easily found. These solutions permit investigations of in-vogue scalar-Gauss-Bonnet corrections to the gravitational action without the need of resorting to approximations or numerical methods. We discuss on the connection to the recently formulated 4D Einstein-Gauss-Bonnet theory of gravity.


title: Assessing the Syntactic Capabilities of Transformer-based Multilingual Language Models

abstract: Multilingual Transformer-based language models, usually pretrained on more than 100 languages, have been shown to achieve outstanding results in a wide range of cross-lingual transfer tasks. However, it remains unknown whether the optimization for different languages conditions the capacity of the models to generalize over syntactic structures, and how languages with syntactic phenomena of different complexity are affected. In this work, we explore the syntactic generalization capabilities of the monolingual and multilingual versions of BERT and RoBERTa. More specifically, we evaluate the syntactic generalization potential of the models on English and Spanish tests, comparing the syntactic abilities of monolingual and multilingual models on the same language (English), and of multilingual models on two different languages (English and Spanish). For English, we use the available SyntaxGym test suite; for Spanish, we introduce SyntaxGymES, a novel ensemble of targeted syntactic tests in Spanish, designed to evaluate the syntactic generalization capabilities of language models through the SyntaxGym online platform.


title: Joule-Thomson expansion for quantum corrected AdS-Reissner-Nordstrom black holes in Kiselev spacetime

abstract: In this work we study the inversion temperature associated with the Joule-Thomson expansion from the thermodynamics of AdS-Reissner-N\"ordstrom black holes. We include quantum corrections in a cosmological fluid that can describe phantom dark matter or quintessence, both in a Kiselev scenario. The description of such physical systems involves numerical solutions and the results are presented as temperature-pressure plots for various values of the parameters of our model. We find non-zero minimum inversion temperatures as well as non-zero minimum pressures depending on the values of those parameters. Completing our study, we also find isenthalpic curves associated with black hole fixed mass processes.


title: Automated quantitative analysis of first-pass myocardial perfusion magnetic resonance imaging data

abstract: Coronary artery disease (CAD) remains the world's leading cause of mortality and the disease burden is continually expanding as the population ages. Recently, the MR-INFORM randomised trial has demonstrated that the management of patients with stable CAD can be guided by stress perfusion cardiovascular magnetic resonance (CMR) imaging and it is non-inferior to the using the invasive reference standard of fractional flow reserve. The benefits of using stress perfusion CMR include that it is non-invasive and significantly reduces the number of unnecessary coronary revascularisations. As compared to other ischaemia tests, it boasts a high spatial resolution and does not expose the patient to ionising radiation. However, the main limitation of stress perfusion CMR is that the diagnostic accuracy is highly dependent on the level of training of the operator, resulting in the test only being performed routinely in experienced tertiary centres. The clinical translation of stress perfusion CMR would be greatly aided by a fully automated, user-independent, quantitative evaluation of myocardial blood flow. This thesis presents major steps towards this goal: robust motion correction, automated image processing, reliable quantitative modelling, and thorough validation. The motion correction scheme makes use of data decomposition techniques, such as robust principal component analysis, to mitigate the difficulties in image registration caused by the dynamic contrast enhancement. The motion corrected image series are input to a processing pipeline which leverages the recent advances in image processing facilitated by deep learning. The pipeline utilises convolutional neural networks to perform a series of computer vision tasks including myocardial segmentation and right ventricular insertion point detection...


title: Steps minimize dissipation in rapidly driven stochastic systems

abstract: Micro- and nano-scale systems driven by rapid changes in control parameters (control protocols) dissipate significant energy. In the fast-protocol limit, we find that protocols that minimize dissipation at fixed duration are universally given by a two-step process, jumping to and from a point that balances jump size with fast relaxation. Jump protocols could be exploited by molecular machines or thermodynamic computing to improve energetic efficiency, and implemented in nonequilibrium free-energy estimation to improve accuracy.


title: Budget-Constrained Coalition Strategies with Discounting

abstract: Discounting future costs and rewards is a common practice in accounting, game theory, and machine learning. In spite of this, existing logics for reasoning about strategies with cost and resource constraints do not account for discounting. The paper proposes a sound and complete logical system for reasoning about budget-constrained strategic abilities that incorporates discounting into its semantics.


title: Emergence of Structural Bias in Differential Evolution

abstract: Heuristic optimisation algorithms are in high demand due to the overwhelming amount of complex optimisation problems that need to be solved. The complexity of these problems is well beyond the boundaries of applicability of exact optimisation algorithms and therefore require modern heuristics to find feasible solutions quickly. These heuristics and their effects are almost always evaluated and explained by particular problem instances. In previous works, it has been shown that many such algorithms show structural bias, by either being attracted to a certain region of the search space or by consistently avoiding regions of the search space, on a special test function designed to ensure uniform 'exploration' of the domain. In this paper, we analyse the emergence of such structural bias for Differential Evolution (DE) configurations and, specifically, the effect of different mutation, crossover and correction strategies. We also analyse the emergence of the structural bias during the run-time of each algorithm. We conclude with recommendations of which configurations should be avoided in order to run DE unbiased.


title: Chaotic switching in driven-dissipative Bose-Hubbard dimers: when a flip bifurcation meets a T-point in $R^4$

abstract: The Bose--Hubbard dimer model is a celebrated fundamental quantum mechanical model that accounts for the dynamics of bosons at two interacting sites. It has been realized experimentally by two coupled, driven and lossy photonic crystal nanocavities, which are optical devices that operate with only a few hundred photons due to their extremely small size. Our work focuses on characterizing the different dynamics that arise in the semiclassical approximation of such driven-dissipative photonic Bose--Hubbard dimers. Mathematically, this system is a four-dimensional autonomous vector field that describes two specific coupled oscillators, where both the amplitude and the phase are important. We perform a bifurcation analysis of this system to identify regions of different behavior as the pump power $f$ and the detuning $\delta$ of the driving signal are varied, for the case of fixed positive coupling. The bifurcation diagram in the $(f,\delta)$-plane is organized by two points of codimension-two bifurcations -- a $Z_2$-equivariant homoclinic flip bifurcation and a Bykov T-point -- and provides a roadmap for the observable dynamics, including different types of chaotic behavior. To illustrate the overall structure and different accumulation processes of bifurcation curves and associated regions, our bifurcation analysis is complemented by the computation of kneading invariants and of maximum Lyapunov exponents in the $(f,\delta)$-plane. The bifurcation diagram displays a menagerie of dynamical behavior and offers insights into the theory of global bifurcations in a four-dimensional phase space, including novel bifurcation phenomena such as degenerate singular heteroclinic cycles.


title: Known cases of the Hodge conjecture

abstract: The Hodge conjecture is a major open problem in complex algebraic geometry. In this survey, we discuss the main cases where the conjecture is known, and also explain an approach by Griffiths-Green to solve the problem.


title: Antiferromagnetic spatial photonic Ising machine through optoelectronic correlation computing

abstract: Recently, spatial photonic Ising machines (SPIM) have been demonstrated to compute the minima of Hamiltonians for large-scale spin systems. Here we propose to implement an antiferromagnetic model through optoelectronic correlation computing with SPIM. Also we exploit the gauge transformation which enables encoding the spins and the interaction strengths in a single phase-only spatial light modulator. With a simple setup, we experimentally show the ground state search of an antiferromagnetic model with $40000$ spins in number-partitioning problem. Thus such an optoelectronic computing exhibits great programmability and scalability for the practical applications of studying statistical systems and combinatorial optimization problems.

categories: cs.ET cond-mat.dis-nn physics.app-ph physics.optics


title: Correlation-Robust Optimal Auctions

abstract: I study the design of auctions in which the auctioneer is assumed to have information only about the marginal distribution of a generic bidder's valuation, but does not know the correlation structure of the joint distribution of bidders' valuations. I assume that a generic bidder's valuation is bounded and $\bar{v}$ is the maximum valuation of a generic bidder. The performance of a mechanism is evaluated in the worst case over the uncertainty of joint distributions that are consistent with the marginal distribution. For the two-bidder case, the second-price auction with the uniformly distributed random reserve maximizes the worst-case expected revenue across all dominant-strategy mechanisms under certain regularity conditions. For the $N$-bidder ($N\ge3$) case, the second-price auction with the $\bar{v}-$scaled $Beta (\frac{1}{N-1},1)$ distributed random reserve maximizes the worst-case expected revenue across standard (a bidder whose bid is not the highest will never be allocated) dominant-strategy mechanisms under certain regularity conditions. When the probability mass condition (part of the regularity conditions) does not hold, the second-price auction with the $s^*-$scaled $Beta (\frac{1}{N-1},1)$ distributed random reserve maximizes the worst-case expected revenue across standard dominant-strategy mechanisms, where $s^*\in (0,\bar{v})$.


title: The Coupling of Galactic Dark Matter Halos with Stellar Bars

abstract: Resonant torques couple stellar bars to dark matter halos. Here we use high-resolution numerical simulations to demonstrate long-term angular momentum transfer between stellar bars and dark matter orbits of varying orientation. We show that bar-driven reversals of dark matter orbit orientations can play a surprisingly large role in the evolution of the bar pattern speed. In predominantly prograde (co-rotating) halos, dark matter orbits become trapped in the stellar bar forming a parallel dark matter bar. This dark matter bar reaches more than double the vertical height of the stellar bar. In halos dominated by retrograde orbits, a dark matter wake forms oriented perpendicular to the stellar bar. These dark matter over-densities provide a novel space to look for dark matter annihilation or decay signals. % We predict that the Milky Way hosts a dark matter bar aligned with the stellar bar as well as a dark matter wake the near-side of which should extend from Galactic center to a galactic longitude of $l \approx 323^\circ$.


title: Adaptive Policy Transfer in Reinforcement Learning

abstract: Efficient and robust policy transfer remains a key challenge for reinforcement learning to become viable for real-wold robotics. Policy transfer through warm initialization, imitation, or interacting over a large set of agents with randomized instances, have been commonly applied to solve a variety of Reinforcement Learning tasks. However, this seems far from how skill transfer happens in the biological world: Humans and animals are able to quickly adapt the learned behaviors between similar tasks and learn new skills when presented with new situations. Here we seek to answer the question: Will learning to combine adaptation and exploration lead to a more efficient transfer of policies between domains? We introduce a principled mechanism that can "Adapt-to-Learn", that is adapt the source policy to learn to solve a target task with significant transition differences and uncertainties. We show that the presented method learns to seamlessly combine learning from adaptation and exploration and leads to a robust policy transfer algorithm with significantly reduced sample complexity in transferring skills between related tasks.


title: Ultrafast Distributed Coloring of High Degree Graphs

abstract: We give a new randomized distributed algorithm for the $\Delta+1$-list coloring problem. The algorithm and its analysis dramatically simplify the previous best result known of Chang, Li, and Pettie [SICOMP 2020]. This allows for numerous refinements, and in particular, we can color all $n$-node graphs of maximum degree $\Delta \ge \log^{2+\Omega(1)} n$ in $O(\log^* n)$ rounds. The algorithm works in the CONGEST model, i.e., it uses only $O(\log n)$ bits per message for communication. On low-degree graphs, the algorithm shatters the graph into components of size $\operatorname{poly}(\log n)$ in $O(\log^* \Delta)$ rounds, showing that the randomized complexity of $\Delta+1$-list coloring in CONGEST depends inherently on the deterministic complexity of related coloring problems.


title: Service Function Chaining in MEC: A Mean-Field Game and Reinforcement Learning Approach

abstract: Multi-access edge computing (MEC) and network virtualization technologies are important enablers for fifth-generation (5G) networks to deliver diverse applications and services. Services are often provided as fully connected virtual network functions (VNF)s, through service function chaining (SFC). However, the problem of allocating SFC resources at the network edge still faces many challenges related to the way VNFs are placed, chained and scheduled. In this paper, to solve these problems, we propose a game theory-based approach with the objective to reduce service latency in the context of SFC at the network edge. The problem of allocating SFC resources can be divided into two subproblems. 1) The VNF placement and routing subproblem, and 2) the VNF scheduling subproblem. For the former subproblem, we formulate it as a mean field game (MFG) in which VNFs are modeled as entities contending over edge resources with the goal of reducing the resource consumption of MEC nodes and reducing latency for users. We propose a on a reinforcement learning-based technique, where the Ishikawa-Mann learning algorithm (IMLA) is used. For the later subproblem we formulate it as a matching game between the VFNs and an edge resources in order to find the execution order of the VNFs while reducing the latency. To efficiently solve it, we propose a modified version of the many-to-one deferred acceptance algorithm (DAA), called the enhanced multi-step deferred acceptance algorithm (eMSDA). To illustrate the performance of the proposed approaches, we perform extensive simulations. The obtained results show that the proposed approaches outperform the benchmarks other state-of-the-art methods.


title: ppsim: A software package for efficiently simulating and visualizing population protocols

abstract: We introduce ppsim, a software package for efficiently simulating population protocols, a widely-studied subclass of chemical reaction networks (CRNs) in which all reactions have two reactants and two products. Each step in the dynamics involves picking a uniform random pair from a population of $n$ molecules to collide and have a (potentially null) reaction. In a recent breakthrough, Berenbrink, Hammer, Kaaser, Meyer, Penschuck, and Tran [ESA 2020] discovered a population protocol simulation algorithm quadratically faster than the naive algorithm, simulating $\Theta(\sqrt{n})$ reactions in *constant* time (independently of $n$, though the time scales with the number of species), while preserving the *exact* stochastic dynamics. ppsim implements this algorithm, with a tightly optimized Cython implementation that can exactly simulate hundreds of billions of reactions in seconds. It dynamically switches to the CRN Gillespie algorithm for efficiency gains when the number of applicable reactions in a configuration becomes small. As a Python library, ppsim also includes many useful tools for data visualization in Jupyter notebooks, allowing robust visualization of time dynamics such as histogram plots at time snapshots and averaging repeated trials. Finally, we give a framework that takes any CRN with only bimolecular (2 reactant, 2 product) or unimolecular (1 reactant, 1 product) reactions, with arbitrary rate constants, and compiles it into a continuous-time population protocol. This lets ppsim exactly sample from the chemical master equation (unlike approximate heuristics such as tau-leaping or LNA), while achieving asymptotic gains in running time. In linked Jupyter notebooks, we demonstrate the efficacy of the tool on some protocols of interest in molecular programming, including the approximate majority CRN and CRN models of DNA strand displacement reactions.


title: Friends and Foes: Conditional Occurrence Rates of Exoplanet Companions and their Impact on Radial Velocity Follow-up Surveys

abstract: Population studies of Kepler's multi-planet systems have revealed a surprising degree of structure in their underlying architectures. Information from a detected transiting planet can be combined with a population model to make predictions about the presence and properties of additional planets in the system. Using a statistical model for the distribution of planetary systems (He et al. 2020; arXiv:2007.14473), we compute the conditional occurrence of planets as a function of the period and radius of Kepler-detectable planets. About half ($0.52 \pm 0.03$) of the time, the detected planet is not the planet with the largest semi-amplitude $K$ in the system, so efforts to measure the mass of the transiting planet with radial velocity (RV) follow-up will have to contend with additional planetary signals in the data. We simulate RV observations to show that assuming a single-planet model to measure the $K$ of the transiting planet often requires significantly more observations than in the ideal case with no additional planets, due to the systematic errors from unseen planet companions. Our results show that planets around 10-day periods with $K$ close to the single-measurement RV precision ($\sigma_{1,\rm obs}$) typically require $\sim 100$ observations to measure their $K$ to within 20% error. For a next generation RV instrument achieving $\sigma_{1,\rm obs} = 10$ cm/s, about $\sim 200$ ($600$) observations are needed to measure the $K$ of a transiting Venus in a Kepler-like system to better than 20% (10%) error, which is $\sim 2.3$ times as many as what would be necessary for a Venus without any planetary companions.


title: Lecture notes on descriptional complexity and randomness

abstract: A didactical survey of the foundations of Algorithmic Information Theory. These notes are short on motivation, history and background but introduce some of the main techniques and concepts of the field. The "manuscript" has been evolving over the years. Please, look at "Version history" below to see what has changed when.


title: Crowding breaks the forward/backward symmetry of transition times in biased random walks

abstract: Microscopic mechanisms of natural processes are frequently understood in terms of random walk models by analyzing local particle transitions. This is because these models properly account for dynamic processes at the molecular level and provide a clear physical picture. Recent theoretical studies made a surprising discovery that in complex systems the symmetry of molecular forward/backward transition times with respect to local bias in the dynamics may be broken and it may take longer to go downhill than uphill. The physical origins of these phenomena remain not fully understood. Here we explore in more detail the microscopic features of the symmetry breaking in the forward/backward transition times by analyzing exactly solvable discrete-state stochastic models. In particular, we consider a specific case of two random walkers on four-sites periodic lattice as the way to represent the general systems with multiple pathways. It is found that the asymmetry in transition times depends on several factors that include the degree of deviation from equilibrium, the particle crowding and methods of measurements of dynamic properties. Our theoretical analysis suggests that the asymmetry in transition times can be explored experimentally for determining the important microscopic features of natural processes by quantitative measuring the local deviations from equilibrium and the degrees of crowding.


title: Latest results on quiescent and post-disruption runaway electron mitigation experiments at Frascati Tokamak Upgrade

abstract: Results from the last FTU campaigns on the deuterium large (wrt FTU volume) pellet REs suppression capability, mainly due to the induced burst MHD activity expelling REs seed are presented for discharges with 0.5 MA and 5.3T. Clear indications of avalanche multiplication of REs following single pellet injection on 0.36 MA flat-top discharges is shown together with quantitative indications of dissipative effects in terms of critical electrical field increase due to fan-like instabilities. Analysis of large fan-like instabilities on post-disruption RE beams, that seem to be correlated with low electrical field and background density drops, reveal their strong RE energy suppression capability suggesting a new strategy for RE energy suppression controlling large fan instabilities. We demonstrate how such density drops can be induced using modulated ECRH power on post-disruption beams.


title: Accountable Error Characterization

abstract: Customers of machine learning systems demand accountability from the companies employing these algorithms for various prediction tasks. Accountability requires understanding of system limit and condition of erroneous predictions, as customers are often interested in understanding the incorrect predictions, and model developers are absorbed in finding methods that can be used to get incremental improvements to an existing system. Therefore, we propose an accountable error characterization method, AEC, to understand when and where errors occur within the existing black-box models. AEC, as constructed with human-understandable linguistic features, allows the model developers to automatically identify the main sources of errors for a given classification system. It can also be used to sample for the set of most informative input points for a next round of training. We perform error detection for a sentiment analysis task using AEC as a case study. Our results on the sample sentiment task show that AEC is able to characterize erroneous predictions into human understandable categories and also achieves promising results on selecting erroneous samples when compared with the uncertainty-based sampling.


title: Word-level Human Interpretable Scoring Mechanism for Novel Text Detection Using Tsetlin Machines

abstract: Recent research in novelty detection focuses mainly on document-level classification, employing deep neural networks (DNN). However, the black-box nature of DNNs makes it difficult to extract an exact explanation of why a document is considered novel. In addition, dealing with novelty at the word-level is crucial to provide a more fine-grained analysis than what is available at the document level. In this work, we propose a Tsetlin machine (TM)-based architecture for scoring individual words according to their contribution to novelty. Our approach encodes a description of the novel documents using the linguistic patterns captured by TM clauses. We then adopt this description to measure how much a word contributes to making documents novel. Our experimental results demonstrate how our approach breaks down novelty into interpretable phrases, successfully measuring novelty.


title: Personalized Popular Music Generation Using Imitation and Structure

abstract: Many practices have been presented in music generation recently. While stylistic music generation using deep learning techniques has became the main stream, these models still struggle to generate music with high musicality, different levels of music structure, and controllability. In addition, more application scenarios such as music therapy require imitating more specific musical styles from a few given music examples, rather than capturing the overall genre style of a large data corpus. To address requirements that challenge current deep learning methods, we propose a statistical machine learning model that is able to capture and imitate the structure, melody, chord, and bass style from a given example seed song. An evaluation using 10 pop songs shows that our new representations and methods are able to create high-quality stylistic music that is similar to a given input song. We also discuss potential uses of our approach in music evaluation and music therapy.


title: Collision Between Molecular Clouds-III: The effects of cloud initial density profile on head-on collisions

abstract: In this third paper in our cloud-collision series, we present the results from simulations of head-on collisions with a strongly centrally-condensed initial density profile of $\rm \rho \propto R^{-2}$. We investigate the impact of these density profiles on the overall evolution of the simulations: the structures formed, their dynamical evolution, and their star formation activity. We consider clouds which are globally bound and globally unbound, leading to three different scenarios -- the collision of a bound cloud with a bound cloud, the collision of two unbound clouds, or the collision of one cloud of each type. In all the simulations dense star clusters form before the collisions occur, and we find that star formation remains confined to these systems and is little affected by the collisions. If the clouds' are both initially bound, the collision forms a filamentary structure, but otherwise, this does not occur. We observe that rotating structures form around the clusters, but they also form in our non-colliding control simulations, so are not a consequence of the collisions. Dissipation of kinetic energy in these simulations is inefficient because of the substructure created in the clouds by turbulence before the collisions. As a result, although some gas is left bound in the COM frame, the star clusters formed in the simulations do not become bound to each other.


title: High-dimensional percolation criticality and hints of mean-field-like caging of the random Lorentz gas

abstract: The random Lorentz gas (RLG) is a minimal model for transport in disordered media. Despite the broad relevance of the model, theoretical grasp over its properties remains weak. For instance, the scaling with dimension d of its localization transition at the void percolation threshold is not well controlled analytically nor computationally. A recent study [Biroli et al. Phys. Rev. E L030104 (2021)] of the caging behavior of the RLG motivated by the mean-field theory of glasses has uncovered physical inconsistencies in that scaling that heighten the need for guidance. Here, we first extend analytical expectations for asymptotic high-d bounds on the void percolation threshold, and then computationally evaluate both the threshold and its criticality in various d. In high-d systems, we observe that the standard percolation physics is complemented by a dynamical slowdown of the tracer dynamics reminiscent of mean-field caging. A simple modification of the RLG is found to bring the interplay between percolation and mean-field-like caging down to d=3.


title: Near Neighbor Search via Efficient Average Distortion Embeddings

abstract: A recent series of papers by Andoni, Naor, Nikolov, Razenshteyn, and Waingarten (STOC 2018, FOCS 2018) has given approximate near neighbour search (NNS) data structures for a wide class of distance metrics, including all norms. In particular, these data structures achieve approximation on the order of $p$ for $\ell_p^d$ norms with space complexity nearly linear in the dataset size $n$ and polynomial in the dimension $d$, and query time sub-linear in $n$ and polynomial in $d$. The main shortcoming is the exponential in $d$ pre-processing time required for their construction. In this paper, we describe a more direct framework for constructing NNS data structures for general norms. More specifically, we show via an algorithmic reduction that an efficient NNS data structure for a given metric is implied by an efficient average distortion embedding of it into $\ell_1$ or into Euclidean space. In particular, the resulting data structures require only polynomial pre-processing time, as long as the embedding can be computed in polynomial time. As a concrete instantiation of this framework, we give an NNS data structure for $\ell_p$ with efficient pre-processing that matches the approximation factor, space and query complexity of the aforementioned data structure of Andoni et al. On the way, we resolve a question of Naor (Analysis and Geometry in Metric Spaces, 2014) and provide an explicit, efficiently computable embedding of $\ell_p$, for $p \ge 2$, into $\ell_2$ with (quadratic) average distortion on the order of $p$. We expect our approach to pave the way for constructing efficient NNS data structures for all norms.


title: Lessons Learned from Virtual Transitions of Undergraduate Student Support Programs due to the COVID-19 Pandemic

abstract: The COVID-19 global pandemic spurred sudden and dramatic changes in the way universities and research programs operate, often with negative consequences that have disproportionately affected minorities, women, and people with low income. These consequences compound the difficulties these groups faced prior to the pandemic | in particular, the increased obstacles marginalized students encounter while pursuing STEM career paths. The National Astronomy Consortium (NAC) was founded to support underrepresented physics and astronomy students by providing research opportunities and long-term mentoring. In response to the pandemic and nation-wide shutdowns, the NAC shifted both its summer research program and its annual fall conference to a fully-remote format. Here, we discuss the changes the NAC made to achieve a successful virtual transition, and provide temporary and long-term recommendations for programs seeking to provide support to marginalized undergraduate and graduate researchers moving forward.


title: Sample and Computation Redistribution for Efficient Face Detection

abstract: Although tremendous strides have been made in uncontrolled face detection, efficient face detection with a low computation cost as well as high precision remains an open challenge. In this paper, we point out that training data sampling and computation distribution strategies are the keys to efficient and accurate face detection. Motivated by these observations, we introduce two simple but effective methods (1) Sample Redistribution (SR), which augments training samples for the most needed stages, based on the statistics of benchmark datasets; and (2) Computation Redistribution (CR), which reallocates the computation between the backbone, neck and head of the model, based on a meticulously defined search methodology. Extensive experiments conducted on WIDER FACE demonstrate the state-of-the-art efficiency-accuracy trade-off for the proposed \scrfd family across a wide range of compute regimes. In particular, \scrfdf{34} outperforms the best competitor, TinaFace, by $3.86\%$ (AP at hard set) while being more than \emph{3$\times$ faster} on GPUs with VGA-resolution images. We also release our code to facilitate future research.


title: A gauged non-universal $U(1)_{X}$ model to study muon $g-2$ and $B$ meson anomalies

abstract: We study a non-universal $U(1)_{X}$ extension of the Standard Model with an extended scalar sector of two doublets and one singlet plus three additional exotic quarks and two exotic charged leptons on the fermionic sector. In order to obtain the observed fermion mass hierarchy, an additional $\mathbb{Z}_{2}$ discrete symmetry is imposed, where the heaviest fermions acquire their masses from three different scales determined by two Higgs doublets and one singlet, whereas the lightest fermions obtain their masses from effective operators up to dimension seven. From chiral anomalies cancellation, the model also includes heavy right-handed neutrinos which get massive via an inverse see-saw mechanism, reproducing the observed mass differences for the active neutrinos. We analyze phenomenological consequences of the model in view of the so-called flavour anomalies, namely the latest measurement made by Fermilab of the anomalous magnetic moment of the muon $g-2$, and additionally, the fit to semi-leptonic $B$ meson decays made by different flavour groups, dominated mainly by the LHCb 2020 data. We obtain that the model can explain the former at the $1\sigma$ level by means of contributions coming from charged $W^{+}$ bosons interacting with exotic Majorana neutrinos at one-loop level, with the $Z'$ boson contribution itself coming from the $U(1)_X$ symmetry being negligible. However, we find that the model is able to accommodate the $b\to s\ell^{+}\ell^{-}$ transitions associated to the $B$ meson anomalies only at the $2\sigma$ level via tree-level $Z'$ boson exchange, while simultaneously respecting various constraints from the recent $R_{K^{(*)}}$ measurements made by the LHCb, neutrino trident production and $B-\bar B$ oscillations.


title: $\mathbb{i}$CAS: Imposed Automatic Selection and Localization of Complete Active Spaces

abstract: It is shown that a prechosen set of occupied/virtual valence/core atomic/fragmental orbitals can be transformed to an equivalent set of localized occupied/virtual pre-molecular orbitals (pre-LMO), which can then be taken as probes to select the same number of maximally matching localized occupied/virtual Hartree-Fock (HF) or restricted open-shell Hartree-Fock (ROHF) molecular orbitals as the initial local orbitals spanning the desired complete active space (CAS). In each cycle of the self-consistent field (SCF) calculation, the CASSCF orbitals can be localized by means of the noniterative ``top-down least-change'' algorithm for localizing ROHF orbitals, such that the maximum matching between the orbitals of two adjacent iterations can readily be monitored, leading finally to converged localized CASSCF orbitals that overlap most the guess orbitals. Such an approach is to be dubbed as ``imposed CASSCF'' ($\mathbb{i}$CASSCF or simply $\mathbb{i}$CAS in short) for good reasons: (1) it has been assumed that only those electronic states that have largest projections onto the active space defined by the prechosen atomic/fragmental orbitals are to be targeted. This is certainly an imposed constraint but has wide applications in organic and transition metal chemistry where valence (or core) atomic/fragmental orbitals can readily be identified. (2) The selection of both initial and optimized local active orbitals is imposed from the very beginning by the pre-LMOs (which span the same space as the prechosen atomic/fragmental orbitals). Apart from the (imposed) automation and localization, $\mathbb{i}$CAS has an additional merit: the CAS is guaranteed to be the same for all geometries for the pre-LMOs do not change in character with geometry. Both organic molecules and transition metal complexes are taken as showcases to reveal the efficacy of $\mathbb{i}$CAS.


title: Refining the mass estimate for the intermediate-mass black hole candidate in NGC 3319

abstract: Recent X-ray observations by Jiang et al. have identified an active galactic nucleus (AGN) in the bulgeless spiral galaxy NGC 3319, located just $14.3\pm1.1\,$Mpc away, and suggest the presence of an intermediate-mass black hole (IMBH; $10^2\leq M_\bullet/\mathrm{M_{\odot}}\leq10^5$) if the Eddington ratios are as high as 3 to $3\times10^{-3}$. In an effort to refine the black hole mass for this (currently) rare class of object, we have explored multiple black hole mass scaling relations, such as those involving the (not previously used) velocity dispersion, logarithmic spiral-arm pitch angle, total galaxy stellar mass, nuclear star cluster mass, rotational velocity, and colour of NGC 3319, to obtain ten mass estimates, of differing accuracy. We have calculated a mass of $3.14_{-2.20}^{+7.02}\times10^4\,\mathrm{M_\odot}$, with a confidence of 84% that it is $\leq$$10^5\,\mathrm{M_\odot}$, based on the combined probability density function from seven of these individual estimates. Our conservative approach excluded two black hole mass estimates (via the nuclear star cluster mass, and the fundamental plane of black hole activity $\unicode{x2014}$ which only applies to black holes with low accretion rates) that were upper limits of $\sim$$10^5\,{\rm M}_{\odot}$, and it did not use the $M_\bullet\unicode{x2013}L_{\rm 2-10\,keV}$ relation's prediction of $\sim$$10^5\,{\rm M}_{\odot}$. This target provides an exceptional opportunity to study an IMBH in AGN mode and advance our demographic knowledge of black holes. Furthermore, we introduce our novel method of meta-analysis as a beneficial technique for identifying new IMBH candidates by quantifying the probability that a galaxy possesses an IMBH.


title: Economic analysis of tidal stream turbine arrays: a review

abstract: This tidal stream energy industry has to date been comprised of small demonstrator projects made up of one to a four turbines. However, there are currently plans to expand to commercially sized projects with tens of turbines or more. As the industry moves to large-scale arrays for the first time, there has been a push to develop tools to optimise the array design and help bring down the costs. This review investigates different methods of modelling the economic performance of tidal-stream arrays, for use within these optimisation tools. The different cost reduction pathways are discussed from costs falling as the global installed capacity increases, due to greater experience, improved power curves through larger-diameter higher-rated turbines, to economic efficiencies that can be found by moving to large-scale arrays. A literature review is conducted to establish the most appropriate input values for use in economic models. This includes finding a best case, worst case and typical values for costs and other related parameters. The information collated in this review can provide a useful steering for the many optimisation tools that have been developed, especially when cost information is commercially sensitive and a realistic parameter range is difficult to obtain.


title: Speech2Slot: An End-to-End Knowledge-based Slot Filling from Speech

abstract: In contrast to conventional pipeline Spoken Language Understanding (SLU) which consists of automatic speech recognition (ASR) and natural language understanding (NLU), end-to-end SLU infers the semantic meaning directly from speech and overcomes the error propagation caused by ASR. End-to-end slot filling (SF) from speech is an essential component of end-to-end SLU, and is usually regarded as a sequence-to-sequence generation problem, heavily relied on the performance of language model of ASR. However, it is hard to generate a correct slot when the slot is out-of-vovabulary (OOV) in training data, especially when a slot is an anti-linguistic entity without grammatical rule. Inspired by object detection in computer vision that is to detect the object from an image, we consider SF as the task of slot detection from speech. In this paper, we formulate the SF task as a matching task and propose an end-to-end knowledge-based SF model, named Speech-to-Slot (Speech2Slot), to leverage knowledge to detect the boundary of a slot from the speech. We also release a large-scale dataset of Chinese speech for slot filling, containing more than 830,000 samples. The experiments show that our approach is markedly superior to the conventional pipeline SLU approach, and outperforms the state-of-the-art end-to-end SF approach with 12.51% accuracy improvement.


title: Distributed In-memory Data Management for Workflow Executions

abstract: Complex scientific experiments from various domains are typically modeled as workflows and executed on large-scale machines using a Parallel Workflow Management System (WMS). Since such executions usually last for hours or days, some WMSs provide user steering support, i.e., they allow users to run data analyses and, depending on the results, adapt the workflows at runtime. A challenge in the parallel execution control design is to manage workflow data for efficient executions while enabling user steering support. Data access for high scalability is typically transaction-oriented, while for data analysis, it is online analytical-oriented so that managing such hybrid workloads makes the challenge even harder. In this work, we present SchalaDB, an architecture with a set of design principles and techniques based on distributed in-memory data management for efficient workflow execution control and user steering. We propose a distributed data design for scalable workflow task scheduling and high availability driven by a parallel and distributed in-memory DBMS. To evaluate our proposal, we develop d-Chiron, a WMS designed according to SchalaDB's principles. We carry out an extensive experimental evaluation on an HPC cluster with up to 960 computing cores. Among other analyses, we show that even when running data analyses for user steering, SchalaDB's overhead is negligible for workloads composed of hundreds of concurrent tasks on shared data. Our results encourage workflow engine developers to follow a parallel and distributed data-oriented approach not only for scheduling and monitoring but also for user steering.


title: Zero-power calibration of photonic circuits at cryogenic temperatures

abstract: The continual success of superconducting photon-detection technologies in quantum photonics asserts cryogenic-compatible systems as a cornerstone of full quantum photonic integration. Here, we present a way to reversibly fine-tune the optical properties of individual waveguide structures through local changes to their geometry using solidified xenon. Essentially, we remove the need for additional on-chip calibration elements, effectively zeroing the power consumption tied to reconfigurable elements, with virtually no detriment to photonic device performance. We enable passive circuit tuning in pressure-controlled environments, locally manipulating the cladding thickness over portions of optical waveguides. We realize this in a cryogenic environment, through controlled deposition of xenon gas and precise tuning of its thickness using sublimation, triggered by on-chip resistive heaters. $\pi$ phase shifts occur over a calculated length of just $L_{\pi}$ = 12.3$\pm$0.3 $\mu m$. This work paves the way towards the integration of compact, reconfigurable photonic circuits alongside superconducting detectors, devices, or otherwise.


title: On the Electrostatic Interaction between Point Charges due to Dielectrical Shielding

abstract: How will the electrostatic interaction between two point charges change if they are shielded from the other by a dielectrical slab? While the physical setting of this electromagnetic problem is relatively simple, it is easy to be wronged and the correct solution is surprisingly complicated. Here we will show a general answer using the method of images, in which the electrical field are not found by solving the Poisson's equation but by superposing an infinite number of image charges to recurrently satisfy all interfaces' boundary conditions. We also obtain analytical and algebraic results in some special cases.


title: Nanowire networks: how does small-world character evolve with dimensionality?

abstract: Networks of nanowires are currently under consideration for a wide range of electronic and optoelectronic applications. Nanowire devices are usually made by sequential deposition, which inevitably leads to stacking of the wires on top of one another. Here we demonstrate the effect of stacking on the topology of the resulting networks. We compare perfectly 2D networks with quasi3D networks, and compare both nanowire networks to the corresponding Watts Strogatz networks, which are standard benchmark systems. By investigating quantities such as clustering, path length, modularity, and small world propensity we show that the connectivity of the quasi-3D networks is significantly different to that of the 2D networks, a result which may have important implications for applications of nanowire networks.


title: Vibronic coupling in energy transfer dynamics and two-dimensional electronic-vibrational spectra

abstract: We introduce a heterodimer model in which multiple mechanisms of vibronic coupling and their impact on energy transfer can be explicitly studied. We consider vibronic coupling that arises through either Franck-Condon activity in which each site in the heterodimer has a local electron-phonon coupling and as Herzberg-Teller activity in which the transition dipole moment coupling the sites has an explicit vibrational mode-dependence. We have computed two-dimensional electronic-vibrational (2DEV) spectra for this model while varying the magnitude of these two effects and find that 2DEV spectra contain static and dynamic signatures of both types of vibronic coupling. Franck-Condon activity emerges through a change in the observed excitonic structure while Herzberg-Teller activity is evident in the appearance of significant side-band transitions that mimic the lower-energy excitonic structure. A comparison of quantum beating patterns obtained from analysis of the simulated 2DEV spectra shows that this technique can report on the mechanism of energy transfer, elucidating a means of experimentally determining the role of specific vibronic coupling mechanisms in such processes.


title: Note: Vibronic coupling in light-harvesting complex II revisited

abstract: We discuss our recent theoretical work on vibronic coupling mechanisms in a model energy transfer system in the context of previous 2DEV experiments on a natural light-harvesting system, light-harvesting complex II (LHCII), where vibronic signatures were suggested to be involved in energy transfer. In this comparison, we directly assign the vibronic coupling mechanism in LHCII as arising from Herzberg-Teller activity and show how this coupling modulates the energy transfer dynamics in this photosynthetic system.


title: On Conformal Transformation with Multiple Scalar Fields and Geometric Property of Field Space with Einstein-like Solutions

abstract: Multiple scalar fields appear in vast modern particle physics and gravity models. When they couple to gravity non-minimally, conformal transformation is utilized to bring the theory into Einstein frame. However, the kinetic terms of scalar fields are usually not canonical, which makes analytic treatment difficult. Here we investigate under what conditions the theories can be transformed to the quasi-canonical form, in which case the effective metric tensor in field space is conformally flat. We solve the relevant nonlinear partial differential equations for arbitrary number of scalar fields and present several solutions that may be useful for future phenomenological model building, including the $\sigma$-model with a particular non-minimal coupling. We also find conformal flatness can always be achieved in some modified gravity theories, for example, Starobinsky model.


title: Separate but Together: Unsupervised Federated Learning for Speech Enhancement from Non-IID Data

abstract: We propose FEDENHANCE, an unsupervised federated learning (FL) approach for speech enhancement and separation with non-IID distributed data across multiple clients. We simulate a real-world scenario where each client only has access to a few noisy recordings from a limited and disjoint number of speakers (hence non-IID). Each client trains their model in isolation using mixture invariant training while periodically providing updates to a central server. Our experiments show that our approach achieves competitive enhancement performance compared to IID training on a single device and that we can further facilitate the convergence speed and the overall performance using transfer learning on the server-side. Moreover, we show that we can effectively combine updates from clients trained locally with supervised and unsupervised losses. We also release a new dataset LibriFSD50K and its creation recipe in order to facilitate FL research for source separation problems.


title: Optimal Online Algorithms for Peak-Demand Reduction Maximization with Energy Storage

abstract: The high proportions of demand charges in electric bills motivate large-power customers to leverage energy storage for reducing the peak procurement from the outer grid. Given limited energy storage, we expect to maximize the peak-demand reduction in an online fashion, challenged by the highly uncertain demands and renewable injections, the non-cumulative nature of peak consumption, and the coupling of online decisions. In this paper, we propose an optimal online algorithm that achieves the best competitive ratio, following the idea of maintaining a constant ratio between the online and the optimal offline peak-reduction performance. We further show that the optimal competitive ratio can be computed by solving a linear number of linear-fractional programs. Moreover, we extend the algorithm to adaptively maintain the best competitive ratio given the revealed inputs and actions at each decision-making round. The adaptive algorithm retains the optimal worst-case guarantee and attains improved average-case performance. We evaluate our proposed algorithms using real-world traces and show that they obtain up to 81% peak reduction of the optimal offline benchmark. Additionally, the adaptive algorithm achieves at least 20% more peak reduction against baseline alternatives.


title: Unsupervised domain adaptation via double classifiers based on high confidence pseudo label

abstract: Unsupervised domain adaptation (UDA) aims to solve the problem of knowledge transfer from labeled source domain to unlabeled target domain. Recently, many domain adaptation (DA) methods use centroid to align the local distribution of different domains, that is, to align different classes. This improves the effect of domain adaptation, but domain differences exist not only between classes, but also between samples. This work rethinks what is the alignment between different domains, and studies how to achieve the real alignment between different domains. Previous DA methods only considered one distribution feature of aligned samples, such as full distribution or local distribution. In addition to aligning the global distribution, the real domain adaptation should also align the meso distribution and the micro distribution. Therefore, this study propose a double classifier method based on high confidence label (DCP). By aligning the centroid and the distribution between centroid and sample of different classifiers, the meso and micro distribution alignment of different domains is realized. In addition, in order to reduce the chain error caused by error marking, This study propose a high confidence marking method to reduce the marking error. To verify its versatility, this study evaluates DCP on digital recognition and target recognition data sets. The results show that our method achieves state-of-the-art results on most of the current domain adaptation benchmark datasets.


title: Travelling wave analysis of cellular invasion into surrounding tissues

abstract: Single-species reaction-diffusion equations, such as the Fisher-KPP and Porous-Fisher equations, support travelling wave solutions that are often interpreted as simple mathematical models of biological invasion. Such travelling wave solutions are thought to play a role in various applications including development, wound healing and malignant invasion. One criticism of these single-species equations is that they do not explicitly describe interactions between the invading population and the surrounding environment. In this work we study a reaction-diffusion equation that describes malignant invasion which has been used to interpret experimental measurements describing the invasion of malignant melanoma cells into surrounding human skin tissues. This model explicitly describes how the population of cancer cells degrade the surrounding tissues, thereby creating free space into which the cancer cells migrate and proliferate to form an invasion wave of malignant tissue that is coupled to a retreating wave of skin tissue. We analyse travelling wave solutions of this model using a combination of numerical simulation, phase plane analysis and perturbation techniques. Our analysis shows that the travelling wave solutions involve a range of very interesting properties that resemble certain well-established features of both the Fisher-KPP and Porous-Fisher equations, as well as a range of novel properties that can be thought of as extensions of these well-studied single-species equations. Of particular interest is that travelling wave solutions of the invasion model are very well approximated by trajectories in the Fisher-KPP phase plane that are normally disregarded. This observation establishes a previously unnoticed link between coupled multi-species reaction diffusion models of invasion and a different class of models of invasion that involve moving boundary problems.


title: An improved tile-based scalable distributed management model of massive high-resolution satellite images

abstract: The amount of remote sensing (RS) data has increased at an unexpected scale, due to the rapid progress of earth-observation and the growth of satellite RS and sensor technologies. Traditional relational databases attend their limit to meet the needs of high-resolution and large-scale RS Big Data management. As a result, massive RS data management is currently one of the most imperative topics. To address this problem, this paper describes a distributed architecture for big RS data storage based on a unified metadata file, pyramid model, and Hilbert curve for data composition and indexing using NoSQL databases (i.e, Apache Hbase). In this paper, a Hadoop-based framework in AzureInsight cloud platform is designed to manage massive RS data in a parallel and distributed way. Experimental results prove that our method has the potential to overcome the weakness of traditional methods. The proposed model is suitable for massive high-resolution image data management.


title: Recursive sequences attached to modular representations of finite groups

abstract: The core of a finite-dimensional modular representation $M$ of a finite group $G$ is its largest non-projective summand. We prove that the dimensions of the cores of $M^{\otimes n}$ have algebraic Hilbert series when $M$ is Omega-algebraic, in the sense that the non-projective summands of $M^{\otimes n}$ fall into finitely many orbits under the action of the syzygy operator $\Omega$. Similarly, we prove that these dimension sequences are eventually linearly recursive when $M$ is what we term $\Omega^{+}$-algebraic. This partially answers a conjecture by Benson and Symonds. Along the way, we also prove a number of auxiliary permanence results for linear recurrence under operations on multi-variable sequences.


title: High-dimensional coherent one-way quantum key distribution

abstract: High-dimensional quantum key distribution (QKD) offers secure communication, with secure key rates that surpass those achievable by QKD protocols utilizing two-dimensional encoding. However, existing high-dimensional QKD protocols require additional experimental resources, such as multiport interferometers and multiple detectors, thus raising the cost of practical high-dimensional systems and limiting their use. Here, we present and analyze a novel protocol for arbitrary-dimensional QKD, that requires only the hardware of a standard two-dimensional system. We provide security proofs against individual attacks and coherent attacks, setting an upper and lower bound on the secure key rates. Then, we test the new high-dimensional protocol in a standard two-dimensional QKD system over a 40 km fiber link. The new protocol yields a two-fold enhancement of the secure key rate compared to the standard two-dimensional coherent one-way protocol, without introducing any hardware modifications to the system. This work, therefore, holds great potential to enhance the performance of already deployed time-bin QKD systems through a software update alone. Furthermore, its applications extend across different encoding schemes of QKD qudits.

categories: quant-ph cs.CR cs.IT math.IT physics.optics


title: Proof of a conjecture of Dahmen and Beukers on counting integral Lam\'{e} equations with finite monodromy

abstract: In this paper, we prove Dahmen and Beukers' conjecture that the number of integral Lam\'{e} equations with index $n$ modulo scalar equivalence with the monodromy group dihedral $D_{N}$ of order $2N$ is given by \[L_{n}(N)=\frac{1}{2}\left( \frac{n(n+1)\Psi(N)}{24}-\left( a_{n}% \phi(N)+b_{n}\phi(\tfrac{N}{2}) \right) \right) +\frac{2}% {3}\varepsilon_{n}(N).\] Our main tool is the new pre-modular form $Z_{r,s}^{(n)}(\tau)$ of weight $n(n+1)/2$ introduced by Lin and Wang \cite{LW2} and the associated modular form $M_{n,N}(\tau):=\prod_{(r,s)}Z_{r,s}^{(n)}(\tau)$ of weight $\Psi(N)n(n+1)/{2}$, where the product runs over all $N$-torsion points $(r,s)$ of exact order $N$. We show that this conjecture is equivalent to the precise formula of the vanishing order of $M_{n,N}(\tau)$ at infinity: \[v_{\infty}(M_{n,N}(\tau))=a_{n}\phi(N)+b_{n}\phi( N/2).\] This formula is extremely hard to prove because the explicit expression of $Z_{r,s}^{(n)}(\tau)$ is not known for general $n$. Here we succeed to prove it by using certain Painlev\'{e} VI equations. Our result also indicates that this conjecture is intimately connected with the problem of counting pole numbers of algebraic solutions of certain Painlev\'{e} VI equations. The main results of this paper has been announced in \cite{Lin-CDM}.


title: A 3-approximation list scheduling algorithm for a single-machine scheduling problem with a non-renewable resource and total weighted completion time criterion

abstract: This paper studies a single-machine scheduling problem with a non-renewable resource (NR-SSP) and total weighted completion time criterion. The non-renewable resource is consumed when the machine starts processing a job. We consider the case where each job's weight in the objective function is proportional to its resource consumption amount. The problem is known to be NP-hard in this case. We propose a 3-approximation list scheduling algorithm for this problem. Besides, we show that the approximation ratio 3 is tight for the algorithm.


title: Quantum Embedding Theories to Simulate Condensed Systems on Quantum Computers

abstract: Quantum computers hold promise to improve the efficiency of quantum simulations of materials and to enable the investigation of systems and properties more complex than tractable at present on classical architectures. Here, we discuss computational frameworks to carry out electronic structure calculations of solids on noisy intermediate scale quantum computers using embedding theories, and we give examples for a specific class of materials, i.e., spin defects in solids. These are promising systems to build future quantum technologies, e.g., computers, sensors and devices for quantum communications. Although quantum simulations on quantum architectures are in their infancy, promising results for realistic systems appear to be within reach.


title: Wavelength conversion for single-photon polarization qubits through continuous variable quantum teleportation

abstract: A quantum internet connects remote quantum processors that need interact and exchange quantum signals over a long distance through photonic channels. However, these quantum nodes are usually composed of quantum systems with emitted photons unsuitable for long-distance transmission. Therefore, quantum wavelength conversion to telecom is crucial for long-distance quantum networks based on optical fiber. Here we propose wavelength conversion devices for single-photon polarization qubits using continuous variable quantum teleportation, which can efficiently convert qubits between near-infrared (780/795 nm suitable for interacting with atomic quantum nodes) and telecom wavelength (1300-1500 nm suitable for long-distance transmission). The teleportation uses entangled photon sources (i.e., non-degenerate two-mode squeezed state) that can be generated by four-wave mixing in rubidium atomic vapor cells, with a diamond configuration of atomic transitions. The entangled fields can be emitted in two orthogonal polarizations with locked relative phase, making them especially suitable for interfacing with single-photon polarization qubits. Our work paves the way for the realization of long-distance quantum networks.


title: Lightweight Distributed Gaussian Process Regression for Online Machine Learning

abstract: In this paper, we study the problem where a group of agents aim to collaboratively learn a common static latent function through streaming data. We propose a lightweight distributed Gaussian process regression (GPR) algorithm that is cognizant of agents' limited capabilities in communication, computation and memory. Each agent independently runs agent-based GPR using local streaming data to predict test points of interest; then the agents collaboratively execute distributed GPR to obtain global predictions over a common sparse set of test points; finally, each agent fuses results from distributed GPR with agent-based GPR to refine its predictions. By quantifying the transient and steady-state performances in predictive variance and error, we show that limited inter-agent communication improves learning performances in the sense of Pareto. Monte Carlo simulation is conducted to evaluate the developed algorithm.