title: Two's Company, Three's a Crowd: Consensus-Halving for a Constant Number of Agents

abstract: We consider the $\varepsilon$-Consensus-Halving problem, in which a set of heterogeneous agents aim at dividing a continuous resource into two (not necessarily contiguous) portions that all of them simultaneously consider to be of approximately the same value (up to $\varepsilon$). This problem was recently shown to be PPA-complete, for $n$ agents and $n$ cuts, even for very simple valuation functions. In a quest to understand the root of the complexity of the problem, we consider the setting where there is only a constant number of agents, and we consider both the computational complexity and the query complexity of the problem. For agents with monotone valuation functions, we show a dichotomy: for two agents the problem is polynomial-time solvable, whereas for three or more agents it becomes PPA-complete. Similarly, we show that for two monotone agents the problem can be solved with polynomially-many queries, whereas for three or more agents, we provide exponential query complexity lower bounds. These results are enabled via an interesting connection to a monotone Borsuk-Ulam problem, which may be of independent interest. For agents with general valuations, we show that the problem is PPA-complete and admits exponential query complexity lower bounds, even for two agents.


title: Towards a Formal Foundation of Intermittent Computing

abstract: Intermittently powered devices enable new applications in harsh or inaccessible environments, such as space or in-body implants, but also introduce problems in programmability and correctness. Researchers have developed programming models to ensure that programs make progress and do not produce erroneous results due to memory inconsistencies caused by intermittent executions. As the technology has matured, more and more features are added to intermittently powered devices, such as I/O. Prior work has shown that all existing intermittent execution models have problems with repeated device or sensor inputs (RIO). RIOs could leave intermittent executions in an inconsistent state. Such problems and the proliferation of existing intermittent execution models necessitate a formal foundation for intermittent computing. In this paper, we formalize intermittent execution models, their correctness properties with respect to memory consistency and inputs, and identify the invariants needed to prove systems correct. We prove equivalence between several existing intermittent systems. To address RIO problems, we define an algorithm for identifying variables affected by RIOs that need to be restored after reboot and prove the algorithm correct. Finally, we implement the algorithm in a novel intermittent runtime system that is correct with respect to input operations and evaluate its performance.


title: On the connectivity of the disjointness graph of segments of point sets in general position in the plane

abstract: Let $P$ be a set of $n\geq 3$ points in general position in the plane. The edge disjointness graph $D(P)$ of $P$ is the graph whose vertices are all the closed straight line segments with endpoints in $P$, two of which are adjacent in $D(P)$ if and only if they are disjoint. We show that the connectivity of $D(P)$ is at least $\binom{\lfloor\frac{n-2}{2}\rfloor}{2}+\binom{\lceil\frac{n-2}{2}\rceil}{2}$, and that this bound is tight for each $n\geq 3$.


title: Deep Hedging of Long-Term Financial Derivatives

abstract: This study presents a deep reinforcement learning approach for global hedging of long-term financial derivatives. A similar setup as in Coleman et al. (2007) is considered with the risk management of lookback options embedded in guarantees of variable annuities with ratchet features. The deep hedging algorithm of Buehler et al. (2019a) is applied to optimize neural networks representing global hedging policies with both quadratic and non-quadratic penalties. To the best of the author's knowledge, this is the first paper that presents an extensive benchmarking of global policies for long-term contingent claims with the use of various hedging instruments (e.g. underlying and standard options) and with the presence of jump risk for equity. Monte Carlo experiments demonstrate the vast superiority of non-quadratic global hedging as it results simultaneously in downside risk metrics two to three times smaller than best benchmarks and in significant hedging gains. Analyses show that the neural networks are able to effectively adapt their hedging decisions to different penalties and stylized facts of risky asset dynamics only by experiencing simulations of the financial market exhibiting these features. Numerical results also indicate that non-quadratic global policies are significantly more geared towards being long equity risk which entails earning the equity risk premium.


title: Integrating Machine Learning for Planetary Science: Perspectives for the Next Decade

abstract: Machine learning (ML) methods can expand our ability to construct, and draw insight from large datasets. Despite the increasing volume of planetary observations, our field has seen few applications of ML in comparison to other sciences. To support these methods, we propose ten recommendations for bolstering a data-rich future in planetary science.


title: Unnormalized Variational Bayes

abstract: We unify empirical Bayes and variational Bayes for approximating unnormalized densities. This framework, named unnormalized variational Bayes (UVB), is based on formulating a latent variable model for the random variable $Y=X+N(0,\sigma^2 I_d)$ and using the evidence lower bound (ELBO), computed by a variational autoencoder, as a parametrization of the energy function of $Y$ which is then used to estimate $X$ with the empirical Bayes least-squares estimator. In this intriguing setup, the $\textit{gradient}$ of the ELBO with respect to noisy inputs plays the central role in learning the energy function. Empirically, we demonstrate that UVB has a higher capacity to approximate energy functions than the parametrization with MLPs as done in neural empirical Bayes (DEEN). We especially showcase $\sigma=1$, where the differences between UVB and DEEN become visible and qualitative in the denoising experiments. For this high level of noise, the distribution of $Y$ is very smoothed and we demonstrate that one can traverse in a single run $-$ without a restart $-$ all MNIST classes in a variety of styles via walk-jump sampling with a fast-mixing Langevin MCMC sampler. We finish by probing the encoder/decoder of the trained models and confirm UVB $\neq$ VAE.


title: Learning To Pay Attention To Mistakes

abstract: In convolutional neural network based medical image segmentation, the periphery of foreground regions representing malignant tissues may be disproportionately assigned as belonging to the background class of healthy tissues \cite{attenUnet}\cite{AttenUnet2018}\cite{InterSeg}\cite{UnetFrontNeuro}\cite{LearnActiveContour}. This leads to high false negative detection rates. In this paper, we propose a novel attention mechanism to directly address such high false negative rates, called Paying Attention to Mistakes. Our attention mechanism steers the models towards false positive identification, which counters the existing bias towards false negatives. The proposed mechanism has two complementary implementations: (a) "explicit" steering of the model to attend to a larger Effective Receptive Field on the foreground areas; (b) "implicit" steering towards false positives, by attending to a smaller Effective Receptive Field on the background areas. We validated our methods on three tasks: 1) binary dense prediction between vehicles and the background using CityScapes; 2) Enhanced Tumour Core segmentation with multi-modal MRI scans in BRATS2018; 3) segmenting stroke lesions using ultrasound images in ISLES2018. We compared our methods with state-of-the-art attention mechanisms in medical imaging, including self-attention, spatial-attention and spatial-channel mixed attention. Across all of the three different tasks, our models consistently outperform the baseline models in Intersection over Union (IoU) and/or Hausdorff Distance (HD). For instance, in the second task, the "explicit" implementation of our mechanism reduces the HD of the best baseline by more than $26\%$, whilst improving the IoU by more than $3\%$. We believe our proposed attention mechanism can benefit a wide range of medical and computer vision tasks, which suffer from over-detection of background.


title: Coherently amplifying photon production from vacuum with a dense cloud of accelerating photodetectors

abstract: An accelerating photodetector is predicted to see photons in the electromagnetic vacuum. However, the extreme accelerations required have prevented the direct experimental verification of this quantum vacuum effect. In this work, we consider many accelerating photodetectors that are contained within an electromagnetic cavity. We show that the resulting photon production from the cavity vacuum can be collectively enhanced such as to be measurable. The combined cavity-photodetectors system maps onto a parametrically driven Dicke-type model; when the detector number exceeds a certain critical value, the vacuum photon production undergoes a phase transition from a normal phase to an enhanced superradiant-like, inverted lasing phase. Such a model may be realized as a mechanical membrane with a dense concentration of optically active defects undergoing gigahertz flexural motion within a superconducting microwave cavity. We provide estimates suggesting that recent related experimental devices are close to demonstrating this inverted, vacuum photon lasing phase.


title: Algebraic 3D Graphic Statics: Constrained Areas

abstract: This research provides algorithms and numerical methods to geometrically control the magnitude of the internal and external forces in the reciprocal diagrams of 3D/Polyhedral Graphic statics (3DGS). In 3DGS, the form of the structure and its equilibrium of forces is represented by two polyhedral diagrams that are geometrically and topologically related. The areas of the faces of the force diagram represent the magnitude of the internal and external forces in the system. For the first time, the methods of this research allow the user to control and constrain the areas and edge lengths of the faces of general polyhedrons that can be convex, self-intersecting, or concave. As a result, a designer can explicitly control the force magnitudes in the force diagram and explore the equilibrium of a variety of compression and tension-combined funicular structural forms. In this method, a quadratic formulation is used to compute the area of a single face based on its edge lengths. The approach is applied to manipulating the face geometry with a predefined area and the edge lengths. Subsequently, the geometry of the polyhedron is updated with newly changed faces. This approach is a multi-step algorithm where each step includes computing the geometry of a single face and updating the polyhedral geometry. One of the unique results of this framework is the construction of the zero-area, self-intersecting faces, where the sum of the signed areas of a self-intersecting face is zero, representing a member with zero force in the form diagram. The methodology of this research can clarify the equilibrium of some systems that could not be previously justified using reciprocal polyhedral diagrams. Therefore, it generalizes the principle of the equilibrium of polyhedral frames and opens a completely new horizon in the design of highly-sophisticated funicular polyhedral structures beyond compression-only systems.


title: Cosmological constraints with the sub-millimetre galaxies Magnification Bias after large scale bias corrections

