title: On Cross-validation for Sparse Reduced Rank Regression

abstract: In high-dimensional data analysis, regularization methods pursuing sparsity and/or low rank have received a lot of attention recently. To provide a proper amount of shrinkage, it is typical to use a grid search and a model comparison criterion to find the optimal regularization parameters. However, we show that fixing the parameters across all folds may result in an inconsistency issue, and it is more appropriate to cross-validate projection-selection patterns to obtain the best coefficient estimate. Our in-sample error studies in jointly sparse and rank-deficient models lead to a new class of information criteria with four scale-free forms to bypass the estimation of the noise level. By use of an identity, we propose a novel scale-free calibration to help cross-validation achieve the minimax optimal error rate non-asymptotically. Experiments support the efficacy of the proposed methods.


title: On the structure of distance sets over prime fields

abstract: Let $\mathbb{F}_q$ be a finite field of order $q$ and $\mathcal{E}$ be a set in $\mathbb{F}_q^d$. The distance set of $\mathcal{E}$, denoted by $\Delta(\mathcal{E})$, is the set of distinct distances determined by the pairs of points in $\mathcal{E}$. Very recently, Iosevich, Koh, and Parshall (2018) proved that if $|\mathcal{E}|\gg q^{d/2}$, then the quotient set of $\Delta(\mathcal{E})$ satisfies \[\left\vert\frac{\Delta(\mathcal{E})}{\Delta(\mathcal{E})}\right\vert=\left\vert \left\lbrace\frac{a}{b}\colon a, b\in \Delta(\mathcal{E}), b\ne 0\right\rbrace\right\vert\gg q.\] In this paper, we break the exponent $d/2$ when $\mathcal{E}$ is a Cartesian product of sets over a prime field. More precisely, let $p$ be a prime and $A\subset \mathbb{F}_p$. If $\mathcal{E}=A^d\subset \mathbb{F}_p^d$ and $|\mathcal{E}|\gg p^{\frac{d}{2}-\varepsilon}$ for some $\varepsilon>0$, then we have \[\left\vert\frac{\Delta(\mathcal{E})}{\Delta(\mathcal{E})}\right\vert, ~\left\vert \Delta(\mathcal{E})\cdot \Delta(\mathcal{E})\right\vert \gg p.\] Such improvements are not possible over arbitrary finite fields. These results give us a better understanding about the structure of distance sets and the Erd\H{o}s-Falconer distance conjecture over finite fields.


title: Weak field limit and gravitational waves in higher-order gravity

abstract: We derive the weak field limit for a gravitational Lagrangian density $L_{g}=(R+a_{0}R^{2}+\sum_{k=1}^{p} a_{k}R\Box^{k}R)\sqrt{-g}$ where higher-order derivative terms in the Ricci scalar $R$ are taken into account. The interest for this kind of effective theories comes out from the consideration of the infrared and ultraviolet behaviors of gravitational field and, in general, from the formulation of quantum field theory in curved spacetimes. Here, we obtain solutions in weak field regime both in vacuum and in the presence of matter and derive gravitational waves considering the contribution of $R\Box^{k}R$ terms. By using a suitable set of coefficients $a_{k}$, it is possible to find up to $(p+2)$ normal modes of oscillation with six polarization states with helicity 0 or 2. Here $p$ is the higher order term in the $\Box$ operator appearing in the gravitational Lagrangian. More specifically: the mode $\omega_{1}$, with $k^{2}=0$, has transverse polarizations $\epsilon_{\mu\nu}^{\left(+\right)}$ and $\epsilon_{\mu\nu}^{\left(\times\right)}$ with helicity 2; the $(p+1)$ modes $\omega_{m}$, with $k^{2}\neq0$, have transverse polarizations $\epsilon_{\mu\nu}^{\left(1\right)}$ and non-transverse ones $\epsilon_{\mu\nu}^{\left(\text{TT}\right)}$, $\epsilon_{\mu\nu}^{\left(\text{TS}\right)}$, $\epsilon_{\mu\nu}^{\left(L\right)}$ with helicity 0.


title: Expander Graphs -- Both Local and Global

abstract: Let $G=(V,E)$ be a finite graph. For $v\in V$ we denote by $G_v$ the subgraph of $G$ that is induced by $v$'s neighbor set. We say that $G$ is $(a,b)$-regular for $a>b>0$ integers, if $G$ is $a$-regular and $G_v$ is $b$-regular for every $v\in V$. Recent advances in PCP theory call for the construction of infinitely many $(a,b)$-regular expander graphs $G$ that are expanders also locally. Namely, all the graphs $\{G_v|v\in V\}$ should be expanders as well. While random regular graphs are expanders with high probability, they almost surely fail to expand locally. Here we construct two families of $(a,b)$-regular graphs that expand both locally and globally. We also analyze the possible local and global spectral gaps of $(a,b)$-regular graphs. In addition, we examine our constructions vis-a-vis properties which are considered characteristic of high-dimensional expanders.


title: Variational Self-attention Model for Sentence Representation

abstract: This paper proposes a variational self-attention model (VSAM) that employs variational inference to derive self-attention. We model the self-attention vector as random variables by imposing a probabilistic distribution. The self-attention mechanism summarizes source information as an attention vector by weighted sum, where the weights are a learned probabilistic distribution. Compared with conventional deterministic counterpart, the stochastic units incorporated by VSAM allow multi-modal attention distributions. Furthermore, by marginalizing over the latent variables, VSAM is more robust against overfitting. Experiments on the stance detection task demonstrate the superiority of our method.


title: Monte-Carlo Sampling applied to Multiple Instance Learning for Histological Image Classification

abstract: We propose a patch sampling strategy based on a sequential Monte-Carlo method for high resolution image classification in the context of Multiple Instance Learning. When compared with grid sampling and uniform sampling techniques, it achieves higher generalization performance. We validate the strategy on two artificial datasets and two histological datasets for breast cancer and sun exposure classification.


title: Learning to Selectively Transfer: Reinforced Transfer Learning for Deep Text Matching

abstract: Deep text matching approaches have been widely studied for many applications including question answering and information retrieval systems. To deal with a domain that has insufficient labeled data, these approaches can be used in a Transfer Learning (TL) setting to leverage labeled data from a resource-rich source domain. To achieve better performance, source domain data selection is essential in this process to prevent the "negative transfer" problem. However, the emerging deep transfer models do not fit well with most existing data selection methods, because the data selection policy and the transfer learning model are not jointly trained, leading to sub-optimal training efficiency. In this paper, we propose a novel reinforced data selector to select high-quality source domain data to help the TL model. Specifically, the data selector "acts" on the source domain data to find a subset for optimization of the TL model, and the performance of the TL model can provide "rewards" in turn to update the selector. We build the reinforced data selector based on the actor-critic framework and integrate it to a DNN based transfer learning model, resulting in a Reinforced Transfer Learning (RTL) method. We perform a thorough experimental evaluation on two major tasks for text matching, namely, paraphrase identification and natural language inference. Experimental results show the proposed RTL can significantly improve the performance of the TL model. We further investigate different settings of states, rewards, and policy optimization methods to examine the robustness of our method. Last, we conduct a case study on the selected data and find our method is able to select source domain data whose Wasserstein distance is close to the target domain data. This is reasonable and intuitive as such source domain data can provide more transferability power to the model.


title: Expanders - how to find them, and what to find in them

abstract: A graph $G=(V,E)$ is called an expander if every vertex subset $U$ of size up to $|V|/2$ has an external neighborhood whose size is comparable to $|U|$. Expanders have been a subject of intensive research for more than three decades and have become one of the central notions of modern graph theory. We first discuss the above definition of an expander and its alternatives. Then we present examples of families of expanding graphs and state basic properties of expanders. Next, we introduce a way to argue that a given graph contains a large expanding subgraph. Finally we research properties of expanding graphs, such as existence of small separators, of cycles (including cycle lengths), and embedding of large minors.


title: Thermodynamic equilibrium of binary mixtures on curved surfaces

abstract: We study the global influence of curvature on the free energy landscape of two-dimensional binary mixtures confined on closed surfaces. Starting from a generic effective free energy, constructed on the basis of symmetry considerations and conservation laws, we identify several model-independent phenomena, such as a curvature-dependent line tension and local shifts in the binodal concentrations. To shed light on the origin of the phenomenological parameters appearing in the effective free energy, we further construct a lattice-gas model of binary mixtures on non-trivial substrates, based on the curved-space generalization of the two-dimensional Ising model. This allows us to decompose the interaction between the local concentration of the mixture and the substrate curvature into four distinct contributions, as a result of which the phase diagram splits into critical sub-diagrams. The resulting free energy landscape can admit, as stable equilibria, strongly inhomogeneous mixed phases, which we refer to as antimixed states below the critical temperature. We corroborate our semi-analytical findings with phase-field numerical simulations on realistic curved lattices. Despite this work being primarily motivated by recent experimental observations of multi-component lipid vesicles supported by colloidal scaffolds, our results are applicable to any binary mixture confined on closed surfaces of arbitrary geometry.


title: Spectral methods for testing cluster structure of graphs

abstract: In the framework of graph property testing, we study the problem of determining if a graph admits a cluster structure. We say that a graph is $(k, \phi)$-clusterable if it can be partitioned into at most $k$ parts such that each part has conductance at least $\phi$. We present an algorithm that accepts all graphs that are $(2, \phi)$-clusterable with probability at least $\frac{2}3$ and rejects all graphs that are $\epsilon$-far from $(2, \phi^*)$-clusterable for $\phi^* \le \mu \phi^2 \epsilon^2$ with probability at least $\frac{2}3$ where $\mu > 0$ is a parameter that affects the query complexity. This improves upon the work of Czumaj, Peng, and Sohler by removing a $\log n$ factor from the denominator of the bound on $\phi^*$ for the case of $k=2$. Our work was concurrent with the work of Chiplunkar et al.\@ who achieved the same improvement for all values of $k$. Our approach for the case $k=2$ relies on the geometric structure of the eigenvectors of the graph Laplacian and results in an algorithm with query complexity $O(n^{1/2+O(1)\mu} \cdot \text{poly}(1/\epsilon, 1/\phi,\log n))$.


title: Augmented data for the Backus problem

