title: Leveraging the Template and Anchor Framework for Safe, Online Robotic Gait Design

abstract: Online control design using a high-fidelity, full-order model for a bipedal robot can be challenging due to the size of the state space of the model. A commonly adopted solution to overcome this challenge is to approximate the full-order model (anchor) with a simplified, reduced-order model (template), while performing control synthesis. Unfortunately it is challenging to make formal guarantees about the safety of an anchor model using a controller designed in an online fashion using a template model. To address this problem, this paper proposes a method to generate safety-preserving controllers for anchor models by performing reachability analysis on template models while bounding the modeling error. This paper describes how this reachable set can be incorporated into a Model Predictive Control framework to select controllers that result in safe walking on the anchor model in an online fashion. The method is illustrated on a 5-link RABBIT model, and is shown to allow the robot to walk safely while utilizing controllers designed in an online fashion.


title: Projective representations of the inhomogeneous symplectic group: Quantum symmetry origins of the Heisenberg commutation relations

abstract: Quantum symmetries that leave invariant physical transition probabilities are described by projective representations of Lie groups. The mathematical theory of projected representations for topologically connected Lie groups is reviewed and applied to the inhomogeneous symplectic group ISp(2n). The projective representations are given in terms of the ordinary unitary representations of the central extension of the group with respect to its first homotopy group direct product with its second cohomology group. The second cohomology group of ISp(2n) is the one dimensional abelian group and the central extension turns the 2n dimension abelian normal subgroup of translations of ISp(2n) into the Weyl-Heisenberg group. This is the quantum symmetry origin of the Heisenberg commutation relations that are the Hermitian representation of the Lie algebra of the Weyl-Heisenberg group.


title: Holographic line defects in F(4) gauged supergravity

abstract: In this note we construct a solution of six-dimensional $F(4)$ gauged supergravity using $AdS_2\times S^3$ warped over an interval as an ansatz. The solution is completely regular, preserves eight of the sixteen supersymmetries of the $AdS_6$ vacuum and is a holographic realization of a line defect in a dual five-dimensional theory. We calculate the expectation value of the defect and the one-point function of the stress tensor in the presence of the defect using holographic renormalization.


title: Accept Synthetic Objects as Real: End-to-End Training of Attentive Deep Visuomotor Policies for Manipulation in Clutter

abstract: Recent research demonstrated that it is feasible to end-to-end train multi-task deep visuomotor policies for robotic manipulation using variations of learning from demonstration (LfD) and reinforcement learning (RL). In this paper, we extend the capabilities of end-to-end LfD architectures to object manipulation in clutter. We start by introducing a data augmentation procedure called Accept Synthetic Objects as Real (ASOR). Using ASOR we develop two network architectures: implicit attention ASOR-IA and explicit attention ASOR-EA. Both architectures use the same training data (demonstrations in uncluttered environments) as previous approaches. Experimental results show that ASOR-IA and ASOR-EA succeed ina significant fraction of trials in cluttered environments where previous approaches never succeed. In addition, we find that both ASOR-IA and ASOR-EA outperform previous approaches even in uncluttered environments, with ASOR-EA performing better even in clutter compared to the previous best baseline in an uncluttered environment.


title: The suspension of a 4-manifold and its applications

abstract: Let $M$ be a smooth, orientable, closed, connected $4$-manifold and suppose that $H_1(M;\mathbb{Z})$ is finitely generated and has no $2$-torsion. We give a homotopy decomposition of the suspension of $M$ in terms of spheres, Moore spaces and $\Sigma\mathbb{C}P^{2}$. This is used to calculate any reduced generalized cohomology theory of $M$ as a group and to determine the homotopy types of certain current groups and gauge groups.


title: Simple Woods-Saxon type form for $\Omega \alpha $ and $\Xi \alpha$ interactions using Folding Model

abstract: We derive a simple Woods-Saxon type form for the potentials between $Y=\Xi, \Omega$ and $\alpha$ by using a single-folding potential method, based on a separable $Y$-nucleon potential. Accordingly, the potentials $\Xi+\alpha$ and $\Omega+\alpha$ are obtained using the ESC08c Nijmegens $\Xi N$ potential (in $^{3}S_{1}$ channel) and HAL QCD Collaboration $\Omega N$ interactions (in lattice QCD), respectively. In deriving the potential between $Y$ and $\alpha$, the same potential between $Y$ and $N$ is used. Binding energy, scattering length and effective range of $Y$ particle on the alpha particle are approximated by the resultant potentials. The depths of the potentials in $\Omega \alpha $ and $\Xi \alpha $ systems are obtained $-61$ and $-24.4$ MeV, respectively. In the case of $\Xi \alpha$ potential, a fairly good agreement is observed between the single-folding potential method and the phenomenological potential of Dover-Gal model. These potentials can be used in 3-,4- and 5-body cluster structures of $ \Omega$ and $\Xi$ hypernuclei.


title: Metrics and pseudometrics on unitary groups with applications in quantum information processing

abstract: Metrics and pseudometrics are defined on the group of unitary operators in a Hilbert space. Several explicit formulas are derived. A special feature of the work is investigation of pseudometrics in unitary groups. The rich classes of pseudometrics have many interesting applications. Three such applications, distinguishibility of unitary operators, quantum coding and quantum search problems are discussed.


title: Quasi-Assouad dimensions for random measures supported on $[0,1]^d$

abstract: We introduce a probability distribution on $\mathcal{P}([0,1]^d)$, the space of all Borel probability measures on $[0,1]^d$. Under this distribution, almost all measures are shown to have infinite upper quasi-Assouad dimension and zero lower quasi-Assouad dimension (hence the upper and lower Assouad dimensions are almost surely infinite or zero). We also indicate how the results extend to other Assouad-like dimensions.


title: Inverse problems for Schr\"odinger equations with unbounded potentials

abstract: We summarize in these notes the course given at the Summer School of AIP 2019 held in Grenoble from July 1st to July 5th. This course was mainly devoted to the determination of the unbounded potential in a Schr\"odinger equation from the associated Dirichlet-to-Neumann map (abbreviated to DN map in this text). We establish a stability inequality for potentials belonging to $L^n$, where $n\ge 3$ is the dimension of the space. Next, we prove a uniqueness result for potentials in $L^{n/2}$, $n\ge 3$, and apply this uniqueness result to demonstrate a Borg-Levinson type theorem. We use a classical approach which is essentially based on the construction of the so-called complex geometric optic solutions (abbreviated to CGO solutions in this text).


title: Structural phase transitions in VSe2: energetics, electronic structure and magnetism

abstract: First principles calculations of magnetic and electronic properties of VSe2 describing the transition between two structural phases(H,T) were performed. Results of the calculations evidence rather low energy barrier ( 0.60 eV for monolayer) for transition betweenthe phases. The energy required for the deviation of Se atom or whole layer of selenium atoms on a small angle up to 10 degrees from initial positions is also rather low, 0.32 and 0.19 eV/Se, respectively. The changes in band structure of VSe2 caused by these motions of Se atoms should be taken into account for analysis of the experimental data. Simulations of the strain effects suggest that the experimentally observed T phase of VSe2 monolayer is the ground state due a substrate-induced strain. Calculations of the difference in total energies of ferromagnetic and antiferromagnetic configurations evidence that the ferromagnetic configuration is the ground state of the system for all stable and intermediate atomic structures. Calculated phonon dispersions suggest visible influence of magnetic configurations on vibrational properties.


title: Precision Small Scattering Angle Measurements of Elastic Proton-Proton Single and Double Spin Analyzing Powers at the RHIC Hydrogen Jet Polarimeter

abstract: The Polarized Atomic Hydrogen Gas Jet Target polarimeter is employed by the Relativistic Heavy Ion Collider (RHIC) to measure the absolute polarization of each colliding proton beam. Polarimeter detectors and data acquisition were upgraded in 2015 to increase solid angle, energy range and energy resolution. These upgrades and advanced systematic error analysis along with improved beam intensity and polarization in RHIC runs 2015 ($E_\textrm{beam}=100\,\text{GeV}$) and 2017 ($255\,\text{GeV}$) allowed us to greatly reduce the statistical and systematic uncertainties for elastic spin asymmetries, $A_N(t)$ and $A_NN(t)$, in the Coulomb-nuclear interference momentum transfer range $0.0013<-t<0.018\,\text{GeV}^2$. For the first time hadronic single spin-flip $r_5$ and double spin-flip $r_2$ amplitude parameters were reliably isolated at these energies and momentum transfers. Measurements at two beam energies enable a separation of Pomeron and Regge pole contributions to $r_5(s)$ and $r_2(s)$, indicating that the spin component may persist at high energies.


title: Implementation of on-chip multi-channel focusing wavelength demultiplexing with regularized digital metamaterials

abstract: Adiabatic waveguide taper and on-chip wavelength demultiplexer are the key components of photonic integrated circuits. However, these two kinds of devices which designed by traditional semi-analytic methods or brute-force search methods usually have large size. Here, based on regularized digital metamaterials, we have designed, fabricated and characterized a two-channel focused wavelength demultiplexer with a footprint of 2.4 x 10 um2. The designed demultiplexer can directly connect to a grating coupler under the absence of an adiabatic waveguide taper. The objective first method and modified steepest descent method are used to design the demultiplexer which splits 1520 nm and 1580 nm light from a 10-um-wide input waveguide into two 0.48-um-wide output waveguides. Experimental results show that the insertion loss of the upper (lower) channel of the demultiplexer is -1.77 dB (-2.10 dB) and the crosstalk is -25.17 dB (-12.14 dB). Besides, the simulation results indicate that the fabrication tolerance of our devices can reach 20 nm in etching depth and 10 nm in plane size changing. Benefit From the extensibility of our design method, we can design other types of ultra-compact 'focused' devices, like mode splitters, mode converters and power splitters, and we can also design devices with more complicated functionalities, for example, we have designed a three-channel focused wavelength demultiplexer.


title: On decay of shock-like waves into compact oscillons

