title: A Euclidean Signature Semi-Classical Program

abstract: In this article we discuss our ongoing program to extend the scope of certain, well-developed microlocal methods for the asymptotic solution of Schr\"{o}dinger's equation (for suitable `nonlinear oscillatory' quantum mechanical systems) to the treatment of several physically significant, interacting quantum field theories. Our main focus is on applying these `Euclidean-signature semi-classical' methods to self-interacting (real) scalar fields of renormalizable type in 2, 3 and 4 spacetime dimensions and to Yang-Mills fields in 3 and 4 spacetime dimensions. A central argument in favor of our program is that the asymptotic methods for Schr\"{o}dinger operators developed in the microlocal literature are far superior, for the quantum mechanical systems to which they naturally apply, to the conventional WKB methods of the physics literature and that these methods can be modified, by techniques drawn from the calculus of variations and the analysis of elliptic boundary value problems, to apply to certain (bosonic) quantum field theories. Unlike conventional (Rayleigh/ Schr\"{o}dinger) perturbation theory these methods allow one to avoid the artificial decomposition of an interacting system into an approximating `unperturbed' system and its perturbation and instead to keep the nonlinearities (and, if present gauge invariances) of an interacting system intact at every level of the analysis.


title: Generalization of the small-world effect on a model approaching the Erd\H{o}s-R\'enyi random graph

abstract: The famous Watts-Strogatz (WS) small-world network model does not approach the Erd\H{o}s-R\'enyi (ER) random graph model in the limit of total randomization which can lead to confusion and complicates certain analyses. In this paper we discuss a simple alternative which was first introduced by Song and Wang, where instead of rewiring, edges are drawn between pairs of nodes with a distance-based connection probability. We show that this model is simpler to analyze, approaches the true ER random graph model in the completely randomized limit, and demonstrate that the WS model and the alternative model may yield different quantitative results using the example of a random walk temporal observable. An efficient sampling algorithm for the alternative model is proposed. Analytic results regarding the degree distribution, degree variance, number of two-stars per node, number of triangles per node, clustering coefficient, and random walk mixing time are presented. Subsequently, the small-world effect is illustrated by showing that the clustering coefficient decreases much slower than an upper bound on the message delivery time with increasing long-range connection probability which generalizes the small-world effect from informed searches to random search strategies. Due to its accessibility for analytic evaluations, we propose that this modified model should be used as an alternative reference model for studying the influence of small-world topologies on dynamic systems as well as a simple model to introduce numerous topics when teaching network science.


title: Description of shape coexistence in $^{96}$Zr based on the collective quadrupole Bohr Hamiltonian

abstract: Experimental data on $^{96}$Zr indicate coexisting spherical and deformed structures with small mixing amplitudes. We investigate the properties of the low-lying collective states of $^{96}$Zr based on the collective quadrupole Bohr Hamiltonian. The $\beta$-dependent collective potential having two minima -- spherical and deformed, is fixed so to describe experimental data in the best way.Good agreement with the experimental data on the excitation energies, $B(E2)$ and $B(M1)$ reduced transition probabilities is obtained. It is shown that the low-energy structure of $^{96}$Zr can be reproduced in a satisfactory way in the geometrical model with a potential function supporting shape coexistence. However, the excitation energy of the $2^+_2$ state can be reproduced only if the rotation inertia coefficient is taken five times smaller then the vibrational one in the region of the deformed well. It is shown also that shell effects are important for the description of the $B(M1;2^+_2 \rightarrow 2^+_1)$ value. An indication on the influence of the pairing vibrational mode on the $\rho^2 (0^+_2 \rightarrow 0^+_1)$ value is obtained.


title: The secondary transit of the hot Jupiter WASP-121b at 2 $\mu$m

