title: An efficient Monte Carlo interior penalty discontinuous Galerkin method for the time-harmonic Maxwell's equations with random coefficients

abstract: This paper develops an efficient Monte Carlo interior penalty discontinuous Galerkin method for electromagnetic wave propagation in random media. This method is based on a multi-modes expansion of the solution to the time-harmonic random Maxwell equations. It is shown that each mode function satisfies a Maxwell system with random sources defined recursively. An unconditionally stable IP-DG method is employed to discretize the nearly deterministic Maxwell system and the Monte Carlo method combined with an efficient acceleration strategy is proposed for computing the mode functions and the statistics of the electromagnetic wave. A complete error analysis is established for the proposed multi-modes Monte Carlo IP-DG method. It is proved that the proposed method converges with an optimal order for each of three levels of approximations. Numerical experiments are provided to validate the theoretical results and to gauge the performance of the proposed numerical method and approach.


title: Leading Coefficients and the Multiplicity of Known Roots

abstract: We show that a monic univariate polynomial over a field of characteristic zero, with $k$ distinct non-zero known roots, is determined by its $k$ proper leading coefficients by providing an explicit algorithm for computing the multiplicities of each root. We provide a version of the result and accompanying algorithm when the field is not algebraically closed by considering the minimal polynomials of the roots. Furthermore, we show how to perform the aforementioned algorithm in a numerically stable manner over $\mathbb{C}$, and then apply it to obtain new characteristic polynomials of hypergraphs.


title: Control of Magnetic Dipole Emission with Surface Plasmon Polaritons

abstract: The possibility to use surface plasmon polaritons for enhancement of weak magnetic dipole transitions is analyzed theoretically and demonstrated experimentally for simple flat geometry and sine-wave profile modulated plasmonic films. Spontaneous emission of Eu3+ in organic matrices deposited onto plasmonic surfaces demonstrates specific angular and polarization patterns at both electric and magnetic dipole transitions with a well-defined maximum at the plasmon decoupling conditions manifesting the character and magnitude of the effect.


title: Modeling the Ferroelectric Phase Transition in Perovskite Relaxors and New Corresponding Solving Method

abstract: Three-dimensional-random-site-Ising-model (3D-RSIM) along with Glauber dynamics of pseudo-spins (PSs) is applied for the first time to study the relaxion-ferroelectric-phase-transition (RFPT) of perovskite-relaxor-ferroelectrics (PRFEs) in detail. To solve this model, we proposed a new method, mean-field of PS-strings (MF-PSSs). We find: 1) 3D-RSIM is a mixed system consisting of permanent-paraelectric-subsystem (PPSS), low-transition-temperature ferroelectric-subsystem (LTT-FSS), and high-transition-temperature ferroelectric-subsystem (HTT-FSS). The contents of these three subsystems change with the PS concentration (phi) of the model; 2) When phi < phi_p (phi_p is the percolation threshold of the 3D-RSIM), HTT-FSS dominates the whole system. On cooling, the system undergoes an inhomogeneous-diffuse-ferroelectric-phase-transition (IDFPT), and as phi increases, it shows a crossover from the critical, gradually to Vogel-Fulcher, and to Arrehnius type slowdown. Specifically, we observe this crossover from the relationship between the peak temperature (Tm) of the real susceptibility of the correlated relaxation of PSs and angular frequency (omega); 3) When phi = phi_p, the three subsystems have comparable contents. However, as temperature decreases, the spontaneous polarization and specific heat of the whole system show an IDFPT, but Tm shifts to lower temperature with omega until 0K. Since the relaxation time of PPSS rapidly increases at lower temperature and diverges at 0K, during a cooling process at a limited rate, PPSS must be frozen at a certain temperature, i.e. PS glass transition occurs. In short, with the increase of phi, 3D-RSIGM gradually evolves from the IDFPT system to PS glass system, and phi_p is the characteristic concentration of this evolution. Moreover, based on the coupling between PSs and crystal lattice, we also give a new possible mechanism of Burns transformation.


title: Planar embeddings of chainable continua

abstract: We prove that for a chainable continuum $X$ and every non-zigzag $x\in X$ there exists a planar embedding $\phi:X\to \phi(X)\subset\mathbb R^2$ such that $\phi(x)$ is accessible, partially answering the question of Nadler and Quinn from 1972. Two embeddings $\phi,\psi:X \to \mathbb R^2$ are called strongly equivalent if $\phi \circ \psi^{-1}: \psi(X) \to \phi(X)$ can be extended to a homeomorphism of $\mathbb R^2$. We also prove that every indecomposable chainable continuum can be embedded in the plane in uncountably many strongly non-equivalent ways.


title: Human Activity Recognition Based on Wearable Sensor Data: A Standardization of the State-of-the-Art

abstract: Human activity recognition based on wearable sensor data has been an attractive research topic due to its application in areas such as healthcare and smart environments. In this context, many works have presented remarkable results using accelerometer, gyroscope and magnetometer data to represent the activities categories. However, current studies do not consider important issues that lead to skewed results, making it hard to assess the quality of sensor-based human activity recognition and preventing a direct comparison of previous works. These issues include the samples generation processes and the validation protocols used. We emphasize that in other research areas, such as image classification and object detection, these issues are already well-defined, which brings more efforts towards the application. Inspired by this, we conduct an extensive set of experiments that analyze different sample generation processes and validation protocols to indicate the vulnerable points in human activity recognition based on wearable sensor data. For this purpose, we implement and evaluate several top-performance methods, ranging from handcrafted-based approaches to convolutional neural networks. According to our study, most of the experimental evaluations that are currently employed are not adequate to perform the activity recognition in the context of wearable sensor data, in which the recognition accuracy drops considerably when compared to an appropriate evaluation approach. To the best of our knowledge, this is the first study that tackles essential issues that compromise the understanding of the performance in human activity recognition based on wearable sensor data.


title: Well-posedness and numerical treatment of the Blackstock equation in nonlinear acoustics

abstract: We study the Blackstock equation which models the propagation of nonlinear sound waves through dissipative fluids. Global well-posedness of the model with homogeneous Dirichlet boundary conditions is shown for small initial data. To this end, we employ a fixed-point technique coupled with well-posedness results for a linearized model and appropriate energy estimates. Furthermore, we obtain exponential decay for the energy of the solution. We present additionally a finite element-based method for solving the Blackstock equation and illustrate the behavior of solutions through several numerical experiments.


title: 3D-CODED : 3D Correspondences by Deep Deformation

abstract: We present a new deep learning approach for matching deformable shapes by introducing {\it Shape Deformation Networks} which jointly encode 3D shapes and correspondences. This is achieved by factoring the surface representation into (i) a template, that parameterizes the surface, and (ii) a learnt global feature vector that parameterizes the transformation of the template into the input surface. By predicting this feature for a new shape, we implicitly predict correspondences between this shape and the template. We show that these correspondences can be improved by an additional step which improves the shape feature by minimizing the Chamfer distance between the input and transformed template. We demonstrate that our simple approach improves on state-of-the-art results on the difficult FAUST-inter challenge, with an average correspondence error of 2.88cm. We show, on the TOSCA dataset, that our method is robust to many types of perturbations, and generalizes to non-human shapes. This robustness allows it to perform well on real unclean, meshes from the the SCAPE dataset.


title: Identifying Recurring Patterns with Deep Neural Networks for Natural Image Denoising

abstract: Image denoising methods must effectively model, implicitly or explicitly, the vast diversity of patterns and textures that occur in natural images. This is challenging, even for modern methods that leverage deep neural networks trained to regress to clean images from noisy inputs. One recourse is to rely on "internal" image statistics, by searching for similar patterns within the input image itself. In this work, we propose a new method for natural image denoising that trains a deep neural network to determine whether patches in a noisy image input share common underlying patterns. Given a pair of noisy patches, our network predicts whether different sub-band coefficients of the original noise-free patches are similar. The denoising algorithm then aggregates matched coefficients to obtain an initial estimate of the clean image. Finally, this estimate is provided as input, along with the original noisy image, to a standard regression-based denoising network. Experiments show that our method achieves state-of-the-art color image denoising performance, including with a blind version that trains a common model for a range of noise levels, and does not require knowledge of level of noise in an input image. Our approach also has a distinct advantage when training with limited amounts of training data.


title: Dependently Typed Folds for Nested Data Types

abstract: We present an approach to develop folds for nested data types using dependent types. We call such folds $\textit{dependently typed folds}$, they have the following properties. (1) Dependently typed folds are defined by well-founded recursion and they can be defined in a total dependently typed language. (2) Dependently typed folds do not depend on maps, map functions and many terminating functions can be defined using dependently typed folds. (3) The induction principles for nested data types follow from the definitions of dependently typed folds and the programs defined by dependently typed folds can be formally verified. (4) Dependently typed folds exist for any nested data types and they can be specialized to the traditional $\textit{higher-order folds}$. Using various of examples, we show how to program and reason about dependently typed folds. We also show how to obtain dependently typed folds in general and how to specialize them to the corresponding higher-order folds.


title: Beyond Bags of Words: Inferring Systemic Nets

abstract: Textual analytics based on representations of documents as bags of words have been reasonably successful. However, analysis that requires deeper insight into language, into author properties, or into the contexts in which documents were created requires a richer representation. Systemic nets are one such representation. They have not been extensively used because they required human effort to construct. We show that systemic nets can be algorithmically inferred from corpora, that the resulting nets are plausible, and that they can provide practical benefits for knowledge discovery problems. This opens up a new class of practical analysis techniques for textual analytics.


title: A latent spatial factor approach for synthesizing opioid associated deaths and treatment admissions in Ohio counties

