title: Magnetic field-tuned superconductor/insulator transition in TiN nanostrips

abstract: We have measured the electric transport properties of TiN nanostrips with different widths. At zero magnetic field the temperature dependent resistance R(T) saturates at a finite resistance towards low temperatures, which results from quantum phase slips in the narrower strips. We find that the current-voltage (I-V) characteristics of the narrowest strips are equivalent to those of small Josephson junctions. Applying a transverse magnetic field drives the devices into a reentrant insulating phase, with I-V-characteristics dual to those in the superconducting regime. The results evidence that our critically disordered superconducting nanostrips behave like small self-organized random Josephson networks.


title: What powers Hyperluminous Infrared galaxies at z~1-2?

abstract: We investigate what powers hyperluminous infrared galaxies (HyLIRGs; LIR(8-1000um)>10^13 Lsun) at z~1-2, by examining the behaviour of the infrared AGN luminosity function in relation to the infrared galaxy luminosity function. The former corresponds to emission from AGN-heated dust only, whereas the latter includes emission from dust heated by stars and AGN. Our results show that the two luminosity functions are substantially different below 10^13 Lsun but converge in the HyLIRG regime. We find that the fraction of AGN dominated sources increases with total infrared luminosity and at LIR >10^13.5 Lsun AGN can account for the entire infrared emission. We conclude that the bright end of the 1<z<2 infrared galaxy luminosity function is shaped by AGN rather than star-forming galaxies.


title: Backdrop: Stochastic Backpropagation

abstract: We introduce backdrop, a flexible and simple-to-implement method, intuitively described as dropout acting only along the backpropagation pipeline. Backdrop is implemented via one or more masking layers which are inserted at specific points along the network. Each backdrop masking layer acts as the identity in the forward pass, but randomly masks parts of the backward gradient propagation. Intuitively, inserting a backdrop layer after any convolutional layer leads to stochastic gradients corresponding to features of that scale. Therefore, backdrop is well suited for problems in which the data have a multi-scale, hierarchical structure. Backdrop can also be applied to problems with non-decomposable loss functions where standard SGD methods are not well suited. We perform a number of experiments and demonstrate that backdrop leads to significant improvements in generalization.


title: First-principles study of CO and OH adsorption on In-doped ZnO surfaces

abstract: We present a first-principles computational study of non-polar ZnO (10$\bar{1}$0) surfaces doped with indium. The calculations were performed using a model ZnO slab. The position of the In dopants was varied from deep bulk-like layers to the surface layers. It was established that the preferential location of the In atoms is at the surface by examining the dependence of the defect formation energy as well as the surface energy on In location. The adsorption sites on the surface of ZnO and the energy of adsorption of CO molecules and OH-species were determined in connection to In doping. It was found that OH has higher bonding energy to the surface than CO. The presence of In atoms at the surface of ZnO is favorable for CO adsorption, resulting in an elongation of the C-O bond and in charge transfer to the surface. The effect of CO and OH adsorption on the electronic and conduction properties of surfaces was assessed. We conclude that In-doped ZnO surfaces should present a higher electronic response upon adsorption of CO.


title: Computing the Spatial Probability of Inclusion inside Partial Contours for Computer Vision Applications

abstract: In Computer Vision, edge detection is one of the favored approaches for feature and object detection in images since it provides information about their objects boundaries. Other region-based approaches use probabilistic analysis such as clustering and Markov random fields, but those methods cannot be used to analyze edges and their interaction. In fact, only image segmentation can produce regions based on edges, but it requires thresholding by simply separating the regions into binary in-out information. Hence, there is currently a gap between edge-based and region-based algorithms, since edges cannot be used to study the properties of a region and vice versa. The objective of this paper is to present a novel spatial probability analysis that allows determining the probability of inclusion inside a set of partial contours (strokes). To answer this objective, we developed a new approach that uses electromagnetic convolutions and repulsion optimization to compute the required probabilities. Hence, it becomes possible to generate a continuous space of probability based only on the edge information, thus bridging the gap between the edge-based methods and the region-based methods. The developed method is consistent with the fundamental properties of inclusion probabilities and its results are validated by comparing an image with the probability-based estimation given by our algorithm. The method can also be generalized to take into consideration the intensity of the edges or to be used for 3D shapes. This is the first documented method that allows computing a space of probability based on interacting edges, which opens the path to broader applications such as image segmentation and contour completion.


title: Design of optimal illumination patterns in single-pixel imaging using image dictionaries

abstract: Single-pixel imaging (SPI) has a major drawback that many sequential illuminations are required for capturing one single image with long acquisition time. Basis illumination patterns such as Fourier patterns and Hadamard patterns can achieve much better imaging efficiency than random patterns. But the performance is still sub-optimal since the basis patterns are fixed and non-adaptive for varying object images. This Letter proposes a novel scheme for designing and optimizing the illumination patterns adaptively from an image dictionary by extracting the common image features using principal component analysis (PCA). Simulation and experimental results reveal that our proposed scheme outperforms conventional Fourier SPI in terms of imaging efficiency.


title: Driven power-law oscillator

abstract: We explore the nonlinear dynamics of a driven power law oscillator whose shape varies periodically in time covering a broad spectrum of anharmonicities. Combining weak and strong confinement of different geometry within a single driving period the phase space allows not only for regular and chaotic bounded motion but in particular also for an unbounded motion which exhibits an exponential net growth of the corresponding energies. Our computational study shows that phases of motion with energy gain and loss as well as approximate energy conservation alternate within a single period of the oscillator and can be assigned to the change of the underlying confinement geometry. We demonstrate how the crossover from a single- to a two-component phase space takes place with varying frequency and amplitude and analyze the corresponding volumes in phase space. In the high frequency regime an effective potential is derived that combines the different features of the driven power-law oscillator. Possible experimental realizations are discussed.


title: Asymmetry Helps: Improved Private Information Retrieval Protocols for Distributed Storage

abstract: We consider private information retrieval (PIR) for distributed storage systems (DSSs) with noncolluding nodes where data is stored using a non maximum distance separable (MDS) linear code. It was recently shown that if data is stored using a particular class of non-MDS linear codes, the MDS-PIR capacity, i.e., the maximum possible PIR rate for MDS-coded DSSs, can be achieved. For this class of codes, we prove that the PIR capacity is indeed equal to the MDS-PIR capacity, giving the first family of non-MDS codes for which the PIR capacity is known. For other codes, we provide asymmetric PIR protocols that achieve a strictly larger PIR rate compared to existing symmetric PIR protocols.


title: Contributions to ${b \rightarrow s \ell \ell}$ Anomalies from ${R}$-Parity Violating Interactions

abstract: We examine the parameter space of supersymmetric models with $R$-parity violating interactions of the form $\lambda' L Q D^c$ to explain the various anomalies observed in $b \rightarrow s \ell \ell$ transitions. To generate the appropriate operator in the low energy theory, we are led to a region of parameter space where loop contributions dominate. In particular, we concentrate on parameters for which diagrams involving winos, which have not been previously considered, give large contributions. Many different potentially constraining processes are analyzed, including $\tau \rightarrow \mu \mu \mu$, $B_s-\bar{B}_s$ mixing, $B \rightarrow K^{(*)} \nu \bar{\nu}$, $Z$ decays to charged leptons, and direct LHC searches. We find that it is possible to explain the anomalies, but it requires large values of $\lambda'$, which lead to relatively low Landau poles.


title: Data-Driven Participation Factors for Nonlinear Systems Based on Koopman Mode Decomposition

abstract: This paper develops a novel data-driven technique to compute the participation factors for nonlinear systems based on the Koopman mode decomposition. Provided that certain conditions are satisfied, it is shown that the proposed technique generalizes the original definition of the linear mode-in-state participation factors. Two numerical examples are provided to demonstrate the performance of our approach: one relying on a canonical nonlinear dynamical system, and the other based on the two-area four-machine power system. The Koopman mode decomposition is capable of coping with a large class of nonlinearity, thereby making our technique able to deal with oscillations arising in practice due to nonlinearities while being fast to compute and compatible with real-time applications.


title: Extracting relevant structures from self-determination theory questionnaires via Information Bottleneck method

abstract: In this paper, we introduce the application of Information Bottleneck as a method to investigate properties of questionnaires developed for the study of motivational profiles based on self-determination theory. Founded on information theory, the Information Bottleneck method compresses multidimensional categorical data into clusters by minimizing the loss of relevant information. It does not require linearity such as Pearson correlation. For exploratory data analysis, it allows investigating the structure of questionnaire items grouping with a hierarchical partition as perceived by the sample, based only on conditional probability distributions calculated from answers to the questionnaire items. We applied the approach to an instrument adapted from the academic motivation scale. The Information Bottleneck approach showed a partition of items aligned with the self-determination continuum. The results of this first application of Information Bottleneck to psychometric studies suggested that this procedure can could be useful as a consistent, complementary analysis for traditional exploratory procedures such as factor analysis.


title: A first-principles global multiphase equation of state for hydrogen

abstract: We present and discuss a wide-range hydrogen equation of state model based on a consistent set of ab initio simulations including quantum protons and electrons. Both the process of constructing this model and its predictions are discussed in detail. The cornerstones of this work are the specification of simple physically motivated free energy models, a general multiparameter/multiderivative fitting method, and the use of the most accurate simulation methods to date. The resulting equation of state aims for a global range of validity ($T = 1-10^9 K$ and $V_m = 10^{-9}-1 m^3/mol$), as the models are specifically constructed to reproduce exact thermodynamic and mechanical limits. Our model is for the most part analytic or semianalytic and is thermodynamically consistent by construction; the problem of interpolating between distinctly different models -often a cause for thermodynamic inconsistencies and spurious discontinuities- is avoided entirely.


title: Importance Sampling Policy Evaluation with an Estimated Behavior Policy

