title: Structures of Spurious Local Minima in $k$-means

abstract: $k$-means clustering is a fundamental problem in unsupervised learning. The problem concerns finding a partition of the data points into $k$ clusters such that the within-cluster variation is minimized. Despite its importance and wide applicability, a theoretical understanding of the $k$-means problem has not been completely satisfactory. Existing algorithms with theoretical performance guarantees often rely on sophisticated (sometimes artificial) algorithmic techniques and restricted assumptions on the data. The main challenge lies in the non-convex nature of the problem; in particular, there exist additional local solutions other than the global optimum. Moreover, the simplest and most popular algorithm for $k$-means, namely Lloyd's algorithm, generally converges to such spurious local solutions both in theory and in practice. In this paper, we approach the $k$-means problem from a new perspective, by investigating the structures of these spurious local solutions under a probabilistic generative model with $k$ ground truth clusters. As soon as $k=3$, spurious local minima provably exist, even for well-separated and balanced clusters. One such local minimum puts two centers at one true cluster, and the third center in the middle of the other two true clusters. For general $k$, one local minimum puts multiple centers at a true cluster, and one center in the middle of multiple true clusters. Perhaps surprisingly, we prove that this is essentially the only type of spurious local minima under a separation condition. Our results pertain to the $k$-means formulation for mixtures of Gaussians or bounded distributions. Our theoretical results corroborate existing empirical observations and provide justification for several improved algorithms for $k$-means clustering.


title: A Novel Numerical Index for Assessing Results of Frequency Response Analysis (FRA): an Experimental Study on Electrical Machines

abstract: This paper is an effort to evaluate the numerical indices for the assessment of the frequency response analysis (FRA) to provide a deeper understanding of their characteristics. The study introduces the indices in the literature and categorizes them into two groups. The results of an actual FRA experimental setup on an electrical machine are employed to examine the indices from the various aspects. The main features of the indices are extracted, and the better ones in various aspects are marked. Moreover, based on the observations from the experimental studies, a new numerical index for comparing FRA results is proposed in this paper to enhance FRA results assessment and interpretation in both electrical machines and transformers. The advantages of the proposed index are also shown by comparing it with the existing ones.


title: Unitarization of the Horocyclic Radon Transform on Homogeneous Trees

abstract: Following previous work in the continuous setup, we construct the unitarization of the horocyclic Radon transform on a homogeneous tree X and we show that it intertwines the quasi regular representations of the group of isometries of X on the tree itself and on the space of horocycles.


title: A posteriori error estimates of finite element methods by preconditioning

abstract: We present a framework that relates preconditioning with a posteriori error estimates in finite element methods. In particular, we use standard tools in subspace correction methods to obtain reliable and efficient error estimators. As a simple example, we recover the classical residual error estimators for the second order elliptic equations.


title: How massive are the superfluid cores in the Crab and Vela pulsars and why their glitch-events are accompanied with under and overshootings?

abstract: The Crab and Vela are well-studied glitching pulsars and the data obtained sofar should enable us to test the reliability of models of their internal structures. Very recently it was proposed that glitching pulsars are embedded in bimetric spacetime: their incompressible superfluid cores (SuSu-cores) are embedded in flat spacetime, whereas the ambient compressible and dissipative media are enclosed in Schwarzschild spacetime. In this letter we apply this model to the Crab and Vela pulsars and show that a newly born pulsar initially of $1.25\text{M}_{\odot}$ and an embryonic SuSu-core of $0.029\text{M}_{\odot}$ could evolve into a \ARZL{Crab-like}\ARZR pulsar after 1000 years and into a \ARZL{Vela-like}\ARZR pulsar 10000 years later to finally fade away as an invisible dark energy object after roughly 10 Myr. Based thereon we infer that the Crab and the Vela pulsars should have SuSu-cores of $0.15\text{M}_{\odot}$ and $0.55\text{M}_{\odot}$, respectively. Futhermore, the under- and overshootings phenomena observed to accompany the glitch events of the Vela pulsar are rather a common phenomenon of glitching pulsars that can be well-explained within the framework of bimetric spacetime.


title: The nucleon isovector tensor charge from lattice QCD using chiral fermions

abstract: In this work we present the isovector flavor combination for the nucleon tensor charge extracted from lattice QCD simulations using overlap fermions on $N_f=2+1$ domain-wall configurations. The pion mass dependence is studied using six valence quark masses, each reproducing a value for the pion mass in the valence sector between 147 and 330 MeV. We investigate and eliminate systematic uncertainties due to contamination by excited states, by employing several values for the source-sink separation that span from 1 fm to 1.6 fm. We apply a chiral extrapolation in the valence sector using a quadratic and a logarithmic term to fit the pion mass dependence, which describes well the lattice data. The lattice matrix element is renormalized non-perturbatively, and the final result is $g_T=1.096(30)$ in the $\overline{\rm MS}$ scheme at a renormalization scale of 2 GeV.


title: Non-existence of dead cores in fully nonlinear elliptic models

abstract: We investigate non-existence of nonnegative dead-core solutions for the problem $$|Du|^\gamma F(x, D^2u)+a(x)u^q = 0 \quad \mbox{in} \quad \Omega, \quad u=0 \quad \mbox{ on } \quad \partial\Omega.$$ Here $\Omega \subset \mathbb{R}^N$ is a bounded smooth domain, $F$ is a fully nonlinear elliptic operator, $a: \Omega \to \mathbb{R}$ is a sign-changing weight, $\gamma \geq 0$, and $0<q<\gamma+1$. We show that this problem has no non-trivial dead core solutions if either $q$ is close enough to $\gamma+1$ or the negative part of $a$ is sufficiently small. In addition, we obtain the existence and uniqueness of a positive solution under these conditions on $q$ and $a$. Our results extend previous ones established in the semilinear case, and are new even for the simple model $|D u(x)|^{\gamma} \mathrm{Tr}(\mathrm{A}(x) D^2 u(x)) + a(x)u^{q}(x) = 0$, where $\mathrm{A} \in C^0(\Omega;Sym(N))$ is a uniformly elliptic and non-negative matrix.


title: Gaussian Smoothen Semantic Features (GSSF) -- Exploring the Linguistic Aspects of Visual Captioning in Indian Languages (Bengali) Using MSCOCO Framework

abstract: In this work, we have introduced Gaussian Smoothen Semantic Features (GSSF) for Better Semantic Selection for Indian regional language-based image captioning and introduced a procedure where we used the existing translation and English crowd-sourced sentences for training. We have shown that this architecture is a promising alternative source, where there is a crunch in resources. Our main contribution of this work is the development of deep learning architectures for the Bengali language (is the fifth widely spoken language in the world) with a completely different grammar and language attributes. We have shown that these are working well for complex applications like language generation from image contexts and can diversify the representation through introducing constraints, more extensive features, and unique feature spaces. We also established that we could achieve absolute precision and diversity when we use smoothened semantic tensor with the traditional LSTM and feature decomposition networks. With better learning architecture, we succeeded in establishing an automated algorithm and assessment procedure that can help in the evaluation of competent applications without the requirement for expertise and human intervention.


title: Multi-item Non-truthful Auctions Achieve Good Revenue

abstract: We present a general framework for designing approximately revenue-optimal mechanisms for multi-item additive auctions, which applies to both truthful and non-truthful auctions. Given a (not necessarily truthful) single-item auction format $A$ satisfying certain technical conditions, we run simultaneous item auctions augmented with a personalized entry fee for each bidder that must be paid before the auction can be accessed. These entry fees depend only on the prior distribution of bidder types, and in particular are independent of realized bids. We bound the revenue of the resulting two-part tariff mechanism using a novel geometric technique that enables revenue guarantees for many common non-truthful auctions that previously had none. Our approach adapts and extends the duality framework of Cai et al [CDW16] beyond truthful auctions. Our framework can be used with many common auction formats, such as simultaneous first-price, simultaneous second-price, and simultaneous all-pay auctions. Our results for first price and all-pay are the first revenue guarantees of non-truthful mechanisms in multi-dimensional environments, addressing an open question in the literature [RST17]. If all-pay auctions are used, we prove that the resulting mechanism is also credible in the sense that the auctioneer cannot benefit by deviating from the stated mechanism after observing agent bids. This is the first static credible mechanism for multi-item additive auctions that achieves a constant factor of the optimal revenue. If second-price auctions are used, we obtain a truthful $O(1)$-approximate mechanism with fixed entry fees that are amenable to tuning via online learning techniques.


title: Investigating Simple Object Representations in Model-Free Deep Reinforcement Learning

abstract: We explore the benefits of augmenting state-of-the-art model-free deep reinforcement algorithms with simple object representations. Following the Frostbite challenge posited by Lake et al. (2017), we identify object representations as a critical cognitive capacity lacking from current reinforcement learning agents. We discover that providing the Rainbow model (Hessel et al.,2018) with simple, feature-engineered object representations substantially boosts its performance on the Frostbite game from Atari 2600. We then analyze the relative contributions of the representations of different types of objects, identify environment states where these representations are most impactful, and examine how these representations aid in generalizing to novel situations.


title: Performance Limits of Differential Power Processing

abstract: This paper investigates the performance limits of differential power processing (DPP) and presents quantitative and systematic design guidelines for the selection and comparison of DPP topologies. A stochastic model is developed to evaluate the expected power losses of a variety of DPP topologies with probabilistic load distribution. The expected losses of several DPP topologies are derived and compared against traditional dc-dc converters to reveal their performance limits. The impacts of the load distribution and load scale on the expected losses are investigated. The theoretical models are verified with SPICE simulations and experimental results.


title: Strong colorings over partitions