abstract: Background: Opioid misuse is a major public health issue in the United States and in particular Ohio. However, the burden of the epidemic is challenging to quantify as public health surveillance measures capture different aspects of the problem. Here we synthesize county-level death and treatment counts to compare the relative burden across counties and assess associations with social environmental covariates. Methods: We construct a generalized spatial factor model to jointly model death and treatment rates for each county. For each outcome, we specify a spatial rates parameterization for a Poisson regression model with spatially varying factor loadings. We use a conditional autoregressive model to account for spatial dependence within a Bayesian framework. Results: The estimated spatial factor was highest in the southern and southwestern counties of the state, representing a higher burden of the opioid epidemic. We found that relatively high rates of treatment contributed to the factor in the southern part of the state; whereas, relatively higher rates of death contributed in the southwest. The estimated factor was also positively associated with the proportion of residents aged 18-64 on disability and negatively associated with the proportion of residents reporting white race. Conclusions: We synthesized the information in the opioid associated death and treatment counts through a spatial factor model to estimate a latent factor representing the consensus between the two surveillance measures. We believe this framework provides a coherent approach to describe the epidemic while leveraging information from multiple surveillance measures.


title: End-to-End Parkinson Disease Diagnosis using Brain MR-Images by 3D-CNN

abstract: In this work, we use a deep learning framework for simultaneous classification and regression of Parkinson disease diagnosis based on MR-Images and personal information (i.e. age, gender). We intend to facilitate and increase the confidence in Parkinson disease diagnosis through our deep learning framework.


title: Understanding the Meaning of Understanding

abstract: Can we train a machine to detect if another machine has understood a concept? In principle, this is possible by conducting tests on the subject of that concept. However we want this procedure to be done by avoiding direct questions. In other words, we would like to isolate the absolute meaning of an abstract idea by putting it into a class of equivalence, hence without adopting straight definitions or showing how this idea "works" in practice. We discuss the metaphysical implications hidden in the above question, with the aim of providing a plausible reference framework.


title: ${\rm SDiff}(S^2)$ and the orbit method

abstract: The group of area preserving diffeomorphisms of the two sphere, ${\rm SDiff}(S^2)$, is one of the simplest examples of an infinite dimensional Lie group. It plays a key role in incompressible hydrodynamics and it recently appeared in general relativity as a subgroup of two closely related, newly defined symmetry groups. We investigate its representation theory using the method of coadjoint orbits. We describe the Casimir functions and the Cartan algebra. Then we evaluate the trace of a simple ${\rm SDiff}(S^2)$ operator using the Atiyah-Bott fixed point formula. The trace is divergent but we show that it has well-defined truncations related to the structure of ${\rm SDiff}(S^2)$. Finally, we relate our results back to the recent appearances of ${\rm SDiff}(S^2)$ in black hole physics.


title: Manifold Mixup: Better Representations by Interpolating Hidden States

abstract: Deep neural networks excel at learning the training data, but often provide incorrect and confident predictions when evaluated on slightly different test examples. This includes distribution shifts, outliers, and adversarial examples. To address these issues, we propose Manifold Mixup, a simple regularizer that encourages neural networks to predict less confidently on interpolations of hidden representations. Manifold Mixup leverages semantic interpolations as additional training signal, obtaining neural networks with smoother decision boundaries at multiple levels of representation. As a result, neural networks trained with Manifold Mixup learn class-representations with fewer directions of variance. We prove theory on why this flattening happens under ideal conditions, validate it on practical situations, and connect it to previous works on information theory and generalization. In spite of incurring no significant computation and being implemented in a few lines of code, Manifold Mixup improves strong baselines in supervised learning, robustness to single-step adversarial attacks, and test log-likelihood.


title: Weighted Tanimoto Coefficient for 3D Molecule Structure Similarity Measurement

abstract: Similarity searching of molecular structure has been an important application in the Chemoinformatics, especially in drug discovery. Similarity searching is a common method used for identification of molecular structure. It involve three main principal component of similarity searching: structure representation; weighting scheme; and similarity coefficient. In this paper, we introduces Weighted Tanimoto Coefficient based on weighted Euclidean distance in order to investigate the effect of weight function on the result for similarity searching. The Tanimoto coefficient is one of the popular similarity coefficients used to measure the similarity between pairs of the molecule. The most of research area found that the similarity searching is based on binary or fingerprint data. Meanwhile, we used non-binary data and was set amphetamine structure as a reference or targeted structure and the rest of the dataset becomes a database structure. Throughout this study, it showed that there is definitely gives a different result between a similarity searching with and without weight.


title: The ORNL Analysis Technique for Extracting $\beta$-Delayed Multi-Neutron Branching Ratios with BRIKEN

abstract: Many choices are available in order to evaluate large radioactive decay networks. %multi-particle decay data. There are many parameters that influence the calculated $\beta$-decay delayed single and multi-neutron emission branching fractions. We describe assumptions about the decay model, background, and other parameters and their influence on $\beta$-decay delayed multi-neutron emission analysis. An analysis technique, the ORNL BRIKEN analysis procedure, for determining $\beta$-delayed multi-neutron branching ratios in $\beta$-neutron precursors produced by means of heavy-ion fragmentation is presented. The technique is based on estimating the initial activities of zero, one, and two neutrons occurring in coincidence with an ion-implant and $\beta$ trigger. The technique allows one to extract $\beta$-delayed multi-neutron decay branching ratios measured with the hybrid \textsuperscript{3}He BRIKEN neutron counter. As an example, two analyses of the $\beta$-neutron emitter \textsuperscript{77}Cu based on different {\it a priori} assumptions are presented along with comparisons to literature values.


title: The $e$-vector of a simplicial complex

abstract: We study the exponential Hilbert series (both coarsely- and finely-graded) of the Stanley-Reisner ring of an abstract simplicial complex, $\Delta$, and we introduce the $e$-vector of $\Delta$, which relates to the coefficients of the exponential Hilbert series. We explore the relationship of the $e$-vector with the classical $f$-vector and $h$-vector of $\Delta$ while simultaneously investigating the geometric information that the $e$-vector encodes about $\Delta$. We then prove a simple combinatorial identity for the $e$-vector in the case where $\Delta$ is an Eulerian manifold.


title: The mechanism of nanoparticle precipitation induced by electron irradiation in transmission electron microscopy

abstract: Ru nanoparticles (NPs) can precipitate in Ru doped SiO2 amorphous thin films, triggered by electron irradiation in (scanning) transmission electron microscope ((S)TEM). A new mechanism was introduced to interpret the formation of metal NPs in (S)TEM. The induced electric field by electron irradiation, which originates from charging due to ionizations and excitations of atom electrons, can reduce the Gibbs free energy barrier for nucleation of metal particles. Furthermore, the directional ion drifting driven by the electric forces may accelerate the kinetic process of metal particle precipitation.


title: Quasi-tight Framelets with Directionality or High Vanishing Moments Derived from Arbitrary Refinable Functions

abstract: Construction of multivariate tight framelets is known to be a challenging problem. Multivariate dual framelets with vanishing moments generalize tight framelets and are not easy to be constructed either. Compactly supported multivariate framelets with directionality or high vanishing moments are of interest and importance in both theory and applications. In this paper we introduce the notion of a quasi-tight framelet, which is a dual framelet, but behaves almost like a tight framelet. Let $\phi\in L_2(R^d)$ be an arbitrary compactly supported $M$-refinable function such that its underlying low-pass filter satisfies the basic sum rule. We first constructively prove by a step-by-step algorithm that we can always easily derive from the arbitrary $M$-refinable function $\phi$ a directional compactly supported quasi-tight $M$-framelet in $L_2(R^d)$ associated with a directional quasi-tight $M$-framelet filter bank, each of whose high-pass filters has only two nonzero coefficients with opposite signs. If in addition all the coefficients of its low-pass filter are nonnegative, such a quasi-tight $M$-framelet becomes a directional tight $M$-framelet in $L_2(R^d)$. Furthermore, we show by a constructive algorithm that we can always derive from the arbitrary $M$-refinable function $\phi$ a compactly supported quasi-tight $M$-framelet in $L_2(R^d)$ with the highest possible order of vanishing moments. We shall also present a result on quasi-tight framelets whose associated high-pass filters are purely differencing filters with the highest order of vanishing moments. Several examples will be provided to illustrate our main theoretical results and algorithms in this paper.


title: Blockchain Enabled Enhanced IoT Ecosystem Security

abstract: Blockchain (BC), the technology behind the Bitcoin cryptocurrency system, is starting to be adopted for ensuring enhanced security and privacy in the Internet of Things (IoT) ecosystem. Fervent research is currently being focused in both academia and industry in this domain. Proof of Work (PoW), a cryptographic puzzle, plays a vital role in ensuring BC security by maintaining a digital ledger of transactions, which are considered to be incorruptible. Furthermore, BC uses a changeable Public Key (PK) to record the identity of users, thus providing an extra layer of privacy. Not only in cryptocurrency has the successful adoption of the BC been implemented, but also in multifaceted non-monetary systems, such as in: distributed storage systems, proof of location and healthcare. Recent research articles and projects or applications were surveyed to assess the implementation of the BC for IoT Security and identify associated challenges and propose solutions for BC enabled enhanced security for the IoT ecosystem.


title: Spectral Test of the MIXMAX Random Number Generators

abstract: An important statistical test on the pseudo-random number generators is called the spectral test. The test is aimed at answering the question of distribution of the generated pseudo-random vectors in dimensions $d$ that are larger than the genuine dimension of a generator $N$. In particular, the default MIXMAX generators have various dimensions: $N=8,17,240$ and higher. Therefore the spectral test is important to perform in dimensions $d > 8$ for $N=8$ generator, $d> 17$ for $N=17$ and $d> 240$ for $N=240$ generator. These tests have been performed by L'Ecuyer and collaborators. When $d > N$ the vectors of the generated numbers fall into the parallel hyperplanes and the distances between them can be larger than the genuine "resolution" of the MIXMAX generators, which is $ l=2^{-61}$. The aim of this article is to further study the spectral properties of the MIXMAX generators, to investigate the dependence of the spectral properties of the MIXMAX generators as a function of their internal parameters and in particular their dependence on the parameter $m$. We found that the best spectral properties are realized when $m$ is between $2^{24}$ and $2^{36}$, a range which is inclusive of the value of the $N=17$ generator. We also provide the alternative parameters for the generators, $N=8$ and $N=240$ with $m$ in this optimised range.


title: On a scalar field of time and electromagnetism