abstract: The signum-Gordon model in 1+1 dimensions possesses the exact shockwave solution with discontinuity of the field at the light cone and infinite gradient energy. The energy of a regular part of the wave inside the light cone is finite and it grows linearly with time. The initial data for such waves contain a field configuration which is null in the space and has time derivative proportional to the Dirac delta. We study regularized initial data that lead to shock-like waves with finite gradient energy. We found that such waves exist in the finite time intervals and finally they decay and produce a cascade of oscillon-like structures. A pattern of the decay is very similar to the one observed in process of scattering of compact oscillons.


title: Quantization of spin Hall conductivity in two-dimensional topological insulators versus symmetry and spin-orbit interaction

abstract: The third-rank tensor of the static spin Hall conductivity is investigated for two-dimensional (2D) topological insulators by electronic structure calculations. Its seeming quantization is numerically demonstrated for highly symmetric systems independent of the gap size. 2D crystals with hexagonal and square Bravais lattice show similar effects, while true rectangular translational symmetry yields conductivity values much below the quantum $e^2/h$. Field-induced lifting the inversion symmetry does not influence the quantum spin Hall state up to band inversion but the conductivity quantization. Weak symmetry-conserving biaxial but also uniaxial strain has a minor influence as long as inverted gaps dictate the topological character. The results are discussed in terms of the atomic geometry and the Rashba contribution to the spin-orbit interaction (SOI). Translational and point-group symmetry as well as SOI rule the deviation from the quantization of the spin Hall conductance.


title: Thin Loop Groups

abstract: We verify that for a finite simplicial complex $X$ and for piecewise linear loops on $X$, the "thin" loop space is a topological group of the same homotopy type as the space of continuous loops. This turns out not to be the case for the higher loops.


title: Lick Observatory Supernova Search Follow-Up Program: Photometry Data Release of 93 Type Ia Supernovae

abstract: We present BVRI and unfiltered light curves of 93 Type Ia supernovae (SNe Ia) from the Lick Observatory Supernova Search (LOSS) follow-up program conducted between 2005 and 2018. Our sample consists of 78 spectroscopically normal SNe Ia, with the remainder divided between distinct subclasses (three SN 1991bg-like, three SN 1991T-like, four SNe Iax, two peculiar, and three super-Chandrasekhar events), and has a median redshift of 0.0192. The SNe in our sample have a median coverage of 16 photometric epochs at a cadence of 5.4 days, and the median first observed epoch is ~4.6 days before maximum B-band light. We describe how the SNe in our sample are discovered, observed, and processed, and we compare the results from our newly developed automated photometry pipeline to those from the previous processing pipeline used by LOSS. After investigating potential biases, we derive a final systematic uncertainty of 0.03 mag in BVRI for our dataset. We perform an analysis of our light curves with particular focus on using template fitting to measure the parameters that are useful in standardising SNe Ia as distance indicators. All of the data are available to the community, and we encourage future studies to incorporate our light curves in their analyses.


title: Sleep Stage Classification Using Bidirectional LSTM in Wearable Multi-sensor Systems

abstract: Understanding the sleep quality and architecture is essential to human being's health, which is usually represented using multiple sleep stages. A standard sleep stage determination requires Electroencephalography (EEG) signals during the expensive and labor-intensive Polysomnography (PSG) test. To overcome this inconvenience, cardiorespiratory signals are proposed for the same purpose because of the easy and comfortable acquisition by simplified devices. In this paper, we leverage our low-cost wearable multi-sensor system to acquire the cardiorespiratory signals from subjects. Three novel features are designed during the feature extraction. We then apply a Bi-directional Recurrent Neural Network architecture with Long Short-term Memory (BLSTM) to predict the four-class sleep stages. Our prediction accuracy is 80.25% on a large public dataset (417 subjects), and 80.75% on our 32 enrolled subjects, respectively. Our results outperform the previous works which either used small data sets and had the potential over-fitting issues, or used the conventional machine learning methods on large data sets.


title: CAGE: Context-Aware Grasping Engine

abstract: Semantic grasping is the problem of selecting stable grasps that are functionally suitable for specific object manipulation tasks. In order for robots to effectively perform object manipulation, a broad sense of contexts, including object and task constraints, needs to be accounted for. We introduce the Context-Aware Grasping Engine, which combines a novel semantic representation of grasp contexts with a neural network structure based on the Wide & Deep model, capable of capturing complex reasoning patterns. We quantitatively validate our approach against three prior methods on a novel dataset consisting of 14,000 semantic grasps for 44 objects, 7 tasks, and 6 different object states. Our approach outperformed all baselines by statistically significant margins, producing new insights into the importance of balancing memorization and generalization of contexts for semantic grasping. We further demonstrate the effectiveness of our approach on robot experiments in which the presented model successfully achieved 31 of 32 suitable grasps. The code and data are available at: https://github.com/wliu88/rail_semantic_grasping


title: Spontaneous Fruit Fly Optimisation for truss weight minimisation: Performance evaluation based on the no free lunch theorem

abstract: Over the past decade, several researchers have presented various optimisation algorithms for use in truss design. The no free lunch theorem implies that no optimisation algorithm fits all problems; therefore, the interest is not only in the accuracy and convergence rate of the algorithm but also the tuning effort and population size required for achieving the optimal result. The latter is particularly crucial for computationally intensive or high-dimensional problems. Contrast-based Fruit-fly Optimisation Algorithm (c-FOA) proposed by Kanarachos et al. in 2017 is based on the efficiency of fruit flies in food foraging by olfaction and visual contrast. The proposed Spontaneous Fruit Fly Optimisation (s-FOA) enhances c-FOA and addresses the difficulty in solving nonlinear optimisation algorithms by presenting standard parameters and lean population size for use on all optimisation problems. Six benchmark problems were studied to assess the performance of s-FOA. A comparison of the results obtained from documented literature and other investigated techniques demonstrates the competence and robustness of the algorithm in truss optimisation.


title: Ab Initio Linear and Pump-Probe Spectroscopy of Excitons in Molecular Crystals

abstract: Linear and non-linear spectroscopies are powerful tools used to investigate the energetics and dynamics of electronic excited states of both molecules and crystals. While highly accurate \emph{ab initio} calculations of molecular spectra can be performed relatively routinely, extending these calculations to periodic systems is challenging. Here, we present calculations of the linear absorption spectrum and pump-probe two-photon photoemission spectra of the naphthalene crystal using equation-of-motion coupled-cluster theory with single and double excitations (EOM-CCSD). Molecular acene crystals are of interest due to the low-energy multi-exciton singlet states they exhibit, which have been studied extensively as intermediates involved in singlet fission. Our linear absorption spectrum is in good agreement with experiment, predicting a first exciton absorption peak at 4.4 eV, and our two-photon photoemission spectra capture the behavior of multi-exciton states, whose double-excitation character cannot be captured by current methods. The simulated pump-probe spectra provide support for existing interpretations of two-photon photoemission in closely-related acene crystals such as tetracene and pentacene.


title: Brain-Inspired Hardware for Artificial Intelligence: Accelerated Learning in a Physical-Model Spiking Neural Network

abstract: Future developments in artificial intelligence will profit from the existence of novel, non-traditional substrates for brain-inspired computing. Neuromorphic computers aim to provide such a substrate that reproduces the brain's capabilities in terms of adaptive, low-power information processing. We present results from a prototype chip of the BrainScaleS-2 mixed-signal neuromorphic system that adopts a physical-model approach with a 1000-fold acceleration of spiking neural network dynamics relative to biological real time. Using the embedded plasticity processor, we both simulate the Pong arcade video game and implement a local plasticity rule that enables reinforcement learning, allowing the on-chip neural network to learn to play the game. The experiment demonstrates key aspects of the employed approach, such as accelerated and flexible learning, high energy efficiency and resilience to noise.


title: Computational Chemistry on Quantum Computers

abstract: The purpose of this experiment was to use the known analytical techniques to study the creation, simulation, and measurements of molecular Hamiltonians. The techniques used consisted of the Linear Combination of Atomic Orbitals (LCAO), the Linear Combination of Unitaries (LCU), and the Phase Estimation Algorithm (PEA). The molecules studied were $H_2$ with and without spin, as well as $He_2$ without spin. Hamiltonians were created under the LCAO basis, and reconstructed using the Jordan-Winger transform in order to create a linear combination of Pauli spin operators. The lengths of each molecular Hamiltonian greatly increased from the $H_2$ without spin, to $He_2$. This resulted in a reduced ability to simulate the Hamiltonians under ideal conditions. Thus, only low orders of l = 1 and l = 2 were used when expanding the Hamiltonian in accordance to the LCU method of simulation. The resulting Hamiltonians were measured using PEA, and plotted against function of $\frac{2\pi(K)}{N}$ and the probability distribution of each register. The resolution of the graph was dependent on the amount of registers, N, being used. However, the reduction of order hardly changed the image of the $H_2$ graphs. Qualitative comparisons between the three molecules were drawn.


title: Random $k$-out subgraph leaves only $O(n/k)$ inter-component edges

abstract: Each vertex of an arbitrary simple graph on $n$ vertices chooses $k$ random incident edges. What is the expected number of edges in the original graph that connect different connected components of the sampled subgraph? We prove that the answer is $O(n/k)$, when $k\ge c\log n$, for some large enough $c$. We conjecture that the same holds for smaller values of $k$, possibly for any $k\ge 2$. Such a result is best possible for any $k\ge 2$. As an application, we use this sampling result to obtain a one-way communication protocol with \emph{private} randomness for finding a spanning forest of a graph in which each vertex sends only ${O}(\sqrt{n}\log n)$ bits to a referee.


title: Modelling Connective Spectra via Multicategories

abstract: We put a model structure on a full subcategory of based multicategories in which the weak equivalences are created by the K-theory functor of Elmendorf-Mandell, providing a model categorical lift of Thomason's theorem on the modeling of connective spectra by symmetric monoidal categories. We note that this lifts to a semi-modelstructure on based multicategories itself. As a corollary we show that to model connective spectra up to stable equivalence it suffices to restrict to symmetric monoidal groupoids.


title: Can each number be specified by a finite text?