abstract: We consider the problem of off-policy evaluation in Markov decision processes. Off-policy evaluation is the task of evaluating the expected return of one policy with data generated by a different, behavior policy. Importance sampling is a technique for off-policy evaluation that re-weights off-policy returns to account for differences in the likelihood of the returns between the two policies. In this paper, we study importance sampling with an estimated behavior policy where the behavior policy estimate comes from the same set of data used to compute the importance sampling estimate. We find that this estimator often lowers the mean squared error of off-policy evaluation compared to importance sampling with the true behavior policy or using a behavior policy that is estimated from a separate data set. Intuitively, estimating the behavior policy in this way corrects for error due to sampling in the action-space. Our empirical results also extend to other popular variants of importance sampling and show that estimating a non-Markovian behavior policy can further lower large-sample mean squared error even when the true behavior policy is Markovian.


title: Gas-phase microresonator-based comb spectroscopy without an external pump laser

abstract: We present a novel approach to realize microresonator-comb-based high resolution spectroscopy that combines a fiber-laser cavity with a microresonator. Although the spectral resolution of a chip-based comb source is typically limited by the free spectral range (FSR) of the microresonator, we overcome this limit by tuning the 200-GHz repetition-rate comb over one FSR via control of an integrated heater. Our dual-cavity scheme allows for self-starting comb generation without the need for conventional pump-cavity detuning while achieving a spectral resolution equal to the comb linewidth. We measure broadband molecular absorption spectra of acetylene by interleaving 800 spectra taken at 250-MHz per spectral step using a 60-GHz-coarse-resolution spectrometer and exploits advances of integrated heater which can locally and rapidly change the refractive index of a microresonator with low electrical consumption (0.9 GHz/mW), which is orders of magnitude lower than a fiber-based comb. This approach offers a path towards a simple, robust and low-power consumption CMOS-compatible platform capable of remote sensing.


title: gprHOG and the popularity of Histogram of Oriented Gradients (HOG) for Buried Threat Detection in Ground-Penetrating Radar

abstract: Substantial research has been devoted to the development of algorithms that automate buried threat detection (BTD) with ground penetrating radar (GPR) data, resulting in a large number of proposed algorithms. One popular algorithm GPR-based BTD, originally applied by Torrione et al., 2012, is the Histogram of Oriented Gradients (HOG) feature. In a recent large-scale comparison among five veteran institutions, a modified version of HOG referred to here as "gprHOG", performed poorly compared to other modern algorithms. In this paper, we provide experimental evidence demonstrating that the modifications to HOG that comprise gprHOG result in a substantially better-performing algorithm. The results here, in conjunction with the large-scale algorithm comparison, suggest that HOG is not competitive with modern GPR-based BTD algorithms. Given HOG's popularity, these results raise some questions about many existing studies, and suggest gprHOG (and especially HOG) should be employed with caution in future studies.


title: Trigger and Timing Distributions using the TTC-PON and GBT Bridge Connection in ALICE for the LHC Run 3 Upgrade

abstract: The ALICE experiment at CERN is preparing for a major upgrade for the third phase of data taking run (Run 3), when the high luminosity phase of the Large Hadron Collider (LHC) starts. The increase in the beam luminosity will result in high interaction rate causing the data acquisition rate to exceed 3 TB/sec. In order to acquire data for all the events and to handle the increased data rate, a transition in the readout electronics architecture from the triggered to the trigger-less acquisition mode is required. In this new architecture, a dedicated electronics block called the Common Readout Unit (CRU) is defined to act as a nodal communication point for detector data aggregation and as a distribution point for timing, trigger and control (TTC) information. TTC information in the upgraded triggerless readout architecture uses two asynchronous high-speed serial links connections: the TTC-PON and the GBT. We have carried out a study to evaluate the quality of the embedded timing signals forwarded by the CRU to the connected electronics using the TTC-PON and GBT bridge connection. We have used four performance metrics to characterize the communication bridge: (a)the latency added by the firmware logic, (b)the jitter cleaning effect of the PLL on the timing signal, (c)BER analysis for quantitative measurement of signal quality, and (d)the effect of optical transceivers parameter settings on the signal strength. Reliability study of the bridge connection in maintaining the phase consistency of timing signals is conducted by performing multiple iterations of power on/off cycle, firmware upgrade and reset assertion/de-assertion cycle (PFR cycle). The test results are presented and discussed concerning the performance of the TTC-PON and GBT bridge communication chain using the CRU prototype and its compliance with the ALICE timing requirements.


title: Document Chunking and Learning Objective Generation for Instruction Design

abstract: Instructional Systems Design is the practice of creating of instructional experiences that make the acquisition of knowledge and skill more efficient, effective, and appealing. Specifically in designing courses, an hour of training material can require between 30 to 500 hours of effort in sourcing and organizing reference data for use in just the preparation of course material. In this paper, we present the first system of its kind that helps reduce the effort associated with sourcing reference material and course creation. We present algorithms for document chunking and automatic generation of learning objectives from content, creating descriptive content metadata to improve content-discoverability. Unlike existing methods, the learning objectives generated by our system incorporate pedagogically motivated Bloom's verbs. We demonstrate the usefulness of our methods using real world data from the banking industry and through a live deployment at a large pharmaceutical company.


title: Modeling Impact of Human Errors on the Data Unavailability and Data Loss of Storage Systems

abstract: Data storage systems and their availability play a crucial role in contemporary datacenters. Despite using mechanisms such as automatic fail-over in datacenters, the role of human agents and consequently their destructive errors is inevitable. Due to very large number of disk drives used in exascale datacenters and their high failure rates, the disk subsystem in storage systems has become a major source of Data Unavailability (DU) and Data Loss (DL) initiated by human errors. In this paper, we investigate the effect of Incorrect Disk Replacement Service (IDRS) on the availability and reliability of data storage systems. To this end, we analyze the consequences of IDRS in a disk array, and conduct Monte Carlo simulations to evaluate DU and DL during mission time. The proposed modeling framework can cope with a) different storage array configurations and b) Data Object Survivability (DOS), representing the effect of system level redundancies such as remote backups and mirrors. In the proposed framework, the model parameters are obtained from industrial and scientific reports alongside field data which have been extracted from a datacenter operating with 70 storage racks. The results show that ignoring the impact of IDRS leads to unavailability underestimation by up to three orders of magnitude. Moreover, our study suggests that by considering the effect of human errors, the conventional beliefs about the dependability of different Redundant Array of Independent Disks (RAID) mechanisms should be revised. The results show that RAID1 can result in lower availability compared to RAID5 in the presence of human errors. The results also show that employing automatic fail-over policy (using hot spare disks) can reduce the drastic impacts of human errors by two orders of magnitude.


title: Natural Language Generation for Electronic Health Records

abstract: A variety of methods existing for generating synthetic electronic health records (EHRs), but they are not capable of generating unstructured text, like emergency department (ED) chief complaints, history of present illness or progress notes. Here, we use the encoder-decoder model, a deep learning algorithm that features in many contemporary machine translation systems, to generate synthetic chief complaints from discrete variables in EHRs, like age group, gender, and discharge diagnosis. After being trained end-to-end on authentic records, the model can generate realistic chief complaint text that preserves much of the epidemiological information in the original data. As a side effect of the model's optimization goal, these synthetic chief complaints are also free of relatively uncommon abbreviation and misspellings, and they include none of the personally-identifiable information (PII) that was in the training data, suggesting it may be used to support the de-identification of text in EHRs. When combined with algorithms like generative adversarial networks (GANs), our model could be used to generate fully-synthetic EHRs, facilitating data sharing between healthcare providers and researchers and improving our ability to develop machine learning methods tailored to the information in healthcare data.


title: Long time behavior of random and nonautonomous Fisher-KPP equations. Part I. Stability of equilibria and spreading speeds

abstract: In the current series of two papers, we study the long time behavior of the following random Fisher-KPP equation $$ u_t =u_{xx}+a(\theta_t\omega)u(1-u),\quad x\in\mathbb{R} $$ where $\omega\in\Omega$, $(\Omega, \mathcal{F},\mathbb{P})$ is a given probability space, $\theta_t$ is an ergodic metric dynamical system on $\Omega$, and $a(\omega)>0$ for every $\omega\in\Omega$. We also study the long time behavior of the following nonautonomous Fisher-KPP equation, $$ u_t=u_{xx}+a_0(t)u(1-u),\quad x\in\mathbb{R}$$ where $a_0(t)$ is a positive locally H\"older continuous function. In this first part of the series, we investigate the stability of positive equilibria and the spreading speeds. Under some proper assumption on $a(\omega)$, we show that the constant solution $u=1$ of (1) is asymptotically stable with respect to strictly positive perturbations and show that (1) has a deterministic spreading speed interval $[2\sqrt{\underline a}, 2\sqrt{\bar a}]$, where $\underline{a}$ and $\bar a$ are the least and the greatest means of $a(\cdot)$, respectively, and hence the spreading speed interval is linearly determinant. It is shown that the solution of (1) with the initial function which is bounded away from $0$ for $x\ll -1$ and is $0$ for $x\gg 1$ propagates at the speed $2\sqrt {\hat a}$, where $\hat a$ is the average of $a(\cdot)$. Under some assumption on $a_0(\cdot)$, we also show that the constant solution $u=1$ of (2) is asymptotically stably and (2) admits a bounded spreading speed interval. It is not assumed that $a(\omega)$ and $a_0(t)$ are bounded above and below by some positive constants. The results obtained in this part are new and extend the existing results in literature on spreading speeds of Fisher-KPP equations. In the second part of the series, we will study the existence and stability of transition fronts of (1) and (2).


title: Steady-state generation of Wigner negative states in 1D resonance fluorescence