abstract: A strong coloring on a cardinal $\kappa$ is a function $f:[\kappa]^2\to \kappa$ such that for every $A\subseteq \kappa$ of full size $\kappa$, every color $\gamma<\kappa$ is attained by $f\upharpoonright[A]^2$. The symbol $\kappa\nrightarrow [\kappa]^2_\kappa$ asserts the existence of a strong coloring on $\kappa$. We introduce the symbol $\kappa\nrightarrow_p[\kappa]^2_\kappa$ which asserts the existence of a coloring $f:[\kappa]^2\to \kappa$ which is strong over a partition $p:[\kappa]^2\to\theta$. A coloring $f$ is strong over $p$ if for every $A\in [\kappa]^\kappa$ there is $i<\theta$ so that every color $\gamma<\kappa$ is attained by $f\upharpoonright ([A]^2\cap p^{-1}(i))$. We prove that whenever $\kappa\nrightarrow[\kappa]^2_\kappa$ holds, also $\kappa\nrightarrow_p[\kappa]^2_\kappa$ holds for an arbitrary finite partition $p$. Similarly, arbitrary finite $p$-s can be added to stronger symbols which hold in any model of ZFC. If $\kappa^\theta=\kappa$, then $\kappa\nrightarrow_p[\kappa]^2_\kappa$ and stronger symbols, like $\mathrm{Pr}_1(\kappa,\kappa,\kappa,\chi)$ or $\mathrm{Pr}_0(\kappa,\kappa,\kappa,\aleph_0)$, hold also for an arbitrary partition $p$ to $\theta$ parts.


title: On certain extensions of vector bundles in p-adic geometry

abstract: Given two arbitrary vector bundles on the Fargues-Fontaine curve, we give an explicit criterion in terms of Harder-Narasimhan polygons on whether they realize a semistable vector bundle as their extensions. Our argument is largely combinatorial and builds upon the dimension analysis of certain moduli spaces of bundle maps developed in [BFH+22].


title: Stochastic Normalizing Flows

abstract: The sampling of probability distributions specified up to a normalization constant is an important problem in both machine learning and statistical mechanics. While classical stochastic sampling methods such as Markov Chain Monte Carlo (MCMC) or Langevin Dynamics (LD) can suffer from slow mixing times there is a growing interest in using normalizing flows in order to learn the transformation of a simple prior distribution to the given target distribution. Here we propose a generalized and combined approach to sample target densities: Stochastic Normalizing Flows (SNF) -- an arbitrary sequence of deterministic invertible functions and stochastic sampling blocks. We show that stochasticity overcomes expressivity limitations of normalizing flows resulting from the invertibility constraint, whereas trainable transformations between sampling steps improve efficiency of pure MCMC/LD along the flow. By invoking ideas from non-equilibrium statistical mechanics we derive an efficient training procedure by which both the sampler's and the flow's parameters can be optimized end-to-end, and by which we can compute exact importance weights without having to marginalize out the randomness of the stochastic blocks. We illustrate the representational power, sampling efficiency and asymptotic correctness of SNFs on several benchmarks including applications to sampling molecular systems in equilibrium.

categories: stat.ML cs.LG physics.chem-ph physics.data-an


title: Combining Observational and Experimental Datasets Using Shrinkage Estimators

abstract: We consider the problem of combining data from observational and experimental sources to make causal conclusions. This problem is increasingly relevant, as the modern era has yielded passive collection of massive observational datasets in areas such as e-commerce and electronic health. These data may be used to supplement experimental data, which is frequently expensive to obtain. In Rosenman et al. (2018), we considered this problem under the assumption that all confounders were measured. Here, we relax the assumption of unconfoundedness. To derive combined estimators with desirable properties, we make use of results from the Stein Shrinkage literature. Our contributions are threefold. First, we propose a generic procedure for deriving shrinkage estimators in this setting, making use of a generalized unbiased risk estimate. Second, we develop two new estimators, prove finite sample conditions under which they have lower risk than an estimator using only experimental data, and show that each achieves a notion of asymptotic optimality. Third, we draw connections between our approach and results in sensitivity analysis, including proposing a method for evaluating the feasibility of our estimators.


title: Relativity of Depth and Sophistication

abstract: Logical depth and sophistication are two quantitative measures of the non-trivial organization of an object. Although apparently different, these measures have been proven equivalent, when the logical depth is renormalized by the busy beaver function. In this article, the measures are relativized to auxiliary information and re-compared to one another. The ability of auxiliary information to solve the halting problem introduces a distortion between the measures. Finally, similar to algorithmic complexity, sophistication and logical depth (renormalized) each offer a relation between their expression of $(x, y)$, $(x)$ and $(y | x)$.


title: Understanding the spatial burden of gender-based violence: Modelling patterns of violence in Nairobi, Kenya through geospatial information

abstract: We present statistical techniques for analyzing global positioning system (GPS) data in order to understand, communicate about, and prevent patterns of violence. In this pilot study, participants in Nairobi, Kenya were asked to rate their safety at several locations, with the goal of predicting safety and learning important patterns. These approaches are meant to help articulate differences in experiences, fostering a discussion that will help communities identify issues and policymakers develop safer communities. A generalized linear mixed model incorporating spatial information taken from existing maps of Kibera showed significant predictors of perceived lack of safety included being alone and time of day; in debrief interviews, participants described feeling unsafe in spaces with hiding places, disease carrying animals, and dangerous individuals. This pilot study demonstrates promise for detecting spatial patterns of violence, which appear to be confirmed by actual rates of measured violence at schools. Several factors relevant to community building consistently predict perceived safety and emerge in participants' qualitative descriptions, telling a cohesive story about perceived safety and empowering communication to community stakeholders.


title: A method of enhancing the detection sensitivity of transient sources in time series with Gaussian stationary noise

abstract: The Gaussian phase noise of intensity time series is demonstrated to be drastically reduced when the raw voltage data are digitally filtered through an arbitrarily large number $n$ of orthornormal bandpass profiles (eigen-filters) sharing the same intensity bandwidth, and the resulting intensity series are co-added. Specifically, the relative noise variance of the summed series at the resolution of one coherence time or less, goes down with increasing $n$ as $1/n$, although (consistent with the radiometer equation) the advantage gradually disappears when the series is bin averaged to lower resolution. Thus the algorithm is designed to enhance the sensitivity of detecting transients that are smoothed out by time averaging and too faint to be visible in the noisy unaveraged time series, as demonstrated by the simulation of a weak embedded time varying signal of either a periodic nature or a fast and unrepeated pulse. The algorithm is then applied to a 10 minute observation of the pulsar PSR 1937+21 by the VLA, where the theoretical predictions were verified by the data. Moreover, it is shown that microstructures within the time profile are better defined as the number $n$ of filters used increases, and a periodic signal of period $1.86 \times 10^{-5}$~s ($53.9$~kHz) is discovered in the pulse profile. Lastly, we apply the algorithm to the first binary black hole merger detected by LIGO, GW150914. We find the SNR of the mean peak intensity increases as $\sqrt{n}$ and cross correlation of the event between the LIGO-Hanford-Livingston detector pair increases with filter order $n$.


title: Computing Boundary Cycle of a Pseudo-Triangle Polygon from its Visibility Graph

abstract: Visibility graph of a simple polygon is a graph with the same vertex set in which there is an edge between a pair of vertices if and only if the segment through them lies completely inside the polygon. Each pair of adjacent vertices on the boundary of the polygon are assumed to be visible. Therefore, the visibility graph of each polygon always contains its boundary edges. This implies that we have always a Hamiltonian cycle in a visibility graph which determines the order of vertices on the boundary of the corresponding polygon. In this paper, we propose a polynomial time algorithm for determining such a Hamiltonian cycle for a pseudo-triangle polygon from its visibility graph.


title: AMOUN: Asymmetric lightweight cryptographic scheme for wireless group communication

abstract: Multi-recipient cryptographic schemes provide secure communication, between one sender and multiple recipients, in a multi-party group. Providing secure multi-party communication is very challenging, especially in dynamic networks. Existing multi-recipient cryptographic schemes pose a variety of limitations. These include high computational overhead for both encryption and decryption, additional communication overhead and high setup cost due to change in membership, and collusion among recipients. In order to overcome these limitations, this paper introduces a novel asymmetric multi-recipient cryptographic scheme, AMOUN. In the proposed scheme, to better utilize network resources, the sender transmits a ciphertext containing different messages to multiple recipients, where each recipient only allowed to retrieve its own designated message. Security analysis demonstrates that the proposed scheme is indistinguishable under adaptive chosen plaintext attack. Quantitative analysis reveals that lightweight AMOUN shows lower average computational cost than both RSA and Multi-RSA, for both encryption and decryption, even when the key sizes are four times larger. For a given prime size, in case of encryption, AMOUN shows 98% and 99% lower average computational cost than RSA and Multi-RSA, respectively. For decryption, AMOUN shows a performance improvement of 99% compared to RSA and Multi-RSA.


title: Multi-layer Representation Fusion for Neural Machine Translation

abstract: Neural machine translation systems require a number of stacked layers for deep models. But the prediction depends on the sentence representation of the top-most layer with no access to low-level representations. This makes it more difficult to train the model and poses a risk of information loss to prediction. In this paper, we propose a multi-layer representation fusion (MLRF) approach to fusing stacked layers. In particular, we design three fusion functions to learn a better representation from the stack. Experimental results show that our approach yields improvements of 0.92 and 0.56 BLEU points over the strong Transformer baseline on IWSLT German-English and NIST Chinese-English MT tasks respectively. The result is new state-of-the-art in German-English translation.


title: BatchEnsemble: An Alternative Approach to Efficient Ensemble and Lifelong Learning

abstract: Ensembles, where multiple neural networks are trained individually and their predictions are averaged, have been shown to be widely successful for improving both the accuracy and predictive uncertainty of single neural networks. However, an ensemble's cost for both training and testing increases linearly with the number of networks, which quickly becomes untenable. In this paper, we propose BatchEnsemble, an ensemble method whose computational and memory costs are significantly lower than typical ensembles. BatchEnsemble achieves this by defining each weight matrix to be the Hadamard product of a shared weight among all ensemble members and a rank-one matrix per member. Unlike ensembles, BatchEnsemble is not only parallelizable across devices, where one device trains one member, but also parallelizable within a device, where multiple ensemble members are updated simultaneously for a given mini-batch. Across CIFAR-10, CIFAR-100, WMT14 EN-DE/EN-FR translation, and out-of-distribution tasks, BatchEnsemble yields competitive accuracy and uncertainties as typical ensembles; the speedup at test time is 3X and memory reduction is 3X at an ensemble of size 4. We also apply BatchEnsemble to lifelong learning, where on Split-CIFAR-100, BatchEnsemble yields comparable performance to progressive neural networks while having a much lower computational and memory costs. We further show that BatchEnsemble can easily scale up to lifelong learning on Split-ImageNet which involves 100 sequential learning tasks.