abstract: Ground-based observations of the secondary transit in the 2MASS K band are presented for the hot Jupiter WASP-121b. These are the first occultation observations of an extrasolar planet carried out with an instrument attached to a 1m-class telescope (SMARTS' 1.3 m). We find a highly significant transit depth of (0.228 +/- 0.023)%. Together with the Hubble Space Telescope near infrared emission spectrum, current data support more involved atmosphere models with species producing emission/absorption features, rather than simple smooth black body emission. Analysis of the time difference between the primary and secondary transits and the durations of these events yield an eccentricity of e=0.0207 +/- 0.0153, which is consistent with the earlier estimates of low/zero eccentricity, but with a smaller error. Together with the existing K-band data on other systems, WASP-121b lends further support to the lack of efficient heat transport between the day and night sides for nearly all Hot Jupiters.


title: Public Health and access to medicine. Pharmaceutical industry's role

abstract: Every year, 10 million people die from lack of access to treatment for curable diseases, specially in developing countries. Meanwhile, legal but unsafe drugs cause 130 thousand deaths per year. How can this be happening in 21st Century? What role does the pharmaceutical industry play in this tragedy? In this research, WHO reports are analyzed and primary information gathered so as to answer this questions.


title: Stochastic analysis of emergence of evolutionary cyclic behavior in population dynamics with transfer

abstract: Horizontal gene transfer consists in exchanging genetic materials between microorganisms during their lives. This is a major mechanism of bacterial evolution and is believed to be of main importance in antibiotics resistance. We consider a stochastic model for the evolution of a discrete population structured by a trait taking finitely many values, with density-dependent competition. Traits are vertically inherited unless a mutation occurs, and can also be horizontally transferred by unilateral conjugation with frequency dependent rate. Our goal is to analyze the trade-off between natural evolution to higher birth rates on one side, and transfer which drives the population towards lower birth rates on the other side. Simulations show that evolutionary outcomes include evolutionary suicide or cyclic re-emergence of small populations with well-adapted traits. We focus on a parameter scaling where individual mutations are rare but the global mutation rate tends to infinity. This implies that negligible sub-populations may have a strong contribution to evolution. Our main result quantifies the asymptotic dynamics of subpopulation sizes on a logarithmic scale. We characterize the possible evolutionary outcomes with explicit criteria on the model parameters. An important ingredient for the proofs lies in comparisons of the stochastic population process with linear or logistic birth-death processes with immigration. For the latter processes, we derive several results of independent interest.


title: A mathematical model of contact inhibition of locomotion: coupling contractility and focal adhesions

abstract: Cell migration is often accompanied by collisions with other cells, which can lead to cessation of movement, repolarization, and migration away from the contact site - a process termed contact inhibition of locomotion (CIL). During CIL, the coupling between actomyosin contractilityand cell-substrate adhesions is modified. However, mathematical models describing stochastic cell migration and collision outcomes as a result of the coupling remain elusive. Here, we extend our previously developed stochastic model of single cell migration to include CIL. Our simulation results explain, in terms of the modified contractility and adhesion dynamics, several experimentally observed findings regarding CIL. These include response modulation in the presence of an external cue and alterations of group migration in the absence of CIL. Together with our previous findings, our work is able to explain a wide range of observations about single and collective cell migration.


title: Mitigation of Cosmic Ray Effect on Microwave Kinetic Inductance Detector Arrays

abstract: For space observatories, the glitches caused by high energy phonons created by the interaction of cosmic ray particles with the detector substrate lead to dead time during observation. Mitigating the impact of cosmic rays is therefore an important requirement for detectors to be used in future space missions. In order to investigate possible solutions, we carry out a systematic study by testing four large arrays of Microwave Kinetic Inductance Detectors (MKIDs), each consisting of $\sim$960 pixels and fabricated on monolithic 55 mm $\times$ 55 mm $\times$ 0.35 mm Si substrates. We compare the response to cosmic ray interactions in our laboratory for different detector arrays: A standard array with only the MKID array as reference; an array with a low $T_c$ superconducting film as phonon absorber on the opposite side of the substrate; and arrays with MKIDs on membranes. The idea is that the low $T_c$ layer down-converts the phonon energy to values below the pair breaking threshold of the MKIDs, and the membranes isolate the sensitive part of the MKIDs from phonons created in the substrate. We find that the dead time can be reduced up to a factor of 40 when compared to the reference array. Simulations show that the dead time can be reduced to below 1 % for the tested detector arrays when operated in a spacecraft in an L2 or a similar far-Earth orbit. The technique described here is also applicable and important for large superconducting qubit arrays for future quantum computers.


title: A focusing and defocusing semi-discrete complex short pulse equation and its varioius soliton solutions

abstract: In this paper, we are concerned with a a semi-discrete complex short pulse (CSP) equation of both focusing and defocusing types, which can be viewed as an analogue to the Ablowitz-Ladik (AL) lattice in the ultra-short pulse regime. By using a generalized Darboux transformation method, various solutions to this newly integrable semi-discrete equation are studied with both zero and nonzero boundary conditions. To be specific, for the focusing CSP equation, the multi-bright solution (zero boundary condition), multi-breather and high-order rogue wave solutions (nonzero boudanry conditions) are derived, while for the defocusing CSP equation with nonzero boundary condition, the multi-dark soliton solution is constructed. We further show that, in the continuous limit, all the solutions obtained converge to the ones for its original CSP equation (see Physica D, 327 13-29 and Phys. Rev. E 93 052227)


title: Non-ergodic extended states in the SYK model

abstract: We analytically study spectral correlations and many body wave functions of an SYK-model deformed by a one body contribution to the Hamiltonian. Our main result is the identification of a wide range of intermediate coupling strengths where the spectral statistics is of Wigner-Dyson type, while wave functions are non-uniformly distributed over Fock space and show fractal behavior. The structure of the theory suggests that such manifestations of non-ergodic extendedness may be a prevalent phenomenon in strongly interacting chaotic quantum systems.


title: Energy Crowdsourcing and Peer-to-Peer Energy Trading in Blockchain-Enabled Smart Grids

abstract: The power grid is rapidly transforming, and while recent grid innovations increased the utilization of advanced control methods, the next-generation grid demands technologies that enable the integration of distributed energy resources (DERs)---and consumers that both seamlessly buy and sell electricity. This paper develops an optimization model and blockchain-based architecture to manage the operation of crowdsourced energy systems (CES), with peer-to-peer (P2P) energy trading transactions. An operational model of CESs in distribution networks is presented considering various types of energy trading transactions and crowdsourcees. Then, a two-phase operation algorithm is presented: Phase I focuses on the day-ahead scheduling of generation and controllable DERs, whereas Phase II is developed for hour-ahead or real-time operation of distribution networks. The developed approach supports seamless P2P energy trading between individual prosumers and/or the utility. The presented operational model can also be used to operate islanded microgrids. The CES framework and the operation algorithm are then prototyped through an efficient blockchain implementation, namely the IBM Hyperledger Fabric. This implementation allows the system operator to manage the network users to seamlessly trade energy. Case studies and prototype illustration are provided.


title: Modeling tax distribution in metropolitan regions with PolicySpace

abstract: Brazilian executive body has consistently vetoed legislative initiatives easing creation and emancipation of municipalities. The literature lists evidence of the negative results of municipal fragmentation, especially so for metropolitan regions. In order to provide evidences for the argument of metropolitan union, this paper quantifies the quality of life of metropolitan citizens in the face of four alternative rules of distribution of municipal tax collection. Methodologically, a validated agent-based spatial model is simulated. On top of that, econometric models are tested using real exogenous variables and simulated data. Results suggest two central conclusions. First, the progressiveness of the Municipal Participation Fund and its relevance to a better quality of life in metropolitan municipalities is confirmed. Second, municipal financial merging would improve citizens' quality of life, compared to the status quo for 23 Brazilian metropolises. Further, the paper presents quantitative evidence that allows comparing alternative tax distributions for each of the 40 simulated metropolises, identifying more efficient forms of fiscal distribution and contributing to the literature and to contemporary parliamentary debate.


title: Forced magnetic reconnection and plasmoid coalescence: I - MHD Simulations

abstract: Forced magnetic reconnection, a reconnection event triggered by external perturbation, should be ubiquitous in the solar corona. Energy released during such cases can be much greater than that which was introduced by the perturbation. It is unclear how the properties of the external perturbation and the initial current sheet affect the reconnection region properties, and thereby the reconnection dynamics and energy release profile. We investigate the effect of the form of the external perturbation and initial current sheet on the evolution of the reconnection region and the energy release process. Chiefly we explore the non-linear interactions between multiple, simultaneous perturbations, which represent more realistic scenarios. Future work will use these results in test particle simulations to investigate particle acceleration over multiple reconnection events. Simulations are performed using Lare2d, a 2.5D Lagrangian-remap solver for the visco-resistive MHD equations. The model of forced reconnection is extended to include superpositions of sinusoidal driving disturbances, including localised Gaussian perturbations. A transient perturbation is applied to the boundaries of a region containing a force-free current sheet. The simulation domain is sufficiently wide to allow multiple magnetic islands to form and coalesce. Island coalescence contributes significantly to energy release and involves rapid reconnection. Long wavelength modes in perturbations dominate the evolution, without the presence of which reconnection is either slow, as in the case of short wavelength modes, or the initial current sheet remains stable, as in the case of noise perturbations. Multiple perturbations combine in a highly non-linear manner: reconnection is typically faster than when either disturbance is applied individually, with multiple low-energy events contributing to the same total energy release.


title: Fair Algorithms for Clustering

abstract: We study the problem of finding low-cost Fair Clusterings in data where each data point may belong to many protected groups. Our work significantly generalizes the seminal work of Chierichetti et.al. (NIPS 2017) as follows. - We allow the user to specify the parameters that define fair representation. More precisely, these parameters define the maximum over- and minimum under-representation of any group in any cluster. - Our clustering algorithm works on any $\ell_p$-norm objective (e.g. $k$-means, $k$-median, and $k$-center). Indeed, our algorithm transforms any vanilla clustering solution into a fair one incurring only a slight loss in quality. - Our algorithm also allows individuals to lie in multiple protected groups. In other words, we do not need the protected groups to partition the data and we can maintain fairness across different groups simultaneously. Our experiments show that on established data sets, our algorithm performs much better in practice than what our theoretical results suggest.


title: On Completion of \textit{C}*-algebra-valued metric spaces

abstract: The concept of a $ C $*-algebra-valued metric space was introduced in 2014. It is a generalization of a metric space by replacing the set of real numbers by a $ C $*-algebra. In this paper, we show that $ C $*-algebra-valued metric spaces are cone metric spaces in some point of view which is useful to extend results of the cone case to $ C $*-algebra-valued metric spaces. Then the completion theorem of $ C $*-algebra-valued metric spaces is obtained. Moreover, the completion theorem of $ C $*-algebra-valued normed spaces is verified and the connection with Hilbert $ C $*-modules, generalized inner product spaces, is also provided.


title: High entropy oxides: An emerging prospect for magnetic rare earth - transition metal perovskites

abstract: It has been shown that oxide ceramics containing multiple transition and/or rare-earth elements in equimolar ratios have a strong tendency to crystallize in simple single phase structures, stabilized by the high configurational entropy. In analogy to the metallic alloy systems, these oxides are denoted high entropy oxides (HEOs). The HEO concept allows to access hitherto uncharted areas in the multi-element phase diagram. Among the already realized structures there is the highly complex class of rare earth - transition element perovskites. This fascinating class of materials generated by applying the innovative concept of high entropy stabilization provides a new and manyfold research space with promise of discoveries of unprecedented properties and phenomena. The present study provides a first investigation of the magnetic properties of selected compounds of this novel class of materials. Comprehensive studies by DC and AC magnetometry are combined with element specific spectroscopy in order to understand the interplay between magnetic exchange and the high degree of chemical disorder in the systems. We observe a predominant antiferromagnetic behavior in the single phase materials, combined with a small ferromagnetic contribution possibly stemming from small ferromagnetic clusters or configurations in the antiferromagnetic matrix. In the long term perspective it is proposed to screen the properties of this family of compounds with high throughput methods, including combined experimental and theoretical approaches.


title: Thermal bridging of graphene nanosheets via covalent molecular junctions. a Non-Equilibrium Green Functions Density Functional Tight-Binding study

abstract: Despite the uniquely high thermal conductivity of graphene is well known, the exploitation of graphene into thermally conductive nanomaterials and devices is limited by the inefficiency of thermal contacts between the individual nanosheets. A fascinating yet experimentally challenging route to enhance thermal conductance at contacts between graphene nanosheets is through molecular junctions, allowing covalently connecting nanosheets otherwise interacting only via weak Van der Waals forces. Beside the bare existence of covalent connections, the choice of molecular structures to be used as thermal junctions should be guided by their vibrational properties, in terms of phonon transfer through the molecular junction. In this paper, density functional tight-binding combined with Green functions formalism was applied for the calculation of thermal conductance and phonon spectra of several different aliphatic and aromatic molecular junctions between graphene nanosheets. Effect of molecular junction length, conformation, and aromaticity were studied in detail and correlated with phonon tunnelling spectra. The theoretical insight provided by this work can guide future experimental studies to select suitable molecular junctions in order to enhance the thermal transport by suppressing the interfacial thermal resistances, particularly attractive for various systems, including graphene nanopapers and graphene polymer nanocomposites, as well as related devices. In a broader view, the possibility to design molecular junctions to control phonon transport currently appears as an efficient way to produce phononic devices and controlling heat management in nanostructures.


title: Ito's conjecture and the coset construction for ${\mathcal W}^k(\mathfrak{sl}(3|2))$

abstract: Many ${\mathcal W}$-(super)algebras which are defined by the generalized Drinfeld-Sokolov reduction are also known or expected to have coset realizations. For example, it was conjectured by Ito that the principal ${\mathcal W}$-superalgebra ${\mathcal W}^k(\mathfrak{sl}(n+1|n))$ is isomorphic to the coset of $V^{l+1}(\mathfrak{gl}_n)$ inside $V^{l}(\mathfrak{sl}_{n+1}) \otimes {\mathcal E}(n)$ for generic values of $l$. Here ${\mathcal E}(n)$ denotes the rank $n$ $bc$-system, which carries an action of $V^1(\mathfrak{gl}_n)$, and $k$ and $l$ are related by $(k + 1) (l + n + 1) = 1$. This conjecture is known in the case $n=1$, which is somewhat degenerate, and we shall prove it in the first nontrivial case $n=2$. As a consequence, we show that the simple quotient ${\mathcal W}_k(\mathfrak{sl}(3|2))$ is lisse and rational for all positive integers $l>1$. These are new examples of rational ${\mathcal W}$-superalgebras.


title: Monotone Least Squares and Isotonic Quantiles

abstract: We consider bivariate observations $(X_1,Y_1), \ldots, (X_n,Y_n)$ such that, conditional on the $X_i$, the $Y_i$ are independent random variables with distribution functions $F_{X_i}$, where $(F_x)_x$ is an unknown family of distribution functions. Under the sole assumption that $x \mapsto F_x$ is isotonic with respect to stochastic order, one can estimate $(F_x)_x$ in two ways: (i) For any fixed $y$ one estimates the antitonic function $x \mapsto F_x(y)$ via nonparametric monotone least squares, replacing the responses $Y_i$ with the indicators $1_{[Y_i \le y]}$. (ii) For any fixed $\beta \in (0,1)$ one estimates the isotonic quantile function $x \mapsto F_x^{-1}(\beta)$ via a nonparametric version of regression quantiles. We show that these two approaches are closely related, with (i) being more flexible than (ii). Then, under mild regularity conditions, we establish rates of convergence for the resulting estimators $\hat{F}_x(y)$ and $\hat{F}_x^{-1}(\beta)$, uniformly over $(x,y)$ and $(x,\beta)$ in certain rectangles as well as uniformly in $y$ or $\beta$ for a fixed $x$.


title: Service Rate Region of Content Access from Erasure Coded Storage

abstract: We consider storage systems in which $K$ files are stored over $N$ nodes. A node may be systematic for a particular file in the sense that access to it gives access to the file. Alternatively, a node may be coded, meaning that it gives access to a particular file only when combined with other nodes (which may be coded or systematic). Requests for file $f_k$ arrive at rate $\lambda_k$, and we are interested in the rate that can be served by a particular system. In this paper, we determine the set of request arrival rates for the a $3$-file coded storage system. We also provide an algorithm to maximize the rate of requests served for file $K$ given $\lambda_1,\dots, \lambda_{K-1}$ in a general $K$-file case.


title: Can the homogeneity scale be used as a standard ruler?

abstract: The Universe on comoving scales larger than 100 Mpc/h is assumed to be statistically homogeneous, with the transition scale where we have $1\%$ deviation from homogeneity being known as the homogeneity scale $R_H$. The latter was recently proposed in Ntelis et al. (2018) [arXiv:1810.09362] to be a standard ruler. In this paper we perform a comparison between the baryon acoustic oscillations and the homogeneity scales $R_H$, assuming a spatially flat universe, linear cosmological perturbation theory and the cosmological constant cold dark matter $\Lambda$CDM model. By direct theoretical calculations, we demonstrate that the answer to the question of whether the homogeneity scale can be used as a standard ruler, is clearly negative due to the non-monotonicity of the homogeneity scale $R_H$ with the matter density parameter $\Omega_{m0}$ of the $\Lambda$CDM model. Furthermore, we also consider the effect of redshift space distortions and we find that they do not break the degeneracies but instead only change the value of $R_H$ by a few percent. Finally, we also confirm our findings with the help of an N-Body simulation.


title: The Buzzard Flock: Dark Energy Survey Synthetic Sky Catalogs

abstract: We present a suite of 18 synthetic sky catalogs designed to support science analysis of galaxies in the Dark Energy Survey Year 1 (DES Y1) data. For each catalog, we use a computationally efficient empirical approach, ADDGALS, to embed galaxies within light-cone outputs of three dark matter simulations that resolve halos with masses above ~5x10^12 h^-1 m_sun at z <= 0.32 and 10^13 h^-1 m_sun at z~2. The embedding method is tuned to match the observed evolution of galaxy counts at different luminosities as well as the spatial clustering of the galaxy population. Galaxies are lensed by matter along the line of sight --- including magnification, shear, and multiple images --- using CALCLENS, an algorithm that calculates shear with 0.42 arcmin resolution at galaxy positions in the full catalog. The catalogs presented here, each with the same LCDM cosmology (denoted Buzzard), contain on average 820 million galaxies over an area of 1120 square degrees with positions, magnitudes, shapes, photometric errors, and photometric redshift estimates. We show that the weak-lensing shear catalog, redMaGiC galaxy catalogs and redMaPPer cluster catalogs provide plausible realizations of the same catalogs in the DES Y1 data by comparing their magnitude, color and redshift distributions, angular clustering, and mass-observable relations, making them useful for testing analyses that use these samples. We make public the galaxy samples appropriate for the DES Y1 data, as well as the data vectors used for cosmology analyses on these simulations.


title: Contamination Attacks and Mitigation in Multi-Party Machine Learning

abstract: Machine learning is data hungry; the more data a model has access to in training, the more likely it is to perform well at inference time. Distinct parties may want to combine their local data to gain the benefits of a model trained on a large corpus of data. We consider such a case: parties get access to the model trained on their joint data but do not see each others individual datasets. We show that one needs to be careful when using this multi-party model since a potentially malicious party can taint the model by providing contaminated data. We then show how adversarial training can defend against such attacks by preventing the model from learning trends specific to individual parties data, thereby also guaranteeing party-level membership privacy.


title: Formulation of the laws of conservation and non-conservation

abstract: For each substance-like quantity, a theorem about its conservation or non-conservation can be formulated. For the electric charge e.g. it reads: Electric charge can neither be created nor destroyed. Such a statement is short and easy to understand. For some quantities, however, the proposition about conservation or non-conservation is usually formulated in an unnecessarily complicated way. Sometimes the formulation is not generally valid; in other cases only a consequence of the conservation or non-conservation is pronounced. A clear and unified formulation could improve the comprehensibility and simplify teaching.


title: GILT: Generating Images from Long Text

abstract: Creating an image reflecting the content of a long text is a complex process that requires a sense of creativity. For example, creating a book cover or a movie poster based on their summary or a food image based on its recipe. In this paper we present the new task of generating images from long text that does not describe the visual content of the image directly. For this, we build a system for generating high-resolution 256 $\times$ 256 images of food conditioned on their recipes. The relation between the recipe text (without its title) to the visual content of the image is vague, and the textual structure of recipes is complex, consisting of two sections (ingredients and instructions) both containing multiple sentences. We used the recipe1M dataset to train and evaluate our model that is based on a the StackGAN-v2 architecture.


title: A High Resolution PDE Approach to Quadrilateral Mesh Generation

abstract: We describe a high order technique to generate quadrilateral decompositions and meshes for complex two dimensional domains using spectral elements in a field guided procedure. Inspired by cross field methods, we never actually compute crosses. Instead, we compute a high order accurate guiding field using a continuous Galerkin (CG) or discontinuous Galerkin (DG) spectral element method to solve a Laplace equation for each of the field variables using the open source code Nektar++. The spectral method provides spectral convergence and sub-element resolution of the fields. The DG approximation allows meshing of corners that are not multiples of $\pi/2$ in a discretization consistent manner, when needed. The high order field can then be exploited to accurately find irregular nodes, and can be accurately integrated using a high order separatrix integration method to avoid features like limit cycles. The result is a mesh with naturally curved quadrilateral elements that do not need to be curved a posteriori to eliminate invalid elements. The mesh generation procedure is implemented in the open source mesh generation program NekMesh.


title: Using ZDDs in the Mapping of Quantum Circuits

abstract: A critical step in quantum compilation is the transformation of a technology-independent quantum circuit into a technology-dependent form for a targeted device. In addition to mapping quantum gates into the supported gate set, it is necessary to map pseudo qubits in the technology-independent circuit into physical qubits of the technology-dependent circuit such that coupling constraints among qubits acting in multiple-qubit gates are satisfied. It is usually not possible to find such a mapping without adding SWAP gates into the circuit. To cope with the technical limitations of NISQ-era quantum devices, it is advantageous to find a mapping that requires as few additional gates as possible. The large search space of possible mappings makes this task a difficult combinatorial optimization problem. In this work, we demonstrate how zero-suppressed decision diagrams (ZDDs) can be used for typical implementation tasks in quantum mapping algorithms. We show how to maximally partition a quantum circuit into blocks of adjacent gates, and if adjacent gates within a circuit do not share common mapping permutations, we attempt to combine them using parallelized SWAP operations represented in a ZDD. Boundaries for the partitions are formed where adjacent gates are unable to be combined. Within each partition block, ZDDs represent all possible mappings of pseudo qubits to physical qubits.


title: The temporal requirements of directly observing self-gravitating spiral waves in protoplanetary discs with ALMA

abstract: We investigate how the detectability of signatures of self-gravity in a protoplanetary disc depends on its temporal evolution. We run a one-dimensional model for secular timescales to follow the disc mass as a function of time. We then combine this with three-dimensional global hydrodynamics simulations that employ a hybrid radiative transfer method to approximate realistic heating and cooling. We simulate ALMA continuum observations of these systems, and find that structures induced by the gravitational instability (GI) are readily detectable when $q=M_\mathrm{disc}/M_*\gtrsim 0.25$ and $R_\mathrm{outer}\lesssim 100$ au. The high accretion rate generated by gravito-turbulence in such a massive disc drains its mass to below the detection threshold in $\sim10^4$ years, or approximately 1 % of the typical disc lifetime. Therefore, discs with spiral arms detected in ALMA dust observations, if generated by self-gravity, must either be still receiving infall to maintain a high $q$ value, or have just emerged from their natal envelope. Detection of substructure in systems with lower $q$ is possible, but would require a specialist integration with the most extended configuration over several days. This disfavours the possibility of GI-caused spiral structure in systems with $q<0.25$ being detected in relatively short integration times, such as those found in the DSHARP ALMA survey (Andrews et al. 2018; Huang et al. 2018). We find no temporal dependence of detectability on dynamical timescales


title: Further results for a subclass of univalent functions related with differential equation

abstract: Let $\Omega$ denote the class of functions $f$ analytic in the open unit disc $\Delta$, normalized by the condition $f(0)=f'(0)-1=0$ and satisfying the inequality \begin{equation*} \left|zf'(z)-f(z)\right|<\frac{1}{2}\quad(z\in\Delta). \end{equation*} The class $\Omega$ was introduced recently by Peng and Zhong (Acta Math Sci {\bf37B(1)}:69--78, 2017). Also let $\mathcal{U}$ denote the class of functions $f$ analytic and normalized in $\Delta$ and satisfying the condition \begin{equation*} \left|\left(\frac{z}{f(z)}\right)^2f'(z)-1\right|<1\quad(z\in\Delta). \end{equation*} In this article, we obtain some further results for the class $\Omega$ including, an extremal function and more examples of $\Omega$, inclusion relation between $\Omega$ and $\mathcal{U}$, the radius of starlikeness, convexity and close--to--convexity and sufficient condition for function $f$ to be in $\Omega$. Furthermore, along with the settlement of the coefficient problem and the Fekete--Szeg\"{o} problem for the elements of $\Omega$, the Toeplitz matrices for $\Omega$ are also discussed in this article.


title: Regularity of stable solutions to quasilinear elliptic equations on Riemannian models

abstract: We investigate the regularity of semi-stable, radially symmetric, and decreasing solutions for a class of quasilinear reaction-diffusion equations in the inhomogeneous context of Riemannian manifolds. We prove uniform boundedness, Lebesgue and Sobolev estimates for this class of solutions for equations involving the p-Laplace Beltrami operator and locally Lipschitz non-linearity. We emphasize that our results do not depend on the boundary conditions and the specific form of the non-linearities and metric. Moreover, as an application, we establish regularity of the extremal solutions for equations involving the p-Laplace Beltrami operator with zero Dirichlet boundary conditions.


title: An Open Source, Versatile, Affordable Waves in Ice Instrument for Scientific Measurements in the Polar Regions

abstract: Sea ice is a major feature of the polar environments. Recent changes in the climate and extent of the sea ice, together with increased economic activity and research interest in these regions, are driving factors for new measurements of sea ice dynamics. Waves in ice are important as they participate in the coupling between the open ocean and the ice-covered regions. Measurements are challenging to perform due to remoteness and harsh environmental conditions. While progress has been made in observing wave propagtion in sea ice using remote methods, these are still relatively new measurements and would benefit from more in situ data for validation. In this article, we present an open source instrument that was developed for performing such measurements. The versatile design includes an ultra-low power unit, a microcontroller-based logger, a small microcomputer for on-board data processing, and an Iridium modem for satellite communications. Virtually any sensor can be used with this design. In the present case, we use an Inertial Motion Unit to record wave motion. High quality results were obtained, which opens new possibilities for in situ measurements in the polar regions. Our instrument can be easily customized to fit many in situ measurement tasks, and we hope that our work will provide a framework for future developments of a variety of such open source instruments.


title: Morphological Networks for Image De-raining

abstract: Mathematical morphological methods have successfully been applied to filter out (emphasize or remove) different structures of an image. However, it is argued that these methods could be suitable for the task only if the type and order of the filter(s) as well as the shape and size of operator kernel are designed properly. Thus the existing filtering operators are problem (instance) specific and are designed by the domain experts. In this work we propose a morphological network that emulates classical morphological filtering consisting of a series of erosion and dilation operators with trainable structuring elements. We evaluate the proposed network for image de-raining task where the SSIM and mean absolute error (MAE) loss corresponding to predicted and ground-truth clean image is back-propagated through the network to train the structuring elements. We observe that a single morphological network can de-rain an image with any arbitrary shaped rain-droplets and achieves similar performance with the contemporary CNNs for this task with a fraction of trainable parameters (network size). The proposed morphological network(MorphoN) is not designed specifically for de-raining and can readily be applied to similar filtering / noise cleaning tasks. The source code can be found here https://github.com/ranjanZ/2D-Morphological-Network


title: Forecasting Granular Audience Size for Online Advertising

abstract: Orchestration of campaigns for online display advertising requires marketers to forecast audience size at the granularity of specific attributes of web traffic, characterized by the categorical nature of all attributes (e.g. {US, Chrome, Mobile}). With each attribute taking many values, the very large attribute combination set makes estimating audience size for any specific attribute combination challenging. We modify Eclat, a frequent itemset mining (FIM) algorithm, to accommodate categorical variables. For consequent frequent and infrequent itemsets, we then provide forecasts using time series analysis with conditional probabilities to aid approximation. An extensive simulation, based on typical characteristics of audience data, is built to stress test our modified-FIM approach. In two real datasets, comparison with baselines including neural network models, shows that our method lowers computation time of FIM for categorical data. On hold out samples we show that the proposed forecasting method outperforms these baselines.


title: Interpretable CNNs for Object Classification

abstract: This paper proposes a generic method to learn interpretable convolutional filters in a deep convolutional neural network (CNN) for object classification, where each interpretable filter encodes features of a specific object part. Our method does not require additional annotations of object parts or textures for supervision. Instead, we use the same training data as traditional CNNs. Our method automatically assigns each interpretable filter in a high conv-layer with an object part of a certain category during the learning process. Such explicit knowledge representations in conv-layers of CNN help people clarify the logic encoded in the CNN, i.e., answering what patterns the CNN extracts from an input image and uses for prediction. We have tested our method using different benchmark CNNs with various structures to demonstrate the broad applicability of our method. Experiments have shown that our interpretable filters are much more semantically meaningful than traditional filters.


title: Resource Allocation in One-dimensional Distributed Service Networks

abstract: We consider assignment policies that allocate resources to users, where both resources and users are located on a one-dimensional line. First, we consider unidirectional assignment policies that allocate resources only to users located to their left. We propose the Move to Right (MTR) policy, which scans from left to right assigning nearest rightmost available resource to a user, and contrast it to the Unidirectional Gale-Shapley (UGS) matching policy. While both these policies are optimal among all unidirectional policies, we show that they are equivalent with respect to the expected distance traveled by a request (request distance), although MTR is fairer. Moreover, we show that when user and resource locations are modeled by statistical point processes, and resources are allowed to satisfy more than one user, the spatial system under unidirectional policies can be mapped into bulk service queuing systems, thus allowing the application of a plethora of queuing theory results that yield closed form expressions. As we consider a case where different resources can satisfy different numbers of users, we also generate new results for bulk service queues. We also consider bidirectional policies where there are no directional restrictions on resource allocation and develop an algorithm for computing the optimal assignment which is more efficient than known algorithms in the literature when there are more resources than users. Finally, numerical evaluation of performance of unidirectional and bidirectional allocation schemes yields design guidelines beneficial for resource placement.


title: SNRA: A Spintronic Neuromorphic Reconfigurable Array for In-Circuit Training and Evaluation of Deep Belief Networks

abstract: In this paper, a spintronic neuromorphic reconfigurable Array (SNRA) is developed to fuse together power-efficient probabilistic and in-field programmable deterministic computing during both training and evaluation phases of restricted Boltzmann machines (RBMs). First, probabilistic spin logic devices are used to develop an RBM realization which is adapted to construct deep belief networks (DBNs) having one to three hidden layers of size 10 to 800 neurons each. Second, we design a hardware implementation for the contrastive divergence (CD) algorithm using a four-state finite state machine capable of unsupervised training in N+3 clocks where N denotes the number of neurons in each RBM. The functionality of our proposed CD hardware implementation is validated using ModelSim simulations. We synthesize the developed Verilog HDL implementation of our proposed test/train control circuitry for various DBN topologies where the maximal RBM dimensions yield resource utilization ranging from 51 to 2,421 lookup tables (LUTs). Next, we leverage spin Hall effect (SHE)-magnetic tunnel junction (MTJ) based non-volatile LUTs circuits as an alternative for static random access memory (SRAM)-based LUTs storing the deterministic logic configuration to form a reconfigurable fabric. Finally, we compare the performance of our proposed SNRA with SRAM-based configurable fabrics focusing on the area and power consumption induced by the LUTs used to implement both CD and evaluation modes. The results obtained indicate more than 80% reduction in combined dynamic and static power dissipation, while achieving at least 50% reduction in device count.


title: Variational analysis of low-lying states in supersymmetric Yang-Mills theory

abstract: We have calculated the masses of bound states numerically in N = 1 supersymmetric Yang-Mills theory with gauge group SU(2). Using the suitably optimised variational method with an operator basis consisting of smeared Wilson loops and mesonic operators, we are able to obtain the masses of the ground states and first excited states in the scalar, pseudoscalar and spin-1/2 sectors. Extrapolated to the continuum limit, the corresponding particles appear to be approximately mass degenerate and to fit into the predicted chiral supermultiplets. The extended operator basis including both glueball-like and mesonic operators leads to improved results compared to earlier studies, and moreover allows us to investigate the mixing content of the physical states, which we compare to predictions in the literature.


title: Reflecting on Inaccessible J\'onsson Cardinals

abstract: Given an inaccessible J\'onsson cardinal $\lambda$, a sequence of results due to Shelah from Cardinal Arithmetic and Sh413 tell us that $\lambda$ must be at least $\lambda\times\omega$-Mahlo. We may then ask ourselves whether we can improve this bound, and unfortunately not much progress has been made on this question since Sh413. Part of the problem may be that the proofs involved are rather complicated and make use of somewhat ad hoc machinery. In order to remedy this, we present a survey and new proofs of these results regarding the Mahlo degree of an inaccessible J\'onsson cardinal with the hope of making this material accessible.


title: Cosmology with kSZ: breaking the optical depth degeneracy with Fast Radio Bursts

abstract: The small-scale cosmic microwave background (CMB) is dominated by anisotropies from the kinematic Sunyaev-Zeldovich (kSZ) effect, and upcoming experiments will measure it very precisely, but the optical depth degeneracy limits the cosmological information that can be extracted. At the same time, fast radio bursts (FRBs) are an exciting new frontier for astrophysics, but their usefulness as cosmological probes is currently unclear. We show that FRBs are uniquely suited for breaking the kSZ optical depth degeneracy. This opens up new possibilities for constraining cosmology with the kSZ effect, and new cosmological applications for FRBs.


title: Dynamic tail inference with log-Laplace volatility

abstract: We propose a family of models that enable predictive estimation of time-varying extreme event probabilities in heavy-tailed and nonlinearly dependent time series. The models are a white noise process with conditionally log-Laplace stochastic volatility. In contrast to other, similar stochastic volatility formalisms, this process has analytic expressions for its conditional probabilistic structure that enable straightforward estimation of dynamically changing extreme event probabilities. The process and volatility are conditionally Pareto-tailed, with tail exponent given by the reciprocal of the log-volatility's mean absolute innovation. This formalism can accommodate a wide variety of nonlinear dependence, as well as conditional power law-tail behavior ranging from weakly non-Gaussian to Cauchy-like tails. We provide a computationally straightforward estimation procedure that uses an asymptotic approximation of the process' dynamic large deviation probabilities. We demonstrate the estimator's utility with a simulation study. We then show the method's predictive capabilities on a simulated nonlinear time series where the volatility is driven by the chaotic Lorenz system. Lastly we provide an empirical application, which shows that this simple modeling method can be effectively used for dynamic and predictive tail inference in financial time series.

categories: stat.ME econ.EM math.ST q-fin.RM q-fin.ST stat.TH


title: Using deceased-donor kidneys to initiate chains of living donor kidney paired donations: algorithms and experimentation

abstract: We design a flexible algorithm that exploits deceased donor kidneys to initiate chains of living donor kidney paired donations, combining deceased and living donor allocation mechanisms to improve the quantity and quality of kidney transplants. The advantages of this approach have been measured using retrospective data on the pool of donor/recipient incompatible and desensitized pairs at the Padua University Hospital, the largest center for living donor kidney transplants in Italy. The experiments show a remarkable improvement on the number of patients with incompatible donor who could be transplanted, a decrease in the number of desensitization procedures, and an increase in the number of UT patients (that is, patients unlikely to be transplanted for immunological reasons) in the waiting list who could receive an organ.


title: Stationary waves with prescribed $L^2$-norm for the planar Schr\"odinger-Poisson system

abstract: The paper deals with the existence of standing wave solutions for the Schr\"odinger-Poisson system with prescribed mass in dimension $N=2$. This leads to investigate the existence of normalized solutions for an integro-differential equation involving a logarithmic convolution potential, namely $$ \left \{ \begin{aligned} - \Delta u & + \lambda u + \gamma \Bigl(\log {| \cdot |} * |u|^2 \Bigr) u =a |u|^{p-2} u \qquad \text{in $\mathbb R^2$,} \\ &\int_{\mathbb R^2} |u|^2 dx = c \end{aligned} \right. $$ where $c>0$ is a given real number. Under different assumptions on $\gamma \in \mathbb R$, $a \in \mathbb R$, $p>2$, we prove several existence and multiplicity results. Here $\lambda \in \mathbb R $ appears as a Lagrange parameter and is part of the unknowns. With respect to the related higher dimensional cases, the presence of the logarithmic kernel, which is unbounded from above and below, makes the structure of the solution set much richer and it forces the implementation of new ideas to catch the normalized solutions.


title: Rank one density property for a class of $M$-bases

abstract: In the early 1990s the works of Larson, Wogen and Argyros, Lambrou, Longstaff disclosed an example of a strong $M$-basis that did not admit a linear summation method. We study a class of $M$-bases $\mathfrak{F}=\{f_n\}_{n=1}^\infty$ in the Hilbert space that generalizes the Larson--Wogen system. We determine the conditions under which $\mathfrak{F}$ admits a linear summation method. In order to do that we employ some of the graph theory techniques.


title: Using think-aloud interviews to characterize model-based reasoning in electronics for a laboratory course assessment

abstract: Models of physical systems are used to explain and predict experimental results and observations. The Modeling Framework for Experimental Physics describes the process by which physicists revise their models to account for the newly acquired observations, or change their apparatus to better represent their models when they encounter discrepancies between actual and expected behavior of a system. While modeling is a nationally recognized learning outcome for undergraduate physics lab courses, no assessments of students' model-based reasoning exist for upper-division labs. As part of a larger effort to create two assessments of students' modeling abilities, we used the Modeling Framework to develop and code think-aloud problem-solving activities centered on investigating an inverting amplifier circuit. This study is the second phase of a multiphase assessment instrument development process. Here, we focus on characterizing the range of modeling pathways students employ while interpreting the output signal of a circuit functioning far outside its recommended operation range. We end by discussing four outcomes of this work: (1) Students engaged in all modeling subtasks, and they spent the most time making measurements, making comparisons, and enacting revisions; (2) Each subtask occurred in close temporal proximity to all over subtasks; (3) Sometimes, students propose causes that do not follow logically from observed discrepancies; (4) Similarly, students often rely on their experiential knowledge and enact revisions that do not follow logically from articulated proposed causes.


title: On the value distribution of the Riemann zeta-function and the Euler gamma-function

abstract: We prove some uniqueness results for the Riemann zeta-function and the Euler gamma-function by virtue of shared values using the value distribution theory.


title: Richer and Deeper Supervision Network for Salient Object Detection

abstract: Recent Salient Object Detection (SOD) systems are mostly based on Convolutional Neural Networks (CNNs). Specifically, Deeply Supervised Saliency (DSS) system has shown it is very useful to add short connections to the network and supervising on the side output. In this work, we propose a new SOD system which aims at designing a more efficient and effective way to pass back global information. Richer and Deeper Supervision (RDS) is applied to better combine features from each side output without demanding much extra computational space. Meanwhile, the backbone network used for SOD is normally pre-trained on the object classification dataset, ImageNet. But the pre-trained model has been trained on cropped images in order to only focus on distinguishing features within the region of the object. But the ignored background information is also significant in the task of SOD. We try to solve this problem by introducing the training data designed for object detection. A coarse global information is learned based on an entire image with its bounding box before training on the SOD dataset. The large-scale of object images can slightly improve the performance of SOD. Our experiment shows the proposed RDS network achieves the state-of-the-art results on five public SOD datasets.


title: Designing Data Protection for GDPR Compliance into IoT Healthcare Systems

abstract: In this paper, we investigate the implications of the General Data Privacy Regulation (GDPR) on the design of an IoT healthcare system. On 25th May 2018, the GDPR has become mandatory within the European Union and hence also for all suppliers of IT products. Infringements on the regulation are now fined with penalties of up 20 Million EUR or 4\% of the annual turnover of a company whichever is higher. This is a clear motivation for system designers to guarantee compliance to the GDPR. We propose a data labeling model to support access control for privacy-critical patient data together with the Fusion/UML process to design GDPR compliant system. We illustrate this design process on the case study of IoT based monitoring of Alzheimer's patients that we work on in the CHIST-ERA project SUCCESS.


title: Adaptive Activity Monitoring with Uncertainty Quantification in Switching Gaussian Process Models

abstract: Emerging wearable sensors have enabled the unprecedented ability to continuously monitor human activities for healthcare purposes. However, with so many ambient sensors collecting different measurements, it becomes important not only to maintain good monitoring accuracy, but also low power consumption to ensure sustainable monitoring. This power-efficient sensing scheme can be achieved by deciding which group of sensors to use at a given time, requiring an accurate characterization of the trade-off between sensor energy usage and the uncertainty in ignoring certain sensor signals while monitoring. To address this challenge in the context of activity monitoring, we have designed an adaptive activity monitoring framework. We first propose a switching Gaussian process to model the observed sensor signals emitting from the underlying activity states. To efficiently compute the Gaussian process model likelihood and quantify the context prediction uncertainty, we propose a block circulant embedding technique and use Fast Fourier Transforms (FFT) for inference. By computing the Bayesian loss function tailored to switching Gaussian processes, an adaptive monitoring procedure is developed to select features from available sensors that optimize the trade-off between sensor power consumption and the prediction performance quantified by state prediction entropy. We demonstrate the effectiveness of our framework on the popular benchmark of UCI Human Activity Recognition using Smartphones.


title: Quantum stability in open string theory with broken supersymmetry

abstract: We consider the 1-loop effective potential in type I string theory compactified on a torus, with supersymmetry broken by the Scherk-Schwarz mechanism. At fixed supersymmetry breaking scale M, and up to exponentially suppressed terms, we show that the potential admits local minima of arbitrary sign, in dimension d<=5. While the open string Wilson lines are massive, the closed string moduli are flat directions. In a T-dual picture, the relevant backgrounds involve isolated 1/2-branes, whose positions are frozen on orientifold planes, thus decreasing the rank of the gauge group, and introducing massless fermions in fundamental representations.


title: Kinetic fluid moments closure for a magnetized plasma with collisions

abstract: A novel method aimed at a kinetic moments closure for a magnetized plasma with arbitrary collisionality is proposed. The intended first application is to a tokamak edge and scrape-off-layer plasma. The velocity distribution function for each species is expanded in 8 Gaussian Radial Basis Functions (GRBFs) which are essentially shifted Maxwellians at eight representative 3D-velocity points of drift. The vector of 8 fluid moments (for particle density, 3 particle fluxes, total energy density, and 3 energy fluxes) has an 8x8 analytic linear matrix relation to the vector of 8 GRBF density weights in 3D real space. The 8 fluid moments with sources for each species are advanced in time while the 8 GRBF weighs are determined from the 8x8 inverse matrix. The two closure moments (for the stress tensor and the energy weighted stress tensor) are linearly determined from the GRBF weights. Most importantly the velocity moments of the nonlinear Coulomb Fokker-Planck collision operator [Rosenbluth et al, Phys. Rev 107, 1957] are evaluated from the GRBF weights. Generalization from 8 to 12, 16, 20 .., in an energy weighted moment hierarchy is straightforward. The electric field follows from a generalized vorticity (quasi-neutral current continuity) equation. A strong drift ordering approximation can be applied to eliminate any spuriously unstable high frequency cyclotron motions. A novel weak drift ordering two-time step scheme avoids the vorticity equation by following ion cyclotron motion in time to get the electric field with ion gyroaveraging. Inclusion of low-beta magnetic perturbations is straightforward.


title: Edge and bulk dissipative solitons in modulated PT-symmetric waveguide arrays

abstract: We address dissipative soliton formation in modulated PT-symmetric continuous waveguide arrays composed from waveguides with amplifying and absorbing sections, whose density gradually increases (due to decreasing waveguide separation) either towards the center of the array or towards its edges. In such a structure the level of gain/loss at which PT-symmetry gets broken depends on the direction of increase of the waveguide density. Breakup of the PT-symmetry occurs when eigenvalues of modes localized in the region, where waveguide density is largest, collide and move into complex plane. In this regime of broken symmetry the inclusion of focusing Kerr-type nonlinearity of the material and weak two-photon absorption allows to arrest the growth of amplitude of amplified modes and may lead to the appearance of stable attractors either in the center or at the edge of the waveguide array, depending on the type of array modulation. Such solitons can be stable, they acquire specific triangular shapes and notably broaden with increase of gain-loss level. Our results illustrate how spatial array modulation that breaks PT-symmetry locally can be used to control specific location of dissipative solitons forming in the array.


title: Regularity of spectral stacks and discreteness of weight-hearts

abstract: We study regularity in the context of ring spectra and spectral stacks. Parallel to that, we construct a weight structure on the category of compact quasi-coherent sheaves on spectral quotient stacks of the form $X=[\operatorname{Spec} R/G]$ defined over a field, where $R$ is a connective ${\mathcal{E}_{\infty}}$-$k$-algebra and $G$ is a linearly reductive group acting on $R$. Under reasonable assumptions we show that regularity of $X$ is equivalent to regularity of $R$. We also show that if $R$ is bounded, such a stack is discrete. This result can be interpreted in terms of weight structures and suggests a general phenomenon: for a symmetric monoidal stable $\infty$-category with a compatible bounded weight structure, the existence of an adjacent t-structure satisfying a strong boundedness condition should imply discreteness of the weight-heart. We also prove a gluing result for weight structures and adjacent t-structures, in the setting of a semi-orthogonal decomposition of stable $\infty$-categories.


title: On integrands and loop momentum in string and field theory

abstract: The notion of a unique integrand does not a priori makes sense in field theory: different Feynman diagrams have different loop momenta and there should be no reason to compare them. In string theory, however, a global integrand is natural and allows, for instance, to make explicit the separation between left and right-moving degrees of freedom. However, the significance of this integrand had not really been investigated so far. It is even more important in view of the recently discovered loop monodromies that are related to the duality between color and kinematics in gauge and gravity loop amplitudes. This paper intends to start filling this gap, by presenting a careful definition of the loop momentum in string theory, and describing precisely the resulting global integrand obtained in the field theory limit. We will then apply this technology to write down some monodromy relations at two and three loops, and make contact with the color/kinematics duality.


title: Learning with Collaborative Neural Network Group by Reflection

abstract: For the present engineering of neural systems, the preparing of extensive scale learning undertakings generally not just requires a huge neural system with a mind boggling preparing process yet additionally troublesome discover a clarification for genuine applications. In this paper, we might want to present the Collaborative Neural Network Group (CNNG). CNNG is a progression of neural systems that work cooperatively to deal with various errands independently in a similar learning framework. It is advanced from a solitary neural system by reflection. Along these lines, in light of various circumstances removed by the calculation, the CNNG can perform diverse techniques when handling the information. The examples of chose methodology can be seen by human to make profound adapting more reasonable. In our execution, the CNNG is joined by a few moderately little neural systems. We give a progression of examinations to assess the execution of CNNG contrasted with other learning strategies. The CNNG is able to get a higher accuracy with a much lower training cost. We can reduce the error rate by 74.5% and reached the accuracy of 99.45% in MNIST with three feedforward networks (4 layers) in one training epoch.


title: Sampled-Data Observers for 1-D Parabolic PDEs with Non-Local Outputs

abstract: The present work provides a systematic approach for the design of sampled-data observers to a wide class of 1-D, parabolic PDEs with non-local outputs. The studied class of parabolic PDEs allows the presence of globally Lipschitz nonlinear and non-local terms in the PDE. Two different sampled-data observers are presented: one with an inter-sample predictor for the unavailable continuous measurement signal and one without an inter-sample predictor. Explicit conditions on the upper diameter of the (uncertain) sampling schedule for both designs are derived for exponential convergence of the observer error to zero in the absence of measurement noise and modeling errors. Moreover, explicit estimates of the convergence rate can be deduced based on the knowledge of the upper diameter of the sampling schedule. When measurement noise and/or modeling errors are present, Input-to-Output Stability (IOS) estimates of the observer error hold for both designs with respect to noise and modeling errors. The main results are illustrated by two examples which show how the proposed methodology can be extended to other cases (e.g., boundary point measurements).


title: Finite Rogers--Ramanujan type identities

abstract: Polynomial generalizations of all 130 of the identities in Slater's list of identities of the Rogers-Ramanujan type are presented. Furthermore, duality relationships among many of the identities are derived. Some of the these polynomial identities were previously known but many are new. The author has implemented much of the finitization process in a Maple package which is available for free download from the author's website.


title: Analysis of Distributed ADMM Algorithm for Consensus Optimization in Presence of Node Error

abstract: Alternating Direction Method of Multipliers (ADMM) is a popular convex optimization algorithm, which can be employed for solving distributed consensus optimization problems. In this setting agents locally estimate the optimal solution of an optimization problem and exchange messages with their neighbors over a connected network. The distributed algorithms are typically exposed to different types of errors in practice, e.g., due to quantization or communication noise or loss. We here focus on analyzing the convergence of distributed ADMM for consensus optimization in presence of additive random node error, in which case, the nodes communicate a noisy version of their latest estimate of the solution to their neighbors in each iteration. We present analytical upper and lower bounds on the mean squared steady state error of the algorithm in case that the local objective functions are strongly convex and have Lipschitz continuous gradients. In addition we show that, when the local objective functions are convex and the additive node error is bounded, the estimation error of the noisy ADMM for consensus optimization is also bounded. Numerical results are provided which demonstrate the effectiveness of the presented analyses and shed light on the role of the system and network parameters on performance.


title: Photospheric plasma and magnetic field dynamics during the formation of solar AR 11190

abstract: The Sun features on its surface typical flow patterns called the granulation, mesogranulation, and supergranulation. These patterns arise due to convective flows transporting energy from the interior of the Sun to its surface. In this paper we will shed light on the interaction between the convective flows in large-scale cells as well as the large-scale magnetic fields in active regions, and investigate in detail the statistical distribution of flow velocities during the evolution and formation of National Oceanic and Atmospheric Administration (NOAA) active region 11190. To do so, we employed local correlation tracking methods on data obtained by the Solar Dynamics Observatory (SDO) spacecraft in the continuum as well as on processed line-of-sight (LOS) magnetograms. We find that the flow fields in an active region can be modelled by a two-component distribution. One component is very stable, follows a Rayleigh distribution, and can be assigned to the background flows, whilst the other component is variable in strength and velocity range and can be attributed to the flux emergence visible both in the continuum maps as well as magnetograms. Generally, the plasma flows, as seen by the distribution of the magnitude of the velocity, follow a Rayleigh distribution even through the time of formation of active regions. However, at certain moments of large-scale fast flux emergence, a second component featuring higher velocities is formed in the velocity magnitudes distribution. The plasma flows are generally highly correlated to the motion of magnetic elements and vice versa except during the times of fast magnetic flux emergence as observed by rising magnetic elements. At these times, the magnetic fields are found to move faster than the corresponding plasma.


title: Using fuzzy bits and neural networks to partially invert few rounds of some cryptographic hash functions

abstract: We consider fuzzy, or continuous, bits, which take values in [0;1] and (-1;1] instead of {0;1}, and operations on them (NOT, XOR etc.) and on their sequences (ADD), to obtain the generalization of cryptographic hash functions, CHFs, for the messages consisting of fuzzy bits, so that CHFs become smooth and non-constant functions of each bit of the message. We then train the neural networks to predict the message that has a given hash, where the loss function for the hash of predicted message and given true hash is backpropagatable. The results of the trainings for the standard CHFs - MD5, SHA1, SHA2-256, and SHA3/Keccak - with small number of (optionally weakened) rounds are presented and compared.


title: Counting twisted Higgs bundles

abstract: We count invariants of the moduli spaces of twisted Higgs bundles on a smooth projective curve.


title: Microscopic theory of spin transport at the interface between the superconductor and a ferromagnetic insulator

abstract: We theoretically investigate spin transport at the interface between the ferromagnetic insulator(FI) and a superconductor(SC). Considering a simple FI-SC interface model, we derive formulas for the spin current and spin-current noise induced by microwave irradiation (spin pumping) or the temperature gradient (the spin Seebeck effect). We show how the superconducting coherence factor affects the temperature dependence of the spin current. We also calculate the spin-current noise in thermal equilibrium and in non-equilibrium states induced by the spin pumping, and compare them quantitatively for an yttrium-iron-garnet-NbN interface.


title: Lower bounds for maximal matchings and maximal independent sets

abstract: There are distributed graph algorithms for finding maximal matchings and maximal independent sets in $O(\Delta + \log^* n)$ communication rounds; here $n$ is the number of nodes and $\Delta$ is the maximum degree. The lower bound by Linial (1987, 1992) shows that the dependency on $n$ is optimal: these problems cannot be solved in $o(\log^* n)$ rounds even if $\Delta = 2$. However, the dependency on $\Delta$ is a long-standing open question, and there is currently an exponential gap between the upper and lower bounds. We prove that the upper bounds are tight. We show that any algorithm that finds a maximal matching or maximal independent set with probability at least $1-1/n$ requires $\Omega(\min\{\Delta,\log \log n / \log \log \log n\})$ rounds in the LOCAL model of distributed computing. As a corollary, it follows that any deterministic algorithm that finds a maximal matching or maximal independent set requires $\Omega(\min\{\Delta, \log n / \log \log n\})$ rounds; this is an improvement over prior lower bounds also as a function of $n$.


title: Stable Electromyographic Sequence Prediction During Movement Transitions using Temporal Convolutional Networks

abstract: Transient muscle movements influence the temporal structure of myoelectric signal patterns, often leading to unstable prediction behavior from movement-pattern classification methods. We show that temporal convolutional network sequential models leverage the myoelectric signal's history to discover contextual temporal features that aid in correctly predicting movement intentions, especially during interclass transitions. We demonstrate myoelectric classification using temporal convolutional networks to effect 3 simultaneous hand and wrist degrees-of-freedom in an experiment involving nine human-subjects. Temporal convolutional networks yield significant $(p<0.001)$ performance improvements over other state-of-the-art methods in terms of both classification accuracy and stability.


title: Production of rotationally cold methyl radicals in pulsed supersonic beams

abstract: We present a comparison of two technically distinct methods for the generation of rotationally cold, pulsed supersonic beams of methyl radicals (CH3): a plate discharge source operating in the glow regime, and a dielectric barrier discharge source (DBD). The results imply that the efficiency of both sources is comparable, and that molecular beams with similar translational and rotational temperatures are formed. Methane (CH4) proved to be the most suitable radical precursor species.


title: Unseen Object Segmentation in Videos via Transferable Representations

abstract: In order to learn object segmentation models in videos, conventional methods require a large amount of pixel-wise ground truth annotations. However, collecting such supervised data is time-consuming and labor-intensive. In this paper, we exploit existing annotations in source images and transfer such visual information to segment videos with unseen object categories. Without using any annotations in the target video, we propose a method to jointly mine useful segments and learn feature representations that better adapt to the target frames. The entire process is decomposed into two tasks: 1) solving a submodular function for selecting object-like segments, and 2) learning a CNN model with a transferable module for adapting seen categories in the source domain to the unseen target video. We present an iterative update scheme between two tasks to self-learn the final solution for object segmentation. Experimental results on numerous benchmark datasets show that the proposed method performs favorably against the state-of-the-art algorithms.


