title: Networks of Necessity: Simulating COVID-19 Mitigation Strategies for Disabled People and Their Caregivers

abstract: A major strategy to prevent the spread of COVID-19 is the limiting of in-person contacts. However, this is impractical or impossible for the many disabled people who do not live in care facilities, but still require caregivers. We seek to determine which interventions can prevent infections among disabled people and their caregivers. We simulate transmission with a model that includes susceptible, exposed, asymptomatic, symptomatically ill, hospitalized, and removed individuals. The networks on which we simulate disease spread incorporate heterogeneity in the risks of different types of interactions, time-dependent lockdown and reopening measures, and contact distributions for four different groups (caregivers, disabled people, essential workers, and the general population). We find the probability of becoming infected is largest for caregivers and second largest for disabled people. Our analysis of network structure illustrates that caregivers have the largest modal eigenvector centrality. We find that two interventions -- contact-limiting by all groups and mask-wearing by disabled people and caregivers -- most reduce the cases among disabled people and caregivers. We also test which group spreads COVID-19 most readily by seeding infections in a subset of each group. We find caregivers are the most potent spreaders of COVID-19, particularly to other caregivers and to disabled people. We test where to use limited vaccine doses most effectively and find (1) vaccinating caregivers better protects disabled people than vaccinating the general population or essential workers and (2) vaccinating caregivers protects disabled people about as much as vaccinating disabled people themselves. Our results highlight the potential effectiveness of mask-wearing, contact-limiting throughout society, and strategic vaccination for limiting the exposure of disabled people and their caregivers to COVID-19.

categories: cs.SI math.DS nlin.AO physics.soc-ph q-bio.PE


title: The Query Complexity of Local Search and Brouwer in Rounds

abstract: We consider the query complexity of finding a local minimum of a function defined on a graph, where at most $k$ rounds of interaction with the oracle are allowed. Rounds model parallel settings, where each query takes resources to complete and is executed on a separate processor. Thus the query complexity in $k$ rounds informs how many processors are needed to achieve a parallel time of $k$. We focus on the d-dimensional grid $[n]^d$, where the dimension $d$ is a constant, and consider two regimes for the number of rounds: constant and polynomial in n. We give algorithms and lower bounds that characterize the trade-off between the number of rounds of adaptivity and the query complexity of local search. When the number of rounds $k$ is constant, we show that the query complexity of local search in $k$ rounds is $\Theta\bigl(n^{\frac{d^{k+1} - d^k}{d^k - 1}}\bigl)$, for both deterministic and randomized algorithms. When the number of rounds is polynomial, i.e. $k = n^{\alpha}$ for $0 < \alpha < d/2$, the randomized query complexity is $\Theta\left(n^{d-1 - \frac{d-2}{d}\alpha}\right)$ for all $d \geq 5$. For $d=3$ and $d=4$, we show the same upper bound expression holds and give almost matching lower bounds. The local search analysis also enables us to characterize the query complexity of computing a Brouwer fixed point in rounds. Our proof technique for lower bounding the query complexity in rounds may be of independent interest as an alternative to the classical relational adversary method of Aaronson from the fully adaptive setting.


title: FGF-GAN: A Lightweight Generative Adversarial Network for Pansharpening via Fast Guided Filter

abstract: Pansharpening is a widely used image enhancement technique for remote sensing. Its principle is to fuse the input high-resolution single-channel panchromatic (PAN) image and low-resolution multi-spectral image and to obtain a high-resolution multi-spectral (HRMS) image. The existing deep learning pansharpening method has two shortcomings. First, features of two input images need to be concatenated along the channel dimension to reconstruct the HRMS image, which makes the importance of PAN images not prominent, and also leads to high computational cost. Second, the implicit information of features is difficult to extract through the manually designed loss function. To this end, we propose a generative adversarial network via the fast guided filter (FGF) for pansharpening. In generator, traditional channel concatenation is replaced by FGF to better retain the spatial information while reducing the number of parameters. Meanwhile, the fusion objects can be highlighted by the spatial attention module. In addition, the latent information of features can be preserved effectively through adversarial training. Numerous experiments illustrate that our network generates high-quality HRMS images that can surpass existing methods, and with fewer parameters.


title: EarlyBERT: Efficient BERT Training via Early-bird Lottery Tickets

abstract: Heavily overparameterized language models such as BERT, XLNet and T5 have achieved impressive success in many NLP tasks. However, their high model complexity requires enormous computation resources and extremely long training time for both pre-training and fine-tuning. Many works have studied model compression on large NLP models, but only focusing on reducing inference time while still requiring an expensive training process. Other works use extremely large batch sizes to shorten the pre-training time, at the expense of higher computational resource demands. In this paper, inspired by the Early-Bird Lottery Tickets recently studied for computer vision tasks, we propose EarlyBERT, a general computationally-efficient training algorithm applicable to both pre-training and fine-tuning of large-scale language models. By slimming the self-attention and fully-connected sub-layers inside a transformer, we are the first to identify structured winning tickets in the early stage of BERT training. We apply those tickets towards efficient BERT training, and conduct comprehensive pre-training and fine-tuning experiments on GLUE and SQuAD downstream tasks. Our results show that EarlyBERT achieves comparable performance to standard BERT, with 35~45% less training time. Code is available at https://github.com/VITA-Group/EarlyBERT.


title: Adaptive Quantile Computation for Brownian Bridge in Change-Point Analysis

abstract: As an example for the fast calculation of distributional parameters of Gaussian processes, we propose a new Monte Carlo algorithm for the computation of quantiles of the supremum norm of weighted Brownian bridges. As it is known, the corresponding distributions arise asymptotically for weighted CUSUM statistics for change-point detection. The new algorithm employs an adaptive (sequential) time discretization for the trajectories of the Brownian bridge. A simulation study shows that the new algorithm by far outperforms the standard approach, which employs a uniform time discretization.


title: The growing topology of the C. elegans connectome

abstract: Probing the developing neural circuitry in Caenorhabditis elegans has enhanced our understanding of nervous systems. The C. elegans connectome, like those of other species, is characterized by a rich club of densely connected neurons embedded within a small-world architecture. This organization of neuronal connections, captured by quantitative network statistics, provides insight into the system's capacity to perform integrative computations. Yet these network measures are limited in their ability to detect weakly connected motifs, such as topological cavities, that may support the systems capacity to perform segregated computations. We address this limitation by using persistent homology to track the evolution of topological cavities in the growing C. elegans connectome throughout neural development, and assess the degree to which the growing connectomes topology is resistant to biological noise. We show that the developing connectome topology is both relatively robust to changes in neuron birth times and not captured by similar growth models. Additionally, we quantify the consequence of a neurons specific birth time and ask if this metric tracks other biological properties of neurons. Our results suggest that the connectomes growing topology is a robust feature of the developing connectome that is distinct from other network properties, and that the growing topology is particularly sensitive to the exact birth times of a small set of predominantly motor neurons. By utilizing novel measurements that track biological features, we anticipate that our study will be helpful in the construction of more accurate models of neuronal development in C. elegans


title: Radio frequency mixing modules for superconducting qubit room temperature control systems

abstract: As the number of qubits in nascent quantum processing units increases, the connectorized RF (radio frequency) analog circuits used in first generation experiments become exceedingly complex. The physical size, cost and electrical failure rate all become limiting factors in the extensibility of control systems. We have developed a series of compact RF mixing boards to address this challenge by integrating I/Q quadrature mixing, IF(intermediate frequency)/LO(local oscillator)/RF power level adjustments, and DC (direct current) bias fine tuning on a 40 mm $\times $ 80 mm 4-layer PCB (printed circuit board) board with EMI (electromagnetic interference) shielding. The RF mixing module is designed to work with RF and LO frequencies between 2.5 and 8.5 GHz. The typical image rejection and adjacent channel isolation are measured to be $\sim$27 dBc and $\sim$50 dB. By scanning the drive phase in a loopback test, the module short-term amplitude and phase linearity are typically measured to be 5$\times$10$^{-4}$ (V$_{\mathrm{pp}}$/V$_{\mathrm{mean}}$) and 1$\times$10$^{-3}$ radian (pk-pk). The operation of RF mixing board was validated by integrating it into the room temperature control system of a superconducting quantum processor and executing randomized benchmarking characterization of single and two qubit gates. We measured a single-qubit process infidelity of $9.3(3) \times 10^{-4}$ and a two-qubit process infidelity of $2.7(1) \times 10^{-2}$.


title: Screening and understanding Li adsorption on 2-dimensional metallic materials by learning physics

abstract: Two-dimensional (2D) materials have received considerable attention as possible electrodes in Li-ion batteries (LIBs), although a deeper understanding of the Li adsorption behavior as well as broad screening of the materials space is still needed. In this work, we build a high-throughput screening scheme that incorporates a learned interaction. First, density functional theory and graph convolution networks are utilized to calculate minimum Li adsorption energies for a small set of 2D metallic materials. The data is then used to find a dependence of the minimum Li adsorption energies on the sum of ionization potential, work function of the 2D metal, and coupling energy between Li+ and substrate. Our results show that variances of elemental properties and density are the most correlated features with coupling. To illustrate the applicability of this approach, the model is employed to show that some fluorides and chromium oxides are potential high-voltage materials with adsorption energies < -7 eV, and the found physics is used as the design principle to enhance the Li adsorption ability of graphene. This physics-driven approach shows higher accuracy and transferability compared with purely data-driven models.


title: Toward Reliable Designs of Data-Driven Reinforcement Learning Tracking Control for Euler-Lagrange Systems

abstract: This paper addresses reinforcement learning based, direct signal tracking control with an objective of developing mathematically suitable and practically useful design approaches. Specifically, we aim to provide reliable and easy to implement designs in order to reach reproducible neural network-based solutions. Our proposed new design takes advantage of two control design frameworks: a reinforcement learning based, data-driven approach to provide the needed adaptation and (sub)optimality, and a backstepping based approach to provide closed-loop system stability framework. We develop this work based on an established direct heuristic dynamic programming (dHDP) learning paradigm to perform online learning and adaptation and a backstepping design for a class of important nonlinear dynamics described as Euler-Lagrange systems. We provide a theoretical guarantee for the stability of the overall dynamic system, weight convergence of the approximating nonlinear neural networks, and the Bellman (sub)optimality of the resulted control policy. We use simulations to demonstrate significantly improved design performance of the proposed approach over the original dHDP.