abstract: Contrary to popular misconception, the question in the title is far from simple. It involves sets of numbers on the first level, sets of sets of numbers on the second level, and so on, endlessly. The infinite hierarchy of the levels involved distinguishes the concept of "definable number" from such notions as "natural number", "rational number", "algebraic number", "computable number" etc. (Explanatory essay for non-experts.)


title: Exascale Deep Learning for Scientific Inverse Problems

abstract: We introduce novel communication strategies in synchronous distributed Deep Learning consisting of decentralized gradient reduction orchestration and computational graph-aware grouping of gradient tensors. These new techniques produce an optimal overlap between computation and communication and result in near-linear scaling (0.93) of distributed training up to 27,600 NVIDIA V100 GPUs on the Summit Supercomputer. We demonstrate our gradient reduction techniques in the context of training a Fully Convolutional Neural Network to approximate the solution of a longstanding scientific inverse problem in materials imaging. The efficient distributed training on a dataset size of 0.5 PB, produces a model capable of an atomically-accurate reconstruction of materials, and in the process reaching a peak performance of 2.15(4) EFLOPS$_{16}$.

categories: cs.LG cond-mat.mtrl-sci cs.DC physics.comp-ph stat.ML


title: K-Motives and Koszul Duality

abstract: We construct an ungraded version of Beilinson-Ginzburg-Soergel's Koszul duality for Langlands dual flag varieties, inspired by Beilinson's construction of rational motivic cohomology in terms of $K$-theory. For this, we introduce and study categories $\operatorname{DK}_{\mathcal{S}}(X)$ of $\mathcal{S}$-constructible $K$-motivic sheaves on varieties $X$ with an affine stratification $\mathcal{S}$. We show that there is a natural and geometric functor, called Beilinson realisation, from $\mathcal{S}$-constructible mixed sheaves $\operatorname{D}^{mix}_{\mathcal{S}}(X)$ to $\operatorname{DK}_{\mathcal{S}}(X)$. We then show that Koszul duality intertwines the Betti realisation and Beilinson realisation functors and descends to an equivalence of constructible sheaves and constructible $K$-motivic sheaves on Langlands dual flag varieties.


title: On the edge-length ratio of 2-trees

abstract: We study planar straight-line drawings of graphs that minimize the ratio between the length of the longest and the shortest edge. We answer a question of Lazard et al. [Theor. Comput. Sci. 770 (2019), 88--94] and, for any given constant $r$, we provide a $2$-tree which does not admit a planar straight-line drawing with a ratio bounded by $r$. When the ratio is restricted to adjacent edges only, we prove that any $2$-tree admits a planar straight-line drawing whose edge-length ratio is at most $4 + \varepsilon$ for any arbitrarily small $\varepsilon > 0$, hence the upper bound on the local edge-length ratio of partial $2$-trees is $4$.


title: Dimension-free estimates for Riesz transforms related to the harmonic oscillator

abstract: We study $L^p$ bounds for two kinds of Riesz transforms on $\mathbb{R}^d$ related to the harmonic oscillator. We pursue an explicit estimate of their $L^p$ norms that is independent of the dimension $d$ and linear in $\max(p, p/(p-1))$.


title: Social contagion models on hypergraphs

abstract: Our understanding of the dynamics of complex networked systems has increased significantly in the last two decades. However, most of our knowledge is built upon assuming pairwise relations among the system's components. This is often an oversimplification, for instance, in social interactions that occur frequently within groups. To overcome this limitation, here we study the dynamics of social contagion on hypergraphs. We develop an analytical framework and provide numerical results for arbitrary hypergraphs, which we also support with Monte Carlo simulations. Our analyses show that the model has a vast parameter space, with first and second-order transitions, bi-stability, and hysteresis. Phenomenologically, we also extend the concept of latent heat to social contexts, which might help understanding oscillatory social behaviors. Our work unfolds the research line of higher-order models and the analytical treatment of hypergraphs, posing new questions and paving the way for modeling dynamical processes on these networks.


title: Anchor Loss: Modulating Loss Scale based on Prediction Difficulty

abstract: We propose a novel loss function that dynamically rescales the cross entropy based on prediction difficulty regarding a sample. Deep neural network architectures in image classification tasks struggle to disambiguate visually similar objects. Likewise, in human pose estimation symmetric body parts often confuse the network with assigning indiscriminative scores to them. This is due to the output prediction, in which only the highest confidence label is selected without taking into consideration a measure of uncertainty. In this work, we define the prediction difficulty as a relative property coming from the confidence score gap between positive and negative labels. More precisely, the proposed loss function penalizes the network to avoid the score of a false prediction being significant. To demonstrate the efficacy of our loss function, we evaluate it on two different domains: image classification and human pose estimation. We find improvements in both applications by achieving higher accuracy compared to the baseline methods.


title: Completely uniformly distributed sequences based on de Bruijn sequences

abstract: We study a construction published by Donald Knuth in 1965 yielding a completely uniformly distributed sequence of real numbers. Knuth's work is based on de Bruijn sequences of increasing orders and alphabet sizes, which grow exponentially in each of the successive segments composing the generated sequence. In this work we present a similar albeit simpler construction using linearly increasing alphabet sizes, and give an elementary proof showing that the sequence it yields is also completely uniformly distributed. In addition, we present an alternative proof of the same result based on Weyl's criterion.


title: A Robust and Unified Solution for Choosing the Phases of Adiabatic States as a Function of Geometry: Extending Parallel Transport Concepts to the cases of Trivial & Near Trivial Crossings

abstract: We investigate a simple and robust scheme for choosing the phases of adiabatic electronic states smoothly (as a function of geometry) so as to maximize the performance of ab initio non-adiabatic dynamics methods. Our approach is based upon consideration of the overlap matrix ($\mathbf{U}$) between basis functions at successive points in time and selecting the phases so as to minimize the matrix norm of $\log(\mathbf{U})$. In so doing, one can extend the concept of parallel transport to cases with sharp curve crossings. We demonstrate that this algorithm performs well under extreme situations where dozens of states cross each other either through trivial crossings (where there is zero effective diabatic coupling), or through nontrivial crossings (when there is a nonzero diabatic coupling), or through a combination of both. In all cases, we compute the time-derivative coupling matrix elements (or equivalently non-adiabatic derivative coupling matrix elements) that are as smooth as possible. Our results should be of interest to all who are interested in either non-adiabatic dynamics, or more generally, parallel transport in large systems.


title: The Imprint of Spiral Arms on the Galactic Rotation Curve

abstract: We discuss a model for the Milky Way obtained by fitting the observed terminal velocities with the radial acceleration relation. The resulting stellar surface density profile departs from a smooth exponential disk, having bumps and wiggles that correspond to massive spiral arms. These features are used to estimate the term for the logarithmic density gradient in the Jeans equation, which turn out to have exactly the right location and amplitude to reconcile the apparent discrepancy between the stellar rotation curve and that of the interstellar gas. This model also predicts a gradually declining rotation curve outside the solar circle with slope $-1.7\;\mathrm{km}\,\mathrm{s}^{-1}\,\mathrm{kpc}^{-1}$, as subsequently observed.


title: Efficient Automata-based Planning and Control under Spatio-Temporal Logic Specifications

abstract: The use of spatio-temporal logics in control is motivated by the need to impose complex spatial and temporal behavior on dynamical systems, and to control these systems accordingly. Synthesizing correct-by-design control laws is a challenging task resulting in computationally demanding methods. We consider efficient automata-based planning for continuous-time systems under signal interval temporal logic specifications, an expressive fragment of signal temporal logic. The planning is based on recent results for automata-based verification of metric interval temporal logic. A timed signal transducer is obtained accepting all Boolean signals that satisfy a metric interval temporal logic specification, which is abstracted from the signal interval temporal logic specification at hand. This transducer is modified to account for the spatial properties of the signal interval temporal logic specification, characterizing all real-valued signals that satisfy this specification. Using logic-based feedback control laws, such as the ones we have presented in earlier works, we then provide an abstraction of the system that, in a suitable way, aligns with the modified timed signal transducer. This allows to avoid the state space explosion that is typically induced by forming a product automaton between an abstraction of the system and the specification.


title: Energy Efficiency Maximization Via Joint Sub-Carrier Assignment and Power Control for OFDMA Full Duplex Networks

abstract: In this paper, we develop an energy efficient resource allocation scheme for orthogonal frequency division multiple access (OFDMA) networks with in-band full-duplex (IBFD) communication between the base station and user equipments (UEs) considering a realistic self-interference (SI) model. Our primary aim is to maximize the system energy efficiency (EE) through a joint power control and sub-carrier assignment in both the downlink (DL) and uplink (UL), where the quality of service requirements of the UEs in DL and UL are guaranteed. The formulated problem is non-convex due to the non-linear fractional objective function and the non-convex feasible set which is generally intractable. In order to handle this difficulty, we first use fractional programming to transform the fractional objective function to the subtractive form. Then, by employing Dinkelbach method, we propose an iterative algorithm in which an inner problem is solved in each iteration. Applying majorization-minimization approximation, we make the inner problem convex. Also, by introducing a penalty function to handle integer sub-carrier assignment variables, we propose an iterative algorithm for addressing the inner problem. We show that our proposed algorithm converges to the locally optimal solution which is also demonstrated by our simulation results. In addition, simulation results show that by applying the IBFD capability in OFDMA networks with efficient SI cancellation techniques, our proposed resource allocation algorithm attains a 75% increase in the EE as compared to the half-duplex system.


title: Selecting a Scale for Spatial Confounding Adjustment

abstract: Unmeasured, spatially-structured factors can confound associations between spatial environmental exposures and health outcomes. Adding flexible splines to a regression model is a simple approach for spatial confounding adjustment, but the spline degrees of freedom do not provide an easily interpretable spatial scale. We describe a method for quantifying the extent of spatial confounding adjustment in terms of the Euclidean distance at which variation is removed. We develop this approach for confounding adjustment with splines and using Fourier and wavelet filtering. We demonstrate differences in the spatial scales these bases can represent and provide a comparison of methods for selecting the amount of confounding adjustment. We find the best performance for selecting the amount of adjustment using an information criterion evaluated on an outcome model without exposure. We apply this method to spatial adjustment in an analysis of particulate matter and blood pressure in a cohort of United States women.