title: On the nonexistence of F{\o}lner sets

abstract: We show that there is $n\in \mathbf N$, a finite system $\Sigma(\vec x,\vec y)$ of equations and inequations having a solution in some group, where $\vec x$ has length $n$, and $\epsilon>0$ such that: for any group $G$ and any $\vec a\in G^n$, if the system $\Sigma(\vec a,\vec y)$ has a solution in $G$, then there is no $(\vec a,\epsilon)$-F{\o} lner set in $G$. The proof uses ideas from model-theoretic forcing together with the observation that no amenable group can be existentially closed. Along the way, we also observe that no existentially closed group can be exact, have the Haagerup property, or have property (T). Finally, we show that, for $n$ large enough and for $\epsilon$ small enough, the existence of $(F,\epsilon)$-F{\o} lner sets, where $F$ has size at most $n$, cannot be expressed in a first-order way uniformly in all groups.


title: Panoptic Feature Pyramid Networks

abstract: The recently introduced panoptic segmentation task has renewed our community's interest in unifying the tasks of instance segmentation (for thing classes) and semantic segmentation (for stuff classes). However, current state-of-the-art methods for this joint task use separate and dissimilar networks for instance and semantic segmentation, without performing any shared computation. In this work, we aim to unify these methods at the architectural level, designing a single network for both tasks. Our approach is to endow Mask R-CNN, a popular instance segmentation method, with a semantic segmentation branch using a shared Feature Pyramid Network (FPN) backbone. Surprisingly, this simple baseline not only remains effective for instance segmentation, but also yields a lightweight, top-performing method for semantic segmentation. In this work, we perform a detailed study of this minimally extended version of Mask R-CNN with FPN, which we refer to as Panoptic FPN, and show it is a robust and accurate baseline for both tasks. Given its effectiveness and conceptual simplicity, we hope our method can serve as a strong baseline and aid future research in panoptic segmentation.