title: Predicting trends in the quality of state-of-the-art neural networks without access to training or testing data

abstract: In many applications, one works with neural network models trained by someone else. For such pretrained models, one may not have access to training data or test data. Moreover, one may not know details about the model, e.g., the specifics of the training data, the loss function, the hyperparameter values, etc. Given one or many pretrained models, it is a challenge to say anything about the expected performance or quality of the models. Here, we address this challenge by providing a detailed meta-analysis of hundreds of publicly-available pretrained models. We examine norm based capacity control metrics as well as power law based metrics from the recently-developed Theory of Heavy-Tailed Self Regularization. We find that norm based metrics correlate well with reported test accuracies for well-trained models, but that they often cannot distinguish well-trained versus poorly-trained models. We also find that power law based metrics can do much better -- quantitatively better at discriminating among series of well-trained models with a given architecture; and qualitatively better at discriminating well-trained versus poorly-trained models. These methods can be used to identify when a pretrained neural network has problems that cannot be detected simply by examining training/test accuracies.


title: Static adhesion hysteresis in elastic structures

abstract: Adhesive interactions between elastic structures such as graphene sheets, carbon nanotubes, and microtubules have been shown to exhibit hysteresis due to irrecoverable energy loss associated with bond breakage, even in static (rate-independent) experiments. To understand this phenomenon, we start with a minimal theory for the peeling of a thin sheet from a substrate, coupling the local event of bond breaking to the nonlocal elastic relaxation of the sheet and show that this can drive static adhesion hysteresis over a bonding/debonding cycle. Using this model we quantify hysteresis in terms of the adhesion and elasticity parameters of the system. This allows us to derive a scaling relation that preserves hysteresis at different levels of granularity while resolving a seeming paradox of lattice trapping in the continuum limit of a discrete fracture process. Finally, to verify our theory, we use new experiments to demonstrate and measure adhesion hysteresis in bundled microtubules.

categories: physics.bio-ph cond-mat.mtrl-sci cond-mat.soft physics.app-ph q-bio.SC


title: Discrete Approximation of Topologically Protected Modes in Magneto-Optical Media

abstract: Topologically protected waves in the linearly polarized Maxwell's equations with gyrotropic, magneto-optic media were studied a decade ago both computationally and experimentally. This paper develops a robust tight-binding model for this system that makes careful use of Wannier function representations. The model provides very good approximations to the underlying band structure. When solved on a semi-infinite strip, it produces exponentially localized edge modes whose corresponding eigenvalues span the frequency band gaps. A set of coupled differential equations are derived which allows one to find how the electromagnetic field propagates unidirectionally, without backscatter from defects. Furthermore, the discrete model predicts topologically protected edge modes with nontrivial Chern number which are consistent with direct simulation.


title: Reliable Trajectories for Dynamic Quadrupeds using Analytical Costs and Learned Initializations

abstract: Dynamic traversal of uneven terrain is a major objective in the field of legged robotics. The most recent model predictive control approaches for these systems can generate robust dynamic motion of short duration; however, planning over a longer time horizon may be necessary when navigating complex terrain. A recently-developed framework, Trajectory Optimization for Walking Robots (TOWR), computes such plans but does not guarantee their reliability on real platforms, under uncertainty and perturbations. We extend TOWR with analytical costs to generate trajectories that a state-of-the-art whole-body tracking controller can successfully execute. To reduce online computation time, we implement a learning-based scheme for initialization of the nonlinear program based on offline experience. The execution of trajectories as long as 16 footsteps and 5.5 s over different terrains by a real quadruped demonstrates the effectiveness of the approach on hardware. This work builds toward an online system which can efficiently and robustly replan dynamic trajectories.


title: Color, Composition, and Thermal Environment of Kuiper Belt Object (486958) Arrokoth

abstract: The outer Solar System object (486958) Arrokoth (provisional designation 2014 MU$_{69}$) has been largely undisturbed since its formation. We study its surface composition using data collected by the New Horizons spacecraft. Methanol ice is present along with organic material, which may have formed through radiation of simple molecules. Water ice was not detected. This composition indicates hydrogenation of carbon monoxide-rich ice and/ or energetic processing of methane condensed on water ice grains in the cold, outer edge of the early Solar System. There are only small regional variations in color and spectra across the surface, suggesting Arrokoth formed from a homogeneous or well-mixed reservoir of solids. Microwave thermal emission from the winter night side is consistent with a mean brightness temperature of 29$\pm$5 K.


title: How many surface modes does one see on the boundary of a Dirac material?

abstract: We present full expressions for the surface part of polarization tensor of a Dirac fermion confined in a half-space in $3+1$ dimensions. We compare this tensor to the polarization tensor of eventual surface mode (which is a $2+1$ dimensional Dirac fermion) and find essential differences in the conductivities in both Hall and normal sectors. Thus, the interaction with electromagnetic field near the boundary differs significantly in the full model and in the effective theory for the surface mode.


title: Two Component Jets of GRB160623A as Shocked Jet cocoon afterglow

abstract: Two components of jets associated with the afterglow of the gamma-ray burst GRB 160623A were observed with multi-frequency observations including long-term monitoring in a sub-millimetre range (230 GHz) using the SMA. The observed light curves with temporal breaks suggests on the basis of the standard forward-shock synchrotron radiation model that the X-ray radiation is narrowly collimated with an opening angle $\theta_{n,j}<\sim6^{\circ}$ whereas the radio radiation originated from wider jets ($\sim27^{\circ}$). The temporal and spectral evolutions of the radio afterglow agree with those expected from a synchrotron radiation modelling with typical physical parameters except for the fact that the observed wide jet opening angle for the radio emission is significantly larger than the theoretical maximum opening angle. By contrast, the opening angle of the X-ray afterglow is consistent with the typical value of GRB jets. Since the theory of the relativistic cocoon afterglow emission is similar to that of a regular afterglow with an opening angle of $\sim30^{\circ}$, the observed radio emission can be interpreted as the shocked jet cocoon emission. This result therefore indicates that the two components of the jets observed in the GRB 160623A afterglow is caused by the jet and the shocked jet cocoon afterglows.


title: Reward Design for Driver Repositioning Using Multi-Agent Reinforcement Learning

abstract: A large portion of passenger requests is reportedly unserviced, partially due to vacant for-hire drivers' cruising behavior during the passenger seeking process. This paper aims to model the multi-driver repositioning task through a mean field multi-agent reinforcement learning (MARL) approach that captures competition among multiple agents. Because the direct application of MARL to the multi-driver system under a given reward mechanism will likely yield a suboptimal equilibrium due to the selfishness of drivers, this study proposes a reward design scheme with which a more desired equilibrium can be reached. To effectively solve the bilevel optimization problem with upper level as the reward design and the lower level as a multi-agent system, a Bayesian optimization (BO) algorithm is adopted to speed up the learning process. We then apply the bilevel optimization model to two case studies, namely, e-hailing driver repositioning under service charge and multiclass taxi driver repositioning under NYC congestion pricing. In the first case study, the model is validated by the agreement between the derived optimal control from BO and that from an analytical solution. With a simple piecewise linear service charge, the objective of the e-hailing platform can be increased by 8.4%. In the second case study, an optimal toll charge of $5.1 is solved using BO, which improves the objective of city planners by 7.9%, compared to that without any toll charge. Under this optimal toll charge, the number of taxis in the NYC central business district is decreased, indicating a better traffic condition, without substantially increasing the crowdedness of the subway system.


title: The first $p$-widths of the unit disk

abstract: On compact 2-manifolds with non-empty convex boundary, we prove a regularity result for integral 1-varifolds $V$ that are stationary with free boundary and $\mathbb{Z}_2$-almost minimizing in small annuli. That regularity says that $V$ is a free boundary finite geodesic network. Next, using that regularity, we compute the first $p$-widths of the unit closed ball $B^2,$ for $p=1, . . . , 4.$


title: Jet quenching parameters in strongly coupled anisotropic plasmas in the presence of magnetic fields

abstract: I use the holographic gauge/gravity duality to systematically calculate the jet quenching parameters in strongly coupled anisotropic plasmas in the presence of external magnetic fields. The magnetic field breaks down spatial rotation symmetry from $SO(3)$ to $SO(2)$, leading to the presence of multiple anisotropic jet quenching parameters, which are evaluated here in two quite different holographic settings. One of them corresponds to a top-down deformation of the strongly coupled $\mathcal{N} = 4$ Super Yang-Mills plasma triggered by an external magnetic field, while the other one is a bottom-up Einstein-Maxwell-Dilaton model of phenomenological relevance for high energy peripheral heavy ion collisions, since it is able to provide a quantitative description of $(2+1)$-flavors lattice QCD thermodynamics with physical quark masses at zero and nonzero magnetic fields. I find for both models an overall enhancement of all the anisotropic jet quenching parameters with increasing magnetic fields. Moreover, I also conclude that for both models transverse momentum broadening is larger in transverse directions than in the direction of the magnetic field. Since these conclusions are shown to hold for two rather different holographic setups at finite temperature and magnetic fields, they are suggested as fairly robust features of strongly coupled anisotropic magnetized plasmas.


title: On the Approximability of Weighted Model Integration on DNF Structures