abstract: The study of the magnification bias produced on high-redshift sub-millimetre galaxies by foreground galaxies through the analysis of the cross-correlation function was recently demonstrated as an interesting independent alternative to the weak-lensing shear as a cosmological probe. In the case of the proposed observable, most of the cosmological constraints mainly depend on the largest angular separation measurements. Therefore, we aim to study and correct the main large-scale biases that affect foreground and background galaxy samples to produce a robust estimation of the cross-correlation function. Then we analyse the corrected signal to derive updated cosmological constraintsWe measured the large-scale, bias-corrected cross-correlation functions using a background sample of H-ATLAS galaxies with photometric redshifts > 1.2 and two different foreground samples (GAMA galaxies with spectroscopic redshifts or SDSS galaxies with photometric ones, both in the range 0.2 < z < 0.8). These measurements are modelled using the traditional halo model description that depends on both halo occupation distribution and cosmological parameters. We then estimated these parameters by performing a Markov chain Monte Carlo under multiple scenarios to study the performance of this observable and how to improve its results. After the large-scale bias corrections, we obtain only minor improvements with respect to the previous magnification bias results, mainly confirming their conclusions: a lower bound on $\Omega_m > 0.22$ at $95\%$ C.L. and an upper bound $\sigma_8 < 0.97$ at $95\%$ C.L. (results from the $z_{spec}$ sample). However, by combining both foreground samples into a simplified tomographic analysis, we were able to obtain interesting constraints on the $\Omega_m$-$\sigma_8$ plane as follows: $\Omega_m= 0.50_{- 0.20}^{+ 0.14}$ and $\sigma_8= 0.75_{- 0.10}^{+ 0.07}$ at 68\% CL.


title: The Return of Lexical Dependencies: Neural Lexicalized PCFGs

abstract: In this paper we demonstrate that $\textit{context free grammar (CFG) based methods for grammar induction benefit from modeling lexical dependencies}$. This contrasts to the most popular current methods for grammar induction, which focus on discovering $\textit{either}$ constituents $\textit{or}$ dependencies. Previous approaches to marry these two disparate syntactic formalisms (e.g. lexicalized PCFGs) have been plagued by sparsity, making them unsuitable for unsupervised grammar induction. However, in this work, we present novel neural models of lexicalized PCFGs which allow us to overcome sparsity problems and effectively induce both constituents and dependencies within a single model. Experiments demonstrate that this unified framework results in stronger results on both representations than achieved when modeling either formalism alone. Code is available at https://github.com/neulab/neural-lpcfg.


title: Order-Disorder Behavior at Thin Film Oxide Interfaces

abstract: Order-disorder processes fundamentally determine the structure and properties of many important oxide systems for energy and computing applications. While these processes have been intensively studied in bulk materials, they are less investigated and understood for nanostructured oxides in highly non-equilibrium conditions. These systems can now be realized through a range of deposition techniques and probed at exceptional spatial and chemical resolution, leading to a greater focus on interface dynamics. Here we survey a selection of recent studies of order-disorder behavior at thin film oxide interfaces, with a particular emphasis on the emergence of order during synthesis and disorder in extreme irradiation environments. We summarize key trends and identify directions for future study in this growing research area.


title: Real-time detection of a change-point in a linear expectile model

abstract: In the present paper we address the real-time detection problem of a change-point in the coefficients of a linear model with the possibility that the model errors are asymmetrical and that the explanatory variables number is large. We build test statistics based on the cumulative sum (CUSUM) of the expectile function derivatives calculated on the residuals obtained by the expectile and adaptive LASSO expectile estimation methods. The asymptotic distribution of these statistics are obtained under the hypothesis that the model does not change. Moreover, we prove that they diverge when the model changes at an unknown observation. The asymptotic study of the test statistics under these two hypotheses allows us to find the asymptotic critical region and the stopping time, that is the observation where the model will change. The empirical performance is investigated by a comparative simulation study with other statistics of CUSUM type. Two examples on real data are also presented to demonstrate its interest in practice.


title: Sufficient conditions for adiabaticity in open quantum systems

abstract: The adiabatic approximation exhibits wide applicability in quantum mechanics, providing a simple approach for nontransitional dynamics in quantum systems governed by slowly varying time-dependent Hamiltonians. However, the standard adiabatic theorem is specifically derived for closed quantum systems. In a realistic open system scenario, the inevitable system-reservoir interaction must be taken into account, which strongly impacts the generalization of the adiabatic behavior. In this paper, we introduce sufficient conditions for the adiabatic approximation in open quantum systems. These conditions are simple yet general, providing a suitable instrument to investigate adiabaticity for arbitrary initial mixed states evolving under time local master equations. We first illustrate our results by showing that the adiabatic approximation for open systems is compatible with the description of quantum thermodynamics at thermal equilibrium, where irreversible entropy production is vanishing. We also apply our sufficient conditions as a tool in quantum control, evaluating the adiabatic behavior for the Hamiltonians of both the Deutsch algorithm and the Landau-Zener model under decoherence.


title: Deriving Differential Target Propagation from Iterating Approximate Inverses

abstract: We show that a particular form of target propagation, i.e., relying on learned inverses of each layer, which is differential, i.e., where the target is a small perturbation of the forward propagation, gives rise to an update rule which corresponds to an approximate Gauss-Newton gradient-based optimization, without requiring the manipulation or inversion of large matrices. What is interesting is that this is more biologically plausible than back-propagation yet may turn out to implicitly provide a stronger optimization procedure. Extending difference target propagation, we consider several iterative calculations based on local auto-encoders at each layer in order to achieve more precise inversions for more accurate target propagation and we show that these iterative procedures converge exponentially fast if the auto-encoding function minus the identity function has a Lipschitz constant smaller than one, i.e., the auto-encoder is coarsely succeeding at performing an inversion. We also propose a way to normalize the changes at each layer to take into account the relative influence of each layer on the output, so that larger weight changes are done on more influential layers, like would happen in ordinary back-propagation with gradient descent.


title: Computing Optimal Decision Sets with SAT

abstract: As machine learning is increasingly used to help make decisions, there is a demand for these decisions to be explainable. Arguably, the most explainable machine learning models use decision rules. This paper focuses on decision sets, a type of model with unordered rules, which explains each prediction with a single rule. In order to be easy for humans to understand, these rules must be concise. Earlier work on generating optimal decision sets first minimizes the number of rules, and then minimizes the number of literals, but the resulting rules can often be very large. Here we consider a better measure, namely the total size of the decision set in terms of literals. So we are not driven to a small set of rules which require a large number of literals. We provide the first approach to determine minimum-size decision sets that achieve minimum empirical risk and then investigate sparse alternatives where we trade accuracy for size. By finding optimal solutions we show we can build decision set classifiers that are almost as accurate as the best heuristic methods, but far more concise, and hence more explainable.


title: Pairing strategies for the Maker-Breaker game on the hypercube with subcubes as winning sets

abstract: We consider the Maker-Breaker positional game on the vertices of the $n$-dimensional hypercube $\{0,1\}^n$ with $k$-dimensional subcubes as winning sets. We describe a pairing strategy which allows Breaker to win if $n$ is a power of 4 and $k \ge n/4 +1$. Our results also imply that for all $n \geq 3$ there is a Breaker's win pairing strategy if $k \ge \left\lfloor\frac{3}{7}n\right\rfloor +1$.


title: The Eichler integral of $E_2$ and $q$-brackets of $t$-hook functions

abstract: For functions $f: \mathcal{P}\rightarrow \mathbb{C}$ on partitions, Bloch and Okounkov defined a power series $\langle f\rangle_q$ that is the "weighted average" of $f$. As Fourier series in $q=e^{2\pi i z}$, such $q$-brackets generate the ring of quasimodular forms, and the modular forms that are powers of Dedekind's eta-function. Using work of Berndt and Han, we build modular objects from $$ f_t(\lambda):= t\sum_{h\in \mathcal{H}_t(\lambda)}\frac{1}{h^2}, $$ weighted sums over partition hook numbers that are multiples of $t$. We find that $\langle f_t \rangle_q$ is the Eichler integral of $(1-E_2(tz))/24,$ which we modify to construct a function $M_t(z)$ that enjoys weight 0 modularity properties. As a consequence, the non-modular Fourier series $$H_t^*(z):=\sum_{\lambda \in \mathcal{P}} f_t(\lambda)q^{|\lambda|-\frac{1}{24}} $$ inherits weight $-1/2$ modularity properties. These are sufficient to imply a Chowla-Selberg type result, generalizing the fact that weight $k$ algebraic modular forms evaluated at discriminant $D<0$ points $\tau$ are algebraic multiples of $\Omega_D^k,$ the $k$th power of the canonical period. If we let $\Psi(\tau):=-\pi i \left(\frac{\tau^2-3\tau+1}{12\tau}\right)-\frac{\log(\tau)}{2},$ then for $t=1$ we prove that $$ H_1^*(-1/\tau)-\frac{1}{\sqrt{-i\tau}}\cdot H_1^*(\tau)\in \overline{\mathbb{Q}}\cdot \frac{\Psi(\tau)}{\sqrt{\Omega_D}}.$$


title: A splitting theorem for capillary graphs under Ricci lower bounds

abstract: In this paper, we study capillary graphs defined on a domain $\Omega$ of a complete Riemannian manifold $M$, where a graph is said to be capillary if it has constant mean curvature and locally constant Dirichlet and Neumann conditions on $\partial \Omega$. Our main result is a splitting theorem both for $\Omega$ and for the graph function on a class of manifolds with nonnegative Ricci curvature. As a corollary, we classify capillary graphs over domains that are globally Lipschitz epigraphs or slabs in a product space $M = N \times \mathbb{R}$, where $N$ has slow volume growth and non-negative Ricci curvature, including the case $M = \mathbb{R}^2,\mathbb{R}^3$. A technical core of the paper is a new gradient estimate for positive CMC graphs on manifolds with Ricci lower bounds.


title: Single Image Cloud Detection via Multi-Image Fusion

abstract: Artifacts in imagery captured by remote sensing, such as clouds, snow, and shadows, present challenges for various tasks, including semantic segmentation and object detection. A primary challenge in developing algorithms for identifying such artifacts is the cost of collecting annotated training data. In this work, we explore how recent advances in multi-image fusion can be leveraged to bootstrap single image cloud detection. We demonstrate that a network optimized to estimate image quality also implicitly learns to detect clouds. To support the training and evaluation of our approach, we collect a large dataset of Sentinel-2 images along with a per-pixel semantic labelling for land cover. Through various experiments, we demonstrate that our method reduces the need for annotated training data and improves cloud detection performance.


title: Proof of Learning (PoLe): Empowering Machine Learning with Consensus Building on Blockchains

