title: Statistical inference for heavy tailed series with extremal independence

abstract: We consider stationary time series $\{X_j, j \in Z\} whose finite dimensional distributions are regularly varying with extremal independence. We assume that for each $h \geq 1$, conditionally on $X_0$ to exceed a threshold tending to infinity, the conditional distribution of $X_h$ suitably normalized converges weakly to a non degenerate distribution. We consider in this paper the estimation of the normalization and of the limiting distribution.


title: On Optimally Partitioning Variable-Byte Codes

abstract: The ubiquitous Variable-Byte encoding is one of the fastest compressed representation for integer sequences. However, its compression ratio is usually not competitive with other more sophisticated encoders, especially when the integers to be compressed are small that is the typical case for inverted indexes. This paper shows that the compression ratio of Variable-Byte can be improved by 2x by adopting a partitioned representation of the inverted lists. This makes Variable-Byte surprisingly competitive in space with the best bit-aligned encoders, hence disproving the folklore belief that Variable-Byte is space-inefficient for inverted index compression. Despite the significant space savings, we show that our optimization almost comes for free, given that: we introduce an optimal partitioning algorithm that does not affect indexing time because of its linear-time complexity; we show that the query processing speed of Variable-Byte is preserved, with an extensive experimental analysis and comparison with several other state-of-the-art encoders.


title: A Robust Wald-type Test for Testing the Equality of Two Means from Log-Normal Samples

abstract: The log-normal distribution is one of the most common distributions used for modeling skewed and positive data. It frequently arises in many disciplines of science, specially in the biological and medical sciences. The statistical analysis for comparing the means of two independent log-normal distributions is an issue of significant interest. In this paper we present a robust test for this problem. The unknown parameters of the model are estimated by minimum density power divergence estimators (Basu et al 1998, Biometrika, 85(3), 549-559). The robustness as well as the asymptotic properties of the proposed test statistics are rigorously established. The performance of the test is explored through simulations and real data analysis. The test is compared with some existing methods, and it is demonstrated that the proposed test outperforms the others in the presence of outliers.


title: Hairy black holes of Lovelock-Born-Infeld-Scalar gravity

abstract: In this paper we construct new hairy black holes of Lovelock-Born-Infeld-Scalar gravity in diverse dimensions. We study horizon structures of the black holes and find the black holes may possess three horizons, two horizons or one horizon. We give the thermodynamic quantities and check the first laws. We study the thermal phase transition behaviors of the black holes in $T$-$S$ plane, and find for even-dimensional black holes, in some ranges of the parameters Van der Waals-like phase transitions occur, in some ranges of the parameters, reentrant phase transitions (RPTs) appear. For odd-dimensional black holes no such phase transitions are found.


title: Density of uniqueness triples from the diamond axiom

abstract: We work with a pre-$\lambda$-frame, which is an abstract elementary class (AEC) endowed with a collection of basic types and a non-forking relation satisfying certain natural properties with respect to models of cardinality $\lambda$. We investigate the density of uniqueness triples in a given pre-$\lambda$-frame $\mathfrak s$, that is, under what circumstances every basic triple admits a non-forking extension that is a uniqueness triple. Prior results in this direction required strong hypotheses on $\mathfrak s$. Our main result is an improvement, in that we assume far fewer hypotheses on $\mathfrak s$. In particular, we do not require $\mathfrak s$ to satisfy the extension, uniqueness, stability, or symmetry properties, or any form of local character, though we do impose the amalgamation and stability properties in $\lambda^+$, and we do assume $\diamondsuit(\lambda^+)$. As a corollary, by applying our main result to the trivial $\lambda$-frame, it follows that in any AEC $\mathbf K$ satisfying modest hypotheses on $\mathbf K_\lambda$ and $\mathbf K_{\lambda^+}$, the set of $*$-domination triples in $\mathbf K_\lambda$ is dense among the non-algebraic triples. We also apply our main result to the non-splitting relation, obtaining the density of uniqueness triples from very few hypotheses.


title: Embedding with a Rigid Substructure

abstract: This paper presents a new distance geometry algorithm for calculating atomic coordinates from estimates of the interatomic distances, which maintains the positions of the atoms in a known rigid substructure. Given an $M \times 3$ matrix of coordinates for the rigid substructure $\mathbf X$, this problem consists of finding the $N \times 3$ matrix $\mathbf Y$ that yields of global minimum of the so-called STRAIN, i.e. \[ \min_{\mathbf Y} \left\| \begin{bmatrix} \mathbf{XX}^\top & \mathbf{XY}^\top \\ \mathbf{YX}^\top & \mathbf{YY}^\top \end{bmatrix} \,-\, \begin{bmatrix} \mathbf A & \mathbf B \\ \mathbf B^\top & \mathbf C \end{bmatrix} \right\|_{\mathsf F}^2 ~, \] where $\mathbf A = \mathbf{XX}^\top$ , and $\mathbf B, \mathbf C$ are matrices of inner products calculated from the estimated distances. The vanishing of the gradient of the STRAIN is shown to be equivalent to a system of only six nonlinear equations in six unknowns for the inertial tensor associated with the solution Y . The entire solution space is characterized in terms of the geometry of the intersection curves between the unit sphere and certain variable ellipsoids. Upon deriving tight bilateral bounds on the moments of inertia of any possible solution, we construct a search procedure that reliably locates the global minimum. The effectiveness of this method is demonstrated on realistic simulated and chemical test problems.


title: Chaotic blowup in the 3D incompressible Euler equations on a logarithmic lattice

abstract: The dispute on whether the three-dimensional (3D) incompressible Euler equations develop an infinitely large vorticity in a finite time (blowup) keeps increasing due to ambiguous results from state-of-the-art direct numerical simulations (DNS), while the available simplified models fail to explain the intrinsic complexity and variety of observed structures. Here, we propose a new model formally identical to the Euler equations, by imitating the calculus on a 3D logarithmic lattice. This model clarifies the present controversy at the scales of existing DNS and provides the unambiguous evidence of the following transition to the blowup, explained as a chaotic attractor in a renormalized system. The chaotic attractor spans over the anomalously large six-decade interval of spatial scales. For the original Euler system, our results suggest that the existing DNS strategies at the resolution accessible now (and presumably rather long into the future) are unsuitable, by far, for the blowup analysis, and establish new fundamental requirements for the approach to this long-standing problem.


title: Two new Super Li-rich core He-burning giants: A new twist to the long tale of Li enhancement in K giants

abstract: In this paper we report two new super Li-rich K giants: KIC 2305930 and KIC 12645107 with Li abundances exceeding that of ISM (A(Li) $\geq$ 3.2 dex). Importantly, both the giants have been classified as core He-burning red clump stars based on asteroseismic data from Kepler mission. Also, both the stars are found to be low mass (M $\approx$ 1.0 M$_{\odot}$) which, together with an evidence of their evolutionary status of being red clump imply that the stars have gone through both the luminosity bump and He-flash during their RGB evolution. Stars' large Li abundance and evolutionary phase suggest that Li enrichment occurred very recently probably at the tip of RGB either during He-flash, an immediate preceding event on RGB, or by some kind of external event such as merger of RGB star with white dwarf. The findings will provide critical constraints to theoretical models for understanding of Li enhancement origin in RGB stars.


title: The Regularity Problem for Lie Groups with Asymptotic Estimate Lie Algebras

abstract: We solve the regularity problem for Milnor's infinite dimensional Lie groups in the asymptotic estimate context. Specifically, let $G$ be a Lie group with asymptotic estimate Lie algebra $\mathfrak{g}$, and denote its evolution map by $\mathrm{evol}\colon \mathrm{D}\equiv \mathrm{dom}[\mathrm{evol}]\rightarrow G$, i.e., $\mathrm{D}\subseteq C^0([0,1],\mathfrak{g})$. We show that $\mathrm{evol}$ is $C^\infty$-continuous on $\mathrm{D}\cap C^\infty([0,1],\mathfrak{g})$ if and only if $\mathrm{evol}$ is $C^0$-continuous on $\mathrm{D}\cap C^0([0,1],\mathfrak{g})$. We furthermore show that $G$ is k-confined for $k\in \mathbb{N}\sqcup\{\mathrm{lip},\infty\}$ if $G$ is constricted. (The latter condition is slightly less restrictive than to be asymptotic estimate.) Results obtained in a previous paper then imply that an asymptotic estimate Lie group $G$ is $C^\infty$-regular if and only if it is Mackey complete, locally $\mu$-convex, and has Mackey complete Lie algebra - In this case, $G$ is $C^k$-regular for each $k\in \mathbb{N}_{\geq 1}\sqcup\{\mathrm{lip},\infty\}$ (with ``smoothness restrictions'' for $k\equiv\mathrm{lip}$), as well as $C^0$-regular if $G$ is even sequentially complete with integral complete Lie algebra.


title: Interpreting Quantile Independence

abstract: How should one assess the credibility of assumptions weaker than statistical independence, like quantile independence? In the context of identifying causal effects of a treatment variable, we argue that such deviations should be chosen based on the form of selection on unobservables they allow. For quantile independence, we characterize this form of treatment selection. Specifically, we show that quantile independence is equivalent to a constraint on the average value of either a latent propensity score (for a binary treatment) or the cdf of treatment given the unobservables (for a continuous treatment). In both cases, this average value constraint requires a kind of non-monotonic treatment selection. Using these results, we show that several common treatment selection models are incompatible with quantile independence. We introduce a class of assumptions which weakens quantile independence by removing the average value constraint, and therefore allows for monotonic treatment selection. In a potential outcomes model with a binary treatment, we derive identified sets for the ATT and QTT under both classes of assumptions. In a numerical example we show that the average value constraint inherent in quantile independence has substantial identifying power. Our results suggest that researchers should carefully consider the credibility of this non-monotonicity property when using quantile independence to weaken full independence.


title: Non-commutative integrability, exact solvability and the Hamilton-Jacobi theory