abstract: Weighted model counting (WMC) consists of computing the weighted sum of all satisfying assignments of a propositional formula. WMC is well-known to be #P-hard for exact solving, but admits a fully polynomial randomized approximation scheme (FPRAS) when restricted to DNF structures. In this work, we study weighted model integration, a generalization of weighted model counting which involves real variables in addition to propositional variables, and pose the following question: Does weighted model integration on DNF structures admit an FPRAS? Building on classical results from approximate volume computation and approximate weighted model counting, we show that weighted model integration on DNF structures can indeed be approximated for a class of weight functions. Our approximation algorithm is based on three subroutines, each of which can be a weak (i.e., approximate), or a strong (i.e., exact) oracle, and in all cases, comes along with accuracy guarantees. We experimentally verify our approach over randomly generated DNF instances of varying sizes, and show that our algorithm scales to large problem instances, involving up to 1K variables, which are currently out of reach for existing, general-purpose weighted model integration solvers.


title: Generating clause sequences of a CNF formula

abstract: Given a CNF formula $\Phi$ with clauses $C_1,\ldots,C_m$ and variables $V=\{x_1,\ldots,x_n\}$, a truth assignment $a:V\rightarrow\{0,1\}$ of $\Phi$ leads to a clause sequence $\sigma_\Phi(a)=(C_1(a),\ldots,C_m(a))\in\{0,1\}^m$ where $C_i(a) = 1$ if clause $C_i$ evaluates to $1$ under assignment $a$, otherwise $C_i(a) = 0$. The set of all possible clause sequences carries a lot of information on the formula, e.g. SAT, MAX-SAT and MIN-SAT can be encoded in terms of finding a clause sequence with extremal properties. We consider a problem posed at Dagstuhl Seminar 19211 "Enumeration in Data Management" (2019) about the generation of all possible clause sequences of a given CNF with bounded dimension. We prove that the problem can be solved in incremental polynomial time. We further give an algorithm with polynomial delay for the class of tractable CNF formulas. We also consider the generation of maximal and minimal clause sequences, and show that generating maximal clause sequences is NP-hard, while minimal clause sequences can be generated with polynomial delay.


title: Spontaneous emergence of Josephson junctions in homogeneous rings of single-crystal Sr2RuO4

abstract: The chiral p-wave order parameter in Sr$_\mathsf{2}$RuO$_\mathsf{4}$ would make it a special case amongst the unconventional superconductors. A consequence of this symmetry is the possible existence of superconducting domains of opposite chirality. At the boundary of such domains, the locally supressed condensate can produce an intrinsic Josephson junction. Here, we provide evidence of such junctions using mesoscopic rings, structured from Sr$_\mathsf{2}$RuO$_\mathsf{4}$ single crystals. Our order parameter simulations predict such rings to host stable domain walls across their arms. Thisis verified with transport experiments on loops, with a sharp transition at 1.5 K, which show distinct critical current oscillations with periodicity corresponding to the flux quantum. In contrast, loops with broadened transitions at around 3 K are void of such junctions and show standard Little-Parks oscillations. Our analysis demonstrates the junctions are of intrinsic origin and makes a compelling case for the existence of superconducting domains.


title: Vanishing of the Euler class in Power subgroups of the punctured mapping class group

abstract: In this paper, we use a result of Dahmani to show that the Euler class of some power subgroup (the subgroup normally generated by a fixed power of Dehn twist about a non-separating curve) is trivial inside the mapping class group of once punctured surface.


title: SED modelling of broadband emission in the pulsar wind nebula 3C 58

abstract: We investigate broadband emission properties of the pulsar wind nebula (PWN) 3C 58 using a spectral energy distribution (SED) model. We attempt to match simultaneously the broadband SED and spatial variations of X-ray emission in the PWN. We further the model to explain a possible far-IR feature of which a hint is recently suggested in 3C 58: a small bump at $\sim$$10^{11}$ GHz in the PLANCK and Herschel band. While external dust emission may easily explain the observed bump, it may be internal emission of the source implying an additional population of particles. Although significance for the bump is not high, here we explore possible origins of the IR bump using the emission model and find that a population of electrons with GeV energies can explain the bump. If it is produced in the PWN, it may provide new insights into particle acceleration and flows in PWNe.


title: A two-stage algorithm for aircraft conflict resolution with trajectory recovery

abstract: As air traffic volume is continuously increasing, it has become a priority to improve traffic control algorithms to handle future air travel demand and improve airspace capacity. We address the conflict resolution problem in air traffic control using a novel approach for aircraft collision avoidance with trajectory recovery. We present a two-stage algorithm that first solves all initial conflicts by adjusting aircraft headings and speeds, before identifying the optimal time for aircraft to recover towards their target destination. The collision avoidance stage extends an existing mixed-integer programming formulation to heading control. For the trajectory recovery stage, we introduce a novel exact mixed-integer programming formulation as well as a greedy heuristic algorithm. The proposed two-stage approach guarantees that all trajectories during both the collision avoidance and recovery stages are conflict-free. Numerical results on benchmark problems show that the proposed heuristic for trajectory recovery is competitive while also emphasizing the difficulty of this optimization problem. The proposed approach can be used as a decision-support tool for introducing automation in air traffic control.


title: From Fundamental First-Principle Calculations to NanoEngineering Applications: Review of the NESSIE project

abstract: This paper outlines how modern first-principle calculations can adequately address the needs for ever higher levels of numerical accuracy and high-performance in large-scale electronic structure simulations, and pioneer the fundamental study of quantum many-body effects in a large number of emerging nanomaterials.


title: Covert Capacity of Bosonic Channels

abstract: We investigate the quantum-secure covert-communication capabilities of lossy thermal-noise bosonic channels, the quantum-mechanical model for many practical channels. We determine the expressions for the covert capacity of these channels: $L_{\text{no-EA}}$, when Alice and Bob share only a classical secret, and $L_{\text{EA}}$, when they benefit from entanglement assistance. We find that entanglement assistance alters the fundamental scaling law for covert communication. Instead of $L_{\text{no-EA}}\sqrt{n}-r_{\text{no-EA}}(n)$, $r_{\text{no-EA}}(n)=o(\sqrt{n})$, entanglement assistance allows $L_{\text{EA}}\sqrt{n}\log n-r_{\text{EA}}(n)$, $r_{\text{EA}}(n)=o(\sqrt{n}\log n)$, covert bits to be transmitted reliably over $n$ channel uses.


title: Automatic Frame Selection using CNN in Ultrasound Elastography

abstract: Ultrasound elastography is used to estimate the mechanical properties of the tissue by monitoring its response to an internal or external force. Different levels of deformation are obtained from different tissue types depending on their mechanical properties, where stiffer tissues deform less. Given two radio frequency (RF) frames collected before and after some deformation, we estimate displacement and strain images by comparing the RF frames. The quality of the strain image is dependent on the type of motion that occurs during deformation. In-plane axial motion results in high-quality strain images, whereas out-of-plane motion results in low-quality strain images. In this paper, we introduce a new method using a convolutional neural network (CNN) to determine the suitability of a pair of RF frames for elastography in only 5.4 ms. Our method could also be used to automatically choose the best pair of RF frames, yielding a high-quality strain image. The CNN was trained on 3,818 pairs of RF frames, while testing was done on 986 new unseen pairs, achieving an accuracy of more than 91%. The RF frames were collected from both phantom and in vivo data.


title: SpotTheFake: An Initial Report on a New CNN-Enhanced Platform for Counterfeit Goods Detection

abstract: The counterfeit goods trade represents nowadays more than 3.3% of the whole world trade and thus it's a problem that needs now more than ever a lot of attention and a reliable solution that would reduce the negative impact it has over the modern society. This paper presents the design and early stage development of a novel counterfeit goods detection platform that makes use of the outstsanding learning capabilities of the classical VGG16 convolutional model trained through the process of "transfer learning" and a multi-stage fake detection procedure that proved to be not only reliable but also very robust in the experiments we have conducted so far using an image dataset of various goods which we gathered ourselves.


title: Directional Deep Embedding and Appearance Learning for Fast Video Object Segmentation

abstract: Most recent semi-supervised video object segmentation (VOS) methods rely on fine-tuning deep convolutional neural networks online using the given mask of the first frame or predicted masks of subsequent frames. However, the online fine-tuning process is usually time-consuming, limiting the practical use of such methods. We propose a directional deep embedding and appearance learning (DDEAL) method, which is free of the online fine-tuning process, for fast VOS. First, a global directional matching module, which can be efficiently implemented by parallel convolutional operations, is proposed to learn a semantic pixel-wise embedding as an internal guidance. Second, an effective directional appearance model based statistics is proposed to represent the target and background on a spherical embedding space for VOS. Equipped with the global directional matching module and the directional appearance model learning module, DDEAL learns static cues from the labeled first frame and dynamically updates cues of the subsequent frames for object segmentation. Our method exhibits state-of-the-art VOS performance without using online fine-tuning. Specifically, it achieves a J & F mean score of 74.8% on DAVIS 2017 dataset and an overall score G of 71.3% on the large-scale YouTube-VOS dataset, while retaining a speed of 25 fps with a single NVIDIA TITAN Xp GPU. Furthermore, our faster version runs 31 fps with only a little accuracy loss. Our code and trained networks are available at https://github.com/YingjieYin/Directional-Deep-Embedding-and-Appearance-Learning-for-Fast-Video-Object-Segmentation.


title: Four-valued monitorability of $\omega$-regular languages

abstract: Runtime Verification (RV) is a lightweight formal technique in which program or system execution is monitored and analyzed, to check whether certain properties are satisfied or violated after a finite number of steps. The use of RV has led to interest in deciding whether a property is monitorable: whether it is always possible for the satisfaction or violation of the property to be determined after a finite future continuation. However, classical two-valued monitorability suffers from two inherent limitations. First, a property can only be evaluated as monitorable or non-monitorable; no information is available regarding whether only one verdict (satisfaction or violation) can be detected. Second, monitorability is defined at the language-level and does not tell us whether satisfaction or violation can be detected starting from the current monitor state during system execution. To address these limitations, this paper proposes a new notion of four-valued monitorability for $\omega$-languages and applies it at the state-level. Four-valued monitorability is more informative than two-valued monitorability as a property can be evaluated as a four-valued result, denoting that only satisfaction, only violation, or both are active for a monitorable property. We can also compute state-level weak monitorability, i.e., whether satisfaction or violation can be detected starting from a given state in a monitor, which enables state-level optimizations of monitoring algorithms. Based on a new six-valued semantics, we propose procedures for computing four-valued monitorability of $\omega$-regular languages, both at the language-level and at the state-level. We have developed a new tool that implements the proposed procedure for computing monitorability of LTL formulas.