title: Braid group action on projective quantum $\mathfrak{sl}(2)$ modules

abstract: We define a family of the braid group representations via the action of the $R$-matrix (of the quasitriangular extension) of the restricted quantum $\mathfrak{sl}(2)$ on a tensor power of a simple projective module. This family is an extension of the Lawrence representation specialized at roots of unity. Although the center of the braid group has finite order on the specialized Laurence representations, this action is faithful for our extension.


title: Descriptive complexity of subsets of the space of finitely generated groups

abstract: In this paper, we determine the descriptive complexity of subsets of the Polish space of marked groups defined by various group theoretic properties. In particular, using Grigorchuk groups, we establish that the sets of solvable groups, groups of exponential growth and groups with decidable word problem are $\mathbf{\Sigma}^0_2$-complete and that the sets of periodic groups and groups of intermediate growth are $\mathbf{\Pi}^0_2$-complete. We also provide bounds for the descriptive complexity of simplicity, amenability, residually finiteness, Hopficity and co-Hopficity. This paper is intended to serve as a compilation of results on this theme.


title: P2FAAS: Toward Privacy-Preserving Fuzzing as a Service

abstract: Global corporations (e.g., Google and Microsoft) have recently introduced a new model of cloud services, fuzzing-as-a-service (FaaS). Despite effectively alleviating the cost of fuzzing, the model comes with privacy concerns. For example, the end user has to trust both cloud and service providers who have access to the application to be fuzzed. Such concerns are due to the platform is under the control of its provider and the application and the fuzzer are highly coupled. In this paper, we propose P2FaaS, a new ecosystem that preserves end user's privacy while providing FaaS in the cloud. The key idea of P2FaaS is to utilize Intel SGX for preventing cloud and service providers from learning information about the application. Our preliminary evaluation shows that P2FaaS imposes 45% runtime overhead to the fuzzing compared to the baseline. In addition, P2FaaS demonstrates that, with recently introduced hardware, Intel SGX Card, the fuzzing service can be scaled up to multiple servers without native SGX support.


title: Thermal cycling memory in phase separated manganites

abstract: We have studied the irreversibility of the magnetization induced by thermal cycles in La0.5Ca0.5MnO3 manganites, which present a low temperature state characterized by the coexistence of phases. The effect is evidenced by a decrease of the magnetization after cycling the sample between 300 and 50 K. We developed a phenomenological model that allows us to correlate the value of the magnetization with the number of cycles performed. The experimental results show excellent agreement with our model, suggesting that this material could be used for the development of a device to monitor thermal changes. The effect of thermal cycling is towards an increase of the amount of the non ferromagnetic phase in the compounds and it might be directly related with the strain at the contact surface among the coexisting phases.


title: Ethical Hacking for IoT Security: A First Look into Bug Bounty Programs and Responsible Disclosure

abstract: The security of the Internet of Things (IoT) has attracted much attention due to the growing number of IoT-oriented security incidents. IoT hardware and software security vulnerabilities are exploited affecting many companies and persons. Since the causes of vulnerabilities go beyond pure technical measures, there is a pressing demand nowadays to demystify IoT "security complex" and develop practical guidelines for both companies, consumers, and regulators. In this paper, we present an initial study targeting an unexplored sphere in IoT by illuminating the potential of crowdsource ethical hacking approaches for enhancing IoT vulnerability management. We focus on Bug Bounty Programs (BBP) and Responsible Disclosure (RD), which stimulate hackers to report vulnerability in exchange for monetary rewards. We carried out a qualitative investigation supported by literature survey and expert interviews to explore how BBP and RD can facilitate the practice of identifying, classifying, prioritizing, remediating, and mitigating IoT vulnerabilities in an effective and cost-efficient manner. Besides deriving tangible guidelines for IoT stakeholders, our study also sheds light on a systematic integration path to combine BBP and RD with existing security practices (e.g., penetration test) to further boost overall IoT security.


title: Intelligent image synthesis to attack a segmentation CNN using adversarial learning

abstract: Deep learning approaches based on convolutional neural networks (CNNs) have been successful in solving a number of problems in medical imaging, including image segmentation. In recent years, it has been shown that CNNs are vulnerable to attacks in which the input image is perturbed by relatively small amounts of noise so that the CNN is no longer able to perform a segmentation of the perturbed image with sufficient accuracy. Therefore, exploring methods on how to attack CNN-based models as well as how to defend models against attacks have become a popular topic as this also provides insights into the performance and generalization abilities of CNNs. However, most of the existing work assumes unrealistic attack models, i.e. the resulting attacks were specified in advance. In this paper, we propose a novel approach for generating adversarial examples to attack CNN-based segmentation models for medical images. Our approach has three key features: 1) The generated adversarial examples exhibit anatomical variations (in form of deformations) as well as appearance perturbations; 2) The adversarial examples attack segmentation models so that the Dice scores decrease by a pre-specified amount; 3) The attack is not required to be specified beforehand. We have evaluated our approach on CNN-based approaches for the multi-organ segmentation problem in 2D CT images. We show that the proposed approach can be used to attack different CNN-based segmentation models.


title: An interpretation of Temam's extra force in the quasi-incompressible Navier-Stokes system