abstract: The non-commutative integrability (NCI) is a property fulfilled by some Hamiltonian systems that ensures, among other things, the exact solvability of their corresponding equations of motion. The latter means that an "explicit formula" for the trajectories of these systems can be constructed. Such a construction rests mainly on the so-called Lie theorem on integrability by quadratures. It is worth mentioning that, in the context of Hamiltonian systems, the NCI has been for around 40 years, essentially, the unique criterium for exact solvability expressed in the terms of first integrals (containing the usual Liouville-Arnold integrability criterium as a particular case). Concretely, a Hamiltonian system with $n$ degrees of freedom is said to be non-commutative integrable if a set of independent first integrals $F_{1},...,F_{l}$ are known such that: the kernel of the $l\times l$ matrix with coefficients $\left\{ F_{i},F_{j}\right\} $, where $\left\{ \cdot,\cdot\right\} $ denotes the canonical Poisson bracket, has dimension $2n-l$ (isotropy); and each bracket $\left\{ F_{i},F_{j}\right\} $ is functionally dependent on $F_{1},...,F_{l}$ (closure). In this paper, we develop two procedures for constructing the trajectories of a Hamiltonian system which only require isotropic first integrals (closure condition is not needed). One of them is based on an extended version of the geometric Hamilton-Jacobi theory, and does not rely on the above mentioned Lie's theorem. We do all that in the language of functions of several variables.


title: Subword Regularization: Improving Neural Network Translation Models with Multiple Subword Candidates

abstract: Subword units are an effective way to alleviate the open vocabulary problems in neural machine translation (NMT). While sentences are usually converted into unique subword sequences, subword segmentation is potentially ambiguous and multiple segmentations are possible even with the same vocabulary. The question addressed in this paper is whether it is possible to harness the segmentation ambiguity as a noise to improve the robustness of NMT. We present a simple regularization method, subword regularization, which trains the model with multiple subword segmentations probabilistically sampled during training. In addition, for better subword sampling, we propose a new subword segmentation algorithm based on a unigram language model. We experiment with multiple corpora and report consistent improvements especially on low resource and out-of-domain settings.


title: Generating Interpretable Fuzzy Controllers using Particle Swarm Optimization and Genetic Programming

abstract: Autonomously training interpretable control strategies, called policies, using pre-existing plant trajectory data is of great interest in industrial applications. Fuzzy controllers have been used in industry for decades as interpretable and efficient system controllers. In this study, we introduce a fuzzy genetic programming (GP) approach called fuzzy GP reinforcement learning (FGPRL) that can select the relevant state features, determine the size of the required fuzzy rule set, and automatically adjust all the controller parameters simultaneously. Each GP individual's fitness is computed using model-based batch reinforcement learning (RL), which first trains a model using available system samples and subsequently performs Monte Carlo rollouts to predict each policy candidate's performance. We compare FGPRL to an extended version of a related method called fuzzy particle swarm reinforcement learning (FPSRL), which uses swarm intelligence to tune the fuzzy policy parameters. Experiments using an industrial benchmark show that FGPRL is able to autonomously learn interpretable fuzzy policies with high control performance.


title: Simultaneous Parameter Learning and Bi-Clustering for Multi-Response Models

abstract: We consider multi-response and multitask regression models, where the parameter matrix to be estimated is expected to have an unknown grouping structure. The groupings can be along tasks, or features, or both, the last one indicating a bi-cluster or "checkerboard" structure. Discovering this grouping structure along with parameter inference makes sense in several applications, such as multi-response Genome-Wide Association Studies. This additional structure can not only can be leveraged for more accurate parameter estimation, but it also provides valuable information on the underlying data mechanisms (e.g. relationships among genotypes and phenotypes in GWAS). In this paper, we propose two formulations to simultaneously learn the parameter matrix and its group structures, based on convex regularization penalties. We present optimization approaches to solve the resulting problems and provide numerical convergence guarantees. Our approaches are validated on extensive simulations and real datasets concerning phenotypes and genotypes of plant varieties.


title: Stress anisotropy in shear-jammed packings of frictionless disks

abstract: We perform computational studies of repulsive, frictionless disks to investigate the development of stress anisotropy in mechanically stable (MS) packings. We focus on two protocols for generating MS packings: 1) isotropic compression and 2) applied simple or pure shear strain $\gamma$ at fixed packing fraction $\phi$. MS packings of frictionless disks occur as geometric families (i.e. parabolic segments with positive curvature) in the $\phi$-$\gamma$ plane. MS packings from protocol 1 populate parabolic segments with both signs of the slope, $d\phi/d\gamma >0$ and $d\phi/d\gamma <0$. In contrast, MS packings from protocol 2 populate segments with $d\phi/d\gamma <0$ only. For both simple and pure shear, we derive a relationship between the stress anisotropy and dilatancy $d\phi/d\gamma$ obeyed by MS packings along geometrical families. We show that for MS packings prepared using isotropic compression, the stress anisotropy distribution is Gaussian centered at zero with a standard deviation that decreases with increasing system size. For shear jammed MS packings, the stress anisotropy distribution is a convolution of Weibull distributions that depend on strain, which has a nonzero average and standard deviation in the large-system limit. We also develop a framework to calculate the stress anisotropy distribution for packings generated via protocol 2 in terms of the stress anisotropy distribution for packings generated via protocol 1. These results emphasize that for repulsive frictionless disks, different packing-generation protocols give rise to different MS packing probabilities, which lead to differences in macroscopic properties of MS packings.


title: Some $q$-congruences with parameters

abstract: Let $\Phi_n(q)$ be the $n$-th cyclotomic polynomial in $q$. Recently, the author and Zudilin provide a creative microscoping method to prove some $q$-supercongruences mainly modulo $\Phi_n(q)^3$ by introducing an additional parameter $a$. In this paper, we use this creative microscoping method to confirm some conjectures on $q$-supercongruences modulo $\Phi_n(q)^2$. We also give some parameter-generalizations of known $q$-supercongruences. For instance, we present further generalizations of a $q$-analogue of a famous supercongruence of Rodriguez-Villegas: $$ \sum_{k=0}^{p-1}\frac{{2k\choose k}^2}{16^k} \equiv (-1)^{(p-1)/2}\pmod{p^2}\quad\text{for any odd prime $p$.} $$


title: The Bayesian optimist's guide to adaptive immune receptor repertoire analysis

abstract: Probabilistic modeling is fundamental to the statistical analysis of complex data. In addition to forming a coherent description of the data-generating process, probabilistic models enable parameter inference about given data sets. This procedure is well-developed in the Bayesian perspective, in which one infers probability distributions describing to what extent various possible parameters agree with the data. In this paper we motivate and review probabilistic modeling for adaptive immune receptor repertoire data then describe progress and prospects for future work, from germline haplotyping to adaptive immune system deployment across tissues. The relevant quantities in immune sequence analysis include not only continuous parameters such as gene use frequency, but also discrete objects such as B cell clusters and lineages. Throughout this review, we unravel the many opportunities for probabilistic modeling in adaptive immune receptor analysis, including settings for which the Bayesian approach holds substantial promise (especially if one is optimistic about new computational methods). From our perspective the greatest prospects for progress in probabilistic modeling for repertoires concern ancestral sequence estimation for B cell receptor lineages, including uncertainty from germline genotype, rearrangement, and lineage development.


title: How aggressive a driver is? - A quantitative analysis

abstract: Consider a bottleneck in a road through which only one car can pass through. Suppose that at a time the car passing will have the most aggressive driver in queue and that the aggressiveness of an individual is measured by an attribute $A \equiv N \tau^{\sigma}$ where the quantity $N$ varies randomly from person to person in the range 0 to 1, $\tau$ is the time for which the driver is waiting in the bottleneck and the parameter $\sigma$ is the same for all individuals. Thus, we assume that the aggressiveness depends on the nature of the individual and increases with waiting time in a traffic jam. In support of the algebraic form of $A$, we show (numerically and analytically) that our hypothesis implies that the probability of waiting for a time $\tau$ will be $P(\tau) \propto \tau^{\alpha}$ with the value of $\alpha$ fixed by $\sigma$. Empirical studies confirm such variation in $P(\tau)$ with an exponent of 3.0 to 3.5 in two different cities of India and 1.5 for a traffic intersection in Germany. There is a possibility that the parameter $\sigma$ (and hence $\alpha$) is characteristic of a geographical region.


title: Intrinsic energy spread and bunch length growth in plasma-based accelerators due to betatron motion

abstract: Plasma-based accelerators (PBAs), having demonstrated the production of GeV electron beams in only centimetre scales, offer a path towards a new generation of highly compact and cost-effective particle accelerators. However, achieving the required beam quality, particularly on the energy spread for applications such as free-electron lasers, remains a challenge. Here we investigate fundamental sources of energy spread and bunch length in PBAs which arise from the betatron motion of beam electrons. We present an analytical theory, validated against particle-in-cell simulations, which accurately describes these phenomena. Significant impact on the beam quality is predicted for certain configurations, explaining previously observed limitations on the achievable bunch length and energy spread. Guidelines for mitigating these contributions towards high-quality beams are deduced.


title: $T\bar T$-deformations in closed form

abstract: We consider the problem of exact integration of the $T\bar{T}$-deformation of two dimensional quantum field theories, as well as some higher dimensional extensions in the form of $\det T$-deformations. When the action can be shown to only depend algebraically on the background metric the solution of the deformation equation on the Lagrangian can be given in closed form in terms of solutions of the (extended) Burgers' equation. We present such examples in two and higher dimensions.


title: Ramsey's theorem and products in the Weihrauch degrees

abstract: We study the positions in the Weihrauch lattice of parallel products of various combinatorial principles related to Ramsey's theorem. Among other results, we obtain an answer to a question of Brattka, by showing that Ramsey's theorem for pairs ($\mathsf{RT}^2_2$) is strictly Weihrauch below the parallel product of the stable Ramsey's theorem for pairs and the cohesive principle ($\mathsf{SRT}^2_2 \times \mathsf{COH}$).


title: UNIQ: Uniform Noise Injection for Non-Uniform Quantization of Neural Networks