title: Shifted Lanczos method for quadratic forms with Hermitian matrix resolvents

abstract: Quadratic forms of Hermitian matrix resolvents involve the solutions of shifted linear systems. Efficient iterative solutions use the shift-invariance property of Krylov subspaces The Hermitian Lanczos method reduces a given vector and matrix to a Jacobi matrix (real symmetric tridiagonal matrix with positive super and sub-diagonal entries) and approximates the quadratic form using the Jacobi matrix. This study develops a shifted Lanczos method that deals directly with the Hermitian matrix resolvent. We derive a matrix representation of a linear operator that approximates the resolvent by solving a Vorobyev moment problem associated with the shifted Lanczos method. We show that an entry of the Jacobi matrix resolvent can approximate the quadratic form, matching the moments. We give a sufficient condition such that the method does not break down, an error bound, and error estimates. Numerical experiments on matrices drawn from real-world applications compare the proposed method with previous methods and show that the proposed method outperforms well-established methods in solving some problems.


title: Multiple Flat Projections for Cross-manifold Clustering

abstract: Cross-manifold clustering is a hard topic and many traditional clustering methods fail because of the cross-manifold structures. In this paper, we propose a Multiple Flat Projections Clustering (MFPC) to deal with cross-manifold clustering problems. In our MFPC, the given samples are projected into multiple subspaces to discover the global structures of the implicit manifolds. Thus, the cross-manifold clusters are distinguished from the various projections. Further, our MFPC is extended to nonlinear manifold clustering via kernel tricks to deal with more complex cross-manifold clustering. A series of non-convex matrix optimization problems in MFPC are solved by a proposed recursive algorithm. The synthetic tests show that our MFPC works on the cross-manifold structures well. Moreover, experimental results on the benchmark datasets show the excellent performance of our MFPC compared with some state-of-the-art clustering methods.


title: The spectral dimension of human mobility

abstract: Human mobility patterns are surprisingly structured. In spite of many hard to model factors, such as climate, culture, and socioeconomic opportunities, aggregate migration rates obey a universal, parameter-free, `radiation' model. Recent work has further shown that the detailed spectral decomposition of these flows -- defined as the number of individuals that visit a given location with frequency $f$ from a distance $r$ away -- also obeys simple rules, namely, scaling as a universal inverse square law in the combination, $rf$. However, this surprising regularity, derived on general grounds, has not been explained through microscopic mechanisms of individual behavior. Here we confirm this by analyzing large-scale cell phone datasets from three distinct regions and show that a direct consequence of this scaling law is that the average `travel energy' spent by visitors to a given location is constant across space, a finding reminiscent of the well-known travel budget hypothesis of human movement. The attractivity of different locations, which we define by the total number of visits to that location, also admits non-trivial, spatially-clustered structure. The observed pattern is consistent with the well-known central place theory in urban geography, as well as with the notion of Weber optimality in spatial economy, hinting to a collective human capacity of optimizing recurrent movements. We close by proposing a simple, microscopic human mobility model which simultaneously captures all our empirical findings. Our results have relevance for transportation, urban planning, geography, and other disciplines in which a deeper understanding of aggregate human mobility is key.


title: Representations with minimal support for quantized Gieseker varieties

abstract: We study the minimally supported representations of quantizations of Gieseker moduli spaces. We relate them to $\operatorname{SL}_n$-equivariant D-modules on the nilpotent cone of $\mathfrak{sl}_n$ and to minimally supported representations of type A rational Cherednik algebras. Our main result is character formulas for minimally supported representations of quantized Gieseker moduli spaces.


title: Individual Fairness for $k$-Clustering

abstract: We give a local search based algorithm for $k$-median and $k$-means (and more generally for any $k$-clustering with $\ell_p$ norm cost function) from the perspective of individual fairness. More precisely, for a point $x$ in a point set $P$ of size $n$, let $r(x)$ be the minimum radius such that the ball of radius $r(x)$ centered at $x$ has at least $n/k$ points from $P$. Intuitively, if a set of $k$ random points are chosen from $P$ as centers, every point $x\in P$ expects to have a center within radius $r(x)$. An individually fair clustering provides such a guarantee for every point $x\in P$. This notion of fairness was introduced in [Jung et al., 2019] where they showed how to get an approximately feasible $k$-clustering with respect to this fairness condition. In this work, we show how to get a bicriteria approximation for fair $k$-clustering: The $k$-median ($k$-means) cost of our solution is within a constant factor of the cost of an optimal fair $k$-clustering, and our solution approximately satisfies the fairness condition (also within a constant factor). Further, we complement our theoretical bounds with empirical evaluation.


title: Study of polycrystalline bulk Sr$_3$OsO$_6$ double-perovskite insulator: comparison with 1000 K ferromagnetic epitaxial films

abstract: Polycrystalline Sr$_3$OsO$_6$, which is an ordered double-perovskite insulator, is synthesized via solid-state reaction under high-temperature and high-pressure conditions of 1200 $^\circ$C and 6 GPa. The synthesis enables us to conduct a comparative study of the bulk form of Sr$_3$OsO$_6$ toward revealing the driving mechanism of 1000 K ferromagnetism, which has recently been discovered for epitaxially grown Sr$_3$OsO$_6$ films. Unlike the film, the bulk is dominated by antiferromagnetism rather than ferromagnetism. Therefore, robust ferromagnetic order appears only when Sr$_3$OsO$_6$ is under the influence of interfaces. A specific heat capacity of 39.6(9) 10$^{-3}$ J mol$^{-1}$ K$^{-2}$ is found at low temperatures ($<$17 K). This value is remarkably high, suggesting the presence of possible fermionic-like excitations at the magnetic ground state. Although the bulk and film forms of Sr$_3$OsO$_6$ share the same lattice basis and electrically insulating state, the magnetism is entirely different between them.


title: Criterion for fingering instabilities in colloidal gels

abstract: We sandwich a colloidal gel between two parallel plates and induce a radial flow by lifting the upper plate at a constant velocity. Two distinct scenarios result from such a tensile test: ($i$) stable flows during which the gel undergoes a tensile deformation without yielding, and ($ii$) unstable flows characterized by the radial growth of air fingers into the gel. We show that the unstable regime occurs beyond a critical energy input, independent of the gel's macroscopic yield stress. This implies a local fluidization of the gel at the tip of the growing fingers and results in the most unstable wavelength of the patterns exhibiting the characteristic scalings of the classical viscous fingering instability. Our work provides a quantitative criterion for the onset of fingering in colloidal gels based on a local shear-induced yielding, in agreement with the delayed failure framework.


title: Low-PMEPR Preamble Sequence Design for Dynamic Spectrum Allocation in OFDMA Systems

abstract: Orthogonal Frequency Division Multiple Access (OFDMA) with Dynamic spectrum allocation (DSA) is able to provide a wide range of data rate requirements. This paper is focused on the design of preamble sequences in OFDMA systems with low peak-to-mean envelope power ratio (PMEPR) property in the context of DSA. We propose a systematic preamble sequence design which gives rise to low PMEPR for possibly non-contiguous spectrum allocations. With the aid of Golay-Davis-Jedwab (GDJ) sequences, two classes of preamble sequences are presented. We prove that their PMEPRs are upper bounded by 4 for any DSA over a chunk of four contiguous resource blocks.


title: Learning Individually Fair Classifier with Path-Specific Causal-Effect Constraint

abstract: Machine learning is used to make decisions for individuals in various fields, which require us to achieve good prediction accuracy while ensuring fairness with respect to sensitive features (e.g., race and gender). This problem, however, remains difficult in complex real-world scenarios. To quantify unfairness under such situations, existing methods utilize {\it path-specific causal effects}. However, none of them can ensure fairness for each individual without making impractical functional assumptions on the data. In this paper, we propose a far more practical framework for learning an individually fair classifier. To avoid restrictive functional assumptions, we define the {\it probability of individual unfairness} (PIU) and solve an optimization problem where PIU's upper bound, which can be estimated from data, is controlled to be close to zero. We elucidate why our method can guarantee fairness for each individual. Experimental results show that our method can learn an individually fair classifier at a slight cost of accuracy.


title: The nonlinear fractional diffusion equations with Nagumo-type sources and perturbed orders

abstract: We consider a class of nonlinear fractional equations having the Caputo fractional derivative of the time variable $t$, the fractional order of the self-adjoint positive definite unbounded operator in a Hilbert space and a singular nonlinear source. These equations are generalizations of some well-known fractional equation such as the fractional Cahn-Allen equation, the fractional Burger equation, the fractional Cahn-Hilliard equation, the fractional Kuramoto-Sivashinsky equation, etc. We study both the initial value and the final value problem. Under some suitable assumptions, we investigate the existence, uniqueness of maximal solution, and stability of solution of the problems with respect to perturbed fractional orders. For $t=0$, we show that the final value problem is instable and deduce that the problem is ill-posed. A regularization method is proposed to recover the initial data from the inexact fractional orders and the final data. By some regularity assumptions of the exact solutions of the problems, we obtain an error estimate of H\"older type.


title: From Matching with Diversity Constraints to Matching with Regional Quotas

abstract: In the past few years, several new matching models have been proposed and studied that take into account complex distributional constraints. Relevant lines of work include (1) school choice with diversity constraints where students have (possibly overlapping) types and (2) hospital-doctor matching where various regional quotas are imposed. In this paper, we present a polynomial-time reduction to transform an instance of (1) to an instance of (2) and we show how the feasibility and stability of corresponding matchings are preserved under the reduction. Our reduction provides a formal connection between two important strands of work on matching with distributional constraints. We then apply the reduction in two ways. Firstly, we show that it is NP-complete to check whether a feasible and stable outcome for (1) exists. Due to our reduction, these NP-completeness results carry over to setting (2). In view of this, we help unify some of the results that have been presented in the literature. Secondly, if we have positive results for (2), then we have corresponding results for (1). One key conclusion of our results is that further developments on axiomatic and algorithmic aspects of hospital-doctor matching with regional quotas will result in corresponding results for school choice with diversity constraints.