title: Thermal Boundary Conductance of Two-Dimensional $MoS_{2}$ Interfaces

abstract: Understanding the thermal properties of two-dimensional (2D) materials and devices is essential for thermal management of 2D applications. Here we perform molecular dynamics simulations to evaluate both the specific heat of $MoS_{2}$ as well as the thermal boundary conductance (TBC) between one to five layers of $MoS_{2}$ with amorphous $SiO_{2}$ and between single-layer $MoS_{2}$ and crystalline $AlN$. The results of all calculations are compared to existing experimental data. In general, the TBC of such 2D interfaces is low, below ~20 $MWm^{-2}K^{-1}$, due to the weak van der Waals (vdW) coupling and mismatch of phonon density of states (PDOS) between materials. However, the TBC increases with vdW coupling strength, with temperature, and with the number of $MoS_{2}$ layers (which introduce additional phonon modes). These findings suggest that the TBC of 2D materials is tunable by modulating their interface interaction, the number of layers, and finding a PDOS-matched substrate, with important implications for future energy-efficient 2D electronics, photonics, and thermoelectrics.


title: Magnetic island merger as a mechanism for inverse magnetic energy transfer

abstract: Magnetic energy transfer from small to large scales due to successive magnetic island coalescence is investigated. A solvable analytical model is introduced and shown to correctly capture the evolution of the main quantities of interest, as borne out by numerical simulations. Magnetic reconnection is identified as the key mechanism enabling the inverse transfer, and setting its properties: magnetic energy decays as $\tilde t^{-1}$, where $\tilde t$ is time normalized to the (appropriately defined) reconnection timescale; and the correlation length of the field grows as $\tilde t^{1/2}$. The magnetic energy spectrum is self-similar, and evolves as $\propto \tilde t ^{-3/2}k^{-2}$, where the $k$-dependence is imparted by the formation of thin current sheets.