abstract: We present a novel method for neural network quantization that emulates a non-uniform $k$-quantile quantizer, which adapts to the distribution of the quantized parameters. Our approach provides a novel alternative to the existing uniform quantization techniques for neural networks. We suggest to compare the results as a function of the bit-operations (BOPS) performed, assuming a look-up table availability for the non-uniform case. In this setup, we show the advantages of our strategy in the low computational budget regime. While the proposed solution is harder to implement in hardware, we believe it sets a basis for new alternatives to neural networks quantization.


title: The method of decoupling fields generalized to higher spatial derivatives

abstract: This work studies the spatial derivatives of decoupling fields to strongly coupled forward-backward stochastic differential equations in a Brownian setting. We formally deduce the backward dynamics of the first and higher spatial derivatives. In addition, we study necessary conditions under which singularities in either one of them can occur while moving backwards in time.


title: Some remarks on the volume of log varieties

abstract: In this note, using methods introduced by Hacon, McKernan and Xu, we study the accumulation points of volumes of varieties of log general type. First, we show that, if the set of boundary coefficients $\Lambda$ is DCC, closed under limits and contains 1, then also the corresponding set of volumes is DCC and closed under limits. Then, we consider the case of $\epsilon$-log canonical varieties, for $0 < \epsilon < 1$. In this situation, we prove that, if $\Lambda$ is finite, then the corresponding set of volumes is discrete.


title: An adaptive local discrete convolution method for the numerical solution of Maxwell's equations

abstract: We present a numerical method for solving the free-space Maxwell's equations in three dimensions using compact convolution kernels on a rectangular grid. We first rewrite Maxwell's Equations as a system of wave equations with auxiliary variables and discretize its solution from the method of spherical means. The algorithm has been extended to be used on a locally-refined nested hierarchy of rectangular grids.


title: A Basic Result on the Superposition of Arrival Processes in Deterministic Networks

abstract: Time-Sensitive Networking (TSN) and Deterministic Networking (DetNet) are emerging standards to enable deterministic, delay-critical communication in such networks. This naturally (re-)calls attention to the network calculus theory (NC), since a rich set of results for delay guarantee analysis have already been developed there. One could anticipate an immediate adoption of those existing network calculus results to TSN and DetNet. However, the fundamental difference between the traffic specification adopted in TSN and DetNet and those traffic models in NC makes this difficult, let alone that there is a long-standing open challenge in NC. To address them, this paper considers an arrival time function based max-plus NC traffic model. In particular, for the former, the mapping between the TSN / DetNet and the NC traffic model is proved. For the latter, the superposition property of the arrival time function based NC traffic model is found and proved. Appealingly, the proved superposition property shows a clear analogy with that of a well-known counterpart traffic model in NC. These results help make an important step towards the development of a system theory for delay guarantee analysis of TSN / DetNet networks.


title: From Credit Assignment to Entropy Regularization: Two New Algorithms for Neural Sequence Prediction

abstract: In this work, we study the credit assignment problem in reward augmented maximum likelihood (RAML) learning, and establish a theoretical equivalence between the token-level counterpart of RAML and the entropy regularized reinforcement learning. Inspired by the connection, we propose two sequence prediction algorithms, one extending RAML with fine-grained credit assignment and the other improving Actor-Critic with a systematic entropy regularization. On two benchmark datasets, we show the proposed algorithms outperform RAML and Actor-Critic respectively, providing new alternatives to sequence prediction.


title: Matryoshka Networks: Predicting 3D Geometry via Nested Shape Layers

abstract: In this paper, we develop novel, efficient 2D encodings for 3D geometry, which enable reconstructing full 3D shapes from a single image at high resolution. The key idea is to pose 3D shape reconstruction as a 2D prediction problem. To that end, we first develop a simple baseline network that predicts entire voxel tubes at each pixel of a reference view. By leveraging well-proven architectures for 2D pixel-prediction tasks, we attain state-of-the-art results, clearly outperforming purely voxel-based approaches. We scale this baseline to higher resolutions by proposing a memory-efficient shape encoding, which recursively decomposes a 3D shape into nested shape layers, similar to the pieces of a Matryoshka doll. This allows reconstructing highly detailed shapes with complex topology, as demonstrated in extensive experiments; we clearly outperform previous octree-based approaches despite having a much simpler architecture using standard network components. Our Matryoshka networks further enable reconstructing shapes from IDs or shape similarity, as well as shape sampling.


title: Connection formulas between Coulomb wave functions

abstract: The mathematical relations between the regular Coulomb function $F_{\eta\ell}(\rho)$ and the irregular Coulomb functions $H^\pm_{\eta\ell}(\rho)$ and $G_{\eta\ell}(\rho)$ are obtained in the complex plane of the variables $\eta$ and $\rho$ for integer or half-integer values of $\ell$. These relations, referred to as "connection formulas", form the basis of the theory of Coulomb wave functions, and play an important role in many fields of physics, especially in the quantum theory of charged particle scattering. As a first step, the symmetry properties of the regular function $F_{\eta\ell}(\rho)$ are studied, in particular under the transformation $\ell\mapsto-\ell-1$, by means of the modified Coulomb function $\Phi_{\eta\ell}(\rho)$, which is entire in the dimensionless energy $\eta^{-2}$ and the angular momentum $\ell$. Then, it is shown that, for integer or half-integer $\ell$, the irregular functions $H^\pm_{\eta\ell}(\rho)$ and $G_{\eta\ell}(\rho)$ can be expressed in terms of the derivatives of $\Phi_{\eta,\ell}(\rho)$ and $\Phi_{\eta,-\ell-1}(\rho)$ with respect to $\ell$. As a consequence, the connection formulas directly lead to the description of the singular structures of $H^\pm_{\eta\ell}(\rho)$ and $G_{\eta\ell}(\rho)$ at complex energies in their whole Riemann surface. The analysis of the functions is supplemented by novel graphical representations in the complex plane of $\eta^{-1}$.


title: Block Sparse Multi-lead ECG Compression Exploiting between-lead Collaboration

abstract: Multilead ECG compression (MlEC) has attracted tremendous attention in long-term monitoring of the patients heart behavior. This paper proposes a method denoted by block sparse MlEC (BlS MlEC) in order to exploit between-lead correlations to compress the signals in a more efficient way. This is due to the fact that multi-lead ECG signals are multiple observations of the same source (heart) from different locations. Consequently, they have high correlation in terms of the support set of their sparse models which leads them to share dominant common structure. In order to obtain the block sparse model, the collaborative version of lasso estimator is applied. In addition, we have shown that raised cosine kernel has advantages over conventional Gaussian and wavelet (Daubechies family) due to its specific properties. It is demonstrated that using raised cosine kernel in constructing the sparsifying basis matrix gives a sparser model which results in higher compression ratio and lower reconstruction error. The simulation results show the average improvement of 37%, 88% and 90-97% for BlS M-lEC compared to the non-collaborative case with raised cosine kernel, Gaussian kernel and collaborative case with Daubechies wavelet kernels, respectively, in terms of reconstruction error while the compression ratio is considered fixed.


title: Optimization based evaluation of grating interferometric phase stepping series and analysis of mechanical setup instabilities

abstract: The diffraction contrast modalities accessible by X-ray grating interferometers are not imaged directly but have to be inferred from sine like signal variations occurring in a series of images acquired at varying relative positions of the interferometer's gratings. The absolute spatial translations involved in the acquisition of these phase stepping series usually lie in the range of only a few hundred nanometers, wherefore positioning errors as small as 10nm will already translate into signal uncertainties of one to ten percent in the final images if not accounted for. Classically, the relative grating positions in the phase stepping series are considered input parameters to the analysis and are, for the Fast Fourier Transform that is typically employed, required to be equidistantly distributed over multiples of the gratings' period. More general optimization based approaches with relaxed requirements on the sampling positions are often deemed too time consuming in contrast. In the following, a fast converging optimization scheme is presented simultaneously determining the phase stepping curves' parameters as well as the actually performed motions of the stepped grating, including also erroneous rotational motions which are commonly neglected. While the correction of solely the translational errors along the stepping direction is found to be sufficient with regard to the reduction of image artifacts, the possibility to also detect minute rotations about all axes proves to be a valuable tool for system calibration and monitoring. The simplicity of the provided algorithm, in particular when only considering translational errors, makes it well suitable as a standard evaluation procedure also for large image series.


title: Memory effect, conformal symmetry and gravitational plane waves

abstract: We discuss in some detail the interaction of classical particles, including the scattering and memory effect, with a pulse of gravitational plane wave. The key point is the conformal symmetry of gravitational plane waves. In particular, we obtain, in the limit of short pulse, some results for impulsive gravitational waves. Furthermore, in the general case, we give certain conditions which allow us to completely describe the interaction in terms of the singular Baldwin-Jeffery-Rosen coordinates.


title: Neural-networks-based Photon-Counting Data Correction: Pulse Pileup Effect

abstract: Compared with the start-of-art energy integration detectors (EIDs), photon-counting detectors (PCDs) with energy discrimination capabilities have demonstrated great potentials in various applications of medical x-ray radiography and computed tomography (CT). However, the advantages of PCDs may be compromised by various hardware-related degradation factors. For example, at high count rate, quasi-coincident photons may be piled up and result in not only lost counts but also distorted spectrums, which is called the pulse pileup effects. Considering the relative slow detection speed of current PCDs and high x-ray flux for clinical imaging, it is vital to develop an effective approach to compensate or correct for the pileup effect. The aim of this paper was to develop a model-free, end-to-end, and data-driven approach in the neural network framework. We first introduce the trigger threshold concept to regularize the correction process. Then, we design a fully-connected cascade artificial neural network and train it with measurements and true counts. After training, the network is used to recover raw photon-counting data for higher fidelity. To evaluate the proposed approach, Monte Carlo simulations are conducted. We simulate four types of PCDs combining the unipolar or bipolar pulse shape with the nonparalyzable or paralyzable detection mode. The simulation results demonstrate that the trained neural network decreases the relative error of count averages of raw data by an order of magnitude while maintaining a reasonable variance over a wide range of attenuation paths.


title: On the Enumeration of Maximal $(\Delta, \gamma)$-Cliques of a Temporal Network