title: Equivariant K-theory and Resolution II: Non-Abelian actions

abstract: The smooth action of a compact Lie group on a compact manifold can be resolved to an iterated space, as made explicit by Pierre Albin and the second author. On the resolution the lifted action has fixed isotropy type corresponding to the open stratum and also in an iterated sense, with connecting equivariant fibrations over the boundary hypersurfaces covering the resolutions of the other strata. This structure descends to a resolution of the quotient as a stratified space. For an Abelian group action the equivariant K-theory can then be described in terms of bundles over the bases `dressed' by the representations of the isotropy types with morphisms covering the connecting maps. A similar model is given here covering the non-Abelian case. Now the reduced objects are torsion-twisted bundles over finite covers of the bases, corresponding to the projective action of the normalizers on the representations of the isotropy groups, again with morphisms over all the boundaries. This leads to a closely related iterated deRham model for equivariant cohomology and, now with values in forms twisted by flat bundles of representation rings over the bases, for delocalized equivariant cohomology. We show, as envisioned by Baum, Brylinksi and MacPherson, that the usual equivariant Chern character, mapping to equivariant cohomology, factors through a natural Chern character from equivariant K-theory to delocalized equivariant cohomology with the latter giving an Atiyah-Hirzebruch isomorphism.


title: Homological bulk-edge correspondence for Weyl semimetal

abstract: For a certain translation invariant tight-binding model of three-dimensional Weyl semimetals, we establish a bulk-edge correspondence as an equality of two relative homology classes, based on an idea of Mathai and Thiang: From spectral information on the edge Hamiltonian, we construct a relative homology class on the surface momentum space. This class agrees with the image under the surface projection of a homology class on the bulk momentum space relative to the Weyl points, constructed from the bulk Hamiltonian. Furthermore, the relative homology class on the surface momentum space can be represented by homology cycles whose images constitute the Fermi arc, the locus where the edge Hamiltonian admits zero spectrum.


title: QubiC: An open source FPGA-based control and measurement system for superconducting quantum information processors

abstract: As quantum information processors grow in quantum bit (qubit) count and functionality, the control and measurement system becomes a limiting factor to large scale extensibility. To tackle this challenge and keep pace with rapidly evolving classical control requirements, full control stack access is essential to system level optimization. We design a modular FPGA (field-programmable gate array) based system called QubiC to control and measure a superconducting quantum processing unit. The system includes room temperature electronics hardware, FPGA gateware, and engineering software. A prototype hardware module is assembled from several commercial off-the-shelf evaluation boards and in-house developed circuit boards. Gateware and software are designed to implement basic qubit control and measurement protocols. System functionality and performance are demonstrated by performing qubit chip characterization, gate optimization, and randomized benchmarking sequences on a superconducting quantum processor operating at the Advanced Quantum Testbed at Lawrence Berkeley National Laboratory. The single-qubit and two-qubit process fidelities are measured to be 0.9980$\pm$0.0001 and 0.948$\pm$0.004 by randomized benchmarking. With fast circuit sequence loading capability, the QubiC performs randomized compiling experiments efficiently and improves the feasibility of executing more complex algorithms.


title: Explicit regularization and implicit bias in deep network classifiers trained with the square loss

abstract: Deep ReLU networks trained with the square loss have been observed to perform well in classification tasks. We provide here a theoretical justification based on analysis of the associated gradient flow. We show that convergence to a solution with the absolute minimum norm is expected when normalization techniques such as Batch Normalization (BN) or Weight Normalization (WN) are used together with Weight Decay (WD). The main property of the minimizers that bounds their expected error is the norm: we prove that among all the close-to-interpolating solutions, the ones associated with smaller Frobenius norms of the unnormalized weight matrices have better margin and better bounds on the expected classification error. With BN but in the absence of WD, the dynamical system is singular. Implicit dynamical regularization -- that is zero-initial conditions biasing the dynamics towards high margin solutions -- is also possible in the no-BN and no-WD case. The theory yields several predictions, including the role of BN and weight decay, aspects of Papyan, Han and Donoho's Neural Collapse and the constraints induced by BN on the network weights.


title: A Multi-modal Deep Learning Model for Video Thumbnail Selection

abstract: Thumbnail is the face of online videos. The explosive growth of videos both in number and variety underpins the importance of a good thumbnail because it saves potential viewers time to choose videos and even entice them to click on them. A good thumbnail should be a frame that best represents the content of a video while at the same time capturing viewers' attention. However, the techniques and models in the past only focus on frames within a video, and we believe such narrowed focus leave out much useful information that are part of a video. In this paper, we expand the definition of content to include title, description, and audio of a video and utilize information provided by these modalities in our selection model. Specifically, our model will first sample frames uniformly in time and return the top 1,000 frames in this subset with the highest aesthetic scores by a Double-column Convolutional Neural Network, to avoid the computational burden of processing all frames in downstream task. Then, the model incorporates frame features extracted from VGG16, text features from ELECTRA, and audio features from TRILL. These models were selected because of their results on popular datasets as well as their competitive performances. After feature extraction, the time-series features, frames and audio, will be fed into Transformer encoder layers to return a vector representing their corresponding modality. Each of the four features (frames, title, description, audios) will pass through a context gating layer before concatenation. Finally, our model will generate a vector in the latent space and select the frame that is most similar to this vector in the latent space. To the best of our knowledge, we are the first to propose a multi-modal deep learning model to select video thumbnail, which beats the result from the previous State-of-The-Art models.


title: Three-quarter Sibling Regression for Denoising Observational Data

abstract: Many ecological studies and conservation policies are based on field observations of species, which can be affected by systematic variability introduced by the observation process. A recently introduced causal modeling technique called 'half-sibling regression' can detect and correct for systematic errors in measurements of multiple independent random variables. However, it will remove intrinsic variability if the variables are dependent, and therefore does not apply to many situations, including modeling of species counts that are controlled by common causes. We present a technique called 'three-quarter sibling regression' to partially overcome this limitation. It can filter the effect of systematic noise when the latent variables have observed common causes. We provide theoretical justification of this approach, demonstrate its effectiveness on synthetic data, and show that it reduces systematic detection variability due to moon brightness in moth surveys.


title: Investigation of thrust augmentation and nozzle exit pressure in starting jets from a tube nozzle

abstract: It has been widely accepted that nozzle exit over-pressure is responsible for thrust augmentation in starting jets over steady jets and a few pulsed jet propulsion systems have been developed based on this finding. However, no other study has been conducted to confirm the nozzle exit over-pressure effect. In this paper, thrust and nozzle exit pressure are numerically investigated under different jet velocity programs (negative slope (NS), positive slope (PS) and impulsive). Starting jets are generated for fluid slug length to diameter (L/D) ratios of 2-5 at each jet velocity program. In contradiction to the findings presented in [P. Krueger and M. Gharib, J. Phys. Fluids, 15, 1271 (2003)] the results in this paper reveal that the time-averaged thrust due to nozzle exit pressure in starting jets from a tube nozzle is negligible or negative resulting in no thrust augmentation when compared to steady jets with hypothetical identical velocity programs. Also, nozzle exit experiences both over-pressure and under-pressure in accelerating and decelerating portions of jet velocity programs respectively in both NS and PS velocity programs. Furthermore, it is found that secondary counter-rotating vortex is generated during decelerating portion of jet velocity programs which may be linked to the nozzle exit underpressure explaining the reduced thrust during starting jets generated by different velocity programs.


title: Independence and almost automorphy of high order

abstract: In this paper, it is shown that for a minimal system $(X,T)$ and $d,k\in \mathbb{N}$, if $(x,x_i)$ is regionally proximal of order $d$ for $1\leq i\leq k$, then $(x,x_1,\ldots,x_k)$ is $(k+1)$-regionally proximal of order $d$. Meanwhile, we introduce the notion of $\mathrm{IN}^{[d]}$-pair: for a dynamical system $(X,T)$ and $d\in \mathbb{N}$, a pair $(x_0,x_1)\in X\times X$ is called an $\mathrm{IN}^{[d]}$-pair if for any $k\in \mathbb{N}$ and any neighborhoods $U_0 ,U_1 $ of $x_0$ and $x_1$ respectively, there exist integers $p_j^{(i)},1\leq i\leq k,$ $1\leq j\leq d$ such that $$ \bigcup_{i=1}^k\{ p_1^{(i)}\epsilon(1)+\ldots+p_d^{(i)} \epsilon(d):\epsilon(j)\in \{0,1\},1\leq j\leq d\}\backslash \{0\}\subset \mathrm{Ind}(U_0,U_1), $$ where $\mathrm{Ind}(U_0,U_1)$ denotes the collection of all independence sets for $(U_0,U_1)$. It turns out that for a minimal system, if it dose not contain any nontrivial $\mathrm{IN}^{[d]}$-pair, then it is an almost one-to-one extension of its maximal factor of order $d$.


title: Presentations for Globular Operads

abstract: In this paper we develop the theory of presentations for globular operads and construct presentations for the globular operads corresponding to several key theories of $n$-category for $n \leqslant 4$.


title: Controlled Analyses of Social Biases in Wikipedia Bios

abstract: Social biases on Wikipedia, a widely-read global platform, could greatly influence public opinion. While prior research has examined man/woman gender bias in biography articles, possible influences of other demographic attributes limit conclusions. In this work, we present a methodology for analyzing Wikipedia pages about people that isolates dimensions of interest (e.g., gender), from other attributes (e.g., occupation). Given a target corpus for analysis (e.g.~biographies about women), we present a method for constructing a comparison corpus that matches the target corpus in as many attributes as possible, except the target one. We develop evaluation metrics to measure how well the comparison corpus aligns with the target corpus and then examine how articles about gender and racial minorities (cis. women, non-binary people, transgender women, and transgender men; African American, Asian American, and Hispanic/Latinx American people) differ from other articles. In addition to identifying suspect social biases, our results show that failing to control for covariates can result in different conclusions and veil biases. Our contributions include methodology that facilitates further analyses of bias in Wikipedia articles, findings that can aid Wikipedia editors in reducing biases, and a framework and evaluation metrics to guide future work in this area.


title: Graph Networks with Spectral Message Passing

abstract: Graph Neural Networks (GNNs) are the subject of intense focus by the machine learning community for problems involving relational reasoning. GNNs can be broadly divided into spatial and spectral approaches. Spatial approaches use a form of learned message-passing, in which interactions among vertices are computed locally, and information propagates over longer distances on the graph with greater numbers of message-passing steps. Spectral approaches use eigendecompositions of the graph Laplacian to produce a generalization of spatial convolutions to graph structured data which access information over short and long time scales simultaneously. Here we introduce the Spectral Graph Network, which applies message passing to both the spatial and spectral domains. Our model projects vertices of the spatial graph onto the Laplacian eigenvectors, which are each represented as vertices in a fully connected "spectral graph", and then applies learned message passing to them. We apply this model to various benchmark tasks including a graph-based variant of MNIST classification, molecular property prediction on MoleculeNet and QM9, and shortest path problems on random graphs. Our results show that the Spectral GN promotes efficient training, reaching high performance with fewer training iterations despite having more parameters. The model also provides robustness to edge dropout and outperforms baselines for the classification tasks. We also explore how these performance benefits depend on properties of the dataset.


title: Elastocaloric signature of nematic fluctuations

abstract: The elastocaloric effect (ECE) is a thermodynamic quantity relating changes in entropy to changes in strain experienced by a material. As such, ECE measurements can provide valuable information about the entropy landscape proximate to strain-tuned phase transitions. For ordered states that break only point symmetries, bilinear coupling of the order parameter with strain implies that the ECE can also provide a window on fluctuations above the critical temperature, and hence, in principle, can also provide a thermodynamic measure of the associated susceptibility. To demonstrate this, we use the ECE to sensitively reveal the presence of nematic fluctuations in the archetypal Fe-based superconductor Ba(Fe$_{1-x}$Co$_{x}$)$_2$As$_2$. By performing these measurements simultaneously with elastoresistivity in a multimodal fashion, we are able to make a direct and unambiguous comparison of these closely related thermodynamic and transport properties, both of which are sensitive to nematic fluctuations. As a result, we have uncovered an unanticipated doping-dependence of the nemato-elastic coupling and of the magnitude of the scattering of low energy quasi-particles by nematic fluctuations -- while the former weakens, the latter increases dramatically with increasing doping.


title: Detection in Molecular Communications with Ligand Receptors under Molecular Interference

abstract: Molecular Communications (MC) is a bio-inspired communication technique that uses molecules to transfer information among bio-nano devices. In this paper, we focus on the detection problem for biological MC receivers employing ligand receptors to infer the transmitted messages encoded into the concentration of molecules, i.e., ligands. In practice, receptors are not ideally selective against target ligands, and in physiological environments, they can interact with multiple types of ligands at different reaction rates depending on their binding affinity. This molecular cross-talk can cause a substantial interference on MC. Here we consider a particular scenario, where there is non-negligible concentration of interferer molecules in the channel, which have similar receptor-binding characteristics with the information molecules, and the receiver employs single type of receptors. We investigate the performance of four different detection methods, which make use of different statistics of the ligand-receptor binding reactions: instantaneous number of bound receptors, unbound time durations of receptors, bound time durations of receptors, and combination of unbound and bound time durations of receptors within a sampling time interval. The performance of the introduced detection methods are evaluated in terms of bit error probability for varying strength of molecular interference, similarity between information and interferer molecules, number of receptors, and received concentration difference between bit-0 and bit-1 transmissions. We propose synthetic receptor designs that can convert the required receptor statistics to the concentration of intracellular molecules, and chemical reaction networks that can chemically perform the computations required for detection.


title: Bosonic Random Walk Networks for Graph Learning

abstract: The development of Graph Neural Networks (GNNs) has led to great progress in machine learning on graph-structured data. These networks operate via diffusing information across the graph nodes while capturing the structure of the graph. Recently there has also seen tremendous progress in quantum computing techniques. In this work, we explore applications of multi-particle quantum walks on diffusing information across graphs. Our model is based on learning the operators that govern the dynamics of quantum random walkers on graphs. We demonstrate the effectiveness of our method on classification and regression tasks.


title: Gauge Principle and Gauge Invariance in Two-Level Systems

abstract: The quantum Rabi model is a widespread description of the coupling between a two-level system and a quantized single mode of an electromagnetic resonator. Issues about this model's gauge invariance have been raised. These issues become evident when the light-matter interaction reaches the so-called ultrastrong coupling regime. Recently, a modified quantum Rabi model able to provide gauge-invariant physical results (e.g., energy levels, expectation values of observables, quantum probabilities) in any interaction regime was introduced [Nature Physics 15, 803 (2019)]. Here we provide an alternative derivation of this result, based on the implementation in two-state systems of the gauge principle, which is the principle from which all the fundamental interactions in quantum field theory are derived. The adopted procedure can be regarded as the two-site version of the general method used to implement the gauge principle in lattice gauge theories. Applying this method, we also obtain the gauge-invariant quantum Rabi model for asymmetric two-state systems, and the multi-mode gauge-invariant quantum Rabi model beyond the dipole approximation.


title: Towards Threshold Key Exchange Protocols

abstract: Threshold schemes exist for many cryptographic primitives like signatures, key derivation functions, and ciphers. At the same time, practical key exchange protocols based on Diffie-Hellman (DH) or ECDSA primitives are not designed or implemented in a threshold setting. In this paper, we implement popular key exchange protocols in a threshold manner and show that this approach can be used in practice. First, we introduce two basic threshold DH key agreement schemes that provide enhanced security features in comparison with the classic DH primitive: dealerless distributed key generation, threshold shared key computation, and private key shares refreshing. We implemented the proposed DH schemes within WireGuard protocol to demonstrate its effectiveness, efficiency, and usability in practice. The open question is the security of the proposed schemes and their instantiation from the elliptic curves used in key agreement protocols: NIST curves, Russian GOST curves, and Curve25519. Second, we propose an idea of implementing TLS in a threshold setting that can be used instead of Keyless SSL/TLS technology, and provide the measurements of TLS key exchanges based on threshold ECDSA. Even if we don't provide any formal definitions, security analysis, and mathematical proofs, we believe that the ideas and mechanisms suggested in this paper can be interesting and useful. The main intention of the paper is to start discussions and raise awareness of the challenges and problems arising when moving to threshold key exchange protocols.


title: Moderate deviations for systems of slow-fast stochastic reaction-diffusion equations

abstract: The goal of this paper is to study the Moderate Deviation Principle (MDP) for a system of stochastic reaction-diffusion equations with a time-scale separation in slow and fast components and small noise in the slow component. Based on weak convergence methods in infinite dimensions and related stochastic control arguments, we obtain an exact form for the moderate deviations rate function in different regimes as the small noise and time-scale separation parameters vanish. Many issues that appear due to the infinite dimensionality of the problem are completely absent in their finite-dimensional counterpart. In comparison to corresponding Large Deviation Principles, the moderate deviation scaling necessitates a more delicate approach to establishing tightness and properly identifying the limiting behavior of the underlying controlled problem. The latter involves regularity properties of a solution of an associated elliptic Kolmogorov equation on Hilbert space along with a finite-dimensional approximation argument.


title: calculus: High Dimensional Numerical and Symbolic Calculus in R

abstract: The R package calculus implements C++ optimized functions for numerical and symbolic calculus, such as the Einstein summing convention, fast computation of the Levi-Civita symbol and generalized Kronecker delta, Taylor series expansion, multivariate Hermite polynomials, high-order derivatives, ordinary differential equations, differential operators and numerical integration in arbitrary orthogonal coordinate systems. The library applies numerical methods when working with R functions or symbolic programming when working with characters or expressions. The package handles multivariate numerical calculus in arbitrary dimensions and coordinates and implements the symbolic counterpart of the numerical methods whenever possible, without depending on external computer algebra systems. Except for Rcpp, the package has no strict dependencies in order to provide a stable self-contained toolbox that invites re-use.


title: Theory of Faradaically Modulated Redox Active Electrodes for Electrochemically Mediated Selective Adsorption Processes

abstract: Electrochemically mediated selective adsorption is an emerging electrosorption technique that utilizes Faradaically enhanced redox active electrodes, which can adsorb ions not only electrostatically, but also electrochemically. The superb selectivity (>100) of this technique enables selective removal of toxic or high-value target ions under low energy consumption. Here, we develop a general theoretical framework to describe the competitive electrosorption phenomena involving multiple ions and surface-bound redox species. The model couples diffusion, convection and electromigration with competitive surface adsorption reaction kinetics, consistently derived from non-equilibrium thermodynamics. To optimize the selective removal of the target ions, design criteria were derived analytically from physically relevant dimensionless groups and time scales, where the propagation of the target anions concentration front is the limiting step. Detailed computational studies are reported for three case studies that cover a wide range of inlet concentration ratios between the competing ions. And in all three cases, target anions in the electrosorption cell forms a self-sharpening reaction-diffusion wave front. Based on the model, a three-step stop-flow operation scheme with a pure stripping solution of target anions is proposed that optimizes the ion adsorption performance and increases the purity of the regeneration stream to almost 100%, which is beneficial for downstream processing.


title: Canonical embeddings of pairs of arcs

abstract: We show that for given four points in the Riemann sphere and a given isotopy class of two disjoint arcs connecting these points in two pairs, there exists a unique configuration with the property that each arc is a hyperbolic geodesic segment in the complement of the other arc.


abstract: Asymptotic expansion of a variation with anticipative weights is derived by the theory of asymptotic expansion for Skorohod integrals having a mixed normal limit. The expansion formula is expressed with the quasi-torsion, quasi-tangent and other random symbols. To specify these random symbols, it is necessary to classify the level of the effect of each term appearing in the stochastic expansion of the variable in question. To solve this problem, we consider a class ${\cal L}$ of certain sequences $({\cal I}_n)_{n\in{\mathbb N}}$ of Wiener functionals and we give a systematic way of estimation of the order of $({\cal I}_n)_{n\in{\mathbb N}}$. Based on this method, we introduce a notion of exponent of the sequence $({\cal I}_n)_{n\in{\mathbb N}}$, and investigate the stability and contraction effect of the operators $D_{u_n}$ and $D$ on ${\cal L}$, where $u_n$ is the integrand of a Skorohod integral. After constructed these machineries, we derive asymptotic expansion of the variation having anticipative weights. An application to robust volatility estimation is mentioned.


title: PHP code smells in web apps: survival and anomalies

abstract: Context: Code smells are considered symptoms of poor design, leading to future problems, such as reduced maintainability. Except for anecdotal cases (e. g. code dropout), a code smell survives until it gets explicitly refactored or removed. This paper presents a longitudinal study on the survival of code smells for web apps built with PHP. Objectives: RQ: (i) code smells survival depends on their scope? (ii) practitioners attitudes towards code smells removal in web apps have changed throughout time? (iii) how long code smells survive in web applications? (iv) are there sudden variations (anomalies) in the density of code smells through the evolution of web apps? Method: We analyze the evolution of 6 code smells in 8 web applications written in PHP at the server side, across several years, using the survival analysis technique. We classify code smells according to scope in two categories: scattered and localized. Scattered code smells are expected to be more harmful since their influence is not circumscribed as in localized code smells. We split the observations for each web app into two equal and consecutive timeframes, to test the hypothesis that code smells awareness has increased throughout time. As for the anomalies, we standardize their detection criteria. Results: We present some evidence that code smells survival depends on their scope: the average survival rate decreases in some of them, while the opposite is observed for the remainder. The survival of localized code smells is around 4 years, while the scattered ones live around 5 years. Around 60% of the smells are removed, and some live through all the application life. We also show how a graphical representation of anomalies found in the evolution of code smells allows unveiling the story of a development project and make managers aware of the need for enforcing regular refactoring practices.


title: Equivalent circuit model for electrosorption with redox active materials

abstract: Electrosorption is a promising technique for brackish water deionization and waste water remediation. Faradaic materials with redox activity have recently been shown to enhance both the adsorption capacity and the selectivity of electrosorption processes. Development of the theory of electrosorption with redox active materials can provide a fundamental understanding of the electrosorption mechanism and a means to extract material properties from small-scale experiments for process optimization and scale-up. Here, we present an intuitive, physics-based equivalent circuit model to describe the electrosorption performance of redox active materials, which is able to accurately fit experimental cyclic voltammetry measurements. The model can serve as an efficient and easy-to-implement tool to evaluate properties of redox active materials and help to distinguish between the transport-limited and reaction-limited regimes in electrosorption processes. And the extracted intrinsic material properties can be further incorporated into process models under lower supporting electrolyte concentrations for realistic electrosorption applications.


title: Fuzzy Rate Analysis of Operators and its Applications in Linear Spaces

abstract: In this paper, a new concept, the fuzzy rate of an operator in linear spaces is proposed for the very first time. Some properties and basic principles of it are studied. Fuzzy rate of an operator $B$ which is specific in a plane is discussed. As its application, a new fixed point existence theorem is proved.


title: A representation type of compression space of rank 2

abstract: In this paper, we proved that a compression space of rank $2$ is equivalent to an irreducible representation over a Lie algebra.


title: Family of prime-representing constants: use of the ceiling function

abstract: The analysis of regularities and randomness in the distribution of prime numbers remains at the research frontiers for many generations of mathematicians from different groups and topical fields. In 2019 D. Fridman et al. (Am. Math. Mon. 2019, 126:1, 70-73) have suggested the constant $f_1 = 2.9200509773...$ for generation of the complete sequence of primes with using of a recursive relation for $f_n$ such that the floor function $\lfloor f_n \rfloor = p_n$, where $p_n$ is the nth prime. Here I present the family of constants $h_n (h_1 = 1.2148208055...)$ such that the ceiling function $\lceil h_n \rceil = p_n$. The proposed recursive relation $h_n=\lceil h_n \rceil(h_{n-1}-\lceil h_{n-1} \rceil+2)$ generates the sequence of all known prime numbers. I also show that constants $h_n$ are irrational.


title: Analysis of a new chaotic system, electronic realization and use in navigation of differential drive mobile robot

abstract: This paper presents a new chaotic system having four attractors, including two fixed point attractors and two symmetrical chaotic strange attractors. Dynamical properties of the system, viz. sensitive dependence on initial conditions, Lyapunov spectrum, strangeness measure, attraction basin, including the class and size of it, existence of strange attractor, bifurcation analysis, multistability, electronic circuit design, and hardware implementation, are rigorously treated. Numerical computations are used to compute the basin of attraction and show that the system has a far-reaching composite basin of attraction. Such a basin of attraction is vital for engineering applications. Moreover, a circuit model of the system is realized using analog electronic components. A procedure is detailed for converting the system parameters into corresponding electronic component values such as the circuital resistances while ensuring the dynamic ranges are bounded. Besides, the system is used as the source of control inputs for independent navigation of a differential drive mobile robot, which is subject to the Pfaffian velocity constraint. Due to the properties of sensitivity on initial conditions and topological mixing, the robot's path becomes unpredictable and guaranteed to scan the workspace, respectively.


title: Generalized Hardy Type and Caffarelli-Kohn-Nirenberg Type Inequalities on Finsler Manifolds

abstract: In this paper we derive both local and global geometric inequalities on general Riemannnian and Finsler manifolds and prove generalized Caffarelli-Kohn-Nirenberg type and Hardy type inequalities on Finsler manifolds, illuminating curvatures of both Riemannian and Finsler manifolds influence geometric inequalities.


title: Adaptive Surgical Robotic Training Using Real-Time Stylistic Behavior Feedback Through Haptic Cues

abstract: Surgical skill directly affects surgical procedure outcomes; thus, effective training is needed to ensure satisfactory results. Many objective assessment metrics have been developed and some are widely used in surgical training simulators. These objective metrics provide the trainee with descriptive feedback about their performance however, often lack feedback on how to proceed to improve performance. The most effective training method is one that is intuitive, easy to understand, personalized to the user and provided in a timely manner. We propose a framework to enable user-adaptive training using near-real-time detection of performance, based on intuitive styles of surgical movements (e.g., fluidity, smoothness, crispness, etc.), and propose a haptic feedback framework to assist with correcting styles of movement. We evaluate the ability of three types of force feedback (spring, damping, and spring plus damping feedback), computed based on prior user positions, to improve different stylistic behaviors of the user during kinematically constrained reaching movement tasks. The results indicate that four out of the six styles studied here were statistically significantly improved (p<0.05) using spring guidance force feedback and a significant reduction in task time was also found using spring feedback. The path straightness and targeting error in the task were other task performance metrics studied which were improved significantly using the spring-damping feedback. This study presents a groundwork for adaptive training in robotic surgery based on near-real-time human-centric models of surgical behavior.


title: OralViewer: 3D Demonstration of Dental Surgeries for Patient Education with Oral Cavity Reconstruction from a 2D Panoramic X-ray

abstract: Patient's understanding on forthcoming dental surgeries is required by patient-centered care and helps reduce fear and anxiety. Due to the gap of expertise between patients and dentists, conventional techniques of patient education are usually not effective for explaining surgical steps. In this paper, we present \textit{OralViewer} -- the first interactive application that enables dentist's demonstration of dental surgeries in 3D to promote patients' understanding. \textit{OralViewer} takes a single 2D panoramic dental X-ray to reconstruct patient-specific 3D teeth structures, which are then assembled with registered gum and jaw bone models for complete oral cavity modeling. During the demonstration, \textit{OralViewer} enables dentists to show surgery steps with virtual dental instruments that can animate effects on a 3D model in real-time. A technical evaluation shows our deep learning based model achieves a mean Intersection over Union (IoU) of 0.771 for 3D teeth reconstruction. A patient study with 12 participants shows \textit{OralViewer} can improve patients' understanding of surgeries. An expert study with 3 board-certified dentists further verifies the clinical validity of our system.


title: Deep learning for low frequency extrapolation of multicomponent data in elastic full waveform inversion

abstract: Full waveform inversion (FWI) strongly depends on an accurate starting model to succeed. This is particularly true in the elastic regime: The cycle-skipping phenomenon is more severe in elastic FWI compared to acoustic FWI, due to the short S-wave wavelength. In this paper, we extend our work on extrapolated FWI (EFWI) by proposing to synthesize the low frequencies of multi-component elastic seismic records, and use those "artificial" low frequencies to seed the frequency sweep of elastic FWI. Our solution involves deep learning: we separately train the same convolutional neural network (CNN) on two training datasets, one with vertical components and one with horizontal components of particle velocities, to extrapolate the low frequencies of elastic data. The architecture of this CNN is designed with a large receptive field, by either large convolutional kernels or dilated convolution. Numerical examples on the Marmousi2 model show that the 2-4Hz low frequency data extrapolated from band-limited data above 4Hz provide good starting models for elastic FWI of P-wave and S-wave velocities. Additionally, we study the generalization ability of the proposed neural network over different physical models. For elastic test data, collecting the training dataset by elastic simulation shows better extrapolation accuracy than acoustic simulation, i.e., a smaller generalization gap.


title: The Mystery of Bloom-Gilman Duality: A Light Front Holographic QCD Perspective

abstract: Light front wave functions motivated by holographic constructions are used to study Bloom-Gilman duality of deep inelastic scattering. Separate expressions for structure functions in terms of quark and hadronic degrees of freedom are presented, with a goal of relating the two expressions. A two-parton model is defined and resonance transition form factors are computed using previously derived light front wave functions. A new form of global duality is derived from the valence quark-number sum rule. Using a complete set of hadronic states is necessary for this new global duality to be achieved. Previous original work does not provide such a set. This is remedied by amending the model to include a longitudinal confining potential, and the resulting complete set is sufficient to carry out the study of Bloom-Gilman duality. Expressions for transition form factors are obtained and all are shown to fall asymptotically as 1/Q2. The Feynman mechanism dominates the asymptotic behavior of the model. These transition form factors are used to assess the validity of the global and local duality sum rules, with the result that both neither are satisfied. Evaluations of the hadronic expression for q(x,Q2) provide more details about this lack. This result shows that the observed validity of both global and local forms of duality for deep inelastic scattering must be related to a feature of QCD that is deeper than completeness. Our simple present model suggests a prediction that Bloom-Gilman duality would not be observed if deep inelastic scattering experiments were to be made on the pion. The underlying origin of the duality phenomenon in deep inelastic scattering is deeply buried within the confinement aspects of QCD, and remains a mystery.


title: Climbing LP Algorithms

abstract: NP (search) problems allow easy correctness tests for solutions. Climbing algorithms allow also easy assessment of how close to yielding the correct answer is the configuration at any stage of their run. This offers a great flexibility, as how sensible is any deviation from the standard procedures can be instantly assessed. An example is the Dual Matrix Algorithm (DMA) for linear programming, variations of which were considered by A.Y. Levin in 1965 and by Yamnitsky and myself in 1982. It has little sensitivity to numerical errors and to the number of inequalities. It offers substantial flexibility and, thus, potential for further developments.


title: Verifying a Cruise Control System using Simulink and SpaceEx

abstract: This article aims to provide a simple step-by-step guide highlighting the steps needed to verify a control system with formal verification tools. Starting from a description of the physical system and a control objective in natural language, we design the plant and the controller, we use Simulink for simulation and we employ a reachability analysis tool, SpaceEx, for formal verification.


title: On the Structure of the Generalized Group of Units

abstract: Let $R$ be a finite commutative ring with identity and $U(R)$ be its group of units. In 2005, El-Kassar and Chehade presented a ring structure for $U(R)$ and as a consequence they generalized this group of units to the generalized group of units $U^{k}\left( R\right) $ defined iteratively as the group of the units of $U^{k-1}(R)$, with $U^{1}\left( R\right) =U(R) $. In this paper, we examine the structure of this group, when $R=\mathbb{Z}_{n}.$ We find a decomposition of $U^{k}\left(\mathbb{Z}_{n}\right)$ as a direct product of cyclic groups for the general case of any $k$, and we study when these groups are boolean and trivial. We also show that this decomposition structure is directly related to the Pratt Tree primes.


title: Indefinite Sturm-Liouville operators in polar form

abstract: We consider the indefinite Sturm-Liouville differential expression \[\mathfrak{a}(f) := - \frac{1}{w}\left( \frac{1}{r} f' \right)',\] where $\mathfrak{a}$ is defined on a finite or infinite open interval $I$ with $0\in I$ and the coefficients $r$ and $w$ are locally summable and such that $r(x)$ and $(\operatorname{sgn} x) w(x)$ are positive a.e. on $I$. With the differential expression $\mathfrak{a}$ we associate a nonnegative self-adjoint operator $A$ in the Krein space $L^2_w(I)$, which is viewed as a coupling of symmetric operators in Hilbert spaces related to the intersections of $I$ with the positive and the negative semi-axis. For the operator $A$ we derive conditions in terms of the coefficients $w$ and $r$ for the existence of a Riesz basis consisting of generalized eigenfunctions of $A$ and for the similarity of $A$ to a self-adjoint operator in a Hilbert space $L^2_{|w|}(I)$. These results are obtained as consequences of abstract results about the regularity of critical points of nonnegative self-adjoint operators in Krein spaces, which are couplings of two symmetric operators acting in Hilbert spaces.


title: A selective review on calibration information from similar studies based on parametric likelihood or empirical likelihood

abstract: In multi-center clinical trials, due to various reasons, the individual-level data are strictly restricted to be assessed publicly. Instead, the summarized information is widely available from published results. With the advance of computational technology, it has become very common in data analyses to run on hundreds or thousands of machines simultaneous, with the data distributed across those machines and no longer available in a single central location. How to effectively assemble the summarized clinical data information or information from each machine in parallel computation has become a challenging task for statisticians and computer scientists. In this paper, we selectively review some recently-developed statistical methods, including communication efficient distributed statistical inference, and renewal estimation and incremental inference, which can be regarded as the latest development of calibration information methods in the era of big data. Even though those methods were developed in different fields and in different statistical frameworks, in principle, they are asymptotically equivalent to those well known methods developed in meta analysis. Almost no or little information is lost compared with the case when full data are available. As a general tool to integrate information, we also review the generalized method of moments and estimating equations approach by using empirical likelihood method.


title: A Zonal Volt/VAR Control Mechanism for High PV Penetration Distribution Systems

abstract: This paper presents a zonal Volt/VAR control scheme that coordinates Photovoltaic (PV) inverters for providing voltage regulation on 3-phase unbalanced distribution feeders. Voltage sensitivity studies are conducted to uncover the dependency between nodal voltage changes and the reactive power injections at nodes with smart PV inverters. A fast incremental clustering method is developed to divide the distribution feeder into weakly coupled zones based on correlations of nodal voltage sensitivities. Because each zone is weakly coupled, voltage of each zone can be controlled independently. Thus, in each zone, a rule-based voltage controller will dispatch PV smart inverters to provide reactive power control for correcting the over/under voltages. An actual distribution feeder in North Carolina is used as a test bed. Simulation results show that the proposed zonal based Volt/VAR control mechanism can maintain the voltage in the distribution system within limits and solves faster than the centralized controller.


title: Rank of Near Uniform Matrices

abstract: A central question in random matrix theory is universality. When an emergent phenomena is observed from a large collection of chosen random variables it is natural to ask if this behavior is specific to the chosen random variable or if the behavior occurs for a larger class of random variables. The rank statistics of random matrices chosen uniformly from $\operatorname{Mat}(\mathbf{F}_q)$ over a finite field are well understood. The universality properties of these statistics are not yet fully understood however. Recently Wood [39] and Maples [26] considered a natural requirement where the random variables are not allowed to be too close to constant and they showed that the rank statistics match with the uniform model up to an error of type $e^{-cn}$. In this paper we explore a condition called near uniform, under which we are able to prove tighter bounds $q^{-cn}$ on the asymptotic convergence of the rank statistics. Our method is completely elementary, and allows for a small number of the entries to be deterministic, and for the entries to not be identically distributed so long as they are independent. More importantly, the method also extends to near uniform symmetric, alternating matrices. Our method also applies to two models of perturbations of random matrices sampled uniformly over $\operatorname{GL}_n(\mathbf{F}_q)$: subtracting the identity or taking a minor of a uniformly sampled invertible matrix.


title: The wake of a transversely oscillating circular cylinder in a flowing soap film at low Reynolds number

abstract: An inclined gravity-driven soap film channel was used to study the wake patterns formed behind a transversely oscillating cylinder at $Re =235 \pm 14$. The natural frequency of vortex shedding from a stationary cylinder, $f_{\text{St}}$, was used to identify the oscillation frequencies of interest. The (dimensionless) frequency, $f^*=f/{f_{\text{St}}}$, and amplitude, $A^*=A/D$, of the cylinder's motion was varied over a large portion of the fundamental synchronization region (i.e., for $f^* \approx 1$), and a `map' of wake patterns was constructed in ${(f^*, A^*)}$ space. Lock-on between the frequency of the cylinder's motion and the dominant frequency of the resulting vortex wake was observed for a large range of this parameter space, predominantly manifested as synchronized `2S' and `2P' wake modes. Synchronized `P+S', `2T', and `transitional' wakes were also found in smaller regions of parameter space. Unsynchronized `coalescing' and `perturbed von~\karman' wakes were observed as the oscillation frequency became sufficiently different from $f_{\text_{St}}$. The wake patterns and vortex formation processes found in this study, particularly for the 2P mode wakes, are more similar to those observed by \cite{Williamson1988} in three-dimensional experiments than those found by \cite{Leontini2006} in two-dimensional simulations, despite the physical constraint from the soap film that prevents three-dimensional effects in the wake.


title: Achievable Rates for Binary Two-hop Channel with Energy Harvesting Relay and Finite Battery

abstract: We study the problem of joint information and energy transfer in a binary two-hop channel with an energy harvesting relay. We consider a finite battery size at the relay and energy loss in transmitting energy. In other words, to be able to send an energy-contained symbol, the relay must receive multiple energy-contained symbols. Thus, we face a kind of channel with memory. We model the energy saved in the battery as the channel state with the challenge that the receiver does not know the channel state. We propose two different achievable schemes, the first one is based on state-dependent superposition coding and the second one is based on the equivalent timing channel approach. Both of our schemes are based on block Markov coding and backward decoding techniques. Due to these two approaches, we find achievable rates with a single-letter expression for the model.


title: Task Offloading and Resource Allocation with Multiple CAPs and Selfish Users

abstract: In this work, we consider a multi-user mobile edge computing system with multiple computing access points (CAPs). Each mobile user has multiple dependent tasks that must be processed in a round-by-round schedule. In every round, a user may process their individual task locally, or choose to offload their task to one of the $M$ CAPs or the remote cloud server, in order to possibly reduce their processing cost. We aim to jointly optimize the offloading decisions of the users and the resource allocation decisions for each CAP over a global objective function, defined as a weighted sum of total energy consumption and the round time. We first present a centralized heuristic solution, termed MCAP, where the original problem is relaxed to a semi-definite program (SDP) to probabilistically generate the offloading decision. Then, recognizing that the users often exhibit selfish behavior to reduce their individual cost, we propose a game-theoretical approach, termed MCAP-NE, which allows us to compute a Nash Equilibrium (NE) through a finite improvement method starting from the previous SDP solution. This approach leads to a solution from which the users have no incentive to deviate, with substantially reduced NE computation time. In simulation, we compare the system cost of the NE solution with those of MCAP, MCAP-NE, a random mapping, and the optimal solution, showing that our NE solution attains near optimal performance under a wide set of parameter settings, as well as demonstrating the advantages of using MCAP to produce the initial point for MCAP-NE.


title: Simulating Effective QED on Quantum Computers

abstract: In recent years simulations of chemistry and condensed materials has emerged as one of the preeminent applications of quantum computing, offering an exponential speedup for the solution of the electronic structure for certain strongly correlated electronic systems. To date, most treatments have ignored the question of whether relativistic effects, which are described most generally by quantum electrodynamics (QED), can also be simulated on a quantum computer in polynomial time. Here we show that effective QED, which is equivalent to QED to second order in perturbation theory, can be simulated in polynomial time under reasonable assumptions while properly treating all four components of the wavefunction of the fermionic field. In particular, we provide a detailed analysis of such simulations in position and momentum basis using Trotter-Suzuki formulas. We find that the number of $T$-gates needed to perform such simulations on a $3D$ lattice of $n_s$ sites scales at worst as $O(n_s^3/\epsilon)^{1+o(1)}$ in the thermodynamic limit for position basis simulations and $O(n_s^{4+2/3}/\epsilon)^{1+o(1)}$ in momentum basis. We also find that qubitization scales slightly better with a worst case scaling of $\widetilde{O}(n_s^{2+2/3}/\epsilon)$ for lattice eQED and complications in the prepare circuit leads to a slightly worse scaling in momentum basis of $\widetilde{O}(n_s^{5+2/3}/\epsilon)$. We further provide concrete gate counts for simulating a relativistic version of the uniform electron gas that show challenging problems can be simulated using fewer than $10^{13}$ non-Clifford operations and also provide a detailed discussion of how to prepare multi-reference configuration interaction states in effective QED which can provide a reasonable initial guess for the ground state. Finally, we estimate the planewave cutoffs needed to accurately simulate heavy elements such as gold.


title: Space-Time Analyticity of Weak Solutions to Semilinear Parabolic Systems with Variable Coefficients

abstract: Analytic smooth solutions of a general, strongly parabolic semi-linear Cauchy problem of $2m$-th order in $\mathbb{R}^N\times (0,T)$ with analytic coefficients (in space and time variables) and analytic initial data (in space variables) are investigated. They are expressed in terms of holomorphic continuation of global (weak) solutions to the system valued in a suitable Besov interpolation space of $B^{s;p,p}$-type at every time moment $t\in [0,T]$. Given $0 < T'< T\leq \infty$, it is proved that any $B^{s;p,p}$-type solution $u\colon \mathbb{R}^N\times (0,T)\to \mathbb{C}^M$ with analytic initial data possesses a bounded holomorphic continuation $u(x + \mathrm{i}y, \sigma + \mathrm{i}\tau)$ into a complex domain in $\mathbb{C}^N\times \mathbb{C}$ defined by $(x,\sigma)\in \mathbb{R}^N\times (T',T)$, $|y| < A'$ and $|\tau | < B'$, where $A', B'> 0$ are constants depending upon~$T'$. The proof uses the extension of a weak solution to a $B^{s;p,p}$-type solution in a domain in $\mathbb{C}^N\times \mathbb{C}$, such that this extension satisfies the Cauchy-Riemann equations. The holomorphic extension is obtained with a help from holomorphic semigroups and maximal regularity theory for parabolic problems in Besov interpolation spaces of $B^{s;p,p}$-type imbedded (densely and continuously) into an $L^p$-type Lebesgue space. Applications include risk models for European options in Mathematical Finance.


title: Cosmological inference from emulator based halo model I: Validation tests with HSC and SDSS mock catalogs

abstract: We present validation tests of emulator-based halo model method for cosmological parameter inference, assuming hypothetical measurements of the projected correlation function of galaxies, $w_{\rm p}(R)$, and the galaxy-galaxy weak lensing, $\Delta\!\Sigma(R)$, from the spectroscopic SDSS galaxies and the Hyper Suprime-Cam Year1 (HSC-Y1) galaxies. To do this, we use \textsc{Dark Emulator} developed in Nishimichi et al. based on an ensemble of $N$-body simulations, which is an emulation package enabling a fast, accurate computation of halo clustering quantities for flat-geometry $w$CDM cosmologies. Adopting the halo occupation distribution, the emulator allows us to obtain model predictions of $\Delta\!\Sigma$ and $w_{\rm p}$ for the SDSS-like galaxies at a few CPU seconds for an input set of parameters. We present performance and validation of the method by carrying out Markov Chain Monte Carlo analyses of the mock signals measured from a variety of mock catalogs that mimic the SDSS and HSC-Y1 galaxies. We show that the halo model method can recover the underlying true cosmological parameters to within the 68\% credible interval, except for the mocks including the assembly bias effect (although we consider the unrealistically-large amplitude of assembly bias effect). Even for the assembly bias mock, we demonstrate that the cosmological parameters can be recovered {\it if} the analysis is restricted to scales $R\gtrsim 10~h^{-1}{\rm Mpc}$. We also show that, by using a single population of source galaxies to infer the relative strengths of $\Delta\!\Sigma$ for multiple lens samples at different redshifts, the joint probes method allows for self-calibration of photometric redshift errors and multiplicative shear bias. Thus we conclude that the emulator-based halo model method can be safely applied to the HSC-Y1 dataset, achieving a precision of $\sigma(S_8)\simeq 0.04$.


title: Entanglement harvesting in the presence of a reflecting boundary

abstract: We study, in the framework of the entanglement harvesting protocol, the entanglement harvesting of both a pair of inertial and uniformly accelerated detectors locally interacting with vacuum massless scalar fields subjected to a perfectly reflecting plane boundary. We find that the presence of the boundary generally degrades the harvested entanglement when two detectors are very close to the boundary. However, when the distance between detectors and the boundary becomes comparable to the interaction duration parameter, the amount of the harvested entanglement approaches a peak, which even goes beyond that without a boundary. Remarkably, the parameter space of the detectors' separation and the magnitude of acceleration that allows entanglement harvesting to occur is enlarged due to the presence of the boundary. In this sense, the boundary plays a double-edged role on entanglement harvesting, degrading in general the harvested entanglement while enlarging the entanglement harvesting-achievable parameter space. A comparison of three different acceleration scenarios of the detectors with respect to the boundary, i.e., parallel, anti-parallel and mutually perpendicular acceleration, shows that the phenomenon of entanglement harvesting crucially depends on the acceleration, the separation between two detectors and the detectors' distance from the boundary.


title: Energy Performance Analysis of Distributed Renewables: Pacific Northwest Smart Grid Demonstration

abstract: The Pacific Northwest Smart Grid Demonstration was an electricity grid modernization project conducted in the Northwest U.S. This paper presents the analysis of renewable generation at the Renewable Energy Park located in the City of Ellensburg, WA. The community energy park concept is an intriguing model for community investment in renewable resources,but the lessons in this paper should be considered.


title: Dynamic system optimal traffic assignment with atomic users: Convergence and stability

abstract: In this study, we analyse the convergence and stability of dynamic system optimal (DSO) traffic assignment with fixed departure times. We first formulate the DSO traffic assignment problem as a strategic game wherein atomic users select routes that minimise their marginal social costs, called a 'DSO game'. By utilising the fact that the DSO game is a potential game, we prove that a globally optimal state is a stochastically stable state under the logit response dynamics, and the better/best response dynamics converges to a locally optimal state. Furthermore, as an application of DSO assignment, we examine characteristics of the evolutionary implementation scheme of marginal cost pricing. Through theoretical comparison with a fixed pricing scheme, we found the following properties of the evolutionary implementation scheme: (i) the total travel time decreases smoother to an efficient traffic state as congestion externalities are perfectly internalised; (ii) a traffic state would reach a more efficient state as the globally optimal state is stabilised. Numerical experiments also suggest that these properties make the evolutionary scheme robust in the sense that they prevent a traffic state from going to worse traffic states with high total travel times.


title: Multi-task Retrieval for Knowledge-Intensive Tasks

abstract: Retrieving relevant contexts from a large corpus is a crucial step for tasks such as open-domain question answering and fact checking. Although neural retrieval outperforms traditional methods like tf-idf and BM25, its performance degrades considerably when applied to out-of-domain data. Driven by the question of whether a neural retrieval model can be universal and perform robustly on a wide variety of problems, we propose a multi-task trained model. Our approach not only outperforms previous methods in the few-shot setting, but also rivals specialised neural retrievers, even when in-domain training data is abundant. With the help of our retriever, we improve existing models for downstream tasks and closely match or improve the state of the art on multiple benchmarks.


title: A Two Stage Adaptive Metropolis Algorithm

abstract: We propose a new sampling algorithm combining two quite powerful ideas in the Markov chain Monte Carlo literature -- adaptive Metropolis sampler and two-stage Metropolis-Hastings sampler. The proposed sampling method will be particularly very useful for high-dimensional posterior sampling in Bayesian models with expensive likelihoods. In the first stage of the proposed algorithm, an adaptive proposal is used based on the previously sampled states and the corresponding acceptance probability is computed based on an approximated inexpensive target density. The true expensive target density is evaluated while computing the second stage acceptance probability only if the proposal is accepted in the first stage. The adaptive nature of the algorithm guarantees faster convergence of the chain and very good mixing properties. On the other hand, the two-stage approach helps in rejecting the bad proposals in the inexpensive first stage, making the algorithm computationally efficient. As the proposals are dependent on the previous states the chain loses its Markov property, but we prove that it retains the desired ergodicity property. The performance of the proposed algorithm is compared with the existing algorithms in two simulated and two real data examples.


title: On the heat equation with drift in $L_{d+1}$

abstract: In this paper we deal with the heat equation with drift in $L_{d+1}$. Basically, we prove that, if the free term is in $L_{q}$ with high enough $q$, then the equation is uniquely solvable in a rather unusual class of functions such that $\partial_{t}u, D^{2}u\in L_{p}$ with $p<d+1$ and $Du\in L_{q}$.


title: Simple close curve magnetization and application to Bellman's lost in the forest problem

abstract: In this paper we introduce and develop the notion of simple close curve magnetization. We provide an application to Bellman's lost in the forest problem assuming special geometric conditions between the hiker and the boundary of the forest.


title: WARP: Word-level Adversarial ReProgramming

abstract: Transfer learning from pretrained language models recently became the dominant approach for solving many NLP tasks. A common approach to transfer learning for multiple tasks that maximize parameter sharing trains one or more task-specific layers on top of the language model. In this paper, we present an alternative approach based on adversarial reprogramming, which extends earlier work on automatic prompt generation. Adversarial reprogramming attempts to learn task-specific word embeddings that, when concatenated to the input text, instruct the language model to solve the specified task. Using up to 25K trainable parameters per task, this approach outperforms all existing methods with up to 25M trainable parameters on the public leaderboard of the GLUE benchmark. Our method, initialized with task-specific human-readable prompts, also works in a few-shot setting, outperforming GPT-3 on two SuperGLUE tasks with just 32 training samples.


title: Generative Max-Mahalanobis Classifiers for Image Classification, Generation and More

abstract: Joint Energy-based Model (JEM) of Grathwohl et al. shows that a standard softmax classifier can be reinterpreted as an energy-based model (EBM) for the joint distribution p(x,y); the resulting model can be optimized to improve calibration, robustness, and out-of-distribution detection, while generating samples rivaling the quality of recent GAN-based approaches. However, the softmax classifier that JEM exploits is inherently discriminative and its latent feature space is not well formulated as probabilistic distributions, which may hinder its potential for image generation and incur training instability. We hypothesize that generative classifiers, such as Linear Discriminant Analysis (LDA), might be more suitable for image generation since generative classifiers model the data generation process explicitly. This paper therefore investigates an LDA classifier for image classification and generation. In particular, the Max-Mahalanobis Classifier (MMC), a special case of LDA, fits our goal very well. We show that our Generative MMC (GMMC) can be trained discriminatively, generatively, or jointly for image classification and generation. Extensive experiments on multiple datasets show that GMMC achieves state-of-the-art discriminative and generative performances, while outperforming JEM in calibration, adversarial robustness, and out-of-distribution detection by a significant margin. Our source code is available at https://github.com/sndnyang/GMMC.


title: Intent Classification and Slot Filling for Privacy Policies

abstract: Understanding privacy policies is crucial for users as it empowers them to learn about the information that matters to them. Sentences written in a privacy policy document explain privacy practices, and the constituent text spans convey further specific information about that practice. We refer to predicting the privacy practice explained in a sentence as intent classification and identifying the text spans sharing specific information as slot filling. In this work, we propose PolicyIE, an English corpus consisting of 5,250 intent and 11,788 slot annotations spanning 31 privacy policies of websites and mobile applications. PolicyIE corpus is a challenging real-world benchmark with limited labeled examples reflecting the cost of collecting large-scale annotations from domain experts. We present two alternative neural approaches as baselines, (1) intent classification and slot filling as a joint sequence tagging and (2) modeling them as a sequence-to-sequence (Seq2Seq) learning task. The experiment results show that both approaches perform comparably in intent classification, while the Seq2Seq method outperforms the sequence tagging approach in slot filling by a large margin. We perform a detailed error analysis to reveal the challenges of the proposed corpus.


title: Discourse-level Relation Extraction via Graph Pooling

abstract: The ability to capture complex linguistic structures and long-term dependencies among words in the passage is essential for discourse-level relation extraction (DRE) tasks. Graph neural networks (GNNs), one of the methods to encode dependency graphs, have been shown effective in prior works for DRE. However, relatively little attention has been paid to receptive fields of GNNs, which can be crucial for cases with extremely long text that requires discourse understanding. In this work, we leverage the idea of graph pooling and propose to use pooling-unpooling framework on DRE tasks. The pooling branch reduces the graph size and enables the GNNs to obtain larger receptive fields within fewer layers; the unpooling branch restores the pooled graph to its original resolution so that representations for entity mention can be extracted. We propose Clause Matching (CM), a novel linguistically inspired graph pooling method for NLP tasks. Experiments on two DRE datasets demonstrate that our models significantly improve over baselines when modeling long-term dependencies is required, which shows the effectiveness of the pooling-unpooling framework and our CM pooling method.


title: Self-Adaptive Systems: A Systematic Literature Review Across Categories and Domains

abstract: Context: Championed by IBM's vision of autonomic computing paper in 2003, the autonomic computing research field has seen increased research activity over the last 20 years. Several conferences and workshops have been established and have contributed to the autonomic computing knowledge base in search of a new kind of system -- a self-adaptive system (SAS). These systems are characterized by being context-aware and can act on that awareness. The actions carried out could be on the system or on the context (or environment). The underlying goal of a SAS is the sustained achievement of its goals despite changes in its environment. Objective: Despite a number of literature reviews on specific aspects of SASs ranging from their requirements to quality attributes, we lack a systematic understanding of the current state of the art. Method: This paper contributes a systematic literature review into self-adaptive systems using the dblp computer science bibliography as a database. We filtered the records systematically in successive steps to arrive at 293 relevant papers. Each paper was critically analyzed and categorized into an attribute matrix. This matrix consisted of five categories, with each category having multiple attributes. The attributes of each paper, along with the summary of its contents formed the basis of the literature review that spanned 30 years. Results: We characterize the maturation process of the research area from theoretical papers over practical implementations to more holistic and generic approaches, frameworks, and exemplars, applied to areas such as networking, web services, and robotics, with much of the recent work focusing on IoT and IaaS. Conclusion: While there is an ebb and flow of application domains, domains like bio-inspired approaches, security, and cyber physical systems showed promise to grow heading into the 2020s.


title: Transport in helical Luttinger liquids in the fractional quantum Hall regime

abstract: Domain walls in fractional quantum Hall ferromagnets are gapless helical one-dimensional channels formed at the boundaries of topologically distinct quantum Hall (QH) liquids. Na\"{i}vely, these helical domain walls (hDWs) constitute two counter-propagating chiral states with opposite spins. Coupled to an s-wave superconductor, helical channels are expected to lead to topological superconductivity with high order non-Abelian excitations. Here we investigate transport properties of hDWs in the $\nu=2/3$ fractional QH regime. Experimentally we found that current carried by hDWs is substantially smaller than the prediction of the na\"{i}ve model. Luttinger liquid theory of the system reveals redistribution of currents between quasiparticle charge, spin and neutral modes, and predicts the reduction of the hDW current. Inclusion of spin-non-conserving tunneling processes reconciles theory with experiment. The theory confirms emergence of spin modes required for the formation of fractional topological superconductivity.


title: Formalizing Hall's Marriage Theorem in Lean

abstract: We formalize Hall's Marriage Theorem in the Lean theorem prover for inclusion in mathlib, which is a community-driven effort to build a unified mathematics library for Lean. One goal of the mathlib project is to contain all of the topics of a complete undergraduate mathematics education. We provide three presentations of the main theorem statement: in terms of indexed families of finite sets, of relations on types, and of matchings in bipartite graphs. We also formalize a version of K\H{o}nig's lemma (in terms of inverse limits) to boost the theorem to the case of countably infinite index sets. We give a description of the design of the recent mathlib library for simple graphs, and we also give a necessary and sufficient condition for a simple graph to carry a function.


title: VLBI observations of supernova PTF11qcj: Direct constraints on the size of the radio ejecta

abstract: We present High Sensitivity Array (HSA) and enhanced Multi Element Remotely Linked Interferometer Network (eMERLIN) observations of the radio-loud broad-lined type Ic supernova PTF11qcj obtained $\sim7.5$ years after the explosion. Previous observations of this supernova at 5.5 yrs since explosion showed a double-peaked radio light curve accompanied by a detection in the X-rays, but no evidence for broad H$\alpha$ spectral features. The Very Long Baseline Interferometry (VLBI) observations presented here show that the PTF11qcj GHz radio ejecta remains marginally resolved at the sub-milliarcsecond level at $\approx 7.5$ yrs after the explosion, pointing toward a non-relativistic expansion. Our VLBI observations thus favor a scenario in which the second peak of the PTF11qcj radio light curve is related to strong interaction of the supernova ejecta with a circumstellar medium of variable density, rather than to the emergence of an off-axis jet. Continued VLBI monitoring of PTF11qcj in the radio may strengthen further this conclusion.


title: Universality of Weyl Unitaries

abstract: Weyl's unitary matrices, which were introduced in Weyl's 1927 paper on group theory and quantum mechanics, are $p\times p$ unitary matrices given by the diagonal matrix whose entries are the $p$-th roots of unity and the cyclic shift matrix. Weyl's unitaries, which we denote by $\mathfrak u$ and $\mathfrak v$, satisfy $\mathfrak u^p=\mathfrak v^p=1_p$ (the $p\times p$ identity matrix) and the commutation relation $\mathfrak u\mathfrak v=\zeta \mathfrak v\mathfrak u$, where $\zeta$ is a primitive $p$-th root of unity. We prove that Weyl's unitary matrices are universal in the following sense: if $u$ and $v$ are any $d\times d$ unitary matrices such that $u^p= v^p=1_d$ and $ u v=\zeta vu$, then there exists a unital completely positive linear map $\phi:\mathcal M_p(\mathbb C)\rightarrow\mathcal M_d(\mathbb C)$ such that $\phi(\mathfrak u)= u$ and $\phi(\mathfrak v)=v$. We also show, moreover, that any two pairs of $p$-th order unitary matrices that satisfy the Weyl commutation relation are completely order equivalent. When $p=2$, the Weyl matrices are two of the three Pauli matrices from quantum mechanics. It was recently shown that $g$-tuples of Pauli-Weyl-Brauer unitaries are universal for all $g$-tuples of anticommuting selfadjoint unitary matrices; however, we show here that the analogous result fails for positive integers $p>2$. Finally, we show that the Weyl matrices are extremal in their matrix range, using recent ideas from noncommutative convexity theory.


title: Sensei: Self-Supervised Sensor Name Segmentation

abstract: A sensor name, typically an alphanumeric string, encodes the key context (e.g., function and location) of a sensor needed for deploying smart building applications. Sensor names, however, are curated in a building vendor-specific manner using different structures and vocabularies that are often esoteric. They thus require tremendous manual effort to annotate on a per-building basis; even to just segment these sensor names into meaningful chunks. In this paper, we propose a fully automated self-supervised framework, Sensei, which can learn to segment sensor names without any human annotation. Specifically, we employ a neural language model to capture the underlying sensor naming structure and then induce self-supervision based on information from the language model to build the segmentation model. Extensive experiments on five real-world buildings comprising thousands of sensors demonstrate the superiority of Sensei over baseline methods.


title: Formation and evolution of protostellar accretion discs. I. Angular-momentum budget, gravitational self-regulation, and numerical convergence

abstract: We investigate the formation and early evolution of a protostellar disc from a magnetized pre-stellar core using non-ideal magnetohydrodynamic (MHD) simulations including ambipolar diffusion and Ohmic dissipation. The dynamical contraction of the pre-stellar core ultimately leads to the formation of a first hydrostatic core, after ambipolar diffusion decouples the magnetic field from the predominantly neutral gas. The hydrostatic core accumulates angular momentum from the infalling material, evolving into a rotationally supported torus; this `first hydrostatic torus' then forms an accreting protostar and a rotationally supported disc. The disc spreads out by gravitational instability, reaching $\sim$30 au in diameter at $\sim$3 kyr after protostar formation. The total mass and angular momentum of the protostar-disc system are determined mainly by accretion of gas from an infalling pseudo-disc, which has low specific angular momentum because of magnetic braking; their removal from the protostar-disc system by outflow and disc magnetic braking are negligible, in part because the magnetic field is poorly coupled there. The redistribution of angular momentum within the protostar-disc system is facilitated mainly by gravitational instability; this allows formation of relatively large discs even when the specific angular momentum of infalling material is low. We argue that such discs should remain marginally unstable as they grow (with Toomre $Q\sim 1$-$2$), an idea that is broadly consistent with recent observational estimates for Class 0/I discs. We discuss the numerical convergence of our results, and show that properly treating the inner boundary condition is crucial for achieving convergence at an acceptable computational cost.