title: Zero modes and low-energy resolvent expansion for three dimensional Schrodinger operators with point interactions

abstract: We investigate the low-energy behavior of the resolvent of Schrodinger operators with finitely many point interactions in three dimensions. We also discuss the occurrence and the multiplicity of zero energy obstructions.


title: A C-DAG task model for scheduling complex real-time tasks on heterogeneous platforms: preemption matters

abstract: Recent commercial hardware platforms for embedded real-time systems feature heterogeneous processing units and computing accelerators on the same System-on-Chip. When designing complex real-time application for such architectures, the designer needs to make a number of difficult choices: on which processor should a certain task be implemented? Should a component be implemented in parallel or sequentially? These choices may have a great impact on feasibility, as the difference in the processor internal architectures impact on the tasks' execution time and preemption cost. To help the designer explore the wide space of design choices and tune the scheduling parameters, in this paper we propose a novel real-time application model, called C-DAG, specifically conceived for heterogeneous platforms. A C-DAG allows to specify alternative implementations of the same component of an application for different processing engines to be selected off-line, as well as conditional branches to model if-then-else statements to be selected at run-time. We also propose a schedulability analysis for the C-DAG model and a heuristic allocation algorithm so that all deadlines are respected. Our analysis takes into account the cost of preempting a task, which can be non-negligible on certain processors. We demonstrate the effectiveness of our approach on a large set of synthetic experiments by comparing with state of the art algorithms in the literature.