title: Angular Expansion of Nova Shells

abstract: Nova shells can provide us with important information on their distance, their interactions with the circumstellar and interstellar media, and the evolution in morphology of the ejecta. We have obtained narrow-band images of a sample of five nova shells, namely DQHer, FHSer, TAur, V476Cyg, and V533Her, with ages in the range from 50 to 130 years. These images have been compared with suitable available archival images to derive their angular expansion rates. We find that all the nova shells in our sample are still in the free expansion phase, which can be expected, as the mass of the ejecta is 7-45 times larger than the mass of the swept-up circumstellar medium. The nova shells will keep expanding freely for time periods up to a few hundred years, reducing their time dispersal into the interstellar medium.


title: The alpha invariant of complete intersections

abstract: We compute the alpha invariant of any smooth complex projective spin complete intersection of complex dimension $1 \; ({\rm mod} \; 4)$. We prove that the alpha invariant depends only on the total degree and Pontryagin classes. Our findings are consistent with a long-standing conjecture, often called the Sullivan Conjecture, which states that two complete intersections with the same dimensions, total degrees, Pontryagin and Euler classes are diffeomorphic.


title: Second-order Conic Programming Approach for Wasserstein Distributionally Robust Two-stage Linear Programs

abstract: This paper proposes a second-order conic programming (SOCP) approach to solve distributionally robust two-stage stochastic linear programs over 1-Wasserstein balls. We start from the case with distribution uncertainty only in the objective function and exactly reformulate it as an SOCP problem. Then, we study the case with distribution uncertainty only in constraints, and show that such a robust program is generally NP-hard as it involves a norm maximization problem over a polyhedron. However, it is reduced to an SOCP problem if the extreme points of the polyhedron are given as a prior. This motivates to design a constraint generation algorithm with provable convergence to approximately solve the NP-hard problem. In sharp contrast to the exiting literature, the distribution achieving the worst-case cost is given as an "empirical" distribution by simply perturbing each sample for both cases. Finally, experiments illustrate the advantages of the proposed model in terms of the out-of-sample performance and the computational complexity.


title: Overlaid oriented Voronoi diagrams and the 1-Steiner tree problem

abstract: Overlaid oriented Voronoi diagrams (OOVDs) are known to provide useful data for the construction of optimal Euclidean $1$-Steiner trees. The theoretical time complexity of construction methods exploiting the OOVD is $O(n^2)$, but a computational study has never been performed, and robust constructions for OOVDs have not previously been implemented. In this paper, we outline a numerically stable implementation for constructing OOVDs using tools from the Computational Geometry Algorithms Library (CGAL), and test its performance on random point sets. We then study the effect that the OOVD data has in reducing the complexity of $1$-Steiner tree construction when compared to a naive approach. The number of iterations of the main loop of the 1-Steiner algorithm is directly determined by the number of faces in the OOVD, and this appears to be linear for the random inputs we tested. We also discuss methods for processing the OOVD data that lead to a reduction in construction time by roughly a factor of 12.


title: Unraveling Meta-Learning: Understanding Feature Representations for Few-Shot Tasks

abstract: Meta-learning algorithms produce feature extractors which achieve state-of-the-art performance on few-shot classification. While the literature is rich with meta-learning methods, little is known about why the resulting feature extractors perform so well. We develop a better understanding of the underlying mechanics of meta-learning and the difference between models trained using meta-learning and models which are trained classically. In doing so, we introduce and verify several hypotheses for why meta-learned models perform better. Furthermore, we develop a regularizer which boosts the performance of standard training routines for few-shot classification. In many cases, our routine outperforms meta-learning while simultaneously running an order of magnitude faster.


title: Velocity Gradient in the Presence of Self-Gravity: Identifying Gravity-induced Inflow and Determining Collapsing Stage

abstract: Understanding how star formation is regulated requires studying the energy balance between turbulence, magnetic fields, stellar feedback, and gravity within molecular clouds. However, identifying the transition region where the gravity takes over remains elusive. Recent studies of the Velocity Gradient Technique (VGT), which is an advanced tool for magnetic field studies, reveal that the gradients of spectroscopic observables change their directions by 90 degrees with respect to the magnetic fields in the regions of gravitational collapse. In this study, we perform 3D MHD numerical simulations. We observe that star formation successfully proceeds in strongly magnetized and fully ionized media. We confirm that the self-gravity induces the change of gradients' orientation and high gradients' amplitude. We explore two ways of identifying collapsing self-gravitating regions through the double-peak feature in the histogram of gradients' orientation and the curvature of gradients. We show that velocity gradients' morphology and amplitude can be synthetically used to trace the convergent inflows. By comparing with the column density Probability Density Functions (N-PDFs) method, we show that VGT is a powerful new tool for studying the gas dynamics and tracing magnetic field in star-forming regions. By analogy with VGT, we extend the Intensity Gradient Technique (IGT) to locate the gravitational collapsing region and shocks. We demonstrate the synergy of VGT and IGT can determine the collapsing stages in a star-forming region.


title: Unifying Graph Convolutional Neural Networks and Label Propagation

abstract: Label Propagation (LPA) and Graph Convolutional Neural Networks (GCN) are both message passing algorithms on graphs. Both solve the task of node classification but LPA propagates node label information across the edges of the graph, while GCN propagates and transforms node feature information. However, while conceptually similar, theoretical relation between LPA and GCN has not yet been investigated. Here we study the relationship between LPA and GCN in terms of two aspects: (1) feature/label smoothing where we analyze how the feature/label of one node is spread over its neighbors; And, (2) feature/label influence of how much the initial feature/label of one node influences the final feature/label of another node. Based on our theoretical analysis, we propose an end-to-end model that unifies GCN and LPA for node classification. In our unified model, edge weights are learnable, and the LPA serves as regularization to assist the GCN in learning proper edge weights that lead to improved classification performance. Our model can also be seen as learning attention weights based on node labels, which is more task-oriented than existing feature-based attention models. In a number of experiments on real-world graphs, our model shows superiority over state-of-the-art GCN-based methods in terms of node classification accuracy.


title: Strong convergence and asymptotic stability of explicit numerical schemes for nonlinear stochastic differential equations

abstract: In this article we introduce several kinds of easily implementable explicit schemes, which are amenable to Khasminski's techniques and are particularly suitable for highly nonlinear stochastic differential equations (SDEs). We show that without additional restriction conditions except those which guarantee the exact solutions possess their boundedness in expectation with respect to certain Lyapunov functions, the numerical solutions converge strongly to the exact solutions in finite-time. Moreover, based on the nonnegative semimartingale convergence theorem, positive results about the ability of explicit numerical approximation to reproduce the well-known LaSalle-type theorem of SDEs are proved here, from which we deduce the asymptotic stability of numerical solutions. Some examples and simulations are provided to support the theoretical results and to demonstrate the validity of the approach.


title: Relational Message Passing for Knowledge Graph Completion

abstract: Knowledge graph completion aims to predict missing relations between entities in a knowledge graph. In this work, we propose a relational message passing method for knowledge graph completion. Different from existing embedding-based methods, relational message passing only considers edge features (i.e., relation types) without entity IDs in the knowledge graph, and passes relational messages among edges iteratively to aggregate neighborhood information. Specifically, two kinds of neighborhood topology are modeled for a given entity pair under the relational message passing framework: (1) Relational context, which captures the relation types of edges adjacent to the given entity pair; (2) Relational paths, which characterize the relative position between the given two entities in the knowledge graph. The two message passing modules are combined together for relation prediction. Experimental results on knowledge graph benchmarks as well as our newly proposed dataset show that, our method PathCon outperforms state-of-the-art knowledge graph completion methods by a large margin. PathCon is also shown applicable to inductive settings where entities are not seen in training stage, and it is able to provide interpretable explanations for the predicted results. The code and all datasets are available at https://github.com/hwwang55/PathCon.


title: Interactive Text-to-Speech System via Joint Style Analysis

abstract: While modern TTS technologies have made significant advancements in audio quality, there is still a lack of behavior naturalness compared to conversing with people. We propose a style-embedded TTS system that generates styled responses based on the speech query style. To achieve this, the system includes a style extraction model that extracts a style embedding from the speech query, which is then used by the TTS to produce a matching response. We faced two main challenges: 1) only a small portion of the TTS training dataset has style labels, which is needed to train a multi-style TTS that respects different style embeddings during inference. 2) The TTS system and the style extraction model have disjoint training datasets. We need consistent style labels across these two datasets so that the TTS can learn to respect the labels produced by the style extraction model during inference. To solve these, we adopted a semi-supervised approach that uses the style extraction model to create style labels for the TTS dataset and applied transfer learning to learn the style embedding jointly. Our experiment results show user preference for the styled TTS responses and demonstrate the style-embedded TTS system's capability of mimicking the speech query style.


title: Asymptotically flat vacuum initial data sets from a modified parabolic-hyperbolic formulation of the Einstein vacuum constraint equations

abstract: In this paper we continue earlier investigations of evolutionary formulations of the Einstein vacuum constraint equations originally introduced by R\'{a}cz. Motivated by the strong evidence from these works that the resulting vacuum initial data sets are generically not asymptotically flat we analyse the asymptotics of the solutions of a modified formulation by a combination of analytical and numerical techniques. We conclude that the vacuum initial data sets generated with this new formulation are generically asymptotically flat.


title: Precoder Design for mmWave UAV Communications with Physical Layer Security