abstract: It is possible to describe a universal scalar field of time but not a universal coordinate of time and to attribute its non-geodesic alignment to the electromagnetic phenomena. A very surprising outcome is that not only mass generates gravity, but also electric charge does. Charge is, however, coupled to a non-geodesic vector field and thus is not totally equivalent to inertial mass. Only the entire "Energy-Momentum" tensor has a vanishing divergence. The model can be seen as misalignment of physically accessible events in an observer spacetime and of gravity as a controlling response by volumetric contraction of the observer spacetime in the direction where events bend or accelerate to. This non geodesic acceleration is described by a generalization of the Reeb vector. Misalignment of events can be described by 1, 2, and 3 such vectors. The paper presents a term with 4 vectors but does not discuss its physical meaning. The paper also discusses particle mass ratios and the Fine Structure Constant where added or subtracted area in relation to a disk does not involve a ratio 1/24 but 1/96 due to the physical meaning of the orientation of a space foliation which is perpendicular to a time-like vector and due to the orientation of a plane which is perpendicular to a time-like vector. These two orientations mean that only one side of a 3-dimensional foliation has a physical meaning and only one side of a sub-plane of that foliation has a physical meaning then (1/2)*(1/2)*(1/24_=1/96. Another interpretation of the factor 1/4 is the Bekenstein - Hawking entropy to area constant. An additional coefficient 4/{\pi} describes an acceleration field strength and has a compelling source in mainstream physics. Other two field strength coefficients are 95/96 and a critical value due to an imbalance equation between gravity and anti-gravity ~1.556198537190348396563877031439915299.


title: Sinks and sources for C1 dynamics whose Lyapunov exponents have constant sign

abstract: Let $f:M\to M$ be a $C^1$ map of a compact manifold $M$, with dimension at least $2$, admitting some point whose future trajectory has only negative Lyapunov exponents. Then this trajectory converges to a periodic sink. We need only assume that $Df$ is never the null map at any point (in particular, we need no extra smoothness assumption on $Df$), encompassing a wide class of possible critical behavior. Similarly, a trajectory having only positive Lyapunov exponents for a $C^1$ diffeomorphism is itself a periodic repeller (source). Analogously for a $C^1$ open and dense subset of vector field on finite dimensional manifolds: for a flow $\phi_t$ generated by such a vector field, if a trajectory admits weak asymptotic sectional contraction (the extreme rates of expansion of the Linear Poincar\'e Flow are all negative), then this trajectory belongs either to the basin of attraction of a periodic hyperbolic attracting orbit (a periodic sink or an attracting equilibrium); or the trajectory accumulates a codimension one saddle singularity. Similar results hold for weak sectional expanding trajectories. Both results extend part of the non-uniform hyperbolic theory (Pesin's Theory) from the $C^{1+}$ diffeomorphism setting to $C^1$ endomorphisms and $C^1$ flows. Some ergodic theoretical consequences are discussed. The proofs use versions of Pliss' Lemma for maps and flows translated as (reverse) hyperbolic times, and a result ensuring that certain subadditive cocycles over vector fields are in fact additive.


title: Measurements of properties of the Higgs boson decaying to a W boson pair in pp collisions at $\sqrt{s} =$ 13 TeV

abstract: Measurements of the production of the standard model Higgs boson decaying to a W boson pair are reported. The W$^+$W$^-$ candidates are selected in events with an oppositely charged lepton pair, large missing transverse momentum, and various numbers of jets. To select Higgs bosons produced via vector boson fusion and associated production with a W or Z boson, events with two jets or three or four leptons are also selected. The event sample corresponds to an integrated luminosity of 35.9 fb$^{-1}$, collected in pp collisions at $\sqrt{s} =$ 13 TeV by the CMS detector at the LHC during 2016. Combining all channels, the observed cross section times branching fraction is 1.28 $^{+0.18}_{-0.17}$ times the standard model prediction for the Higgs boson with a mass of 125.09 GeV. This is the first observation of the Higgs boson decay to W boson pairs by the CMS experiment.


title: Reduced words for clans

abstract: Clans are combinatorial objects indexing the orbits of $GL(\mathbb{C}^p) \times GL(\mathbb{C}^q)$ on the variety of flags in $\mathbb{C}^{p+q}$. This geometry leads to a partial order on the set of clans analogous to weak Bruhat order on the symmetric group, and we study the saturated chains in this order. We prove an analogue of the Matsumoto-Tits theorem on reduced words in a Coxeter group. We also obtain enumerations of reduced word sets for particular clans in terms of standard tableaux and shifted standard tableaux.


title: Cover-Encodings of Fitness Landscapes

abstract: The traditional way of tackling discrete optimization problems is by using local search on suitably defined cost or fitness landscapes. Such approaches are however limited by the slowing down that occurs when the local minima that are a feature of the typically rugged landscapes encountered arrest the progress of the search process. Another way of tackling optimization problems is by the use of heuristic approximations to estimate a global cost minimum. Here we present a combination of these two approaches by using cover-encoding maps which map processes from a larger search space to subsets of the original search space. The key idea is to construct cover-encoding maps with the help of suitable heuristics that single out near-optimal solutions and result in landscapes on the larger search space that no longer exhibit trapping local minima. We present cover-encoding maps for the problems of the traveling salesman, number partitioning, maximum matching and maximum clique; the practical feasibility of our method is demonstrated by simulations of adaptive walks on the corresponding encoded landscapes which find the global minima for these problems.


title: Observational signature of near-extremal Kerr-like black holes in a modified gravity theory at the Event Horizon Telescope

abstract: We study the shadows cast by near-extremal Kerr-MOG black holes for different values of the parameter in modified gravity (MOG). In particular, we consider an isotropic emitter orbiting near such black holes and analytically compute the positions, fluxes and redshift factors of their images. The size of the shadow decreases when the modified parameter is increased. For each shadow, the images of the emitter appear on a special part of the shadow which has a rich structure. The primary image and secondary images are similar to those produced for the near-extremal (high spin) Kerr black hole, but the near-extremal Kerr-MOG black hole can have a spin ($\hat{J}/M^2_{\alpha}$) which is finitely lower than 1. When the modified parameter is varied, the typical positions of the corresponding images do not change, nor does the typical redshift factor associated with the primary image. However, another typical redshift factor associated with the secondary image increases when the modified parameter is increased. We also find that the fluxes increase in that case. These images appear periodically with period greater than that of Kerr. This provides an alternative signature away from the Kerr case which may be tested by the Event Horizon Telescope.


title: What About Applied Fairness?

abstract: Machine learning practitioners are often ambivalent about the ethical aspects of their products. We believe anything that gets us from that current state to one in which our systems are achieving some degree of fairness is an improvement that should be welcomed. This is true even when that progress does not get us 100% of the way to the goal of "complete" fairness or perfectly align with our personal belief on which measure of fairness is used. Some measure of fairness being built would still put us in a better position than the status quo. Impediments to getting fairness and ethical concerns applied in real applications, whether they are abstruse philosophical debates or technical overhead such as the introduction of ever more hyper-parameters, should be avoided. In this paper we further elaborate on our argument for this viewpoint and its importance.


title: A tight-coupling scheme sharing minimum information across a spatial interface between gyrokinetic turbulence codes

abstract: A new scheme that tightly couples kinetic turbulence codes across a spatial interface is introduced. This scheme evolves from considerations of competing strategies and down-selection. It is found that the use of a composite kinetic distribution function and fields with global boundary conditions as if the coupled code were one, makes the coupling problem tractable. In contrast, coupling the two solutions from each code across the overlap region is found to be more difficult due to numerical dephasing of the turbulent solutions between two solvers. Another advantage of the new scheme is that the data movement can be limited to the 3D fluid quantities, instead of higher dimensional kinetic information, which is computationally more efficient for large scale simulations on leadership class computers.


title: Finding your Lookalike: Measuring Face Similarity Rather than Face Identity

abstract: Face images are one of the main areas of focus for computer vision, receiving on a wide variety of tasks. Although face recognition is probably the most widely researched, many other tasks such as kinship detection, facial expression classification and facial aging have been examined. In this work we propose the new, subjective task of quantifying perceived face similarity between a pair of faces. That is, we predict the perceived similarity between facial images, given that they are not of the same person. Although this task is clearly correlated with face recognition, it is different and therefore justifies a separate investigation. Humans often remark that two persons look alike, even in cases where the persons are not actually confused with one another. In addition, because face similarity is different than traditional image similarity, there are challenges in data collection and labeling, and dealing with diverging subjective opinions between human labelers. We present evidence that finding facial look-alikes and recognizing faces are two distinct tasks. We propose a new dataset for facial similarity and introduce the Lookalike network, directed towards similar face classification, which outperforms the ad hoc usage of a face recognition network directed at the same task.


title: The weak Bernoulli property for matrix Gibbs states

abstract: We study the ergodic properties of a class of measures on $\Sigma^{\mathbb{Z}}$ for which $\mu_{\mathcal{A},t}[x_{0}\cdots x_{n-1}]\approx e^{-nP}\left \|A_{x_{0}}\cdots A_{x_{n-1}}\right \| ^{t}$, where $\mathcal{A}=(A_{0}, \ldots , A_{M-1})$ is a collection of matrices. The measure $\mu_{\mathcal{A},t}$ is called a matrix Gibbs state. In particular we give a sufficient condition for a matrix Gibbs state to have the weak Bernoulli property. We employ a number of techniques to understand these measures including a novel approach based on Perron-Frobenius theory. We find that when $t$ is an even integer the ergodic properties of $\mu_{\mathcal{A} ,t}$ are readily deduced from finite dimensional Perron-Frobenius theory. We then consider an extension of this method to $t>0$ using operators on an infinite dimensional space. Finally we use a general result of Bradley to prove the main theorem.


title: Properness for circle packings and Delaunay circle patterns on complex projective structures

abstract: We consider circle packings and, more generally, Delaunay circle patterns - arrangements of circles arising from a Delaunay decomposition of a finite set of points - on surfaces equipped with a complex projective structure. Motivated by a conjecture of Kojima, Mizushima and Tan, we prove that the forgetful map sending a complex projective structure admitting a circle packing with given nerve (resp. a Delaunay circle pattern with given nerve and intersection angles) to the underlying complex structure is proper.


title: Neutronic Analysis of a Nuclear-Chicago NH3 Neutron Howitzer

abstract: Neutron howitzers are the irradiation instruments, generally used for the chemical analysis in the research laboratories, where the neutron sources are encapsulated within a moderating medium and the target materials are irradiated at the channels inserted through the corresponding neutron sources. In this study, we employ Monte Carlo simulations by using the MCNP5 code in order to perform the neutronic analysis of a Nuclear-Chicago NH3 neutron howitzer with a 5-Ci 239Pu-Be neutron source. This neutron howitzer is a cylindrical drum that is covered by aluminum at the level of outermost layer, and the cylindrical body is installed with paraffin for the neutron moderation. The drum consists of two vertical and two horizontal irradiation channels and a vertical source channel in which the 5-Ci 239Pu-Be neutron source is located and moved upward or downward by the aid of a plexiglass rod. We focus our study on the two horizontal irradiation channels, separated by an angle of 120o, which are installed with a cylindrical plexiglass bar. We construct our geometry by respecting the dimensional properties of the NH3 howitzer components and we define a biased cylindrical source of 1.02 in.x 4.425 in. by taking into account the neutron energy spectrum and source strength of the 5-Ci 239Pu-Be neutron source. We determine the neutron flux within the plexiglass bar for various energy bins in the energy range between 0 and 10.5 MeV, and also, we investigate the axial variation of the neutron flux along the plexiglass region in the horizontal irradiation channels. Our simulation results show that a significant number of thermal neutrons as well as a non-negligible population of fast neutrons propagate through the plexiglass bar, thereby demonstrating the opportunity to product from a wide range of energy spectrum for the neutron activation analysis.


title: Room Temperature Terahertz Spectrometer with Quantum-Level Sensitivity

abstract: Photon detection with quantum-level sensitivity is particularly challenging in the terahertz regime (0.1-10 THz), which contains ~98% of all the photons existing in the universe. Near-quantum-limited terahertz spectrometry has so far only been possible through the use of cryogenically cooled superconducting mixers as frequency downconverters. Here we introduce a spectrometry scheme that uses plasmonic photomixing for frequency downconversion to offer quantum-level sensitivities at room temperature for the first time. Frequency downconversion is achieved by mixing terahertz radiation and a heterodyning optical beam with a terahertz beat frequency in a plasmonics-enhanced semiconductor active region. We demonstrate spectrometer sensitivities down to 3 times the quantum-limit at room temperature. With a versatile design capable of broadband spectrometry, this plasmonic photomixer has broad applicability to quantum optics, chemical sensing, biological studies, medical diagnosis, high data-rate communication, as well as astronomy and atmospheric studies.


title: Distance difference representations of Riemannian manifolds

abstract: Let $M$ be a complete Riemannian manifold and $F\subset M$ a set with a nonempty interior. For every $x\in M$, let $D_x$ denote the function on $F\times F$ defined by $D_x(y,z)=d(x,y)-d(x,z)$ where $d$ is the geodesic distance in $M$. The map $x\mapsto D_x$ from $M$ to the space of continuous functions on $F\times F$, denoted by $\mathcal D_F$, is called a distance difference representation of $M$. This representation, introduced recently by M. Lassas and T. Saksala, is motivated by geophysical imaging among other things. We prove that the distance difference representation $\mathcal D_F$ is a locally bi-Lipschitz homeomorphism onto its image $\mathcal D_F(M)$ and that for every open set $U\subset M$ the set $\mathcal D_F(U)$ uniquely determines the Riemannian metric on $U$. Furthermore the determination of $M$ from $\mathcal D_F(M)$ is stable if $M$ has a priori bounds on its diameter, curvature, and injectivity radius. This extends and strengthens earlier results by M. Lassas and T. Saksala.


title: SMHD: A Large-Scale Resource for Exploring Online Language Usage for Multiple Mental Health Conditions

abstract: Mental health is a significant and growing public health concern. As language usage can be leveraged to obtain crucial insights into mental health conditions, there is a need for large-scale, labeled, mental health-related datasets of users who have been diagnosed with one or more of such conditions. In this paper, we investigate the creation of high-precision patterns to identify self-reported diagnoses of nine different mental health conditions, and obtain high-quality labeled data without the need for manual labelling. We introduce the SMHD (Self-reported Mental Health Diagnoses) dataset and make it available. SMHD is a novel large dataset of social media posts from users with one or multiple mental health conditions along with matched control users. We examine distinctions in users' language, as measured by linguistic and psychological variables. We further explore text classification methods to identify individuals with mental conditions through their language.


title: Analysis of Search Stratagem Utilisation

abstract: In Interactive IR, researchers consider the user behaviour towards systems and search tasks in order to adapt search results and to improve the search experience of users. Analysing the users' past interactions with the system is one typical approach. In this paper, we analyse the user behaviour in retrieval sessions towards Marcia Bates' search stratagems such as Footnote Chasing, Citation Searching, Keyword Searching, Author Searching and Journal Run in a real-life academic search engine. In fact, search stratagems represent high-level search behaviour as the users go beyond simple execution of queries and investigate more of the system functionalities. We performed analyses of these five search stratagems using two datasets extracted from the social sciences search engine sowiport. A specific focus was the detection of the search phase and frequency of the usage of these stratagems. In addition, we explored the impact of these stratagems on the whole search process performance. We addressed mainly the usage patterns' observation of the stratagems, their impact on the conduct of retrieval sessions and explore whether they are used similarly in both datasets. From the observation and metrics proposed, we can conclude that the utilisation of search stratagems in real retrieval sessions leads to an improvement of the precision in terms of positive interactions. However, the difference is that Footnote Chasing, Citation Searching and Journal Run appear mostly at the end of a session while Keyword and Author Searching appear typically at the beginning. Thus, we can conclude from the log analysis that the improvement of search functionalities including personalisation and/or recommendation could be achieved by considering references, citations, and journals in the ranking process.


title: The fibering method approach for a non-linear Schr\"odinger equation coupled with the electromagnetic field

abstract: We study, with respect to the parameter $q\neq0$, the following Schr\"odinger-Bopp-Podolsky system in $\mathbb R^{3}$ \begin{equation*} \left\{ \begin{aligned} -&\Delta u+\omega u+q^2\phi u=|u|^{p-2}u, \\ &-\Delta \phi+a^2\Delta^2 \phi = 4\pi u^2, \end{aligned} \right. \end{equation*} where $p\in(2,3], \omega>0, a\geq0$ are fixed. We prove, by means of the fibering approach, that the system has no solutions at all for large values of $q's$, and has two radial solutions for small $q's$. We give also qualitative properties about the energy level of the solutions and a variational characterization of these extremals values of $q$. Our results recover and improve some results in the literature.


title: Robert Gibson and teaching philosophy of nature in 18th century Ireland

abstract: The dissemination of natural philosophy in the 18th-century, which was based primarily on Newton's pioneering work in mechanics, optics and astrophysics, is presented as seen through a remarkable textbook written by a little known Irish mathematics teacher, Robert Gibson. Later, he became the deputy surveyor general of Ireland (1752-1760).


title: How much income inequality is fair? Nash bargaining solution and its connection to entropy

abstract: The question about fair income inequality has been an important open question in economics and in political philosophy for over two centuries with only qualitative answers such as the ones suggested by Rawls, Nozick, and Dworkin. We provided a quantitative answer recently, for an ideal free-market society, by developing a game-theoretic framework that proved that the ideal inequality is a lognormal distribution of income at equilibrium. In this paper, we develop another approach, using the Nash Bargaining Solution (NBS) framework, which also leads to the same conclusion. Even though the conclusion is the same, the new approach, however, reveals the true nature of NBS, which has been of considerable interest for several decades. Economists have wondered about the economic meaning or purpose of the NBS. While some have alluded to its fairness property, we show more conclusively that it is all about fairness. Since the essence of entropy is also fairness, we see an interesting connection between the Nash product and entropy for a large population of rational economic agents.


title: Four-dimensional topological lattices through connectivity

abstract: Thanks to recent advances, the 4D quantum Hall (QH) effect is becoming experimentally accessible in various engineered set-ups. In this paper, we propose a new type of 4D topological system that, unlike other 2D and 4D QH models, does not require complicated (artificial) gauge fields and/or time-reversal symmetry breaking. Instead, we show that there are 4D QH systems that can be engineered for spinless particles by designing the lattice connectivity with real-valued hopping amplitudes, and we explain how this physics can be intuitively understood in analogy with the 2D Haldane model. We illustrate our discussion with a specific 4D lattice proposal, inspired by the widely-studied 2D honeycomb and brickwall lattice geometries. This also provides a minimal model for a topological system in Class AI, which supports nontrivial topological band invariants only in four spatial dimensions or higher.


title: Search for supersymmetric partners of electrons and muons in proton-proton collisions at $\sqrt{s} =$ 13 TeV

abstract: A search for direct production of the supersymmetric (SUSY) partners of electrons or muons is presented in final states with two opposite-charge, same-flavor leptons (electrons and muons), no jets, and large missing transverse momentum. The data sample corresponds to an integrated luminosity of 35.9 fb$^{-1}$ of proton-proton collisions at $\sqrt{s} =$ 13 TeV, collected with the CMS detector at the LHC in 2016. The search uses the $m_{\mathrm{T}2}$ variable, which generalizes the transverse mass for systems with two invisible objects and provides a discrimination against standard model backgrounds containing W bosons. The observed yields are consistent with the expectations from the standard model. The search is interpreted in the context of simplified SUSY models and probes slepton masses up to approximately 290, 400, and 450 GeV, assuming right-handed only, left-handed only, and both right- and left-handed sleptons (mass degenerate selectrons and smuons), and a massless lightest supersymmetric particle. Limits are also set on selectrons and smuons separately. These limits show an improvement on the existing limits of approximately 150 GeV.


title: A Hybrid RF-VLC System for Energy Efficient Wireless Access

abstract: In this paper, we propose a new paradigm in designing and realizing energy efficient wireless indoor access networks, namely, a hybrid system enabled by traditional RF access, such as WiFi, as well as the emerging visible light communication (VLC). VLC facilitates the great advantage of being able to jointly perform illumination and communications, and little extra power beyond illumination is required to empower communications, thus rendering wireless access with almost zero power consumption. On the other hand, when illumination is not required from the light source, the energy consumed by VLC could be more than that consumed by the RF. By capitalizing on the above properties, the proposed hybrid RF-VLC system is more energy efficient and more adaptive to the illumination conditions than the individual VLC or RF systems. To demonstrate the viability of the proposed system, we first formulate the problem of minimizing the power consumption of the hybrid RF-VLC system while satisfying the users requests and maintaining acceptable level of illumination, which is NP-complete. Therefore, we divide the problem into two subproblems. In the first subproblems, we determine the set of VLC access points (AP) that needs to be turned on to satisfy the illumination requirements. Given this set, we turn our attention to satisfying the users' requests for real-time communications, and we propose a randomized online algorithm that, against an oblivious adversary, achieves a competitive ratio of $\log(N)\log(M)$ with probability of success $1 - \frac{1}{N}$, where $N$ is the number of users and $M$ is the number of VLC and RF APs. We also show that the best online algorithm to solve this problem can achieve a competitive ratio of $\log(M)$. Simulation results further demonstrate the advantages of the hybrid system.


title: Statistical Significance of CP Violation in Long Baseline Neutrino Experiments

abstract: The p-value or statistical significance of a CP conservation null hypothesis test is determined from counting electron neutrino and antineutrino appearance oscillation events. The statistical estimates include cases with background events and different data sample sizes, graphical plots to interpret results and methods to combine p-values from different experiments. These estimates are useful for optimizing the search for CP violation with different amounts of neutrino and antineutrino beam running, comparing results from different experiments and for simple cross checks of more elaborate statistical estimates that use likelihood fitting of neutrino parameters.


title: A group with at least subexponential hyperlinear profile

abstract: The hyperlinear profile of a group measures the growth rate of the dimension of unitary approximations to the group. We construct a finitely-presented group whose hyperlinear profile is at least subexponential, i.e. at least $\exp(1/\epsilon^{k})$ for some $0 < k < 1$. We use this group to give an example of a two-player non-local game requiring subexponential Hilbert space dimension to play near-perfectly.


title: Factorization in weak products of complete Pick spaces

abstract: Let $\mathcal H$ be a reproducing kernel Hilbert space with a normalized complete Nevanlinna-Pick (CNP) kernel. We prove that if $(f_n)$ is a sequence of functions in $\mathcal H$ with $\sum\|f_n\|^2<\infty$, then there exists a contractive column multiplier $(\varphi_n)$ of $\mathcal H$ and a cyclic vector $F\in \mathcal H$ so that $\varphi_ n F=f_n$ for all $n$. The space of weak products $\mathcal H\odot\mathcal H$ is the set of functions of the form $h=\sum_{i=1}^\infty f_ig_i$ with $f_i, g_i\in\mathcal H$ and $\sum_{i=1}^\infty \|f_i\|\|g_i\|<\infty$. Using the above result, in combination with a recent result of Aleman, Hartz, McCarthy, and Richter, we show that for a large class of CNP spaces (including the Drury-Arveson spaces $H^2_d$ and the Dirichlet space in the unit disk) every $h\in\mathcal H\odot\mathcal H$ can be factored as a single product $h=fg$ with $f,g\in\mathcal H$.


title: Online Self-supervised Scene Segmentation for Micro Aerial Vehicles

abstract: Recently, there have been numerous advances in the development of payload and power constrained lightweight Micro Aerial Vehicles (MAVs). As these robots aspire for high-speed autonomous flights in complex dynamic environments, robust scene understanding at long-range becomes critical. The problem is heavily characterized by either the limitations imposed by sensor capabilities for geometry-based methods, or the need for large-amounts of manually annotated training data required by data-driven methods. This motivates the need to build systems that have the capability to alleviate these problems by exploiting the complimentary strengths of both geometry and data-driven methods. In this paper, we take a step in this direction and propose a generic framework for adaptive scene segmentation using self-supervised online learning. We present this in the context of vision-based autonomous MAV flight, and demonstrate the efficacy of our proposed system through extensive experiments on benchmark datasets and real-world field tests.


title: The Smirnov classes for the Fock space and complete Pick spaces

abstract: For a Hilbert function space $\mathcal H$ the Smirnov class $\mathcal N^+(\mathcal H)$ is defined to be the set of functions expressible as a ratio of bounded multipliers of $\mathcal H$, whose denominator is cyclic for the action of $Mult(\mathcal H)$. It is known that for spaces $\mathcal H$ with complete Nevanlinna-Pick (CNP) kernel, the inclusion $\mathcal H\subset \mathcal N^+(\mathcal H)$ holds. We give a new proof of this fact, which includes the new conclusion that every $h\in\mathcal H$ can be expressed as a ratio $b/a\in\mathcal N^+(\mathcal H)$ with $1/a$ already belonging to $\mathcal H$. The proof for CNP kernels is based on another Smirnov-type result of independent interest. We consider the Fock space $\mathfrak F^2_d$ of free (non-commutative) holomorphic functions and its algebra of bounded (left) multipliers $\mathfrak F^\infty_d$. We introduce the (left) {\em free Smirnov class} $\mathcal N^+_{left}$ and show that every $H \in \mathfrak F^2_d$ belongs to it. The proof of the Smirnov theorem for CNP kernels is then obtained by lifting holomorphic functions on the ball to free holomorphic functions, and applying the free Smirnov theorem.


title: Large monochromatic components in multicolored bipartite graphs

abstract: It is well-known that in every $r$-coloring of the edges of the complete bipartite graph $K_{m,n}$ there is a monochromatic connected component with at least ${m+n\over r}$ vertices. In this paper we study an extension of this problem by replacing complete bipartite graphs by bipartite graphs of large minimum degree. We conjecture that in every $r$-coloring of the edges of an $(X,Y)$-bipartite graph with $|X|=m$, $|Y|=n$, $\delta(X,Y) > \left( 1 - \frac{1}{r+1}\right) n$ and $\delta(Y,X) > \left( 1 - \frac{1}{r+1}\right) m$, there exists a monochromatic component on at least $\frac{m+n}{r}$ vertices (as in the complete bipartite graph). If true, the minimum degree condition is sharp (in that both inequalities cannot be made weak when $m$ and $n$ are divisible by $r+1$). We prove the conjecture for $r=2$ and we prove a weaker bound for all $r\geq 3$. As a corollary, we obtain a result about the existence of monochromatic components with at least $\frac{n}{r-1}$ vertices in $r$-colored graphs with large minimum degree.


title: Benchmarks for Image Classification and Other High-dimensional Pattern Recognition Problems

abstract: A good classification method should yield more accurate results than simple heuristics. But there are classification problems, especially high-dimensional ones like the ones based on image/video data, for which simple heuristics can work quite accurately; the structure of the data in such problems is easy to uncover without any sophisticated or computationally expensive method. On the other hand, some problems have a structure that can only be found with sophisticated pattern recognition methods. We are interested in quantifying the difficulty of a given high-dimensional pattern recognition problem. We consider the case where the patterns come from two pre-determined classes and where the objects are represented by points in a high-dimensional vector space. However, the framework we propose is extendable to an arbitrarily large number of classes. We propose classification benchmarks based on simple random projection heuristics. Our benchmarks are 2D curves parameterized by the classification error and computational cost of these simple heuristics. Each curve divides the plane into a "positive- gain" and a "negative-gain" region. The latter contains methods that are ill-suited for the given classification problem. The former is divided into two by the curve asymptote; methods that lie in the small region under the curve but right of the asymptote merely provide a computational gain but no structural advantage over the random heuristics. We prove that the curve asymptotes are optimal (i.e. at Bayes error) in some cases, and thus no sophisticated method can provide a structural advantage over the random heuristics. Such classification problems, an example of which we present in our numerical experiments, provide poor ground for testing new pattern classification methods.


title: A theory of maximum likelihood for weighted infection graphs

abstract: We study the problem of parameter estimation based on infection data from an epidemic outbreak on a graph. We assume that successive infections occur via contagion; i.e., transmissions can only spread across existing directed edges in the graph. Our stochastic spreading model allows individual nodes to be infected more than once, and the probability of the transmission spreading across a particular edge is proportional to both the cumulative number of times the source nodes has been infected in previous stages of the epidemic and the weight parameter of the edge. We propose a maximum likelihood estimator for inferring the unknown edge weights when full information is available concerning the order and identity of successive edge transmissions. When the weights take a particular form as exponential functions of a linear combination of known edge covariates, we show that maximum likelihood estimation amounts to optimizing a convex function, and produces a solution that is both consistent and asymptotically normal. Our proofs are based on martingale convergence theorems and the theory of weighted P\'{o}lya urns. We also show how our theory may be generalized to settings where the weights are not exponential. Finally, we analyze the case where the available infection data comes in the form of an unordered set of edge transmissions. We propose two algorithms for weight parameter estimation in this setting and derive corresponding theoretical guarantees. Our methods are validated using both synthetic data and real-world data from the Ebola spread in West Africa.


title: Influence Spatial Degeneracy on Rotational Spectroscopy Three Wave Mixing and enantiomeric state separation of Chiral Molecules

abstract: Pulse flip angles are calculated for three wave mixing, three state cycles of chiral molecules to produce optimized free induction decay amplitudes proportional to the enantiomeric excess of a sample, and to produced optimized degree of state specific enantiomeric separation. These calculations account for the spatial degeneracy of the levels involved and the resulting inhomogeneous distribution of transition dipole moments. It is found that cycles of transitions that include R followed by Q followed by P branch transitions display only modest reductions of the calculated optimal signals if spatial degeneracy is ignored. Transitions cycles P - Q - R are only slightly worse, while the Q - Q - Q cycles are much worse, increasingly so as the rotational total quantum number increases.


title: The "Hot Spots" Conjecture on the Vicsek Set

abstract: We prove the Hot Spot conjecture on the Vicsek set. Specifically, we show that every eigenfunction of the second smallest eigenvalue of the Neumann Laplacian on the Vicsek set attains its maximum and minimum on the boundary.


title: High-efficiency measurement of an artificial atom embedded in a parametric amplifier

abstract: A crucial limit to measurement efficiencies of superconducting circuits comes from losses involved when coupling to an external quantum amplifier. Here, we realize a device circumventing this problem by directly embedding a two-level artificial atom, comprised of a transmon qubit, within a flux-pumped Josephson parametric amplifier. Surprisingly, this configuration is able to enhance dispersive measurement without exposing the qubit to appreciable excess backaction. This is accomplished by engineering the circuit to permit high-power operation that reduces information loss to unmonitored channels associated with the amplification and squeezing of quantum noise. By mitigating the effects of off-chip losses downstream, the on-chip gain of this device produces end-to-end measurement efficiencies of up to 80 percent. Our theoretical model accurately describes the observed interplay of gain and measurement backaction, and delineates the parameter space for future improvement. The device is compatible with standard fabrication and measurement techniques, and thus provides a route for definitive investigations of fundamental quantum effects and quantum control protocols.


title: The phase-polarization curve of asteroid (3200) Phaethon

abstract: A multi-colour phase-polarization curve of asteroid (3200)~Phaethon has been obtained during the December 2017 apparition by merging measurements taken at the observing station of Calern (France) and at the Rhozen observatory (Bulgaria). All the observations were obtained in the positive polarization branch, the phase angle ranging from 36$^\circ$ to 116$^\circ$. The measured values of linear polarization are among the highest ever observed for a Solar system body. The covered interval of phase angle was not sufficiently extended to derive a firm determination of the $P_{\rm max}$ parameter, but this appears to occur at a phase angle around 130$^\circ$ and reaches more than 45\% of linear polarization. Phaethon is the parent body of the Geminid meteor shower, and the real physical nature of this object (asteroid or comet) has been a long-debated subject. Our polarimetric measurements seem to support the asteroid hypothesis with a phase-polarization curve similar to the asteroid (2)~Pallas, but further observations at smaller phase angles are needed to draw definitive conclusions.


title: Observation of localized modes at effective gauge field interface in synthetic photonic lattice

abstract: We predict a generic mechanism of wave localization at an interface between uniform gauge fields, arising due to propagation-dependent phase accumulation similar to Aharonov-Bohm phenomenon. We realize experimentally a photonic mesh lattice with real-time control over the vector gauge field, and observe robust localization under a broad variation of gauge strength and direction, as well as structural lattice parameters. This suggests new possibilities for confining and guiding waves in diverse physical systems through the synthetic gauge fields.


title: Kuramoto model for excitation-inhibition-based oscillations

abstract: The Kuramoto model (KM) is a theoretical paradigm for investigating the emergence of rhythmic activity in large populations of oscillators. A remarkable example of rhythmogenesis is the feedback loop between excitatory (E) and inhibitory (I) cells in large neuronal networks. Yet, although the EI-feedback mechanism plays a central role in the generation of brain oscillations, it remains unexplored whether the KM has enough biological realism to describe it. Here we derive a two-population KM that fully accounts for the onset of EI-based neuronal rhythms and that, as the original KM, is analytically solvable to a large extent. Our results provide a powerful theoretical tool for the analysis of large-scale neuronal oscillations.


title: Augmented Lagrangian-Based Decomposition Methods with Non-Ergodic Optimal Rates

abstract: We develop two new variants of alternating direction methods of multipliers (ADMM) and two parallel primal-dual decomposition algorithms to solve a wide range class of constrained convex optimization problems. Our approach relies on a novel combination of the augmented Lagrangian framework, partial alternating/linearization scheme, Nesterov's acceleration technique, and adaptive strategy. The proposed algorithms have the following new features compared to existing ADMM variants. First, they have a Nesterov's acceleration step on the primal variables instead of the dual ones as in several ADMM variants. Second, they possess an optimal $O(\frac{1}{k})$-convergence rate guarantee in a non-ergodic sense without any smoothness or strong convexity-type assumption, where $k$ is the iteration counter. When one objective term is strongly convex, our algorithm achieves an optimal $O(\frac{1}{k^2})$-non-ergodic rate. Third, our methods have better per-iteration complexity than standard ADMM due to the linearization step in the second subproblem. Fourth, we provide a set of conditions to derive update rules for algorithmic parameters, and give a concrete update for these parameters as an example. Finally, when the objective function is separable, our methods can naturally be implemented in a parallel fashion. We also study two extensions of our methods and a connection to existing primal-dual methods. We verify our theoretical development via different numerical examples and compare our methods with some existing state-of-the-art algorithms.


title: Full Bayesian Modeling for fMRI Group Analysis

abstract: Functional magnetic resonance imaging or functional MRI (fMRI) is a non-invasive way to assess brain activity by detecting changes associated with blood flow. In this work, we propose a full Bayesian procedure to analyze fMRI data for individual and group stages. For the individual stage we use a multivariate dynamic linear model (MDLM), where the temporal dependence is modeled through the state parameters and the spatial dependence is modeled only locally, taking the nearest neighbors of each voxel location. For the group stage we take advantage of the posterior distribution of the state parameters obtained in the individual stage and create a new posterior distribution that represents the updated beliefs for the group analysis. Since the posterior distribution for the state parameters is indexed by the time $t$, we propose an algorithm that allows on-line estimated curves of the state parameters to be drawn and posterior probabilities computed in order to assess brain activation for both individual and group analysis. We propose an alternative analysis for the group stage using a Gaussian process ANOVA model, where the on-line estimated curves obtained in the individual stage are modeled as a functional response. Finally, we assess our proposed modeling procedure using real resting-state data and computing empirical false-positive brain activation rates.


title: Limiting Behaviors of High Dimensional Stochastic Spin Ensembles

abstract: Lattice spin models in statistical physics are used to understand magnetism. Their Hamiltonians are a discrete form of a version of a Dirichlet energy, signifying a relationship to the Harmonic map heat flow equation. The Gibbs distribution, defined with this Hamiltonian, is used in the Metropolis-Hastings (M-H) algorithm to generate dynamics tending towards an equilibrium state. In the limiting situation when the inverse temperature is large, we establish the relationship between the discrete M-H dynamics and the continuous Harmonic map heat flow associated with the Hamiltonian. We show the convergence of the M-H dynamics to the Harmonic map heat flow equation in two steps: First, with fixed lattice size and proper choice of proposal size in one M-H step, the M-H dynamics acts as gradient descent and will be shown to converge to a system of Langevin stochastic differential equations (SDE). Second, with proper scaling of the inverse temperature in the Gibbs distribution and taking the lattice size to infinity, it will be shown that this SDE system converges to the deterministic Harmonic map heat flow equation. Our results are not unexpected, but show remarkable connections between the M-H steps and the SDE Stratonovich formulation, as well as reveal trajectory-wise out of equilibrium dynamics to be related to a canonical PDE system with geometric constraints.


title: An unbiased approach to compressed sensing

abstract: In compressed sensing a sparse vector is approximately retrieved from an under-determined equation system $Ax=b$. Exact retrieval would mean solving a large combinatorial problem which is well known to be NP-hard. For $b$ of the form $Ax_0+\epsilon$ where $x_0$ and $\epsilon$ is noise, the `oracle solution' is the one you get if you a priori know the support of $x_0$, and is the best solution one could hope for. We provide a non-convex functional whose global minimum is the oracle solution, with the property that any other local minimizer necessarily has high cardinality. We provide estimates of the type $\|\hat x-x_0\|_2\leq C\|\epsilon\|_2$ with constants $C$ that are significantly lower than for competing methods or theorems, and our theory relies on soft assumptions on the matrix $A$, in comparison with standard results in the field. The framework also allows to incorporate a priori information on the cardinality of the sought vector. In this case we show that despite being non-convex, our cost functional has no spurious local minima and the global minima is again the `oracle solution', thereby providing the first method which is guaranteed to find this point for reasonable levels of noise, without resorting to combinatorial methods.


title: How Predictable is Your State? Leveraging Lexical and Contextual Information for Predicting Legislative Floor Action at the State Level

abstract: Modeling U.S. Congressional legislation and roll-call votes has received significant attention in previous literature. However, while legislators across 50 state governments and D.C. propose over 100,000 bills each year, and on average enact over 30% of them, state level analysis has received relatively less attention due in part to the difficulty in obtaining the necessary data. Since each state legislature is guided by their own procedures, politics and issues, however, it is difficult to qualitatively asses the factors that affect the likelihood of a legislative initiative succeeding. Herein, we present several methods for modeling the likelihood of a bill receiving floor action across all 50 states and D.C. We utilize the lexical content of over 1 million bills, along with contextual legislature and legislator derived features to build our predictive models, allowing a comparison of the factors that are important to the lawmaking process. Furthermore, we show that these signals hold complementary predictive power, together achieving an average improvement in accuracy of 18% over state specific baselines.


title: A Flexible Convolutional Solver with Application to Photorealistic Style Transfer

abstract: We propose a new flexible deep convolutional neural network (convnet) to perform fast visual style transfer. In contrast to existing convnets that address the same task, our architecture derives directly from the structure of the gradient descent originally used to solve the style transfer problem [Gatys et al., 2016]. Like existing convnets, ours approximately solves the original problem much faster than the gradient descent. However, our network is uniquely flexible by design: it can be manipulated at runtime to enforce new constraints on the final solution. In particular, we show how to modify it to obtain a photorealistic result with no retraining. We study the modifications made by [Luan et al., 2017] to the original cost function of [Gatys et al., 2016] to achieve photorealistic style transfer. These modifications affect directly the gradient descent and can be reported on-the-fly in our network. These modifications are possible as the proposed architecture stems from unrolling the gradient descent.


title: Bounds on the localization number

abstract: We consider the localization game played on graphs, wherein a set of cops attempt to determine the exact location of an invisible robber by exploiting distance probes. The corresponding optimization parameter for a graph $G$ is called the localization number and is written $\zeta (G)$. We settle a conjecture of \cite{nisse1} by providing an upper bound on the localization number as a function of the chromatic number. In particular, we show that every graph with $\zeta (G) \le k$ has degeneracy less than $3^k$ and, consequently, satisfies $\chi(G) \le 3^{\zeta (G)}$. We show further that this degeneracy bound is tight. We also prove that the localization number is at most 2 in outerplanar graphs, and we determine, up to an additive constant, the localization number of hypercubes.


title: Asymptotic distribution of least square estimators for linear models with dependent errors

abstract: In this paper, we consider the usual linear regression model in the case where the error process is assumed strictly stationary. We use a result from Hannan (1973), who proved a Central Limit Theorem for the usual least square estimator under general conditions on the design and on the error process. Whatever the design satisfying Hannan's conditions, we define an estimator of the covariance matrix and we prove its consistency under very mild conditions. As an application, we show how to modify the usual tests on the linear model in this dependent context, in such a way that the type-I error rate remains asymptotically correct, and we illustrate the performance of this procedure through different sets of simulations.


title: Charge-dependent Flow Induced by Magnetic and Electric Fields in Heavy Ion Collisions

abstract: We investigate the charge-dependent flow induced by magnetic and electric fields in heavy ion collisions. We simulate the evolution of the expanding cooling droplet of strongly coupled plasma hydrodynamically, using the iEBE-VISHNU framework, and add the magnetic and electric fields as well as the electric currents they generate in a perturbative fashion. We confirm the previously reported effect of the electromagnetically induced currents, that is a charge-odd directed flow $\Delta v_1$ that is odd in rapidity, noting that it is induced by magnetic fields (\`a la Faraday and Lorentz) and by electric fields (the Coulomb field from the charged spectators). In addition, we find a charge-odd $\Delta v_3$ that is also odd in rapidity and that has a similar physical origin. We furthermore show that the electric field produced by the net charge density of the plasma drives rapidity-even charge-dependent contributions to the radial flow $\langle p_T \rangle$ and the elliptic flow $\Delta v_2$. Although their magnitudes are comparable to the charge-odd $\Delta v_1$ and $\Delta v_3$, they have a different physical origin, namely the Coulomb forces within the plasma.


title: Regimes of magnetic reconnection in colliding laser-produced magnetized plasma bubbles

abstract: We conduct a multiparametric study of driven magnetic reconnection relevant to recent experiments on colliding magnetized laser produced plasmas using particle-in-cell simulations. Varying the background plasma density, plasma resistivity, and plasma bubble geometry, the 2D simulations demonstrate a rich variety of reconnection behavior and show the coupling between magnetic reconnection and the global hydrodynamical evolution of the system. We consider both the collision between two radially expanding bubbles where reconnection is seeded by the pre-existing X-point, and the collision between two flows in a quasi-1D geometry with initially anti-parallel fields where reconnection must be initiated by the tearing instability. In both geometries, at a baseline case of low-collisionality and low background density, the current sheet is strongly compressed to below scale of the ion-skin-depth scale, and rapid, multi-plasmoid reconnection results. Increasing the plasma resistivity, we observe a collisional slow-down of reconnection and stabilization of plasmoid instability for Lundquist numbers less than approximately $S \sim 10^3$. Secondly, increasing the background plasma density modifies the compressibility of the plasma and can also slow-down or even prevent reconnection, even in completely collisionless regimes, by preventing the current sheet from thinning down to the scale of the ion-skin depth. These results have implications for understanding recent and future experiments, and signatures for these processes for proton-radiography diagnostics of these experiments are discussed.


title: Metal-adsorbed graphene nanoribbons

abstract: The metal atoms, the alkali ones excepted, might provide the multiple outermost orbitals for the multi-orbital hybridizations with the out-of-plane $\pi$ bondings on the honeycomb lattice. This will dominate the fundamental properties of Al-, Ti- and Bi-adsorbed graphene nanoribbons, in which they are explored thoroughly by using the first-principles calculations. The principle focuses are the adatom-dependent binding energies, the adatom-carbon lengths, the optimal position, the maximum adatom concentrations, the free electron density transferred per adatom, the adatom-related valence and conduction bands, the various van Hove singularities in DOSs, the transition-metal-induced magnetic properties, and the significant competitions of the zigzag edge carbons and the metal/transition metal adatoms in spin configurations. The distinct chemical bondings are clearly identified from three kinds of metal adatoms under the delicate physical quantities. The important differences between Al-/Ti-/Bi- and alkali-adsorbed graphene nanoribbons will be discussed in detail, covering band structures, relation of conduction electron density and adatom concentration, spatial charge distributions, orbital-decomposed DOSs, and magnetic configurations $\&$ moments.


title: On the variability of the solar mean magnetic field: contributions from various magnetic features on the surface of the Sun

abstract: The solar mean magnetic field (SMMF) is referred to as the disc-averaged line-of-sight (LOS) magnetic field that also reflects the polarity imbalance of the magnetic field on the Sun. The origin of the SMMF has been debated over the past few decades, with one school of thought suggesting that the contribution to the SMMF is mostly due to the large scale magnetic field structure, also termed as the background magnetic field, whereas other and more recent studies have indicated that active regions have a major contribution to the observed SMMF. In this paper, we re-investigate the issue of the origin of the SMMF by decomposing the solar disk into plages, networks, sunspots and background regions, thereby calculating the variation in the observed SMMF due to each of these features. We have used full disk images from SDO/{\it AIA} recorded at 1600~\AA\, for earmarking plages, networks and background regions and 4500~\AA\, images for separating the sunspots. The LOS fields corresponding to each of these regions are estimated from the co-temporal SDO/{\it HMI} full disk magnetograms. The temporal variation of the SMMF shows a near one-to-one correspondence with that of the background field regions suggesting that they constitute the major component of the observed SMMF. A linear regression analysis based on the coefficient of determination shows that the background field dominates, and accounts for 89\% of the variation in the SMMF whereas the magnetic field from the other features accounts for the rest 11\%.


title: Automatic formation of the structure of abstract machines in hierarchical reinforcement learning with state clustering

abstract: We introduce a new approach to hierarchy formation and task decomposition in hierarchical reinforcement learning. Our method is based on the Hierarchy Of Abstract Machines (HAM) framework because HAM approach is able to design efficient controllers that will realize specific behaviors in real robots. The key to our algorithm is the introduction of the internal or "mental" environment in which the state represents the structure of the HAM hierarchy. The internal action in this environment leads to changes the hierarchy of HAMs. We propose the classical Q-learning procedure in the internal environment which allows the agent to obtain an optimal hierarchy. We extends the HAM framework by adding on-model approach to select the appropriate sub-machine to execute action sequences for certain class of external environment states. Preliminary experiments demonstrated the prospects of the method.


title: Generalized framework for applying the Kelly criterion to stock markets

abstract: We develop a general framework for applying the Kelly criterion to stock markets. By supplying an arbitrary probability distribution modeling the future price movement of a set of stocks, the Kelly fraction for investing each stock can be calculated by inverting a matrix involving only first and second moments. The framework works for one or a portfolio of stocks and the Kelly fractions can be efficiently calculated. For a simple model of geometric Brownian motion of a single stock we show that our calculated Kelly fraction agrees with existing results. We demonstrate that the Kelly fractions can be calculated easily for other types of probabilities such as the Gaussian distribution and correlated multivariate assets.


title: Hessian spectrum at the global minimum of high-dimensional random landscapes

abstract: Using the replica method we calculate the mean spectral density of the Hessian matrix at the global minimum of a random $N \gg 1$ dimensional isotropic, translationally invariant Gaussian random landscape confined by a parabolic potential with fixed curvature $\mu>0$. Simple landscapes with generically a single minimum are typical for $\mu>\mu_{c}$, and we show that the Hessian at the global minimum is always {\it gapped}, with the low spectral edge being strictly positive. When approaching from above the transitional point $\mu= \mu_{c}$ separating simple landscapes from 'glassy' ones, with exponentially abundant minima, the spectral gap vanishes as $(\mu-\mu_c)^2$. For $\mu<\mu_c$ the Hessian spectrum is qualitatively different for 'moderately complex' and 'genuinely complex' landscapes. The former are typical for short-range correlated random potentials and correspond to 1-step replica-symmetry breaking mechanism. Their Hessian spectra turn out to be again gapped, with the gap vanishing on approaching $\mu_c$ from below with a larger critical exponent, as $(\mu_c-\mu)^4$. At the same time in the 'most complex' landscapes with long-ranged power-law correlations the replica symmetry is completely broken. We show that in that case the Hessian remains gapless for all values of $\mu<\mu_c$, indicating the presence of 'marginally stable' spatial directions. Finally, the potentials with {\it logarithmic} correlations share both 1RSB nature and gapless spectrum. The spectral density of the Hessian always takes the semi-circular form, up to a shift and an amplitude that we explicitly calculate.

categories: cond-mat.dis-nn cond-mat.stat-mech math-ph math.MP


title: A homological characterization for freeness of multi-arrangements

abstract: Building on work of Brandt and Terao in their study of $k$-formality, we introduce a co-chain complex associated to a multi-arrangement and prove that its cohomologies determine freeness of the associated module of multi-derivations. This provides a new homological method for determining freeness of arrangements and multi-arrangements. We work out many applications of this homological method. For instance, we prove that if a multi-arrangement is free then the underlying arrangement is $k$-formal for all $k\ge 2$. We also use this method to completely characterize freeness of certain families of multi-arrangements in moduli, showcasing how the geometry of multi-arrangements with the same intersection lattice may have considerable impact on freeness. New counter-examples to Orlik's conjecture also arise in connection to this latter analysis.


title: Multi-View Networks for Denoising of Arbitrary Numbers of Channels

abstract: We propose a set of denoising neural networks capable of operating on an arbitrary number of channels at runtime, irrespective of how many channels they were trained on. We coin the proposed models multi-view networks since they operate using multiple views of the same data. We explore two such architectures and show how they outperform traditional denoising models in multi-channel scenarios. Additionally, we demonstrate how multi-view networks can leverage information provided by additional recordings to make better predictions, and how they are able to generalize to a number of recordings not seen in training.


title: Pattern Dependence Detection using n-TARP Clustering

abstract: Consider an experiment involving a potentially small number of subjects. Some random variables are observed on each subject: a high-dimensional one called the "observed" random variable, and a one-dimensional one called the "outcome" random variable. We are interested in the dependencies between the observed random variable and the outcome random variable. We propose a method to quantify and validate the dependencies of the outcome random variable on the various patterns contained in the observed random variable. Different degrees of relationship are explored (linear, quadratic, cubic, ...). This work is motivated by the need to analyze educational data, which often involves high-dimensional data representing a small number of students. Thus our implementation is designed for a small number of subjects; however, it can be easily modified to handle a very large dataset. As an illustration, the proposed method is used to study the influence of certain skills on the course grade of students in a signal processing class. A valid dependency of the grade on the different skill patterns is observed in the data.


title: Apuntes de Redes Neuronales Artificiales

abstract: These handouts are designed for people who is just starting involved with the topic artificial neural networks. We show how it works a single artificial neuron (McCulloch & Pitt model), mathematically and graphically. We do explain the delta rule, a learning algorithm to find the neuron weights. We also present some examples in MATLAB/Octave. There are examples for classification task for lineal and non-lineal problems. At the end, we present an artificial neural network, a feed-forward neural network along its learning algorithm backpropagation. ----- Estos apuntes est\'an dise\~nados para personas que por primera vez se introducen en el tema de las redes neuronales artificiales. Se muestra el funcionamiento b\'asico de una neurona, matem\'aticamente y gr\'aficamente. Se explica la Regla Delta, algoritmo deaprendizaje para encontrar los pesos de una neurona. Tambi\'en se muestran ejemplos en MATLAB/Octave. Hay ejemplos para problemas de clasificaci\'on, para problemas lineales y no-lineales. En la parte final se muestra la arquitectura de red neuronal artificial conocida como backpropagation.


title: Geometric Shape Features Extraction Using a Steady State Partial Differential Equation System

abstract: A unified method for extracting geometric shape features from binary image data using a steady state partial differential equation (PDE) system as a boundary value problem is presented in this paper. The PDE and functions are formulated to extract the thickness, orientation, and skeleton simultaneously. The main advantages of the proposed method is that the orientation is defined without derivatives and thickness computation is not imposed a topological constraint on the target shape. A one-dimensional analytical solution is provided to validate the proposed method. In addition, two-dimensional numerical examples are presented to confirm the usefulness of the proposed method.


title: A Graphical Interactive Debugger for Distributed Systems

abstract: Designing and debugging distributed systems is notoriously difficult. The correctness of a distributed system is largely determined by its handling of failure scenarios. The sequence of events leading to a bug can be long and complex, and it is likely to include message reorderings and failures. On single-node systems, interactive debuggers enable stepping through an execution of the program, but they lack the ability to easily simulate failure scenarios and control the order in which messages are delivered. Oddity is a graphical, interactive debugger for distributed systems. It brings the power of traditional step-through debugging---fine-grained control and observation of a program as it executes---to distributed systems. It also enables exploratory testing, in which an engineer examines and perturbs the behavior of a system in order to better understand it, perhaps without a specific bug in mind. A programmer can directly control message and failure interleaving. Oddity supports time travel, allowing a developer to explore multiple branching executions of a system within a single debugging session. Above all, Oddity encourages distributed systems thinking: rather than assuming the normal case and attaching failure handling as an afterthought, distributed systems should be developed around the certainty of message loss and node failure. Graduate and undergraduate students used Oddity in two distributed systems classes. Usage tracking and qualitative surveys showed that students found Oddity useful for both debugging and exploratory testing.


title: Localization of topological charge density near $T_c$ in quenched QCD with Wilson flow

abstract: We smear quenched lattice QCD ensembles with lattice volume $32^3\times 8$ by using Wilson flow. Six ensembles at temperature near the critical temperature $T_c$ corresponding to the critical inverse coupling $\beta_c=6.06173(49)$ are used to investigate the localization of topological charge density. If the effective smearing radius of Wilson flow is large enough, the density, size and peak of Harrington-Shepard (HS) caloron-like topological lumps of ensembles are stable when $\beta\leq 6.050$, but start to change significantly when $\beta\geq6.055$. The inverse participation ratio (IPR) of topological charge density shows similar results, it begins to increase when $\beta\geq 6.055$ and is stable when $\beta\leq 6.050$. The pseudoscalar glueball mass is extracted from the topological charge density correlator (TCDC) of ensembles at $T=1.19T_c,~\textrm{and }1.36T_c$, the masses are $1.915(98)\textrm{ GeV}$ and $1.829(123)\textrm{ GeV}$ respectively, they are consistent with results from conventional methods.


title: VERA monitoring of the radio jet 3C 84 during 2007--2013: detection of non-linear motion

abstract: We present a kinematic study of the subparsec-scale radio jet of the radio galaxy 3C 84/NGC 1275 with the VLBI Exploration of Radio Astrometry (VERA) array at 22 GHz for 80 epochs from 2007 October to 2013 December. The averaged radial velocity of the bright component "C3" with reference to the radio core is found to be $0.27 pm 0.02c$ between 2007 October and 2013 December. This constant velocity of C3 is naturally explained by the advancing motion of the head of the mini-radio lobe. We also find a non-linear component in the motion of C3 with respect to the radio core. We briefly discuss possible origins of this non-linear motion.


title: Bounds on sizes of general caps in $AG(n,q)$ via the Croot-Lev-Pach polynomial method

abstract: In 2016, Ellenberg and Gijswijt employed a method of Croot, Lev, and Pach to show that a maximal cap in $AG(n, q)$ has size $O(q^{cn})$ for some $c < 1$. In this paper, we show more generally that if $S$ is a subset of $AG(n, q)$ containing no $m$ points on any $(m - 2)$- flat, then $|S| < q^{c_mn}$ for some $c_m < 1$, as long as $q$ is odd or $m$ is even.


title: RPM-Drive: A robust, safe, and reversible gene drive system that remains functional after 200+ generations

abstract: Despite the advent of several novel, synthetic gene drive mechanisms and their potential to one-day control a number of devastating diseases, among other applications, practical use of these systems remains contentious and risky. In particular, there is little in the way of empirical evidence of the long-term robustness of these synthetic systems against mutational breakdown. Rather, most existing systems are either known or predicted to be susceptible to rapid inactivation, though methodological designs continue to be refined. Here we evaluate a currently existing synthetic, underdominance-based gene drive system 200+ generations after it was first established in a laboratory colony of Drosophila melanogaster. Not only do we find that the system is still functioning as designed, we also show evidence that disruptions to the genetic construct are highly likely to be removed by natural selection, contributing to the system's robust, long-term stability. This stability appears to be a result of a fundamental relationship between ribosomal proteins (a novel target of the system) and natural cellular defenses that protect against cancer development. As far as we are aware, this is the longest continually functioning synthetic gene drive system thus verified, making it highly appropriate for additional research into its eventual suitability for field trials. Due to inherent properties of this gene drive, it is also likely to be adaptable for use in many different species. The insect lines established and used to test this system have been deposited at a Drosophila stock center, and are available to labs for further, independent testing.


title: Origins of the unidirectional spin Hall magnetoresistance in metallic bilayers

abstract: Recent studies evidence the emergence of asymmetric electron transport in layered conductors owing to the interplay between electrical conductivity, magnetization, and the spin Hall or Rashba- Edelstein effects. Here, we investigate the unidirectional magnetoresistance (UMR) caused by the current-induced spin accumulation in Co/Pt and CoCr/Pt bilayers. We identify three competing mechanisms underpinning the resistance asymmetry, namely interface and bulk spin-dependent electron scattering and electron-magnon scattering. Our measurements provide a consistent description of the current, magnetic field, and temperature dependence of the UMR and show that both positive and negative UMR can be obtained by tuning the interface and bulk spin-dependent scattering terms relative to the magnon population.


title: A Bound on the Cohomology of Quasiregularly Elliptic Manifolds

abstract: We show that a closed, connected and orientable Riemannian manifold of dimension $d$ that admits a quasiregular mapping from $\mathbb R^d$ must have bounded cohomological dimension independent of the distortion of the map. The dimension of the degree $l$ de Rham cohomology of $M$ is bounded above by $\binom{d}{l}$. This is a sharp upper bound that proves the Bonk-Heinonen conjecture. A corollary of this theorem answers an open problem posed by Gromov in 1981. He asked whether there exists a $d$-dimensional, simply connected manifold that does not admit a quasiregular map from $\mathbb R^d$. Our result gives an affirmative answer to this question.


title: Positive Grassmannian and polyhedral subdivisions

abstract: The nonnegative Grassmannian is a cell complex with rich geometric, algebraic, and combinatorial structures. Its study involves interesting combinatorial objects, such as positroids and plabic graphs. Remarkably, the same combinatorial structures appeared in many other areas of mathematics and physics, e.g., in the study of cluster algebras, scattering amplitudes, and solitons. We discuss new ways to think about these structures. In particular, we identify plabic graphs and more general Grassmannian graphs with polyhedral subdivisions induced by 2-dimensional projections of hypersimplices. This implies a close relationship between the positive Grassmannian and the theory of fiber polytopes and the generalized Baues problem. This suggests natural extensions of objects related to the positive Grassmannian.


title: Demonstration of nonpolar m-plane vertical GaN-on-GaN p-n power diodes grown on free-standing GaN substrates

abstract: This work demonstrates the first nonpolar vertical GaN on GaN pn power diodes grown on m-plane free standing substrates by MOCVD. The SEM and HRXRD results showed the good crystal quality of the homoepitaxial nonpolar structure with low defect densities. The CL result confirmed the nonpolar p GaN was of high quality with considerably reduced deep level states. At forward bias, the device showed good rectifying behaviors with a turn-on voltage of 4.0 V, an on-resistance of 2.3 mohmcm2, and a high on off ratio of 1e10. At reverse bias, the current leakage and breakdown were described by the trap assisted space charge limited current conduction mechanism, where I was proportional to V power 4.5. The critical electrical field was calculated to be 2.0 MV per cm without field plates or edge termination, which is the highest value reported on nonpolar power devices. The high performance m-plane p-n diodes can serve as key building blocks to further develop nonpolar GaN power electronics and polarization-engineering-related advanced power device structures for power conversion applications.