title: Investigation of dissipative particle dynamics with colored noise

abstract: We investigate the behavior of dissipative particle dynamics (DPD) with time-correlated random noise. A new stochastic force for DPD is proposed which consists of a random force whose noise has an algebraic correlation proportional to 1/t and is generated by the so called Kangaroo process. We stress the benefits of a time correlated noise in stochastic systems. We show that the system exhibits significantly different properties from classical DPD, driven by Wiener noise. While the probability distribution function of the velocity is Gaussian, the acceleration develops a bi-modal character. Although the fluctuation dissipation theorem may not strictly hold, we demonstrate that the system reaches equilibrium states with fluctuation-dissipation balance. We believe that our explorative research on the DPD model may stimulate the application of modified DPD to unconventional problems beyond molecular modeling.


title: Face Recognition System

abstract: Deep learning is one of the new and important branches in machine learning. Deep learning refers to a set of algorithms that solve various problems such as images and texts by using various machine learning algorithms in multi-layer neural networks. Deep learning can be classified as a neural network from the general category, but there are many changes in the concrete realization. At the core of deep learning is feature learning, which is designed to obtain hierarchical information through hierarchical networks, so as to solve the important problems that previously required artificial design features. Deep Learning is a framework that contains several important algorithms. For different applications (images, voice, text), you need to use different network models to achieve better results. With the development of deep learning and the introduction of deep convolutional neural networks, the accuracy and speed of face recognition have made great strides. However, as we said above, the results from different networks and models are very different. In this paper, facial features are extracted by merging and comparing multiple models, and then a deep neural network is constructed to train and construct the combined features. In this way, the advantages of multiple models can be combined to mention the recognition accuracy. After getting a model with high accuracy, we build a product model. This article compares the pure-client model with the server-client model, analyzes the pros and cons of the two models, and analyzes the various commercial products that are required for the server-client model.