abstract: We discuss the role of the extra force $$ \mathbf f_{\rm e}=-\frac 1 2(\nabla\cdot\mathbf v)\mathbf v $$ in the system of partial differential equations \begin{equation*} \left\{ \begin{aligned} &\frac{\partial\mathbf v}{\partial \rm t}+(\mathbf v\cdot\nabla)\mathbf v+\nabla \mathrm p-\frac 1 {\rm Re}\Delta\mathbf v=\mathbf f+\mathbf f_{\rm e},\\ &\frac 1 {\mathrm K}\frac{\partial \mathrm p}{\partial \rm t}+\nabla\cdot\mathbf v=0,\qquad \mathrm K>>1, \end{aligned}\right. \end{equation*} whose weak solution has been proved in [Arch. Rat. Mech. Analysis, 32:135-153] to approximate, in the limit $\mathrm K\to\infty$, the weak solution of the incompressible Navier-Stokes system. Taking the cue from [Ann. Mat. Pura Appl. 172:103-124], we provide a physical interpretation of the extra force $\mathbf f_{\rm e}$, showing that it is a manifestation of inertia.


title: Intense white luminescence in ZnTe embedded porous silicon

abstract: Porous silicon layers were embedded with ZnTe using the isothermal close space sublimation technique. The presence of ZnTe was demonstrated using cross-sectional energy dispersive spectroscopy maps. ZnTe embedded samples present intense room temperature photoluminescence along the whole visible range. We ascribe this PL to ZnTe nanocrystals of different sizes grown on the internal pore surface. Such crystals, with different orientations and sizes, were observed in transmission electron microscopy images, while transmission electron diffraction images of the same regions reveal ZnTe characteristic patterns.


title: Labels of real projective varieties

abstract: Let $X$ be a complex projective variety defined over $\mathbb R$. Recently, Bernardi and the first author introduced the notion of admissible rank with respect to $X$. This rank takes into account only decompositions that are stable under complex conjugation. Such a decomposition carries a label, i.e., a pair of integers recording the cardinality of its totally real part. We study basic properties of admissible ranks for varieties, along with special examples of curves; for instance, for rational normal curves admissible and complex ranks coincide. Along the way, we introduce the scheme theoretic version of admissible rank. Finally, analogously to the situation of real ranks, we analyze typical labels, i.e., those arising as labels of a full-dimensional Euclidean open set. We highlight similarities and differences with typical ranks.


title: Survival analysis as a classification problem

abstract: In this paper, we explore a method for treating survival analysis as a classification problem. The method uses a "stacking" idea that collects the features and outcomes of the survival data in a large data frame, and then treats it as a classification problem. In this framework, various statistical learning algorithms (including logistic regression, random forests, gradient boosting machines and neural networks) can be applied to estimate the parameters and make predictions. For stacking with logistic regression, we show that this approach is approximately equivalent to the Cox proportional hazards model with both theoretical analysis and simulation studies. For stacking with other machine learning algorithms, we show through simulation studies that our method can outperform Cox proportional hazards model in terms of estimated survival curves. This idea is not new, but we believe that it should be better known by statistiicians and other data scientists.


title: Generic uniqueness and stability for the mixed ray transform

abstract: We consider the mixed ray transform of tensor fields on a three-dimensional compact simple Riemannian manifold with boundary. We prove the injectivity of the transform, up to natural obstructions, and establish stability estimates for the normal operator on generic three dimensional simple manifold in the case of 1+1 and 2+2 tensors fields. We show how the anisotropic perturbations of averaged isotopic travel-times of qS-polarized elastic waves provide partial information about the mixed ray transform of 2+2 tensors fields. If in addition we include the measurement of the shear wave amplitude, the complete mixed ray transform can be recovered. We also show how one can obtain the mixed ray transform from an anisotropic perturbation of the Dirichlet-to-Neumann map of an isotropic elastic wave equation on a smooth and bounded domain in three dimensional Euclidean space.


title: Active Goal Recognition

abstract: To coordinate with other systems, agents must be able to determine what the systems are currently doing and predict what they will be doing in the future---plan and goal recognition. There are many methods for plan and goal recognition, but they assume a passive observer that continually monitors the target system. Real-world domains, where information gathering has a cost (e.g., moving a camera or a robot, or time taken away from another task), will often require a more active observer. We propose to combine goal recognition with other observer tasks in order to obtain \emph{active goal recognition} (AGR). We discuss this problem and provide a model and preliminary experimental results for one form of this composite problem. As expected, the results show that optimal behavior in AGR problems balance information gathering with other actions (e.g., task completion) such as to achieve all tasks jointly and efficiently. We hope that our formulation opens the door for extensive further research on this interesting and realistic problem.


title: Coherent WDM transmission using quantum-dash mode-locked laser diodes as multi-wavelength source and local oscillator

abstract: Quantum-dash (QD) mode-locked laser diodes (MLLD) lend themselves as chip-scale frequency comb generators for highly scalable wavelength-division multiplexing (WDM) links in future data-center, campus-area, or metropolitan networks. Driven by a simple DC current, the devices generate flat broadband frequency combs, containing tens of equidistant optical tones with line spacings of tens of GHz. Here we show that QD-MLLDs can not only be used as multi-wavelength light sources at a WDM transmitter, but also as multi-wavelength local oscillators (LO) for parallel coherent reception. In our experiments, we demonstrate transmission of an aggregate data rate of 4.1 Tbit/s (23x45 GBd PDM-QPSK) over 75 km standard single-mode fiber (SSMF). To the best of our knowledge, this represents the first demonstration of a coherent WDM link that relies on QD-MLLD both at the transmitter and the receiver.


title: Advances in Microwave Near-Field Imaging: Prototypes, Systems, and Applications

abstract: A near-field microwave imaging system attempts to reveal the presence of an object and/or an electrical property distribution by measuring the scattered field from many positions surrounding the object. Over the past few decades, both the hardware and software components of a near-field microwave imaging system technology have attracted interest throughout the world. Due to limitations of hardware technology (unavailability of data acquisition apparatus), experimental microwave imaging is very challenging for the pioneers. However, Probably due to the hardware cost, most of the studies (operating at a few GHz) were still focused on software only. The feasibility of using microwave approaches to image different types of objects have been tested and verified by simulations in a variety of applications. Further, work has been conducted on improving both quantitative and qualitative algorithms to improve simulated reconstruction results. Nowadays, benefitting from the hardware progress and reduction of their cost, researchers are eager to pursue real experimental validations instead of simulations, and, more unique prototypes and commercial systems have been built for various applications. These prototypes and systems are a result of years of dedicated work and it is important to review the advancements in developed prototype systems. The article will provide an overview of the many of the systems designed from different research groups throughout the world, for applications of near-field microwave imaging. The article further outlines challenges faced in current microwave near-field imaging, developmental tendencies of engineers and scientists, and the future outlook.


title: Neuromorphic Metasurface

abstract: Metasurfaces have been used to realize optical functions such as focusing and beam steering. They use sub-wavelength nanostructures to control the local amplitude and phase of light. Here we show that such control could also enable a new function of artificial neural inference. We demonstrate that metasurfaces can directly recognize objects by focusing light from an object to different spatial locations that correspond to the class of the object.


title: Using infrared gas sensors in an in-vitro dynamic gut model for detecting short-chain fatty-acids: Technical Report

abstract: Short-chain fatty acids (SCFAs), including acetate, propionate and butyrate, are organic fatty-acids that are produced when indigested carbohydrates are fermented in the colon by gut-bacteria. Butyrate is especially considered a beneficial compound in relation to gut health and the maintenance of colonic homeostasis. Little is known about butyrate production and the current measurement methods of fecal samples are not representative and there is a strong unmet need to measure SCFAs in-vivo to further understand their role and effect on colonic function in human gut health and disease. The general aim of the project is to develop a novel sensor based on infrared technology that is capable of online detecting short chain fatty acids (SCFA; e.g. acetate, propionate and butyrate) and other small biomolecules (e.g. ammonium) produced in bioreactors. Whereas currently such levels of small biomolecules are usually not detected online due to practical constraints (sampling-analysis-processing of data), reliable online detection could significantly decrease analysis costs. Further, it would allow to perform online quality control of the bioreactors and allow to e.g. detect deviating values at an early stage. Hence, one can react much faster tackling certain problems which will further improve the quality of each experiment. Using machine learning algorithms to process the sensor signals, the results from testing in a so-called artificial gut setup are described in this technical report.


title: Is the surface fluid during surface melting, a study of the surface shear modulus of solid Gallium

abstract: Melting of 3D solids is often preceded by a melting of their surface, a distinct process which begins at a temperature lower than the melting temperature. Until now, surface melting was investigated mostly by diffraction or other non-contact techniques, revealing how the surface becomes progressively disordered with temperature. We designed a method to measure an effective shear modulus of the surface, another property which should change at melting. This property was not measured before. We applied this method to surface melting of Gallium, and found that this surface shear modulus vanishes abruptly near the onset of premelting, about 9K below melting temperature.


title: Radial dose distribution and effective delta ray radius (Penumbra radius): Determination for some ions passing through water

abstract: An analytical equation for calculating radial dose of heavy ions in water is introduced by Awad et al. (Applied Radiation and Isotopes, 142 (2018) 135-142. It is simple alternative to Monte Carlo code and is a promising code, however, still needs refinement. Refinement was added through adjusting radial dose integration upper limit which gives the effective delta-ray range, rmax and ions penumbra radius in water as well. Radial dose distributions for 85 ions forming fifteen energy groups from 0.25 to 24 MeV/n were studied. By employing the effective delta-ray range, it was possible to get more consistent radial dose distribution in comparison to experimental and Monte Carlo simulation data. The corresponding LET values of those ions were estimated and compared with SRIM program. Penumbra radii for 85 ions were determined. Good description for the penumbra radii was obtained using a proposed new equation which fits experimental data as well.


title: Assessment of an Isogeometric Approach with Catmull-Clark Subdivision Surfaces using the Laplace-Beltrami Problems

abstract: An isogeometric approach for solving the Laplace-Beltrami equation on a two-dimensional manifold embedded in three-dimensional space using a Galerkin method based on Catmull-Clark subdivision surfaces is presented and assessed. The scalar-valued Laplace-Beltrami equation requires only C0 continuity and is adopted to elucidate key features and properties of the isogeometric method using Catmull-Clark subdivision surfaces. Catmull-Clark subdivision bases are used to discretise both the geometry and the physical field. A fitting method generates control meshes to approximate any given geometry with Catmull-Clark subdivision surfaces. The performance of the Catmull-Clark subdivision method is compared to the conventional finite element method. Subdivision surfaces without extraordinary vertices show the optimal convergence rate. However, extraordinary vertices introduce error, which decreases the convergence rate. A comparative study shows the effect of the number and valences of the extraordinary vertices on accuracy and convergence. An adaptive quadrature scheme is shown to reduce the error.


title: Black-body Thermal Radiative and Thermodynamic Functions of 1-Dimensional Self-Assembly Nanotubes

abstract: Assuming that 1-dimensional self-assembly nanotubes can be represented as a 1-dimensional cavity-type black-body radiator of length L and radius r in thermal equilibrium at a temperature T, analytical expressions for the thermal radiative and thermodynamic functions of the emitted black-body radiation are obtained in the finite spectral range of the electromagnetic spectrum. The total energy density, the Stefan-Boltzmann law, the number density of the photons, the Helmholtz free energy density, the entropy density, the heat capacity at constant volume, and the pressure are expressed in terms of the polylogarithm functions. In the frequency range 0.05-0.19 PHz, the thermal radiative and thermodynamic functions of the 1-dimensional black-body radiation at different temperatures are calculated. To confirm our assumption, it is necessary to measure the total energy density in the direction L in the finite frequency range 0.05-0.19 PHz, using a single nanoparticle mass spectrometer or other optical devices.


title: On speeding up an asymptotic-analysis-based homogenisation scheme for designing gradient porous structured materials using a zoning strategy

abstract: Gradient porous structured materials possess significant potential of being applied in many engineering fields. To accelerate the design process of infill graded microstructures, a novel asymptotic homogenisation topology optimisation method was proposed by Zhu et al.(2019), aiming for 1) significantly enriching the pool of representable graded microstructures; 2) deriving an homogenised formulation for stress analysis in consistency with fine-scale results. But the work is severely confined from being widely applied, mainly due to the following two reasons. Firstly, to circumvent the macroscopically pointwise computation for solving various microscopic cell problems, linearisation had to be adopted for its numerical implementation, and this significantly reduces the design freedom. Secondly, lacking of sensitive analysis, genetic algorithm was chosen for optimisation, inevitably decreasing the computational efficiency. To address these bottleneck challenging issues, a zoning scheme empowered by computational parallelism is introduced, and the sensitivity analysis associated with the new asymptotic framework is conducted. Through comparisons with fine-scale simulation results, the proposed algorithm is shown to be an effective tool for evaluating the mechanical behaviour of graded microstructures. As an optimisation tool, the mapping function takes a concise and explicit form. But its parameterisation still needs further investigation, so as to improve the solution optimality of the present approach, especially in comparison with another recently proposed method by Groen et al (2018). The high computational efficiency shown by the proposed scheme is also demonstrated with optimisation results for three-dimensional graded microstructures, which are not frequently discussed in literature, possibly because of the high computational cost generated.


title: The geometry of convergence in numerical analysis

abstract: The domains of mesh functions are strict subsets of the underlying space of continuous independent variables. Spaces of partial maps between topological spaces admit topologies which do not depend on any metric. Such topologies geometrically generalize the usual numerical analysis definitions of convergence.


title: A Note on Fuzzy Automorphism and Inner Automorphism of Groups

abstract: The fuzzification of classical set theory came into existence when Zadeh [1] laid down the concept of a fuzzy set as a generalization of a crisp set. The objective of this paper is to extend the concept of fuzzy endomorphism to fuzzy automorphism. Notions of fuzzy inner automorphism and fuzzy inner automorphism induced by a fuzzy subgroup are introduced. Finally, we obtain the fuzzy analogues of well-known results of classical group theory.


title: Level-set topology optimization considering nonlinear thermoelasticity

abstract: At elevated temperature environments, elastic structures experience a change of the stress-free state of the body that can strongly influence the optimal topology of the structure. This work presents level-set based topology optimization of structures undergoing large deformations due to thermal and mechanical loads. The nonlinear analysis model is constructed by multiplicatively decomposing thermal and mechanical effects and introducing an intermediate stress-free state between the undeformed and deformed coordinates. By incorporating the thermoelastic nonlinearity into the level-set topology optimization scheme, wider design spaces can be explored with the consideration of both mechanical and thermal loads. Four numerical examples are presented that demonstrate how temperature changes affect the optimal design of large-deforming structures. In particular, we show how optimization can manipulate the material layout in order to create a counteracting effect between thermal and mechanical loads, even up to a degree that buckling and snap-through are suppressed. Hence the consideration of large deformations in conjunction with thermoelasticity opens many new possibilities for controlling and manipulating the thermo-mechanical response via topology optimization.