abstract: Backus considered a boundary value problem for the Laplace equation with the non-linear data in the form of the magnitude $|Du|$ of the gradient of the solution u. We consider this problem with the data expanded by $(\partial/\partial\nu)|Du|$ given on the boundary of the domain. To justify the requirement for additional data, we use them to estimate the number of sources for the related inverse source problem in the plane. We show that, for an arbitrary dimension, a harmonic function satisfies a quasi-linear equation on the boundary of the domain with the coefficients involving the augmented data. We use the finite element method to recover the harmonic function on the boundary by solving numerically the derived equation.


title: Finite-dimensional pointed Hopf algebras over finite simple groups of Lie type V. Mixed classes in Chevalley and Steinberg groups

abstract: We show that all classes that are neither semisimple nor unipotent in finite simple Chevalley or Steinberg groups different from $PSL_n(q)$ collapse (i.e. are never the support of a finite-dimensional Nichols algebra). As a consequence, we prove that the only finite-dimensional pointed Hopf algebra whose group of group-like elements is $PSp_{2n}(q)$, $P\Omega^+_{4n}(q)$, $P\Omega^-_{4n}(q)$, $^3D_4(q)$, $E_7(q)$, $E_8(q)$, $F_4(q)$, or $G_2(q)$ with q even is the group algebra.


title: Metric Regularity of Quasidifferentiable Mappings and Optimality Conditions for Nonsmooth Mathematical Programming Problems

abstract: This article is devoted to the analysis of necessary and/or sufficient conditions for metric regularity in terms of Demyanov-Rubinov-Polyakova quasidifferentials. We obtain new necessary and sufficient conditions for the local metric regularity of a multifunction in terms of quasidifferentials of the distance function to this multifunction. We also propose a new MFCQ-type constraint qualification for a parametric system of quasidifferentiable equality and inequality constraints and prove that it ensures the metric regularity of a multifunction associated with this system. As an application, we utilize our constraint qualification to strengthen existing optimality conditions for quasidifferentiable programming problems with equality and inequality constraints. We also prove the independence of the optimality conditions of the choice of quasidifferentials and present a simple example in which the optimality conditions in terms of quasidifferentials detect the non-optimality of a given point, while optimality conditions in terms of various subdifferentials fail to disqualify this point as non-optimal.


title: Search for top-quark decays $t \rightarrow Hq$ with 36 fb$^{-1}$ of $pp$ collision data at $\sqrt{s}=13$ TeV with the ATLAS detector

abstract: A search for flavour-changing neutral current decays of a top quark into an up-type quark ($q=u, c$) and the Standard Model Higgs boson, $t\rightarrow Hq$, is presented. The search is based on a dataset of $pp$ collisions at $\sqrt{s}=13$ TeV recorded in 2015 and 2016 with the ATLAS detector at the CERN Large Hadron Collider and corresponds to an integrated luminosity of 36.1 fb$^{-1}$. Two complementary analyses are performed that search for top-quark pair events in which one top quark decays into $Wb$ and the other top quark decays into $Hq$, and target the $H\rightarrow b\bar{b}$ and $H\rightarrow\tau^+\tau^-$ decay modes, respectively. The high multiplicity of $b$-quark jets, or the presence of hadronically decaying $\tau$-leptons, are exploited in the two analyses respectively. Multivariate techniques are used to separate the signal from the background, which is dominated by top-quark pair production. No significant excess of events above the background expectation is found, and 95% CL upper limits on the $t\rightarrow Hq$ branching ratios are derived. The combination of these searches with ATLAS searches in diphoton and multilepton final states yields observed (expected) 95% CL upper limits on the $t\rightarrow Hc$ and $t\rightarrow Hu$ branching ratios of $1.1 \times 10^{-3}$ ($8.3 \times 10^{-4}$) and $1.2 \times 10^{-3}$ ($8.3 \times 10^{-4}$), respectively. The corresponding combined observed (expected) upper limits on the $|\lambda_{tcH}|$ and $|\lambda_{tuH}|$ couplings are 0.064 (0.055) and 0.066 (0.055), respectively.


title: Numerical simulation of compressible flows using two-phase observable Euler equations

abstract: Most fluid flow problems that are vital in engineering applications involve at least one of the following features: turbulence, shocks, and/or material interfaces. While seemingly different phenomena, these flows all share continuous generation of high wavenumber modes, which we term the '$k_\infty$ irregularity.' In this work, an inviscid regularization technique called 'observable regularization' is proposed for the simulation of two-phase compressible flows. The proposed approach regularizes the equations at the level of the partial differential equation and as a result, any numerical method can be used to solve the system of equations. The regularization is accomplished by introducing an 'observability limit' that represents the length scale below which one cannot properly model or continue to resolve flow structures. An observable volume fraction equation is derived for capturing the material interface, which satisfies the pressure equilibrium at the interface. The efficacy of the observable regularization method is demonstrated using several test cases, including a one-dimensional material interface tracking, one-dimensional shock-tube and shock-bubble problems, and two-dimensional simulations of a shock interacting with a cylindrical bubble. The results show favorable agreement, both qualitatively and quantitatively, with available exact solutions or numerical and experimental data from the literature. The computational saving by using the current method is estimated to be about one order of magnitude in two-dimensional computations and significantly higher in three-dimensional computations. Lastly, the effect of the observability limit and best practices to choose its value are discussed.


title: Effects of Violation of Equivalence Principle on UHE Neutrinos at IceCube in 4 Flavour Scenario

abstract: If weak equivalence principle is violated then different types of neutrinos would couple differently with gravity and that may produce a gravity induced oscillation for the neutrinos of different flavour. We explore here the possibility that very small violation of the principle of weak equivalence (VEP) can be probed by ultra high energy neutrinos from distant astrophysical sources. The very long baseline length and the ultra high energies of such neutrinos could be helpful to probe very small VEP. We consider a 4-flavour neutrino scenario (3 active + 1 sterile) with both mass-flavour and gravity induced oscillations and compare the detection signatures for these neutrinos (muon tracks and shower events) with and without gravity induced oscillations at a kilometer scale detector such as IceCube. We find that even very small VEP ($\sim 10^{-42}$) can considerably affect the detected muon yield produced by UHE neutrinos from distant Gamma Ray Bursts (GRBs)


title: A Formal Separation Between Strategic and Nonstrategic Behavior

abstract: It is common to make a distinction between "strategic" behavior and other forms of intentional but "nonstrategic" behavior: typically, that strategic agents model other agents while nonstrategic agents do not. However, a crisp boundary between these concepts has proven elusive. This problem is pervasive throughout the game theoretic literature on bounded rationality and particularly critical in parts of the behavioral game theory literature that make an explicit distinction between the behavior of "nonstrategic" level-0 agents and "strategic" higher-level agents (e.g., the level-k and cognitive hierarchy models). Overall, work discussing bounded rationality rarely gives clear guidance on how the rationality of nonstrategic agents must be bounded, instead typically just singling out specific decision rules (e.g., randomizing uniformly, playing toward the best case, optimizing the worst case) and informally asserting that they are nonstrategic. In this work, we propose a new, formal characterization of nonstrategic behavior. Our main contribution is to show that it satisfies two properties: (1) it is general enough to capture all purportedly "nonstrategic" decision rules of which we are aware in the behavioral game theory literature; (2) behavior that obeys our characterization is distinct from strategic behavior in a precise sense.


title: Van der Waals equation of state and PVT properties of real fluid

abstract: It is shown that: in the case when two parameters of the Van der Waals equation of state are defined from the critical temperature and pressure the exact parametrical solution of the equations of the liquid-vapor phase equilibrium of the Van der Waals fluid quantitatively describes the experimental dependencies of the saturated pressure of argon on the temperature and reduced vapor density, and it gives the quantitative description of the temperature dependencies of the reduced densities near critical point. When the parameters are defined from the critical pressure and density the parametric solution describes quantitatively the experimental dependencies of the saturated pressure of argon on the density and reduced temperature, it can describe qualitatively the dependencies of the vapor and liquid densities on the reduced temperature, and it gives the quantitative description of the dependencies of the densities on the reduced temperature near critical point. If the parameters are defined from the critical temperature and density then the exact solution describes quantitatively the experimental dependencies of the reduced saturated pressure of argon on the density and temperature, it describes qualitatively the temperature dependencies of the vapor and liquid densities of argon, and it gives the quantitative description of the temperature dependencies of the vapor and liquid densities near critical point. It is also shown that the Van der Waals equation of state describes quantitatively the reference experimental PVT- data for the gas and supercritical fluid states for the under-critical densities of argon, the dependencies of the saturation pressure on the temperature and vapor density, and the dependence of the vapor density of argon on temperature if the parameters are defined from the critical pressure and temperature.


title: A Probabilistic and Non-Deterministic Call-by-Push-Value Language

abstract: There is no known way of giving a domain-theoretic semantics to higher-order probabilistic languages, in such a way that the involved domains are continuous or quasi-continuous - the latter is required to do any serious mathematics. We argue that the problem naturally disappears for languages with two kinds of types, where one kind is interpreted in a Cartesian-closed category of continuous dcpos, and the other is interpreted in a category that is closed under the probabilistic powerdomain functor. Such a setting is provided by Paul B. Levy's call-by-push-value paradigm. Following this insight, we define a call-by-push-value language, with probabilistic choice sitting inside the value types, and where conversion from a value type to a computation type involves demonic non-determinism. We give both a domain-theoretic semantics and an operational semantics for the resulting language, and we show that they are sound and adequate. With the addition of statistical termination testers and parallel if, we show that the language is even fully abstract - and those two primitives are required for that.


title: Fingerprint Presentation Attack Detection: Generalization and Efficiency

abstract: We study the problem of fingerprint presentation attack detection (PAD) under unknown PA materials not seen during PAD training. A dataset of 5,743 bonafide and 4,912 PA images of 12 different materials is used to evaluate a state-of-the-art PAD, namely Fingerprint Spoof Buster. We utilize 3D t-SNE visualization and clustering of material characteristics to identify a representative set of PA materials that cover most of PA feature space. We observe that a set of six PA materials, namely Silicone, 2D Paper, Play Doh, Gelatin, Latex Body Paint and Monster Liquid Latex provide a good representative set that should be included in training to achieve generalization of PAD. We also propose an optimized Android app of Fingerprint Spoof Buster that can run on a commodity smartphone (Xiaomi Redmi Note 4) without a significant drop in PAD performance (from TDR = 95.7% to 95.3% @ FDR = 0.2%) which can make a PA prediction in less than 300ms.