title: Audio Content Analysis

abstract: Preprint for a book chapter introducing Audio Content Analysis. With a focus on Music Information Retrieval systems, this chapter defines musical audio content, introduces the general process of audio content analysis, and surveys basic approaches to audio content analysis. The various tasks in Audio Content Analysis are categorized into three classes: music transcription, music performance analysis, and music identification and categorization. The examples for music transcription systems include music key detection, fundamental frequency detection, and music structure detection. Music performance analysis systems feature an overview of beat and tempo detection approaches as well as music performance assessment. The covered music classification systems are audio fingerprinting, music genre classification, and music emotion recognition. The chapter concludes with a discussion and current challenges in the field and a speculation on future perspectives.


title: NeurIPS 2020 EfficientQA Competition: Systems, Analyses and Lessons Learned

abstract: We review the EfficientQA competition from NeurIPS 2020. The competition focused on open-domain question answering (QA), where systems take natural language questions as input and return natural language answers. The aim of the competition was to build systems that can predict correct answers while also satisfying strict on-disk memory budgets. These memory budgets were designed to encourage contestants to explore the trade-off between storing retrieval corpora or the parameters of learned models. In this report, we describe the motivation and organization of the competition, review the best submissions, and analyze system predictions to inform a discussion of evaluation for open-domain QA.