abstract: In this work we demonstrate numerically that the nonlinearity provided by a continuously driven two-level system (TLS) allows for the generation of Wigner-negative states of the electromagnetic field confined in one spatial dimension. Wigner-negative states, a.k.a. Wigner nonclassical states, are desirable for quantum information protocols beyond the scope of classical computers. Focusing on the steady-state emission from the TLS, we find the largest negativity at the drive strength where the coherent reflection vanishes.


title: New 3d $\mathcal{N}=2$ Dualities from Quadratic Monopoles

abstract: Aspects of three dimensional $\mathcal{N}=2$ gauge theories with monopole superpotentials and their dualities are investigated. The moduli spaces of a number of such theories are studied using Hilbert series. Moreover, we propose new dualities involving quadratic powers for the monopole superpotentials, for unitary, symplectic and orthogonal gauge groups. These dualities are then tested using the three sphere partition function and matching of the Hilbert series. We also provide an argument for the obstruction to the duality for theories with quartic monopole superpotentials.


title: Adversarial Domain Adaptation for Classification of Prostate Histopathology Whole-Slide Images

abstract: Automatic and accurate Gleason grading of histopathology tissue slides is crucial for prostate cancer diagnosis, treatment, and prognosis. Usually, histopathology tissue slides from different institutions show heterogeneous appearances because of different tissue preparation and staining procedures, thus the predictable model learned from one domain may not be applicable to a new domain directly. Here we propose to adopt unsupervised domain adaptation to transfer the discriminative knowledge obtained from the source domain to the target domain without requiring labeling of images at the target domain. The adaptation is achieved through adversarial training to find an invariant feature space along with the proposed Siamese architecture on the target domain to add a regularization that is appropriate for the whole-slide images. We validate the method on two prostate cancer datasets and obtain significant classification improvement of Gleason scores as compared with the baseline models.


title: Collision strengths for FIR and UV transtions in PIII and the phosphorus abundance

abstract: Phosphorus abundance is crucial for DNA-based extraterrestrial life in exoplanets. Atomic data for observed spectral lines of P-ions are needed for its accurate determination. We present the first calculations for collision strengths for the forbidden PIII fine structure transition $3s^23p (^2P^o_{1/2-3/2})$ within the ground state at 17.9 $\mu$m, as well as allowed UV transitions in the $3s^23p (^2P^o_{1/2,3/2}) \rightarrow 3s3p^2 (^2D_{3/2,5/2}, ^2S_{1/2}, ^2P_{1/2,3/2})$ multiplets between 915-1345 $\AA$. Collision strengths are computed using the Breit-Pauli R-Matrix method including the first 18 levels, and they exhibit extensive auto-ionizing resonance structures. In particular, the Maxwellian averaged effective collision strength for the FIR 17.9 $\mu$m transition shows a factor 3 temperature variation broadly peaking at typical nebular temperatures. Its theoretical emissivity with solar phosphorus abundance is computed relative to H$\beta$ and found to be similar to observed intensties from planetary nebulae; the abundances derived in earlier works are 3-5 times sub-solar. The results pertain to the reported paucity of phosphorus from preferred production sites in supernovae, and abundances in planetary nebulae and supernova remnants.


title: Sums of squares in pseudoconvex hypersurfaces and torsion phenomena for Catlin's boundary systems

abstract: Given a pseudoconvex hypersurface in C^n and an arbitrary weight, we show the existence of local coordinates in which the polynomial model contains a particularly simple sum of squares of monomials. Our second main result provides a normalization of a part of any Catlin boundary system. We illustrate by an example that this normalization cannot be extended to the rest of the boundary system due to the existence of what we refer to as torsion.


title: Evaluating Impact of Human Errors on the Availability of Data Storage Systems

abstract: In this paper, we investigate the effect of incorrect disk replacement service on the availability of data storage systems. To this end, we first conduct Monte Carlo simulations to evaluate the availability of disk subsystem by considering disk failures and incorrect disk replacement service. We also propose a Markov model that corroborates the Monte Carlo simulation results. We further extend the proposed model to consider the effect of automatic disk fail-over policy. The results obtained by the proposed model show that overlooking the impact of incorrect disk replacement can result up to three orders of magnitude unavailability underestimation. Moreover, this study suggests that by considering the effect of human errors, the conventional believes about the dependability of different RAID mechanisms should be revised. The results show that in the presence of human errors, RAID1 can result in lower availability compared to RAID5.


title: The Double Pentaladder Integral to All Orders

abstract: We compute dual-conformally invariant ladder integrals that are capped off by pentagons at each end of the ladder. Such integrals appear in six-point amplitudes in planar N=4 super-Yang-Mills theory. We provide exact, finite-coupling formulas for the basic double pentaladder integrals as a single Mellin integral over hypergeometric functions. For particular choices of the dual conformal cross ratios, we can evaluate the integral at weak coupling to high loop orders in terms of multiple polylogarithms. We argue that the integrals are exponentially suppressed at strong coupling. We describe the space of functions that contains all such double pentaladder integrals and their derivatives, or coproducts. This space, a prototype for the space of Steinmann hexagon functions, has a simple algebraic structure, which we elucidate by considering a particular discontinuity of the functions that localizes the Mellin integral and collapses the relevant symbol alphabet. This function space is endowed with a coaction, both perturbatively and at finite coupling, which mixes the independent solutions of the hypergeometric differential equation and constructively realizes a coaction principle of the type believed to hold in the full Steinmann hexagon function space.


title: Metallic and antiferromagnetic fixed points from gravity

abstract: We consider SU(2)$\times$ U(1) gauge theory coupled to matter field in adjoints and study RG group flow. We constructed Callan Symanzik equation and subsequent $\beta$ functions and study the fixed points. We find there are two fixed points, showing metallic and antiferromagnetic behaviour. We have shown that metallic phase develops an instability if certain parametric conditions are satisfied.


title: Playing Atari with Six Neurons