abstract: The integration of unmanned aerial vehicles (UAVs) into the terrestrial cellular networks is envisioned as one key technology for next-generation wireless communications. In this work, we consider the physical layer security of the communications links in the millimeter-wave (mmWave) spectrum which are maintained by UAVs functioning as base stations (BS). In particular, we propose a new precoding strategy which incorporates the channel state information (CSI) of the eavesdropper (Eve) compromising link security. We show that our proposed precoder strategy eliminates any need for artificial noise (AN) transmission in underloaded scenarios (fewer users than number of antennas). In addition, we demonstrate that our nonlinear precoding scheme provides promising secrecy-rate performance even for overloaded scenarios at the expense of transmitting low-power AN.


title: Hybrid Embedded Deep Stacked Sparse Autoencoder with w_LPPD SVM Ensemble

abstract: Deep learning is a kind of feature learning method with strong nonliear feature transformation and becomes more and more important in many fields of artificial intelligence. Deep autoencoder is one representative method of the deep learning methods, and can effectively extract abstract the information of datasets. However, it does not consider the complementarity between the deep features and original features during deep feature transformation. Besides, it suffers from small sample problem. In order to solve these problems, a novel deep autoencoder - hybrid feature embedded stacked sparse autoencoder(HESSAE) has been proposed in this paper. HFESAE is capable to learn discriminant deep features with the help of embedding original features to filter weak hidden-layer outputs during training. For the issue that class representation ability of abstract information is limited by small sample problem, a feature fusion strategy has been designed aiming to combining abstract information learned by HFESAE with original feature and obtain hybrid features for feature reduction. The strategy is hybrid feature selection strategy based on L1 regularization followed by an support vector machine(SVM) ensemble model, in which weighted local discriminant preservation projection (w_LPPD), is designed and employed on each base classifier. At the end of this paper, several representative public datasets are used to verify the effectiveness of the proposed algorithm. The experimental results demonstrated that, the proposed feature learning method yields superior performance compared to other existing and state of art feature learning algorithms including some representative deep autoencoder methods.


title: How fast can you update your MST? (Dynamic algorithms for cluster computing)

abstract: Imagine a large graph that is being processed by a cluster of computers, e.g., described by the $k$-machine model or the Massively Parallel Computation Model. The graph, however, is not static; instead it is receiving a constant stream of updates. How fast can the cluster process the stream of updates? The fundamental question we want to ask in this paper is whether we can update the graph fast enough to keep up with the stream. We focus specifically on the problem of maintaining a minimum spanning tree (MST), and we give an algorithm for the $k$-machine model that can process $O(k)$ graph updates per $O(1)$ rounds with high probability. (And these results carry over to the Massively Parallel Computation (MPC) model.) We also show a lower bound, i.e., it is impossible to process $k^{1+\epsilon}$ updates in $O(1)$ rounds. Thus we provide a nearly tight answer to the question of how fast a cluster can respond to a stream of graph modifications while maintaining an MST.


title: Optimization and Validation of a Deep Learning CuZr Atomistic Potential: Robust Applications for Crystalline and Amorphous Phases with near-DFT Accuracy

abstract: We show that a deep-learning neural network potential (DP) based on density functional theory (DFT) calculations can well describe Cu-Zr materials, an example of a binary alloy system that can coexist in several ordered intermetallics and as an amorphous phase. The complex phase diagram for Cu-Zr makes it a challenging system for traditional atomistic force-fields that fail to describe well the different properties and phases. Instead, we show that a DP approach using a large database with ~300k configurations can render results generally on par with DFT. The training set includes configurations of pristine and bulk elementary metals and intermetallics in the liquid and solid phases in addition to slab and amorphous configurations. The DP model was validated by comparing bulk properties such as lattice constants, elastic constants, bulk moduli, phonon spectra, surface energies to DFT values for identical structures. Further, we contrast the DP results with values obtained using well-established two embedded atom method potentials. Overall, our DP potential provides near DFT accuracy for the different Cu-Zr phases but with a fraction of its computational cost, thus enabling accurate computations of realistic atomistic models especially for the amorphous phase.


title: Computing Covers under Substring Consistent Equivalence Relations

abstract: Covers are a kind of quasiperiodicity in strings. A string $C$ is a cover of another string $T$ if any position of $T$ is inside some occurrence of $C$ in $T$. The shortest and longest cover arrays of $T$ have the lengths of the shortest and longest covers of each prefix of $T$, respectively. The literature has proposed linear-time algorithms computing longest and shortest cover arrays taking border arrays as input. An equivalence relation $\approx$ over strings is called a substring consistent equivalence relation (SCER) iff $X \approx Y$ implies (1) $|X| = |Y|$ and (2) $X[i:j] \approx Y[i:j]$ for all $1 \le i \le j \le |X|$. In this paper, we generalize the notion of covers for SCERs and prove that existing algorithms to compute the shortest cover array and the longest cover array of a string $T$ under the identity relation will work for any SCERs taking the accordingly generalized border arrays.


title: Superpixel Segmentation via Convolutional Neural Networks with Regularized Information Maximization

abstract: We propose an unsupervised superpixel segmentation method by optimizing a randomly-initialized convolutional neural network (CNN) in inference time. Our method generates superpixels via CNN from a single image without any labels by minimizing a proposed objective function for superpixel segmentation in inference time. There are three advantages to our method compared with many of existing methods: (i) leverages an image prior of CNN for superpixel segmentation, (ii) adaptively changes the number of superpixels according to the given images, and (iii) controls the property of superpixels by adding an auxiliary cost to the objective function. We verify the advantages of our method quantitatively and qualitatively on BSDS500 and SBD datasets.


title: Evolutionary Bi-objective Optimization for the Dynamic Chance-Constrained Knapsack Problem Based on Tail Bound Objectives

abstract: Real-world combinatorial optimization problems are often stochastic and dynamic. Therefore, it is essential to make optimal and reliable decisions with a holistic approach. In this paper, we consider the dynamic chance-constrained knapsack problem where the weight of each item is stochastic, the capacity constraint changes dynamically over time, and the objective is to maximize the total profit subject to the probability that total weight exceeds the capacity. We make use of prominent tail inequalities such as Chebyshev's inequality, and Chernoff bound to approximate the probabilistic constraint. Our key contribution is to introduce an additional objective which estimates the minimal capacity bound for a given stochastic solution that still meets the chance constraint. This objective helps to cater for dynamic changes to the stochastic problem. We apply single- and multi-objective evolutionary algorithms to the problem and show how bi-objective optimization can help to deal with dynamic chance-constrained problems.


title: Building Flat-Band Lattice Models from Gram Matrices

abstract: We propose a powerful and convenient method to systematically design flat-band lattice models, which overcomes the difficulties underlying the previous method. Especially, our method requires no elaborate calculations, applies to arbitrary spatial dimensions, and guarantees to result in a completely flat ground band. We use this method to generate several classes of lattice models, including models with both short- and long-range hoppings, both topologically trivial and non-trivial flat bands. Some of these models were previously known. Our method, however, provides crucial new insights. For example, we have reproduced and generalized the Kapit-Mueller model [Kapit and Mueller, Phys. Rev. Lett. \textbf{105}, 215303 (2010)] and demonstrated a universal scaling rule between the flat band degeneracy and the magnetic flux that was not noticed in previous studies. We show that the flat band of this model results from the (over-)completeness properties of coherent states.


title: Last iterate convergence in no-regret learning: constrained min-max optimization for convex-concave landscapes

abstract: In a recent series of papers it has been established that variants of Gradient Descent/Ascent and Mirror Descent exhibit last iterate convergence in convex-concave zero-sum games. Specifically, \cite{DISZ17, LiangS18} show last iterate convergence of the so called "Optimistic Gradient Descent/Ascent" for the case of \textit{unconstrained} min-max optimization. Moreover, in \cite{Metal} the authors show that Mirror Descent with an extra gradient step displays last iterate convergence for convex-concave problems (both constrained and unconstrained), though their algorithm does not follow the online learning framework; it uses extra information rather than \textit{only} the history to compute the next iteration. In this work, we show that "Optimistic Multiplicative-Weights Update (OMWU)" which follows the no-regret online learning framework, exhibits last iterate convergence locally for convex-concave games, generalizing the results of \cite{DP19} where last iterate convergence of OMWU was shown only for the \textit{bilinear case}. We complement our results with experiments that indicate fast convergence of the method.


title: Schmidt's subspace theorem for moving hypersurface targets in subgeneral position with index in algebraic variety

abstract: Recently, Xie-Cao [15] obtained a Second Main Theorem for moving hypersurfaces located in subgeneral position with index which is extended the result of Ru [11]. By using some methods due to Son-Tan-Thin [13], Quang [9] and Xie-Cao [15], we shall give a Schmidt's Subspace Theorem for moving hypersurface targets in subgeneral position with index intersecting algebraic variety. Our result is a extension the Schmidt's Subspace Theorem due to Son-Tan-Thin [13] and Quang [9].


title: Unsupervised Image-generation Enhanced Adaptation for Object Detection in Thermal images

abstract: Object detection in thermal images is an important computer vision task and has many applications such as unmanned vehicles, robotics, surveillance and night vision. Deep learning based detectors have achieved major progress, which usually need large amount of labelled training data. However, labelled data for object detection in thermal images is scarce and expensive to collect. How to take advantage of the large number labelled visible images and adapt them into thermal image domain, is expected to solve. This paper proposes an unsupervised image-generation enhanced adaptation method for object detection in thermal images. To reduce the gap between visible domain and thermal domain, the proposed method manages to generate simulated fake thermal images that are similar to the target images, and preserves the annotation information of the visible source domain. The image generation includes a CycleGAN based image-to-image translation and an intensity inversion transformation. Generated fake thermal images are used as renewed source domain. And then the off-the-shelf Domain Adaptive Faster RCNN is utilized to reduce the gap between generated intermediate domain and the thermal target domain. Experiments demonstrate the effectiveness and superiority of the proposed method.


title: Toeplitz operators on weighted Bergman spaces induced by a class of radial weights