title: Neural Inverse Rendering of an Indoor Scene from a Single Image

abstract: Inverse rendering aims to estimate physical attributes of a scene, e.g., reflectance, geometry, and lighting, from image(s). Inverse rendering has been studied primarily for single objects or with methods that solve for only one of the scene attributes. We propose the first learning-based approach that jointly estimates albedo, normals, and lighting of an indoor scene from a single image. Our key contribution is the Residual Appearance Renderer (RAR), which can be trained to synthesize complex appearance effects (e.g., inter-reflection, cast shadows, near-field illumination, and realistic shading), which would be neglected otherwise. This enables us to perform self-supervised learning on real data using a reconstruction loss, based on re-synthesizing the input image from the estimated components. We finetune with real data after pretraining with synthetic data. To this end, we use physically-based rendering to create a large-scale synthetic dataset, which is a significant improvement over prior datasets. Experimental results show that our approach outperforms state-of-the-art methods that estimate one or more scene attributes.


title: Reconcilement of VHE $\gamma$-ray/X-ray correlation studies in Mrk 421 and break-down at high fluxes

abstract: Multi-wavelength campaigns have been carried out to study the correlation between the very high energy (VHE) $\gamma$-ray and the X-ray emissions in blazars but, no conclusive results have been achieved yet. In this paper, we add Milagro data to the existing VHE $\gamma$-ray data from HEGRA-CT1 and Whipple and test the consistency and robustness of the reported correlation between VHE $\gamma$-ray and X-ray fluxes in Mrk 421. We found that at monthly time scale the correlation is robust, consistent between instruments and described as a linear function. Furthermore, most of the fluxes on shorter time scales are consistent with the correlation within 3 $\sigma_A$ even, where $\sigma_A$ is an estimated intrinsic scatter. However, a break-down of the correlation becomes clearly evident at high states of activity with fluxes $\rm \gtrsim 2.5\times 10^{-10}\, cm^{-2}s^{-1}$ at energies above 400 GeV independently of the time scale, observational period or instrument, even for single flares, the X-ray and VHE $\gamma$-ray emissions lie on the correlation until the VHE $\gamma$-ray flux reaches values higher than the one mentioned above. We have interpreted our results within the one-zone synchrotron self-Compton model. We found that describing a single and unique $\gamma$-ray/X-ray correlation strongly narrows the range of possible values of the magnetic field $B$ when a constant value of the spectral index along the correlation is assumed.