title: Revisiting neutron propagation-based phase contrast imaging and tomography: use of phase retrieval to amplify the effective degree of brilliance

abstract: Propagation-based neutron phase-contrast tomography was demonstrated using the ISIS pulsed spallation source. The proof-of-concept tomogram with Paganin-type phase-retrieval filter applied exhibited an effective net boost of $23\pm 1$ in the signal-to-noise ratio as compared to an attenuation-based tomogram, implying a boost in the effective degree of neutron brilliance of over two orders of magnitude. This comparison is for phase retrieval versus conventional absorption with no additional collimation in place. Expressions are provided for the optimal phase-contrast geometry as well as conditions for the validity of the method. The underpinning theory is derived under the assumption of the sample being composed of a single material. The effective boost in brilliance may be employed to give reduced acquisition time, or may instead be used to keep exposure times fixed while improving the measured contrast.


title: ProvMark: A Provenance Expressiveness Benchmarking System

abstract: System level provenance is of widespread interest for applications such as security enforcement and information protection. However, testing the correctness or completeness of provenance capture tools is challenging and currently done manually. In some cases there is not even a clear consensus about what behavior is correct. We present an automated tool, ProvMark, that uses an existing provenance system as a black box and reliably identifies the provenance graph structure recorded for a given activity, by a reduction to subgraph isomorphism problems handled by an external solver. ProvMark is a beginning step in the much needed area of testing and comparing the expressiveness of provenance systems. We demonstrate ProvMark's usefuless in comparing three capture systems with different architectures and distinct design philosophies.


title: Towards Variable Assistance for Lower Body Exoskeletons

abstract: This paper presents and experimentally demonstrates a novel framework for variable assistance on lower body exoskeletons, based upon safety-critical control methods. Existing work has shown that providing some freedom of movement around a nominal gait, instead of rigidly following it, accelerates the spinal learning process of people with a walking impediment when using a lower body exoskeleton. With this as motivation, we present a method to accurately control how much a subject is allowed to deviate from a given gait while ensuring robustness to patient perturbation. This method leverages control barrier functions to force certain joints to remain inside predefined trajectory tubes in a minimally invasive way. The effectiveness of the method is demonstrated experimentally with able-bodied subjects and the Atalante lower body exoskeleton.


title: Diachronic Topics in New High German Poetry

abstract: Statistical topic models are increasingly and popularly used by Digital Humanities scholars to perform distant reading tasks on literary data. It allows us to estimate what people talk about. Especially Latent Dirichlet Allocation (LDA) has shown its usefulness, as it is unsupervised, robust, easy to use, scalable, and it offers interpretable results. In a preliminary study, we apply LDA to a corpus of New High German poetry (textgrid, with 51k poems, 8m token), and use the distribution of topics over documents for a classification of poems into time periods and for authorship attribution.


title: Mobile Phone Data for Children on the Move: Challenges and Opportunities

abstract: Today, 95% of the global population has 2G mobile phone coverage and the number of individuals who own a mobile phone is at an all time high. Mobile phones generate rich data on billions of people across different societal contexts and have in the last decade helped redefine how we do research and build tools to understand society. As such, mobile phone data has the potential to revolutionize how we tackle humanitarian problems, such as the many suffered by refugees all over the world. While promising, mobile phone data and the new computational approaches bring both opportunities and challenges. Mobile phone traces contain detailed information regarding people's whereabouts, social life, and even financial standing. Therefore, developing and adopting strategies that open data up to the wider humanitarian and international development community for analysis and research while simultaneously protecting the privacy of individuals is of paramount importance. Here we outline the challenging situation of children on the move and actions UNICEF is pushing in helping displaced children and youth globally, and discuss opportunities where mobile phone data can be used. We identify three key challenges: data access, data and algorithmic bias, and operationalization of research, which need to be addressed if mobile phone data is to be successfully applied in humanitarian contexts.


title: IoT Technologies for Augmented Human: a Survey

abstract: Internet of Things (IoT) technology has delivered new enablers for improving human abilities. These enablers promise an enhanced quality of life and professional efficiency; however, the synthesis of IoT and human augmentation technologies has also extended IoT-related challenges far beyond the current scope. These potential challenges associated with IoT-empowered Augmented Human (AH) have so far not been well-investigated. Thus, this article attempts to introduce readers to AH concept as well as summarize notable research challenges raised by such systems, in order to facilitate reader's further interest in this topic. The article considers emerging IoT applications for human augmentation, devices and design principles, connectivity demands, and security aspects.


title: The Bouncing Penny and Nonholonomic Impacts

abstract: The evolution of a Lagrangian mechanical system is variational. Likewise, when dealing with a hybrid Lagrangian system (a system with discontinuous impacts), the impacts can also be described by variations. These variational impacts are given by the so-called Weierstrass-Erdmann corner conditions. Therefore, hybrid Lagrangian systems can be completely understood by variational principles. Unlike typical (unconstrained / holonomic) Lagrangian systems, nonholonomically constrained Lagrangian systems are not variational. However, by using the Lagrange-d'Alembert principle, nonholonomic systems can be described as projections of variational systems. This paper works out the analogous version of the Weierstrass-Erdmann corner conditions for nonholonomic systems and examines the billiard problem with a rolling disk.


title: Scaling-Translation-Equivariant Networks with Decomposed Convolutional Filters

abstract: Encoding the scale information explicitly into the representation learned by a convolutional neural network (CNN) is beneficial for many computer vision tasks especially when dealing with multiscale inputs. We study, in this paper, a scaling-translation-equivariant (ST-equivariant) CNN with joint convolutions across the space and the scaling group, which is shown to be both sufficient and necessary to achieve equivariance for the regular representation of the scaling-translation group ST . To reduce the model complexity and computational burden, we decompose the convolutional filters under two pre-fixed separable bases and truncate the expansion to low-frequency components. A further benefit of the truncated filter expansion is the improved deformation robustness of the equivariant representation, a property which is theoretically analyzed and empirically verified. Numerical experiments demonstrate that the proposed scaling-translation-equivariant network with decomposed convolutional filters (ScDCFNet) achieves significantly improved performance in multiscale image classification and better interpretability than regular CNNs at a reduced model size.


title: Data-driven ambiguity sets with probabilistic guarantees for dynamic processes

abstract: Distributional ambiguity sets provide quantifiable ways to characterize the uncertainty about the true probability distribution of random variables of interest. This makes them a key element in data-driven robust optimization by exploiting high-confidence guarantees to hedge against uncertainty. This paper explores the construction of Wasserstein ambiguity sets in dynamic scenarios where data is collected progressively and may only reveal partial information about the unknown random variable. For random variables evolving according to known dynamics, we leverage assimilated samples to make inferences about their unknown distribution at the end of the sampling horizon. Under exact knowledge of the flow map, we provide sufficient conditions that relate the growth of the trajectories with the sampling rate to establish a reduction of the ambiguity set size as the horizon increases. Further, we characterize the exploitable sample history that results in a guaranteed reduction of ambiguity sets under errors in the computation of the flow and when the dynamics is subject to bounded unknown disturbances. Our treatment deals with both full- and partial-state measurements and, in the latter case, exploits the sampled-data observability properties of linear time-varying systems under irregular sampling. Simulations on a UAV detection application show the superior performance resulting from the proposed dynamic ambiguity sets.


title: Unconventional Hund Metal in a Weak Itinerant Ferromagnet

abstract: The physics of weak itinerant ferromagnets is challenging due to their small magnetic moments and the ambiguous role of local interactions governing their electronic properties, many of which violate Fermi liquid theory. While magnetic fluctuations play an important role in the materials' unusual electronic states, the nature of these fluctuations and the paradigms through which they arise remain debated. Here we use inelastic neutron scattering to study magnetic fluctuations in the canonical weak itinerant ferromagnet MnSi. Data reveal that short-wavelength magnons continue to propagate until a mode crossing predicted for strongly interacting quasiparticles is reached, and the local susceptibility peaks at a coherence energy predicted for a correlated Hund metal by first-principles many-body theory. Scattering between electrons and orbital and spin fluctuations in MnSi can be understood at the local level to generate non-Fermi liquid character. These results provide crucial insight into the role of interorbital Hund's exchange within the broader class of enigmatic multiband itinerant, weak ferromagnets.


title: Occurrence of semiclassical vacuum decay

abstract: Recently, a novel phenomenon is observed for vacuum decay to proceed via classically allowed dynamical evolution from initial configurations of false vacuum fluctuations. With the help of some occasionally developed large fluctuations in time derivative of a homogeneous scalar field that is initially sitting at the false vacuum, the flyover vacuum decay is proposed if the initial Gaussian profile of field velocity contains a spatial region where the field velocity is large enough to flyover the potential barrier. In this paper, we point out that, even if the initial profile of field velocity is nowhere to be large enough to classically overcome the barrier, the semiclassical vacuum decay could still be possible if the initial profile of field velocity extends over large enough region.


title: Graph-Partitioning-Based Diffusion Convolutional Recurrent Neural Network for Large-Scale Traffic Forecasting

abstract: Traffic forecasting approaches are critical to developing adaptive strategies for mobility. Traffic patterns have complex spatial and temporal dependencies that make accurate forecasting on large highway networks a challenging task. Recently, diffusion convolutional recurrent neural networks (DCRNNs) have achieved state-of-the-art results in traffic forecasting by capturing the spatiotemporal dynamics of the traffic. Despite the promising results, however, applying DCRNNs for large highway networks still remains elusive because of computational and memory bottlenecks. We present an approach for implementing a DCRNN for a large highway network that overcomes these limitations. Our approach uses a graph-partitioning method to decompose a large highway network into smaller networks and trains them independently. We demonstrate the efficacy of the graph-partitioning-based DCRNN approach to model the traffic on a large California highway network with 11,160 sensor locations. We develop an overlapping nodes approach for the graph-partitioning-based DCRNN to include sensor locations from partitions that are geographically close to a given partition. Furthermore, we demonstrate that the DCRNN model can be used to forecast the speed and flow simultaneously and that the forecasted values preserve fundamental traffic flow dynamics. Our approach to developing DCRNN models that represent large highway networks can be a potential core capability in advanced highway traffic monitoring systems, where a trained DCRNN model forecasting traffic at all sensor locations can be used to adjust traffic management strategies proactively based on anticipated future conditions.


