title: Counting and Enumerating Independent Sets with Applications to Knapsack Problems

abstract: We introduce methods to count and enumerate all maximal independent, all maximum independent sets, and all independent sets in threshold graphs and k-threshold graphs. Within threshold graphs and k-threshold graphs independent sets correspond to feasible solutions in related knapsack instances. We give several characterizations for knapsack instances and multidimensional knapsack instances which allow an equivalent graph. This allows us to solve special knapsack instances as well as special multidimensional knapsack instances for fixed number of dimensions in polynomial time. We also conclude lower bounds on the number of necessary bins within several bin packing problems.


title: Sieve-SDP: a simple facial reduction algorithm to preprocess semidefinite programs

abstract: We introduce Sieve-SDP, a simple facial reduction algorithm to preprocess semidefinite programs (SDPs). Sieve-SDP inspects the constraints of the problem to detect lack of strict feasibility, deletes redundant rows and columns, and reduces the size of the variable matrix. It often detects infeasibility. It does not rely on any optimization solver: the only subroutine it needs is Cholesky factorization, hence it can be implemented in a few lines of code in machine precision. We present extensive computational results on several problem collections from the literature, with many SDPs coming from polynomial optimization.


title: Sign patterns that require $\mathbb{H}_n$ exist for each $n\geq 4$

abstract: The refined inertia of a square real matrix $A$ is the ordered $4$-tuple $(n_+, n_-, n_z, 2n_p)$, where $n_+$ (resp., $n_-$) is the number of eigenvalues of $A$ with positive (resp., negative) real part, $n_z$ is the number of zero eigenvalues of $A$, and $2n_p$ is the number of nonzero pure imaginary eigenvalues of $A$. For $n \geq 3$, the set of refined inertias $\mathbb{H}_n=\{(0, n, 0, 0), (0, n-2, 0, 2), (2, n-2, 0, 0)\}$ is important for the onset of Hopf bifurcation in dynamical systems. We say that an $n\times n$ sign pattern ${\cal A}$ requires $\mathbb{H}_n$ if $\mathbb{H}_n=\{\text{ri}(B) | B \in Q({\cal A})\}$. Bodine et al. conjectured that no $n\times n$ irreducible sign pattern that requires $\mathbb{H}_n$ exists for $n$ sufficiently large, possibly $n\ge 8$. However, for each $n \geq 4$, we identify three $n\times n$ irreducible sign patterns that require $\mathbb{H}_n$, which resolves this conjecture.


title: Equilibrium properties of magnetic filament suspensions

abstract: Langevin dynamics is used to study equilibrium properties of the suspension of magnetic filaments (chains of nanoparticles permanently crosslinked with polymers). It is shown that the filament suspension generally has a larger magnetic susceptibility than the system of unlinked nanoparticles with the same average particle concentration. However, actual susceptibility gain strongly depends on length and flexibility of filaments. It is also shown that in a strong gravitational (centrifugal) field sedimentation profiles of filaments are less homogeneous than that of unlinked particles. The spatial distribution of filaments weakly depends on the intensity of interparticle dipole-dipole interactions.


title: Properties of Short Gamma-Ray Burst Pulses from A BATSE TTE GRB Pulse Catalog

abstract: We analyze pulse properties of Short gamma-ray bursts (GRBs) from a new catalog containing 434 pulses from 387 BATSE Time-Tagged Event (TTE) GRBs. Short GRB pulses exhibit correlated properties of duration, fluence, hardness, and amplitude, and they evolve hard-to-soft while undergoing similar triple-peaked light curves similar to those found in Long/Intermediate bursts. We classify pulse light curves using their temporal complexities, demonstrating that Short GRB pulses exhibit a range of complexities from smooth to highly variable. Most of the bright, hard, chaotic emission seen in complex pulses seems to represent a separate highly-variable emission component. Unlike Long/Intermediate bursts, as many as 90\% of Short GRBs are single-pulsed. However, emission in Short multi-pulsed bursts is coupled such that the first pulse's duration is a predictor of both the interpulse separation and subsequent pulse durations. These results strongly support the idea that external shocks produce the prompt emission seen in Short GRBs. The similarities between the triple-peaked structures and spectral evolution of Long, Short, and Intermediate GRBs then suggests that external shocks are responsible for the prompt emission observed in all GRB classes. In addition to these findings, we identify a new type of gamma-ray transient in which peak amplitudes occur at the end of the burst rather than at earlier times. Some of these "Crescendo" bursts are preceded by rapid-fire "Staccato" pulses, whereas the remaining are preceded by a variable episode that could be unresolved staccato pulses.


title: Tubes, Mono Jets, Squeeze Out and CME

abstract: Glasma Flux Tubes, Mono Jets with squeeze out flow around them plus the Chiral Magnetic Effect(CME) are physical phenomenon that generate two particle correlation with respect to the reaction plane in mid-central 20\% to 30\% Au-Au collision $\sqrt{s_{NN}}$ = 200.0 GeV measured at RHIC.


title: Decay Estimates for Solutions of Porous Medium Equations with Advection

abstract: In this paper, we show that bounded weak solutions of the Cauchy problem for general degenerate parabolic equations of the form \begin{equation} \notag u_t \,+\; \mbox{div}\,f(x,t,u) \;=\; \mbox{div}\,(\;\!|\,u\,|^{\alpha} \, \nabla u \;\!), \quad \;\; x \in \mathbb{R}^{n}\!\:\!, \; t > 0, \end{equation} where $ \alpha > 0 \, $ is constant, decrease to zero, under fairly broad conditions on the advection flux $f$. Besides that, we derive a time decay rate for these solutions.


title: Quantum-electrodynamical approach to the exciton spectrum in Transition-Metal Dichalcogenides