abstract: A temporal network is a mathematical way of precisely representing a time varying relationship among a group of agents. In this paper, we introduce the notion of $(\Delta, \gamma)$-Cliques of a temporal network, where every pair of vertices present in the clique communicates atleast $\gamma$ times in each $\Delta$ period within a given time duration. We present an algorithm for enumerating all such maximal cliques present in the network. We also implement the proposed algorithm with three human contact network data sets. Based on the obtained results, we analyze the data set on multiple values of $\Delta$ and $\gamma$, which helps in finding out contact groups with different frequencies.


title: Of Wines and Reviews: Measuring and Modeling the Vivino Wine Social Network

abstract: This paper presents an analysis of social experiences around wine consumption through the lens of Vivino, a social network for wine enthusiasts with over 26 million users worldwide. We compare users' perceptions of various wine types and regional styles across both New and Old World wines, examining them across price ranges, vintages, regions, varietals, and blends. Among other things, we find that ratings provided by Vivino users are not biased by cost. We then study how wine characteristics, language in wine reviews, and the distribution of wine ratings can be combined to develop prediction models. More specifically, we model user behavior to develop a regression model for predicting wine ratings, and a classifier for determining user review preferences.


title: Universal, high-fidelity quantum gates based on superadiabatic, geometric phases on a solid-state spin-qubit at room temperature

abstract: Geometric phases and holonomies (their non-commuting generalizations) are a promising resource for the realization of high-fidelity quantum operations in noisy devices, due to their intrinsic fault-tolerance against noise and experimental imperfections. Despite their conceptual appeal and proven fault-tolerance, for a long time their practical use in quantum computing was limited to proof of principle demonstrations. Only in 2012 Sj\"oqvist et al. formulated a strategy to generate non-Abelian (i.e. holonomic) quantum gates through non-adiabatic transformation. Successful experimental demonstrations of this concept followed on various physical qubit systems and proved the feasibility of this fast, holonomic quantum gate concept. Despite these successes, the experimental implementation of such non-Abelian quantum gates remains experimentally challenging since in general the emergence of a suitable holonomy requires encoding of the logical qubit within a three (or higher) level system being driven by two (or more) control fields. A very recent proposal by Liang et al. offers an elegant solution generating a non-Abelian, geometric quantum gate on a simple, two-level system driven by one control field. Exploiting the concept of transitionless quantum driving it allows the generation of universal geometric quantum gates through superadiabatic evolution. This concept thus generates fast and robust phase-based quantum gates on the basis of minimal experimental resources. Here, we report on the first such implementation of a set of non-commuting single-qubit superadiabatic geometric quantum gates on the electron spin of the negatively charged nitrogen vacancy center in diamond. The realized quantum gates combine high-fidelity and fast quantum gate performance. This provides a promising and powerful tool for large-scale quantum computing under realistic, noisy experimental conditions.


title: Riesz bases from orthonormal bases by replacement

abstract: Given an orthonormal basis $ {\mathcal V}= \{v_j\} _{j\in N}$ in a separable Hilbert space $H$ and a set of unit vectors $ {\mathcal B}=\{w_j\}_{j\in N}$, we consider the sets $ {\mathcal B}_N$ obtained by replacing the vectors $v_1, ...,\, v_N$ with vectors $w_1,\, ...,\, w_N$. We show necessary and sufficient conditions that ensure that the sets $ {\mathcal B}_N$ are Riesz bases of $H$ and we estimate the frame constants of the $ {\mathcal B}_N$. Then, we prove conditions that ensure that $ {\mathcal B}$ is a Riesz basis. Applications to the construction of exponential bases on domains of $ R^d$ are also presented.


title: Efficient Algorithms for Asymptotic Bounds on Termination Time in VASS

abstract: Vector Addition Systems with States (VASS) provide a well-known and fundamental model for the analysis of concurrent processes, parameterized systems, and are also used as abstract models of programs in resource bound analysis. In this paper we study the problem of obtaining asymptotic bounds on the termination time of a given VASS. In particular, we focus on the practically important case of obtaining polynomial bounds on termination time. Our main contributions are as follows: First, we present a polynomial-time algorithm for deciding whether a given VASS has a linear asymptotic complexity. We also show that if the complexity of a VASS is not linear, it is at least quadratic. Second, we classify VASS according to quantitative properties of their cycles. We show that certain singularities in these properties are the key reason for non-polynomial asymptotic complexity of VASS. In absence of singularities, we show that the asymptotic complexity is always polynomial and of the form $\Theta(n^k)$, for some integer $k\leq d$, where $d$ is the dimension of the VASS. We present a polynomial-time algorithm computing the optimal $k$. For general VASS, the same algorithm, which is based on a complete technique for the construction of ranking functions in VASS, produces a valid lower bound, i.e., a $k$ such that the termination complexity is $\Omega(n^k)$. Our results are based on new insights into the geometry of VASS dynamics, which hold the potential for further applicability to VASS analysis.


title: Sparse Grid Approximation Spaces for Space-Time Boundary Integral Formulations of the Heat Equation

abstract: The aim of this paper is to develop stable and accurate numerical schemes for boundary integral formulations of the heat equation with Dirichlet boundary conditions. The accuracy of Galerkin discretisations for the resulting boundary integral formulations depends mainly on the choice of discretisation space. We develop a-priori error analysis utilising a proof technique that involves norm equivalences in hierarchical wavelet subspace decompositions. We apply this to a full tensor product discretisation, showing improvements over existing results, particularly for discretisation spaces having low polynomial degrees. We then use the norm equivalences to show that an anisotropic sparse grid discretisation yields even higher convergence rates. Finally, a simple adaptive scheme is proposed to suggest an optimal shape for the sparse grid index sets.


title: Feedforward Architectures for Decentralized Precoding in Massive MU-MIMO Systems

abstract: Massive multi-user multiple-input multiple-output (MU-MIMO) enables significant gains in spectral efficiency and link reliability compared to conventional small-scale MIMO technology. Furthermore, linear precoders, e.g., using zero forcing or Wiener filter (WF) precoding, are sufficient to achieve excellent error-rate performance in the massive MU-MIMO downlink. However, these methods necessitate centralized processing at the base-station (BS), which causes (i) excessively high interconnect and chip input/output data rates, and (ii) high implementation complexity. We propose two feedforward architectures and corresponding decentralized WF precoders that parallelize precoding across multiple computing fabrics, effectively mitigating the issues of centralized approaches. To demonstrate the efficacy of our decentralized precoders, we provide implementation results on a multi-GPU system, which show that our solutions achieve throughputs in the Gbit/s regime while achieving (near-)optimal error-rate performance in the massive MU-MIMO downlink.


title: SHADE: Information Based Regularization for Deep Learning

abstract: Regularization is a big issue for training deep neural networks. In this paper, we propose a new information-theory-based regularization scheme named SHADE for SHAnnon DEcay. The originality of the approach is to define a prior based on conditional entropy, which explicitly decouples the learning of invariant representations in the regularizer and the learning of correlations between inputs and labels in the data fitting term. Our second contribution is to derive a stochastic version of the regularizer compatible with deep learning, resulting in a tractable training scheme. We empirically validate the efficiency of our approach to improve classification performances compared to common regularization schemes on several standard architectures.


title: Mapping the emission line strengths and kinematics of Supernova Remnant S147 with extensive LAMOST spectroscopic observations

abstract: We present extensive spectroscopic observations of supernova remnant (SNR) S147 collected with the Large sky Area Multi-Object fiber Spectroscopic Telescope (LAMOST). The spectra were carefully sky-subtracted taking into account the complex filamentary structure of S147. We have utilized all available LAMOST spectra toward S147, including sky and stellar spectra. By measuring the prominent optical emission lines including H{\alpha}, [NII]{\lambda}6584, and [SII]{\lambda}{\lambda}6717,6731, we present maps of radial velocity and line intensity ratio covering the whole nebula of S147 with unprecedented detail. The maps spatially correlated well with the complex filamentary structure of S147. For the central 2 deg of S147, the radial velocity varies from -100 to 100 km/s and peaks between ~ 0 and 10 km/s. The intensity ratios of H{\alpha}/[SII]{\lambda}{\lambda}6717,6731, [SII]{\lambda}6717/{\lambda}6731 and H{\alpha}/[NII]{\lambda}6584 peak at about 0.77, 1.35 and 1.48, respectively, with a scatter of 0.17, 0.19 and 0.37, respectively. The intensity ratios are consistent with the literature values. However, the range of variations of line intensity ratios estimated here and representative of the whole nebula, are larger than previously estimated.


title: On Obtaining Stable Rankings

abstract: Decision making is challenging when there is more than one criterion to consider. In such cases, it is common to assign a goodness score to each item as a weighted sum of its attribute values and rank them accordingly. Clearly, the ranking obtained depends on the weights used for this summation. Ideally, one would want the ranked order not to change if the weights are changed slightly. We call this property {\em stability} of the ranking. A consumer of a ranked list may trust the ranking more if it has high stability. A producer of a ranked list prefers to choose weights that result in a stable ranking, both to earn the trust of potential consumers and because a stable ranking is intrinsically likely to be more meaningful. In this paper, we develop a framework that can be used to assess the stability of a provided ranking and to obtain a stable ranking within an "acceptable" range of weight values (called "the region of interest"). We address the case where the user cares about the rank order of the entire set of items, and also the case where the user cares only about the top-$k$ items. Using a geometric interpretation, we propose algorithms that produce stable rankings. In addition to theoretical analyses, we conduct extensive experiments on real datasets that validate our proposal.


title: Effective throughput of MISO systems over $\kappa-\mu$ shadowed fading channels: MGF based analysis

abstract: The effective throughput of multiple-input single-output (MISO) wireless communication systems over $\kappa-\mu$ shadowed fading channels is analysed. To obtain exact closed-from expression, the moment generating function (MGF) of the instantaneous signal to noise ratio (SNR) with independent and identically distributed (i.i.d.) transmit antennas is employed. In addition, a special case of integer values for the fading parameter is performed to simplify the derived expression. The asymptotic behaviour of the effective throughput at high values of SNR is also studied. The Monte Carlo simulations and the numerical results are presented to verify the validation of our analysis.