title: Constraints on Sterile Neutrinos in the MeV to GeV Mass Range

abstract: A detailed discussion is given of the analysis of recent data to obtain improved upper bounds on the couplings $|U_{e4}|^2$ and $|U_{\mu 4}|^2$ for a mainly sterile neutrino mass eigenstate $\nu_4$. Using the excellent agreement among ${\cal F}t$ values for superallowed nuclear beta decay, an improved upper limit is derived for emission of a $\nu_4$. The agreement of the ratios of branching ratios $R^{(\pi)}_{e/\mu}=BR(\pi^+ \to e^+ \nu_e)/BR(\pi^+ \to \mu^+ \nu_\mu)$, $R^{(K)}_{e/\mu}$, $R^{(D_s)}_{e/\tau}$, $R^{(D_s)}_{\mu/\tau}$, and $R^{(D)}_{e/\tau}$, and the branching ratios $BR(B^+\rightarrow e^+\nu_e)$ and $BR(B^+\rightarrow \mu^+\nu_\mu)$ decays with predictions of the Standard Model, is utilized to derive new constraints on $\nu_4$ emission covering the $\nu_4$ mass range from MeV to GeV. We also discuss constraints from peak search experiments probing for emission of a $\nu_4$ via lepton mixing, as well as constraints from pion beta decay, CKM unitarity, $\mu$ decay, leptonic $\tau$ decay, and other experimental inputs.


title: Security Risk Analysis of the Shorter-Queue Routing Policy for Two Symmetric Servers