title: Topological Quantum Field Theories for character varieties

abstract: This PhD Thesis is devoted to the study of Hodge structures on a special type of complex algebraic varieties, the so-called character varieties. For this purpose, we propose to use a powerful algebro-geometric tool coming from theoretical physics, known as Topological Quantum Field Theory (TQFT). With this idea in mind, in the present Thesis we develop a formalism that allows us to construct TQFTs from two simpler pieces of data: a field theory (geometric data) and a quantisation (algebraic data). As an application, we construct a TQFT computing Hodge structures on representation varieties and we use it for computing explicity the Deligne-Hodge polynomials of parabolic $\mathrm{SL}_2(\mathbb{C})$-character varieties.


title: Lattice of the dual of an Anderson $t$-motive as an element of the maximal Schubert cycle of its flag variety

abstract: Let $M$ be an uniformizable Anderson $t$-motive. We can describe its lattice as an element of the maximal Schubert cycle of its flag variety. This element (called a Siegel element) is a set of matrices parametrized by integer points of a tetrahedron. We give explicit formulas for a Siegel element for $M'$ -- the dual of $M$. As a by-product, we get another proof of the theorem that the lattice of $M'$ is the dual of the lattice of $M$, independent of the proof obtained by U. Hartl and A.-K. Juschka. Generalizations of this result to non-maximal Schubert cycles, other tensor operations (tensor product and Hom) and t-motives having non-trivial endomorphism rings, are subjects of further research.


title: Robust random search with scale-free stochastic resetting

abstract: A new model of search based on stochastic resetting is introduced, wherein rate of resets depends explicitly on time elapsed since the beginning of the process. It is shown that rate inversely proportional to time leads to paradoxical diffusion which mixes self-similarity and linear growth of the mean square displacement with non-locality and non-Gaussian propagator. It is argued that such resetting protocol offers a general and efficient search-boosting method that does not need to be optimized with respect to the scale of the underlying search problem (e.g., distance to the goal) and is not very sensitive to other search parameters. Both subdiffusive and superdiffusive regimes of the mean squared displacement scaling are demonstrated with more general rate functions.


title: Space Telescope and Optical Reverberation Mapping Project. X. Understanding the Absorption-Line Holiday in NGC 5548

abstract: The flux variations in the emission lines in active galactic nuclei (AGNs) are driven by variations in the ionizing continuum flux --which are usually reflected in the observable UV-optical continuum. The "Reverberation mapping" technique measures the delay between line and continuum variations to determine the size of the line emitting region, this is the basis for measurements of the central black hole mass in AGNs. The Space Telescope and Optical Reverberation Mapping Project (AGN STORM) on NGC 5548 in 2014 is the most intensive multi-wavelength AGN monitoring campaign ever. For most of the campaign, the emission-line variations followed changes in the continuum with a time lag, as expected. However, the lines varied independently of the observed UV-optical continuum during a 60 -- 70 day "holiday." To understand this remarkable phenomenon, we study the intrinsic absorption lines present in NGC 5548. We identify a novel cycle that reproduces the absorption line variability and thus identify the physics that allows the holiday to occur. In our model, variations in this obscurer's line-of-sight covering factor modify the soft X-ray continuum. This leads to changes in the ionization of helium gas in the broad-line region. Ionizing radiation produced by recombining helium then affects the ionization of other species as observed during the AGN STORM holiday. It is likely that any other model which selectively changes the soft X-ray part of the continuum during the holiday can also explain the anomalous emission line behavior observed.


title: Dynamics of Magnetoelectric Reversal of Antiferromagnetic Domain

abstract: When electric and magnetic fields are applied together on a magnetoelectric antiferromagnet, the domain state is subject to reversal. Although the initial and final conditions are saturated single-domain states, the process of reversal may decompose into local multi-domain switching events. In thin films of Cr2O3, the magnetoelectric coercivity and the switching speed found from experiments are considerably lower than expected from magnetic anisotropy, similar to Brown's paradox in ferromagnetic materials. Multi-domain effects originate because antiferromagnetic domain walls are metastably pinned by lattice defects, not due to reduction of magnetostatic energy, which is negligible. This paper theoretically analyzes domain reversal in thin-film magnetoelectric antiferromagnets in the form of nucleation, domain wall propagation, and coherent rotation. The timescales of reversal mechanisms are modeled as a function of applied magnetoelectric pressure. The theory is assessed with reference to latest experimental works on magnetoelectric switching of thin-film Cr2O3: domain wall propagation is found to be dominant and responsible for switching in the experiments. The results bear implications in the energy-delay performance of ME memory devices utilizing antiferromagnetic insulators, which are prospective for nonvolatile technology.


title: Relationship between quandle shadow cocycle invariants and Vassiliev invariants of links

abstract: We prove that there exist links and quandle (shadow) cocycle invariants by using Alexander quandles such that quandle (shadow) cocycle invariants by using Alexander quandles induce Vassiliev invariants. It is natural to try to relate quandle (shadow) cocycle invariants to quantum invariants.


title: Unsupervised learning by a nonlinear network with Hebbian excitatory and anti-Hebbian inhibitory neurons

abstract: This paper introduces a rate-based nonlinear neural network in which excitatory (E) neurons receive feedforward excitation from sensory (S) neurons, and inhibit each other through disynaptic pathways mediated by inhibitory (I) interneurons. Correlation-based plasticity of disynaptic inhibition serves to incompletely decorrelate E neuron activity, pushing the E neurons to learn distinct sensory features. The plasticity equations additionally contain "extra" terms fostering competition between excitatory synapses converging onto the same postsynaptic neuron and inhibitory synapses diverging from the same presynaptic neuron. The parameters of competition between S$\to$E connections can be adjusted to make learned features look more like "parts" or "wholes." The parameters of competition between I-E connections can be adjusted to set the typical decorrelatedness and sparsity of E neuron activity. Numerical simulations of unsupervised learning show that relatively few I neurons can be sufficient for achieving good decorrelation, and increasing the number of I neurons makes decorrelation more complete. Excitatory and inhibitory inputs to active E neurons are approximately balanced as a result of learning.


title: Classical ladder functions for Rosen-Morse and curved Kepler-Coulomb systems