title: Semi-parametric Image Synthesis

abstract: We present a semi-parametric approach to photographic image synthesis from semantic layouts. The approach combines the complementary strengths of parametric and nonparametric techniques. The nonparametric component is a memory bank of image segments constructed from a training set of images. Given a novel semantic layout at test time, the memory bank is used to retrieve photographic references that are provided as source material to a deep network. The synthesis is performed by a deep network that draws on the provided photographic material. Experiments on multiple semantic segmentation datasets show that the presented approach yields considerably more realistic images than recent purely parametric techniques. The results are shown in the supplementary video at https://youtu.be/U4Q98lenGLQ


title: On the Iwasawa main conjectures for modular forms at non-ordinary primes

abstract: In this paper, we prove under mild hypotheses the Iwasawa main conjectures of Lei--Loeffler--Zerbes for modular forms of weight $2$ at non-ordinary primes. Our proof is based on the study of the two-variable analogues of these conjectures formulated by B\"uy\"ukboduk--Lei for imaginary quadratic fields in which $p$ splits, and on anticyclotomic Iwasawa theory. As application of our results, we deduce the $p$-part of the Birch and Swinnerton-Dyer formula in analytic ranks $0$ or $1$ for abelian varieties over $\mathbb{Q}$ of ${\rm GL}_2$-type for non-ordinary primes $p>2$.


title: Transmission Capacity of Full-Duplex MIMO Ad-Hoc Network with Limited Self-Interference Cancellation

abstract: In this paper, we propose a joint transceiver beamforming design to simultaneously mitigate self-interference (SI) and partial inter-node interference for full-duplex multiple-input and multiple-output ad-hoc network, and then derive the transmission capacity upper bound (TC-UB) for the corresponding network. Condition on a specified transceiver antenna's configuration, we allow the SI effect to be cancelled at transmitter side, and offer an additional degree-of-freedom at receiver side for more inter-node interference cancellation. In addition, due to the proposed beamforming design and imperfect SI channel estimation, the conventional method to obtain the TC-UB is not applicable. This motivates us to exploit the dominating interferer region plus Newton-Raphson method to iteratively formulate the TC-UB. The results show that the derived TC-UB is quite close to the actual one especially when the number of receive-antenna is small. Moreover, our proposed beamforming design outperforms the existing beamforming strategies, and FD mode works better than HD mode in low signal-to-noise ratio region.


title: Perspectives on Nuclear Structure and Scattering with the Ab Initio No-Core Shell Model

abstract: Nuclear structure and reaction theory are undergoing a major renaissance with advances in many-body methods, strong interactions with greatly improved links to Quantum Chromodynamics (QCD), the advent of high performance computing, and improved computational algorithms. Predictive power, with well-quantified uncertainty, is emerging from non-perturbative approaches along with the potential for new discoveries such as predicting nuclear phenomena before they are measured. We present an overview of some recent developments and discuss challenges that lie ahead. Our focus is on explorations of alternative truncation schemes in the harmonic oscillator basis, of which our Japanese--United States collaborative work on the No-Core Monte-Carlo Shell Model is an example. Collaborations with Professor Takaharu Otsuka and his group have been instrumental in these developments.


title: Probing spin-phonon interactions in silicon carbide with Gaussian acoustics

abstract: Hybrid spin-mechanical systems provide a platform for integrating quantum registers and transducers. Efficient creation and control of such systems require a comprehensive understanding of the individual spin and mechanical components as well as their mutual interactions. Point defects in silicon carbide (SiC) offer long-lived, optically addressable spin registers in a wafer-scale material with low acoustic losses, making them natural candidates for integration with high quality factor mechanical resonators. Here, we show Gaussian focusing of a surface acoustic wave in SiC, characterized by a novel stroboscopic X-ray diffraction imaging technique, which delivers direct, strain amplitude information at nanoscale spatial resolution. Using ab initio calculations, we provide a more complete picture of spin-strain coupling for various defects in SiC with C3v symmetry. This reveals the importance of shear for future device engineering and enhanced spin-mechanical coupling. We demonstrate all-optical detection of acoustic paramagnetic resonance without microwave magnetic fields, relevant to sensing applications. Finally, we show mechanically driven Autler-Townes splittings and magnetically forbidden Rabi oscillations. These results offer a basis for full strain control of three-level spin systems.


title: Large Piezoelectric Response of van der Waals Layered Solids

abstract: The bulk piezoelectric response, as measured by the piezoelectric modulus tensor (\textbf{d}), is determined by a combination of charge redistribution due to strain and the amount of strain produced by the application of stress (stiffness). Motivated by the notion that less stiff materials could exhibit large piezoelectric responses, herein we investigate the piezoelectric modulus of van der Waals-bonded quasi-2D ionic compounds using first-principles calculations. From a pool of 869 known binary and ternary quasi-2D materials, we have identified 135 non-centrosymmetric crystals of which 48 systems are found to have \textbf{d} components larger than the longitudinal piezoelectric modulus of AlN (a common piezoelectric for resonators), and three systems with the response greater than that of PbTiO$_3$, which is among the materials with largest known piezoelectric modulus. None of the identified materials have previously been considered for piezoelectric applications. Furthermore, we find that large \textbf{d} components always couple to the deformations (shearing or axial) of van der Waals "gaps" between the layers and are indeed enabled by the weak intra-layer interactions.


title: Optimal Scheduling for Exposed Datapath Architectures with Buffered Processing Units by ASP

abstract: Conventional processor architectures are restricted in exploiting instruction level parallelism (ILP) due to the relatively low number of programmer-visible registers. Therefore, more recent processor architectures expose their datapaths so that the compiler (1) can schedule parallel instructions to different processing units and (2) can make effective use of local storage of the processing units. Among these architectures, the Synchronous Control Asynchronous Dataflow (SCAD) architecture is a new exposed datapath architecture whose processing units are equipped with first-in first-out (FIFO) buffers at their input and output ports. In contrast to register-based machines, the optimal code generation for SCAD is still a matter of research. In particular, SAT and SMT solvers were used to generate optimal resource constrained and optimal time constrained schedules for SCAD, respectively. As Answer Set Programming (ASP) offers better flexibility in handling such scheduling problems, we focus in this paper on using an answer set solver for both resource and time constrained optimal SCAD code generation. As a major benefit of using ASP, we are able to generate \emph{all} optimal schedules for a given program which allows one to study their properties. Furthermore, the experimental results of this paper demonstrate that the answer set solver can compete with SAT solvers and outperforms SMT solvers. \emph{This paper is under consideration for acceptance in TPLP.}


title: But Who Protects the Moderators? The Case of Crowdsourced Image Moderation

abstract: Though detection systems have been developed to identify obscene content such as pornography and violence, artificial intelligence is simply not good enough to fully automate this task yet. Due to the need for manual verification, social media companies may hire internal reviewers, contract specialized workers from third parties, or outsource to online labor markets for the purpose of commercial content moderation. These content moderators are often fully exposed to extreme content and may suffer lasting psychological and emotional damage. In this work, we aim to alleviate this problem by investigating the following question: How can we reveal the minimum amount of information to a human reviewer such that an objectionable image can still be correctly identified? We design and conduct experiments in which blurred graphic and non-graphic images are filtered by human moderators on Amazon Mechanical Turk (AMT). We observe how obfuscation affects the moderation experience with respect to image classification accuracy, interface usability, and worker emotional well-being.


title: Zolotarev Iterations for the Matrix Square Root

abstract: We construct a family of iterations for computing the principal square root of a square matrix $A$ using Zolotarev's rational minimax approximants of the square root function. We show that these rational functions obey a recursion, allowing one to iteratively generate optimal rational approximants of $\sqrt{z}$ of high degree using compositions and products of low-degree rational functions. The corresponding iterations for the matrix square root converge to $A^{1/2}$ for any input matrix $A$ having no nonpositive real eigenvalues. In special limiting cases, these iterations reduce to known iterations for the matrix square root: the lowest-order version is an optimally scaled Newton iteration, and for certain parameter choices, the principal family of Pad\'e iterations is recovered. Theoretical results and numerical experiments indicate that the iterations perform especially well on matrices having eigenvalues with widely varying magnitudes.


title: A Stochastic Model for UAV Networks Positioned Above Demand Hotspots in Urban Environments

abstract: Wireless access points on unmanned aerial vehicles (UAVs) are being considered for mobile service provisioning in commercial networks. To be able to efficiently use these devices in cellular networks it is necessary to first have a qualitative and quantitative understanding of how their design parameters reflect on the service quality experienced by the end user. In this paper we use stochastic geometry to characterise the behaviour of a network of UAV access points that intelligently position themselves above user hotspots, and we evaluate the performance of such a network against cases where the UAVs are positioned in a rectangular grid or according to heuristic positioning algorithms.


title: Precision Medicine as an Accelerator for Next Generation Cognitive Supercomputing

abstract: In the past several years, we have taken advantage of a number of opportunities to advance the intersection of next generation high-performance computing AI and big data technologies through partnerships in precision medicine. Today we are in the throes of piecing together what is likely the most unique convergence of medical data and computer technologies. But more deeply, we observe that the traditional paradigm of computer simulation and prediction needs fundamental revision. This is the time for a number of reasons. We will review what the drivers are, why now, how this has been approached over the past several years, and where we are heading.


title: Gradient Sampling Methods for Nonsmooth Optimization

abstract: This paper reviews the gradient sampling methodology for solving nonsmooth, nonconvex optimization problems. An intuitively straightforward gradient sampling algorithm is stated and its convergence properties are summarized. Throughout this discussion, we emphasize the simplicity of gradient sampling as an extension of the steepest descent method for minimizing smooth objectives. We then provide overviews of various enhancements that have been proposed to improve practical performance, as well as of several extensions that have been made in the literature, such as to solve constrained problems. The paper also includes clarification of certain technical aspects of the analysis of gradient sampling algorithms, most notably related to the assumptions one needs to make about the set of points at which the objective is continuously differentiable. Finally, we discuss possible future research directions.