title: $\mathcal{L}_1$ Adaptive Control for Switching Reference Systems: Application to Flight Control

abstract: This paper presents a framework for the design and analysis of an $\mathcal{L}_1$ adaptive controller with a switching reference system. The use of a switching reference system allows the desired behavior to be scheduled across the operating envelope, which is often required in aerospace applications. The analysis uses a switched reference system that assumes perfect knowledge of uncertainties and uses a corresponding non-adaptive controller. Provided that this switched reference system is stable, it is shown that the closed-loop system with unknown parameters and disturbances and the $\mathcal{L}_1$ adaptive controller can behave arbitrarily close to this reference system. Simulations of the short period dynamics of a transport class aircraft during the approach phase illustrate the theoretical results.


title: Faults in Deep Reinforcement Learning Programs: A Taxonomy and A Detection Approach

abstract: A growing demand is witnessed in both industry and academia for employing Deep Learning (DL) in various domains to solve real-world problems. Deep Reinforcement Learning (DRL) is the application of DL in the domain of Reinforcement Learning (RL). Like any software systems, DRL applications can fail because of faults in their programs. In this paper, we present the first attempt to categorize faults occurring in DRL programs. We manually analyzed 761 artifacts of DRL programs (from Stack Overflow posts and GitHub issues) developed using well-known DRL frameworks (OpenAI Gym, Dopamine, Keras-rl, Tensorforce) and identified faults reported by developers/users. We labeled and taxonomized the identified faults through several rounds of discussions. The resulting taxonomy is validated using an online survey with 19 developers/researchers. To allow for the automatic detection of faults in DRL programs, we have defined a meta-model of DRL programs and developed DRLinter, a model-based fault detection approach that leverages static analysis and graph transformations. The execution flow of DRLinter consists in parsing a DRL program to generate a model conforming to our meta-model and applying detection rules on the model to identify faults occurrences. The effectiveness of DRLinter is evaluated using 15 synthetic DRLprograms in which we injected faults observed in the analyzed artifacts of the taxonomy. The results show that DRLinter can successfully detect faults in all synthetic faulty programs.