title: Computational aberration compensation by coded aperture-based correction of aberration obtained from Fourier ptychography (CACAO-FP)

abstract: We report a novel generalized optical measurement system and computational approach to determine and correct aberrations in optical systems. We developed a computational imaging method capable of reconstructing an optical system's aberration by harnessing Fourier ptychography (FP) with no spatial coherence requirement. It can then recover the high resolution image latent in the aberrated image via deconvolution. Deconvolution is made robust to noise by using coded apertures to capture images. We term this method: coded aperture-based correction of aberration obtained from Fourier ptychography (CACAO-FP). It is well-suited for various imaging scenarios with the presence of aberration where providing a spatially coherent illumination is very challenging or impossible. We report the demonstration of CACAO-FP with a variety of samples including an in-vivo imaging experiment on the eye of a rhesus macaque to correct for its inherent aberration in the rendered retinal images. CACAO-FP ultimately allows for a poorly designed imaging lens to achieve diffraction-limited performance over a wide field of view by converting optical design complexity to computational algorithms in post-processing.


title: Magnetoelectric Coupling by Giant Piezoelectric Tensor Design

abstract: Strain-coupled magnetoelectric (ME) phenomena in piezoelectric / ferromagnetic thin-film bilayers are a promising paradigm for sensors and information storage devices, where strain is utilized to manipulate the magnetization of the ferromagnetic film. In-plane magnetization rotation with an electric field across the film thickness has been challenging due to the virtual elimination of in-plane piezoelectric strain by substrate clamping, and to the requirement of anisotropic in-plane strain in two-terminal devices. We have overcome both of these limitations by fabricating lithographically patterned devices with a piezoelectric membrane on a soft substrate platform, in which in-plane strain is freely generated, and a patterned edge constraint that transforms the nominally isotropic piezoelectric strain into the required uniaxial strain. We fabricated 500 nm thick, (001) oriented [Pb(Mg$_{1/3}$Nb$_{2/3}$)O$_3$]$_{0.7}$-[PbTiO$_3$]$_{0.3}$ (PMN-PT) unclamped piezoelectric membranes with ferromagnetic Ni overlayers. Guided by analytical and numerical continuum elastic calculations, we designed and fabricated two-terminal devices exhibiting Ni magnetization rotation in response to an electric field across the PMN-PT. Similar membrane heterostructures could be used to apply designed strain patterns to many other materials systems to control properties such as superconductivity, band topology, conductivity, and optical response.


title: Are LGRBs biased tracers of star formation? Clues from the host galaxies of the $Swift$/BAT6 complete sample of bright LGRBs III: Stellar masses, star formation rates and metallicities at $z>1$

abstract: (Abridged) Long gamma-ray bursts (LGRB) have been suggested as promising tracers of star formation owing to their association with the core-collapse of massive stars. The goal of this work is to characterise the population of host galaxies of LGRBs at 1 < z < 2, investigate the conditions in which LGRBs form at these redshifts and assess their use as tracers of star formation. We perform a spectro-photometric analysis to determine the stellar mass, star formation rate, specific star formation rate and metallicity of the complete, unbiased host galaxy sample of the Swift/BAT6 LGRB sample at 1 < z < 2. We compare the distribution of these properties to the ones of typical star-forming galaxies from the MOSDEF and COSMOS2015 Ultra Deep surveys, within the same redshift range. We find that, similarly to z < 1, LGRBs do not directly trace star formation at 1 < z < 2, and they tend to avoid high-mass, high-metallicity host galaxies. We also find evidence for an enhanced fraction of starbursts among the LGRB host sample with respect to the star-forming population of galaxies. Nonetheless we demonstrate that the driving factor ruling the LGRB efficiency is metallicity. The LGRB host distributions can be reconciled with the ones expected from galaxy surveys by imposing a metallicity upper limit of 12+logOH ~ 8.55. Metallicity rules the LGRB production efficiency, which is stifled at Z > 0.7 Zsun. Under this hypothesis we can expect LGRBs to trace star formation at z > 3, once the bulk of the star forming galaxy population are characterised by metallicities below this limit. The moderately high metallicity threshold found is in agreement with the conditions necessary to rapidly produce a fast-rotating Wolf-Rayet star a in close binary system, and could be accommodated by single star models under chemically homogeneous mixing with very rapid rotation and weak magnetic coupling.


title: A Possible Relic Star Cluster in the Sextans Dwarf Galaxy

abstract: We report a possible discovery of a relic star cluster in the Sextans dwarf spheroidal galaxy. Using the \textit{hk} index ($\equiv$($Ca-b$)$-$($b-y$)) as a photometric metallicity indicator, we discriminate the metal-poor and metal-rich stars in the galaxy and find unexpected number density excess of metal-poor stars located 7'.7 ($\sim$190 pc in projected distance) away from the known galactic center. The $V-I$ color$-$magnitude diagram (CMD) for stars around the density excess reveals that both the main sequence and the giant branch are considerably narrower and redder than the bulk of field stars in Sextans. Our stellar population models show (a) that the narrow CMD is best reproduced by a simple stellar population with an age of $\sim$13 Gyr and [Fe/H] of $\sim$$-$2.3 dex, and (b) that the redder $V-I$ color of the $hk$-weak population is explained $only$ if it is $\sim$2 Gyr older than the field stars. The results lead us to conclude that the off-centered density peak is likely associated with an old, metal-poor globular cluster. The larger spatial extent ($>$ 80 pc in radius) and the smaller number of stars ($\sim$1000) than typical globular clusters point to a star cluster that is in the process of dissolution. The finding serves as the first detection of a surviving star cluster in Sextans, supporting previous suggestions of the presence of star cluster remnants in the galaxy. If confirmed, the survival of a relic star cluster until now implies a $cored$ dark matter halo profile for this dwarf galaxy.


title: Magnetic Skyrmion Lattice by Fourier Transform Method

abstract: We demonstrate a fast numerical method of theoretical studies of skyrmion lattice or spiral order in magnetic materials with Dzyaloshinsky-Moriya interaction. The method is based on the Fourier expansion of the magnetization combined with a minimization of the free energy functional of the magnetic material in Fourier space, yielding the optimal configuration of the system for any given set of parameters. We employ a Lagrange multiplier technique in order to satisfy micromagnetic constraints. We apply this method to a system that exhibits, depending on the parameter choice, ferromagnetic, skyrmion lattice, or spiral (helical) order. Known critical fields corresponding to the helical-skyrmion as well as the skyrmion-ferromagnet phase transitions are reproduced with high precision. Using this numerical method we predict new types of excited (metastable) states of the skyrmion lattice, which may be stabilized by coupling the skyrmion lattice with a superconducting vortex lattice. The method can be readily adapted to other micromagnetic systems.


title: Handling the Uncertainties in the Galactic Dark Matter Distribution for Particle Dark Matter Searches

abstract: In this work we characterize the distribution of Dark Matter (DM) in the Milky Way (MW), and its uncertainties, adopting the well known "Rotation Curve" method. We perform a full marginalization over the uncertainties of the Galactic Parameters and over the lack of knowledge on the morphology of the baryonic components of the Galaxy. The local DM density rho0 is constrained to the range 0.3 - 0.8 GeV/cm3 at the 2 sigma level, and has a strong positive correlation to R0, the local distance from the Galactic Center (GC). The not well-known value of R0 is thus, at the moment, a major limitation in determining rho0. Similarly, we find that the inner slope of the DM profile, gamma, is very weakly constrained, showing no preference for a cored profile (gamma~0) or a cuspy one (gamma~[1.0,1.4]). Some combination of parameters can be, however, strongly constrained. For example the often used standard rho0=0.3 GeV/cm3, R0=8.5 kpc is excluded at more than 4 sigma. We release the full likelihood of our analysis in a tabular form over a multidimensional grid in the parameters characterizing the DM distribution, namely the scale radius Rs, the scale density rhos, the inner slope of the profile gamma, and R0. The likelihood can be used to include the effect of the DM distribution uncertainty on the results of searches for an indirect DM signal in gamma-rays or neutrinos, from the GC, or the Halo region surrounding it. As one example, we study the case of the GC excess in gamma rays. Further applications of our tabulated uncertainties in the DM distribution involve local DM searches, like direct detection and anti-matter observations, or global fits combining local and GC searches.


title: Electromagnetism and hidden vector fields in modified gravity theories: spontaneous and induced vectorization

abstract: In general relativity, Maxwell's equations are embedded in curved spacetime through the minimal prescription, but this could change if strong-gravity modifications are present. We show that with a nonminimal coupling between gravity and a massless vector field, nonperturbative effects can arise in compact stars. We find solutions describing stars with nontrivial vector field configurations, some of which are associated to an instability, while others are not. The vector field can be interpreted either as the electromagnetic field, or as a hidden vector field weakly coupled with the standard model.


title: Experimental quantum stochastic walks simulating associative memory of Hopfield neural networks

abstract: With the increasing crossover between quantum information and machine learning, quantum simulation of neural networks has drawn unprecedentedly strong attention, especially for the simulation of associative memory in Hopfield neural networks due to their wide applications and relatively simple structures that allow for easier mapping to the quantum regime. Quantum stochastic walk, a strikingly powerful tool to analyze quantum dynamics, has been recently proposed to simulate the firing pattern and associative memory with a dependence on Hamming Distance. We successfully map the theoretical scheme into a three-dimensional photonic quantum chip and realize quantum stochastic walk evolution through well-controlled detunings of the propagation constant. We demonstrate a good match rate of the associative memory between the experimental quantum scheme and the expected result for Hopfield neural networks. The ability of quantum simulation for an important feature of a neural network, combined with the scalability of our approach through low-loss integrated chip and straightforward Hamiltonian engineering, provides a primary but steady step towards photonic artificial intelligence devices for optimization and computation tasks of greatly improved efficiencies.