abstract: In this article, we study the classical shortest queue problem under the influence of malicious attacks, which is relevant to a variety of engineering system including transportation, manufacturing, and communications. We consider a homogeneous Poisson arrival process of jobs and two parallel exponential servers with symmetric service rates. A system operator route incoming jobs to the shorter queue; if the queues are equal, the job is routed randomly. A malicious attacker is able to intercept the operator's routing instruction and overwrite it with a randomly generated one. The operator is able to defend individual jobs to ensure correct routing. Both attacking and defending induce technological costs. The attacker's (resp. operator's) decision is the probability of attacking (resp. defending) the routing of each job. We first quantify the queuing cost for given strategy profiles by deriving a theoretical upper bound for the cost. Then, we formulate a non-zero-sum attacker-defender game, characterize the equilibria in multiple regimes, and quantify the security risk. We find that the attacker's best strategy is either to attack all jobs or not to attack, and the defender's strategy is strongly influenced by the arrival rate of jobs. Finally, as a benchmark, we compare the security risks of the feedback-controlled system to a corresponding open-loop system with Bernoulli routing. We show that the shorter-queue policy has a higher (resp. lower) security risk than the Bernoulli policy if the demand is lower (resp. higher) than the service rate of one server.


title: Improving Noise Robustness In Speaker Identification Using A Two-Stage Attention Model

abstract: While the use of deep neural networks has significantly boosted speaker recognition performance, it is still challenging to separate speakers in poor acoustic environments. To improve robustness of speaker recognition system performance in noise, a novel two-stage attention mechanism which can be used in existing architectures such as Time Delay Neural Networks (TDNNs) and Convolutional Neural Networks (CNNs) is proposed. Noise is known to often mask important information in both time and frequency domain. The proposed mechanism allows the models to concentrate on reliable time/frequency components of the signal. The proposed approach is evaluated using the Voxceleb1 dataset, which aims at assessment of speaker recognition in real world situations. In addition three types of noise at different signal-noise-ratios (SNRs) were added for this work. The proposed mechanism is compared with three strong baselines: X-vectors, Attentive X-vector, and Resnet-34. Results on both identification and verification tasks show that the two-stage attention mechanism consistently improves upon these for all noise conditions.


title: Matrix Sketching for Secure Collaborative Machine Learning

abstract: Collaborative learning allows participants to jointly train a model without data sharing. To update the model parameters, the central server broadcasts model parameters to the clients, and the clients send updating directions such as gradients to the server. While data do not leave a client device, the communicated gradients and parameters will leak a client's privacy. Attacks that infer clients' privacy from gradients and parameters have been developed by prior work. Simple defenses such as dropout and differential privacy either fail to defend the attacks or seriously hurt test accuracy. We propose a practical defense which we call Double-Blind Collaborative Learning (DBCL). The high-level idea is to apply random matrix sketching to the parameters (aka weights) and re-generate random sketching after each iteration. DBCL prevents clients from conducting gradient-based privacy inferences which are the most effective attacks. DBCL works because from the attacker's perspective, sketching is effectively random noise that outweighs the signal. Notably, DBCL does not much increase computation and communication costs and does not hurt test accuracy at all.


title: A Visual Analytics Framework for Adversarial Text Generation

abstract: This paper presents a framework which enables a user to more easily make corrections to adversarial texts. While attack algorithms have been demonstrated to automatically build adversaries, changes made by the algorithms can often have poor semantics or syntax. Our framework is designed to facilitate human intervention by aiding users in making corrections. The framework extends existing attack algorithms to work within an evolutionary attack process paired with a visual analytics loop. Using an interactive dashboard a user is able to review the generation process in real time and receive suggestions from the system for edits to be made. The adversaries can be used to both diagnose robustness issues within a single classifier or to compare various classifier options. With the weaknesses identified, the framework can also be used as a first step in mitigating adversarial threats. The framework can be used as part of further research into defense methods in which the adversarial examples are used to evaluate new countermeasures. We demonstrate the framework with a word swapping attack for the task of sentiment classification.


title: Asynchronous and time-varying proximal type dynamics multi-agent network games

abstract: In this paper, we study proximal type dynamics in the context of noncooperative multi-agent network games. These dynamics arise in different applications, since they describe distributed decision making in multi-agent networks, e.g., in opinion dynamics, distributed model fitting and network information fusion, where the goal of each agent is to seek an equilibrium using local information only. We analyse several conjugations of this class of games, providing convergence results, or designing equilibrium seeking algorithms when the original dynamics fail to converge. For the games subject only to local constraints we look into both synchronous/asynchronous dynamics and time-varying communication networks. For games subject in addition to coupling constraints, we design an equilibrium seeking algorithm converging to a special class of game equilibria. Finally, we validate the theoretical results via numerical simulations on opinion dynamics and distributed model fitting.


title: Automatic Snake Gait Generation Using Model Predictive Control

abstract: In this paper, we propose a method for generating undulatory gaits for snake robots. Instead of starting from a pre-defined movement pattern such as a serpenoid curve, we use a Model Predictive Control approach to automatically generate effective locomotion gaits via trajectory optimization. An important advantage of this approach is that the resulting gaits are automatically adapted to the environment that is being modeled as part of the snake dynamics. To illustrate this, we use a novel model for anisotropic dry friction, along with existing models for viscous friction and fluid dynamic effects such as drag and added mass. For each of these models, gaits generated without any change in the method or its parameters are as efficient as Pareto-optimal serpenoid gaits tuned individually for each environment. Furthermore, the proposed method can also produce more complex or irregular gaits, e.g. for obstacle avoidance or executing sharp turns.


title: Photon--Matter Quantum Correlations in Spontaneous Raman Scattering

abstract: We develop a Hamiltonian formalism to study energy and position/momentum correlations between a single Stokes photon and a single material excitation that are created as a pair in the spontaneous Raman scattering process. Our approach allows for intuitive separation of the effects of spectral linewidth, chromatic dispersion, and collection angle on these correlations, and we compare the predictions of the model to experiment. These results have important implications for the use of Raman scattering in quantum protocols that rely on spectrally unentangled photons and collective excitations.


title: Using human-in-the-loop synthesis to author functional reactive programs

abstract: Programs that respond to asynchronous events are challenging to write; they are difficult to reason about and tricky to test and debug. Because these programs can have a huge space of possible input timings and interleaving, the programmer may easily miss corner cases. We propose applying synthesis to aid programmers in creating programs more easily and with a higher degree of confidence in their correctness. We have written an efficient encoding of functional reactive programming (FRP) semantics based on functional programming over lists lifted in Rosette. We demonstrate that this technique is state-of-the-art by first comparing its performance against two existing synthesis tools that produce list manipulation programs, and then by synthesizing a suite of benchmarks given complete specifications. We also propose an interactive tool in which a programmer provides some initial partial specification in the form of input/output examples or invariants; the tool finds ambiguity in the specification by synthesizing two candidate programs and gives the user an input that distinguishes them; the user updates the specification and continues iterating until the correct program is found. As evaluation, we demonstrate the use of the tool on a suite of benchmarks from the web programming and Internet of Things domains and walk through a sample interaction on a realistic web benchmark, showing that we can converge on the target program with a tractable number of interactions. As future work, we discuss encoding additional FRP languages to in order to explore metalinguistic features, strategies for decomposition that would allow the synthesis of larger programs, and improved programmer tools such as a GUI to more easily elicit specifications.


title: Simple and Almost Assumption-Free Out-of-Sample Bound for Random Feature Mapping

abstract: Random feature mapping (RFM) is a popular method for speeding up kernel methods at the cost of losing a little accuracy. We study kernel ridge regression with random feature mapping (RFM-KRR) and establish novel out-of-sample error upper and lower bounds. While out-of-sample bounds for RFM-KRR have been established by prior work, this paper's theories are highly interesting for two reasons. On the one hand, our theories are based on weak and valid assumptions. In contrast, the existing theories are based on various uncheckable assumptions, which makes it unclear whether their bounds are the nature of RFM-KRR or simply the consequence of strong assumptions. On the other hand, our analysis is completely based on elementary linear algebra and thereby easy to read and verify. Finally, our experiments lend empirical supports to the theories.


title: The Kauffman skein algebra of the torus

abstract: We give a presentation of the Kauffman (BMW) skein algebra of the torus, which is the "type BCD" analogue of the Homflypt skein algebra of torus which was computed by the first and third authors. In the appendix we show this presentation is compatible with the Frohman-Gelca description of the Kauffman bracket (Temperley-Lieb) skein algebra of the torus [FG00].


title: Vertical transport and tunnelling in rare-earth nitride heterostructures

abstract: We report an investigation of the ferromagnetic semiconductor rare earth nitrides (RENs) for their potential for cryogenic-temperature electronics and spintronics application. We have indentified ohmic contacts suitable for the device structures that demand electron transport through interface layers, and grown REN/insulator/REN heterostructures that display tunnelling characteristics, an enormous 400% tunneling magnetoresistance and a hysteresis promising their exploitation in non-volatile magnetic random access memory.


title: Gene ologs: a categorical framework for Gene Ontology

abstract: Gene Ontology (GO) is the most important resource for gene function annotation. It provides a way to unify biological knowledge across different species via a dynamic and controlled vocabulary. GO is now widely represented in the Semantic Web standard Web Ontology Language (OWL). OWL renders a rich logic constructs to GO but also has its limitations. On the other hand, olog is a different language for ontology and it is based on category theory. Due to its solid mathematical background, olog can be rigorously formulated yet considerably expressive. These advantages make ologs complementary, if not better than, OWL. We therefore think it worthwhile to adopt ologs for GO and took an initial step in this work.


title: Characterizing physiological and symptomatic variation in menstrual cycles using self-tracked mobile health data

abstract: The menstrual cycle is a key indicator of overall health for women of reproductive age. Previously, menstruation was primarily studied through survey results; however, as menstrual tracking mobile apps become more widely adopted, they provide an increasingly large, content-rich source of menstrual health experiences and behaviors over time. By exploring a database of user-tracked observations from the Clue app by BioWink of over 378,000 users and 4.9 million natural cycles, we show that self-reported menstrual tracker data can reveal statistically significant relationships between per-person cycle length variability and self-reported qualitative symptoms. A concern for self-tracked data is that they reflect not only physiological behaviors, but also the engagement dynamics of app users. To mitigate such potential artifacts, we develop a procedure to exclude cycles lacking user engagement, thereby allowing us to better distinguish true menstrual patterns from tracking anomalies. We uncover that women located at different ends of the menstrual variability spectrum, based on the consistency of their cycle length statistics, exhibit statistically significant differences in their cycle characteristics and symptom tracking patterns. We also find that cycle and period length statistics are stationary over the app usage timeline across the variability spectrum. The symptoms that we identify as showing statistically significant association with timing data can be useful to clinicians and users for predicting cycle variability from symptoms or as potential health indicators for conditions like endometriosis. Our findings showcase the potential of longitudinal, high-resolution self-tracked data to improve understanding of menstruation and women's health as a whole.


title: Augmenting the Pathology Lab: An Intelligent Whole Slide Image Classification System for the Real World

abstract: Standard of care diagnostic procedure for suspected skin cancer is microscopic examination of hematoxylin \& eosin stained tissue by a pathologist. Areas of high inter-pathologist discordance and rising biopsy rates necessitate higher efficiency and diagnostic reproducibility. We present and validate a deep learning system which classifies digitized dermatopathology slides into 4 categories. The system is developed using 5,070 images from a single lab, and tested on an uncurated set of 13,537 images from 3 test labs, using whole slide scanners manufactured by 3 different vendors. The system's use of deep-learning-based confidence scoring as a criterion to consider the result as accurate yields an accuracy of up to 98\%, and makes it adoptable in a real-world setting. Without confidence scoring, the system achieved an accuracy of 78\%. We anticipate that our deep learning system will serve as a foundation enabling faster diagnosis of skin cancer, identification of cases for specialist review, and targeted diagnostic classifications.


title: Probabilistic Data Association via Mixture Models for Robust Semantic SLAM

abstract: Modern robotic systems sense the environment geometrically, through sensors like cameras, lidar, and sonar, as well as semantically, often through visual models learned from data, such as object detectors. We aim to develop robots that can use all of these sources of information for reliable navigation, but each is corrupted by noise. Rather than assume that object detection will eventually achieve near perfect performance across the lifetime of a robot, in this work we represent and cope with the semantic and geometric uncertainty inherent in methods like object detection. Specifically, we model data association ambiguity, which is typically non-Gaussian, in a way that is amenable to solution within the common nonlinear Gaussian formulation of simultaneous localization and mapping (SLAM). We do so by eliminating data association variables from the inference process through max-marginalization, preserving standard Gaussian posterior assumptions. The result is a max-mixture-type model that accounts for multiple data association hypotheses as well as incorrect loop closures. We provide experimental results on indoor and outdoor semantic navigation tasks with noisy odometry and object detection and find that the ability of the proposed approach to represent multiple hypotheses, including the "null" hypothesis, gives substantial robustness advantages in comparison to alternative semantic SLAM approaches.


title: Periodic continued fractions over $S$-integers in number fields and Skolem's $p$-adic method

abstract: We generalize the classical theory of periodic continued fractions (PCFs) over ${\mathbf Z}$ to rings ${\mathcal O}$ of $S$-integers in a number field. Let ${\mathcal B}=\{\beta, {\beta^*}\}$ be the multi-set of roots of a quadratic polynomial in ${\mathcal O}[x]$. We show that PCFs $P=[b_1,\ldots,b_N,\bar{a_1\ldots ,a_k}]$ of type $(N,k)$ potentially converging to a limit in ${\mathcal B}$ are given by ${\mathcal O}$-points on an affine variety $V:=V({\mathcal B})_{N,k}$ generically of dimension $N+k-2$. We give the equations of $V$ in terms of the continuant polynomials of Wallis and Euler. The integral points $V({\mathcal O})$ are related to writing matrices in $\textrm{SL}_2({\mathcal O})$ as products of elementary matrices. We give an algorithm to determine if a PCF converges and, if so, to compute its limit. Our standard example generalizes the PCF $\sqrt{2}=[1,\bar{2}]$ to the ${\mathbf Z}_2$-extension of ${\mathbf Q}$: $F_n={\mathbf Q}(\alpha_n)$, $\alpha_{n}:=2\cos(2\pi/2^{n+2})$, with integers ${\mathcal O}_n={\mathbf Z}[\alpha_n]$. We want to find the PCFs of $\alpha_{n+1}$ over ${\mathcal O}_{n}$ of type $(N,k)$ by finding the ${\mathcal O}_{n}$-points on $V({\mathcal B}_{n+1})_{N,k}$ for ${\mathcal B}_{n+1}:=\{\alpha_{n+1}, -\alpha_{n+1}\}$. There are three types $(N,k)=(0,3), (1,2), (2,1)$ such that the associated PCF variety $V({\mathcal B})_{N,k}$ is a curve; we analyze these curves. For generic ${\mathcal B}$, Siegel's theorem implies that each of these three $V({\mathcal B})_{N,k}({\mathcal O})$ is finite. We find all the ${\mathcal O}_n$-points on these PCF curves $V({\mathcal B}_{n+1})_{N,k}$ for $n=0,1$. When $n=1$ we make extensive use of Skolem's $p$-adic method for $p=2$, including its application to Ljunggren's equation $x^2 + 1 =2y^4$.


title: Braid group actions for quantum symmetric pairs of type AIII/AIV

abstract: In the present paper we construct braid group actions on quantum symmetric pair coideal subalgebras of type AIII/AIV. This completes the proof of a conjecture by Kolb and Pellegrini in the case where the underlying Lie algebra is $\mathfrak{sl}_n$. The braid group actions are defined on the generators of the coideal subalgebras and the defining relations and braid relations are verified by explicit calculations.


title: Efficient Algorithms for Distributionally Robust Stochastic Optimization with Discrete Scenario Support

abstract: Recently, there has been a growing interest in distributionally robust optimization (DRO) as a principled approach to data-driven decision making. In this paper, we consider a distributionally robust two-stage stochastic optimization problem with discrete scenario support. While much research effort has been devoted to tractable reformulations for DRO problems, especially those with continuous scenario support, few efficient numerical algorithms were developed, and most of them can neither handle the non-smooth second-stage cost function nor the large number of scenarios $K$ effectively. We fill the gap by reformulating the DRO problem as a trilinear min-max-max saddle point problem and developing novel algorithms that can achieve an $\mathcal{O}(1/\epsilon)$ iteration complexity which only mildly depends on $K$. The major computations involved in each iteration of these algorithms can be conducted in parallel if necessary. Besides, for solving an important class of DRO problems with the Kantorovich ball ambiguity set, we propose a slight modification of our algorithms to avoid the expensive computation of the probability vector projection at the price of an $\mathcal{O}(\sqrt{K})$ times more iterations. Finally, preliminary numerical experiments are conducted to demonstrate the empirical advantages of the proposed algorithms.


title: On Transformations of Markov Chains and Poisson Boundary

abstract: A discrete-time Markov chain can be transformed into a new Markov chain by looking at its states along iterations of an almost surely finite stopping time. By the optional stopping theorem, any bounded harmonic function with respect to the transition function of the original chain is harmonic with respect to the transition function of the transformed chain. The reverse inclusion is in general not true. Our main result provides a sufficient condition on the stopping time which guarantees that the space of bounded harmonic functions for the transformed chain embeds in the space of bounded harmonic sequences for the original chain. We also obtain a similar result on positive unbounded harmonic functions, under some additional conditions. Our work was motivated by and is analogous to Forghani-Kaimanovich, the well-studied case when the Markov chain is a random walk on a discrete group.


title: Attention Interpretability Across NLP Tasks

abstract: The attention layer in a neural network model provides insights into the model's reasoning behind its prediction, which are usually criticized for being opaque. Recently, seemingly contradictory viewpoints have emerged about the interpretability of attention weights (Jain & Wallace, 2019; Vig & Belinkov, 2019). Amid such confusion arises the need to understand attention mechanism more systematically. In this work, we attempt to fill this gap by giving a comprehensive explanation which justifies both kinds of observations (i.e., when is attention interpretable and when it is not). Through a series of experiments on diverse NLP tasks, we validate our observations and reinforce our claim of interpretability of attention through manual evaluation.