abstract: The progress of deep learning (DL), especially the recent development of automatic design of networks, has brought unprecedented performance gains at heavy computational cost. On the other hand, blockchain systems routinely perform a huge amount of computation that does not achieve practical purposes in order to build Proof-of-Work (PoW) consensus from decentralized participants. In this paper, we propose a new consensus mechanism, Proof of Learning (PoLe), which directs the computation spent for consensus toward optimization of neural networks (NN). In our mechanism, the training/testing data are released to the entire blockchain network (BCN) and the consensus nodes train NN models on the data, which serves as the proof of learning. When the consensus on the BCN considers a NN model to be valid, a new block is appended to the blockchain. We experimentally compare the PoLe protocol with Proof of Work (PoW) and show that PoLe can achieve a more stable block generation rate, which leads to more efficient transaction processing. We also introduce a novel cheating prevention mechanism, Secure Mapping Layer (SML), which can be straightforwardly implemented as a linear NN layer. Empirical evaluation shows that SML can detect cheating nodes at small cost to the predictive performance.


title: Exceptional points and enhanced sensitivity in PT-symmetric continuous elastic media

abstract: We investigate non-Hermitian degeneracies, also known as exceptional points, in continous elastic media, and their potential application to the detection of mass and stiffness perturbations. Degenerate states are induced by enforcing parity-time symmetry through tailored balanced gain and loss, introduced in the form of complex stiffnesses and may be implemented through piezoelectric transducers. Breaking of this symmetry caused by external perturbations leads to a splitting of the eigenvalues, which is explored as a sentitive approach to detection of such perturbations. Numerical simulations on one-dimensional waveguides illustrate the presence of several exceptional points in their vibrational spectrum, and conceptually demonstrate their sensitivity to point mass inclusions. Second order exceptional points are shown to exhibit a frequency shift in the spectrum with a square root dependence on the perturbed mass, which is confirmed by a perturbation approach and by frequency response predictions. Elastic domains supporting guided waves are then investigated, where exceptional points are formed by the hybridization of Lamb wave modes. After illustrating a similar sensitivity to point mass inclusions, we also show how these concepts can be applied to surface wave modes for sensing crack-type defects. The presented results describe fundamental vibrational properties of PT-symmetric elastic media supporting exceptional points, whose sensitivity to perturbations goes beyond the linear dependency commonly encountered in Hermitian systems. The findings are thus promising for applications involving sensing of perturbations such as added masses, stiffness discontinuities and surface cracks.


title: A General Framework For Detecting Anomalous Inputs to DNN Classifiers

abstract: Detecting anomalous inputs, such as adversarial and out-of-distribution (OOD) inputs, is critical for classifiers (including deep neural networks or DNNs) deployed in real-world applications. While prior works have proposed various methods to detect such anomalous samples using information from the internal layer representations of a DNN, there is a lack of consensus on a principled approach for the different components of such a detection method. As a result, often heuristic and one-off methods are applied for different aspects of this problem. We propose an unsupervised anomaly detection framework based on the internal DNN layer representations in the form of a meta-algorithm with configurable components. We proceed to propose specific instantiations for each component of the meta-algorithm based on ideas grounded in statistical testing and anomaly detection. We evaluate the proposed methods on well-known image classification datasets with strong adversarial attacks and OOD inputs, including an adaptive attack that uses the internal layer representations of the DNN (often not considered in prior work). Comparisons with five recently-proposed competing detection methods demonstrates the effectiveness of our method in detecting adversarial and OOD inputs.


title: Quantitative normal approximations for the stochastic fractional heat equation

abstract: In this article we present a {\it quantitative} central limit theorem for the stochastic fractional heat equation driven by a a general Gaussian multiplicative noise, including the cases of space-time white noise and the white-colored noise with spatial covariance given by the Riesz kernel or a bounded integrable function. We show that the spatial average over a ball of radius $R$ converges, as $R$ tends to infinity, after suitable renormalization, towards a Gaussian limit in the total variation distance. We also provide a functional central limit theorem. As such, we extend recently proved similar results for stochastic heat equation to the case of the fractional Laplacian and to the case of general noise.


title: The $L^p$ Teichm\"uller theory: Existence and regularity of critical points

abstract: We study minimisers of the $p$-conformal energy functionals, \[ \mathsf{E}_p(f):=\int_\ID \IK^p(z,f)\,dz,\quad f|_\IS=f_0|_\IS, \] defined for self mappings $f:\ID\to\ID$ with finite distortion and prescribed boundary values $f_0$. Here \[ \IK(z,f) = \frac{\|Df(z)\|^2}{J(z,f)} = \frac{1+|\mu_f(z)|^2}{1-|\mu_f(z)|^2}\] is the pointwise distortion functional and $\mu_f(z)$ is the Beltrami coefficient of $f$. We show that for quasisymmetric boundary data the limiting regimes $p\to\infty$ recover the classical Teichm\"uller theory of extremal quasiconformal mappings (in part a result of Ahlfors), and for $p\to1$ recovers the harmonic mapping theory. Critical points of $\mathsf{E}_p$ always satisfy the inner-variational distributional equation \[ 2p\int_\ID \IK^p\;\frac{\overline{\mu_f}}{1+|\mu_f|^2}\varphi_\zbar \; dz=\int_\ID \IK^p \; \varphi_z\; dz,\quad\forall\varphi\in C_0^\infty(\ID ). \] We establish the existence of minimisers in the {\em a priori} regularity class $W^{1,\frac{2p}{p+1}}(\ID)$ and show these minimisers have a pseudo-inverse - a continuous $W^{1,2}(\ID)$ surjection of $\ID$ with $(h\circ f)(z)=z$ almost everywhere. We then give a sufficient condition to ensure $C^{\infty}(\ID)$ smoothness of solutions to the distributional equation. For instance $\IK(z,f)\in L^r_{loc}(\ID)$ for any $r>p+1$ is enough to imply the solutions to the distributional equation are local diffeomorphisms. Further $\IK(w,h)\in L^1(\ID)$ will imply $h$ is a homeomorphism, and together these results yield a diffeomorphic minimiser. We show such higher regularity assumptions to be necessary for critical points of the inner variational equation.


title: Higher regularity and uniqueness for inner variational equations