title: Sub-Gaussian Error Bounds for Hypothesis Testing

abstract: We interpret likelihood-based test functions from a geometric perspective where the Kullback-Leibler (KL) divergence is adopted to quantify the distance from a distribution to another. Such a test function can be seen as a sub-Gaussian random variable, and we propose a principled way to calculate its corresponding sub-Gaussian norm. Then an error bound for binary hypothesis testing can be obtained in terms of the sub-Gaussian norm and the KL divergence, which is more informative than Pinsker's bound when the significance level is prescribed. For $M$-ary hypothesis testing, we also derive an error bound which is complementary to Fano's inequality by being more informative when the number of hypotheses or the sample size is not large.


title: Coherent optical communications using coherence-cloned Kerr soliton microcombs

abstract: Dissipative Kerr soliton microcomb has been recognized as a promising on-chip multi-wavelength laser source for fiber optical communications, as its comb lines possess frequency and phase stability far beyond independent lasers. In the scenarios of coherent optical transmission and interconnect, a highly beneficial but rarely explored target is to re-generate a Kerr soliton microcomb at the receiver side as local oscillators that conserve the frequency and phase property of the incoming data carriers, so that to enable coherent detection with minimized optical and electrical compensations. Here, by using the techniques of pump laser conveying and two-point locking, we implement re-generation of a Kerr soliton microcomb that faithfully clones the frequency and phase coherence of another microcomb sent from 50 km away. Moreover, leveraging the coherence-cloned soliton microcombs as carriers and local oscillators, we demonstrate terabit coherent data interconnect, wherein traditional digital processes for frequency offset estimation is totally dispensed with, and carrier phase estimation is substantially simplified via slowed-down phase estimation rate per channel and joint phase estimation among multiple channels. Our work reveals that, in addition to providing a multitude of laser tones, regulating the frequency and phase of Kerr soliton microcombs among transmitters and receivers can significantly improve coherent communication in terms of performance, power consumption, and simplicity.