title: Graphene Based Adaptive Thermal Camouflage

abstract: In nature adaptive coloration has been effectively utilized for concealment and signaling. Various biological mechanisms have evolved that can tune the reflectivity for visible and ultraviolet light. These examples inspire many artificial systems for mimicking adaptive coloration to match the visual appearance to their surroundings. Thermal camouflage, however, has been an outstanding challenge which requires an ability to control of the emitted thermal radiation from the surface. Here we report a new class of active thermal surfaces capable of efficient real-time electrical-control of thermal emission over the full infrared (IR) spectrum without changing the temperature of the surface. Our approach relies on electro-modulation of IR absorptivity and emissivity of multilayer graphene via reversible intercalation of nonvolatile ionic liquids. The demonstrated devices are light (30 g/m2), thin (<50 um) and ultra-flexible which can conformably coat their environment. In addition, combining active thermal surfaces with a feedback mechanism, we demonstrate realization of an adaptive thermal camouflage system which can reconfigure its thermal appearance and blend itself with the thermal background in a few seconds. Furthermore, we show that these devices can disguise hot objects as cold and cold ones as hot in a thermal imaging system. We anticipate that, the electrical control of thermal radiation will impact on a variety of new technologies ranging from adaptive IR optics to heat management for outer space applications.


title: Building Models for Biopathway Dynamics Using Intrinsic Dimensionality Analysis

abstract: An important task for many if not all the scientific domains is efficient knowledge integration, testing and codification. It is often solved with model construction in a controllable computational environment. In spite of that, the throughput of in-silico simulation-based observations become similarly intractable for thorough analysis. This is especially the case in molecular biology, which served as a subject for this study. In this project, we aimed to test some approaches developed to deal with the curse of dimensionality. Among these we found dimension reduction techniques especially appealing. They can be used to identify irrelevant variability and help to understand critical processes underlying high-dimensional datasets. Additionally, we subjected our data sets to nonlinear time series analysis, as those are well established methods for results comparison. To investigate the usefulness of dimension reduction methods, we decided to base our study on a concrete sample set. The example was taken from the domain of systems biology concerning dynamic evolution of sub-cellular signaling. Particularly, the dataset relates to the yeast pheromone pathway and is studied in-silico with a stochastic model. The model reconstructs signal propagation stimulated by a mating pheromone. In the paper, we elaborate on the reason of multidimensional analysis problem in the context of molecular signaling, and next, we introduce the model of choice, simulation details and obtained time series dynamics. A description of used methods followed by a discussion of results and their biological interpretation finalize the paper.

categories: q-bio.MN cs.IT math.IT q-bio.QM stat.AP stat.ML


title: The detectability of radio emission from exoplanets

abstract: Like the magnetised planets in our Solar System, magnetised exoplanets should emit strongly at radio wavelengths. Radio emission directly traces the planetary magnetic fields and radio detections can place constraints on the physical parameters of these features. Large comparative studies of predicted radio emission characteristics for the known population of exoplanets help to identify what physical parameters could be key for producing bright, observable radio emission. Since the last comparative study, many thousands of exoplanets have been discovered. We report new estimates for the radio flux densities and maximum emission frequencies for the current population of known exoplanets orbiting pre-main sequence and main-sequence stars with spectral types F-M. The set of exoplanets predicted to produce observable radio emission are Hot Jupiters orbiting young stars. The youth of these system predicts strong stellar magnetic fields and/or dense winds, which are key for producing bright, observable radio emission. We use a new all-sky circular polarisation Murchison Widefield Array survey to place sensitive limits on 200 MHz emission from exoplanets, with $3\sigma$ values ranging from 4.0 - 45.0 mJy. Using a targeted Giant Metre Wave Radio Telescope observing campaign, we also report a $3\sigma$ upper limit of 4.5 mJy on the radio emission from V830 Tau b, the first Hot Jupiter to be discovered orbiting a pre-main sequence star. Our limit is the first to be reported for the low-frequency radio emission from this source.


title: Triangle Inscribed-Triangle Picking

abstract: Given a triangle ABC, we derive the probability distribution function and the moments of the area of an inscribed triangle RST whose vertices are uniformly distributed on AB, BC, and CA. The theoretical results are confirmed by a Monte Carlo simulation.


title: Vertex-minors and the Erd\H{o}s-Hajnal conjecture

abstract: We prove that for every graph $H$, there exists $\varepsilon>0$ such that every $n$-vertex graph with no vertex-minors isomorphic to $H$ has a pair of disjoint sets $A$, $B$ of vertices such that $|A|, |B|\ge \varepsilon n$ and $A$ is complete or anticomplete to $B$. We deduce this from recent work of Chudnovsky, Scott, Seymour, and Spirkl (2018). This proves the analog of the Erd\H{o}s-Hajnal conjecture for vertex-minors.


title: A Compendium of Hopf-Like Bifurcations in Piecewise-Smooth Dynamical Systems

abstract: This Letter outlines 20 geometric mechanisms by which limit cycles are created locally in two-dimensional piecewise-smooth systems of ODEs. These include boundary equilibrium bifurcations of hybrid systems, Filippov systems, and continuous systems, and limit cycles created from folds and by the addition of hysteresis or time-delay. Scaling laws for the amplitude and period of the limit cycles are compared to (classical) Hopf bifurcations.


title: A Matrix Gaussian Distribution

abstract: In this note, we define a Gaussian probability distribution over matrices. We prove some useful properties of this distribution, namely, the fact that marginalization, conditioning, and affine transformations preserve the matrix Gaussian distribution. We also derive useful results regarding the expected value of certain quadratic forms based solely on covariances between matrices. Previous definitions of matrix normal distributions are severely under-parameterized, assuming unrealistic structure on the covariance (see Section 2). We believe that our generalization is better equipped for use in practice.


title: Joint Analysis of Individual-level and Summary-level GWAS Data by Leveraging Pleiotropy

abstract: A large number of recent genome-wide association studies (GWASs) for complex phenotypes confirm the early conjecture for polygenicity, suggesting the presence of large number of variants with only tiny or moderate effects. However, due to the limited sample size of a single GWAS, many associated genetic variants are too weak to achieve the genome-wide significance. These undiscovered variants further limit the prediction capability of GWAS. Restricted access to the individual-level data and the increasing availability of the published GWAS results motivate the development of methods integrating both the individual-level and summary-level data. How to build the connection between the individual-level and summary-level data determines the efficiency of using the existing abundant summary-level resources with limited individual-level data, and this issue inspires more efforts in the existing area. In this study, we propose a novel statistical approach, LEP, which provides a novel way of modeling the connection between the individual-level data and summary-level data. LEP integrates both types of data by \underline{LE}veraing \underline{P}leiotropy to increase the statistical power of risk variants identification and the accuracy of risk prediction. The algorithm for parameter estimation is developed to handle genome-wide-scale data. Through comprehensive simulation studies, we demonstrated the advantages of LEP over the existing methods. We further applied LEP to perform integrative analysis of Crohn's disease from WTCCC and summary statistics from GWAS of some other diseases, such as Type 1 diabetes, Ulcerative colitis and Primary biliary cirrhosis. LEP was able to significantly increase the statistical power of identifying risk variants and improve the risk prediction accuracy from 63.39\% ($\pm$ 0.58\%) to 68.33\% ($\pm$ 0.32\%) using about 195,000 variants.


title: Stable and contact-free time stepping for dense rigid particle suspensions

abstract: We consider suspensions of rigid bodies in a two-dimensional viscous fluid. Even with high-fidelity numerical methods, unphysical contact between particles occurs because of spatial and temporal discretization errors. We apply the method of Lu et al. [Journal of Computational Physics, 347:160-182, 2017] where overlap is avoided by imposing a minimum separation distance. In its original form, the method discretizes interactions between different particles explicitly. Therefore, to avoid stiffness, a large minimum separation distance is used. In this paper, we extend the method of Lu et al. by treating all interactions implicitly. This new time stepping method is able to simulate dense suspensions with large time step sizes and a small minimum separation distance. The method is tested on various unbounded and bounded flows, and rheological properties of the resulting suspensions are computed.


title: Cycle-Consistent Deep Generative Hashing for Cross-Modal Retrieval

abstract: In this paper, we propose a novel deep generative approach to cross-modal retrieval to learn hash functions in the absence of paired training samples through the cycle consistency loss. Our proposed approach employs adversarial training scheme to lean a couple of hash functions enabling translation between modalities while assuming the underlying semantic relationship. To induce the hash codes with semantics to the input-output pair, cycle consistency loss is further proposed upon the adversarial training to strengthen the correlations between inputs and corresponding outputs. Our approach is generative to learn hash functions such that the learned hash codes can maximally correlate each input-output correspondence, meanwhile can also regenerate the inputs so as to minimize the information loss. The learning to hash embedding is thus performed to jointly optimize the parameters of the hash functions across modalities as well as the associated generative models. Extensive experiments on a variety of large-scale cross-modal data sets demonstrate that our proposed method achieves better retrieval results than the state-of-the-arts.


title: A stellar fly-by close to the Galactic centre: Can we detect stars on highly-relativistic orbits?

abstract: The Galactic centre Nuclear Star Cluster is one of the densest stellar clusters in the Galaxy. The stars in its inner portions orbit the supermassive black hole associated with compact radio source Sgr~A* at the orbital speeds of several thousand km/s. The B-type star S2 is currently the best case to test the general relativity as well as other theories of gravity based on its stellar orbit. Yet its orbital period of $\sim 16\,{\rm yr}$ and the eccentricity of $\sim 0.88$ yields the relativistic pericentre shift of $\sim 11'$, which is observationally still difficult to reliably measure due to possible Newtonian perturbations as well as reference-frame uncertainties. A naive way to solve this problem is to find stars with smaller pericentre distances, $r_{\rm p}\lesssim 1529$ Schwarzschild radii ($120\,{\rm AU}$), and thus more prominent relativistic effects. In this contribution, we show that to detect stars on relativistic orbits is progressively less likely given the volume shrinkage and the expected stellar density distributions. Finally, one arrives to a sparse region, where the total number of bright stars is expected to fall below one. One can, however, still potentially detect stars crossing this region. In this contribution, we provide a simple formula for the detection probability of a star crossing a sparse region. We also examine an approximate time-scale on which the star reappears in the sparse region, i.e. a "waiting" time-scale for observers.