abstract: Ladder functions in classical mechanics are defined in a similar way as ladder operators in the context of quantum mechanics. In the present paper, we develop a new method for obtaining ladder functions of one dimensional systems by means of a product of two `factor functions'. We apply this method to the curved Kepler-Coulomb and Rosen-Morse II systems whose ladder functions were not found yet. The ladder functions here obtained are applied to get the motion of the system.


title: A Gramian Description of the Degree 4 Generalized Elliptope

abstract: One of the most widely studied convex relaxations in combinatorial optimization is the relaxation of the cut polytope $\mathscr C^N$ to the elliptope $\mathscr E^N$, which corresponds to the degree 2 sum-of-squares (SOS) relaxation of optimizing a quadratic form over the hypercube $\{\pm 1\}^N$. We study the extension of this classical idea to degree 4 SOS, which gives an intermediate relaxation we call the degree 4 generalized elliptope $\mathscr E_4^N$. Our main result is a necessary and sufficient condition for the Gram matrix of a collection of vectors to belong to $\mathscr E_4^N$. Consequences include a tight rank inequality between degree 2 and degree 4 pseudomoment matrices, and a guarantee that the only extreme points of $\mathscr E^N$ also in $\mathscr E_4^N$ are the cut matrices; that is, $\mathscr E^N$ and $\mathscr E_4^N$ share no "spurious" extreme point. For Gram matrices of equiangular tight frames, we give a simple criterion for membership in $\mathscr{E}_4^N$. This yields new inequalities satisfied in $\mathscr{E}_4^N$ but not $\mathscr{E}^N$ whose structure is related to the Schl\"{a}fli graph and which cannot be obtained as linear combinations of triangle inequalities. We also give a new proof of the restriction to degree 4 of a result of Laurent showing that $\mathscr{E}_4^N$ does not satisfy certain cut polytope inequalities capturing parity constraints. Though limited to this special case, our proof of the positive semidefiniteness of Laurent's pseudomoment matrix is short and elementary. Our techniques also suggest that membership in $\mathscr{E}_4^N$ is closely related to the partial transpose operation on block matrices, which has previously played an important role in the study of quantum entanglement. To illustrate, we present a correspondence between certain entangled bipartite quantum states and the matrices of $\mathscr{E}_4^N\setminus\mathscr{C}^N$.


title: Successive Wyner-Ziv Coding for the Binary CEO Problem under Logarithmic Loss

abstract: The $L$-link binary Chief Executive Officer (CEO) problem under logarithmic loss is investigated in this paper. A quantization splitting technique is applied to convert the problem under consideration to a $(2L-1)$-step successive Wyner-Ziv (WZ) problem, for which a practical coding scheme is proposed. In the proposed scheme, low-density generator-matrix (LDGM) codes are used for binary quantization while low-density parity-check (LDPC) codes are used for syndrome generation; the decoder performs successive decoding based on the received syndromes and produces a soft reconstruction of the remote source. The simulation results indicate that the rate-distortion performance of the proposed scheme can approach the theoretical inner bound based on binary-symmetric test-channel models.


title: The NANOGrav 11-Year Data Set: Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries

abstract: Observations indicate that nearly all galaxies contain supermassive black holes (SMBHs) at their centers. When galaxies merge, their component black holes form SMBH binaries (SMBHBs), which emit low-frequency gravitational waves (GWs) that can be detected by pulsar timing arrays (PTAs). We have searched the recently-released North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11-year data set for GWs from individual SMBHBs in circular orbits. As we did not find strong evidence for GWs in our data, we placed 95\% upper limits on the strength of GWs from such sources as a function of GW frequency and sky location. We placed a sky-averaged upper limit on the GW strain of $h_0 < 7.3(3) \times 10^{-15}$ at $f_\mathrm{gw}= 8$ nHz. We also developed a technique to determine the significance of a particular signal in each pulsar using ``dropout' parameters as a way of identifying spurious signals in measurements from individual pulsars. We used our upper limits on the GW strain to place lower limits on the distances to individual SMBHBs. At the most-sensitive sky location, we ruled out SMBHBs emitting GWs with $f_\mathrm{gw}= 8$ nHz within 120 Mpc for $\mathcal{M} = 10^9 \, M_\odot$, and within 5.5 Gpc for $\mathcal{M} = 10^{10} \, M_\odot$. We also determined that there are no SMBHBs with $\mathcal{M} > 1.6 \times 10^9 \, M_\odot$ emitting GWs in the Virgo Cluster. Finally, we estimated the number of potentially detectable sources given our current strain upper limits based on galaxies in Two Micron All-Sky Survey (2MASS) and merger rates from the Illustris cosmological simulation project. Only 34 out of 75,000 realizations of the local Universe contained a detectable source, from which we concluded it was unsurprising that we did not detect any individual sources given our current sensitivity to GWs.


title: Leishmaniasis Parasite Segmentation and Classification using Deep Learning

abstract: Leishmaniasis is considered a neglected disease that causes thousands of deaths annually in some tropical and subtropical countries. There are various techniques to diagnose leishmaniasis of which manual microscopy is considered to be the gold standard. There is a need for the development of automatic techniques that are able to detect parasites in a robust and unsupervised manner. In this paper we present a procedure for automatizing the detection process based on a deep learning approach. We train a U-net model that successfully segments leismania parasites and classifies them into promastigotes, amastigotes and adhered parasites.


title: Sentiment Classification of Customer Reviews about Automobiles in Roman Urdu

abstract: Text mining is a broad field having sentiment mining as its important constituent in which we try to deduce the behavior of people towards a specific item, merchandise, politics, sports, social media comments, review sites etc. Out of many issues in sentiment mining, analysis and classification, one major issue is that the reviews and comments can be in different languages like English, Arabic, Urdu etc. Handling each language according to its rules is a difficult task. A lot of research work has been done in English Language for sentiment analysis and classification but limited sentiment analysis work is being carried out on other regional languages like Arabic, Urdu and Hindi. In this paper, Waikato Environment for Knowledge Analysis (WEKA) is used as a platform to execute different classification models for text classification of Roman Urdu text. Reviews dataset has been scrapped from different automobiles sites. These extracted Roman Urdu reviews, containing 1000 positive and 1000 negative reviews, are then saved in WEKA attribute-relation file format (arff) as labeled examples. Training is done on 80% of this data and rest of it is used for testing purpose which is done using different models and results are analyzed in each case. The results show that Multinomial Naive Bayes outperformed Bagging, Deep Neural Network, Decision Tree, Random Forest, AdaBoost, k-NN and SVM Classifiers in terms of more accuracy, precision, recall and F-measure.


title: Cascaded V-Net using ROI masks for brain tumor segmentation

abstract: In this work we approach the brain tumor segmentation problem with a cascade of two CNNs inspired in the V-Net architecture \cite{VNet}, reformulating residual connections and making use of ROI masks to constrain the networks to train only on relevant voxels. This architecture allows dense training on problems with highly skewed class distributions, such as brain tumor segmentation, by focusing training only on the vecinity of the tumor area. We report results on BraTS2017 Training and Validation sets.


title: Relatively complicated? Using models to teach general relativity at different levels

abstract: This review presents an overview of various kinds of models -- physical, abstract, mathematical, visual -- that can be used to present the concepts and applications of Einstein's general theory of relativity at the level of undergraduate and even high-school teaching. After a general introduction dealing with various kinds of models and their properties, specific areas of general relativity are addressed: the elastic sheet model and other models for the fundamental geometric properties of gravity, models for black holes including the river model, cosmological models for an expanding universe, and models for gravitational waves as well as for interferometric gravitational wave detectors.


title: A classification of finite simple amenable ${\cal Z}$-stable $C^*$-algebras, I: $C^*$-algebras with generalized tracial rank one

abstract: A class of $C^*$-algebras, to be called those of generalized tracial rank one, is introduced, and classified by the Elliott invariant. A second class of unital simple separable amenable $C^*$-algebras, those whose tensor products with UHF-algebras of infinite type are in the first class, to be referred to as those of rational generalized tracial rank one, is proved to exhaust all possible values of the Elliott invariant for unital finite simple separable amenable ${\cal Z}$-stable $C^*$-algebras. An isomorphism theorem for a special sub-class of those $C^*$-algebras are presented. This provides the basis for the classification of $C^*$-algebras with rational generalized tracial rank one in Part II.


title: Continuous quantum measurement and It\^o formalism

abstract: A new quantum-stochastic differential calculus is derived for representing continuous quantum measurement of the position operator. Closed nonlinear quantum-stochastic differential equation is given for the quantum state of the observed particle. Possible physical realization of continuous position measurement is constructed.


title: A Geometric Theory of Higher-Order Automatic Differentiation

abstract: First-order automatic differentiation is a ubiquitous tool across statistics, machine learning, and computer science. Higher-order implementations of automatic differentiation, however, have yet to realize the same utility. In this paper I derive a comprehensive, differential geometric treatment of automatic differentiation that naturally identifies the higher-order differential operators amenable to automatic differentiation as well as explicit procedures that provide a scaffolding for high-performance implementations.


title: Simple bounded highest weight modules of basic classical Lie superalgebras

abstract: We classify all simple bounded highest weight modules of a basic classical Lie superalgebra $\mathfrak g$. In particular, our classification leads to the classification of the simple weight modules with finite weight multiplicities over all classical Lie superalgebras. We also obtain some character formulas of strongly typical bounded highest weight modules of $\mathfrak g$.


title: Joint Overlap Analysis of Multiple Genomic Interval Sets

abstract: Next-generation sequencing (NGS) technologies have produced large volumes of genomic data. One common operation on heterogeneous genomic data is genomic interval intersection. Most of the existing tools impose restrictions such as not allowing nested intervals or requiring intervals to be sorted when finding overlaps in two or more interval sets. We proposed segment tree (ST) and indexed segment tree forest (ISTF) based solutions for intersection of multiple genomic interval sets in parallel. We developed these methods as a tool, Joint Overlap Analysis (JOA), which takes n interval sets and finds overlapping intervals with no constraints on the given intervals. The proposed indexed segment tree forest is a novel composite data structure, which leverages on indexing and natural binning of a segment tree. We also presented construction and search algorithms for this novel data structure. We compared JOA ST and JOA ISTF with each other, and with other interval intersection tools for verification of its correctness and for showing that it attains comparable execution times. We implemented JOA in Java using the fork/join framework which speeds up parallel processing by taking advantage of all available processor cores. We compared JOA ST with JOA ISTF and showed that segment tree and indexed segment tree forest methods are comparable with each other in terms of execution time and memory usage. We also carried out execution time comparison analysis for JOA and other tools and demonstrated that JOA has comparable execution time and is able to further reduce its running time by using more processors per node. JOA can be run using its GUI or as a command line tool. JOA is available with source code at https://github.com/burcakotlu/JOA/. A user manual is provided at https://joa.readthedocs.org


title: Julia sets of Newton maps of real quadratic polynomial maps on the plane

abstract: We study numerically the $\alpha$- and $\omega$-limits of the Newton maps of two of the most elementary families of polynomial transformations on the plane: those with a linear component and those with both components of degree two. Our results are fully consistent with some conjectures we posed in a recent work about the dynamics of Newton maps.


title: Towards a CAN IDS based on a neural-network data field predictor

abstract: Modern vehicles contain a few controller area networks (CANs), which allow scores of on-board electronic control units (ECUs) to communicate messages critical to vehicle functions and driver safety. CAN provide a lightweight and reliable broadcast protocol but is bereft of security features. As evidenced by many recent research works, CAN exploits are possible both remotely and with direct access, fueling a growing CAN intrusion detection system (IDS) body of research. A challenge for pioneering vehicle-agnostic IDSs is that passenger vehicles' CAN message encodings are proprietary, defined and held secret by original equipment manufacturers (OEMs). Targeting detection of next-generation attacks, in which messages are sent from the expected ECU at the expected time but with malicious content, researchers are now seeking to leverage "CAN data models", which predict future CAN message contents and use prediction error to identify anomalous, hopefully malicious CAN messages. Yet, current works model CAN signals post-translation, i.e., after applying OEM-donated or reverse-engineered translations from raw data. In this work, we present initial IDS results testing deep neural networks used to predict CAN data at the bit level, thereby providing IDS capabilities but avoiding reverse engineering proprietary encodings. Our results suggest the method is promising for continuous signals in CAN data, but struggles for discrete, e.g., binary, signals.


title: Spacetime Symmetries, Invariant Sets, and Additive Subdynamics of Cellular Automata

abstract: Cellular automata are fully-discrete, spatially-extended dynamical systems that evolve by simultaneously applying a local update function. Despite their simplicity, the induced global dynamic produces a stunning array of richly-structured, complex behaviors. These behaviors present a challenge to traditional closed-form analytic methods. In certain cases, specifically when the local update is additive, powerful techniques may be brought to bear, including characteristic polynomials, the ergodic theorem with Fourier analysis, and endomorphisms of compact Abelian groups. For general dynamics, though, where such analytics generically do not apply, behavior-driven analysis shows great promise in directly monitoring the emergence of structure and complexity in cellular automata. Here we detail a surprising connection between generalized symmetries in the spacetime fields of configuration orbits as revealed by the behavior-driven local causal states, invariant sets of spatial configurations, and additive subdynamics which allow for closed-form analytic methods.

categories: cond-mat.stat-mech cs.FL nlin.CG nlin.PS


title: Salvaging Falsified Instrumental Variable Models

abstract: What should researchers do when their baseline model is refuted? We provide four constructive answers. First, researchers can measure the extent of falsification. To do this, we consider continuous relaxations of the baseline assumptions of concern. We then define the falsification frontier: The smallest relaxations of the baseline model which are not refuted. This frontier provides a quantitative measure of the extent of falsification. Second, researchers can present the identified set for the parameter of interest under the assumption that the true model lies somewhere on this frontier. We call this the falsification adaptive set. This set generalizes the standard baseline estimand to account for possible falsification. Third, researchers can present the identified set for a specific point on this frontier. Finally, as a sensitivity analysis, researchers can present identified sets for points beyond the frontier. To illustrate these four ways of salvaging falsified models, we study overidentifying restrictions in two instrumental variable models: a homogeneous effects linear model, and heterogeneous effect models with either binary or continuous outcomes. In the linear model, we consider the classical overidentifying restrictions implied when multiple instruments are observed. We generalize these conditions by considering continuous relaxations of the classical exclusion restrictions. By sufficiently weakening the assumptions, a falsified baseline model becomes non-falsified. We obtain analogous results in the heterogeneous effect models, where we derive identified sets for marginal distributions of potential outcomes, falsification frontiers, and falsification adaptive sets under continuous relaxations of the instrument exogeneity assumptions. We illustrate our results in four different empirical applications.


title: Counting integers representable as images of polynomials modulo $n$

abstract: Given a polynomial $f(x_1,x_2,\ldots, x_t)$ in $t$ variables with integer coefficients and a positive integer $n$, let $\alpha(n)$ be the number of integers $0\leq a<n$ such that the polynomial congruence $f(x_1, x_2, \ldots, x_t)\equiv a\ (mod\ n)$ is solvable. We describe a method that allows to determine the function $\alpha$ associated to polynomials of the form $c_1x_1^k+c_2x_2^k+\cdots+c_tx_t^k$. Then we apply this method to polynomials that involve sums and differences of squares, mainly to the polynomials $x^2+y^2, x^2-y^2$ and $x^2+y^2+z^2$.


title: Constrained Inverse Optimal Control with Application to a Human Manipulation Task

abstract: This paper presents an inverse optimal control methodology and its application to training a predictive model of human motor control from a manipulation task. It introduces a convex formulation for learning both objective function and constraints of an infinite-horizon constrained optimal control problem with nonlinear system dynamics. The inverse approach utilizes Bellman's principle of optimality to formulate the infinite-horizon optimal control problem as a shortest path problem and Lagrange multipliers to identify constraints. We highlight the key benefit of using the shortest path formulation, i.e., the possibility of training the predictive model with short and selected trajectory segments. The method is applied to training a predictive model of movements of a human subject from a manipulation task. The study indicates that individual human movements can be predicted with low error using an infinite-horizon optimal control problem with constraints on shoulder movement.


title: Allocation strategies for high fidelity models in the multifidelity regime

abstract: We propose a novel approach to allocating resources for expensive simulations of high fidelity models when used in a multifidelity framework. Allocation decisions that distribute computational resources across several simulation models become extremely important in situations where only a small number of expensive high fidelity simulations can be run. We identify this allocation decision as a problem in optimal subset selection, and subsequently regularize this problem so that solutions can be computed. Our regularized formulation yields a type of group lasso problem that has been studied in the literature to accomplish subset selection. Our numerical results compare performance of algorithms that solve the group lasso problem for algorithmic allocation against a variety of other strategies, including those based on classical linear algebraic pivoting routines and those derived from more modern machine learning-based methods. We demonstrate on well known synthetic problems and more difficult real-world simulations that this group lasso solution to the relaxed optimal subset selection problem performs better than the alternatives.


title: Circuit optimization for IBM processors: A way to get higher fidelity and higher values of nonclassicality witnesses

abstract: Recently, various quantum computing and communication tasks have been implemented using IBM's superconductivity-based quantum computers which are available on the cloud. Here, we show that the circuits used in most of those works were not optimized and the use of the optimized circuits can considerably improve the possibility of observing unique features of quantum mechanics. Specifically, a systematic procedure is used here to obtain optimized circuits (circuits having reduced gate count and number of levels) for a large number of Clifford+T circuits which have already been implemented in the IBM quantum computers. Optimized circuits implementable in IBM quantum computers are also obtained for a set of reversible benchmark circuits. With a clear example, it is shown that the reduction in circuit costs enhances the fidelity of the output state (with respect to the theoretically expected state in the absence of noise) as lesser number of gates and levels introduce lesser amount of errors during evolution of the state. Further, considering Mermin inequality as an example, it's shown that the violation of classical limit is enhanced when we use an optimized circuit. Thus, the approach adopted here can be used to identify relatively weaker signature of quantumness and also to establish quantum supremacy in a stronger manner.


title: First order perturbations of hypersurfaces of arbitrary causal character

abstract: In this work we study the problem of first order perturbations of a general hypersurface, i.e. with arbitrary causal character at each point. We extend the framework by Mars (Class. Quantum Grav. 22 3325 (2005)) where this problem was studied to second order for everywhere timelike or spacelike hypersurfaces, and we adapt it to cover the general case. We apply the formalism to the matching of spacetimes across a general hypersurface to first order in perturbation theory.


title: Neutrinos Angra experiment: commissioning and first operational measurements

abstract: The Neutrinos Angra Experiment has completed a major step by finishing the comissioning of the detector and the data acquisition system at the experimental site located in the Angra dos Reis nuclear power plant. The experiment consists of a water-based detector and associated electronics, both designed with the goal of detecting the electron antineutrinos produced by the nuclear reactor. The detection is possible due to the Inverse Beta Decay, where the final products in the water are photons in the UV-to-visible range of the spectrum. The assembled detector comprises three active volumes filled with water: (i) a cubic target detector for electron antineutrinos, covered by 32 8-inches PMTs, (ii) a lateral layer surrounding the target (veto) equipped with 4 PMTs and (iii) a third volume covering the top of both, also equipped with 4~PMTs. In the present document the main features of the detector assembly as well as the integration of the readout electronics on-site are reported. Finally, some operational characteristics are shown based on straightforward analysis of the first measurements performed during the last months with the fully working detector.


title: Geometrical and statistical properties of M-estimates of scatter on Grassmann manifolds

abstract: We consider data from the Grassmann manifold $G(m,r)$ of all vector subspaces of dimension $r$ of $\mathbb{R}^m$, and focus on the Grassmannian statistical model which is of common use in signal processing and statistics. Canonical Grassmannian distributions $\mathbb{G}_{\Sigma}$ on $G(m,r)$ are indexed by parameters $\Sigma$ from the manifold $\mathcal{M}= Pos_{sym}^{1}(m)$ of positive definite symmetric matrices of determinant $1$. Robust M-estimates of scatter (GE) for general probability measures $\mathcal{P}$ on $G(m,r)$ are studied. Such estimators are defined to be the maximizers of the Grassmannian log-likelihood $-\ell_{\mathcal{P}}(\Sigma)$ as function of $\Sigma$. One of the novel features of this work is a strong use of the fact that $\mathcal{M}$ is a CAT(0) space with known visual boundary at infinity $\partial \mathcal{M}$. We also recall that the sample space $G(m,r)$ is a part of $\partial \mathcal{M}$, show the distributions $\mathbb{G}_{\Sigma}$ are $SL(m,\mathbb{R})$--quasi-invariant, and that $\ell_{\mathcal{P}}(\Sigma)$ is a weighted Busemann function. Let $\mathcal{P}_n =(\delta_{U_1}+\cdots+\delta_{U_n})/n$ be the empirical probability measure for $n$-samples of random i.i.d. subspaces $U_i\in G(m,r)$ of common distribution $\mathcal{P}$, whose support spans $\mathbb{R}^m$. For $\Sigma_n$ and $\Sigma_{\mathcal{P}}$ the GEs of $\mathcal{P}_n$ and $\mathcal{P}$, we show the almost sure convergence of $\Sigma_n$ towards $\Sigma$ as $n\to\infty$ using methods from geometry, and provide a central limit theorem for the rescaled process $C_n = \frac{m}{tr(\Sigma_{\mathcal{P}}^{-1} \Sigma_n)}g^{-1} \Sigma_n g^{-1}$, where $\Sigma =gg$ with $g\in SL(m,\mathbb{R})$ the unique symmetric positive-definite square root of $\Sigma$.


title: Solar Potential Analysis of Rooftops Using Satellite Imagery

abstract: Solar energy is one of the most important sources of renewable energy and the cleanest form of energy. In India, where solar energy could produce power around trillion kilowatt-hours in a year, our country is only able to produce power of around in gigawatts only. Many people are not aware of the solar potential of their rooftop, and hence they always think that installing solar panels is very much expensive. In this work, we introduce an approach through which we can generate a report remotely that provides the amount of solar potential of a building using only its latitude and longitude. We further evaluated various types of rooftops to make our solution more robust. We also provide an approximate area of rooftop that can be used for solar panels placement and a visual analysis of how solar panels can be placed to maximize the output of solar power at a location.


title: Ellipsoids are the only local maximizers of the volume product

abstract: Using previous results about shadow systems and Steiner symmetrization, we prove that the local maximizers of the volume product of convex bodies are actually the global maximizers, that is: ellipsoids.


title: Relaxation of structures resulting from head-on mergers of ultralight bosonic dark matter cores

abstract: In this work we study some features of head-on mergers of equilibrium solutions of the Gross-Pitaevskii-Poisson system that rules the dynamics of the ultralight bosonic dark matter model. The importance of equilibrium solutions is that they play the role of halo cores in structure formation simulations. We consider a given range of initial conditions in order to sample the parameter space in terms of mass ratio and head-on momentum. In each case we analyze the relaxation process induced by gravitational cooling in the high and low momentum regimes and estimate the relaxation time scales in each case. We detect a low frequency mode in the whole parameters space and it was found that the resulting configuration oscillates under this mode with amplitude that depends on the mass ratio and head-on momentum. In some cases the resulting configuration oscillates with changes in density of two orders of magnitude and with a matter distribution that is far from isotropic. These results could contribute to the collection of possible mass distributions considered in the reconstruction of mass profiles obtained in structure formation simulations.


title: Spin chains for two-qubit teleportation

abstract: Generating high-quality multi-particle entanglement between communicating parties is the primary resource in quantum teleportation protocols. To this aim, we show that the natural dynamics of a single spin chain is able to sustain the generation of two pairs of Bell states - possibly shared between a sender and a distant receiver - which can in turn enable two-qubit teleportation. In particular, we address a spin-1/2 chain with XX interactions, connecting two pairs of spins located at its boundaries, playing the roles of sender and receiver. In the regime where both end pairs are weakly coupled to the spin chain, it is possible to generate at predefinite times a state that has vanishing infidelity with the product state of two Bell pairs, thereby providing nearly unit fidelity of teleportation. We also derive an effective Hamiltonian via a second-order perturbation approach that faithfully reproduces the dynamics of the full system.


title: State-of-the-Art Economic Load Dispatch of Power Systems Using Particle Swarm Optimization

abstract: Metaheuristic particle swarm optimization (PSO) algorithm has emerged as one of the most promising optimization techniques in solving highly constrained non-linear and non-convex optimization problems in different areas of electrical engineering. Economic operation of the power system is one of the most important areas of electrical engineering where PSO has been used efficiently in solving various issues of practical systems. In this paper, a comprehensive survey of research works in solving various aspects of economic load dispatch (ELD) problems of power system engineering using different types of PSO algorithms is presented. Five important areas of ELD problems have been identified, and the papers published in the general area of ELD using PSO have been classified into these five sections. These five areas are (i) single objective economic load dispatch, (ii) dynamic economic load dispatch, (iii) economic load dispatch with non-conventional sources, (iv) multi-objective environmental/economic dispatch, and (v) economic load dispatch of microgrids. At the end of each category, a table is provided which describes the main features of the papers in brief. The promising future works are given at the conclusion of the review.


title: Large dynamic range autorefraction with a low-cost diffuser wavefront sensor

abstract: Wavefront sensing with a thin diffuser has emerged as a potential low-cost alternative to a lenslet array for aberrometry. Diffuser wavefront sensors (DWS) have previously relied on tracking speckle displacement and consequently require coherent illumination. Here we show that displacement of caustic patterns can be tracked for estimating wavefront gradient, enabling the use of incoherent light sources and large dynamic-range wavefront measurements. We compare the precision of a DWS to a Shack-Hartmann wavefront sensor (SHWS) when using coherent, partially coherent, and incoherent illumination, in the application of autorefraction. We induce spherical and cylindrical errors in a model eye and use a multi-level Demon's non-rigid registration algorithm to estimate caustic displacements relative to an emmetropic model eye. When compared to spherical error measurements with the SHWS using partially coherent illumination, the DWS demonstrates a $\sim$5-fold improvement in dynamic range (-4.0 to +4.5 D vs. -22.0 to +19.5 D) with less than half the reduction in resolution (0.072 vs. 0.116 D), enabling a $\sim$3-fold increase in the number of resolvable prescriptions (118 vs. 358). In addition to being 40x lower-cost, the unique, non-periodic nature of the caustic pattern formed by a diffuser enables a larger dynamic range of aberration measurements compared to a lenslet array.


title: Oscillation modes of ultralight BEC dark matter cores

abstract: Structure formation simulations of ultralight bosonic dark matter, ruled by the Gross-Pitaevskii-Poisson (GPP) system of equations, show that dark matter clumps to form structures with density profiles consisting of a core surrounded by a power-law distribution. The core has a density profile similar to that of spherically symmetric equilibrium configurations of the GPP system. These configurations have been shown to be stable under a variety of perturbations and to have attractor properties. It is interesting to know the dominant frequencies of oscillation of these configurations. One reason is that in galaxies the effects of perturbations can trigger observable effects for specific frequency modes. Another reason is that during the process of structure formation, the oscillation modes can help to characterize a core. Based on these motivations, in this manuscript we present a systematic numerical analysis of the reaction of equilibrium configurations to various axi-symmetric perturbations with modes $l=0,1,2,3,4$. We then calculate the first few oscillation frequencies of equilibrium configurations for each mode.


title: The Device War - The War Between IOT Brands In A Household

abstract: Users buy compatible IOT devices from different brands with an expectation that their cooperation is smooth, but while function may superficially look friendly, cohabitation can subversively cause early battery depletion in competitor devices. The Wi-Fi Direct standard was introduced with the intention of simplifying peer-to-peer connections in home applications while helping devices to save power through centralization of effort into a single group owner device negotiated on start-up. Attacks on the group formation stage can be based on manipulating a victim device to frequently end up being assigned the group owner function, thereby depleting its batteries at faster rates than its peer devices. This manipulation is made easy by the group formation mechanism adopted by the standard. We show that group formation procedures could be better secured with features ensuring fairness by relying on commitments and by learning about the behavior observed for peer devices in the past. Simulations are used to quantify the resistance achieved against several attack strategies.


title: Work exchange, geometric magnetization and fluctuation-dissipation relations in a quantum dot under adiabatic magnetoelectric driving

abstract: We study the adiabatic dynamics of the charge, spin and energy of a quantum dot with a Coulomb interaction under two-parameter driving, associated to time-dependent gate voltage and magnetic field. The quantum dot is coupled to a single reservoir at temperature $T=0$, hence, the net energy dynamics is fully dissipative. However, in the presence of many-body interactions, other interesting mechanisms take place, like the net exchange of work between the two types of forces and the non-equilibrium accumulation of charge with different spin orientations. The latter has a geometric nature. The dissipation takes place in the form of an instantaneous Joule law with the universal resistance $R_0=h/2e^2$. We show the relation between this Joule law and instantaneous fluctuation-dissipation relations. The latter lead to generalized Korringa-Shiba relations, valid in the Fermi-liquid regime.


title: Metaplex networks: influence of the exo-endo structure of complex systems on diffusion

abstract: In a complex system the interplay between the internal structure of its entities and their interconnection may play a fundamental role in the global functioning of the system. Here, we define the concept of metaplex, which describes such trade-off between internal structure of entities and their interconnections. We then define a dynamical system on a metaplex and study diffusive processes on them. We provide analytical and computational evidences about the role played by the size of the nodes, the location of the internal coupling areas, and the strength and range of the coupling between the nodes on the global dynamics of metaplexes. Finally, we extend our analysis to two real-world metaplexes: a landscape and a brain metaplex. We corroborate that the internal structure of the nodes in a metaplex may dominate the global dynamics (brain metaplex) or play a regulatory role (landscape metaplex) to the influence of the interconnection between nodes.


title: Sharp local smoothing estimates for Fourier integral operators

abstract: The theory of Fourier integral operators is surveyed, with an emphasis on local smoothing estimates and their applications. After reviewing the classical background, we describe some recent work of the authors which established sharp local smoothing estimates for a natural class of Fourier integral operators. We also show how local smoothing estimates imply oscillatory integral estimates and obtain a maximal variant of an oscillatory integral estimate of Stein. Together with an oscillatory integral counterexample of Bourgain, this shows that our local smoothing estimates are sharp in odd spatial dimensions. Motivated by related counterexamples, we formulate local smoothing conjectures which take into account natural geometric assumptions arising from the structure of the Fourier integrals.


title: On the third-order Horadam and geometric mean sequences

abstract: This paper, in considering aspects of the geometric mean sequence, offers new results connecting generalized Tribonacci and third-order Horadam numbers which are established and then proved independently.


title: Generalized Hirano inverses in Banach Algebra

abstract: Let $A$ be a Banach Algebra, we say that $a\in A$ has generalized Hirano inverse if there exists some $b \in A$ such that $b=b.a.b, a.b=b.a$ and $a-a^2.b$ is quasi nil potent, if and only if there exists some $p=p^3\in A$ such that $a-p$ is quasi nil potent.


title: Using Reinforcement Learning to find Efficient Qubit Routing Policies for Deployment in Near-term Quantum Computers

abstract: This paper addresses the problem of qubit routing in first-generation and other near-term quantum computers. In particular, it is asserted that the qubit routing problem can be formulated as a reinforcement learning (RL) problem, and that this is sufficient, in principle, to discover the optimal qubit routing policy for any given quantum computer architecture. In order to achieve this, it is necessary to alter the conventional RL framework to allow combinatorial action space, and this represents a second contribution of this paper, which is expected to find additional application, beyond the qubit routing problem addressed herein. Numerical results are included demonstrating the advantage of the RL-trained qubit routing policy over using a sorting network.


title: Nonlinear coherent states associated with a measure on the positive real half line

abstract: We construct a class of generalized nonlinear coherent states by means of a newly obtained class of 2D complex orthogonal polynomials. The associated coherent states transform is discussed. A polynomials realization of the basis of the quantum states Hilbert space is also obtained. Here, the entire structure owes its existence to a certain measure on the positive real half line, of finite total mass, together with all its moments. We illustrate this construction with the example of the measure $r^\beta e^{-r}dr$, which leads to a new generalization of the true-polyanalytic Bargmann transform.


title: Measurement and modeling of electron-cloud-induced betatron tune shifts at the Cornell Electron-positron Storage Ring test accelerator

abstract: We report on extensive measurements at the Cornell Electron-positron Storage Ring of electron-cloud-induced betatron tune shifts for trains of positron bunches at 2.1 and 5.3 GeV with bunch populations ranging between 0.64x10^10 and 9.6x10^10. Measurements using a witness bunch with variable distance from the end of the train and variable bunch population provide information on cloud decay and cloud pinching during the bunch passage. We employ Monte Carlo simulations of the reflection and absorption of synchrotron radiation photons to determine the pattern of absorption sites around the circumference of the storage ring. The Geant4 simulation toolkit is used to model the interactions of the photons with the beampipe wall and determine the production energy and location distributions of the photoelectrons which seed the electron cloud. An electron cloud buildup model based on fitted ring-averaged secondary-yield properties of the vacuum chamber predicts tune shifts in good agreement with the measurements.


title: Ease on Down the Code: Complex Collaborative Qualitative Coding Simplified with 'Code Wizard'

abstract: This paper describes the design and development of a preliminary qualitative coding tool as well as a method to improve the process of achieving inter-coder reliability (ICR) in small teams. Software applications that support qualitative coding do not sufficiently assist collaboration among coders and overlook some fundamental issues related to ICR. We propose a new dimension of collaborative coding called "coders' certainty" and demonstrate its ability to illustrate valuable code disagreements that are missing from existing approaches. Through a case study, we describe the utility of our tool, Code Wizard, and how it helped a group of researchers effectively collaborate to code naturalistic observation data. We report the valuable lessons we learned from the development of our tool and method: (1) identifying coders' certainty constitutes an important part of determining the quality of data analysis and facilitates identifying overlapping and ambiguous codes, (2) making the details of coding process visible helps streamline the coding process and leads to a sense of ownership of the research results, and (3) there is valuable information hidden in coding disagreements that can be leveraged for improving the process of data analysis.


title: A Model for Elastic Evolution on Foliated Shapes

abstract: We study a shape evolution framework in which the deformation of shapes from time t to t + dt is governed by a regularized anisotropic elasticity model. More precisely, we assume that at each time shapes are infinitesimally deformed from a stress-free state to an elastic equilibrium as a result of the application of a small force. The configuration of equilibrium then becomes the new resting state for subsequent evolution. The primary motivation of this work is the modeling of slow changes in biological shapes like atrophy, where a body force applied to the volume represents the location and impact of the disease. Our model uses an optimal control viewpoint with the time derivative of force interpreted as a control, deforming a shape gradually from its observed initial state to an observed final state. Furthermore, inspired by the layered organization of cortical volumes, we consider a special case of our model in which shapes can be decomposed into a family of layers (forming a "foliation"). Preliminary experiments on synthetic layered shapes in two and three dimensions are presented to demonstrate the effect of elasticity.


title: Homogenization of Non-symmetric Jump Processes

abstract: We study the homogenization for a class of non-symmetric pure jump Feller processes. The jump intensity involves periodic and aperiodic constituents, as well as oscillating and non-oscillating constituents. This means that the noise can come both from the underlying periodic medium and from external environments, and is allowed to have different scales. It turns out that the Feller process converges in distribution, as the scaling parameter goes to zero, to a L\'evy process. As special cases of our result, some homogenization problems studied in previous works can be recovered. We also generalize the approach to the homogenization of symmetric stable-like processes with variable order. Moreover, we present some numerical experiments to demonstrate the usage of our homogenization results in the numerical approximation of first exit times.


title: On the Constraint Structure of Vacuum Energy Sequestering

abstract: We carry out the Hamiltonian analysis of the local vacuum energy sequestering model - a manifestly local and diffeomorphism invariant extension of general relativity which has been shown to remove the radiatively unstable contribution to the vacuum energy generated by matter loops. We find that the degravitation of this UV sensitive quantity is enforced via global relations that are a consequence of the model's peculiar constraint structure. We also show that the model propagates the proper number of degrees of freedom and thus locally reduces to general relativity on-shell.


title: Unfolding quantum master equation into a system of real-valued equations: computationally effective expansion over the basis of $SU(N)$ generators

abstract: Dynamics of an open $N$-state quantum system is typically modeled with a Markovian master equation describing the evolution of the system's density operator. By using generators of $SU(N)$ group as a basis, the density operator can be transformed into a real-valued 'Bloch vector'. The Lindbladian, a super-operator which serves a generator of the evolution, %in the master equation, can be expanded over the same basis and recast in the form of a real matrix. Together, these expansions result is a non-homogeneous system of $N^2-1$ real-valued linear differential equations for the Bloch vector. Now one can, e.g., implement a high-performance parallel simplex algorithm to find a solution of this system which guarantees exact preservation of the norm and Hermiticity of the density matrix. However, when performed in a straightforward way, the expansion turns to be an operation of the time complexity $\mathcal{O}(N^{10})$. The complexity can be reduced when the number of dissipative operators is independent of $N$, which is often the case for physically meaningful models. Here we present an algorithm to transform quantum master equation into a system of real-valued differential equations and propagate it forward in time. By using a scalable model, we evaluate computational efficiency of the algorithm and demonstrate that it is possible to handle the model system with $N = 10^3$ states on a single node of a computer cluster.


title: Figure 1 Theory Meets Figure 2 Experiments in the Study of Gene Expression

abstract: It is tempting to believe that we now own the genome. The ability to read and re-write it at will has ushered in a stunning period in the history of science. Nonetheless, there is an Achilles heel exposed by all of the genomic data that has accrued: we still don't know how to interpret it. Many genes are subject to sophisticated programs of transcriptional regulation, mediated by DNA sequences that harbor binding sites for transcription factors which can up- or down-regulate gene expression depending upon environmental conditions. This gives rise to an input-output function describing how the level of expression depends upon the parameters of the regulated gene { for instance, on the number and type of binding sites in its regulatory sequence. In recent years, the ability to make precision measurements of expression, coupled with the ability to make increasingly sophisticated theoretical predictions, have enabled an explicit dialogue between theory and experiment that holds the promise of covering this genomic Achilles heel. The goal is to reach a predictive understanding of transcriptional regulation that makes it possible to calculate gene expression levels from DNA regulatory sequence. This review focuses on the canonical simple repression motif to ask how well the models that have been used to characterize it actually work. We consider a hierarchy of increasingly sophisticated experiments in which the minimal parameter set learned at one level is applied to make quantitative predictions at the next. We show that these careful quantitative dissections provide a template for a predictive understanding of the many more complex regulatory arrangements found across all domains of life.


title: ${\rm SL}_2$ quantum trace in quantum Teichm\"uller theory via writhe

abstract: Quantization of the Teichm\"uller space of a punctured Riemann surface $S$ is an approach to $3$-dimensional quantum gravity, and is a prototypical example of quantization of cluster varieties. Any simple loop $\gamma$ in $S$ gives rise to a natural trace-of-monodromy function $\mathbb{I}(\gamma)$ on the Teichm\"uller space. For any ideal triangulation $\Delta$ of $S$, this function $\mathbb{I}(\gamma)$ is a Laurent polynomial in the square-roots of the exponentiated shear coordinates for the arcs of $\Delta$. An important problem was to construct a quantization of this function $\mathbb{I}(\gamma)$, namely to replace it by a noncommutative Laurent polynomial in the quantum variables. This problem, which is closely related to the framed protected spin characters in physics, has been solved by Allegretti and Kim using Bonahon and Wong's ${\rm SL}_2$ quantum trace for skein algebras, and by Gabella using Gaiotto, Moore and Neitzke's Seiberg-Witten curves, spectral networks, and writhe of links. We show that these two solutions to the quantization problem coincide. We enhance Gabella's solution and show that it is a twist of the Bonahon-Wong quantum trace.


title: Neutrino, parity violaton, V-A: a historical survey

abstract: A concise story of the rise of the four fermion theory of the universal weak interaction and its experimental confirmation, with a special emphasis on the problems related to parity violation.


title: A Comparison of Nonsmooth, Nonconvex, Constrained Optimization Solvers for the Design of Time-Delay Compensators

abstract: We present a detailed set of performance comparisons of two state-of-the-art solvers for the application of designing time-delay compensators, an important problem in the field of robust control. Formulating such robust control mechanics as constrained optimization problems often involves objective and constraint functions that are both nonconvex and nonsmooth, both of which present significant challenges to many solvers and their end-users hoping to obtain good solutions to these problems. In our particular engineering task, the main difficulty in the optimization arises in a nonsmooth and nonconvex stability constraint, which states that the infinite spectrum of zeros of the so-called shaper should remain in the open left half-plane. To perform our evaluation, we make use $\beta$-relative minimization profiles, recently introduced visualization tools that are particularly suited for benchmarking solvers on nonsmooth, nonconvex, constrained optimization problems. Furthermore, we also introduce new visualization tools, called Global-Local Profiles, which for a given problem and a fixed computational budget, assess the tradeoffs of distributing the budget over few or many starting points, with the former getting more budget per point and latter less.


title: Actor Conditioned Attention Maps for Video Action Detection

abstract: While observing complex events with multiple actors, humans do not assess each actor separately, but infer from the context. The surrounding context provides essential information for understanding actions. To this end, we propose to replace region of interest(RoI) pooling with an attention module, which ranks each spatio-temporal region's relevance to a detected actor instead of cropping. We refer to these as Actor-Conditioned Attention Maps (ACAM), which amplify/dampen the features extracted from the entire scene. The resulting actor-conditioned features focus the model on regions that are relevant to the conditioned actor. For actor localization, we leverage pre-trained object detectors, which transfer better. The proposed model is efficient and our action detection pipeline achieves near real-time performance. Experimental results on AVA 2.1 and JHMDB demonstrate the effectiveness of attention maps, with improvements of 7 mAP on AVA and 4 mAP on JHMDB.


title: Calculating Bardeen-Cooper-Schrieffer and Magnetic Superstructure Electronic States with $\Theta\Phi$

abstract: We propose the $\Theta\Phi$ (Theta-Phi) package which addresses two of the most important extensions of the essentially single-particle mean-field paradigm of the computational solid state physics: the admission of the Bardeen-Cooper-Schrieffer electronic ground state and allowance of the magnetically ordered states with an arbitrary superstructure (pitch) wave vector. Both features are implemented in the context of multi-band systems which paves the way to an interplay with the solid state quantum physics packages eventually providing access to the first-principles estimates of the relevant matrix elements of the model Hamiltonians derived from the standard DFT calculations. Several examples showing the workability of the proposed code are given.


title: Vector Boson production with heavy flavor quarks with the CMS experiment at LHC

abstract: In this paper, a brief outlook on the measurements of the total and differential cross-sections of vector bosons production associated with heavy flavor quarks performed using CMS proton-proton collision data taken at center of mass energies of 7, 8, and 13 TeV is presented. Obtained results are compared with predictions of various state-of-the-art Monte-Carlo generators.


title: Adaptation in multivariate log-concave density estimation

abstract: We study the adaptation properties of the multivariate log-concave maximum likelihood estimator over three subclasses of log-concave densities. The first consists of densities with polyhedral support whose logarithms are piecewise affine. The complexity of such densities~$f$ can be measured in terms of the sum $\Gamma(f)$ of the numbers of facets of the subdomains in the polyhedral subdivision of the support induced by $f$. Given $n$ independent observations from a $d$-dimensional log-concave density with $d \in \{2,3\}$, we prove a sharp oracle inequality, which in particular implies that the Kullback--Leibler risk of the log-concave maximum likelihood estimator for such densities is bounded above by $\Gamma(f)/n$, up to a polylogarithmic factor. Thus, the rate can be essentially parametric, even in this multivariate setting. For the second type of adaptation, we consider densities that are bounded away from zero on a polytopal support; we show that up to polylogarithmic factors, the log-concave maximum likelihood estimator attains the rate $n^{-4/7}$ when $d=3$, which is faster than the worst-case rate of $n^{-1/2}$. Finally, our third type of subclass consists of densities whose contours are well-separated; these new classes are constructed to be affine invariant and turn out to contain a wide variety of densities, including those that satisfy H\"older regularity conditions. Here, we prove another sharp oracle inequality, which reveals in particular that the log-concave maximum likelihood estimator attains a risk bound of order $n^{-\min\bigl(\frac{\beta+3}{\beta+7},\frac{4}{7}\bigr)}$ when $d=3$ over the class of $\beta$-H\"older log-concave densities with $\beta\in (1,3]$, again up to a polylogarithmic factor.


title: An explicit Waldspurger formula for Hilbert modular forms II

abstract: We describe a construction of preimages for the Shimura map on Hilbert modular forms using generalized theta series, and give an explicit Waldspurger type formula relating their Fourier coefficients to central values of twisted $L$-functions. This formula extends our previous work, allowing to compute these central values when the main central value vanishes.


title: On the Outage Performance of SWIPT Based Three-step Two-way DF Relay Networks

abstract: In this paper, we study the outage performance of simultaneous wireless information and power transfer (SWIP- T) based three-step two-way decode-and-forward (DF) relay networks, where both power-splitting (PS) and harvest-then-forward are employed. In particular, we derive the expressions of terminal-to-terminal (T2T) and system outage probabilities based on a Gaussian-Chebyshev quadrature approximation, and obtain the T2T and system outage capacities. The effects of various system parameters, e.g., the static power allocation ratio at the relay, symmetric PS, as well as asymmetric PS, on the outage performance of the investigated network are examined. It is shown that our derived expression for T2T outage capacity is more accurate than existing analytical results, and that the asymmetric PS achieves a higher system outage capacity than the symmetric one when the channels between the relay node and the terminal nodes have different statistic gains.


title: Extended \'etale homotopy groups from profinite Galois categories

abstract: In this note we show that the protruncated shape of a spectral $\infty$-topos is a delocalization of its profinite stratified shape. This gives a way to reconstruct the extended \'etale homotopy groups (i.e., the non-profinitely complete \'etale homotopy groups) of a coherent scheme from its profinite Galois category.


title: Thermalization in Solid-State NMR Controlled by Quantum Chaos in Spin Bath

abstract: We theoretically investigate thermalization and spin diffusion driven by a quantum spin bath for a realistic solid-state NMR experiment. We consider polycrystalline L-alanine, and investigate how the spin polarization spreads among several $^{13}$C nuclear spins, which interact via dipole-dipole coupling with the bath of strongly dipolar-coupled $^1$H nuclear (proton) spins. We do this by using direct numerical simulation of the many-spin time-dependent Schr\"odinger equation. We find that, although the proton spins located near the carbon sites interact most strongly with the $^{13}$C spins, this interaction alone is not enough to drive spin diffusion and thermalize the $^{13}$C nuclear spins. We demonstrate that the thermalization within the $^{13}$C subsystem is driven by the collective many-body dynamics of the proton spin bath, and specifically, that the onset of thermalization among the $^{13}$C spins is directly related to the onset of chaotic behavior in the proton spin bath. Therefore, thermalization and spin diffusion within the $^{13}$C subsystem is controlled by the proton spins located far from the C sites. In spite of their weak coupling to the $^{13}$C spins, these far-away protons help produce a network of strongly coupled proton spins with collective dynamics, that drives thermalization.


title: Improving the Kinetic Couplings in Lattice Non-Relativistic QCD

abstract: We improve the non-relativistic QCD (NRQCD) action by comparing the dispersion relation to that of the continuum through $\mathcal{O}(p^6)$ in perturbation theory. The one-loop matching coefficients of the $\mathcal{O}(p^4)$ kinetic operators are determined, as well as the scale at which to evaluate $\alpha_s$ in the $V$-scheme for each quantity. We utilise automated lattice perturbation theory using twisted boundary conditions as an infrared regulator. The one-loop radiative corrections to the mass renormalisation, zero-point energy and overall energy-shift of an NRQCD $b$-quark are also found. We also explore how a Fat$3$-smeared NRQCD action and changes of the stability parameter $n$ affect the coefficients. Finally, we use gluon field ensembles at multiple lattice spacing values, all of which include $u$, $d$, $s$ and $c$ quark vacuum polarisation, to test how the improvements affect the non-perturbatively determined $\Upsilon(1S)$ and $\eta_b(1S)$ kinetic masses, and the tuning of the $b$ quark mass.


title: Factors and connected factors in tough graphs with high isolated toughness

abstract: Let $G$ be a graph and let $f$ be a positive integer-valued function on $V(G)$. Assume that for all $S\subseteq V(G)$, $$\sum_{v\in I(G\setminus S)}f(v)(f(v)+1)\le |S|,$$ where $I(G\setminus S)$ denotes the set of isolated vertices of $G\setminus S$. In this paper, we show that if for all $S\subseteq V(G)$, $$\omega(G\setminus S)\le \sum_{v\in S}(f(v)-1)+1,$$ and $\sum_{v\in V(G)}f(v)$ is even, then $G$ has a factor $F$ such that for each vertex $v$, $d_F(v)=f(v)$, where $\omega(G\setminus S)$ denotes the number of components of $G\setminus S$. Moreover, we show that if for all $S\subseteq V(G)$, $$\omega(G\setminus S)\le \frac{1}{4}|S|+1,$$ and $f\ge 2$, then $G$ has a connected factor $H$ such that for each vertex $v$, $d_H(v)\in \{f(v),f(v)+1\}$.


title: What does the information paradox say about the universe?

abstract: The black hole information paradox is resolved in string theory by a radical change in the picture of the hole: black hole microstates are horizon sized quantum gravity objects called `fuzzballs' instead of vacuum regions with a central singularity. The requirement of causality implies that the quantum gravity wavefunctional $\Psi$ has an important component not present in the semiclassical picture: virtual fuzzballs. The large mass $M$ of the fuzzballs would suppress their virtual fluctuations, but this suppression is compensated by the large number -- $Exp[S_{bek}(M)]$ -- of possible fuzzballs. These fuzzballs are extended compression-resistant objects. The presence of these objects in the vacuum wavefunctional alters the physics of collapse when a horizon is about to form; this resolves the information paradox. We argue that these virtual fuzzballs also resist the curving of spacetime, and so cancel out the large cosmological constant created by the vacuum energy of local quantum fields. Assuming that the Birkoff theorem holds to leading order, we can map the black hole information problem to a problem in cosmology. Using the virtual fuzzball component of the wavefunctional, we give a qualitative picture of the evolution of $\Psi$ which is consistent with the requirements placed by the information paradox.