title: Divisors computing minimal log discrepancies on lc surfaces

abstract: Let $(X\ni x,B)$ be an lc surface germ. If $X\ni x$ is klt, we show that there exists a divisor computing the minimal log discrepancy of $(X\ni x,B)$ that is a Koll\'ar component of $X\ni x$. If $B\not=0$ or $X\ni x$ is not Du Val, we show that any divisor computing the minimal log discrepancy of $(X\ni x,B)$ is a potential lc place of $X\ni x$.


title: Frequency locking, Quasi periodicity and Chaos due to special relativistic effects

abstract: We study quasi periodic and frequency locked states that can occur in a sinusoidally driven linear harmonic oscillator in the special relativistic regime. We show how the shift in natural frequency of the oscillator with increasing relativistic effects leads to frequency locking or quasi periodicity and the chaotic states that arise due to the increasing non linearity. We find the same system can have multi-stable states in the presence of small damping. We also report an enhancement of chaos in the relativistic Henon-Heiles system.


title: An Electronic Textile Embedded Smart Cementitious Composite

abstract: Structural health monitoring (SHM) using self-sensing cement-based materials has been reported before, where nano-fillers have been incorporated in cementitious matrices as functional sensing elements. A percolation threshold is always required in order for conductive nano-fillers modified concrete to be useful for SHM. Nonetheless, the best pressure/strain sensitivity results achieved for any self-sensing cementitious matrix are <0.01 MPa-1. In this work, we introduce novel reduced graphene oxide (RGO) based electronic textile (e-textile) embedded in plain and polymer-binder-modified cementitious matrix for SHM applications. As a proof of concept, it was demonstrated that these coated fabric-based sensors can be successfully embedded within the cement-based structures, which are independent of any percolation threshold due to the interconnected fabric inside the host matrix. The piezo-resistive response was measured by applying direct and cyclic compressive loads (0.1 to 3.9 MPa). A pressure sensitivity of 1.5 MPa-1 and an ultra-high gauge factor of 2000 was obtained for the system of the self-sensing cementitious structure with embedded e-textiles. The sensitivity of this new system with embedded e-textile is many orders of magnitude higher than nanoparticle based self-sensing of cementitious composites. The manufactured e-textile sensors showed mechanical stability and functional durability over long-term cyclic compression tests of 1000 cycles.