abstract: We study local minima of the $p$-conformal energy functionals, \[ \mathsf{E}_{\cal A}^\ast(h):=\int_\ID {\cal A}(\IK(w,h)) \;J(w,h) \; dw,\quad h|_\IS=h_0|_\IS, \] defined for self mappings $h:\ID\to\ID$ with finite distortion of the unit disk with prescribed boundary values $h_0$. Here $\IK(w,h) = \frac{\|Dh(w)\|^2}{J(w,h)} $ is the pointwise distortion functional, and ${\cal A}:[1,\infty)\to [1,\infty)$ is convex and increasing with ${\cal A}(t)\approx t^p$ for some $p\geq 1$, with additional minor technical conditions. Note ${\cal A}(t)=t$ is the Dirichlet energy functional. Critical points of $\mathsf{E}_{\cal A}^\ast$ satisfy the Ahlfors-Hopf inner-variational equation \[ {\cal A}'(\IK(w,h)) h_w \overline{h_\wbar} = \Phi \] where $\Phi$ is a holomorphic function. Iwaniec, Kovalev and Onninen established the Lipschitz regularity of critical points. Here we give a sufficient condition to ensure that a local minimum is a diffeomorphic solution to this equation, and that it is unique. This condition is necessarily satisfied by any locally quasiconformal critical point, and is basically the assumption $\IK(w,h)\in L^1(\ID)\cap L^r_{loc}(\ID)$ for some $r>1$.


title: Fully Dynamic Inference with Deep Neural Networks

abstract: Modern deep neural networks are powerful and widely applicable models that extract task-relevant information through multi-level abstraction. Their cross-domain success, however, is often achieved at the expense of computational cost, high memory bandwidth, and long inference latency, which prevents their deployment in resource-constrained and time-sensitive scenarios, such as edge-side inference and self-driving cars. While recently developed methods for creating efficient deep neural networks are making their real-world deployment more feasible by reducing model size, they do not fully exploit input properties on a per-instance basis to maximize computational efficiency and task accuracy. In particular, most existing methods typically use a one-size-fits-all approach that identically processes all inputs. Motivated by the fact that different images require different feature embeddings to be accurately classified, we propose a fully dynamic paradigm that imparts deep convolutional neural networks with hierarchical inference dynamics at the level of layers and individual convolutional filters/channels. Two compact networks, called Layer-Net (L-Net) and Channel-Net (C-Net), predict on a per-instance basis which layers or filters/channels are redundant and therefore should be skipped. L-Net and C-Net also learn how to scale retained computation outputs to maximize task accuracy. By integrating L-Net and C-Net into a joint design framework, called LC-Net, we consistently outperform state-of-the-art dynamic frameworks with respect to both efficiency and classification accuracy. On the CIFAR-10 dataset, LC-Net results in up to 11.9$\times$ fewer floating-point operations (FLOPs) and up to 3.3% higher accuracy compared to other dynamic inference methods. On the ImageNet dataset, LC-Net achieves up to 1.4$\times$ fewer FLOPs and up to 4.6% higher Top-1 accuracy than the other methods.


title: Accelerating Multi-attribute Unsupervised Seismic Facies Analysis With RAPIDS

abstract: Classification of seismic facies is done by clustering seismic data samples based on their attributes. Year after year, 3D datasets used by exploration geophysics increase in size, complexity, and number of attributes, requiring a continuous rise in the classification performance. In this work, we explore the use of Graphics Processing Units (GPUs) to perform the classification of seismic surveys using the well-established Machine Learning (ML) method k-means. We show that the high-performance distributed implementation of the k-means algorithm available at the RAPIDS library can be used to classify facies in large seismic datasets much faster than a classical parallel CPU implementation (up to 258-fold faster in NVIDIA V100 GPUs), especially for large seismic blocks. We tested the algorithm with different real seismic volumes, including Netherlands, Parihaka, and Kahu (from 12GB to 66GB).


title: Fast, Structured Clinical Documentation via Contextual Autocomplete

abstract: We present a system that uses a learned autocompletion mechanism to facilitate rapid creation of semi-structured clinical documentation. We dynamically suggest relevant clinical concepts as a doctor drafts a note by leveraging features from both unstructured and structured medical data. By constraining our architecture to shallow neural networks, we are able to make these suggestions in real time. Furthermore, as our algorithm is used to write a note, we can automatically annotate the documentation with clean labels of clinical concepts drawn from medical vocabularies, making notes more structured and readable for physicians, patients, and future algorithms. To our knowledge, this system is the only machine learning-based documentation utility for clinical notes deployed in a live hospital setting, and it reduces keystroke burden of clinical concepts by 67% in real environments.


title: Approximate Ridesharing of Personal Vehicles Problem

abstract: The ridesharing problem is that given a set of trips, each trip consists of an individual, a vehicle of the individual and some requirements, select a subset of trips and use the vehicles of selected trips to deliver all individuals to their destinations satisfying the requirements. Requirements of trips are specified by parameters including source, destination, vehicle capacity, preferred paths of a driver, detour distance and number of stops a driver is willing to make, and time constraints. We analyze the relations between the time complexity and parameters for two optimization problems: minimizing the number of selected vehicles and minimizing total travel distance of the vehicles. We consider the following conditions: (1) all trips have the same source or same destination, (2) no detour is allowed, (3) each participant has one preferred path, (4) no limit on the number of stops, and (5) all trips have the same departure and same arrival time. It is known that both minimization problems are NP-hard if one of Conditions (1), (2) and (3) is not satisfied. We prove that both problems are NP-hard and further show that it is NP-hard to approximate both problems within a constant factor if Conditions (4) or (5) is not satisfied. We give $\frac{K+2}{2}$-approximation algorithms for minimizing the number of selected vehicles when condition (4) is not satisfied, where $K$ is the largest capacity of all vehicles.


title: The spectral density of dense random networks and the breakdown of the Wigner semicircle law

abstract: Although the spectra of random networks have been studied for a long time, the influence of network topology on the dense limit of network spectra remains poorly understood. By considering the configuration model of networks with four distinct degree distributions, we show that the spectral density of the adjacency matrices of dense random networks is determined by the strength of the degree fluctuations. In particular, the eigenvalue distribution of dense networks with an exponential degree distribution is governed by a simple equation, from which we uncover a logarithmic singularity in the spectral density. We also derive a relation between the fourth moment of the eigenvalue distribution and the variance of the degree distribution, which leads to a sufficient condition for the breakdown of the Wigner semicircle law for dense random networks. Based on the same relation, we propose a classification scheme of the distinct universal behaviours of the spectral density in the dense limit. Our theoretical findings should lead to important insights on the mean-field behaviour of models defined on graphs.


title: Benchmarking and Comparing Multi-exposure Image Fusion Algorithms

abstract: Multi-exposure image fusion (MEF) is an important area in computer vision and has attracted increasing interests in recent years. Apart from conventional algorithms, deep learning techniques have also been applied to multi-exposure image fusion. However, although much efforts have been made on developing MEF algorithms, the lack of benchmark makes it difficult to perform fair and comprehensive performance comparison among MEF algorithms, thus significantly hindering the development of this field. In this paper, we fill this gap by proposing a benchmark for multi-exposure image fusion (MEFB) which consists of a test set of 100 image pairs, a code library of 16 algorithms, 20 evaluation metrics, 1600 fused images and a software toolkit. To the best of our knowledge, this is the first benchmark in the field of multi-exposure image fusion. Extensive experiments have been conducted using MEFB for comprehensive performance evaluation and for identifying effective algorithms. We expect that MEFB will serve as an effective platform for researchers to compare performances and investigate MEF algorithms.


title: Learning RGB-D Feature Embeddings for Unseen Object Instance Segmentation

abstract: Segmenting unseen objects in cluttered scenes is an important skill that robots need to acquire in order to perform tasks in new environments. In this work, we propose a new method for unseen object instance segmentation by learning RGB-D feature embeddings from synthetic data. A metric learning loss function is utilized to learn to produce pixel-wise feature embeddings such that pixels from the same object are close to each other and pixels from different objects are separated in the embedding space. With the learned feature embeddings, a mean shift clustering algorithm can be applied to discover and segment unseen objects. We further improve the segmentation accuracy with a new two-stage clustering algorithm. Our method demonstrates that non-photorealistic synthetic RGB and depth images can be used to learn feature embeddings that transfer well to real-world images for unseen object instance segmentation.


title: Optimal Local and Remote Controls of Multiple Systems with Multiplicative Noises and Unreliable Uplink Channels

abstract: In this paper, the optimal local and remote linear quadratic (LQ) control problem is studied for a networked control system (NCS) which consists of multiple subsystems and each of which is described by a general multiplicative noise stochastic system with one local controller and one remote controller. Due to the unreliable uplink channels, the remote controller can only access unreliable state information of all subsystems, while the downlink channels from the remote controller to the local controllers are perfect. The difficulties of the LQ control problem for such a system arise from the different information structures of the local controllers and the remote controller. By developing the Pontyagin maximum principle, the necessary and sufficient solvability conditions are derived, which are based on the solution to a group of forward and backward difference equations (G-FBSDEs). Furthermore, by proposing a new method to decouple the G-FBSDEs and introducing new coupled Riccati equations (CREs), the optimal control strategies are derived where we verify that the separation principle holds for the multiplicative noise NCSs with packet dropouts. This paper can be seen as an important contribution to the optimal control problem with asymmetric information structures.


title: Prediction of hierarchical time series using structured regularization and its application to artificial neural networks

abstract: This paper discusses the prediction of hierarchical time series, where each upper-level time series is calculated by summing appropriate lower-level time series. Forecasts for such hierarchical time series should be coherent, meaning that the forecast for an upper-level time series equals the sum of forecasts for corresponding lower-level time series. Previous methods for making coherent forecasts consist of two phases: first computing base (incoherent) forecasts and then reconciling those forecasts based on their inherent hierarchical structure. With the aim of improving time series predictions, we propose a structured regularization method for completing both phases simultaneously. The proposed method is based on a prediction model for bottom-level time series and uses a structured regularization term to incorporate upper-level forecasts into the prediction model. We also develop a backpropagation algorithm specialized for application of our method to artificial neural networks for time series prediction. Experimental results using synthetic and real-world datasets demonstrate the superiority of our method in terms of prediction accuracy and computational efficiency.


title: Asymptotics of sloshing eigenvalues for a triangular prism

abstract: We consider the three-dimensional sloshing problem on a triangular prism whose angles with the sloshing surface are of the form $\frac{\pi}{2q}$, where $q$ is an integer. We are interested in finding a two-term asymptotic expansion of the eigenvalue counting function. When both angles are $\frac{\pi}{4}$, we compute the exact value of the second term. As for the general case, we conjecture an asymptotic expansion by constructing quasimodes for the problem and computing the counting function of the related quasi-eigenvalues. These quasimodes come from solutions of the sloping beach problem and correspond to two kinds of waves, edge waves and surface waves. We show that the quasi-eigenvalues are exponentially close to real eigenvalues of the sloshing problem. The asymptotic expansion of their counting function is closely related to a lattice counting problem inside a perturbed ellipse where the perturbation is in a sense random. The contribution of the angles can then be detected through that perturbation.


title: Rethinking Recurrent Neural Networks and Other Improvements for Image Classification

abstract: Over the long history of machine learning, which dates back several decades, recurrent neural networks (RNNs) have been used mainly for sequential data and time series and generally with 1D information. Even in some rare studies on 2D images, these networks are used merely to learn and generate data sequentially rather than for image recognition tasks. In this study, we propose integrating an RNN as an additional layer when designing image recognition models. We also develop end-to-end multimodel ensembles that produce expert predictions using several models. In addition, we extend the training strategy so that our model performs comparably to leading models and can even match the state-of-the-art models on several challenging datasets (e.g., SVHN (0.99), Cifar-100 (0.9027) and Cifar-10 (0.9852)). Moreover, our model sets a new record on the Surrey dataset (0.949). The source code of the methods provided in this article is available at https://github.com/leonlha/e2e-3m and http://nguyenhuuphong.me.


title: On the spectral asymptotics of waves in periodic media with Dirichlet or Neumann exclusions

abstract: We consider homogenization of the scalar wave equation in periodic media at finite wavenumbers and frequencies, with the focus on continua characterized by: (a) arbitrary Bravais lattice in $\mathbb{R}^d$, $d\!\geqslant\!2$, and (b) exclusions i.e. "voids" that are subject to homogenous (Neumann or Dirichlet) boundary conditions. Making use of the Bloch wave expansion, we pursue this goal via asymptotic ansatz featuring the "spectral distance" from a given wavenumber-eigenfrequency pair (within the first Brillouin zone) as the perturbation parameter. We then introduce the effective wave motion via projection(s) of the scalar wavefield onto the Bloch eigenfunction(s) for the unit cell of periodicity, evaluated at the origin of a spectral neighborhood. For generality, we account for the presence of the source term in the wave equation and we consider -- at a given wavenumber -- generic cases of isolated, repeated, and nearby eigenvalues. In this way we obtain a palette of effective models, featuring both wave- and Dirac-type behaviors, whose applicability is controlled by the local band structure and eigenfunction basis. In all spectral regimes, we pursue the homogenized description up to at least first order of expansion, featuring asymptotic corrections of the homogenized Bloch-wave operator and the homogenized source term. Inherently, such framework provides a convenient platform for the synthesis of a wide range of wave phenomena in metamaterials and phononic crystals. The proposed homogenization framework is illustrated by approximating asymptotically the dispersion relationships for (i) Kagome lattice featuring hexagonal Neumann exclusions, and (ii) "pinned" square lattice with circular Dirichlet exclusions. We complete the numerical portrayal of analytical developments by studying the response of a Kagome lattice due to a dipole-like source term acting near the edge of a band gap.


title: Octonions and the two strictly projective tight 5-designs

abstract: In addition to the vertices of the regular hexagon and icosahedron, there are precisely two strictly projective tight 5-designs: one constructed from the short vectors of the Leech lattice and the other corresponding to a generalized hexagon structure in the octonion projective plane. This paper describes a new connection between these two strictly projective tight 5-designs -- a common construction using octonions. Certain octonion involutionary matrices act on a three-dimensional octonion vector space to produce the first 5-design and these same matrices act on the octonion projective plane to produce the second 5-design. This result uses the octonion construction of the Leech lattice due to Robert Wilson and provides a new link between the generalized hexagon Gh(2,8) and the Leech lattice.


title: Photoinduced dynamics of excitonic order and Rabi oscillation in the two-orbital Hubbard model

abstract: We investigate the condition for the photoinduced enhancement of an excitonic order in a two-orbital Hubbard model, which has been theoretically proposed in our previous work [Phys. Rev. B 97, 115105 (2018)], and analyze it from the viewpoint of the Rabi oscillation. Within the mean-field approximation, we simulate real-time dynamics of an excitonic insulator with a direct gap, where the pair condensation in the initial state is of BEC nature and the photoexcitation is introduced by electric dipole transitions. We first discuss that in the atomic limit our model is reduced to a two-level system that undergoes the Rabi oscillation, so that for single cycle pulses physical quantities after the photoirradiation are essentially determined by the ratio of the Rabi frequency to the pump-light frequency. Then, it is shown that this picture holds even in the case of nonzero transfer integrals where each one-particle state exhibits the Rabi oscillation leading to the enhancement of the excitonic order. We demonstrate that effects of electron-phonon interactions do not alter the results qualitatively. We also examine many-body dynamics by the exact diagonalization method on small clusters, which strongly suggests that our mechanism for the enhancement of the exctionic order survives even when quantum fluctuations are taken into account.


title: Towards a Sustainable Microgrid on Alderney Island Using a Python-based Energy Planning Tool

abstract: In remote or islanded communities, the use of microgrids (MGs) is necessary to ensure electrification and resilience of supply. However, even in small-scale systems, it is computationally and mathematically challenging to design low-cost, optimal, sustainable solutions taking into consideration all the uncertainties of load demands and power generations from renewable energy sources (RESs). This paper uses the open-source Python-based Energy Planning (PyEPLAN) tool, developed for the design of sustainable MGs in remote areas, on the Alderney island, the 3$^{rd}$ largest of the Channel Islands with a population of about 2000 people. A two-stage stochastic model is used to optimally invest in battery storage, solar power, and wind power units. Moreover, the AC power flow equations are modelled by a linearised version of the DistFlow model in PyEPLAN, where the investment variables are here-and-now decisions and not a function of uncertain parameters while the operation variables are wait-and-see decisions and a function of uncertain parameters. The $k$-means clustering technique is used to generate a set of best (risk-seeker), nominal (risk-neutral), and worst (risk-averse) scenarios capturing the uncertainty spectrum using the yearly historical patterns of load demands and solar/wind power generations. The proposed investment planning tool is a mixed-integer linear programming (MILP) model and is coded with Pyomo in PyEPLAN.


title: Symmetry breaking in the double moir\'{e} superlattices of relaxed twisted bilayer graphene on hexagonal boron nitride

abstract: We study the atomic and electronic structures of the commensurate double moir\'{e} superlattices in fully relaxed twisted bilayer graphene (TBG) nearly aligned with the hexagonal boron nitride (BN). The single-particle effective Hamiltonian ($\hat{H}^0$) taking into account the relaxation effect and the full moir\'{e} Hamiltonian introduced by BN has been built for TBG/BN. The mean-field (MF) band structures of the self-consistent Hartree-Fock (SCHF) ground states at different number ($\nu$) of filled flat bands relative to the charge neutrality point (CNP) are obtained based on $\hat{H}^0$ in the plane-wave-like basis. The single-particle flat bands in TBG/BN become separated by the opened gap at CNP due to the symmetry breaking in $\hat{H}^0$. We find that the broken $C_2$ symmetry in $\hat{H}^0$ mainly originates from the intralayer inversion-asymmetric structural deformation in the graphene layer adjacent to BN, which introduces spatially non-uniform modifications of the intralayer Hamiltonian. The gapped flat bands have finite Chern numbers. For TBG/BN with the magic twist angle, the SCHF ground states with $|\nu|$ = 1-3 are all insulating with narrow MF gaps. When the flat conduction bands are filled, the gap at $\nu$ = 1 is smaller than that at $\nu$ = 3, suggesting that the nontrivial topological properties associated with the flat Chern bands are more likely to be observed at $\nu = 3$. This is similar for negative $\nu$ with empty valence bands. The dependence of the electronic structure of TBG/BN on positive $\nu$ is roughly consistent with recent experimental observations.


title: An Improvement for Capsule Networks using Depthwise Separable Convolution

abstract: Capsule Networks face a critical problem in computer vision in the sense that the image background can challenge its performance, although they learn very well on training data. In this work, we propose to improve Capsule Networks' architecture by replacing the Standard Convolution with a Depthwise Separable Convolution. This new design significantly reduces the model's total parameters while increases stability and offers competitive accuracy. In addition, the proposed model on $64\times64$ pixel images outperforms standard models on $32\times32$ and $64\times64$ pixel images. Moreover, we empirically evaluate these models with Deep Learning architectures using state-of-the-art Transfer Learning networks such as Inception V3 and MobileNet V1. The results show that Capsule Networks can perform comparably against Deep Learning models. To the best of our knowledge, we believe that this is the first work on the integration of Depthwise Separable Convolution into Capsule Networks.


title: On fixed point property for $L_p$-representations of Kazhdan groups

abstract: Let $G$ be a topological group with finite Kazhdan set, let $\Omega$ be a standard Borel space and $\mu$ a finite measure on $\Omega$. We prove that for any $p\in [1, \infty)$, any affine isometric action $G \curvearrowright L_p(\Omega, \mu)$ whose linear part arises from an ergodic measure-preserving action $G \curvearrowright (\Omega, \mu)$, has a fixed point.


title: Diving into a simple anguilliform swimmer's sensitivity

abstract: Computational models of aquatic locomotion range from individual modest simple swimmers in 2D to sophisticated 3D individual swimmers to complex multi-swimmer models that attempt to parse collective behavioral dynamics. Each of these models contain a multitude of model input parameters to which the model outputs are inherently dependent, i.e., various swimming performance metrics. In this work, the swimming performance's sensitivity to parameters is investigated for an idealized, simple anguilliform swimming model in 2D. The swimmer considered here propagates forward by dynamically varying its body curvature, similar to motion of a \textit{C. elegan}. The parameter sensitivities were explored with respect to the fluid scale (Reynolds number, Re), stroke (undulation) frequency, as well as a kinematic parameter controlling the velocity and acceleration of each upstroke and downstroke. In total, 5000 fluid-structure interaction simulations were performed, each with a unique parameter combination selected via a Sobol sequence. Thus, global sensitivity analysis was performed using Sobol sensitivity analysis. Results indicate that the swimmer's performance is most sensitive to variation in its stroke frequency. Trends in swimming performance were discovered by projecting the performance metrics from the overall 3D parameter space onto particular 2D subspaces and Pareto-like optimal fronts were identified for swimming efficiency. This work is a natural extension of the parameter explorations of the same model from [1].


title: Counting solvable $S$-unit equations

abstract: We obtain upper bounds on the number of finite sets $\mathcal S$ of primes below a given bound for which various $2$ variable $\mathcal S$-unit equations have a solution.


title: DroneLight: Drone Draws in the Air using Long Exposure Light Painting and ML

abstract: We propose a novel human-drone interaction paradigm where a user directly interacts with a drone to light-paint predefined patterns or letters through hand gestures. The user wears a glove which is equipped with an IMU sensor to draw letters or patterns in the midair. The developed ML algorithm detects the drawn pattern and the drone light-paints each pattern in midair in the real time. The proposed classification model correctly predicts all of the input gestures. The DroneLight system can be applied in drone shows, advertisements, distant communication through text or pattern, rescue, and etc. To our knowledge, it would be the world's first human-centric robotic system that people can use to send messages based on light-painting over distant locations (drone-based instant messaging). Another unique application of the system would be the development of vision-driven rescue system that reads light-painting by person who is in distress and triggers rescue alarm.


title: Skewed link regression models for imbalanced binary response with applications to life insurance

abstract: For a portfolio of life insurance policies observed for a stated period of time, e.g., one year, mortality is typically a rare event. When we examine the outcome of dying or not from such portfolios, we have an imbalanced binary response. The popular logistic and probit regression models can be inappropriate for imbalanced binary response as model estimates may be biased, and if not addressed properly, it can lead to serious adverse predictions. In this paper, we propose the use of skewed link regression models (Generalized Extreme Value, Weibull, and Fre\`chet link models) as more superior models to handle imbalanced binary response. We adopt a fully Bayesian approach for the generalized linear models (GLMs) under the proposed link functions to help better explain the high skewness. To calibrate our proposed Bayesian models, we use a real dataset of death claims experience drawn from a life insurance company's portfolio. Bayesian estimates of parameters were obtained using the Metropolis-Hastings algorithm and for Bayesian model selection and comparison, the Deviance Information Criterion (DIC) statistic has been used. For our mortality dataset, we find that these skewed link models are more superior than the widely used binary models with standard link functions. We evaluate the predictive power of the different underlying models by measuring and comparing aggregated death counts and death benefits.


title: Parametric Dependence of Bound States in the Continuum on Periodic Structures

abstract: Bound states in the continuum (BICs) have some unusual properties and important applications in photonics. A periodic structure sandwiched between two homogeneous media is the most popular platform for observing BICs and realizing their applications. Existing studies on BICs assume the periodic structure has a $C_2$ rotational symmetry about the axis perpendicular to the periodic layer. It is known that all BICs turn to resonant states with finite quality factors if the periodic structure is perturbed by a generic perturbation breaking the $C_2$ symmetry, and a typical BIC continues to exist if the perturbation keeps the $C_2$ symmetry. We study how typical BICs depend on generic structural parameters. For a class of BICs with one opening radiation channel, we show that in the plane of two generic parameters, the BICs exist continuously as a curve. Consequently, BICs can exist on periodic structures without the $C_2$ symmetry, and they can be found by tuning a single structural parameter. The result is established analytically by a perturbation theory with two independent perturbations and validated by numerical examples. Our study reveals a much larger family for BICs on periodic structures, and provides new opportunities for future applications.


title: Learning interaction kernels in stochastic systems of interacting particles from multiple trajectories

abstract: We consider stochastic systems of interacting particles or agents, with dynamics determined by an interaction kernel which only depends on pairwise distances. We study the problem of inferring this interaction kernel from observations of the positions of the particles, in either continuous or discrete time, along multiple independent trajectories. We introduce a nonparametric inference approach to this inverse problem, based on a regularized maximum likelihood estimator constrained to suitable hypothesis spaces adaptive to data. We show that a coercivity condition enables us to control the condition number of this problem and prove the consistency of our estimator, and that in fact it converges at a near-optimal learning rate, equal to the min-max rate of $1$-dimensional non-parametric regression. In particular, this rate is independent of the dimension of the state space, which is typically very high. We also analyze the discretization errors in the case of discrete-time observations, showing that it is of order $1/2$ in terms of the time gaps between observations. This term, when large, dominates the sampling error and the approximation error, preventing convergence of the estimator. Finally, we exhibit an efficient parallel algorithm to construct the estimator from data, and we demonstrate the effectiveness of our algorithm with numerical tests on prototype systems including stochastic opinion dynamics and a Lennard-Jones model.


title: Antiferromagnetic transitions of Dirac fermions in three dimensions

abstract: We use determinant quantum Monte Carlo (DQMC) simulations to study the role of electron-electron interactions on three-dimensional (3D) Dirac fermions based on the $\pi$-flux model on a cubic lattice. We show that the Hubbard interaction drives the 3D Dirac semimetal to an antiferromagnetic (AF) insulator only above a finite critical interaction strength and the long-ranged AF order persists up to a finite temperature. We evaluate the critical interaction strength and temperatures using finite-size scaling of the spin structure factor. The critical behaviors are consistent with the (3+1)d Gross-Neveu universality class for the quantum critical point and 3D Heisenberg universality class for the thermal phase transitions. We further investigate correlation effects in birefringent Dirac fermion system. It is found that the critical interaction strength $U_c$ is decreased by reducing the velocity of the Dirac cone, quantifying the effect of velocity on the critical interaction strength in 3D Dirac fermion systems. Our findings unambiguously uncover correlation effects in 3D Dirac fermions, and may be observed using ultracold atoms in an optical lattice.


title: Domain Adaptive Semantic Segmentation Using Weak Labels

abstract: Learning semantic segmentation models requires a huge amount of pixel-wise labeling. However, labeled data may only be available abundantly in a domain different from the desired target domain, which only has minimal or no annotations. In this work, we propose a novel framework for domain adaptation in semantic segmentation with image-level weak labels in the target domain. The weak labels may be obtained based on a model prediction for unsupervised domain adaptation (UDA), or from a human annotator in a new weakly-supervised domain adaptation (WDA) paradigm for semantic segmentation. Using weak labels is both practical and useful, since (i) collecting image-level target annotations is comparably cheap in WDA and incurs no cost in UDA, and (ii) it opens the opportunity for category-wise domain alignment. Our framework uses weak labels to enable the interplay between feature alignment and pseudo-labeling, improving both in the process of domain adaptation. Specifically, we develop a weak-label classification module to enforce the network to attend to certain categories, and then use such training signals to guide the proposed category-wise alignment method. In experiments, we show considerable improvements with respect to the existing state-of-the-arts in UDA and present a new benchmark in the WDA setting. Project page is at http://www.nec-labs.com/~mas/WeakSegDA.


title: High energy collision without fine tuning: Acceleration and multiple collisions of shells in a bound system

abstract: High energy collision of massive bodies is investigated without fine tuning. We study multiple collisions of two spherical concentric shells in a gravitationally bound system and calculate the center of mass energy between the shells. We solve the equation of motions for two shells without imposing any fine tuning of the initial parameters. In this bound system, the shells collide many times and these motions are highly nontrivial due to chaotic behavior of the shells. Consequently, the center of mass energy for each collision varies nontrivially and even reach almost its theoretical upper limit. We confirm that a significant proportion of the theoretical limit is automatically achieved during multiple collisions without fine tuning. At the same time, we also study shell ejection from the system after some collisions. If the initial shell's energy is large enough, multiple collisions may cause one shell to accumulate energy so that it escapes to infinity, even though two shells are initially confined in the system. The ejection is caused by multiple collisions inducing nontrivial energy transfer between the shells. The relation between the maximum center of mass energy and the energy transfer causing the shell ejection is also discussed.


title: From pure braid groups to hyperbolic groups

abstract: In this note we show that any homomorphism from a pure surface braid group to a torsion-free hyperbolic group either has a cyclic image or factors through a forgetful map. This extends and gives a new proof of an earlier result of the author which works only when the target is a free group or a surface group. We also prove a similar rigidity result for the pure braid group of the disk.


title: Detecting Change Signs with Differential MDL Change Statistics for COVID-19 Pandemic Analysis

abstract: We are concerned with the issue of detecting changes and their signs from a data stream. For example, when given time series of COVID-19 cases in a region, we may raise early warning signals of outbreaks by detecting signs of changes in the cases. We propose a novel methodology to address this issue. The key idea is to employ a new information-theoretic notion, which we call the differential minimum description length change statistics (D-MDL), for measuring the scores of change sign. We first give a fundamental theory for D-MDL. We then demonstrate its effectiveness using synthetic datasets. We apply it to detecting early warning signals of the COVID-19 epidemic. We empirically demonstrate that D-MDL is able to raise early warning signals of events such as significant increase/decrease of cases. Remarkably, for about $64\%$ of the events of significant increase of cases in 37 studied countries, our method can detect warning signals as early as nearly six days on average before the events, buying considerably long time for making responses. We further relate the warning signals to the basic reproduction number $R0$ and the timing of social distancing. The results showed that our method can effectively monitor the dynamics of $R0$, and confirmed the effectiveness of social distancing at containing the epidemic in a region. We conclude that our method is a promising approach to the pandemic analysis from a data science viewpoint. The software for the experiments is available at https://github.com/IbarakikenYukishi/differential-mdl-change-statistics. An online detection system is available at https://ibarakikenyukishi.github.io/d-mdl-html/index.html


title: Zariski density of points with maximal arithmetic degree

abstract: Given a dominant rational self-map on a projective variety over a number field, we can define the arithmetic degree at a rational point. It is known that the arithmetic degree at any point is less than or equal to the first dynamical degree. In this article, we show that there are densely many $\overline{\mathbb Q}$-rational points with maximal arithmetic degree (i.e. whose arithmetic degree is equal to the first dynamical degree) for self-morphisms on projective varieties. For unirational varieties and abelian varieties, we show that there are densely many rational points with maximal arithmetic degree over a sufficiently large number field. We also give a generalization of a result of Kawaguchi and Silverman in the appendix.


title: Model-Based Event-Triggered Control over Lossy Networks

abstract: The event-triggered control problem over lossy communication networks is addressed in this paper. Although packet dropouts have been considered in the implementation of event-triggered controllers, the assumption of protocols that employ acknowledgement messages persists. This paper provides an approach that relaxes such assumption. An event-based controller is implemented at the sensor node of the networked and uncertain system and it transmits feedback measurements to the controller node at asynchronous time instants. The transmitted packets of information are subject to dropouts by the lossy network. We show that the uncertain system can be asymptotically stabilized, that a positive minimum inter-event time exists, and that the proposed approach does not require acknowledgement messages.


title: Visual Analysis of Discrimination in Machine Learning

abstract: The growing use of automated decision-making in critical applications, such as crime prediction and college admission, has raised questions about fairness in machine learning. How can we decide whether different treatments are reasonable or discriminatory? In this paper, we investigate discrimination in machine learning from a visual analytics perspective and propose an interactive visualization tool, DiscriLens, to support a more comprehensive analysis. To reveal detailed information on algorithmic discrimination, DiscriLens identifies a collection of potentially discriminatory itemsets based on causal modeling and classification rules mining. By combining an extended Euler diagram with a matrix-based visualization, we develop a novel set visualization to facilitate the exploration and interpretation of discriminatory itemsets. A user study shows that users can interpret the visually encoded information in DiscriLens quickly and accurately. Use cases demonstrate that DiscriLens provides informative guidance in understanding and reducing algorithmic discrimination.


title: Effect of quantum resonances on local temperature in nonequilibrium open systems

abstract: Measuring local temperatures of open systems out of equilibrium is emerging as a novel approach to study the local thermodynamic properties of nanosystems. An operational protocol has been proposed to determine the local temperature by coupling a probe to the system and then minimizing the perturbation to a certain local observable of the probed system. In this paper, we first show that such a local temperature is unique for a single quantum impurity and the given local observable. We then extend this protocol to open systems consisting of multiple quantum impurities by proposing a local minimal perturbation condition (LMPC). The influence of quantum resonances on the local temperature is elucidated by both analytic and numerical results. In particular, we demonstrate that quantum resonances may give rise to strong oscillations of the local temperature along a multiimpurity chain under a thermal bias.


title: Delta shock solution to the generalized one-dimensional zero-pressure gas dynamics system with linear damping

abstract: In this paper, we propose a time-dependent viscous system and by using the vanishing viscosity method we show the existence of delta shock solution for a particular $2 \times 2$ system of conservation laws with linear damping.


title: Improving probability selecting based weights for Satisfiability Problem

abstract: The Boolean Satisfiability problem (SAT) is important on artificial intelligence community and the impact of its solving on complex problems. Recently, great breakthroughs have been made respectively on stochastic local search (SLS) algorithms for uniform random k-SAT resulting in several state-of-the-art SLS algorithms Score2SAT, YalSAT, ProbSAT, CScoreSAT and on a hybrid algorithm for hard random SAT (HRS) resulting in one state-of-the-art hybrid algorithm SparrowToRiss. However, there is no an algorithm which can effectively solve both uniform random k-SAT and HRS. In this paper, we present a new SLS algorithm named SelectNTS for uniform random k-SAT and HRS. SelectNTS is an improved probability selecting based local search algorithm for SAT problem. The core of SelectNTS relies on new clause and variable selection heuristics. The new clause selection heuristic uses a new clause weighting scheme and a biased random walk. The new variable selection heuristic uses a probability selecting strategy with the variation of CC strategy based on a new variable weighting scheme. Extensive experimental results on the well-known random benchmarks instances from the SAT Competitions in 2017 and 2018, and on randomly generated problems, show that our algorithm outperforms state-of-the-art random SAT algorithms, and our SelectNTS can effectively solve both uniform random k-SAT and HRS.


title: Multiwavelength investigation of extended green object G19.88-0.53: Revealing a protocluster

abstract: A multiwavelength analysis of star formation associated with the extended green object, G19.88-0.53 is presented in this paper. With multiple detected radio and millimetre components, G19.88-0.53 unveils as harbouring a protocluster rather than a single massive young stellar object. We detect an ionized thermal jet using the upgraded Giant Meterwave Radio Telescope, India, which is found to be associated with a massive, dense and hot ALMA 2.7 mm core driving a bipolar CO outflow. Near-infrared spectroscopy with UKIRT-UIST shows the presence of multiple shock-excited H2 lines concurrent with the nature of this region. Detailed investigation of the gas kinematics using ALMA data reveals G19.88-0.53 as an active protocluster with high-mass star forming components spanning a wide evolutionary spectrum from hot cores in accretion phase to cores driving multiple outflows to possible UCHII regions.


title: The Galactic origin for the borders in the Earth history

abstract: The external galactic key factor for developing of massive geochronological boundaries, as well as mass extinctions (ME) is proved based on chronological relationships of the impulses of globally short processes. External galactic key factor is also proved based on the coincidence of the scales of their developing processes, their regular periodicity in the formation of the boundaries of the Phanerozoic epochs and the presence of a unifying complex of cause-effect relationships for boundary processes. The hypotheses of the impact origin of ocean basins and plumes are substantiated, the assumptions about the unsteady nature of the Spiral Arms and the rotational acceleration of the Milky Ways core are confirmed.


title: Developing RNN-T Models Surpassing High-Performance Hybrid Models with Customization Capability

abstract: Because of its streaming nature, recurrent neural network transducer (RNN-T) is a very promising end-to-end (E2E) model that may replace the popular hybrid model for automatic speech recognition. In this paper, we describe our recent development of RNN-T models with reduced GPU memory consumption during training, better initialization strategy, and advanced encoder modeling with future lookahead. When trained with Microsoft's 65 thousand hours of anonymized training data, the developed RNN-T model surpasses a very well trained hybrid model with both better recognition accuracy and lower latency. We further study how to customize RNN-T models to a new domain, which is important for deploying E2E models to practical scenarios. By comparing several methods leveraging text-only data in the new domain, we found that updating RNN-T's prediction and joint networks using text-to-speech generated from domain-specific text is the most effective.


title: Deep Multi-View Spatiotemporal Virtual Graph Neural Network for Significant Citywide Ride-hailing Demand Prediction

abstract: Urban ride-hailing demand prediction is a crucial but challenging task for intelligent transportation system construction. Predictable ride-hailing demand can facilitate more reasonable vehicle scheduling and online car-hailing platform dispatch. Conventional deep learning methods with no external structured data can be accomplished via hybrid models of CNNs and RNNs by meshing plentiful pixel-level labeled data, but spatial data sparsity and limited learning capabilities on temporal long-term dependencies are still two striking bottlenecks. To address these limitations, we propose a new virtual graph modeling method to focus on significant demand regions and a novel Deep Multi-View Spatiotemporal Virtual Graph Neural Network (DMVST-VGNN) to strengthen learning capabilities of spatial dynamics and temporal long-term dependencies. Specifically, DMVST-VGNN integrates the structures of 1D Convolutional Neural Network, Multi Graph Attention Neural Network and Transformer layer, which correspond to short-term temporal dynamics view, spatial dynamics view and long-term temporal dynamics view respectively. In this paper, experiments are conducted on two large-scale New York City datasets in fine-grained prediction scenes. And the experimental results demonstrate effectiveness and superiority of DMVST-VGNN framework in significant citywide ride-hailing demand prediction.


title: Quantitative Understanding of VAE as a Non-linearly Scaled Isometric Embedding

abstract: Variational autoencoder (VAE) estimates the posterior parameters (mean and variance) of latent variables corresponding to each input data. While it is used for many tasks, the transparency of the model is still an underlying issue. This paper provides a quantitative understanding of VAE property through the differential geometric and information-theoretic interpretations of VAE. According to the Rate-distortion theory, the optimal transform coding is achieved by using an orthonormal transform with PCA basis where the transform space is isometric to the input. Considering the analogy of transform coding to VAE, we clarify theoretically and experimentally that VAE can be mapped to an implicit isometric embedding with a scale factor derived from the posterior parameter. As a result, we can estimate the data probabilities in the input space from the prior, loss metrics, and corresponding posterior parameters, and further, the quantitative importance of each latent variable can be evaluated like the eigenvalue of PCA.


title: Plasma dark matter and electronic recoil events in XENON1T

abstract: Dark matter might be in the form of a dark plasma in the Milky Way halo. Specifically, we consider here a hidden sector consisting of a light `dark electron' and a much heavier `dark proton', each charged under an unbroken $U(1)'$ gauge symmetry. These self-interacting dark sector particles can also interact with ordinary matter via the kinetic mixing interaction, and lead to a signal in dark matter direct detection experiments. Indeed, keV electron recoils can arise quite naturally in such models from dark electron scattering off loosely bound atomic electrons. Here we examine the recently reported XENON1T excess in the context of such a plasma dark matter model. We find that the observed excess can be explained if kinetic mixing is in the approximate range: $10^{-12} \lesssim \epsilon \lesssim 10^{-10}$. The allowed parameter space is consistent with astrophysical and cosmological constraints and consistent also with other direct detection experiments.


title: Branch-and-Bound Solves Random Binary IPs in Polytime

abstract: Branch-and-bound is the workhorse of all state-of-the-art mixed integer linear programming (MILP) solvers. These implementations of branch-and-bound typically use variable branching, that is, the child nodes are obtained by fixing some variable to an integer value $v$ in one node and to $v + 1$ in the other node. Even though modern MILP solvers are able to solve very large-scale instances efficiently, relatively little attention has been given to understanding why the underlying branch-and-bound algorithm performs so well. In this paper our goal is to theoretically analyze the performance of the standard variable branching based branch-and-bound algorithm. In order to avoid the exponential worst-case lower bounds, we follow the common idea of considering random instances. More precisely, we consider random integer programs where the entries of the coefficient matrix and the objective function are randomly sampled. Our main result is that with good probability branch-and-bound with variable branching explores only a polynomial number of nodes to solve these instances, for a fixed number of constraints. To the best of our knowledge this is the first known such result for a standard version of branch-and-bound. We believe that this result provides a compelling indication of why branch-and-bound with variable branching works so well in practice.


title: Fermionic Order by Disorder in a van der Waals Antiferromagnet

abstract: CeTe3 is a unique platform to investigate the itinerant magnetism in a van der Waals (vdW) coupled metal. Despite chemical pressure being a promising route to boost quantum fluctuation in this system, a systematic study on the chemical pressure effect on Ce3+(4f1) states is absent. Here, we report on the successful growth of a series of Se doped single crystals of CeTe3. We found a fluctuation driven exotic magnetic rotation from the usual easy-axis ordering to an unusual hard-axis ordering. Unlike in localized magnetic systems, near-critical magnetism can increase itinerancy hand-in-hand with enhancing fluctuation of magnetism. Thus, seemingly unstable hard-axis ordering emerges through kinetic energy gain, with the self-consistent observation of enhanced magnetic fluctuation (disorder). As far as we recognize, this order-by-disorder process in fermionic system is observed for the first time within vdW materials. Our finding opens a unique experimental platform for direct visualization of the rich quasiparticle Fermi surface deformation associated with the Fermionic order-by-disorder process. Also, the search for emergent exotic phases by further tuning of quantum fluctuation is suggested as a promising future challenge.


title: Crowdsampling the Plenoptic Function

abstract: Many popular tourist landmarks are captured in a multitude of online, public photos. These photos represent a sparse and unstructured sampling of the plenoptic function for a particular scene. In this paper,we present a new approach to novel view synthesis under time-varying illumination from such data. Our approach builds on the recent multi-plane image (MPI) format for representing local light fields under fixed viewing conditions. We introduce a new DeepMPI representation, motivated by observations on the sparsity structure of the plenoptic function, that allows for real-time synthesis of photorealistic views that are continuous in both space and across changes in lighting. Our method can synthesize the same compelling parallax and view-dependent effects as previous MPI methods, while simultaneously interpolating along changes in reflectance and illumination with time. We show how to learn a model of these effects in an unsupervised way from an unstructured collection of photos without temporal registration, demonstrating significant improvements over recent work in neural rendering. More information can be found crowdsampling.io.


title: A Vecchia Approximation for High-Dimensional Gaussian Cumulative Distribution Functions Arising from Spatial Data

abstract: We introduce an approach to quickly and accurately approximate the cumulative distribution function of multivariate Gaussian distributions arising from spatial Gaussian processes. This approximation is trivially parallelizable and simple to implement using standard software. We demonstrate its accuracy and computational efficiency in a series of simulation experiments and apply it to analyzing the joint tail of a large precipitation dataset using a recently-proposed scale mixture model for spatial extremes. This dataset is many times larger than what was previously considered possible to fit using preferred inferential techniques.


title: Twist Angle Dependent Interlayer Exciton Lifetimes in van der Waals Heterostructures

abstract: In van der Waals (vdW) heterostructures formed by stacking two monolayers of transition metal dichalcogenides, multiple exciton resonances with highly tunable properties are formed and subject to both vertical and lateral confinement. We investigate how a unique control knob, the twist angle between the two monolayers, can be used to control the exciton dynamics. We observe that the interlayer exciton lifetimes in $\text{MoSe}_{\text{2}}$/$\text{WSe}_{\text{2}}$ twisted bilayers (TBLs) change by one order of magnitude when the twist angle is varied from 1$^\circ$ to 3.5$^\circ$. Using a low-energy continuum model, we theoretically separate two leading mechanisms that influence interlayer exciton radiative lifetimes. The shift to indirect transitions in the momentum space with an increasing twist angle and the energy modulation from the moir\'e potential both have a significant impact on interlayer exciton lifetimes. We further predict distinct temperature dependence of interlayer exciton lifetimes in TBLs with different twist angles, which is partially validated by experiments. While many recent studies have highlighted how the twist angle in a vdW TBL can be used to engineer the ground states and quantum phases due to many-body interaction, our studies explore its role in controlling the dynamics of optically excited states, thus, expanding the conceptual applications of "twistronics".


title: Communication-Efficient Federated Learning via Optimal Client Sampling

abstract: Federated learning (FL) ameliorates privacy concerns in settings where a central server coordinates learning from data distributed across many clients. The clients train locally and communicate the models they learn to the server; aggregation of local models requires frequent communication of large amounts of information between the clients and the central server. We propose a novel, simple and efficient way of updating the central model in communication-constrained settings based on collecting models from clients with informative updates and estimating local updates that were not communicated. In particular, modeling the progression of model's weights by an Ornstein-Uhlenbeck process allows us to derive an optimal sampling strategy for selecting a subset of clients with significant weight updates. The central server collects updated local models from only the selected clients and combines them with estimated model updates of the clients that were not selected for communication. We test this policy on a synthetic dataset for logistic regression and two FL benchmarks, namely, a classification task on EMNIST and a realistic language modeling task using the Shakespeare dataset. The results demonstrate that the proposed framework provides significant reduction in communication while maintaining competitive or achieving superior performance compared to a baseline. Our method represents a new line of strategies for communication-efficient FL that is orthogonal to the existing user-local methods such as quantization or sparsification, thus complementing rather than aiming to replace those existing methods.


title: Formation of "Blanets" from Dust Grains around the Supermassive Black Holes in Galaxies

abstract: In Wada, Tsukamoto, and Kokubo (2019), we proposed for the first time that a new class of planets, "blanets" (i.e., black hole planets), can be formed around supermassive black holes (SMBHs) in the galactic center. Here, we investigate the dust coagulation processes and physical conditions of the blanet formation outside the snow line ($r_{snow} \sim $ several parsecs) in more detail, especially considering the effect of the radial advection of the dust aggregates. We found that a dimensionless parameter $\alpha = v_t^2/c_s^2$, where $v_t$ is the turbulent velocity and $c_s$ is the sound velocity, describing the turbulent viscosity should be smaller than 0.04 in the circumnuclear disk, to prevent the destruction of the aggregates due to collisions. The formation timescale of blanets $\tau_{GI}$ at $r_{snow}$ is, $\tau_{GI} \simeq$ 70-80 Myr for $\alpha = 0.01-0.04$ and $M_{BH} = 10^6 M_\odot$. The mass of the blanets ranges from $\sim 20 M_E$ to $3000 M_E$ in $r < 4$ pc for $\alpha = 0.02$ ($M_E$ is the Earth mass), which is in contrast with $ 4 M_E-6 M_E$ for the case without the radial advection. Our results suggest that \textit{blanets} could be formed around relatively low-luminosity AGNs ($L_{bol} \sim 10^{42}$ erg s$^{-1}$) during their lifetime ($\lesssim 10^8$ yr).


title: Murchison Widefield Array detection of steep-spectrum, diffuse, non-thermal radio emission within Abell 1127

abstract: Diffuse, non-thermal emission in galaxy clusters is increasingly being detected in low-frequency radio surveys and images. We present a new diffuse, steep-spectrum, non-thermal radio source within the cluster Abell 1127 found in survey data from the Murchison Widefield Array (MWA). We perform follow-up observations with the 'extended' configuration MWA Phase II with improved resolution to better resolve the source and measure its low-frequency spectral properties. We use archival Very Large Array S-band data to remove the discrete source contribution from the MWA data, and from a power law model fit we find a spectral index of $-1.83\pm0.29$ broadly consistent with relic-type sources. The source is revealed by the Giant Metrewave Radio Telescope (GMRT) at 150 MHz to have an elongated morphology, with a projected linear size of 850 kpc as measured in the MWA data. Using Chandra observations we derive morphological estimators and confirm quantitatively that the cluster is in a disturbed dynamical state, consistent with the majority of phoenices and relics being hosted by merging clusters. We discuss the implications of relying on morphology and low-resolution imaging alone for the classification of such sources and highlight the usefulness of the MHz to GHz radio spectrum in classifying these types of emission. Finally, we discuss the benefits and limitations of using the MWA Phase II in conjunction with other instruments for detailed studies of diffuse, steep-spectrum, non-thermal radio emission within galaxy clusters.


title: Higher-order Topological Anderson Insulators

abstract: We study disorder effects in a two-dimensional system with chiral symmetry and find that disorder can induce a quadrupole topological insulating phase (a higher-order topological phase with quadrupole moments) from a topologically trivial phase. Their topological properties manifest in a topological invariant defined based on effective boundary Hamiltonians, the quadrupole moment, and zero-energy corner modes. We find gapped and gapless topological phases and a Griffiths regime. In the gapless topological phase, all the states are localized, while in the Griffiths regime, the states at zero energy become multifractal. We further apply the self-consistent Born approximation to show that the induced topological phase arises from disorder renormalized masses. We finally introduce a practical experimental scheme with topoelectrical circuits where the predicted topological phenomena can be observed by impedance measurements. Our work opens the door to studying higher-order topological Anderson insulators and their localization properties.


title: A Finsler type Lipschitz optimal transport metric for a quasilinear wave equation

abstract: We consider the global well-posedness of weak energy conservative solution to a general quasilinear wave equation through variational principle, where the solution may form finite time cusp singularity, when energy concentrates. As a main result in this paper, we construct a Finsler type optimal transport metric, then prove that the solution flow is Lipschitz under this metric. We also prove a generic regularity result by applying Thom's transversality theorem, then find piecewise smooth transportation paths among a dense set of solutions. The results in this paper are for large data solutions, without restriction on the size of solutions.


title: Third-Order Statistics Reconstruction from Compressive Measurements

abstract: Estimation of third-order statistics relies on the availability of a huge amount of data records, which can pose severe challenges on the data collecting hardware in terms of considerable storage costs, overwhelming energy consumption, and unaffordably high sampling rate especially when dealing with high-dimensional data such as wideband signals. To overcome these challenges, this paper focuses on the reconstruction of the third-order cumulants under the compressive sensing framework. Specifically, this paper derives a transformed linear system that directly connects the cross-cumulants of compressive measurements to the desired third-order statistics. We provide sufficient conditions for lossless third-order statistics reconstruction via solving simple least-squares, along with the strongest achievable compression ratio. To reduce the computational burden, we also propose an approach to recover diagonal cumulant slices directly from compressive measurements, which is useful when the cumulant slices are sufficient for the inference task at hand. All the proposed techniques are tested via extensive simulations. The developed joint sampling and reconstruction approach to third-order statistics estimation is able to reduce the required sampling rates significantly by exploiting the cumulant structure resulting from signal stationarity, even in the absence of any sparsity constraints on the signal or cumulants.


title: Algorithmic Stability in Fair Allocation of Indivisible Goods Among Two Agents

abstract: Many allocation problems in multiagent systems rely on agents specifying cardinal preferences. However, allocation mechanisms can be sensitive to small perturbations in cardinal preferences, thus causing agents who make ``small" or ``innocuous" mistakes while reporting their preferences to experience a large change in their utility for the final outcome. To address this, we introduce a notion of algorithmic stability and study it in the context of fair and efficient allocations of indivisible goods among two agents. We show that it is impossible to achieve exact stability along with even a weak notion of fairness and even approximate efficiency. As a result, we propose two relaxations to stability, namely, approximate-stability and weak-approximate-stability, and show how existing algorithms in the fair division literature that guarantee fair and efficient outcomes perform poorly with respect to these relaxations. This leads us to explore the possibility of designing new algorithms that are more stable. Towards this end, we present a general characterization result for pairwise maximin share allocations, and in turn use it to design an algorithm that is approximately-stable and guarantees a pairwise maximin share and Pareto optimal allocation for two agents. Finally, we present a simple framework that can be used to modify existing fair and efficient algorithms in order to ensure that they also achieve weak-approximate-stability.


title: ISS Estimates in the Spatial Sup-Norm for Nonlinear 1-D Parabolic PDEs

abstract: This paper provides novel Input-to-State Stability (ISS)-style maximum principle estimates for classical solutions of highly nonlinear 1-D parabolic Partial Differential Equations (PDEs). The derivation of the ISS-style maximum principle estimates is performed by using an ISS Lyapunov Functional for the sup norm. The estimates provide fading memory ISS estimates in the sup norm of the state with respect to distributed and boundary inputs. The obtained results can handle parabolic PDEs with nonlinear and non-local in-domain terms/boundary conditions. Three illustrative examples show the efficiency of the proposed methodology for the derivation of ISS estimates in the sup norm of the state.


title: Cognition coming about: self-organisation and free-energy

abstract: Wright and Bourkes compelling article rightly points out that existing models of embryogenesis fail to explain the mechanisms and functional significance of the dynamic connections among neurons. We pursue their account of Dynamic Logic by appealing to the Markov blanket formalism that underwrites the Free Energy Principle. We submit that this allows one to model embryogenesis as self-organisation in a dynamical system that minimises free-energy. The ensuing formalism may be extended to also explain the autonomous emergence of cognition, specifically in the brain, as a dynamic self-assembling process.