abstract: Deep reinforcement learning, applied to vision-based problems like Atari games, maps pixels directly to actions; internally, the deep neural network bears the responsibility of both extracting useful information and making decisions based on it. By separating the image processing from decision-making, one could better understand the complexity of each task, as well as potentially find smaller policy representations that are easier for humans to understand and may generalize better. To this end, we propose a new method for learning policies and compact state representations separately but simultaneously for policy approximation in reinforcement learning. State representations are generated by an encoder based on two novel algorithms: Increasing Dictionary Vector Quantization makes the encoder capable of growing its dictionary size over time, to address new observations as they appear in an open-ended online-learning context; Direct Residuals Sparse Coding encodes observations by disregarding reconstruction error minimization, and aiming instead for highest information inclusion. The encoder autonomously selects observations online to train on, in order to maximize code sparsity. As the dictionary size increases, the encoder produces increasingly larger inputs for the neural network: this is addressed by a variation of the Exponential Natural Evolution Strategies algorithm which adapts its probability distribution dimensionality along the run. We test our system on a selection of Atari games using tiny neural networks of only 6 to 18 neurons (depending on the game's controls). These are still capable of achieving results comparable---and occasionally superior---to state-of-the-art techniques which use two orders of magnitude more neurons.


title: The data-driven Schroedinger bridge

abstract: Erwin Schroedinger posed, and to a large extent solved in 1931/32 the problem of finding the most likely random evolution between two continuous probability distributions. This article considers this problem in the case when only samples of the two distributions are available. A novel iterative procedure is proposed, inspired by Fortet-Sinkhorn type algorithms. Since only samples of the marginals are available, the new approach features constrained maximum likelihood estimation in place of the nonlinear boundary couplings, and importance sampling to propagate the functions $\varphi$ and $\hat{\varphi}$ solving the Schroedinger system. This method is well-suited to high-dimensional settings, where introducing grids leads to numerically unfeasible or unreliable methods. The methodology is illustrated in two applications: entropic interpolation of two-dimensional Gaussian mixtures, and the estimation of integrals through a variation of importance sampling.


title: Stochastic and Coarse-Grained Two-Dimensional Modeling of Directional Particle Movement

abstract: We study the evolution of interacting groups of agents in two-dimensional geometries. We introduce a microscopic stochastic model that includes floor fields modeling the global flow of individual groups as well as local interaction rules. From this microscopic model we derive an analytically-tractable system of conservation laws that governs the evolution of the macroscopic densities. Numerical simulations show good agreement between the system of conservation laws and the microscopic model, though the latter is slightly more diffusive. We conclude by deriving second-order corrections to the system of conservation laws.


title: A Systematic Mapping Study on Security in Agile Requirements Engineering

abstract: [Background] The rapidly changing business environments in which many companies operate is challenging traditional Requirements Engineering (RE) approaches. This gave rise to agile approaches for RE. Security, at the same time, is an essential non-functional requirement that still tends to be difficult to address in agile development contexts. Given the fuzzy notion of "agile" in context of RE and the difficulties of appropriately handling security requirements, the overall understanding of how to handle security requirements in agile RE is still vague. [Objective] Our goal is to characterize the publication landscape of approaches that handle security requirements in agile software projects. [Method] We conducted a systematic mapping to outline relevant work and contemporary gaps for future research. [Results] In total, we identified 21 studies that met our inclusion criteria, dated from 2005 to 2017. We found that the approaches typically involve modifying agile methods, introducing new artifacts (e.g., extending the concept of user story to abuser story), or introducing guidelines to handle security issues. We also identified limitations of using these approaches related to environment, people, effort and resources. [Conclusion] Our analysis suggests that more effort needs to be invested into empirically evaluating the existing approaches and that there is an avenue for future research in the direction of mitigating the identified limitations.


title: Understanding diseases as increased heterogeneity: a complex network computational framework

abstract: Due to the complexity of the human body, most diseases present a high inter-personal variability in the way they manifest, i.e. in their phenotype, which has important clinical repercussions - as for instance the difficulty in defining objective diagnostic rules. We here explore the hypothesis that signs and symptoms used to define a disease should be understood in terms of the dispersion (as opposed to the average) of physical observables. To that end, we propose a computational framework, based on complex networks theory, to map groups of subjects to a network structure, based on their pairwise phenotypical similarity. We demonstrate that the resulting structure can be used to improve the performance of classification algorithms, especially in the case of a limited number of instances, both with synthetic and real data sets. Beyond providing an alternative conceptual understanding of diseases, the proposed framework could be of special relevance in the growing field of personalised, or N-to-1, medicine.


title: Adversarial Reinforcement Learning Framework for Benchmarking Collision Avoidance Mechanisms in Autonomous Vehicles

abstract: With the rapidly growing interest in autonomous navigation, the body of research on motion planning and collision avoidance techniques has enjoyed an accelerating rate of novel proposals and developments. However, the complexity of new techniques and their safety requirements render the bulk of current benchmarking frameworks inadequate, thus leaving the need for efficient comparison techniques unanswered. This work proposes a novel framework based on deep reinforcement learning for benchmarking the behavior of collision avoidance mechanisms under the worst-case scenario of dealing with an optimal adversarial agent, trained to drive the system into unsafe states. We describe the architecture and flow of this framework as a benchmarking solution, and demonstrate its efficacy via a practical case study of comparing the reliability of two collision avoidance mechanisms in response to intentional collision attempts.


title: Mode-Coupling Theory of the Glass Transition: A Primer

abstract: Understanding the physics of glass formation remains one of the major unsolved challenges of condensed matter science. As a material solidifies into a glass, it exhibits a spectacular slowdown of the dynamics upon cooling or compression, but at the same time undergoes only minute structural changes. Among the numerous theories put forward to rationalize this complex behavior, Mode-Coupling Theory (MCT) stands out as the only framework that provides a fully first-principles-based description of glass phenomenology. This review outlines the key physical ingredients of MCT, its predictions, successes, and failures, as well as recent improvements of the theory. We also discuss the extension and application of MCT to the emerging field of non-equilibrium active soft matter


title: ${\delta}$ MUS Revisited:Rectifying an 82 Yr-Old Mistake

abstract: The red giant ${\delta}$ Muscae is known since 1919 to be a spectroscopic binary, and the first and only orbit was determined in 1936, claiming the period to be 847 days. This was discrepant with the Hipparcos determined astrometric orbit. Using the latest data available for this object - leading to a 100 yr time-span - we show here that the correct period is 423 d, and are able for the first time to combine the spectroscopic orbit with the Hipparcos orbit. Using all the available information, we find that the $\sim$ 1.2 solar mass red giant must have a $\sim$ 0.3 - 0.4 solar mass M dwarf companion, and that the system will soon evolve towards a He WD binary system. Given its relatively short period, ${\delta}$ Muscae may be an ideal benchmark for testing astrometric orbits obtained by Gaia for very bright stars.


title: Topologically-based fractional diffusion and emergent dynamics with short-range interactions

abstract: We introduce a new class of models for emergent dynamics. It is based on a new communication protocol which incorporates two main features: short-range kernels which restrict the communication to local geometric balls, and anisotropic communication kernels, adapted to the local density in these balls, which form topological neighborhoods. We prove flocking behavior -- the emergence of global alignment for regular, non-vacuous solutions of the $n$-dimensional models based on short-range topological communication. Moreover, global regularity (and hence unconditional flocking) of the one-dimensional model is proved via an application of a De Giorgi-type method. To handle the non-symmetric singular kernels that arise with our topological communication, we develop a new analysis for local fractional elliptic operators, interesting for its own sake, encountered in the construction of our class of models.


title: The local zeta function in enumerating quartic fields

abstract: An exact formula is obtained for the Fourier transform of the local condition of maximality modulo primes $p>3$ in the prehomogeneous vector space $2 \otimes \mathrm{Sym}^2(\mathbb{Z}_p^3)$ parametrizing quartic fields, thus solving the local `quartic case' in enumerating quartic fields.


title: Nonuniqueness of conformal metrics with constant $Q$-curvature

abstract: We establish several nonuniqueness results for the problem of finding complete conformal metrics with constant (fourth-order) $Q$-curvature on compact and noncompact manifolds of dimension $\geq5$. Infinitely many branches of metrics with constant $Q$-curvature, but without constant scalar curvature, are found to bifurcate from Berger metrics on spheres and complex projective spaces. These provide examples of nonisometric metrics with the same constant negative $Q$-curvature in a conformal class with negative Yamabe invariant, echoing the absence of a Maximum Principle. We also discover infinitely many complete metrics with constant $Q$-curvature conformal to $\mathbb S^m\times\mathbb R^d$, $m\geq4$, $d\geq1$, and $\mathbb S^m\times\mathbb H^d$, $2\leq d\leq m-3$; which give infinitely many solutions to the singular constant $Q$-curvature problem on round spheres $\mathbb S^n$ blowing up along a round subsphere $\mathbb S^k$, for all $0\leq k<(n-4)/2$.


title: Improving rewards in overloaded real-time systems

abstract: Competitive analysis of online algorithms has commonly been applied to understand the behaviour of real-time systems during overload conditions. While competitive analysis provides insight into the behaviour of certain algorithms, it is hard to make inferences about the performance of those algorithms in practice. Other approaches to dealing with overload resort to heuristics that seem to perform well but are hard to prove as being good. Further, most work on handling overload in real-time systems does not consider using information regarding the distribution of arrival rates of jobs and execution times to make scheduling decisions. We present an scheduling policy (obtained through stochastic approximation, and using information about the workload) to handle overload in real-time systems and improve the revenue earned when each successful job completion results in revenue accrual. We prove that the policy we outline does lead to increased revenue when compared to a class of scheduling policies that make static resource allocations to different service classes. We also use empirical evidence to underscore the fact that this policy performs better than a variety of other scheduling policies. The ideas presented can be applied to several soft real-time systems, specifically systems with multiple service classes.


title: Atomistic Study of the Electronic Contact Resistivity Between the Half-Heusler Alloys (HfCoSb, HfZrCoSb, HfZrNiSn) and the Metal Ag

abstract: Half-Heusler(HH) alloys have shown promising thermoelectric properties in the medium and high temperature range. To harness these material properties for thermoelectric applications, it is important to realize electrical contacts with low electrical contact resistivity. However, little is known about the detailed structural and electronic properties of such contacts, and the expected values of contact resistivity. Here, we employ atomistic ab initio calculations to study electrical contacts in a subclass of HH alloys consisting of the compounds HfCoSb, HfZrCoSb, and HfZrNiSn. By using Ag as a prototypical metal, we show that the termination of the HH material critically determines the presence or absence of strong deformations at the interface. Our study includes contacts to doped materials, and the results indicate that the p-type materials generally form ohmic contacts while the n-type materials have a small Schottky barrier. We calculate the temperature dependence of the contact resistivity in the low to medium temperature range and provide quantitative values that set lower limits for these systems.


title: Factorized Adversarial Networks for Unsupervised Domain Adaptation

abstract: In this paper, we propose Factorized Adversarial Networks (FAN) to solve unsupervised domain adaptation problems for image classification tasks. Our networks map the data distribution into a latent feature space, which is factorized into a domain-specific subspace that contains domain-specific characteristics and a task-specific subspace that retains category information, for both source and target domains, respectively. Unsupervised domain adaptation is achieved by adversarial training to minimize the discrepancy between the distributions of two task-specific subspaces from source and target domains. We demonstrate that the proposed approach outperforms state-of-the-art methods on multiple benchmark datasets used in the literature for unsupervised domain adaptation. Furthermore, we collect two real-world tagging datasets that are much larger than existing benchmark datasets, and get significant improvement upon baselines, proving the practical value of our approach.


title: Nonlinear mode competition and symmetry-protected power oscillations in topological lasers

abstract: Topological photonics started out as a pursuit to engineer systems that mimic fermionic single-particle Hamiltonians with symmetry-protected modes, whose number can only change in spectral phase transitions such as band inversions. The paradigm of topological lasing, realized in three recent experiments, offers entirely new interpretations of these states, as they can be selectively amplified by distributed gain and loss. A key question is whether such topological mode selection persists when one accounts for the nonlinearities that stabilize such systems at their working point. Here we show that topological defect lasers can indeed stably operate in genuinely topological states. These comprise direct analogues of zero modes from the linear setting, as well as a novel class of states displaying symmetry-protected power oscillations, which appear in a spectral phase transition when the gain is increased. These effects show a remarkable practical resilience against imperfections, even if these break the underlying symmetries, and pave the way to harness the power of topological protection in nonlinear quantum devices.


title: Strong Pseudo Transitivity and Intersection Graphs

abstract: A directed graph $G=(V,E)$ is {\it strongly pseudo transitive} if there is a partition $\{A,E-A\}$ of $E$ so that graphs $G_1=(V,A)$ and $G_2=(V,E-A)$ are transitive, and additionally, if $ab\in A$ and $bc\in E $ implies that $ac\in E$. A strongly pseudo transitive graph $G=(V,E)$ is strongly pseudo transitive of the first type, if $ab\in A$ and $bc\in E$ implies $ac\in A$. An undirected graph is co-strongly pseudo transitive (co-strongly pseudo transitive of the first type) if its complement has an orientation which is strongly pseudo transitive (co-strongly pseudo transitive of the first type). Our purpose is show that the results in computational geometry \cite{CFP, Lu} and intersection graph theory \cite{Ga2, ES} can be unified and extended, using the notion of strong pseudo transitivity. As a consequence the general algorithmic framework in \cite{Sh} is applicable to solve the maximum independent set in $O(n^3)$ time in a variety of problems, thereby, avoiding case by case lengthily arguments for each problem. We show that the intersection graphs of axis parallel rectangles intersecting a diagonal line from bottom, and half segments are co-strongly pseudo transitive. In addition, we show that the class of the interval filament graphs is co-strongly transitive of the first type, and hence the class of polygon circle graphs which is contained in the class of interval filament graphs (but contains the classes of chordal graphs, circular arc, circle, and outer planar graphs), and the class of incomparability graphs are strongly transitive of the first type. For class of chordal graphs we give two different proofs, using two different characterizations, verifying that they are co-strongly transitive of the first type. We present some containment results.


title: Beam Spin Asymmetry in Electroproduction of Pseudoscalar or Scalar Meson Production off the Scalar Target

abstract: We discuss the electroproduction of pseudoscalar ($0^{-+}$) or scalar ($0^{++}$) meson production off the scalar target. The most general formulation of the differential cross section for the $0^{-+}$ or $0^{++}$ meson production process involves only one or two hadronic form factors, respectively, on a scalar target. The Rosenbluth-type separation of the differential cross section provides the explicit relation between the hadronic form factors and the different parts of the differential cross section in a completely model-independent manner. The absence of the beam spin asymmetry for the pseudoscalar meson production provides a benchmark for the experimental data analysis. The measurement of the beam spin asymmetry for the scalar meson production may also provide a unique opportunity not only to explore the imaginary part of the hadronic amplitude in the general formulation but also to examine the significance of the chiral-odd generalized parton distribution (GPD) contribution in the leading-twist GPD formulation.


title: A General Framework for Bandit Problems Beyond Cumulative Objectives

abstract: The stochastic multi-armed bandit (MAB) problem is a common model for sequential decision problems. In the standard setup, a decision maker has to choose at every instant between several competing arms, each of them provides a scalar random variable, referred to as a "reward." Nearly all research on this topic considers the total cumulative reward as the criterion of interest. This work focuses on other natural objectives that cannot be cast as a sum over rewards, but rather more involved functions of the reward stream. Unlike the case of cumulative criteria, in the problems we study here the oracle policy, that knows the problem parameters a priori and is used to "center" the regret, is not trivial. We provide a systematic approach to such problems, and derive general conditions under which the oracle policy is sufficiently tractable to facilitate the design of optimism-based (upper confidence bound) learning policies. These conditions elucidate an interesting interplay between the arm reward distributions and the performance metric. Our main findings are illustrated for several commonly used objectives such as conditional value-at-risk, mean-variance trade-offs, Sharpe-ratio, and more.


title: Wall-crossing for toric mutations

abstract: This note explains how to deduce the wall-crossing formula for toric mutations established by Pascaleff-Tonkonog from the perverse schober of the corresponding local Landau-Ginzburg model. Along the way, we develop a general framework to extract a wall-crossing formula from a perverse schober on the projective line with a single critical point.


title: Invertibility of adjacency matrices for random d-regular directed graphs

abstract: Let $d\geq 3$ be a fixed integer, and a prime number $p$ such that $\gcd(p,d)=1$. Let $A$ be the adjacency matrix of a random $d$-regular directed graph on $n$ vertices. We show that as a random matrix in ${\mathbb F}_p$, \begin{equation} {\mathbb P}(\text{$A$ is singular in ${\mathbb F}_p$})\leq \frac{1+{\mathrm{o}}(1)}{p-1}, \end{equation} as $n$ goes to infinity. As a consequence, as a random matrix in $\mathbb R$, \begin{equation} {\mathbb P}(\text{$A$ is singular in $\mathbb R$})={\mathrm{o}}(1) \end{equation} as $n$ goes to infinity. This answers an open problem by Frieze [12] and Vu [29,30], for random $d$-regular bipartite graphs. The proof combines a local central limit theorem and a large deviation estimate.


title: Towards a systematically improvable many-body description of antiferromagnetic iron oxide

abstract: We report variational and fixed-node diffusion quantum Monte Carlo (QMC) calculations of anti-ferromagnetic iron oxide (FeO) in the ground state B1 crystal structure. The goal of this study was a systematic investigation of the sensitivity of several ground state properties to a variety of QMC wave function generation techniques including advanced wave functions such as multi-determinant expansions and backflow transformations. We found that the predicted lattice distortion was largely controlled by the choice of single particle orbitals used to construct the wave function, rather than by subsequent wave function optimization techniques within QMC. However, the absolute magnetic moment was remarkably insensitive to the method of wave function construction. QMC estimates of total spin density indicate that in addition to strong electronic correlation of the Fe $3d$ states, charge transfer may be an important but challenging piece of physics to accurately capture within existing QMC methods. Finally, we highlight the need for advanced and systematically improvable many-body wave functions suitable for accurately describing challenging real systems.


title: Passive Static Equilibrium with Frictional Contacts and Application to Grasp Stability Analysis

abstract: This paper studies the problem of passive grasp stability under an external disturbance, that is, the ability of a grasp to resist a disturbance through passive responses at the contacts. To obtain physically consistent results, such a model must account for friction phenomena at each contact; the difficulty is that friction forces depend in non-linear fashion on contact behavior (stick or slip). We develop the first polynomial-time algorithm which either solves such complex equilibrium constraints for two-dimensional grasps, or otherwise concludes that no solution exists. To achieve this, we show that the number of possible `slip states' (where each contact is labeled as either sticking or slipping) that must be considered is polynomial (in fact quadratic) in the number of contacts, and not exponential as previously thought. Our algorithm captures passive response behaviors at each contact, while accounting for constraints on friction forces such as the maximum dissipation principle.


title: Nonperturbative Renormalization of Operators in Near-Conformal Systems Using Gradient Flows

abstract: We propose a continuous real space renormalization group transformation based on gradient flow, allowing for a numerical study of renormalization without the need for costly ensemble matching. We apply our technique in a pilot study of SU$(3)$ gauge theory with $N_f = 12$ fermions in the fundamental representation, finding the mass anomalous dimension to be $\gamma_m = 0.23(6)$, consistent with other perturbative and lattice estimates. We also present the first lattice calculation of the nucleon anomalous dimension in this theory, finding $\gamma_N = 0.05(5)$.


title: An empirical characterization of community structures in complex networks using a bivariate map of quality metrics

abstract: Community detection emerges as an important task in the discovery of network mesoscopic structures. However, the concept of a "good" community is very context-dependent and it is relatively complicated to deduce community characteristics using available community detection techniques. In reality, the existence of a gap between structural goodness quality metrics and expected topological patterns creates a confusion in evaluating community structures. In this paper, we introduce an empirical multivariate analysis of different structural goodness properties in order to characterize several detectable community topologies. Specifically, we show that a combination of two representative structural dimensions including community transitivity and hub dominance allows to distinguish different topologies such as star-based, clique-based, string-based and grid-based structures. Additionally, these classes of topology disclose structural proximities with those of graphs created by Erd\H{o}s-R\'{e}nyi, Watts-Strogatz and Barab\'{a}si-Albert generative models. We illustrate popular community topologies identified by different detection methods on a large dataset composing many network categories and associate their structures with the most related graph generative model. Interestingly, this conjunctive representation sheds light on fundamental differences between mesoscopic structures in various network categories including: communication, information, biological, technological, social, ecological, synthetic networks and more.


title: New And Surprising Ways to Be Mean. Adversarial NPCs with Coupled Empowerment Minimisation

abstract: Creating Non-Player Characters (NPCs) that can react robustly to unforeseen player behaviour or novel game content is difficult and time-consuming. This hinders the design of believable characters, and the inclusion of NPCs in games that rely heavily on procedural content generation. We have previously addressed this challenge by means of empowerment, a model of intrinsic motivation, and demonstrated how a coupled empowerment maximisation (CEM) policy can yield generic, companion-like behaviour. In this paper, we extend the CEM framework with a minimisation policy to give rise to adversarial behaviour. We conduct a qualitative, exploratory study in a dungeon-crawler game, demonstrating that CEM can exploit the affordances of different content facets in adaptive adversarial behaviour without modifications to the policy. Changes to the level design, underlying mechanics and our character's actions do not threaten our NPC's robustness, but yield new and surprising ways to be mean.


title: Revisiting black holes and wormholes under Weyl transformations

abstract: The behavior of black holes horizon and wormholes under the Weyl conformal transformation is investigated. First, a shorter, but more general, derivation of the Weyl transformation of the simple prescription for detecting horizons and wormholes given recently in the literature for spherically symmetric spacetimes is provided. The derivation allows for a simple and intuitive way to understand why and when horizons and wormholes might arise in the conformal frame even if they were absent in the original frame. Then, the conformal behavior of black holes horizon and wormholes in more general spacetimes, based on more "sophisticated" definitions, is provided. The study shows that black holes and wormholes might always arise in the new frame even if they were absent in the original frame. Moreover, it is shown that some of the definitions found in the literature might be transformed into one another under such transformations. Worked-out examples are given.


title: An analog magnon adder for all-magnonic neurons

abstract: Spin-waves are excellent data carriers with a perspective use in neuronal networks: Their lifetime gives the spin-wave system an intrinsic memory, they feature strong nonlinearity, and they can be guided and steered through extended magnonic networks. In this work, we present a magnon adder that integrates over incoming spin-wave pulses in an analog fashion. Such an adder is a linear prequel to a magnonic neuron, which would integrate over the incoming pulses until a certain nonlinearity is reached. In this work, the adder is realized by a resonator in combination with a parametric amplifier which is just compensating the resonator losses.


title: Aspects of NMR Reciprocity and Applications in Highly Conductive Media

abstract: In the context of NMR spectroscopy and MRI, the principle of reciprocity provides a convenient method for determining the reception sensitivity from the transmitted rf field pattern. The reciprocity principle for NMR was originally described by Hoult et al. [J. Magn. Reson. (1976),24, 71], and can be seen as being based on the broader Lorentz reciprocity principle, and similar theorems from antenna theory. One frequent application of the reciprocity principle is that for a single coil used for both transmission and detection, the transmit and receive fields can be assumed to be equal. This aspect is also where some of the conceptual difficulty of applying the theorem may be encountered. For example, the questions of whether one should use the complex conjugate field for detection, or whether one should apply the theorem in the rotating frame or the laboratory frame are often where considerable confusion may arise, and incorrect results may be derived. We attempt here to provide a helpful discussion of the application of the reciprocity principle in such a way as to clarify some of the confounding questions. In particular, we avoid the use of the `negatively rotating frame', which is frequently mentioned in this context, since we consider it to unnecessarily complicate the matter. In addition, we also discuss the implications of the theorem for magnetic resonance experiments on conducting samples, and metals, in particular.


title: The Araucaria Project: High-precision Cepheid astrophysics from the analysis of variables in double-lined eclipsing binaries

abstract: Based on new observations and improved modeling techniques, we have reanalyzed seven Cepheids in the Large Magellanic Cloud. Improved physical parameters have been determined for the exotic system OGLE LMC-CEP-1718 composed of two first-overtone Cepheids and a completely new model was obtained for the OGLE LMC-CEP-1812 classical Cepheid. This is now the shortest period Cepheid for which the projection factor is measured. The typical accuracy of our dynamical masses and radii determinations is 1%. The radii of the six classical Cepheids follow period--radius relations in the literature. Our very accurate physical parameter measurements allow us to calculate a purely empirical, tight period--mass--radius relation that agrees well with theoretical relations derived from non-canonical models. This empirical relation is a powerful tool to calculate accurate masses for single Cepheids for which precise radii can be obtained from Baade--Wesselink-type analyses. The mass of the type-II Cepheid $\kappa$ Pav, $0.56 \pm 0.08 M_\odot$, determined using this relation is in a very good agreement with theoretical predictions. We find large differences between the p-factor values derived for the Cepheids in our sample. Evidence is presented that a simple period--p-factor relation shows an intrinsic dispersion, hinting at the relevance of other parameters, such as the masses, radii, and radial velocity variation amplitudes. We also find evidence that the systematic blueshift exhibited by Cepheids, is primarily correlated with their gravity. The companion star of the Cepheid in the OGLE LMC-CEP-4506 system has a very similar temperature and luminosity, and is clearly located inside the Cepheid instability strip, yet it is not pulsating.


title: The shape of a memorised random walk

abstract: We view random walks as the paths of foraging animals, perhaps searching for food or avoiding predators while forming a mental map of their surroundings. The formation of such maps requires them to memorise the locations they have visited. We model memory using a kernel, proportional to the number of locations recalled as a function of the time since they were first observed. We give exact analytic expressions relating the elongation of the memorised walk to the structure of the memory kernel, and confirm these by simulation. We find that more slowly decaying memories lead to less elongated mental maps.


title: REORDER: Securing Dynamic-Priority Real-Time Systems Using Schedule Obfuscation

abstract: The deterministic (timing) behavior of real-time systems (RTS) can be used by adversaries - say, to launch side channel attacks or even destabilize the system by denying access to critical resources. We propose a protocol (named REORDER) to obfuscate this predictable timing behavior of RTS, especially ones designed using dynamic-priority scheduling algorithms (e.g., EDF). We also present a metric (named "schedule entropy") that measures the levels of obfuscation introduced into a given real-time system. The REORDER protocol was integrated into the standard Linux real-time scheduler and evaluated on a realistic embedded platform (Raspberry Pi) running the MiBench automotive benchmark workloads. We also demonstrate how designers of RTS can increase the security of their systems and also quantitatively measure the impact (both in terms of security and performance) of using this protocol.


title: Large spin-mixing conductance in highly Bi-doped Cu thin films

abstract: Spin Hall effect provides an efficient tool for the conversion of a charge current into a spin current, opening the possibility of producing pure spin currents in non-magnetic materials for the next generation of spintronics devices. In this sense, giant Spin Hall Effect has been recently reported in Cu doped with 0.5 % Bi grown by sputtering and larger values are expected for larger Bi doping, according to first principles calculations. In this work we demonstrate the possibility of doping Cu with up to 10 % of Bi atoms without evidences of Bi surface segregation or cluster formation, as studied by different microscopic and spectroscopic techniques. In addition, YIG/BiCu structures have been grown, showing a spin mixing conductance larger that the one shown by similar Pt/YIG structures. These results reflects the potentiality of these new materials in spintronics devices.


title: Quantum gyro-electric effect: Photon spin-1 quantization in continuum topological bosonic phases

abstract: Topological phases of matter arise in distinct fermionic and bosonic flavors. The fundamental differences between them are encapsulated in their rotational symmetries - the spin. Although spin quantization is routinely encountered in fermionic topological edge states, analogous quantization for bosons has proven elusive. To this end, we develop the complete electromagnetic continuum theory characterizing 2+1D topological bosons, taking into account their intrinsic spin and orbital angular momentum degrees of freedom. We demonstrate that spatiotemporal dispersion (momentum and frequency dependence of linear response) captures the matter-mediated interactions between bosons and is a necessary ingredient for topological phases. We prove that the bulk topology of these 2+1D phases is manifested in transverse spin-1 quantization of the photon. From this insight, we predict two unique bosonic phases - one with even parity $C=\pm 2$ and one with odd $C=\pm 1$. To understand the even parity phase $C=\pm 2$, we introduce an exactly solvable model utilizing non-local optical Hall conductivity and reveal a single gapless photon at the edge. This unidirectional photon is spin-1 helically quantized, immune to backscattering, defects, and exists at the boundary of the $C=\pm 2$ bosonic phase and any interface - even vacuum. The contrasting phenomena of transverse quantization in the bulk, but longitudinal (helical) quantization on the edge is addressed as the quantum gyro-electric effect (QGEE). We also validate our bosonic Maxwell theory by direct comparison with the supersymmetric Dirac theory of fermions. To accelerate the discovery of such bosonic phases, we suggest two new probes of topological matter with broken time-reversal symmetry: momentum-resolved electron energy loss spectroscopy and cold atom near-field measurement of non-local optical Hall conductivity.


title: Learning to track on-the-fly using a particle filter with annealed- weighted QPSO modeled after a singular Dirac delta potential

abstract: This paper proposes an evolutionary Particle Filter with a memory guided proposal step size update and an improved, fully-connected Quantum-behaved Particle Swarm Optimization (QPSO) resampling scheme for visual tracking applications. The proposal update step uses importance weights proportional to velocities encountered in recent memory to limit the swarm movement within probable regions of interest. The QPSO resampling scheme uses a fitness weighted mean best update to bias the swarm towards the fittest section of particles while also employing a simulated annealing operator to avoid subpar fine tune during latter course of iterations. By moving particles closer to high likelihood landscapes of the posterior distribution using such constructs, the sample impoverishment problem that plagues the Particle Filter is mitigated to a great extent. Experimental results using benchmark sequences imply that the proposed method outperforms competitive candidate trackers such as the Particle Filter and the traditional Particle Swarm Optimization based Particle Filter on a suite of tracker performance indices.


title: Comprehensive numerical modeling of filamentary RRAM devices including voltage ramp-rate and cycle-to-cycle variations

abstract: The equilibrium ON and OFF states of resistive random access memory (RRAM) are due to formation and destruction of a conducting filament. The laws of thermodynamics dictate that these states correspond to the minimum of free energy. Here, we develop a numerical model that, through the minimization of free energy at a given voltage, determines the filament parameters and thus the electric current. Overall, it simulates the current-voltage (I-V) characteristics of RRAM. The model describes mutual transformations of RRAM states through SET (ON to OFF) and RESET (OFF to ON) processes. From the modeling perspectives, these states and processes constitute four programming modules constructed here in COMSOL Multiphysics software tackling the electrodynamic and heat transfer equations and yielding RRAM energy and I-V. Our modeling uniquely reproduces the observed I-V varying with voltage ramp-rates. That is achieved by accounting for the ramp-rate dependent activation energy of conduction. The underlying mechanism is due to the deformation interaction caused by the double well atomic potentials universally present in amorphous materials and having exponentially broad distribution of relaxation times. As another unique feature, our modeling reproduces the observed cycle-to-cycle variations of RRAM parameters attributed to the lack of self-averaging in small ensembles of double well potentials and electronic states in geometrically small (nano-size) RRAM structures.


title: Pseudofinite H-structures and groups definable in supersimple H-structures

abstract: In this paper we explore some properties of H-structures. We describe a construction of H-structures based on one-dimensional asymptotic classes which preserves pseudo-finiteness. That is, the H-structures we construct are ultraproducts of finite structures. We also prove that under the assumption that the base theory is supersimple of SU-rank one, there are no new definable groups in H-structures.


title: Evolution of the transmission phase through a Coulomb-blockaded Majorana wire

abstract: We present a study of the transmission of electrons through a semiconductor quantum wire with strong spin-orbit coupling in proximity to an s-wave superconductor, which is Coulomb-blockaded. Such a system supports Majorana zero modes in the presence of an external magnetic field. Without superconductivity, phase lapses are expected to occur in the transmission phase, and we find that they disappear when a topological superconducting phase is stabilized. We express tunneling through the nanowire with the help of effective matrix elements, which depend on both the fermion parity of the wire and the overlap with Bogoliubov-de-Gennes wave functions. Using a modified scattering matrix formalism, that allows for including electron-electron interactions, we study the transmission phase in different regimes.


title: Mining large-scale human mobility data for long-term crime prediction

abstract: Traditional crime prediction models based on census data are limited, as they fail to capture the complexity and dynamics of human activity. With the rise of ubiquitous computing, there is the opportunity to improve such models with data that make for better proxies of human presence in cities. In this paper, we leverage large human mobility data to craft an extensive set of features for crime prediction, as informed by theories in criminology and urban studies. We employ averaging and boosting ensemble techniques from machine learning, to investigate their power in predicting yearly counts for different types of crimes occurring in New York City at census tract level. Our study shows that spatial and spatio-temporal features derived from Foursquare venues and checkins, subway rides, and taxi rides, improve the baseline models relying on census and POI data. The proposed models achieve absolute R^2 metrics of up to 65% (on a geographical out-of-sample test set) and up to 89% (on a temporal out-of-sample test set). This proves that, next to the residential population of an area, the ambient population there is strongly predictive of the area's crime levels. We deep-dive into the main crime categories, and find that the predictive gain of the human dynamics features varies across crime types: such features bring the biggest boost in case of grand larcenies, whereas assaults are already well predicted by the census features. Furthermore, we identify and discuss top predictive features for the main crime categories. These results offer valuable insights for those responsible for urban policy or law enforcement.


title: On estimation and inference in latent structure random graphs

abstract: We define a latent structure model (LSM) random graph as a random dot product graph (RDPG) in which the latent position distribution incorporates both probabilistic and geometric constraints, delineated by a family of underlying distributions on some fixed Euclidean space, and a structural support submanifold from which the latent positions for the graph are drawn. For a one-dimensional latent structure model with known structural support, we show how spectral estimates of the latent positions of an RDPG can be used for efficient estimation of the paramaters of the LSM. We describe how to estimate or learn the structural support in cases where it is unknown, with an illustrative focus on graphs with latent positions along the Hardy-Weinberg curve. Finally, we use the latent structure model formulation to test bilateral homology in the Drosophila connectome.


title: Small MAD families whose Isbell-Mr\'owka spaces are pseudocompact

abstract: Given a countable transitive model $M$ for ZFC+CH, we prove that one can produce a maximal almost disjoint family in $M$ whose Vietoris Hyperspace of its Isbell-Mr\'owka space is pseudocompact on every Cohen extension of $M$. We also show that a classical example of $\omega_1$-sized maximal almost disjoint family obtained by a forcing iteration of length $\omega_1$ in a model of non CH is such that the Vietoris Hyperspace of its Isbell-Mr\'owka space is pseudocompact.


title: A Bayesian Penalized Hidden Markov Model for Ant Interactions

abstract: Interactions between social animals provide insights into the exchange and flow of nutrients, disease, and social contacts. We consider a chamber level analysis of trophallaxis interactions between carpenter ants (\textit{Camponotus pennsylvanicus}) over 4 hours of second-by-second observations. The data show clear switches between fast and slow modes of trophallaxis. However, fitting a standard hidden Markov model (HMM) results in an estimated hidden state process that is overfit to this high resolution data, as the state process fluctuates an order of magnitude more quickly than is biologically reasonable. We propose a novel approach for penalized estimation of HMMs through a Bayesian ridge prior on the state transition rates while also incorporating biologically motivated covariates. This penalty induces smoothing, limiting the rate of state switching that combines with appropriate covariates within the colony to ensure more biologically feasible results. We develop a Markov chain Monte Carlo algorithm to perform Bayesian inference based on discretized observations of the contact network.


title: Note on the number of divisors of reducible quadratic polynomials

abstract: In a recent paper, Lapkova uses a Tauberian theorem to derive the asymptotic formula for the divisor sum $\sum_{n \leq x} d( n (n+v))$ where $v$ is a fixed integer and $d(n)$ denotes the number of divisors of $n$. We reprove her result by following a suggestion of Hooley, namely investigating the relationship between this sum and the well-known sum $\sum_{n \leq x} d( n ) d (n+v)$. As such, we are able to furnish additional terms in the asymptotic formula.


title: On the Soundness of Coroutines with Snapshots

abstract: Coroutines are a general control flow construct that can eliminate control flow fragmentation inherent in event-driven programs, and are still missing in many popular languages. Coroutines with snapshots are a first-class, type-safe, stackful coroutine model, which unifies many variants of suspendable computing, and is sufficiently general to express iterators, single-assignment variables, async-await, actors, event streams, backtracking, symmetric coroutines and continuations. In this paper, we develop a formal model called $\lambda_{\rightsquigarrow}$ (lambda-squiggly) that captures the essence of type-safe, stackful, delimited coroutines with snapshots. We prove the standard progress and preservation safety properties. Finally, we show a formal transformation from the $\lambda_{\rightsquigarrow}$ calculus to the simply-typed lambda calculus with references.


title: Spin-Orbit induced phase-shift in Bi$_{2}$Se$_{3}$ Josephson junctions

abstract: The transmission of Cooper pairs between two weakly coupled superconductors produces a superfluid current and a phase difference; the celebrated Josephson effect. Because of time-reversal and parity symmetries, there is no Josephson current without a phase difference between two superconductors. Reciprocally, when those two symmetries are broken, an anomalous supercurrent can exist in the absence of phase bias or, equivalently, an anomalous phase shift $\varphi_0$ can exist in the absence of a superfluid current. We report on the observation of an anomalous phase shift $\varphi_0$ in hybrid Josephson junctions fabricated with the topological insulator Bi$_2$Se$_3$ submitted to an in-plane magnetic field. This anomalous phase shift $\varphi_0$ is observed directly through measurements of the current-phase relationship in a Josephson interferometer. This result provides a direct measurement of the spin-orbit coupling strength and open new possibilities for phase-controlled Josephson devices made from materials with strong spin-orbit coupling.


title: A Solvable Deformation of Quantum Mechanics

abstract: The conventional Hamiltonian $H= p^2+ V_N(x)$, where the potential $V_N(x)$ is a polynomial of degree $N$, has been studied intensively since the birth of quantum mechanics. In some cases, its spectrum can be determined by combining the WKB method with resummation techniques. In this paper we point out that the deformed Hamiltonian $H=2 \cosh(p)+ V_N(x)$ is exactly solvable for any potential: a conjectural exact quantization condition, involving well-defined functions, can be written down in closed form, and determines the spectrum of bound states and resonances. In particular, no resummation techniques are needed. This Hamiltonian is obtained by quantizing the Seiberg-Witten curve of $\mathcal{N}=2$ Yang-Mills theory, and the exact quantization condition follows from the correspondence between spectral theory and topological strings, after taking a suitable four-dimensional limit. In this formulation, conventional quantum mechanics emerges in a scaling limit near the Argyres-Douglas superconformal point in moduli space. Although our deformed version of quantum mechanics is in many respects similar to the conventional version, it also displays new phenomena, like spontaneous parity symmetry breaking.


title: Mixed Ax-Schanuel for the universal abelian varieties and some applications

abstract: In this paper we prove the mixed Ax-Schanuel theorem for the universal abelian varieties (more generally any mixed Shimura variety of Kuga type), and give some simple applications. In particular we present an application to studying the generic rank of the Betti map.


title: Coupling charge-exchange vibrations to nucleons in a relativistic framework: effect on Gamow-Teller transitions and beta-decay half-lives

abstract: The nuclear response theory for isospin-transfer modes in the relativistic particle-vibration coupling framework is extended to include coupling of single nucleons to isospin-flip (charge-exchange) phonons, in addition to the usual neutral vibrations. This new coupling introduces dynamical pion and rho-meson exchange, beyond the Hartree-Fock approximation, up to infinite order. We investigate the impact of this new mechanism on the Gamow-Teller response of a few doubly-magic neutron-rich nuclei, namely $^{48}$Ca, $^{78}$Ni, $^{132}$Sn and $^{208}$Pb. It is found that the coupling to isospin-flip vibrations can have a non negligible impact on the strength distribution and quenching of the giant resonance, globally improving the agreement with the experimental data. The corresponding beta-decay half-lives of $^{78}$Ni and $^{132}$Sn are also calculated, and found to be decreased by a factor $\sim 2$ by the inclusion of the new phonons.


title: A diffuse interface framework for modelling the evolution of multi-cell aggregates as a soft packing problem driven by the growth and division of cells

abstract: We present a model for cell growth, division and packing under soft constraints that arise from the deformability of the cells as well as of a membrane that encloses them. Our treatment falls within the framework of diffuse interface methods, under which each cell is represented by a scalar phase field and the zero level set of the phase field represents the cell membrane. One crucial element in the treatment is the definition of a free energy density function that penalizes cell overlap, thus giving rise to a simple model of cell-cell contact. In order to properly represent cell packing and the associated free energy, we include a simplified representation of the anisotropic mechanical response of the underlying cytoskeleton and cell membrane through appropriate penalization of the cell shape change. Numerical examples are presented to demonstrate the evolution of multi-cell clusters, and the total free energy of the clusters as a consequence of growth, division and packing.


title: FlowNet3D: Learning Scene Flow in 3D Point Clouds

abstract: Many applications in robotics and human-computer interaction can benefit from understanding 3D motion of points in a dynamic environment, widely noted as scene flow. While most previous methods focus on stereo and RGB-D images as input, few try to estimate scene flow directly from point clouds. In this work, we propose a novel deep neural network named $FlowNet3D$ that learns scene flow from point clouds in an end-to-end fashion. Our network simultaneously learns deep hierarchical features of point clouds and flow embeddings that represent point motions, supported by two newly proposed learning layers for point sets. We evaluate the network on both challenging synthetic data from FlyingThings3D and real Lidar scans from KITTI. Trained on synthetic data only, our network successfully generalizes to real scans, outperforming various baselines and showing competitive results to the prior art. We also demonstrate two applications of our scene flow output (scan registration and motion segmentation) to show its potential wide use cases.


title: A fast algorithm for maximum likelihood estimation of mixture proportions using sequential quadratic programming

abstract: Maximum likelihood estimation of mixture proportions has a long history, and continues to play an important role in modern statistics, including in development of nonparametric empirical Bayes methods. Maximum likelihood of mixture proportions has traditionally been solved using the expectation maximization (EM) algorithm, but recent work by Koenker & Mizera shows that modern convex optimization techniques -- in particular, interior point methods -- are substantially faster and more accurate than EM. Here, we develop a new solution based on sequential quadratic programming (SQP). It is substantially faster than the interior point method, and just as accurate. Our approach combines several ideas: first, it solves a reformulation of the original problem; second, it uses an SQP approach to make the best use of the expensive gradient and Hessian computations; third, the SQP iterations are implemented using an active set method to exploit the sparse nature of the quadratic subproblems; fourth, it uses accurate low-rank approximations for more efficient gradient and Hessian computations. We illustrate the benefits of our approach in experiments on synthetic data sets as well as a large genetic association data set. In large data sets (n = 1,000,000 observations, m = 1,000 mixture components), our implementation achieves at least 100-fold reduction in runtime compared with a state-of-the-art interior point solver. Our methods are implemented in Julia, and in an R package available on CRAN (see https://CRAN.R-project.org/package=mixsqp).


title: CFCM: Segmentation via Coarse to Fine Context Memory

abstract: Recent neural-network-based architectures for image segmentation make extensive usage of feature forwarding mechanisms to integrate information from multiple scales. Although yielding good results, even deeper architectures and alternative methods for feature fusion at different resolutions have been scarcely investigated for medical applications. In this work we propose to implement segmentation via an encoder-decoder architecture which differs from any other previously published method since (i) it employs a very deep architecture based on residual learning and (ii) combines features via a convolutional Long Short Term Memory (LSTM), instead of concatenation or summation. The intuition is that the memory mechanism implemented by LSTMs can better integrate features from different scales through a coarse-to-fine strategy; hence the name Coarse-to-Fine Context Memory (CFCM). We demonstrate the remarkable advantages of this approach on two datasets: the Montgomery county lung segmentation dataset, and the EndoVis 2015 challenge dataset for surgical instrument segmentation.


title: The 2016 outburst of PSR J1119-6127: cooling & a spin-down dominated glitch

abstract: We report on the aftermath of a magnetar outburst from the young, high-magnetic-field radio pulsar PSR J1119-6127 that occurred on 2016 July 27. We present the results of a monitoring campaign using the Neil Gehrels Swift X-ray Telescope, NuSTAR, and XMM-Newton. After reaching a peak luminosity of ~300 times the quiescent luminosity, the pulsar's X-ray flux declined by factor of ~50 on a time scale of several months. The X-ray spectra are well described by a blackbody and a hard power-law tail. After an initial rapid decline during the first day of the outburst, we observe the blackbody temperature rising from kT = 0.9 keV to 1.05 keV during the first two weeks of the outburst, before cooling to 0.9 keV. During this time, the blackbody radius decreases monotonically by a factor of ~4 over a span of nearly 200 days. We also report a heretofore unseen highly pulsed hard X-ray emission component, which fades on a similar timescale to the soft X-ray flux, as predicted by models of relaxation of magnetospheric current twists. The previously reported spin-up glitch which accompanied this outburst was followed by a period of enhanced and erratic torque, leading to a net spin-down of $\sim3.5\times10^{-4}$ Hz, a factor of ~24 over-recovery. We suggest that this and other radiatively loud magnetar-type glitch recoveries are dominated by magnetospheric processes, in contrast to conventional radio pulsar glitch recoveries which are dominated by internal physics.


title: Magnetic oscillations in silicene

abstract: In this work the magnetic oscillations (MO) in pristine silicene at $T=0$ K are studied. Considering a constant electron density we obtain analytical expressions for the ground state internal energy and magnetization, under a perpendicular electric and magnetic field, taking in consideration the Zeeman effect. It is found that the MO are sawtooth-like, depending on the change in the last occupied energy level. This leads us to a classification of the MO peaks in terms of the Landau level (LL), valley or spin changes. Using this classification we analyze the MO for different values of the electric field $E_{z}$. When $E_{z}=0$, the energy levels have a valley degeneracy and the MO peaks occur only whenever the last energy level changes its LL and/or spin. When $E_{z}\neq0$, the valley degeneracy is broken and new MO peaks appear, associated with the valley change in the last energy level. By analyzing the MO peaks amplitude it is possible to extract information about the Fermi velocity and the spin-orbit interaction strength. Finally we analyze the MO frequencies, which can also be associated with the change of LL, valley or spin in the last energy level.


title: Phase-type distributions in population genetics

abstract: Probability modelling for DNA sequence evolution is well established and provides a rich framework for understanding genetic variation between samples of individuals from one or more populations. We show that both classical and more recent models for coalescence (with or without recombination) can be described in terms of the so-called phase-type theory, where complicated and tedious calculations are circumvented by the use of matrices. The application of phase-type theory consists of describing the stochastic model as a Markov model by appropriately setting up a state space and calculating the corresponding intensity and reward matrices. Formulae of interest are then expressed in terms of these aforementioned matrices. We illustrate this by a few examples calculating the mean, variance and even higher order moments of the site frequency spectrum in the multiple merger coalescent models, and by analysing the mean and variance for the number of segregating sites for multiple samples in the two-locus ancestral recombination graph. We believe that phase-type theory has great potential as a tool for analysing probability models in population genetics. The compact matrix notation is useful for clarification of current models, in particular their formal manipulation (calculation), but also for further development or extensions.


title: Sensitivity of the PICO-500 Bubble Chamber to Supernova Neutrinos Through Coherent Nuclear Elastic Scattering

abstract: Ton-scale direct dark matter search experiments should be sensitive to neutrino-induced recoil events from either $^8$B solar neutrinos or the brief but intense flux from a core collapse supernova in the Milky Way. These low-threshold detectors are sensitive to the very low recoil energies, of order 10 keV, deposited via coherent elastic scatters between supernova neutrinos and target nuclei. Large superheated fluid detectors like PICO-500, a bubble chamber to be initially filled with an active target of 1 t of C$_3$F$_8$, should see multiple-bubble events from CE$\nu$NS if the detector is live during a neutrino burst from a supernova at a distance up to 10 kpc. This paper discusses conditions under which bubble chambers could be used as an independent measurement in the event of a supernova similar to SN 1987A, with particular sensitivity to the currently less-constrained heavy-lepton $\nu_x$ channel.

categories: astro-ph.HE astro-ph.IM astro-ph.SR hep-ex physics.ins-det


title: Precision ephemerides for gravitational-wave searches -- III. Revised system parameters of Sco X-1

abstract: Neutron stars in low-mass X-ray binaries are considered promising candidate sources of continuous gravitational-waves. These neutron stars are typically rotating many hundreds of times a second. The process of accretion can potentially generate and support non-axisymmetric distortions to the compact object, resulting in persistent emission of gravitational-waves. We present a study of existing optical spectroscopic data for Sco X-1, a prime target for continuous gravitational-wave searches, with the aim of providing revised constraints on key orbital parameters required for a directed search with advanced-LIGO data. From a circular orbit fit to an improved radial velocity curve of the Bowen emission components, we derived an updated orbital period and ephemeris. Centre of symmetry measurements from the Bowen Doppler tomogram yield a centre of the disc component of 90 km $\mathrm{s^{-1}}$, which we interpret as a revised upper limit to the projected orbital velocity of the NS $K_1$. By implementing Monte Carlo binary parameter calculations, and imposing new limits on $K_1$ and the rotational broadening, we obtained a complete set of dynamical system parameter constraints including a new range for $K_1$ of 40--90 km $\mathrm{s^{-1}}$. Finally, we discussed the implications of the updated orbital parameters for future continuous-waves searches.


title: Permanent Magnet Synchronous Motors are Globally Asymptotically Stabilizable with PI Current Control

abstract: This note shows that the industry standard desired equilibrium for permanent magnet synchronous motors (i.e., maximum torque per Ampere) can be globally asymptotically stabilized with a PI control around the current errors, provided some viscous friction (possibly small) is present in the rotor dynamics and the proportional gain of the PI is suitably chosen. Instrumental to establish this surprising result is the proof that the map from voltages to currents of the incremental model of the motor satisfies some passivity properties. The analysis relies on basic Lyapunov theory making the result available to a wide audience.


title: Wolf 1465: Not a Bright Dwarf Carbon Star

abstract: Wolf 1465 has been suggested as a candidate for a bright dwarf carbon star. It is not.


title: Towards Quantum Integrated Information Theory

abstract: Integrated Information Theory (IIT) has emerged as one of the leading research lines in computational neuroscience to provide a mechanistic and mathematically well-defined description of the neural correlates of consciousness. Integrated Information ($\Phi$) quantifies how much the integrated cause/effect structure of the global neural network fails to be accounted for by any partitioned version of it. The holistic IIT approach is in principle applicable to any information-processing dynamical network regardless of its interpretation in the context of consciousness. In this paper we take the first steps towards a formulation of a general and consistent version of IIT for interacting networks of quantum systems. A variety of different phases, from the dis-integrated ($\Phi=0$) to the holistic one (extensive $\log\Phi$), can be identified and their cross-overs studied.


title: A scalable matrix-free spectral element approach for unsteady PDE constrained optimization using PETSc/TAO

abstract: We provide a new approach for the efficient matrix-free application of the transpose of the Jacobian for the spectral element method for the adjoint based solution of partial differential equation (PDE) constrained optimization. This results in optimizations of nonlinear PDEs using explicit integrators where the integration of the adjoint problem is not more expensive than the forward simulation. Solving PDE constrained optimization problems entails combining expertise from multiple areas, including simulation, computation of derivatives, and optimization. The Portable, Extensible Toolkit for Scientific computation (PETSc) together with its companion package, the Toolkit for Advanced Optimization (TAO), is an integrated numerical software library that contains an algorithmic/software stack for solving linear systems, nonlinear systems, ordinary differential equations, differential algebraic equations, and large-scale optimization problems and, as such, is an ideal tool for performing PDE-constrained optimization. This paper describes an efficient approach in which the software stack provided by PETSc/TAO can be used for large-scale nonlinear time-dependent problems. While time integration can involve a range of high-order methods, both implicit and explicit. The PDE-constrained optimization algorithm used is gradient-based and seamlessly integrated with the simulation of the physical problem.