title: Machine Learning Prediction of DNA Charge Transport

abstract: First principle calculations of charge transfer in DNA molecules are computationally expensive given that charge carriers migrate in interaction with intra- and inter-molecular atomic motion. Screening sequences, e.g. to identify excellent electrical conductors is challenging even when adopting coarse-grained models and effective computational schemes that do not explicitly describe atomic dynamics. In this work, we present a machine learning (ML) model that allows the inexpensive prediction of the electrical conductance of millions of {\it long} double-stranded DNA (dsDNA) sequences, reducing computational costs by orders of magnitude. The algorithm is trained on {\it short} DNA nanojunctions with $n=3-7$ base pairs. The electrical conductance of the training set is computed with a quantum scattering method, which captures charge-nuclei scattering processes. We demonstrate that the ML method accurately predicts the electrical conductance of varied dsDNA junctions tracing different transport mechanisms: coherent (short-range) quantum tunneling, on-resonance (ballistic) transport, and incoherent site-to-site hopping. Furthermore, the ML approach supports physical observations that clusters of nucleotides regulate DNA transport behavior. The input features tested in this work could be used in other ML studies of charge transport in complex polymers, in the search for promising electronic and thermoelectric materials.


title: Free boundary problems with nonlocal and local diffusions I: Global solution