abstract: Suppose that $\omega$ is a radial weight on the unit disk that satisfies both forward and reverse doubling conditions. Using Carleson measures and $T1$-type conditions, we obtain necessary and sufficient conditions of the positive Borel measure $\mu$ such that the Toeplitz operator $T_{\mu,\omega}:L^p_a(\omega)\to L_a^1(\omega)$ is bounded and compact for $0<p\leq 1$. In addition, we obtain a bump condition for the bounded Toeplitz operators with $L^1(\omega)$ symbol on $L^1_a(\omega)$. This generalizes a result of Zhu in \cite{zhu1989}.


title: Differentiable Bandit Exploration

abstract: Exploration policies in Bayesian bandits maximize the average reward over problem instances drawn from some distribution $\mathcal{P}$. In this work, we learn such policies for an unknown distribution $\mathcal{P}$ using samples from $\mathcal{P}$. Our approach is a form of meta-learning and exploits properties of $\mathcal{P}$ without making strong assumptions about its form. To do this, we parameterize our policies in a differentiable way and optimize them by policy gradients, an approach that is general and easy to implement. We derive effective gradient estimators and introduce novel variance reduction techniques. We also analyze and experiment with various bandit policy classes, including neural networks and a novel softmax policy. The latter has regret guarantees and is a natural starting point for our optimization. Our experiments show the versatility of our approach. We also observe that neural network policies can learn implicit biases expressed only through the sampled instances.


title: How Do Space-Time Digital Metasurfaces Serve to Perform Analog Signal Processing?

abstract: In the quest to realize analog signal processing using sub-wavelength metasurfaces, in this paper, we demonstrate the first experimental demonstration of programmable time-modulated metasurface processors based on the key properties of spatial Fourier transformation. Exploiting space-time coding strategy enables local, independent, and real-time engineering of not only amplitude but also phase profile of the contributing reflective digital meta-atoms at both central and harmonic frequencies. Several illustrative examples are demonstrated to show that the proposed multifunctional calculus metasurface is capable of implementing a large class of useful mathematical operators, including 1st- and 2nd-order spatial differentiation, 1st-order spatial integration, and integro-differential equation solving accompanied by frequency conversions. Unlike the recent proposals, the designed time-modulated signal processor effectively operates for input signals containing wide spatial frequency bandwidths with an acceptable gain level. Proof-of-principle simulations are also reported along with the successful realization of image processing functions like edge detection. This time-varying wave-based computing system can set the direction for future developments of programmable metasurfaces with highly promising applications in ultrafast equation solving, real-time and continuous signal processing, and imaging.


title: Residual Continual Learning

abstract: We propose a novel continual learning method called Residual Continual Learning (ResCL). Our method can prevent the catastrophic forgetting phenomenon in sequential learning of multiple tasks, without any source task information except the original network. ResCL reparameterizes network parameters by linearly combining each layer of the original network and a fine-tuned network; therefore, the size of the network does not increase at all. To apply the proposed method to general convolutional neural networks, the effects of batch normalization layers are also considered. By utilizing residual-learning-like reparameterization and a special weight decay loss, the trade-off between source and target performance is effectively controlled. The proposed method exhibits state-of-the-art performance in various continual learning scenarios.


title: Micro/Nano Motor Navigation and Localization via Deep Reinforcement Learning

abstract: Efficient navigation and precise localization of Brownian micro/nano self-propelled motor particles within complex landscapes could enable future high-tech applications involving for example drug delivery, precision surgery, oil recovery, and environmental remediation. Here we employ a model-free deep reinforcement learning algorithm based on bio-inspired neural networks to enable different types of micro/nano motors to be continuously controlled to carry out complex navigation and localization tasks. Micro/nano motors with either tunable self-propelling speeds or orientations or both, are found to exhibit strikingly different dynamics. In particular, distinct control strategies are required to achieve effective navigation in free space and obstacle environments, as well as under time constraints. Our findings provide fundamental insights into active dynamics of Brownian particles controlled using artificial intelligence and could guide the design of motor and robot control systems with diverse application requirements.


title: How to 0wn NAS in Your Spare Time

abstract: New data processing pipelines and novel network architectures increasingly drive the success of deep learning. In consequence, the industry considers top-performing architectures as intellectual property and devotes considerable computational resources to discovering such architectures through neural architecture search (NAS). This provides an incentive for adversaries to steal these novel architectures; when used in the cloud, to provide Machine Learning as a Service, the adversaries also have an opportunity to reconstruct the architectures by exploiting a range of hardware side channels. However, it is challenging to reconstruct novel architectures and pipelines without knowing the computational graph (e.g., the layers, branches or skip connections), the architectural parameters (e.g., the number of filters in a convolutional layer) or the specific pre-processing steps (e.g. embeddings). In this paper, we design an algorithm that reconstructs the key components of a novel deep learning system by exploiting a small amount of information leakage from a cache side-channel attack, Flush+Reload. We use Flush+Reload to infer the trace of computations and the timing for each computation. Our algorithm then generates candidate computational graphs from the trace and eliminates incompatible candidates through a parameter estimation process. We implement our algorithm in PyTorch and Tensorflow. We demonstrate experimentally that we can reconstruct MalConv, a novel data pre-processing pipeline for malware detection, and ProxylessNAS- CPU, a novel network architecture for the ImageNet classification optimized to run on CPUs, without knowing the architecture family. In both cases, we achieve 0% error. These results suggest hardware side channels are a practical attack vector against MLaaS, and more efforts should be devoted to understanding their impact on the security of deep learning systems.


title: Fitting ARMA Time Series Models without Identification: A Proximal Approach

abstract: Fitting autoregressive moving average (ARMA) time series models requires model identification before parameter estimation. Model identification involves determining the order of the autoregressive and moving average components which is generally performed by inspection of the autocorrelation and partial autocorrelation functions or other offline methods. In this work, we regularize the parameter estimation optimization problem with a non-smooth hierarchical sparsity-inducing penalty based on two path graphs that allow performing model identification and parameter estimation simultaneously. A proximal block coordinate descent algorithm is then proposed to solve the underlying optimization problem efficiently. The resulting model satisfies the required stationarity and invertibility conditions for ARMA models. Numerical results supporting the proposed method are also presented.


title: Lifter Training and Sub-band Modeling for Computationally Efficient and High-Quality Voice Conversion Using Spectral Differentials

abstract: In this paper, we propose computationally efficient and high-quality methods for statistical voice conversion (VC) with direct waveform modification based on spectral differentials. The conventional method with a minimum-phase filter achieves high-quality conversion but requires heavy computation in filtering. This is because the minimum phase using a fixed lifter of the Hilbert transform often results in a long-tap filter. One of our methods is a data-driven method for lifter training. Since this method takes filter truncation into account in training, it can shorten the tap length of the filter while preserving conversion accuracy. Our other method is sub-band processing for extending the conventional method from narrow-band (16 kHz) to full-band (48 kHz) VC, which can convert a full-band waveform with higher converted-speech quality. Experimental results indicate that 1) the proposed lifter-training method for narrow-band VC can shorten the tap length to 1/16 without degrading the converted-speech quality and 2) the proposed sub-band-processing method for full-band VC can improve the converted-speech quality than the conventional method.


title: Byzantine Lattice Agreement in Asynchronous Systems

abstract: We study the Byzantine lattice agreement (BLA) problem in asynchronous distributed message passing systems. In the BLA problem, each process proposes a value from a join semi-lattice and needs to output a value also in the lattice such that all output values of correct processes lie on a chain despite the presence of Byzantine processes. We present an algorithm for this problem with round complexity of $O(\log f)$ which tolerates $f < \frac{n}{5}$ Byzantine failures in the asynchronous setting without digital signatures, where $n$ is the number of processes. We also show how this algorithm can be modified to work in the authenticated setting (i.e., with digital signatures) to tolerate $f < \frac{n}{3}$ Byzantine failures.


title: Study of silicon sensors for precise timing measurement

abstract: Silicon sensors with high time resolution can help particle identification in the International Linear Collider (ILC). We are studying Low Gain Avalanche Diodes (LGADs) as a high timing resolution sensor. As a step to develop LGADs, we are now focusing to characterize Avalanche Photo Diode (APD)s, because the APDs has the same multiplication structure as LGADs. We studied the characteristics of APDs with particles from radioisotopes.


title: Serpentine optical phased arrays for scalable integrated photonic LIDAR beam steering

abstract: Optical phased arrays (OPAs) implemented in integrated photonic circuits could enable a variety of 3D sensing, imaging, illumination, and ranging applications, and their convergence in new LIDAR technology. However, current integrated OPA approaches do not scale - in control complexity, power consumption, and optical efficiency - to the large aperture sizes needed to support medium to long range LIDAR. We present the serpentine optical phased array (SOPA), a new OPA concept that addresses these fundamental challenges and enables architectures that scale up to large apertures. The SOPA is based on a serially interconnected array of low-loss grating waveguides and supports fully passive, two-dimensional (2D) wavelength-controlled beam steering. A fundamentally space-efficient design that folds the feed network into the aperture also enables scalable tiling of SOPAs into large apertures with a high fill-factor. We experimentally demonstrate the first SOPA, using a 1450 - 1650 nm wavelength sweep to produce 16,500 addressable spots in a 27x610 array. We also demonstrate, for the first time, far-field interference of beams from two separate OPAs on a single silicon photonic chip, as an initial step towards long-range computational imaging LIDAR based on novel active aperture synthesis schemes.


title: Reducing the $H_{0}$ tension with generalized Proca theory

abstract: We investigate the cosmological viability of the generalized proca theory. We first implement the background and linear perturbation equations of motion in the Boltzmann code and then study the constraints on the parameters of the generalized proca theory after running MCMC against the cosmological data set. With Planck + HST data, we obtain the constraint $h=0.7334_{-0.0269}^{+0.0246}$, which indicates that the tension between early universe and late time universe within this theory is removed. By adding other late-time data sets (BAO, RSD, etc.) we show that the tension is reduced, as the 2$\sigma$ allowed region for $h$ in Proca, $h=0.7041_{-0.0087}^{+0.0094}$, overlaps with the 2$\sigma$ region of the HST data.