abstract: Manipulation of intrinsic electron degrees of freedom, such as charge and spin, gives rise to electronics and spintronics, respectively. Electrons in monolayer materials with a honeycomb lattice structure, such as the Transition-Metal Dichalcogenides (TMD's), can be distinguished according to the region (valley) of the Brillouin zone to which they belong. Valleytronics, the manipulation of this electron's property, is expected to set up a new era in the realm of electronic devices. In this work, we accurately determine the energy spectrum and lifetimes of exciton (electron-hole) bound-states for different TMD materials, namely WSe$_2$, WS$_2$ and MoS$_2$. For all of them, we obtain a splitting of the order of 170 meV between the exciton energies from different valleys, corresponding to an effective Zeeman magnetic field of 1400 T. Our approach, which employs quantum-field theory (QFT) techniques based on the Bethe-Salpeter equation and the Schwinger-Dyson formalism, takes into account the full electromagnetic interaction among the electrons. The valley selection mechanism operates through the dynamical breakdown of the time-reversal (TR) symmetry, which originally interconnects the two valleys. This symmetry is spontaneously broken whenever the full electromagnetic interaction vertex is used to probe the response of the system to an external field.


title: Fast and Scalable Distributed Deep Convolutional Autoencoder for fMRI Big Data Analytics

abstract: In recent years, analyzing task-based fMRI (tfMRI) data has become an essential tool for understanding brain function and networks. However, due to the sheer size of tfMRI data, its intrinsic complex structure, and lack of ground truth of underlying neural activities, modeling tfMRI data is hard and challenging. Previously proposed data-modeling methods including Independent Component Analysis (ICA) and Sparse Dictionary Learning only provided a weakly established model based on blind source separation under the strong assumption that original fMRI signals could be linearly decomposed into time series components with corresponding spatial maps. Meanwhile, analyzing and learning a large amount of tfMRI data from a variety of subjects has been shown to be very demanding but yet challenging even with technological advances in computational hardware. Given the Convolutional Neural Network (CNN), a robust method for learning high-level abstractions from low-level data such as tfMRI time series, in this work we propose a fast and scalable novel framework for distributed deep Convolutional Autoencoder model. This model aims to both learn the complex hierarchical structure of the tfMRI data and to leverage the processing power of multiple GPUs in a distributed fashion. To implement such a model, we have created an enhanced processing pipeline on the top of Apache Spark and Tensorflow library, leveraging from a very large cluster of GPU machines. Experimental data from applying the model on the Human Connectome Project (HCP) show that the proposed model is efficient and scalable toward tfMRI big data analytics, thus enabling data-driven extraction of hierarchical neuroscientific information from massive fMRI big data in the future.


title: Some characterizations of the preimage of $A_{\infty}$ for the Hardy-Littlewood maximal operator and consequences

abstract: The purpose of this paper is to give some characterizations of the weight functions $w$ such that $Mw$ is in $A_{\infty}$. We show that for those weights to be in $A_{\infty}$ ensures to be in $A_{1}$. We give a criterion in terms of the local maximal functions $m_{\lambda}$ and we present a pair of applications, among them someone similar to the Coifman-Rochberg characterization of $A_{1}$ but using functions of the form $(f^{\#})^{\delta}$ and $(m_{\lambda}u)^{\delta}$ instead of $(Mf)^{\delta}$.


title: Scaling Text with the Class Affinity Model

abstract: Probabilistic methods for classifying text form a rich tradition in machine learning and natural language processing. For many important problems, however, class prediction is uninteresting because the class is known, and instead the focus shifts to estimating latent quantities related to the text, such as affect or ideology. We focus on one such problem of interest, estimating the ideological positions of 55 Irish legislators in the 1991 D\'ail confidence vote. To solve the D\'ail scaling problem and others like it, we develop a text modeling framework that allows actors to take latent positions on a "gray" spectrum between "black" and "white" polar opposites. We are able to validate results from this model by measuring the influences exhibited by individual words, and we are able to quantify the uncertainty in the scaling estimates by using a sentence-level block bootstrap. Applying our method to the D\'ail debate, we are able to scale the legislators between extreme pro-government and pro-opposition in a way that reveals nuances in their speeches not captured by their votes or party affiliations.


title: Sensitivity of DUNE to long-baseline neutrino oscillation physics

abstract: DUNE is an international project, currently in its design phase, for neutrino physics and proton-decay searches. It will consist of two detectors exposed to a megawatt-scale muon neutrino beam that will be produced at Fermilab (Illinois, USA). One detector will record particle interactions near the source of the beam, while a second, much larger, detector comprising four 10-kilotonne liquid argon TPCs will be installed at a depth of 1.5 km at SURF (South Dakota, USA), about 1300 kilometres away of the neutrino source. Among the primary scientific goals of DUNE is the precision measurement of the parameters that govern neutrino mixing, including those still unknown: the octant in which the $\theta_{23}$ mixing angle lies, the neutrino mass ordering and the value of the CP violation phase. This paper discusses the sensitivity of the DUNE experiment to these parameters.


title: Van der Waals interaction potential between Rydberg atoms near surfaces

abstract: Van der Waals interactions, as a result of the exchange of photons between particles, can be altered by modifying the environment through which these photons propagate. As a consequence, phenomena such as the Rydberg blockade mechanism between highly excited atoms or excitons, can be controlled by adding reflecting surfaces. We provide the quantum electrodynamic framework for the van der Waals interaction in the nonretarded limit that is relevant for long-wavelength transitions such as those between Rydberg systems, and show its intimate connection with common static dipole-dipole interactions, thereby providing a generalization to include macroscopic bodies.


title: Functors and morphisms determined by subcategories

abstract: We study the existence and uniqueness of minimal right determiners in various categories. Particularly in a Hom-finite hereditary abelian category with enough projectives, we prove that the Auslander-Reiten-Smal{\o}-Ringel formula of the minimal right determiner still holds. As an application, we give a formula of minimal right determiners in the category of finitely presented representations of strongly locally finite quivers.


title: Thermal management and non-reciprocal control of phonon flow via optomechanics

abstract: Engineering phonon transport in physical systems is a subject of interest in the study of materials and plays a crucial role in controlling energy and heat transfer. Of particular interest are non-reciprocal phononic systems, which in direct analogy to electric diodes, provide a directional flow of energy. Here, we propose an engineered nanostructured material, in which tunable non-reciprocal phonon transport is achieved through optomechanical coupling. Our scheme relies on breaking time-reversal symmetry by a spatially varying laser drive, which manipulates low-energy acoustic phonons. Furthermore, we take advantage of recent developments in the manipulation of high-energy phonons through controlled scattering mechanisms, such as using alloys and introducing disorder. These combined approaches allow us to design an acoustic isolator and a thermal diode. Our proposed device will have potential impact in phonon-based information processing, and heat management in low temperatures.


title: On the Shapes of Interstellar Grains: Modeling Extinction and Polarization by Spheroids and Continuous Distributions of Ellipsoids

abstract: Although interstellar grains are known to be aspherical, their actual shapes remain poorly constrained. We assess whether three continuous distributions of ellipsoidal shapes from the literature are suitable for describing the shapes of interstellar grains. Randomly-selected shapes from each distribution are shown as illustrations. The often-used Bohren-Huffman CDE includes a very large fraction of extreme shapes: fully 10% of random draws have axial ratio $a_3/a_1 > 19.7$, and 5% have $a_3/a_1 > 33$. The CDE2 distribution includes a much smaller fraction of extreme shapes, and appears to be more realistic. For each of the three CDEs considered, we derive shape-averaged cross sections for extinction and polarization in the Rayleigh limit. Finally, we describe a method for "synthesizing" a dielectric function for an assumed shape or shape distribution if the actual absorption cross sections per grain volume in the Rayleigh limit are known from observations. This synthetic dielectric function predicts the wavelength dependence of polarization, which can then be compared to observations to constrain the grain shape.


title: Efficiently Trainable Text-to-Speech System Based on Deep Convolutional Networks with Guided Attention

abstract: This paper describes a novel text-to-speech (TTS) technique based on deep convolutional neural networks (CNN), without use of any recurrent units. Recurrent neural networks (RNN) have become a standard technique to model sequential data recently, and this technique has been used in some cutting-edge neural TTS techniques. However, training RNN components often requires a very powerful computer, or a very long time, typically several days or weeks. Recent other studies, on the other hand, have shown that CNN-based sequence synthesis can be much faster than RNN-based techniques, because of high parallelizability. The objective of this paper is to show that an alternative neural TTS based only on CNN alleviate these economic costs of training. In our experiment, the proposed Deep Convolutional TTS was sufficiently trained overnight (15 hours), using an ordinary gaming PC equipped with two GPUs, while the quality of the synthesized speech was almost acceptable.


title: Investigating BSM Models with Large Scale Separation

abstract: Mass-split systems based on a conformal infrared fixed point provide a low-energy effective description of beyond the standard model systems with large scale separation. We report results of exploratory investigations with four light and eight heavy flavors using staggered fermions, and up to five different values for the light flavor mass, five different heavy flavor masses, and two values of the bare gauge coupling.


title: A characterization of the $n$-ary many-sorted closure operators and a many-sorted Tarski irredundant basis theorem

abstract: A theorem of single-sorted algebra states that, for a closure space $(A,J)$ and a natural number $n$, the closure operator $J$ on the set $A$ is $n$-ary if, and only if, there exists a single-sorted signature $\Sigma$ and a $\Sigma$-algebra $\mathbf{A}$ such that every operation of $\mathbf{A}$ is of an arity $\leq n$ and $J = \mathrm{Sg}_{\mathbf{A}}$, where $\mathrm{Sg}_{\mathbf{A}}$ is the subalgebra generating operator on $A$ determined by $\mathbf{A}$. On the other hand, a theorem of Tarski asserts that if $J$ is an $n$-ary closure operator on a set $A$ with $n\geq 2$, and if $i<j$ with $i$, $j\in \mathrm{IrB}(A,J)$, where $\mathrm{IrB}(A,J)$ is the set of all natural numbers $n$ such that $(A,J)$ has an irredundant basis ($\equiv$ minimal generating set) of $n$ elements, such that $\{i+1,\ldots, j-1\}\cap \mathrm{IrB}(A,J) = \varnothing$, then $j-i\leq n-1$. In this article we state and prove the many-sorted counterparts of the above theorems. But, we remark, regarding the first one under an additional condition: the uniformity of the many-sorted closure operator.


title: On the Importance of Engaging Students in Crafting Definitions

abstract: In this paper we describe an activity for engaging students in crafting definitions. We explore the strengths of this particular activity as well as the broader implications of engaging students in crafting definitions more generally.


title: A unified treatment of derivative discontinuity, delocalization and static correlation effects in density functional calculations

abstract: We propose a method that incorporates explicit derivative discontinuity of the total energy with respect to the number of electrons and treats both delocalization and static correlation effects in density functional calculations. Our approach is motivated by the exact behavior of the ground state total energy of electrons and involves minimization of the exchange-correlation energy with respect to the Fock space density matrix. The resulting density matrix minimization (DMM) model is simple to implement and can be solved uniquely and efficiently. In a case study of KCuF$_3$, a prototypical Mott-insulator with strong correlation, LDA+DMM correctly reproduced the Mott-Hubbard gap, magnetic ordering and Jahn-Teller distortion.


title: Strict convexity of the free energy of the canonical ensemble under decay of correlations

abstract: We consider a one-dimensional lattice system of unbounded, real-valued spins. We allow arbitrary strong, attractive, nearest-neighbor interaction. We show that the free energy of the canonical ensemble converges uniformly in $C^2$ to the free energy of the grand canonical ensembles. The error estimates are quantitative. A direct consequence is that the free energy of the canonical ensemble is uniformly strictly convex for large systems. Another consequence is a quantitative local Cram\'er theorem which yields the strict convexity of the coarse-grained Hamiltonian. With small adaptations, the argument could be generalized to systems with finite-range interactions on a graph, as long as the degree of the graph is uniformly bounded and the associated grand canonical ensemble has uniform decay of correlations.


title: A Note on Linear Higher Chow Groups

abstract: We give a counterexample to the proof in the literature [K-Theory 25 (2002), 215-231] of the existence of linear representatives of higher Chow groups of number fields.


title: A class of multi-resolution approximations for large spatial datasets

abstract: Gaussian processes are popular and flexible models for spatial, temporal, and functional data, but they are computationally infeasible for large datasets. We discuss Gaussian-process approximations that use basis functions at multiple resolutions to achieve fast inference and that can (approximately) represent any spatial covariance structure. We consider two special cases of this multi-resolution-approximation framework, a taper version and a domain-partitioning (block) version. We describe theoretical properties and inference procedures, and study the computational complexity of the methods. Numerical comparisons and an application to satellite data are also provided.


title: Making Physics Courses Accessible for Blind Students: strategies for course administration, class meetings and course materials

abstract: The Americans with Disabilities Act (ADA) mandates that U.S. institutions of higher education provide "reasonable accommodations" to students with disabilities to ensure equal educational opportunities. However, despite the key role of physics as a gateway to Science, Technology, Engineering and Mathematics (STEM) studies, only limited resources exist for teaching physics to students who are blind or visually impaired. Here we share lessons from our experience creating an accessible physics curriculum for a blind physics major. The authors include the student himself, a blind physics B.S. who graduated from a different institution, a PhD chemist and consultant on STEM accessibility who is himself blind, and several sighted educators and course assistants who worked regularly with the students. Throughout this effort, we learned that many of the principles of universal design described herein enhanced learning for all of our students.


title: Spectral Density Estimation for Random Fields via Periodic Embeddings

abstract: We introduce methods for estimating the spectral density of a random field on a $d$-dimensional lattice from incomplete gridded data. Data are iteratively imputed onto an expanded lattice according to a model with a periodic covariance function. The imputations are convenient computationally, in that circulant embedding and preconditioned conjugate gradient methods can produce imputations in $O(n\log n)$ time and $O(n)$ memory. However, these so-called periodic imputations are motivated mainly by their ability to produce accurace spectral density estimates. In addition, we introduce a parametric filtering method that is designed to reduce periodogram smoothing bias. The paper contains theoretical results studying properties of the imputed data periodogram and numerical and simulation studies comparing the performance of the proposed methods to existing approaches in a number of scenarios. We present an application to a gridded satellite surface temperature dataset with missing values.


title: Intense automorphisms of finite groups

abstract: Let $G$ be a group. An automorphism of $G$ is called intense if it sends each subgroup of $G$ to a conjugate; the collection of such automorphisms is denoted by $\mathrm{Int}(G)$. In the special case in which $p$ is a prime number and $G$ is a finite $p$-group, one can show that $\mathrm{Int}(G)$ is the semidirect product of a normal $p$-Sylow and a cyclic subgroup of order dividing $p-1$. In this thesis we classify the finite $p$-groups whose groups of intense automorphisms are not themselves $p$-groups. It emerges from our investigation that the structure of such groups is almost completely determined by their nilpotency class: for $p>3$, they share a quotient, growing with their class, with a uniquely determined infinite $2$-generated pro-$p$ group.


title: Energy loss and inelastic diffraction of fast atoms at grazing incidence

abstract: The diffraction of fast atoms at grazing incidence on crystal surfaces (GIFAD) was first interpreted only in terms of elastic diffraction from a perfectly periodic rigid surface with atoms fixed at equilibrium position. Recently, a new approach have been proposed, referred here as the quantum binary collision model (QBCM). The QBCM takes into account both the elastic and inelastic momentum transfer via the Lamb-Dicke probability. It suggests that the shape of the inelastic diffraction profiles are log-normal distributions with a variance proportional to the nuclear energy loss deposited on the surface. For keV Neon atoms impinging the LiF surface, the predictions of the QBCM in its analytic version are compared with numerical trajectory simulations. Some of the assumptions such as the planar continuous form, the possibility to neglect the role of lithium atoms and the influence of temperature are investigated. A specific energy loss dependence $\Delta E\propto\theta^7$ is identified in the quasi-elastic regime merging progressively to the classical onset $\Delta E\propto\theta^3$. The ratio of these two predictions highlight the role of quantum effects in the energy loss.


title: Topological origins of bound states in the continuum for systems with conical intersections

abstract: Bound states in the continuum (BSCs) were reported in a linear vibronic coupling model with a conical intersection (CI) [Cederbaum et al. Phys. Rev. Lett. 90, 013001 (2003)]. It was also found that these states are destroyed within the Born-Oppenheimer approximation (BOA). We investigate whether a nontrivial topological or geometric phase (GP) associated with the CI is responsible for BSCs. To address this question we explore modifications of the original two-dimensional two-state linear vibronic coupling model supporting BSCs. These modifications either add GP effects after the BOA or remove the GP within a two-state problem. Using the stabilization graph technique we shown that the GP is crucial for emergence of BSCs.


title: $t$-cores for $(\Delta+t)$-edge-colouring

abstract: We extend the edge-coloring notion of core (subgraph induced by the vertices of maximum degree) to $t$-core (subgraph induced by the vertices $v$ with $d(v)+\mu(v)> \Delta+t$), and find a sufficient condition for $(\Delta+t)$-edge-coloring. In particular, we show that for any $t\geq 0$, if the $t$-core of $G$ has multiplicity at most $t+1$, with its edges of multiplicity $t+1$ inducing a multiforest, then $\chi'(G) \leq \Delta+t$. This extends previous work of Ore, Fournier, and Berge and Fournier. A stronger version of our result (which replaces the multiforest condition with a vertex-ordering condition) generalizes a theorem of Hoffman and Rodger about cores of $\Delta$-edge-colourable simple graphs. In fact, our bounds hold not only for chromatic index, but for the \emph{fan number} of a graph, a parameter introduced by Scheide and Stiebitz as an upper bound on chromatic index. We are able to give an exact characterization of the graphs $H$ such that $\mathrm{Fan}(G) \leq \Delta(G)+t$ whenever $G$ has $H$ as its $t$-core.


title: Torsion in the homology of finite covers of 3-manifolds

abstract: Let $N$ be a prime 3-manifold that is not a closed graph manifold. Building on a result of Hongbin Sun and using a result of Asaf Hadari we show that for every $k\in\Bbb{N}$ there exists a finite cover $\tilde{N}$ of $N$ such that $|\operatorname{Tor} H_1(\tilde{N};\Bbb{Z})|>k$.


title: Higgs decay mediated by top-quark with flavor-changing neutral scalar interactions

abstract: We explore the flavor-changing parameters mediated by a Higgs boson within the THDM-III context. In particular, the $h\to t^* c$ processes, and check the high suppression for the FC in the THDM-III context for the low $t_\beta$ parameters. Our exploration in the $\chi_{ij}^{u}-\chi_{ij}^{d}-$parameter space shows the allowed regions for different $t_\beta^{}$ values. We explored different modes for Higgs decays, considered the experimental constraints to get scattering plots for the FC parameters and some relevant decay modes. We expect future results to figure out the FC and its implications in the scalar sector.


title: Atomic adsorption on pristine graphene along the Periodic Table of Elements - From PBE to non-local functionals

abstract: The understanding of atomic adsorption on graphene is of high importance for many advanced technologies. Here we present a complete database of the atomic adsorption energies for the elements of the Periodic Table up to the atomic number 86 (excluding lanthanides) on pristine graphene. The energies have been calculated using the projector augmented wave (PAW) method with PBE, long-range dispersion interaction corrected PBE (PBE+D2, PBE+D3) as well as non-local vdW-DF2 approach. The inclusion of dispersion interactions leads to an exothermic adsorption for all the investigated elements. Dispersion interactions are found to be of particular importance for the adsorption of low atomic weight earth alkaline metals, coinage and s-metals (11th and 12th groups), high atomic weight p-elements and noble gases. We discuss the observed adsorption trends along the groups and rows of the Periodic Table as well some computational aspects of modelling atomic adsorption on graphene.


title: Multi-Objective Approaches to Markov Decision Processes with Uncertain Transition Parameters

abstract: Markov decision processes (MDPs) are a popular model for performance analysis and optimization of stochastic systems. The parameters of stochastic behavior of MDPs are estimates from empirical observations of a system; their values are not known precisely. Different types of MDPs with uncertain, imprecise or bounded transition rates or probabilities and rewards exist in the literature. Commonly, analysis of models with uncertainties amounts to searching for the most robust policy which means that the goal is to generate a policy with the greatest lower bound on performance (or, symmetrically, the lowest upper bound on costs). However, hedging against an unlikely worst case may lead to losses in other situations. In general, one is interested in policies that behave well in all situations which results in a multi-objective view on decision making. In this paper, we consider policies for the expected discounted reward measure of MDPs with uncertain parameters. In particular, the approach is defined for bounded-parameter MDPs (BMDPs) [8]. In this setting the worst, best and average case performances of a policy are analyzed simultaneously, which yields a multi-scenario multi-objective optimization problem. The paper presents and evaluates approaches to compute the pure Pareto optimal policies in the value vector space.


title: Motivic virtual signed Euler characteristics and applications to Vafa-Witten invariants

abstract: For any scheme $M$ with a perfect obstruction theory, Jiang and Thomas associate a scheme $N$ with symmetric perfect obstruction theory. The scheme $N$ is a cone over $M$ given by the dual of the obstruction sheaf of $M$, and contains $M$ as its zero section. Locally $N$ is the critical locus of a regular function. In this note we prove that $N$ is a $d$-critical scheme in the sense of Joyce. By assuming an orientation on $N$ there exists a global motive for $N$ locally given by the motive of vanishing cycles of the local regular function. We prove a motivic localization formula under the good and circle compact $\C^*$-action for $N$. When taking Euler characteristic the weighted Euler characteristic of $N$ weighted by the Behrend function is the signed Euler characteristic of $M$ by motivic method. As applications we calculate the motivic generating series of the motivic Vafa-Witten invariants for K3 surfaces. This motivic series gives the result of the $\chi_y$-genus for Vafa-Witten invariants of K3 surfaces, which is the same (at instanton branch) as the K-theoretical Vafa-Witten invariants of Thomas.


title: Finding tight Hamilton cycles in random hypergraphs faster

abstract: In an $r$-uniform hypergraph on $n$ vertices a tight Hamilton cycle consists of $n$ edges such that there exists a cyclic ordering of the vertices where the edges correspond to consecutive segments of $r$ vertices. We provide a first deterministic polynomial time algorithm, which finds a.a.s. tight Hamilton cycles in random $r$-uniform hypergraphs with edge probability at least $C \log^3n/n$. Our result partially answers a question of Dudek and Frieze [Random Structures & Algorithms 42 (2013), 374-385] who proved that tight Hamilton cycles exists already for $p=\omega(1/n)$ for $r=3$ and $p=(e + o(1))/n$ for $r\ge 4$ using a second moment argument. Moreover our algorithm is superior to previous results of Allen, B\"ottcher, Kohayakawa and Person [Random Structures & Algorithms 46 (2015), 446-465] and Nenadov and \v{S}kori\'c [arXiv:1601.04034] in various ways: the algorithm of Allen et al. is a randomised polynomial time algorithm working for edge probabilities $p\ge n^{-1+\varepsilon}$, while the algorithm of Nenadov and \v{S}kori\'c is a randomised quasipolynomial time algorithm working for edge probabilities $p\ge C\log^8n/n$.


title: Obliquely Reflected BSDEs

abstract: In this paper, we study existence and uniqueness to multidimensional Reflected Backward Stochastic Differential Equation in an open convex domain, allowing for oblique directions of reflection. In a Markovian framework, combining \emph{a priori} estimates for penalised equations and compactness arguments, we obtain existence results under quite weak assumptions on the driver of the BSDEs and the direction of reflection, which is allowed to depend on both $Y$ and $Z$. In a non Markovian framework, we obtain existence and uniqueness result for direction of reflection depending on time and $Y$. We make use in this case of stability estimates that require some smoothness condition on the domain and the direction of reflection. In a last Section, we illustrate the application of our theoretical results by introducing randomised switching problems.


title: Quadratic Irrationals, Generating Functions, and L\'evy Constants

abstract: We show that the generating function corresponding to the sequence of denominators of the best rational approximants of a quadratic irrational is a rational function with integer coefficients. Consequently we can compute the L\'evy constant of any quadratic irrational explicitly in terms of a finite number of its convergents.


title: An Extended Mean Field Game for Storage in Smart Grids

abstract: We consider a stylized model for a power network with distributed local power generation and storage. This system is modeled as network connection a large number of nodes, where each node is characterized by a local electricity consumption, has a local electricity production (e.g. photovoltaic panels), and manages a local storage device. Depending on its instantaneous consumption and production rates as well as its storage management decision, each node may either buy or sell electricity, impacting the electricity spot price. The objective at each node is to minimize energy and storage costs by optimally controlling the storage device. In a non-cooperative game setting, we are led to the analysis of a non-zero sum stochastic game with $N$ players where the interaction takes place through the spot price mechanism. For an infinite number of agents, our model corresponds to an Extended Mean-Field Game (EMFG). In a linear quadratic setting, we obtain and explicit solution to the EMFG, we show that it provides an approximate Nash-equilibrium for $N$-player game, and we compare this solution to the optimal strategy of a central planner.


title: Gas Eruption Phenomenon Happening from Ga-In Alloy in Electrolyte

abstract: We report a gas eruption phenomenon caused by electrolysis of liquid Ga-In alloy in an electrolyte, especially NaOH solution. A volcanic eruption-like blowout of gas occurred from the orifice on the alloy surface. In addition to gas plume, large gas bubbles were also generated and the total gas yield increased as In ratio was increased. It is found that destructiveness of the passivation layer on the Ga-In alloy is critical to gas generation. The mechanism of gas eruption can be ascribed to a galvanic interaction happens owing to passivation film and alloy with different activity connected as electrode in electrolyte. Further investigation demonstrated that the lattice of the film expands because of the incorporation of indium, which brings about the decrease in band gap and finally enhances more gas generation. These findings regain the basic understanding of room temperature liquid metal inside electrolyte.


title: On Meta-monoids and the Fox-Milnor Condition

abstract: In this paper we introduce an algebraic structure known as meta-monoids which is particularly suited for the study of knot theory. We define a meta-monoid called $\Gamma$-calculus that gives an Alexander invariant of tangles. We believe that $\Gamma$-calculus gives a unifying framework to study the Alexander polynomial. Specifically, we rederive certain important properties of the Alexander polynomial, most notably the Fox-Milnor condition on the Alexander polynomials of ribbon knots. We argue that our proof has some potential for generalization which may help tackle the slice-ribbon conjecture.


title: Variable Partitioning for Distributed Optimization

abstract: This paper is about how to partition decision variables while decomposing a large-scale optimization problem for the best performance of distributed solution methods. Solving a large-scale optimization problem sequen- tially can be computationally challenging. One classic approach is to decompose the problem into smaller sub-problems and solve them in a distributed fashion. However, there is little discussion in the literature on which variables should be grouped together to form the sub-problems, especially when the optimization formulation involves complex constraints. We focus on one of the most popular distributed approaches, dual decomposition and distributed sub-gradient methods. Based on a theoretical guarantee on its convergence rate, we explain that a partition of variables can critically affect the speed of convergence and highlight the importance of the number of dualized constraints. Then, we introduce a novel approach to find a partition that reduces the number of dualized constraints by utilizing a community detection algorithm from physics literature. Roughly speaking, the proposed method groups decision variables that appear together in con- straints and solves the resulting sub-problems with blocks of variables in parallel. Empirical experiments on a real application show that the proposed method significantly accelerates the convergence of the distributed sub-gradient method. The advantage of our approach becomes more significant as the size of the problem increases and each constraint involves more variables.


title: Valence band-anticrossing in GaP$_{1-x}$Bi$_{x}$ dilute bismide alloys: giant bowing of the band gap and spin-orbit splitting energy

abstract: Using spectroscopic ellipsometry measurements on GaP$_{1-x}$Bi$_{x}$/GaP epitaxial layers up to $x = 3.7$% we observe a giant bowing of the direct band gap ($E_{g}^{\Gamma}$) and valence band spin-orbit splitting energy ($\Delta_{\textrm{SO}}$). $E_{g}^{\Gamma}$ ($\Delta_{\textrm{SO}}$) is measured to decrease (increase) by approximately 200 meV (240 meV) with the incorporation of 1% Bi, corresponding to a greater than fourfold increase in $\Delta_{\textrm{SO}}$ in going from GaP to GaP$_{0.99}$Bi$_{0.01}$. The evolution of $E_{g}^{\Gamma}$ and $\Delta_{\textrm{SO}}$ with $x$ is characterised by strong, composition-dependent bowing. We demonstrate that a simple valence band-anticrossing model, parametrised directly from atomistic supercell calculations, quantitatively describes the measured evolution of $E_{g}^{\Gamma}$ and $\Delta_{\textrm{SO}}$ with $x$. In contrast to the well-studied GaAs$_{1-x}$Bi$_{x}$ alloy, in GaP$_{1-x}$Bi$_{x}$ substitutional Bi creates localised impurity states lying energetically within the GaP host matrix band gap. This leads to the emergence of an optically active band of Bi-hybridised states, accounting for the overall large bowing of $E_{g}^{\Gamma}$ and $\Delta_{\textrm{SO}}$ and in particular for the giant bowing observed for $x \lesssim 1$%. Our analysis provides insight into the action of Bi as an isovalent impurity, and constitutes the first detailed experimental and theoretical analysis of the GaP$_{1-x}$Bi$_{x}$ alloy band structure.


title: Permutation Games for the Weakly Aconjunctive $\mu$-Calculus

abstract: We introduce a natural notion of limit-deterministic parity automata and present a method that uses such automata to construct satisfiability games for the weakly aconjunctive fragment of the $\mu$-calculus. To this end we devise a method that determinizes limit-deterministic parity automata of size $n$ with $k$ priorities through limit-deterministic B\"uchi automata to deterministic parity automata of size $\mathcal{O}((nk)!)$ and with $\mathcal{O}(nk)$ priorities. The construction relies on limit-determinism to avoid the full complexity of the Safra/Piterman-construction by using partial permutations of states in place of Safra-Trees. By showing that limit-deterministic parity automata can be used to recognize unsuccessful branches in pre-tableaux for the weakly aconjunctive $\mu$-calculus, we obtain satisfiability games of size $\mathcal{O}((nk)!)$ with $\mathcal{O}(nk)$ priorities for weakly aconjunctive input formulas of size $n$ and alternation-depth $k$. A prototypical implementation that employs a tableau-based global caching algorithm to solve these games on-the-fly shows promising initial results.


title: Multi-Armed Bandits with Metric Movement Costs

abstract: We consider the non-stochastic Multi-Armed Bandit problem in a setting where there is a fixed and known metric on the action space that determines a cost for switching between any pair of actions. The loss of the online learner has two components: the first is the usual loss of the selected actions, and the second is an additional loss due to switching between actions. Our main contribution gives a tight characterization of the expected minimax regret in this setting, in terms of a complexity measure $\mathcal{C}$ of the underlying metric which depends on its covering numbers. In finite metric spaces with $k$ actions, we give an efficient algorithm that achieves regret of the form $\widetilde{O}(\max\{\mathcal{C}^{1/3}T^{2/3},\sqrt{kT}\})$, and show that this is the best possible. Our regret bound generalizes previous known regret bounds for some special cases: (i) the unit-switching cost regret $\widetilde{\Theta}(\max\{k^{1/3}T^{2/3},\sqrt{kT}\})$ where $\mathcal{C}=\Theta(k)$, and (ii) the interval metric with regret $\widetilde{\Theta}(\max\{T^{2/3},\sqrt{kT}\})$ where $\mathcal{C}=\Theta(1)$. For infinite metrics spaces with Lipschitz loss functions, we derive a tight regret bound of $\widetilde{\Theta}(T^{\frac{d+1}{d+2}})$ where $d \ge 1$ is the Minkowski dimension of the space, which is known to be tight even when there are no switching costs.


title: Twisting Functors and Generalized Verma modules

abstract: Let $\mathfrak{g}$ be a reductive Lie algebra. We give a condition that ensures that the character of a generalized Verma module is well-behaved under a twisting functor. We show that a similar result holds for basic classical simple Lie superalgebras, and generalize a result from \cite{CM} about twisting Verma modules.


title: A robust error estimator and a residual-free error indicator for reduced basis methods

abstract: The Reduced Basis Method (RBM) is a rigorous model reduction approach for solving parametrized partial differential equations. It identifies a low-dimensional subspace for approximation of the parametric solution manifold that is embedded in high-dimensional space. A reduced order model is subsequently constructed in this subspace. RBM relies on residual-based error indicators or {\em a posteriori} error bounds to guide construction of the reduced solution subspace, to serve as a stopping criteria, and to certify the resulting surrogate solutions. Unfortunately, it is well-known that the standard algorithm for residual norm computation suffers from premature stagnation at the level of the square root of machine precision. In this paper, we develop two alternatives to the standard offline phase of reduced basis algorithms. First, we design a robust strategy for computation of residual error indicators that allows RBM algorithms to enrich the solution subspace with accuracy beyond root machine precision. Secondly, we propose a new error indicator based on the Lebesgue function in interpolation theory. This error indicator does not require computation of residual norms, and instead only requires the ability to compute the RBM solution. This residual-free indicator is rigorous in that it bounds the error committed by the RBM approximation, but up to an uncomputable multiplicative constant. Because of this, the residual-free indicator is effective in choosing snapshots during the offline RBM phase, but cannot currently be used to certify error that the approximation commits. However, it circumvents the need for \textit{a posteriori} analysis of numerical methods, and therefore can be effective on problems where such a rigorous estimate is hard to derive.


title: Control Problems with Vanishing Lie Bracket Arising from Complete Odd Circulant Evolutionary Games

abstract: We study an optimal control problem arising from a generalization of rock-paper-scissors in which the number of strategies may be selected from any positive odd number greater than 1 and in which the payoff to the winner is controlled by a control variable $\gamma$. Using the replicator dynamics as the equations of motion, we show that a quasi-linearization of the problem admits a special optimal control form in which explicit dynamics for the controller can be identified. We show that all optimal controls must satisfy a specific second order differential equation parameterized by the number of strategies in the game. We show that as the number of strategies increases, a limiting case admits a closed form for the open-loop optimal control. In performing our analysis we show necessary conditions on an optimal control problem that allow this analytic approach to function.


title: A Sequential Approximation Framework for Coded Distributed Optimization

abstract: Building on the previous work of Lee et al. and Ferdinand et al. on coded computation, we propose a sequential approximation framework for solving optimization problems in a distributed manner. In a distributed computation system, latency caused by individual processors ("stragglers") usually causes a significant delay in the overall process. The proposed method is powered by a sequential computation scheme, which is designed specifically for systems with stragglers. This scheme has the desirable property that the user is guaranteed to receive useful (approximate) computation results whenever a processor finishes its subtask, even in the presence of uncertain latency. In this paper, we give a coding theorem for sequentially computing matrix-vector multiplications, and the optimality of this coding scheme is also established. As an application of the results, we demonstrate solving optimization problems using a sequential approximation approach, which accelerates the algorithm in a distributed system with stragglers.


title: Solving Packing and Covering LPs in $\tilde{O}(\frac{1}{\epsilon^2})$ Distributed Iterations with a Single Algorithm and Simpler Analysis

abstract: Packing and covering linear programs belong to the narrow class of linear programs that are efficiently solvable in parallel and distributed models of computation, yet are a powerful modeling tool for a wide range of fundamental problems in theoretical computer science, operations research, and many other areas. Following recent progress in obtaining faster distributed and parallel algorithms for packing and covering linear programs, we present a simple algorithm whose iteration count matches the best known $\tilde{O}(\frac{1}{\epsilon^2})$ for this class of problems. The algorithm is similar to the algorithm of [Allen-Zhu and Orecchia, 2015], it can be interpreted as Nesterov's dual averaging, and it constructs approximate solutions to both primal (packing) and dual (covering) problems. However, the analysis relies on the construction of an approximate optimality gap and a primal-dual view, leading to a more intuitive interpretation. Moreover, our analysis suggests that all existing algorithms for solving packing and covering linear programs in parallel/distributed models of computation are, in fact, unaccelerated, and raises the question of designing accelerated algorithms for this class of problems.


title: The path to high-energy electron-positron colliders: from Wideroe's betatron to Touschek's AdA and to LEP

abstract: We describe the road which led to the construction and exploitation of electron positron colliders, hightlighting how the young physics student Bruno Touschek met the Norwegian engineer Rolf Wideroe in Germany, during WWII, and collaborated in building the 15 MeV betatron, a secret project directed by Wideroe and financed by the Ministry of Aviation of the Reich. This is how Bruno Touschek learnt the science of making particle accelerators and was ready, many years later, to propose and build AdA, the first electron positron collider, in Frascati, Italy, in 1960. We shall then see how AdA was brought from Frascati to Orsay, in France. Taking advantage of the Orsay Linear Accelerator as injector, the Franco-Italian team was able to prove that collisions had taken place, opening the way to the use of particle colliders as a mean to explore high energy physics.


title: Regular Dilation and Nica-covariant Representation on Right LCM Semigroups

abstract: Regular dilation has recently been extended to graph product of $\mathbb{N}$, where having a *-regular dilation is equivalent to having a minimal isometric Nica-covariant dilation. In this paper, we first extend the result to right LCM semigroups, and establish a similar equivalence among *-regular dilation, minimal isometric Nica-covariant dilation, and a Brehmer-type condition. This result can be applied to various semigroups to establish conditions for *-regular dilation.


title: A uniqueness result for 2-soliton solutions of the KdV equation

abstract: Multisoliton solutions of the KdV equation satisfy nonlinear ordinary differential equations which are known as stationary equations for the KdV hierarchy, or sometimes as Lax-Novikov equations. An interesting feature of these equations, known since the 1970's, is that they can be explicitly integrated, by virtue of being finite-dimensional completely integrable Hamiltonian systems. Here we use the integration theory to investigate the question of whether the multisoliton solutions are the only nonsingular solutions of these ordinary differential equations which vanish at infinity. In particular we prove that this is indeed the case for $2$-soliton solutions of the fourth-order stationary equation.


title: Tidal Stresses and Energy Gaps in Microstate Geometries

abstract: We compute energy gaps and study infalling massive geodesic probes in the new families of scaling, microstate geometries that have been constructed recently and for which the holographic duals are known. We find that in the deepest geometries, which have the lowest energy gaps, the geodesic deviation shows that the stress reaches the Planck scale long before the probe reaches the cap of the geometry. Such probes must therefore undergo a stringy transition as they fall into microstate geometry. We discuss the scales associated with this transition and comment on the implications for scrambling in microstate geometries.


title: QCD phase diagram from chiral symmetry restoration: analytic approach at high and low temperature using the Linear Sigma Model with Quarks

abstract: We use the linear sigma model with quarks to study the QCD phase diagram from the point of view of chiral symmetry restoration. We compute the leading order effective potential for high and low temperatures and finite quark chemical potential, up to the contribution of the ring diagrams to account for the plasma screening effects. We fix the values of the model couplings using physical values for the input parameters such as the vacuum pion and sigma masses, the critical temperature at vanishing quark chemical potential and the conjectured end point value of the baryon chemical potential of the transition line at vanishing temperature. We also make the analysis for the same input parameters but with vanishing pion mass. We find that the critical end point (CEP) is located at low temperatures and high quark chemical potentials $(315<\mu^{\text{CEP}}<349\ MeV,18< T^{\text{CEP}}<45\ MeV)$.


title: The Shape of an Image: A Study of Mapper on Images

abstract: We study the topological construction called Mapper in the context of simply connected domains, in particular on images. The Mapper construction can be considered as a generalization for contour, split, and joint trees on simply connected domains. A contour tree on an image domain assumes the height function to be a piecewise linear Morse function. This is a rather restrictive class of functions and does not allow us to explore the topology for most real world images. The Mapper construction avoids this limitation by assuming only continuity on the height function allowing this construction to robustly deal with a significant larger set of images. We provide a customized construction for Mapper on images, give a fast algorithm to compute it, and show how to simplify the Mapper structure in this case. Finally, we provide a simple procedure that guarantees the equivalence of Mapper to contour, join, and split trees on a simply connected domain.


title: Asymptotic Distribution and Simultaneous Confidence Bands for Ratios of Quantile Functions

abstract: Ratio of medians or other suitable quantiles of two distributions is widely used in medical research to compare treatment and control groups or in economics to compare various economic variables when repeated cross-sectional data are available. Inspired by the so-called growth incidence curves introduced in poverty research, we argue that the ratio of quantile functions is a more appropriate and informative tool to compare two distributions. We present an estimator for the ratio of quantile functions and develop corresponding simultaneous confidence bands, which allow to assess significance of certain features of the quantile functions ratio. Derived simultaneous confidence bands rely on the asymptotic distribution of the quantile functions ratio and do not require re-sampling techniques. The performance of the simultaneous confidence bands is demonstrated in simulations. Analysis of the expenditure data from Uganda in years 1999, 2002 and 2005 illustrates the relevance of our approach.


title: Approximate Span Liftings

abstract: We develop new abstractions for reasoning about relaxations of differential privacy: R\'enyi differential privacy, zero-concentrated differential privacy, and truncated concentrated differential privacy, which express different bounds on statistical divergences between two output probability distributions. In order to reason about such properties compositionally, we introduce approximate span-lifting, a novel construction extending the approximate relational lifting approaches previously developed for standard differential privacy to a more general class of divergences, and also to continuous distributions. As an application, we develop a program logic based on approximate span-liftings capable of proving relaxations of differential privacy and other statistical divergence properties.


title: Proton mass decomposition

abstract: We report the results on the proton mass decomposition and also on the related quark and glue momentum fractions. The results are based on overlap valence fermions on four ensembles of $N_f = 2+1$ DWF configurations with three lattice spacings and volumes, and several pion masses including the physical pion mass. With 1-loop perturbative calculation and proper normalization of the glue operator, we find that the $u, d,$ and $s$ quark masses contribute 9(2)\% to the proton mass. The quark energy and glue field energy contribute 31(5)\% and 37(5)\% respectively in the $\overline{MS}$ scheme at $\mu = 2$ GeV. The trace anomaly gives the remaining 23(1)\% contribution. The $u,d,s$ and glue momentum fractions in the $\overline{MS}$ scheme are consistent with the global analysis at $\mu = 2$ GeV.


title: An Energy-Efficient Mixed-Signal Neuron for Inherently Error-Resilient Neuromorphic Systems

abstract: This work presents the design and analysis of a mixed-signal neuron (MS-N) for convolutional neural networks (CNN) and compares its performance with a digital neuron (Dig-N) in terms of operating frequency, power and noise. The circuit-level implementation of the MS-N in 65 nm CMOS technology exhibits 2-3 orders of magnitude better energy-efficiency over Dig-N for neuromorphic computing applications - especially at low frequencies due to the high leakage currents from many transistors in Dig-N. The inherent error-resiliency of CNN is exploited to handle the thermal and flicker noise of MS-N. A system-level analysis using a cohesive circuit-algorithmic framework on MNIST and CIFAR-10 datasets demonstrate an increase of 3% in worst-case classification error for MNIST when the integrated noise power in the bandwidth is ~ 1 {\mu}V2.


title: Central limit theorem for the free energy of the random field Ising model

abstract: A central limit theorem is proved for the free energy of the random field Ising model with all plus or all minus boundary condition, at any temperature (including zero temperature) and any dimension. This solves a problem posed by Wehr and Aizenman in 1990. The proof uses a variant of Stein's method.


title: Hyperfine interactions of $\text{Er}^{3+}$ ions in $\text{Y}_2\text{SiO}_{5}$: electron paramagnetic resonance in a tunable microwave cavity

abstract: The hyperfine structure of the ground state of erbium doped yttrium orthosilicate is analyzed with the use of electron paramagnetic resonance experiments in a tunable microwave resonator. This work was prompted by the disagreement between the measurements made in zero magnetic field [Phys. Rev. B, 94, 075117, (2016)] and a previously published spin Hamiltonian, which is derived from standard EPR measurements at 9.5 GHz [Phys. Rev. B, 74, 214409, (2006)]. The ability to vary magnetic field strength, resonator frequency, and the orientation of our sample enabled us to monitor how the frequencies of hyperfine transitions change as a function of a vector magnetic field. Compared to [Phys. Rev. B, 74, 214409, (2006)], we arrived at a different set of spin Hamiltonian parameters, which are also broadly consistent with their data. We discuss the reliability of our new spin Hamiltonian parameters to make predictions outside the magnetic field and frequency regimes of our data. We also discuss why it proved to be difficult to determine spin Hamiltonian parameters for this material, and present data collection strategies that improve the model reliability.


title: Shadow-free multimers as extreme-performance meta-atoms

abstract: We generalize the concept of parity-time symmetric structures with the goal to create meta-atoms exhibiting extraordinary abilities to overcome the presumed limitations in the scattering of overall lossless particles, such as non-zero forward scattering and the equality of scattering and extinction powers for all lossless particles. Although the forward scattering amplitude and the extinction cross section of our proposed meta-atoms vanish, they scatter incident energy into other directions, with controllable directionality. These meta-atoms possess extreme electromagnetic properties not achievable for passive scatterers. As an example, we study meta-atoms consisting of two or three small dipole scatters. We consider possible microwave realizations in the form of short dipole antennas loaded by lumped elements. The proposed meta-atom empowers extraordinary response of a shadow-free scatterer and theoretically enables most unusual material properties when used as a building block of an artificial medium.


title: Learning Hidden Quantum Markov Models

abstract: Hidden Quantum Markov Models (HQMMs) can be thought of as quantum probabilistic graphical models that can model sequential data. We extend previous work on HQMMs with three contributions: (1) we show how classical hidden Markov models (HMMs) can be simulated on a quantum circuit, (2) we reformulate HQMMs by relaxing the constraints for modeling HMMs on quantum circuits, and (3) we present a learning algorithm to estimate the parameters of an HQMM from data. While our algorithm requires further optimization to handle larger datasets, we are able to evaluate our algorithm using several synthetic datasets. We show that on HQMM generated data, our algorithm learns HQMMs with the same number of hidden states and predictive accuracy as the true HQMMs, while HMMs learned with the Baum-Welch algorithm require more states to match the predictive accuracy.


title: Hard X-ray Quiescent Emission in Magnetars via Resonant Compton Upscattering

abstract: Non-thermal quiescent X-ray emission extending between 10 keV and around 150 keV has been seen in about 10 magnetars by RXTE, INTEGRAL, Suzaku, NuSTAR and Fermi-GBM. For inner magnetospheric models of such hard X-ray signals, inverse Compton scattering is anticipated to be the most efficient process for generating the continuum radiation, because the scattering cross section is resonant at the cyclotron frequency. We present hard X-ray upscattering spectra for uncooled monoenergetic relativistic electrons injected in inner regions of pulsar magnetospheres. These model spectra are integrated over bundles of closed field lines and obtained for different observing perspectives. The spectral turnover energies are critically dependent on the observer viewing angles and electron Lorentz factor. We find that electrons with energies less than around 15 MeV will emit most of their radiation below 250 keV, consistent with the turnovers inferred in magnetar hard X-ray tails. Electrons of higher energy still emit most of the radiation below around 1 MeV, except for quasi-equatorial emission locales for select pulse phases. Our spectral computations use a new state-of-the-art, spin-dependent formalism for the QED Compton scattering cross section in strong magnetic fields.


title: On groups with $S^2$ Bowditch boundary

abstract: We prove that a relatively hyperbolic pair $(G,P)$ has Bowditch boundary a 2-sphere if and only if it is a 3-dimensional Poincare duality pair. We prove this by studying the relationship between the Bowditch and Dahmani boundaries of relatively hyperbolic groups.


title: On generalized quadrangles with a point regular group of automorphisms

abstract: A generalized quadrangle is a point-line incidence geometry such that any two points lie on at most one line and, given a line $\ell$ and a point $P$ not incident with $\ell$, there is a unique point of $\ell$ collinear with $P$. We study the structure of groups acting regularly on the point set of a generalized quadrangle. In particular, we provide a characterization of the generalized quadrangles with a group of automorphisms acting regularly on both the point set and the line set and show that such a thick generalized quadrangle does not admit a polarity. Moreover, we prove that a group $G$ acting regularly on the point set of a generalized quadrangle of order $(u^2, u^3)$ or $(s,s)$, where $s$ is odd and $s+1$ is coprime to $3$, cannot have any nonabelian minimal normal subgroups.


title: Taming heavy-tailed features by shrinkage

abstract: In this work, we focus on a variant of the generalized linear model (GLM) called corrupted GLM (CGLM) with heavy-tailed features and responses. To robustify the statistical inference on this model, we propose to apply $\ell_4$-norm shrinkage to the feature vectors in the low-dimensional regime and apply elementwise shrinkage to them in the high-dimensional regime. Under bounded fourth moment assumptions, we show that the maximum likelihood estimator (MLE) based on the shrunk data enjoys nearly the minimax optimal rate with an exponential deviation bound. Our simulations demonstrate that the proposed feature shrinkage significantly enhances the statistical performance in linear regression and logistic regression on heavy-tailed data. Finally, we apply our shrinkage principle to guard against mislabeling and image noise in the human-written digit recognition problem. We add an $\ell_4$-norm shrinkage layer to the original neural net and reduce the testing misclassification rate by more than $30\%$ relatively in the presence of mislabeling and image noise.


title: Very compact millimeter sizes for composite star-forming/AGN submillimeter galaxies

abstract: We report the study of far-IR sizes of submillimeter galaxies (SMGs) in relation to their dust-obscured star formation rate (SFR) and active galactic nuclei (AGN) presence, determined using mid-IR photometry. We determined the millimeter-wave ($\lambda_{\rm obs}=1100 \mu$m) sizes of 69 ALMA-identified SMGs, selected with $\geq10$$\sigma$ confidence on ALMA images ($F_{\rm 1100 \mu m}=1.7$--7.4 mJy). We found that all the SMGs are located above an avoidance region in the millimeter size-flux plane, as expected by the Eddington limit for star formation. In order to understand what drives the different millimeter-wave sizes in SMGs, we investigated the relation between millimeter-wave size and AGN fraction for 25 of our SMGs at $z=1$--3. We found that the SMGs for which the mid-IR emission is dominated by star formation or AGN have extended millimeter-sizes, with respective median $R_{\rm c,e} = 1.6^{+0.34}_{-0.21}$ and 1.5$^{+0.93}_{-0.24}$ kpc. Instead, the SMGs for which the mid-IR emission corresponds to star-forming/AGN composites have more compact millimeter-wave sizes, with median $R_{\rm c,e}=1.0^{+0.20}_{-0.20}$ kpc. The relation between millimeter-wave size and AGN fraction suggests that this size may be related to the evolutionary stage of the SMG. The very compact sizes for composite star-forming/AGN systems could be explained by supermassive black holes growing rapidly during the SMG coalescing, star-formation phase.


title: Current transport properties and phase diagram of a Kitaev chain with long-range pairing

abstract: We describe a method to probe the quantum phase transition between the short-range topological phase and the long-range topological phase in the superconducting Kitaev chain with long-range pairing, both exhibiting subgap modes localized at the edges. The method relies on the effects of the finite mass of the subgap edge modes in the long-range regime (which survives in the thermodynamic limit) on the single-particle scattering coefficients through the chain connected to two normal leads. Specifically, we show that, when the leads are biased at a voltage V with respect to the superconducting chain, the Fano factor is either zero (in the short-range correlated phase) or 2e (in the long-range correlated phase). As a result, we find that the Fano factor works as a directly measurable quantity to probe the quantum phase transition between the two phases. In addition, we note a remarkable "critical fractionalization effect" in the Fano factor, which is exactly equal to e along the quantum critical line. Finally, we note that a dual implementation of our proposed device makes it suitable as a generator of large-distance entangled two-particle states.


title: Near-infrared laser thermal conjunctivoplasty

abstract: Conjunctivochalasis is a common cause of tear dysfunction due to the conjunctiva becoming loose and wrinkly with age. The current solutions to this disease include either surgical excision in the operating room, or thermoreduction of the loose tissue with hot wire in the clinic. We developed a near-infrared (NIR) laser thermal conjunctivoplasty (LTC) system, which gently shrinks the redundant tissue. The NIR light is mainly absorbed by water, so the heating is even and there is no bleeding. The system utilizes a 1460-nm programmable laser diode system as a light source. A miniaturized handheld probe delivers the laser light and focuses the laser into a 10x1 mm2 line. A foot pedal is used to deliver a preset number of calibrated laser pulses. A fold of loose conjunctiva is grasped by a pair of forceps. The infrared laser light is delivered through an optical fiber and a laser line is focused exactly on the conjunctival fold by a cylindrical lens. Ex vivo experiments using porcine eye were performed with the optimal laser parameters. It was found that up to 50% of conjunctiva shrinkage could be achieved.


title: Polarization-multiplexed, dual-beam swept source optical coherence tomography angiography

abstract: A polarization-multiplexed, dual-beam setup is proposed to expand the field of view for a swept source optical coherence tomography angiography (OCTA) system. This method used a Wollaston prism to split sample path light into two orthogonal-polarized beams. This allowed two beams to shine on the cornea at an angle separation of ~ 14 degrees, which led to a separation of ~ 4.2 mm on the retina. A 3-mm glass plate was inserted into one of the beam paths to set a constant path length difference between the two polarized beams so the interferogram from the two beams are coded at different frequency bands. The resulting OCTA images from the two beams were coded with a depth separation of ~ 2 mm. 5x5 mm2 angiograms from the two beams were obtained simultaneously in 4 seconds. The two angiograms then were montaged to get a wider field of view (FOV) of ~ 5x9.2 mm2.


title: Optimization of population annealing Monte Carlo for large-scale spin-glass simulations

abstract: Population annealing Monte Carlo is an efficient sequential algorithm for simulating k-local Boolean Hamiltonians. Because of its structure, the algorithm is inherently parallel and therefore well suited for large-scale simulations of computationally hard problems. Here we present various ways of optimizing population annealing Monte Carlo using 2-local spin-glass Hamiltonians as a case study. We demonstrate how the algorithm can be optimized from an implementation, algorithmic accelerator, as well as scalable parallelization points of view. This makes population annealing Monte Carlo perfectly suited to study other frustrated problems such as pyrochlore lattices, constraint-satisfaction problems, as well as higher-order Hamiltonians commonly found in, e.g., topological color codes.


title: Trace norm regularization and faster inference for embedded speech recognition RNNs

abstract: We propose and evaluate new techniques for compressing and speeding up dense matrix multiplications as found in the fully connected and recurrent layers of neural networks for embedded large vocabulary continuous speech recognition (LVCSR). For compression, we introduce and study a trace norm regularization technique for training low rank factored versions of matrix multiplications. Compared to standard low rank training, we show that our method leads to good accuracy versus number of parameter trade-offs and can be used to speed up training of large models. For speedup, we enable faster inference on ARM processors through new open sourced kernels optimized for small batch sizes, resulting in 3x to 7x speed ups over the widely used gemmlowp library. Beyond LVCSR, we expect our techniques and kernels to be more generally applicable to embedded neural networks with large fully connected or recurrent layers.


title: User-centric Composable Services: A New Generation of Personal Data Analytics

abstract: Machine Learning (ML) techniques, such as Neural Network, are widely used in today's applications. However, there is still a big gap between the current ML systems and users' requirements. ML systems focus on improving the performance of models in training, while individual users cares more about response time and expressiveness of the tool. Many existing research and product begin to move computation towards edge devices. Based on the numerical computing system Owl, we propose to build the Zoo system to support construction, compose, and deployment of ML models on edge and local devices.


title: Electronic Structure of Ce-Doped and -Undoped Nd$_2$CuO$_4$ Superconducting Thin Films Studied by Hard X-ray Photoemission and Soft X-ray Absorption Spectroscopy

abstract: In order to realize superconductivity in cuprates with the T'-type structure, not only chemical substitution (Ce doping) but also post-growth reduction annealing is necessary. In the case of thin films, however, well-designed reduction annealing alone without Ce doping can induce superconductivity in the T'-type cuprates. In order to unveil the origin of superconductivity in the Ce-undoped T'-type cuprates, we have performed bulk-sensitive hard x-ray photoemission and soft x-ray absorption spectroscopies on superconducting and non-superconducting Nd$_{2-x}$Ce$_x$CuO$_4$ ($x=$ 0, 0.15, and 0.19) thin films. By post-growth annealing, core-level spectra exhibited dramatic changes, which we attributed to the enhancement of core-hole screening in the CuO$_2$ plane and the shift of chemical potential along with changes in the band filling. The result suggests that the superconducting Nd$_2$CuO$_4$ film is doped with electrons and that the electronic states are similar to those of Ce-doped superconductors.


title: What is the Optimal Network Deployment for a Fixed Density of Antennas?

abstract: In this paper, we answer a fundamental question: when the total number of antennas per square kilometer is fixed, what is the optimal network deployment? A denser network with a less number of antennas per base station (BS) or the opposite case. To evaluate network performance, we consider a practical network scenario with a fixed antennas density and multiuser multiple-input-multiple-output (MU-MIMO) operations for single-antenna users. The number of antennas in each BS is calculated by dividing the antenna density by the BS density. With the consideration of several practical network models, i.e., pilot contamination, a limited user equipment (UE) density and probabilistic line-of-sight (LoS)/non-line-of-sight (NLoS) path loss model, we evaluate the area spectral efficiency (ASE) performance. From our simulation results, we conclude that there exists an optimal BS density for a certain UE density to maximize the ASE performance when the antenna density is fixed. The intuition is that (i) by densifying the network with more BSs, we can achieve a receive power gain due to the smaller distance between the typical UE and its serving BS; (ii) by installing more antennas in each BS, we can achieve a beamforming gain for UEs using MU-MIMO, although such beamforming gain is degraded by pilot contamination; (iii) thus, a trade-off exists between the receive power gain and the beamforming gain, if we fix the antenna density in the network.


title: Using Positional Heel-marker Data to More Accurately Calculate Stride Length for Treadmill Walking: A Step Length Approach

abstract: Treadmill walking is a convenient tool for studying the human gait; however, a common gait parameter, stride length, can be difficult to calculate directly because relevant reference points continually move backwards. Although there is no direct calculation of stride length itself, we can use positional heel-marker data to directly determine a similar parameter, step length, and we can sum two step lengths to result in one stride length. This proposed method of calculation is simple but seems to be unexplored in other literature, so this paper displays the details of the calculation. Our experimental results differed from the expected values by 2.2% and had a very low standard deviation, suggesting that this method is viable for practical use. The ability to calculate stride length for treadmill walking using heel-marker data may allow for quick and accurate gait calculations that further contribute to the versatility of heel data as a tool for gait analysis.


title: Massive MIMO Performance Comparison of Beamforming and Multiplexing in the Terahertz Band

abstract: In this paper, we compare the performance of two main MIMO techniques, beamforming and multiplexing, in the Terahertz (THz) band. The main problem with the THz band is its huge propagation loss, which is caused by the tremendous signal attenuation due to molecule absorption of the electro-magnetic wave. To overcome the path loss issue, massive MIMO has been suggested to be employed in the network and is expected to provide Tbps for a distance within a few meters. In this context, beamforming is studied recently as the main technique to take advantage of MIMO in THz and overcome the very high path loss with the assumption that the THz communication channel is Line-of-Sight (LoS) and there are not significant multipath rays. On the other hand, recent studies also showed that the well-known absorbed energy by molecules can be re-radiated immediately in the same frequency. Such re-radiated signal is correlated with the main signal and can provide rich scattering paths for the communication channel. This means that a significant MIMO multiplexing gain can be achieved even in a LoS scenario for the THz band. Our simulation results reveal a surprising observation that the MIMO multiplexing could be a better choice than the MIMO beamforming under certain conditions in THz communications.


title: A New Look at MIMO Capacity in the Millimeter Wave

abstract: In this paper, we present a new theoretical discovery that the multiple-input and multiple-output (MIMO) capacity can be influenced by atmosphere molecules. In more detail, some common atmosphere molecules, such as Oxygen and water, can absorb and re-radiate energy in their natural resonance frequencies, such as 60GHz, 120GHz, and 180GHz, which belong to the millimeter wave (mmWave) spectrum. Such phenomenon can provide equivalent non-line-of-sight (NLoS) paths in an environment that lacks scatterers, and thus greatly improve the spatial multiplexing and diversity of a MIMO system. This kind of performance improvement is particularly useful for most mmWave communications that heavily rely on line-of-sight (LoS) transmissions. To sum up, our study concludes that since the molecular re-radiation happens at certain mmWave frequency bands, the MIMO capacity becomes highly frequency selective and enjoys a considerable boosting at those mmWave frequency bands. The impact of our new discovery is significant, which fundamentally changes our understanding on the relationship between the MIMO capacity and the frequency spectrum. In particular, our results predict that several mmWave bands can serve as valuable spectrum windows for high-efficiency MIMO communications, which in turn may shift the paradigm of research, standardization, and implementation in the field of mmWave communications.


title: Initial Results from the ALFABURST Survey

abstract: Here, we present initial results from the ALFABURST radio transient survey, which is currently running in a commensal mode with the ALFA receiver at the Arecibo telescope. We observed for a total of 1400 hours and have detected single pulses from known pulsars but did not detect any FRBs. The non-detection of FRBs is consistent with the current FRB sky rates.


title: Energy Management for Energy Harvesting Wireless Sensors with Adaptive Retransmission

abstract: This paper analyzes the communication between two energy harvesting wireless sensor nodes. The nodes use automatic repeat request and forward error correction mechanism for the error control. The random nature of available energy and arrivals of harvested energy may induce interruption to the signal sampling and decoding operations. We propose a selective sampling scheme where the length of the transmitted packet to be sampled depends on the available energy at the receiver. The receiver performs the decoding when complete samples of the packet are available. The selective sampling information bits are piggybacked on the automatic repeat request messages for the transmitter use. This way, the receiver node manages more efficiently its energy use. Besides, we present the partially observable Markov decision process formulation, which minimizes the long-term average pairwise error probability and optimizes the transmit power. Optimal and suboptimal power assignment strategies are introduced for retransmissions, which are adapted to the selective sampling and channel state information. With finite battery size and fixed power assignment policy, an analytical expression for the average packet drop probability is derived. Numerical simulations show the performance gain of the proposed scheme with power assignment strategy over the conventional scheme.


title: The Geodesic $2$-center Problem in a Simple Polygon

abstract: The geodesic $k$-center problem in a simple polygon with $n$ vertices consists in the following. Find a set $S$ of $k$ points in the polygon that minimizes the maximum geodesic distance from any point of the polygon to its closest point in $S$. In this paper, we focus on the case where $k=2$ and present an exact algorithm that returns a geodesic $2$-center in $O(n^2\log^2 n)$ time.


title: Hartree-Fock symmetry breaking around conical intersections

abstract: We study the behavior of Hartree-Fock (HF) solutions in the vicinity of conical intersections. These are here understood as regions of a molecular potential energy surface characterized by degenerate or nearly-degenerate eigenfunctions with identical quantum numbers (point group, spin, and electron number). Accidental degeneracies between states with different quantum numbers are known to induce symmetry breaking in HF. The most common closed-shell restricted HF instability is related to singlet-triplet spin degeneracies that lead to collinear unrestricted HF (UHF) solutions. Adding geometric frustration to the mix usually results in noncollinear generalized HF (GHF) solutions, identified by orbitals that are linear combinations of up and down spins. Near conical intersections, we observe the appearance of coplanar GHF solutions that break all symmetries, including complex conjugation and time-reversal, which do not carry good quantum numbers. We discuss several prototypical examples taken from the conical intersection literature. Additionally, we utilize a recently introduced a magnetization diagnostic to characterize these solutions, as well as a solution of a Jahn-Teller active geometry of H$_8^{+2}$.


title: Decay behaviors of the $P_c$ hadronic molecules

abstract: In this proceeding, we present our recent work on decay behaviors of the $P_c$ hadronic molecules, which can help to disentangle the nature of the two $P_c$ pentaquark-like structures. The results turn out that the relative ratio of the decays of $P^+_c(4380)$ to $\bar{D}^* \Lambda_c$ and $J/\psi p$ is very different for $P_c$ being a $\bar D^*\Sigma_c$ or $\bar D\Sigma_c^*$ bound state with $J^P=\frac{3}{2}^-$. And from the total decay width, we find that $P_c(4380)$ being a $\bar D\Sigma_c^*$ molecule state with $J^P=\frac{3}{2}^-$ and $P_c(4450)$ being a $\bar D^*\Sigma_c$ molecule state with $J^P=\frac{5}{2}^+$ is more favorable to the experimental data.


title: The ISLAnds Project III: Variable Stars in Six Andromeda Dwarf Spheroidal Galaxies

abstract: We present a census of variable stars in six M31 dwarf spheroidal satellites observed with the Hubble Space Telescope. We detect 870 RR Lyrae (RRL) stars in the fields of And I (296), II (251), III (111), XV (117), XVI (8), XXVIII (87). We also detect a total of 15 Anomalous Cepheids, three Eclipsing Binaries, and seven field RRL stars compatible with being members of the M31 halo or the Giant Stellar Stream. We derive robust and homogeneous distances to the six galaxies using different methods based on the properties of the RRL stars. Working with the up-to-date set of Period-Wesenheit ($I$, $B$ - $I$) relations published by Marconi et al., we obtain distance moduli of $\mu_0$ = [24.49, 24.16, 24.36, 24.42, 23.70, 24.43] mag (respectively), with systematic uncertainties of 0.08 mag and statistical uncertainties $<$ 0.11 mag. We have considered an enlarged sample of sixteen M31 satellites with published variability studies, and compared their pulsational observables (e.g., periods, amplitudes), with those of fifteen Milky Way satellites for which similar data are available. The properties of the (strictly old) RRL in both satellite systems do not show any significant difference. In particular, we found a strikingly similar correlation between the mean period distribution of the fundamental RRL pulsators (RRab) and the mean metallicities of the galaxies. This indicates that the old RRL progenitors were similar at the early stage in the two environments, suggesting very similar characteristics for the earliest stages of evolution of both satellite systems.


title: Janus: An Uncertain Cache Architecture to Cope with Side Channel Attacks

abstract: Side channel attacks are a major class of attacks to crypto-systems. Attackers collect and analyze timing behavior, I/O data, or power consumption in these systems to undermine their effectiveness in protecting sensitive information. In this work, we propose a new cache architecture, called Janus, to enable crypto-systems to introduce randomization and uncertainty in their runtime timing behavior and power utilization profile. In the proposed cache architecture, each data block is equipped with an on-off flag to enable/disable the data block. The Janus architecture has two special instructions in its instruction set to support the on-off flag. Beside the analytical evaluation of the proposed cache architecture, we deploy it in an ARM-7 processor core to study its feasibility and practicality. Results show a significant variation in the timing behavior across all the benchmarks. The new secure processor architecture has minimal hardware overhead and significant improvement in protecting against power analysis and timing behavior attacks.


title: UGC 4703 interacting pair near to an isolated spiral galaxy NGC 2718: A Milky Way Magellanic Cloud analogue

abstract: We present an analysis of physical and morphological properties of an interacting pair of dwarf galaxies UGC 4703 located at the vicinity of an isolated Milky-Way (MW) type spiral galaxy NGC 2718. Based on the comparison of physical and morphological properties with that of the Large and Small Magellanic Clouds (LMC and SMC), we report that the UGC 4703 pair-NGC 2718 system is probably an LMC-SMC-MW analogue. Located at a sky projected distance of 81 kpc from NGC 2718, we find that UGC 4703 is clearly interacting with its nearby smaller mass companion UGC 4703B forming a bridge of stellar stream between them. Total B-band luminosity of UGC 4703 and its companion is -17.75 and -16.25 mag, respectively. We obtained HI 21-cm line data of UGC 4703 using the GMRT to get a more detailed view of neutral hydrogen (HI) emission. The HI image revealed evidence of interaction between the dwarf galaxy pair but no extended emission, like Magellanic Stream (MS), was detected between NGC 2781 and the UGC 4703 dwarf pair. We also detected star-forming regions along the UGC 4703/ 4703B bridge with stellar mass exceeding 10$^{7}$ M$_{\odot}$. While comparing the optical and HI morphology of the interacting dwarf pairs (UGC 4703-4703B and LMC-SMC) we discuss possible difference in interaction histories of these systems.


title: Generalized Geometric Programming for Rate Allocation in Consensus

abstract: Distributed averaging, or distributed average consensus, is a common method for computing the sample mean of the data dispersed among the nodes of a network in a decentralized manner. By iteratively exchanging messages with neighbors, the nodes of the network can converge to an agreement on the sample mean of their initial states. In real-world scenarios, these messages are subject to bandwidth and power constraints, which motivates the design of a lossy compression strategy. Few prior works consider the rate allocation problem from the perspective of constrained optimization, which provides a principled method for the design of lossy compression schemes, allows for the relaxation of certain assumptions, and offers performance guarantees. We show for Gaussian-distributed initial states with entropy-coded scalar quantization and vector quantization that the coding rates for distributed averaging can be optimized through generalized geometric programming. In the absence of side information from past states, this approach finds a rate allocation over nodes and iterations that minimizes the aggregate coding rate required to achieve a target mean square error within a finite run time. Our rate allocation is compared to some of the prior art through numerical simulations. The results motivate the incorporation of side-information through differential or predictive coding to improve rate-distortion performance.


title: Control Policies Approaching HGI Performance in Heavy Traffic for Resource Sharing Networks

abstract: We consider resource sharing networks of the form introduced in the work of Massouli\'{e} and Roberts(2000) as models for Internet flows. The goal is to study the open problem, formulated in Harrison et al. (2014), of constructing simple form rate allocation policies for broad families of resource sharing networks with associated costs converging to the Hierarchical Greedy Ideal performance in the heavy traffic limit. We consider two types of cost criteria, an infinite horizon discounted cost, and a long time average cost per unit time. We introduce a sequence of rate allocation control policies that are determined in terms of certain thresholds for the scaled queue length processes and prove that, under conditions, both type of costs associated with these policies converge in the heavy traffic limit to the corresponding HGI performance. The conditions needed for these results are satisfied by all the examples considered in Harrison et al. (2014).


title: Construction of Anti-Cyclotomic Euler Systems of Abelian Varieties Associated to $X_1(N)$

abstract: Let $K$ be an imaginary quadratic field, $N$ be a positive integer, $f(z)$ be a newform of level $\Gamma_1(N)$, and $A_f$ be the abelian variety associated to $f$. For each $\tau \in K$ ($\operatorname{Im} \tau >0$), we construct a certain point $P_\tau$ on $A_f$ defined over an extended ring class field of $K$ of level $N$. Our construction generalizes Birch's construction of the Heegner points to the abelian varieties associated to modular forms of level $\Gamma_1(N)$ and nontrivial character. Then, we show that $P_\tau$'s satisfy the distribution and congruence relations of an Euler system, which implies that it should be possible to apply the Euler system techniques to them to show a relation between the non-torsionness of $P_\tau$ and the rank of $A_f(K)$.


title: On the effective cone of higher codimension cycles in $\overline{\mathcal{M}}_{g,n}$

abstract: We exhibit infinitely many extremal effective codimension-$k$ cycles in $\overline{\mathcal{M}}_{g,n}$ in the cases $g\geq 3, n\geq g-1$ and $k=2$, $g\geq 2$, $k\leq n-g,g,$ and $g=1$, $k\leq n-2$. Hence in these cases the effective cone is not rational polyhedral.


title: Kaledin's degeneration theorem and topological Hochschild homology

abstract: We give a short proof of Kaledin's theorem on the degeneration of the noncommutative Hodge-to-de Rham spectral sequence. Our approach is based on topological Hochschild homology and the theory of cyclotomic spectra. As a consequence, we also obtain relative versions of the degeneration theorem, both in characteristic zero and for regular bases in characteristic $p$.


title: Stabilization of Polar Nano Regions in Pb-free ferroelectrics

abstract: Formation of polar nano regions through solid-solution additions are known to enhance significantly the functional properties of ferroelectric materials. Despite considerable progress in characterizing the microscopic behavior of polar nano regions, understanding their real-space atomic structure and dynamics of formation remains a considerable challenge. Here, using the method of dynamic pair distribution function, we provide direct insights into the role of solid-solution additions towards the stabilization of polar nano regions in the Pb-free ferroelectric of Ba(Zr,Ti)O3. It is shown that for an optimum level of substitution of Ti by larger Zr ions, the dynamics of atomic displacements for ferroelectric polarization are slowed sufficiently, which leads to increased local correlation among dipoles below THz frequencies. The dynamic pair distribution function technique demonstrates unique capability to obtain insights into locally correlated atomic dynamics in disordered materials, including new Pb-free ferroelectrics, which is necessary to understand and control their functional properties.


title: Block Coordinate Descent Only Converge to Minimizers

abstract: Given a non-convex twice continuously differentiable cost function with Lipschitz continuous gradient, we prove that all of block coordinate gradient descent, block mirror descent and proximal block coordinate descent converge to a local minimizer, almost surely with random initialization. Furthermore, we show that these results also hold true even for the cost functions with non-isolated critical points.


title: Edge-outer graph embedding and the complexity of the DNA reporter strand problem

abstract: In 2009, Jonoska, Seeman and Wu showed that every graph admits a route for a DNA reporter strand, that is, a closed walk covering every edge either once or twice, in opposite directions if twice, and passing through each vertex in a particular way. This corresponds to showing that every graph has an \emph{edge-outer embedding}, that is, an orientable embedding with some face that is incident with every edge. In the motivating application, the objective is such a closed walk of minimum length. Here we give a short algorithmic proof of the original existence result, and also prove that finding a shortest length solution is NP-hard, even for $3$-connected cubic ($3$-regular) planar graphs. Independent of the motivating application, this problem opens a new direction in the study of graph embeddings, and we suggest new problems emerging from it.


title: Functional analytic approach to Ces\`aro mean

abstract: We study a certain class $\mathcal{P}$ of positive linear functionals $\varphi$ on $L^{\infty}([1,\infty))$ for which $\varphi(f) = \alpha$ if $\lim_{x \to \infty} \frac{1}{x} \int_1^x f(t)dt = \alpha$. It turns out that translations $f(x) \mapsto f(rx)$ on $L^{\infty}([1, \infty))$, where $r \in [1, \infty)$, which are induced by the action of the multiplicative semigroup $[1, \infty)$ on itself, plays an intrinsic role in the study of $\mathcal{P}$. We also deal with an analogue $\mathcal{K}$ of $\mathcal{P}$ of positive linear functionals on $L^{\infty}([0, \infty))$ partaining to the action of the additive semigroup $[0, \infty)$ on itself. In particular, we give some expressions of maximal possible values of $\mathcal{P}$ and $\mathcal{K}$ for a given function respectively.


title: Lattice points in stretched model domains of finite type in $\mathbb{R}^d$

abstract: We study an optimal stretching problem for certain convex domain in $\mathbb{R}^d$ ($d\geq 3$) whose boundary has points of vanishing Gaussian curvature. We prove that the optimal domain which contains the most positive (or least nonnegative) lattice points is asymptotically balanced. This type of problem has its origin in the "eigenvalue optimization among rectangles" problem in spectral geometry. Our proof relies on two-term bounds for lattice counting for general convex domains in $\mathbb{R}^d$ and an explicit estimate of the Fourier transform of the characteristic function associated with the specific domain under consideration.


title: Tuning the Dzyaloshinskii-Moriya Interaction in Pt/Co/MgO heterostructures through MgO thickness

abstract: The interfacial Dzyaloshinskii-Moriya interaction (DMI) in the ferromagnetic/heavy metal ultra-thin film structures , has attracted a lot of attention thanks to its capability to stabilize Neel-type domain walls (DWs) and magnetic skyrmions for the realization of non-volatile memory and logic devices. In this study, we demonstrate that magnetic properties in perpendicularly magnetized Ta/Pt/Co/MgO/Pt heterostructures, such as magnetization and DMI, can be significantly influenced through both the MgO and the Co ultrathin film thickness. By using a field-driven creep regime domain expansion technique, we find that non-monotonic tendencies of DMI field appear when changing the thickness of MgO and the MgO thickness corresponding to the largest DMI field varies as a function of the Co thicknesses. We interpret this efficient control of DMI as subtle changes of both Pt/Co and Co/MgO interfaces, which provide a method to investigate ultra-thin structures design to achieve skyrmion electronics.