abstract: We study a class of free boundary problems of ecological models with nonlocal and local diffusions, which are natural extensions of free boundary problems of reaction diffusion systems in there local diffusions are used to describe the population dispersal, with the free boundary representing the spreading front of the species. We prove that such kind of nonlocal and local diffusion problems has a unique global solution, and then show that a spreading-vanishing dichotomy holds. Moreover, criteria of spreading and vanishing, and long time behavior of solution when spreading happens are established for the classical Lotka-Volterra competition and prey-predator models. Compared with free boundary systems with local diffusions as well as with nonlocal diffusions, the present paper involves some new difficulties, which should be overcome by use of new techniques. This is part I of a two part series, where we prove the existence, uniqueness, regularity and estimates of global solution. The spreading-vanishing dichotomy, criteria of spreading and vanishing, and long-time behavior of solution when spreading happens will be studied in the separate part II.


title: MERA as a holographic strange correlator

abstract: The multi-scale entanglement renormalization ansatz (MERA) is a tensor network that can efficiently parameterize critical ground states on a 1D lattice, and also suggestively implement some aspects of the holographic correspondence of string theory on a lattice. Extending our recent work [S. Singh, Physical Review D 97, 026012 (2018); S. Singh, N. A. McMahon, and G. K. Brennen, Phys. Rev. D 97, 026013 (2018)], we show how the MERA representation of a 1D critical ground state---which has long range entanglement---can be viewed as a strange correlator: the overlap of a 2D state with short range entanglement and a 2D product state. Strange correlators were recently introduced to map 2D symmetry protected or topologically ordered quantum states to critical systems in one lower dimension. The 2D quantum state dual to the input 1D critical state is obtained by lifting the MERA, a procedure which introduces bulk quantum degrees of freedom by inserting intertwiner tensors on each bond of the MERA tensor network. We show how this dual 2D bulk state exhibits several features of holography, for example, appearance of horizon-like holographic screens and bulk gauging of global on-site symmetries at the boundary. We also derive a quantum corrected Ryu-Takayanagi formula relating boundary entanglement entropy to bulk geodesic lengths---as measured by bulk entropy---and numerically test it for ground states of a set of unitary minimal model CFTs, as realized by 1D anyonic Heisenberg models.