title: Effective Bounds on the Dimensions of Jacobians Covering Abelian Varieties

abstract: We show that any polarized abelian variety over a finite field is covered by a Jacobian whose dimension is bounded by an explicit constant. We do this by first proving an effective version of Poonen's Bertini theorem over finite fields, which allows us to show the existence of smooth curves arising as hypersurface sections of bounded degree and genus. Additionally, we show that for simple abelian varieties a better bound is possible. As an application of these results we show that if $E$ is an elliptic curve over a finite field then for any $n\in \mathbb{N}$ there exist smooth curves of bounded genus whose Jacobians have a factor isogenous to $E^n$.


title: Generic H\"older foliations with smooth leaves

abstract: In this work, we consider a specific space of foliations with $C^1$ leaves and H\"older holonomies of the square $M=[0,1]^2$, with some topology and we show that a generic such foliation is non-absolutely continuous, furthermore, the conditional measures defined by Rokhlin disintegration are Dirac measures on the leaves. This space of foliations is motivated by the foliations that appear in hyperbolic systems and partially hyperbolic systems.


title: Site-Directed Insertion: Decision Problems, Maximality and Minimality

abstract: Site-directed insertion is an overlapping insertion operation that can be viewed as analogous to the overlap assembly or chop operations that concatenate strings by overlapping a suffix and a prefix of the argument strings. We consider decision problems and language equations involving site-directed insertion. By relying on the tools provided by semantic shuffle on trajectories we show that one variable equations involving site-directed insertion and regular constants can be solved. We consider also maximal and minimal variants of the site-directed insertion operation.


title: DAMA annual modulation from electron recoils

abstract: Plasma dark matter, which arises in dissipative dark matter models, can give rise to large annual modulation signals from keV electron recoils. Previous work has argued that the DAMA annual modulation signal could have an explanation within such a scenario. However, detailed predictions are difficult due to the inherent complexities involved in modelling the halo plasma interactions with Earth-bound dark matter. Here, we consider a simple phenomenological model for the dark matter velocity function relevant for direct detection experiments, and compare the resulting electron scattering rate with the new DAMA/LIBRA phase 2 data. We also consider the constraints from other experiments, including XENON100 and DarkSide-50.


title: Recurrent Entity Networks with Delayed Memory Update for Targeted Aspect-based Sentiment Analysis

abstract: While neural networks have been shown to achieve impressive results for sentence-level sentiment analysis, targeted aspect-based sentiment analysis (TABSA) --- extraction of fine-grained opinion polarity w.r.t. a pre-defined set of aspects --- remains a difficult task. Motivated by recent advances in memory-augmented models for machine reading, we propose a novel architecture, utilising external "memory chains" with a delayed memory update mechanism to track entities. On a TABSA task, the proposed model demonstrates substantial improvements over state-of-the-art approaches, including those using external knowledge bases.


title: Revisiting Norm Optimization for Multi-Objective Black-Box Problems: A Finite-Time Analysis

abstract: The complexity of Pareto fronts imposes a great challenge on the convergence analysis of multi-objective optimization methods. While most theoretical convergence studies have addressed finite-set and/or discrete problems, others have provided probabilistic guarantees, assumed a total order on the solutions, or studied their asymptotic behaviour. In this paper, we revisit the Tchebycheff weighted method in a hierarchical bandits setting and provide a finite-time bound on the Pareto-compliant additive $\epsilon$-indicator. To the best of our knowledge, this paper is one of few that establish a link between weighted sum methods and quality indicators in finite time.


title: On the Effect of Suboptimal Estimation of Mutual Information in Feature Selection and Classification

abstract: This paper introduces a new property of estimators of the strength of statistical association, which helps characterize how well an estimator will perform in scenarios where dependencies between continuous and discrete random variables need to be rank ordered. The new property, termed the estimator response curve, is easily computable and provides a marginal distribution agnostic way to assess an estimator's performance. It overcomes notable drawbacks of current metrics of assessment, including statistical power, bias, and consistency. We utilize the estimator response curve to test various measures of the strength of association that satisfy the data processing inequality (DPI), and show that the CIM estimator's performance compares favorably to kNN, vME, AP, and H_{MI} estimators of mutual information. The estimators which were identified to be suboptimal, according to the estimator response curve, perform worse than the more optimal estimators when tested with real-world data from four different areas of science, all with varying dimensionalities and sizes.


title: Adversarial Regression for Detecting Attacks in Cyber-Physical Systems

abstract: Attacks in cyber-physical systems (CPS) which manipulate sensor readings can cause enormous physical damage if undetected. Detection of attacks on sensors is crucial to mitigate this issue. We study supervised regression as a means to detect anomalous sensor readings, where each sensor's measurement is predicted as a function of other sensors. We show that several common learning approaches in this context are still vulnerable to \emph{stealthy attacks}, which carefully modify readings of compromised sensors to cause desired damage while remaining undetected. Next, we model the interaction between the CPS defender and attacker as a Stackelberg game in which the defender chooses detection thresholds, while the attacker deploys a stealthy attack in response. We present a heuristic algorithm for finding an approximately optimal threshold for the defender in this game, and show that it increases system resilience to attacks without significantly increasing the false alarm rate.


title: Quantum Rotor Engines

abstract: This chapter presents autonomous quantum engines that generate work in the form of directed motion for a rotor. We first formulate a prototypical clock-driven model in a time-dependent framework and demonstrate how it can be translated into an autonomous engine with the introduction of a planar rotor degree of freedom. The rotor plays both the roles of internal engine clock and of work repository. Using the example of a single-qubit piston engine, the thermodynamic performance is then reviewed. We evaluate the extractable work in terms of ergotropy, the kinetic energy associated to net directed rotation, as well as the intrinsic work based on the exerted torque under autonomous operation; and we compare them with the actual energy output to an external dissipative load. The chapter closes with a quantum-classical comparison of the engine's dynamics. For the single-qubit piston example, we propose two alternative representations of the qubit in an entirely classical framework: (i) a coin flip model and (ii) a classical magnet moment, showing subtle differences between the quantum and classical descriptions.


title: Adversarial Image Registration with Application for MR and TRUS Image Fusion

abstract: Robust and accurate alignment of multimodal medical images is a very challenging task, which however is very useful for many clinical applications. For example, magnetic resonance (MR) and transrectal ultrasound (TRUS) image registration is a critical component in MR-TRUS fusion guided prostate interventions. However, due to the huge difference between the image appearances and the large variation in image correspondence, MR-TRUS image registration is a very challenging problem. In this paper, an adversarial image registration (AIR) framework is proposed. By training two deep neural networks simultaneously, one being a generator and the other being a discriminator, we can obtain not only a network for image registration, but also a metric network which can help evaluate the quality of image registration. The developed AIR-net is then evaluated using clinical datasets acquired through image-fusion guided prostate biopsy procedures and promising results are demonstrated.


title: Field-Free Switching of Perpendicular Magnetic Tunnel Junction via Voltage-Gated Spin Hall Effect for Low-Power Spintronic Memory

abstract: Spin Hall effect (SHE) and voltage-controlled magnetic anisotropy (VCMA) are two promising methods for low-power electrical manipulation of magnetization. Recently, magnetic field-free switching of perpendicular magnetization through SHE has been reported with the aid of an exchange bias from an antiferromagnetic IrMn layer. In this letter, we experimentally demonstrate that the IrMn/CoFeB/MgO structure exhibits a VCMA effect of 39 fJ/Vm, which is comparable to that of the Ta/CoFeB/MgO structure. Magnetization dynamics under a combination of the SHE and VCMA are modeled and simulated. It is found that, by applying a voltage of 1.5 V, the critical SHE switching current can be decreased by 10 times owing to the VCMA effect, leading to low-power operations. Furthermore, a high-density spintronic memory structure can be built with multiple magnetic tunnel junctions (MTJs) located on a single IrMn strip. Through hybrid CMOS/MTJ simulations, we demonstrate that fast-speed write operations can be achieved with power consumption of only 8.5 fJ/bit. These findings reveal the possibility to realize high-density and low-power spintronic memory manipulated by voltage-gated SHE.


title: A unified software framework for solving traffic assignment problems

abstract: We describe a software framework for solving user equilibrium traffic assignment problems. The design is based on the formulation of the problem as a variational inequality. The software implements these as well as several numerical methods for find equilirbria. We compare the solutions obtained under several models: static, Merchant-Nemhauser, `CTM with instantaneous travel time', and `CTM with actual travel time'. Some important differences are demonstrated.


title: Deep Co-attention based Comparators For Relative Representation Learning in Person Re-identification

abstract: Person re-identification (re-ID) requires rapid, flexible yet discriminant representations to quickly generalize to unseen observations on-the-fly and recognize the same identity across disjoint camera views. Recent effective methods are developed in a pair-wise similarity learning system to detect a fixed set of features from distinct regions which are mapped to their vector embeddings for the distance measuring. However, the most relevant and crucial parts of each image are detected independently without referring to the dependency conditioned on one and another. Also, these region based methods rely on spatial manipulation to position the local features in comparable similarity measuring. To combat these limitations, in this paper we introduce the Deep Co-attention based Comparators (DCCs) that fuse the co-dependent representations of the paired images so as to focus on the relevant parts of both images and produce their \textit{relative representations}. Given a pair of pedestrian images to be compared, the proposed model mimics the foveation of human eyes to detect distinct regions concurrent on both images, namely co-dependent features, and alternatively attend to relevant regions to fuse them into the similarity learning. Our comparator is capable of producing dynamic representations relative to a particular sample every time, and thus well-suited to the case of re-identifying pedestrians on-the-fly. We perform extensive experiments to provide the insights and demonstrate the effectiveness of the proposed DCCs in person re-ID. Moreover, our approach has achieved the state-of-the-art performance on three benchmark data sets: DukeMTMC-reID \cite{DukeMTMC}, CUHK03 \cite{FPNN}, and Market-1501 \cite{Market1501}.


title: IMAGINE: Testing a Bayesian pipeline for Galactic Magnetic Field model optimization

abstract: This work contains the details and results of my master's project on testing the IMAGINE pipeline for Galactic magnetic field estimation. The project was carried out from early 2016 to early 2017. For it, an unpublished early development version of the IMAGINE pipeline was tested and debugged. The thesis reports about the kind of difficulties faced when dealing with high dimensional complex parametric Galactic magnetic field models. It was found that such models require extra caution to allow for dependencies between parameters and model implementation errors, which need to be taken into account when performing a Bayesian analysis. These findings, reported here in this thesis, helped to resolve such issues in the later, now published version of the IMAGINE pipeline. The thesis therefore documents the genesis of the pipeline and lessons learned during this process. This document contains original text of the master thesis for reference. Parts of its content therefore do not reflect the current state of the IMAGINE pipeline.


title: A New Nonconvex Strategy to Affine Matrix Rank Minimization Problem

abstract: The affine matrix rank minimization (AMRM) problem is to find a matrix of minimum rank that satisfies a given linear system constraint. It has many applications in some important areas such as control, recommender systems, matrix completion and network localization. However, the problem (AMRM) is NP-hard in general due to the combinational nature of the matrix rank function. There are many alternative functions have been proposed to substitute the matrix rank function, which lead to many corresponding alternative minimization problems solved efficiently by some popular convex or nonconvex optimization algorithms. In this paper, we propose a new nonconvex function, namely, $TL_{\alpha}^{\epsilon}$ function (with $0\leq\alpha<1$ and $\epsilon>0$), to approximate the rank function, and translate the NP-hard problem (AMRM) into the $TL_{p}^{\epsilon}$ function affine matrix rank minimization (TLAMRM) problem. Firstly, we study the equivalence of problem (AMRM) and (TLAMRM), and proved that the uniqueness of global minimizer of the problem (TLAMRM) also solves the NP-hard problem (AMRM) if the linear map $\mathcal{A}$ satisfies a restricted isometry property (RIP). Secondly, an iterative thresholding algorithm is proposed to solve the regularization problem (RTLAMRM) for all $0\leq\alpha<1$ and $\epsilon>0$. At last, some numerical results on low-rank matrix completion problems illustrated that our algorithm is able to recover a low-rank matrix, and the extensive numerical on image inpainting problems shown that our algorithm performs the best in finding a low-rank image compared with some state-of-art methods.


title: Trions and biexcitons in a nanowire

abstract: A theory of the trion and biexciton in a nanowire (NW) in the framework of the effective-mass model using the Born-Oppenheimer approximation is presented. We consider the formation of trions and biexcitons under the action of both the lateral confinement and the localization potential. The analytical expressions for the binding energy and eigenfunctions of the trion and biexciton are obtained and expressed by means of matrix elements of the effective one-dimensional cusp-type Coulomb potentials whose parameters are determined self-consistently by employing the same eigenfunctions of the confined electron and hole states. Our calculations for the ZnO/ZnMgO, CdSe/ZnS and CdSe/CdS core/shell cylindrical shaped NWs show that the trion and biexciton binding energy in NWs are size-dependent and for the same input parameters the biexciton binding energy in NWs is always larger than the binding energy of the trion. The trion and biexciton remain stable in CdSe/ZnS NW with the increase of the dielectric shell, while in ZnO/ZnMgO NW they become unstable when the surrounding dielectric shell exceeds 2.5 nm and 2 nm for each, respectively. The associative ionization of biexciton antibonding states into trion bonding states is studied. Based on the results for size dependence of biexciton binding energy and probability associative ionization an optimal radius for optoelectronic application NW is suggested.


title: Towards Deeper Generative Architectures for GANs using Dense connections

abstract: In this paper, we present the result of adopting skip connections and dense layers, previously used in image classification tasks, in the Fisher GAN implementation. We have experimented with different numbers of layers and inserting these connections in different sections of the network. Our findings suggests that networks implemented with the connections produce better images than the baseline, and the number of connections added has only slight effect on the result.


title: A criterion for the existence of Killing vectors in 3D

abstract: A three-dimensional Riemannian manifold has locally 6, 4, 3, 2, 1 or none independent Killing vectors. We present an explicit algorithm for computing dimension of the infinitesimal isometry algebra. It branches according to the values of curvature invariants. These are relative differential invariants computed via curvature, but they are not scalar polynomial Weyl invariants. We compare our obstructions to the existence of Killing vectors with the known existence criteria due to Singer, Kerr and others.


title: Areas of spherical polyhedral surfaces with regular faces

abstract: For a finite planar graph, it associates with some metric spaces, called (regular) spherical polyhedral surfaces, by replacing faces with regular spherical polygons in the unit sphere and gluing them edge-to-edge. We consider the class of planar graphs which admit spherical polyhedral surfaces with the curvature bounded below by 1 in the sense of Alexandrov, i.e. the total angle at each vertex is at most $2\pi$. We classify all spherical tilings with regular spherical polygons, i.e. total angles at vertices are exactly $2\pi$. We prove that for any graph in this class which does not admit a spherical tiling, the area of the associated spherical polyhedral surface with the curvature bounded below by 1 is at most $4\pi - \epsilon_0$ for some $\epsilon_0 > 0$. That is, we obtain a definite gap between the area of such a surface and that of the unit sphere.


title: Channel thickness optimization for ultra thin and 2D chemically doped TFETs

abstract: 2D material based tunnel FETs are among the most promising candidates for low power electronics applications since they offer ultimate gate control and high current drives that are achievable through small tunneling distances during the device operation. The ideal device is characterized by a minimized tunneling distance. However, devices with the thinnest possible body do not necessarily provide the best performance. For example, reducing the channel thickness increases the depletion width in the source which can be a significant part of the total tunneling distance. Hence, it is important to determine the optimum channel thickness for each channel material individually. In this work, we study the optimum channel thickness for three channel materials: WSe$_{2}$, Black Phosphorus (BP), and InAs using full-band self-consistent quantum transport simulations. To identify the ideal channel thickness for each material at a specific doping density, a new analytic model is proposed and benchmarked against the numerical simulations.


title: Partitions and periodic sequences

abstract: In the first part we associate a periodic sequence to a partition and study the connection the distribution of elements of uniform limit of the sequences. Then some facts of statistical independence of these limits are proved


title: A general framework for boundary equilibrium bifurcations of Filippov systems

abstract: As parameters are varied a boundary equilibrium bifurcation (BEB) occurs when an equilibrium collides with a discontinuity surface in a piecewise-smooth system of ODEs. Under certain genericity conditions, at a BEB the equilibrium either transitions to a pseudo-equilibrium (on the discontinuity surface) or collides and annihilates with a coexisting pseudo-equilibrium. These two scenarios are distinguished by the sign of a certain inner product. Here it is shown that this sign can be determined from the number of unstable directions associated with the two equilibria by using techniques developed by Feigin. A new normal form is proposed for BEBs in systems of any number of dimensions. The normal form involves a companion matrix, as does the leading order sliding dynamics, and so the connection to the stability of the equilibria is explicit. In two dimensions the parameters of the normal form distinguish, in a simple way, the eight topologically distinct cases for the generic local dynamics at a BEB. A numerical exploration in three dimensions reveals that BEBs can create multiple attractors and chaotic attractors, and that the equilibrium at the BEB can be unstable even if both equilibria are stable. The developments presented here stem from seemingly unutilised similarities between BEBs in discontinuous systems (specifically Filippov systems as studied here) and BEBs in continuous systems for which analogous results are, to date, more advanced.


title: Engineering non-equilibrium quantum phase transitions via causally gapped Hamiltonians

abstract: We introduce a phenomenological theory for many-body control of critical phenomena by engineering causally-induced gaps for quantum Hamiltonian systems. The core mechanisms are controlling information flow within and/or between clusters that are created near a quantum critical point. To this end, we construct inhomogeneous quantum phase transitions via designing spatio-temporal quantum fluctuations. We show how non-equilibrium evolution of disordered quantum systems can create new effective correlation length scales and effective dynamical critical exponents. In particular, we construct a class of causally-induced non-adiabatic quantum annealing transitions for strongly disordered quantum Ising chains leading to exponential suppression of topological defects beyond standard Kibble-Zurek predictions. Using exact numerical techniques for 1D quantum Hamiltonian systems, we demonstrate that our approach exponentially outperform adiabatic quantum computing. Using Strong-Disorder Renormalization Group (SDRG), we demonstrate the universality of inhomogeneous quantum critical dynamics and exhibit the causal zones reconstructions during SDRG flow. We derive a scaling relation for minimal causal gaps showing they narrow more slowly than any polynomial with increasing size of system, in contrast to stretched exponential scaling in standard adiabatic evolution. Furthermore, we demonstrate similar scaling behaviour for random cluster-Ising Hamiltonians with higher order interactions.


title: Holistic Management of the GPGPU Memory Hierarchy to Manage Warp-level Latency Tolerance

abstract: In a modern GPU architecture, all threads within a warp execute the same instruction in lockstep. For a memory instruction, this can lead to memory divergence: the memory requests for some threads are serviced early, while the remaining requests incur long latencies. This divergence stalls the warp, as it cannot execute the next instruction until all requests from the current instruction complete. In this work, we make three new observations. First, GPGPU warps exhibit heterogeneous memory divergence behavior at the shared cache: some warps have most of their requests hit in the cache, while other warps see most of their request miss. Second, a warp retains the same divergence behavior for long periods of execution. Third, requests going to the shared cache can incur queuing delays as large as hundreds of cycles, exacerbating the effects of memory divergence. We propose a set of techniques, collectively called Memory Divergence Correction (MeDiC), that reduce the negative performance impact of memory divergence and cache queuing. MeDiC delivers an average speedup of 21.8%, and 20.1% higher energy efficiency, over a state-of-the-art GPU cache management mechanism across 15 different GPGPU applications.