title: Efficient limited-time reachability estimation in temporal networks

abstract: Time-limited states characterise many dynamical processes on networks: disease infected individuals recover after some time, people forget news spreading on social networks, or passengers may not wait forever for a connection. These dynamics can be described as limited waiting-time processes, and they are particularly important for systems modelled as temporal networks. These processes have been studied via simulations, which is equivalent to repeatedly finding all limited-waiting time temporal paths from a source node and time. We propose a method yielding orders of magnitude more efficient way of tracking the reachability of such temporal paths. Our method gives simultaneous estimates of the in- or out-reachability (with any chosen waiting-time limit) from every possible starting point and time. It works on very large temporal networks with hundreds of millions of events on current commodity computing hardware. This opens up the possibility to analyse reachability and dynamics of spreading processes on large temporal networks in completely new ways. For example, one can now compute centralities based on global reachability for all events or can find with high probability the infected node and time, which would lead to the largest epidemic outbreak.


title: Topology via Spectral Projectors with Staggered Fermions

abstract: The spectral projectors method is a way to obtain a theoretically well posed definition of the topological susceptibility on the lattice. Up to now this method has been defined and applied only to Wilson fermions. The goal of this work is to extend the method to staggered fermions, giving a definition for the staggered topological susceptibility and testing it in the pure $SU(3)$ gauge theory. Besides, we also generalize the method to higher-order cumulants of the topological charge distribution.


title: Network Elastic Net for Identifying Smoking specific gene expression for lung cancer

abstract: Survival month for non-small lung cancer patients depend upon which stage of lung cancer is present. Our aim is to identify smoking specific gene expression biomarkers in the prognosis of lung cancer patients. In this paper, we introduce the network elastic net, a generalization of network lasso that allows for simultaneous clustering and regression on graphs. In Network elastic net, we consider similar patients based on smoking cigarettes per year to form the network. We then further find the suitable cluster among patients based on coefficients of genes having different survival month structures and showed the efficacy of the clusters using stage enrichment. This can be used to identify the stage of cancer using gene expression and smoking behavior of patients without doing any tests.


title: Rethinking Irregular Scene Text Recognition

abstract: Reading text from natural images is challenging due to the great variety in text font, color, size, complex background and etc.. The perspective distortion and non-linear spatial arrangement of characters make it further difficult. While rectification based method is intuitively grounded and has pushed the envelope by far, its potential is far from being well exploited. In this paper, we present a bag of tricks that prove to significantly improve the performance of rectification based method. On curved text dataset, our method achieves an accuracy of 89.6% on CUTE-80 and 76.3% on Total-Text, an improvement over previous state-of-the-art by 6.3% and 14.7% respectively. Furthermore, our combination of tricks helps us win the ICDAR 2019 Arbitrary-Shaped Text Challenge (Latin script), achieving an accuracy of 74.3% on the held-out test set. We release our code as well as data samples for further exploration at https://github.com/Jyouhou/ICDAR2019-ArT-Recognition-Alchemy


title: A Decentralized Primal-dual Method for Constrained Minimization of a Strongly Convex Function

abstract: We propose decentralized primal-dual methods for cooperative multi-agent consensus optimization problems over both static and time-varying communication networks, where only local communications are allowed. The objective is to minimize the sum of agent-specific convex functions over conic constraint sets defined by agent-specific nonlinear functions; hence, the optimal consensus decision should lie in the intersection of these private sets. Under the strong convexity assumption, we provide convergence rates for sub-optimality, infeasibility, and consensus violation in terms of the number of communications required; examine the effect of underlying network topology on the convergence rates.


title: VERITAS Observations of 1ES 1218+304 during the 2019 VHE High State

abstract: 1ES 1218+304 is a moderate-redshift (z = 0.182) high-frequency-peaked BL Lac object (HBL). Detected by both the MAGIC and VERITAS gamma-ray observatories, 1ES 1218+304 is frequently monitored by VERITAS as part of its long-term blazar monitoring program. On the 3rd of January 2019, during a regularly scheduled blazar snapshot, VERITAS observed 1ES 1218+304 to have an elevated VHE flux state which continued through January 5 (ATEL #12360). In addition to VERITAS observations, MAGIC also detected elevated VHE emission (ATEL #12354) and Swift-BAT/XRT and Tuorla observations, in X-ray and optical respectively, showed increased multiwavelength activity at the time of the VERITAS observations. In this work, VERITAS observations of 1ES 1218+304 taken as part of a target of opportunity campaign during this elevated state, are presented.


title: Maximum likelihood spectral fitting and its application to EBL constraints

abstract: The extragalactic background light (EBL) is the second-most-intense form of cosmic background light (the first being the cosmic microwave background) and contains the redshifted optical radiation, from infra-red to ultraviolet, emitted across all epochs, making it of great cosmological interest. While direct measurements of the EBL are hampered by foreground contamination, observations of VHE emission from distant sources can be used to obtain indirect measurements of the EBL. In this work a maximum-likelihood fit is applied to the energy spectra of blazars observed by VERITAS, an array of ground-based imaging atmospheric Cherenkov telescopes sensitive to very-high-energy (VHE; E>100 GeV) gamma rays. Using theoretical models of the EBL shape and intensity, the EBL normalization is treated as a free parameter, allowing for model-dependent constraints to be obtained. Details of this maximum-likelihood analysis and preliminary model-dependent constraints on the EBL, are presented.


title: Dark energy and spontaneous mirror symmetry breaking

abstract: Dark energy is interpreted as the leftover of mostly canceled vacuum energy due to the spontaneous mirror symmetry breaking (SMSB) at the electroweak phase transition. Based on the newly proposed mirror-matter model (M$^3$), the extended standard model with mirror matter (SM$^3$) is elaborated to provide a consistent foundation for understanding dark energy, dark matter, baryogenesis, and many other puzzles. New insights of Higgs, top quark, and lepton masses are presented under SM$^3$ using staged quark condensation and four-fermion interactions for SMSB. In particular, the nature and mass scales of neutrinos are naturally explained under the new theory. The new cosmology model based on SM$^3$ could potentially resolve more cosmic enigmas. The possible underlying principles for SMSB and SM$^3$ of a maximally interacting, supersymmetric, and mirrored world are also discussed.


title: Long rotational coherence times of molecules in a magnetic trap

abstract: Polar molecules in superpositions of rotational states exhibit long-range dipolar interactions, but maintaining their coherence in a trapped sample is a challenge. We present calculations that show many laser-coolable molecules have convenient rotational transitions that are exceptionally insensitive to magnetic fields. We verify this experimentally for CaF where we find a transition with sensitivity below 5 Hz G$^{-1}$ and use it to demonstrate a rotational coherence time of 6.4(8) ms in a magnetic trap. Simulations suggest it is feasible to extend this to more than 1 s using a smaller cloud in a biased magnetic trap.


title: Long exit times near a repelling equilibrium

abstract: For a smooth vector field in a neighborhood of a critical point with all positive eigenvalues of the linearization, we consider the associated dynamics perturbed by white noise. Using Malliavin calculus tools, we obtain polynomial asymptotics for probabilities of atypically long exit times in the vanishing noise limit.


title: CodeSwitch-Reddit: Exploration of Written Multilingual Discourse in Online Discussion Forums

abstract: In contrast to many decades of research on oral code-switching, the study of written multilingual productions has only recently enjoyed a surge of interest. Many open questions remain regarding the sociolinguistic underpinnings of written code-switching, and progress has been limited by a lack of suitable resources. We introduce a novel, large, and diverse dataset of written code-switched productions, curated from topical threads of multiple bilingual communities on the Reddit discussion platform, and explore questions that were mainly addressed in the context of spoken language thus far. We investigate whether findings in oral code-switching concerning content and style, as well as speaker proficiency, are carried over into written code-switching in discussion forums. The released dataset can further facilitate a range of research and practical activities.


title: Constraints on Skyrme Equations of State from Doubly Magic Nuclei, Ab-Initio Calculations of Low-Density Neutron Matter, and Neutron Stars

abstract: We use properties of doubly-magic nuclei, ab-initio calculations of low-density neutron matter, and of neutron stars to constrain the parameters of the Skyrme energy-density functional. We find all of these properties can be reproduced within a constrained family of Skyrme parameters. The maximum mass of a neutron star is found to be sensitive to the neutron effective mass. A value of [$ m^{*}_{\rm n}/m](\rho_0) = 0.60-0.65 $ is required to obtain a maximum neutron star mass of 2.1 solar masses. Using the constrained Skyrme functional with the aforementioned effective mass, the predicted radius for a neutron star of 1.4 solar masses is 12.4(1) km and $\Lambda$ = 423(40).


title: Partitioned integrators for thermodynamic parameterization of neural networks

abstract: Traditionally, neural networks are parameterized using optimization procedures such as stochastic gradient descent, RMSProp and ADAM. These procedures tend to drive the parameters of the network toward a local minimum. In this article, we employ alternative "sampling" algorithms (referred to here as "thermodynamic parameterization methods") which rely on discretized stochastic differential equations for a defined target distribution on parameter space. We show that the thermodynamic perspective already improves neural network training. Moreover, by partitioning the parameters based on natural layer structure we obtain schemes with very rapid convergence for data sets with complicated loss landscapes. We describe easy-to-implement hybrid partitioned numerical algorithms, based on discretized stochastic differential equations, which are adapted to feed-forward neural networks, including a multi-layer Langevin algorithm, AdLaLa (combining the adaptive Langevin and Langevin algorithms) and LOL (combining Langevin and Overdamped Langevin); we examine the convergence of these methods using numerical studies and compare their performance among themselves and in relation to standard alternatives such as stochastic gradient descent and ADAM. We present evidence that thermodynamic parameterization methods can be (i) faster, (ii) more accurate, and (iii) more robust than standard algorithms used within machine learning frameworks.


title: Zero index electromagnetic materials

abstract: Here we re-examine one of the most basic quantities in optics: the refractive index. Considering propagation in a plane, we first develop a general formalism for calculating the direction dependent refractive index in a general bi-anisotropic material. From this we derive the general condition for achieving zero refractive index in a given direction. We show that when the zero-index direction is complex valued the material supports waves that can propagate in only one sense, e.g. in only a clockwise direction. We dub such materials complex axis nihility (CAN) media. Our condition shows that there are an infinite family of both time reversible and time irreversible homogeneous electromagnetic media that support unidirectional propagation for a particular polarization. We give examples showing that scattering from such media results in the complete exclusion of partial waves with one sign of the angular momentum, and that interfaces between such media generally support one-way interface states. As well as giving new sets of material parameters, our simple condition reproduces many of the findings derived using topology, such as unidirectional propagation in gyrotropic media, and spin-momentum locking of evanescent waves.


title: On the Hardy number of a domain in terms of harmonic measure and hyperbolic distance

abstract: Let $\psi $ be a conformal map on $\mathbb{D}$ with $ \psi \left( 0 \right)=0$ and let ${F_\alpha }=\left\{ {z \in \mathbb{D}:\left| {\psi \left( z \right)} \right| = \alpha } \right\}$ for $\alpha >0$. Denote by ${H^p}\left( \mathbb{D} \right)$ the classical Hardy space with exponent $p>0$ and by ${\tt h}\left( \psi \right)$ the Hardy number of $\psi$. Consider the limits $$ L:= \lim_{\alpha\to+\infty}\left( \log \omega_{\mathbb D}(0,F_{\alpha})^{-1}/\log \alpha \right), \,\, \mu:= \lim_{\alpha\to+\infty}\left( d_{\mathbb D}(0,F_{\alpha})/\log\alpha \right),$$ where $\omega _\mathbb{D}\left( {0,{F_\alpha }} \right)$ denotes the harmonic measure at $0$ of $F_\alpha $ and $d_\mathbb{D} {\left( {0,{F_\alpha }} \right)}$ denotes the hyperbolic distance between $0$ and $F_\alpha$ in $\mathbb{D}$. We study a problem posed by P. Poggi-Corradini. What is the relation between $L$, $\mu$ and ${\tt h}\left( \psi \right)$? We also provide conditions for the existence of $L$ and $\mu$ and for the equalities $L=\mu={\tt h}\left( \psi \right)$. Poggi-Corradini proved that $\psi \notin {H^{\mu}}\left( \mathbb{D} \right)$ for a wide class of conformal maps $\psi$. We present an example of $\psi$ such that $\psi \in {H^\mu {\left( \mathbb{D} \right)} }$.


title: An Incremental Evaluation Mechanism for the Critical Node Problem

abstract: The Critical Node Problem (CNP) is to identify a subset of nodes in a graph whose removal maximally degrades pairwise connectivity. The CNP is an important variant of the Critical Node Detection Problem (CNDP) with wide applications. Due to its NP-hardness for general graphs, most works focus on local search algorithms that can return a good quality solution in a reasonable time. However, computing the objective function of CNP is a frequent procedure and is time-consuming (with complexity O(|V | + |E|)) during the search, which is a common problem that previous algorithms suffered from. In this paper, we propose a general incremental evaluation mechanism (IEM) to compute the objective function with much lower complexity. In this work, we improved two important greedy operations with IEM, along with experiments. Finally, we evaluate IEM by applying it into an evolutionary algorithm on two popular benchmarks, compared with the state-of-the-art approach. The experimental results showed the significance of IEM.


title: Strong first-order phase transitions in the NMSSM --- a comprehensive survey

abstract: Motivated by the fact that the Next-to-Minimal Supersymmetric Standard Model is one of the most plausible models that can accommodate electroweak baryogenesis, we analyze its phase structure by tracing the temperature dependence of the minima of the effective potential. Our results reveal rich patterns of phase structure that end in the observed electroweak symmetry breaking vacuum. We classify these patterns according to the first transition in their history and show the strong first-order phase transitions that may be possible in each type of pattern. These could allow for the generation of the matter-antimatter asymmetry or potentially observable gravitational waves. For a selection of benchmark points, we checked that the phase transitions completed and calculated the nucleation temperatures. We furthermore present samples that feature strong first-order phase transitions from an extensive scan of the whole parameter space. We highlight common features of our samples, including the fact that the Standard Model like Higgs is often not the lightest Higgs in the model.


title: Dynamic Stale Synchronous Parallel Distributed Training for Deep Learning

abstract: Deep learning is a popular machine learning technique and has been applied to many real-world problems. However, training a deep neural network is very time-consuming, especially on big data. It has become difficult for a single machine to train a large model over large datasets. A popular solution is to distribute and parallelize the training process across multiple machines using the parameter server framework. In this paper, we present a distributed paradigm on the parameter server framework called Dynamic Stale Synchronous Parallel (DSSP) which improves the state-of-the-art Stale Synchronous Parallel (SSP) paradigm by dynamically determining the staleness threshold at the run time. Conventionally to run distributed training in SSP, the user needs to specify a particular staleness threshold as a hyper-parameter. However, a user does not usually know how to set the threshold and thus often finds a threshold value through trial and error, which is time-consuming. Based on workers' recent processing time, our approach DSSP adaptively adjusts the threshold per iteration at running time to reduce the waiting time of faster workers for synchronization of the globally shared parameters, and consequently increases the frequency of parameters updates (increases iteration throughput), which speedups the convergence rate. We compare DSSP with other paradigms such as Bulk Synchronous Parallel (BSP), Asynchronous Parallel (ASP), and SSP by running deep neural networks (DNN) models over GPU clusters in both homogeneous and heterogeneous environments. The results show that in a heterogeneous environment where the cluster consists of mixed models of GPUs, DSSP converges to a higher accuracy much earlier than SSP and BSP and performs similarly to ASP. In a homogeneous distributed cluster, DSSP has more stable and slightly better performance than SSP and ASP, and converges much faster than BSP.


title: Describing correlated observations of neutrino and gamma ray flares from the blazar TXS 0506+056 with proton blazar model

abstract: Recent detection of the neutrino event, IceCube-170922A by IceCube observatory from the Blazar TXS 0506+056 in the state of enhanced gamma ray emission indicates for acceleration of cosmic rays in the blazar jet. The non-detection of the broadline emission in the optical spectrum of TXS 0506+056 and other BL Lac objects suggests that external photons emissions are weak and hence photo-meson (p-gamma) interaction may not be a favored mechanism for high energy neutrino production. The lack of broadline signatures also creates doubt about the presence of a high density cloud in the vicinity of the super-massive black hole (SMBH) of TXS 0506+056 and consequently raised question on hadronuclear (pp) interaction interpretation like relativistic jet meets with high density cloud. Here we demonstrate that non-relativistic protons in the proton blazar model, those come into existence under charge neutrality condition of the blazar jet, offer sufficient target matter for pp-interaction with shock accelerated protons and consequently the model can describe consistently the observed high energy gamma rays and neutrino signal from the blazar TXS 0506+056.


title: MOD: Minimally Ordered Durable Datastructures for Persistent Memory

abstract: Persistent Memory (PM) makes possible recoverable applications that can preserve application progress across system reboots and power failures. Actual recoverability requires careful ordering of cacheline flushes, currently done in two extreme ways. On one hand, expert programmers have reasoned deeply about consistency and durability to create applications centered on a single custom-crafted durable datastructure. On the other hand, less-expert programmers have used software transaction memory (STM) to make atomic one or more updates, albeit at a significant performance cost due largely to ordered log updates. In this work, we propose the middle ground of composable persistent datastructures called Minimally Ordered Durable (MOD) datastructures. MOD is a C++ library of several datastructures---currently, map, set, stack, queue and vector--- that often perform better than STM and yet are relatively easy to use. They allow multiple updates to one or more datastructures to be atomic with respect to failure. Moreover, we provide a recipe to create more recoverable datastructures. MOD is motivated by our analysis of real Intel Optane PM hardware showing that allowing unordered, overlapping flushes significantly improves performance. MOD reduces ordering by adapting existing techniques for out-of-place updates (like shadow paging) with space-reducing structural sharing (from functional programming). MOD exposes a Basic interface for single updates and a Composition interface for atomically performing multiple updates. Relative to the state-of-the-art Intel PMDK v1.5 STM, MOD improves map, set, stack, queue microbenchmark performance by 40%, and speeds up application benchmark performance by 38%.


title: Matrix equation techniques for certain evolutionary partial differential equations

abstract: We show that the discrete operator stemming from the time and space discretization of evolutionary partial differential equations can be represented in terms of a single Sylvester matrix equation. A novel solution strategy that combines projection techniques with the full exploitation of the entry-wise structure of the involved coefficient matrices is proposed. The resulting scheme is able to efficiently solve problems with a tremendous number of degrees of freedom while maintaining a low storage demand as illustrated in several numerical examples.


title: New stable method to solve heat conduction problems in extremely large systems

abstract: We present a new explicit and stable numerical algorithm to solve the homogeneous heat equation. We illustrate the performance of the new method in the cases of two 2D systems with highly inhomogeneous random parameters. Spatial discretization of these problems results in huge and stiff ordinary differential equation systems, which can be solved by our novel method faster than by explicit or the commonly used implicit methods.

categories: cs.CE cs.NA math.NA physics.comp-ph physics.data-an physics.flu-dyn


title: BooVAE: Boosting Approach for Continual Learning of VAE

abstract: Variational autoencoder (VAE) is a deep generative model for unsupervised learning, allowing to encode observations into the meaningful latent space. VAE is prone to catastrophic forgetting when tasks arrive sequentially, and only the data for the current one is available. We address this problem of continual learning for VAEs. It is known that the choice of the prior distribution over the latent space is crucial for VAE in the non-continual setting. We argue that it can also be helpful to avoid catastrophic forgetting. We learn the approximation of the aggregated posterior as a prior for each task. This approximation is parametrised as an additive mixture of distributions induced by encoder evaluated at trainable pseudo-inputs. We use a greedy boosting-like approach with entropy regularisation to learn the components. This method encourages components diversity, which is essential as we aim at memorising the current task with the fewest components possible. Based on the learnable prior, we introduce an end-to-end approach for continual learning of VAEs and provide empirical studies on commonly used benchmarks (MNIST, Fashion MNIST, NotMNIST) and CelebA datasets. For each dataset, the proposed method avoids catastrophic forgetting in a fully automatic way.


title: The influence of impurities on the charge carrier mobility of small molecule organic semiconductors

abstract: Amorphous organic semiconductors based on small molecules and polymers are used in many applications, most prominently organic light emitting diodes (OLEDs) and organic solar cells. Impurities and charge traps are omnipresent in most currently available organic semiconductors and limit charge transport and thus device efficiency. The microscopic cause as well as the chemical nature of these traps are presently not well understood. Using a multiscale model we characterize the influence of impurities on the density of states and charge transport in small-molecule amorphous organic semiconductors. We use the model to quantitatively describe the influence of water molecules and water-oxygen complexes on the electron and hole mobilities. These species are seen to impact the shape of the density of states and to act as explicit charge traps within the energy gap. Our results show that trap states introduced by molecular oxygen can be deep enough to limit the electron mobility in widely used materials.


title: Distinguishing deformation mechanisms in elastocapillary experiments

abstract: Soft materials are known to deform due to a variety of mechanisms, including capillarity, buoyancy, and swelling. In this paper, we present experiments on polyvinylsiloxane gel threads partially-immersed in three liquids with different solubility, wettability, and swellability. Our results demonstrate that deformations due to capillarity, buoyancy, and swelling can be of similar magnitude as such threads come to static equilibrium. To account for all three effects being present in a single system, we derive a model capable of explaining the observed data and use it to determine the force law at the three-phase contact line. The results show that the measured forces are consistent with the expected Young-Dupr\'e equation, and do not require the inclusion of a tangential contact line force.


title: Methods for Measuring Magnetic Flux Crosstalk Between Tunable Transmons

abstract: In the gate model of quantum computing, a program is typically decomposed into a sequence of 1- and 2-qubit gates that are realized as control pulses acting on the system. A key requirement for a scalable control system is that the qubits are addressable - that control pulses act only on the targeted qubits. The presence of control crosstalk makes this addressability requirement difficult to meet. In order to provide metrics that can drive requirements for decreasing crosstalk, we present three measurements that directly quantify the DC and AC flux crosstalk present between tunable transmons, with sensitivities as fine as 0.001%. We develop the theory to connect AC flux crosstalk measures to the infidelity of a parametrically activated two-qubit gate. We employ quantum process tomography in the presence of crosstalk to provide an empirical study of the effects of crosstalk on two-qubit gate fidelity.


title: $O(N^3)$ Measurement Cost for Variational Quantum Eigensolver on Molecular Hamiltonians

abstract: Variational Quantum Eigensolver (VQE) is a promising algorithm for near-term quantum machines. It can be used to estimate the ground state energy of a molecule by performing separate measurements of $O(N^4)$ terms. Several recent papers observed that this scaling may be reducible to $O(N^3)$ by partitioning the terms into linear-sized commuting families that can be measured simultaneously. We confirm these empirical observations by studying the MIN-COMMUTING-PARTITION problem at the level of the fermionic Hamiltonian and its encoding into qubits. Moreover, we provide a fast, pre-computable procedure for creating linearly-sized commuting partitions by solving a round-robin scheduling problem via flow networks.


title: Nash equilibria strategies and equivalent single-objective optimization problems. The case of linear partial differential equations

abstract: In this paper we study the existence and uniqueness of Nash equilibria (solution to competition-wise problems, with several controls trying to reach possibly different goals) associated to linear partial differential equations and show that, in some cases, they are also the solution of suitable single-objective optimization problems (i.e. cooperative-wise problems, where all the controls cooperate to reach a common goal). We use cost functions associated with a particular linear parabolic partial differential equations and distributed controls, but the results are also valid for more general linear differential equations (including elliptic and hyperbolic cases) and controls (e.g. boundary controls, initial value controls,...).


title: KK Higgs produced in association with a top quark pair in the bulk RS Model

abstract: We present a search strategy for the first Kaluza-Klein (KK) mode of the Higgs boson in the framework of the Randall-Sundrum (RS) model with a deformed metric. We study the production of this massive excitation in association with a ttbar pair at the Large Hadron Collider (LHC). The KK Higgs primarily decays into a boosted ttbar final state and we then end up with an interesting four-top final state of which two are boosted. The boosted products in the final state improve the sensitivity for the search of the KK Higgs in this channel whose production cross-section is otherwise rather small. Our results suggest that masses of the KK Higgs resonance upto about 1.2 TeV may be explorable at the highest planned luminosities of the LHC. Beyond this mass, the KK Higgs cross-section is too tiny for it to be explored at the LHC and may be possible only at a future higher energy collider.


title: Adapt or Get Left Behind: Domain Adaptation through BERT Language Model Finetuning for Aspect-Target Sentiment Classification

abstract: Aspect-Target Sentiment Classification (ATSC) is a subtask of Aspect-Based Sentiment Analysis (ABSA), which has many applications e.g. in e-commerce, where data and insights from reviews can be leveraged to create value for businesses and customers. Recently, deep transfer-learning methods have been applied successfully to a myriad of Natural Language Processing (NLP) tasks, including ATSC. Building on top of the prominent BERT language model, we approach ATSC using a two-step procedure: self-supervised domain-specific BERT language model finetuning, followed by supervised task-specific finetuning. Our findings on how to best exploit domain-specific language model finetuning enable us to produce new state-of-the-art performance on the SemEval 2014 Task 4 restaurants dataset. In addition, to explore the real-world robustness of our models, we perform cross-domain evaluation. We show that a cross-domain adapted BERT language model performs significantly better than strong baseline models like vanilla BERT-base and XLNet-base. Finally, we conduct a case study to interpret model prediction errors.


title: The $P{\Phi}$-Compromise Function as a criterion of merit to optimize irreversible thermal engines

abstract: Several authors have proposed out of equilibrium thermal engines models, allowing optimization processes involving a trade off between the power output of the engine and its dissipation. These operating regimes are achieved by using objective functions such as the ecological function ($EF$). In order to measure the quality of the balance between these characteristic functions, it was proposed a relationship where power output and dissipation are evaluated in the above mentioned $EF$-regime and they are compared with respect to its values at the regime of maximum power output. We called this relationship "Compromise Function" and only depends of a parameter that measures the quality of the compromise. Thereafter this function was used to select a value of the mentioned parameter to obtain the generalization of some different objective functions (generalizations of ecological function, omega function and efficient power), by demanding that these generalization parameters maximize the above mentioned functions. In this work we demonstrate that this function can be used directly as an objective function: the "$P{\Phi}$-Compromise Function" ($C_{P\Phi}$), also that the operation modes corresponding to the maximum Generalized Ecological Function, maximum Generalized Omega Function and maximum Efficient power output, are special cases of the operation mode of maximum $C_{P\Phi}$, having the same optimum high reduced temperature, then the characteristic functions will be the same in any of the above three working regimes, independent of the algebraic complexity of each generalized function. These results are presented for two different models of an irreversible energy converter: a non-endoreversible and a totally irreversible, both with heat leakage.


title: Symmetries and Conservation Laws in Quantum Trajectories: Dissipative Freezing

abstract: In driven-dissipative systems, the presence of a strong symmetry guarantees the existence of several steady states belonging to different symmetry sectors. Here we show that, when a system with a strong symmetry is initialized in a quantum superposition involving several of these sectors, each individual stochastic trajectory will randomly select a single one of them and remain there for the rest of the evolution. Since a strong symmetry implies a conservation law for the corresponding symmetry operator on the ensemble level, this selection of a single sector from an initial superposition entails a breakdown of this conservation law at the level of individual realizations. Given that such a superposition is impossible in a classical, stochastic trajectory, this is a a purely quantum effect with no classical analogue. Our results show that a system with a closed Liouvillian gap may exhibit, when monitored over a single run of an experiment, a behaviour completely opposite to the usual notion of dynamical phase coexistence and intermittency, which are typically considered hallmarks of a dissipative phase transition. We discuss our results with a simple, realistic model of squeezed superradiance.


title: Systematic Analysis of Image Generation using GANs

abstract: Generative Adversarial Networks have been crucial in the developments made in unsupervised learning in recent times. Exemplars of image synthesis from text or other images, these networks have shown remarkable improvements over conventional methods in terms of performance. Trained on the adversarial training philosophy, these networks aim to estimate the potential distribution from the real data and then use this as input to generate the synthetic data. Based on this fundamental principle, several frameworks can be generated that are paragon implementations in several real-life applications such as art synthesis, generation of high resolution outputs and synthesis of images from human drawn sketches, to name a few. While theoretically GANs present better results and prove to be an improvement over conventional methods in many factors, the implementation of these frameworks for dedicated applications remains a challenge. This study explores and presents a taxonomy of these frameworks and their use in various image to image synthesis and text to image synthesis applications. The basic GANs, as well as a variety of different niche frameworks, are critically analyzed. The advantages of GANs for image generation over conventional methods as well their disadvantages amongst other frameworks are presented. The future applications of GANs in industries such as healthcare, art and entertainment are also discussed.


title: Throwing away antimatter via neutrino oscillations during the reheating era

abstract: The simplest possibility to explain the baryon asymmetry of the Universe is to assume that radiation is created asymmetrically between baryons and anti-baryons after the inflation. We propose a new mechanism of this kind where CP-violating flavor oscillations of left-handed leptons in the reheating era distribute the lepton asymmetries partially into the right-handed neutrinos while net asymmetry is not created. The asymmetry stored in the right-handed neutrinos is later washed out by the lepton number violating decays, and it ends up with the net lepton asymmetry in the Standard Model particles, which is converted into the baryon asymmetry by the sphaleron process. This scenario works for a range of masses of the right-handed neutrinos while no fine-tuning among the masses is required. The reheating temperature of the Universe can be as low as $O(10)$~TeV if we assume that the decays of inflatons in the perturbative regime are responsible for the reheating. For the case of the reheating via the dissipation effects, the reheating temperature can be as low as $O(100)$~GeV.


title: Network Mechanism for Insect Olfaction

abstract: Early olfactory pathway responses to the presentation of an odor exhibit remarkably similar dynamical behavior across phyla from insects to mammals, and frequently involve transitions among quiescence, collective network oscillations, and asynchronous firing. We hypothesize that the time scales of fast excitation and fast and slow inhibition present in these networks may be the essential element underlying this similar behavior, and design an idealized, conductance-based integrate-and-fire (I&F) model to verify this hypothesis via numerical simulations. To better understand the mathematical structure underlying the common dynamical behavior across species, we derive a firing-rate (FR) model and use it to extract a slow passage through a saddle-node-on-an-invariant-circle (SNIC) bifurcation structure. We expect this bifurcation structure to provide new insights into the understanding of the dynamical behavior of neuronal assemblies and that a similar structure can be found in other sensory systems.


title: Exploring the initial stages in heavy-ion collisions with high-pT RAA and v2 theory and data

abstract: Traditionally, low-pT sector is used to infer the features of initial stages before QGP thermalization. On the other hand, recently acquired wealth of high-pT experimental data paves the way to utilize the high-pT particles energy loss in exploring the initial stages. We here study how four different commonly considered initial-stage scenarios which have the same temperature profile after, but differ in the 'temperature' profile before thermalization affect predictions of high-pT RAA and v2 observables. Contrary to common expectations, we obtain that high-pT v2 is insensitive to the initial stages of medium evolution, being unable to discriminate between different conditions. On the other hand, RAA is sensitive to these conditions, however, within the current errorbars, the sensitivity is not sufficient to distinguish between different initial stages. Moreover, we also reconsider the validity of widely-used procedure of fitting the energy loss parameters, individually for different initial-stage cases, to reproduce the experimentally observed RAA. We here find that previously reported sensitivity of v2 to different initial states is mainly a consequence of the RAA fitting procedure, which may lead to incorrect conclusions. On the other hand, if a global property, in particular the same average temperature, is imposed to tested temperature profiles, high sensitivity of high-pT v2 is again obtained. We however show that this sensitivity would not be a consequence of differences in initial, but rather final, stages. Consequently, the simultaneous study of high-pT RAA and v2, with consistent energy loss parametrization and stringently controlled temperature profiles, is necessary to assess sensitivity of different variables to differences in initial and final stages.


title: Modular Invariant Models of Lepton Masses at Levels 4 and 5

abstract: We explore alternative descriptions of the charged lepton sector in modular invariant models of lepton masses and mixing angles. In addition to the modulus, the symmetry breaking sector of our models includes ordinary flavons. Neutrino mass terms depend only on the modulus and are tailored to minimize the number of free parameters. The charged lepton Yukawa couplings rely upon the flavons alone. We build modular invariant models at levels 4 and 5, where neutrino masses are described both in terms of the Weinberg operator or through a type I seesaw mechanism. At level 4, our models reproduce the hierarchy among electron, muon and tau masses by letting the weights play the role of Froggatt-Nielsen charges. At level 5, our setup allows the treatment of left and right handed charged leptons on the same footing. We have optimized the free parameters of our models in order to match the experimental data, obtaining a good degree of compatibility and predictions for the absolute neutrino masses and the $CP$ violating phases. At a more fundamental level, the whole lepton sector could be correctly described by the simultaneous presence of several moduli. Our examples are meant to make a first step in this direction.


title: Particle re-acceleration and Faraday-complex structures in the RXC J1314.4-2515 galaxy cluster

abstract: Radio relics are sites of electron (re)acceleration in merging galaxy clusters but the mechanism of acceleration and the topology of the magnetic field in and near relics are yet to be understood. We are carrying out an observational campaign on double relic galaxy clusters starting with RXC J1314.4-2515. With $Jansky Very Large Array$ multi-configuration observations in the frequency range 1-4 GHz, we perform both spectral and polarization analyses, using the Rotation Measure synthesis technique. We use archival $XMM-Newton$ observations to constrain the properties of the shocked region. We discover a possible connection between the activity of a radio galaxy and the emission of the eastern radio relic. In the northern elongated arc of the western radio relic, we detect polarized emission with an average polarization fraction of $31 \ \%$ at 3 GHz and we derive the Mach number of the underlying X-ray shock. Our observations reveal low levels of fractional polarization and Faraday-complex structures in the southern region of the relic, which point to the presence of thermal gas and filamentary magnetic field morphology inside the radio emitting volume. We measured largely different Rotation Measure dispersion from the two relics. Finally, we use cosmological magneto-hydrodynamical simulations to constrain the magnetic field, viewing angle, and to derive the acceleration efficiency of the shock. We find that the polarization properties of RXC J1314.4-2515 are consistent with a radio relic observed at $70^{\circ}$ with respect to the line of sight and that efficient re-acceleration of fossil electrons has taken place.


title: LHC and DIS experimental data in the CT18(Z) global QCD analysis

abstract: We discuss implementation of the LHC experimental data sets in the new CT18 global analysis of quantum chromodynamics (QCD) at the next-to-next-leading order of the QCD coupling strength. New methodological developments in the fitting methodology are discussed. Behavior of the CT18 NNLO PDFs for the conventional and "saturation-inspired" factorization scales in deep-inelastic scattering is reviewed. Four new families of (N)NLO CTEQ-TEA PDFs are presented: CT18, A, X, and Z.


title: Minimum $L^q$-distance estimators for non-normalized parametric models

abstract: We propose and investigate a new estimation method for the parameters of models consisting of smooth density functions on the positive half axis. The procedure is based on a recently introduced characterization result for the respective probability distributions, and is to be classified as a minimum distance estimator, incorporating as a distance function the $L^q$-norm. Throughout, we deal rigorously with issues of existence and measurability of these implicitly defined estimators. Moreover, we provide consistency results in a common asymptotic setting, and compare our new method with classical estimators for the exponential-, the Rayleigh-, and the Burr Type XII distribution in Monte Carlo simulation studies. We also assess the performance of different estimators for non-normalized models in the context of an exponential-polynomial family.


title: Properties of the Circumgalactic Medium in Cosmic Ray-Dominated Galaxy Halos

abstract: We investigate the impact of cosmic rays (CRs) on the circumgalactic medium (CGM) in FIRE-2 simulations, for ultra-faint dwarf through Milky Way (MW)-mass halos hosting star-forming (SF) galaxies. Our CR treatment includes injection by supernovae, anisotropic streaming and diffusion along magnetic field lines, collisional and streaming losses, with constant parallel diffusivity $\kappa\sim3\times10^{29}\,\mathrm{cm^2\ s^{-1}}$ chosen to match $\gamma$-ray observations. With this, CRs become more important at larger halo masses and lower redshifts, and dominate the pressure in the CGM in MW-mass halos at $z\lesssim 1-2$. The gas in these "CR-dominated" halos differs significantly from runs without CRs: the gas is primarily cool (a few $\sim10^{4}\,$K), and the cool phase is volume-filling and has a thermal pressure below that needed for virial or local thermal pressure balance. Ionization of the "low" and "mid" ions in this diffuse cool gas is dominated by photo-ionization, with O VI columns $\gtrsim 10^{14.5}\,\mathrm{cm^{-2}}$ at distances $\gtrsim 150\,\mathrm{kpc}$. CR and thermal gas pressure are locally anti-correlated, maintaining total pressure balance, and the CGM gas density profile is determined by the balance of CR pressure gradients and gravity. Neglecting CRs, the same halos are primarily warm/hot ($T\gtrsim 10^{5}\,$K) with thermal pressure balancing gravity, collisional ionization dominates, O VI columns are lower and Ne VIII higher, and the cool phase is confined to dense filaments in local thermal pressure equilibrium with the hot phase.


title: Periodic Fast Radio Bursts from Young Neutron Stars

abstract: Fast radio bursts (FRBs) are highly energetic radio pulses from cosmological origins. Despite an abundance of detections, their nature remains elusive. At least a subset of FRBs is expected to repeat, as the daily FRB rate surpasses that of any known cataclysmic event, which has been confirmed by observations. One of the proposed mechanisms to generate repeating FRBs is supergiant pulses from young and highly spinning NSs, in which case FRBs could inherit the periodicity of their parent NS. Here we examine the consequences of such a population of periodic fast radio bursts (PFRBs). We calculate the rate and lifetime of PFRB progenitors, and find that each newly born highly spinning NS has to emit a number $N_{\rm PFRB}\sim 10^2$ of bursts during its active lifetime of $\tau\sim 100$ years, after which it becomes too dim and crosses a PFRB "death line" analogous to the pulsar one. We propose several tests of this hypothesis. First, the period of PFRBs would increase over time, and their luminosity would decrease, due to the NS spin-down. Second, PFRBs may show modest amounts of rotation measure, given the lack of expelled matter from the pulsar, as opposed to the magnetar-sourced FRBs proposed to explain the first repeater FRB 121102. As an example, we study whether the second confirmed repeater (FRB 180814) is a PFRB, given the preference for an inter-pulse separation of 13 ms within its sub-bursts. We show that, if confirmed, this period would place FRB 180814 in a different category as FRB 121102. We develop tests that would identify---and characterize---the prospective population of PFRBs.


title: Constraining axion-like particles from rare pion decays

abstract: Ultraviolet completions for axion-like particles (ALPs) lighter than the neutral pion generically induce ALP-neutral pion mixing, and are therefore sensitive to direct constraints on the mixing angle. For ALPs below the pion mass, we demonstrate that strong and novel bounds on the ALP-pion mixing angle can be extracted from existing rare pion decay data, measured by the PIENU and PIBETA experiments.


title: Elemental Abundances in M31: A Comparative Analysis of Iron and Alpha Element Abundances in the Outer Disk, Giant Stellar Stream, and Inner Halo of M31

abstract: We measured [Fe/H] and [$\alpha$/Fe] using spectral synthesis of low-resolution stellar spectroscopy for 70 individual red giant branch stars across four fields spanning the outer disk, Giant Stellar Stream (GSS), and inner halo of M31. Fields at M31-centric projected distances of 23 kpc in the halo, 12 kpc in the halo, 22 kpc in the GSS, and 26 kpc in the outer disk are $\alpha$-enhanced, with $\langle$[$\alpha$/Fe]$\rangle$ = 0.43, 0.50, 0.41, and 0.58, respectively. The 23 kpc and 12 kpc halo fields are relatively metal-poor, with $\langle$[Fe/H]$\rangle$ = $-$1.54 and $-$1.30, whereas the 22 kpc GSS and 26 kpc outer disk fields are relatively metal-rich with $\langle$[Fe/H]$\rangle$ = $-$0.84 and $-$0.92, respectively. For fields with substructure, we separated the stellar populations into kinematically hot stellar halo components and kinematically cold components. We did not find any evidence of an [$\alpha$/Fe] gradient along the high surface brightness core of the GSS between $\sim$17$-$22 kpc. However, we found tentative suggestions of a negative [$\alpha$/Fe] gradient in the stellar halo, which may indicate that different progenitor(s) or formation mechanisms contributed to the build up of the inner versus outer halo. Additionally, the [$\alpha$/Fe] distribution of the metal-rich ([Fe/H] $>$ $-$1.5), smooth inner stellar halo (r$_{\rm{proj}}$ $\lesssim$ 26 kpc) is inconsistent with having formed from the disruption of progenitor(s) similar to present-day M31 satellite galaxies. The 26 kpc outer disk is most likely associated with the extended disk of M31, where its high $\alpha$-enhancement provides support for an episode of rapid star formation in M31's disk, possibly induced by a major merger.


title: Probing the $CP$ nature of the top quark Yukawa at hadron colliders

abstract: We analyze the prospects of probing the $CP$-odd $i \tilde \kappa \bar t \gamma^5 t h$ interaction at the LHC and its projected upgrades, the high-luminosity and high-energy LHC, directly using associated on-shell Higgs boson and top quark or top quark pair production. To this end we first construct a $CP$-odd observable based on top quark polarization in $W b \to t h$ scattering with optimal linear sensitivity to $\tilde \kappa$. For the corresponding hadronic process $pp \to t h j$ we present a method of extracting the phase-space dependent weight function that allows to retain close to optimal sensitivity to $\tilde{\kappa}$. We project future sensitivity to the signal in $pp \to t(\to \ell \nu b )h(\to b\bar b) j$. Based on these insights we propose novel $CP$-odd observables for top quark pair production in association with the Higgs, $pp \to t \bar t h$, with semileptonically decaying tops and $h\to \bar b b$, that rely solely on measuring the momenta of leptons and $b$-jets from the decaying tops without having to distinguish the charge of the $b$-jets. Among the many possibilities we single out an observable that can potentially probe $\tilde \kappa \sim 0.5$ at the high-luminosity LHC and $\tilde \kappa \sim 0.1$ at high-energy LHC with $2\sigma$ confidence.


title: Simulating the Complexity of the Dark Matter Sheet I: Numerical Algorithms

abstract: At early times dark matter has a thermal velocity dispersion of unknown amplitude which, for warm dark matter models, can influence the formation of nonlinear structure on observable scales. We propose a new scheme to simulate cosmologies with a small-scale suppression of perturbations that combines two previous methods in a way that avoids the numerical artefacts which have so far prevented either from producing fully reliable results. At low densities and throughout most of the cosmological volume, we represent the dark matter phase-sheet directly using high-accuracy interpolation, thereby avoiding the artificial fragmentation which afflicts particle-based methods in this regime. Such phase-sheet methods are, however, unable to follow the rapidly increasing complexity of the denser regions of dark matter haloes, so for these we switch to an N-body scheme which uses the geodesic deviation equation to track phase-sheet properties local to each particle. In addition, we present a novel high-resolution force calculation scheme based on an oct-tree of cubic force resolution elements which is well suited to approximate the force-field of our combined sheet+particle distribution. Our hybrid simulation scheme enables the first reliable simulations of the internal structure of low-mass haloes in a warm dark matter cosmology.


title: Circumnavigating Collinear Superspace

abstract: In this paper, we extend the collinear superspace formalism to include the full range of $\mathcal{N} = 1$ supersymmetric interactions. Building on the effective field theory rules developed in a companion paper - "Navigating Collinear Superspace" - we construct collinear superspace Lagrangians for theories with non-trivial $F$- and $D$-term auxiliary fields. For (massless) Wess-Zumino models, the key ingredient is a novel type of Grassmann-valued supermultiplet whose lowest component is a (non-propagating) fermionic degree of freedom. For gauge theories coupled to charged chiral matter, the key ingredient is a novel type of vector superfield whose lowest component is a non-propagating gauge potential. This unique vector superfield is used to construct a gauge-covariant derivative; while such an object does not appear in the standard full superspace formalism, it is crucial for modeling gauge interactions when the theory is expressed on a collinear slice. This brings us full circle, by showing that all types of $\mathcal{N} = 1$ theories in four dimensions can be constructed in collinear superspace from purely infrared considerations. We speculate that supersymmetric theories with $\mathcal{N} > 1$ could also be implemented using similar collinear superspace constructions.


title: The Multi-Tracer Optimal Estimator applied to VIPERS

abstract: We use mock galaxy data from the VIMOS Public Extragalactic Redshift Survey (VIPERS) to test the performance of the Multi-Tracer Optimal Estimator (MTOE) of Abramo et al. as a tool to measure the monopoles of the power spectra of multiple tracers of the large-scale structure, $P^{(0)}_\alpha(k)$. We show that MTOE provides more accurate measurements than the standard technique of Feldman, Kaiser & Peacock (FKP), independently of the tracer-selection strategy adopted, on both small and large scales. The largest improvements on individual $P^{(0)}_\alpha(k)$ are obtained using a colour-magnitude selection on small scales, due to MTOE being naturally better equipped to deal with shot noise: we report an average error reduction with respect to FKP of $\sim$ 40$\%$ at $k \, [h$ Mpc$^{-1}]\gtrsim 0.3$. On large scales ($k[h$ Mpc$^{-1}]\lesssim0.1$), the gain in accuracy resulting from cosmic-variance cancellation is $\sim$ 10$\%$ for the ratios of $P^{(0)}_\alpha(k)$. We have carried out a Monte-Carlo Markov Chain analysis to determine the impact of these gains on several quantities derived from $P^{(0)}_\alpha(k)$. If we push the measurement to scales $0.3 < k \, [h$ Mpc$^{-1}]< 0.5$, the average improvements are $\sim$ 30 $\%$ for the amplitudes of the monopoles, $\sim$ 70 $\%$ for the monopole ratios, and $\sim$ 20 $\%$ for the galaxy biases. Our results highlight the potential of MTOE to shed light upon the physics that operate both on large and small cosmological scales. The effect of MTOE on cosmological constraints using VIPERS data will be addressed in a separate paper.


title: HETDEX Pilot Survey. VI. [OIII] Emitters and Expectations for a Local Sample of Star Forming Galaxies in HETDEX

abstract: We assemble an unbiased sample of 29 galaxies with [O II] $\lambda 3727$ and/or [O III] $\lambda 5007$ detections at $z < 0.15$ from the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) Pilot Survey (HPS). HPS finds galaxies without pre-selection based on their detected emission lines via integral field spectroscopy. Sixteen of these objects were followed up with the second-generation, low resolution spectrograph (LRS2) on the upgraded Hobby-Eberly Telescope. Oxygen abundances were then derived via strong emission lines using a Bayesian approach. We find most of the galaxies fall along the mass-metallicity relation derived from photometrically selected star forming galaxies in the Sloan Digital Sky Survey (SDSS). However, two of these galaxies have low metallicity (similar to the very rare green pea galaxies in mass-metallicity space). The star formation rates of this sample fall in an intermediate space between the SDSS star forming main sequence and the extreme green pea galaxies. We conclude that spectroscopic selection fills part of the mass-metallicity-SFR phase space that is missed in photometric surveys with pre-selection like SDSS, i.e., we find galaxies that are actively forming stars but are faint in continuum. We use the results of this pilot investigation to make predictions for the upcoming unbiased, large spectroscopic sample of local line emitters from HETDEX. With the larger HETDEX survey we will determine if galaxies selected spectroscopically without continuum brightness pre-selection have metallicities that fall on a continuum that bridges typical star forming and rarer, more extreme systems like green peas.


title: Brain Energetics, Mitochondria, and Traumatic Brain Injury

abstract: We review current thinking about, and draw connections between, brain energetics and metabolism, mitochondria and traumatic brain injury. In addition to summarizing current thinking in these disciplines, our goal is to suggest a framework for mechanisms and pathways based on optimal energetic decisions.


title: Low-$Q^2$ empirical parametrizations of the $N^\ast$ helicity amplitudes

abstract: The data associated with the electromagnetic excitations of the nucleon ($\gamma^\ast N \to N^\ast$) are usually parametrized by helicity amplitudes at the resonance $N^\ast$ rest frame. The properties of the $\gamma^\ast N \to N^\ast$ transition current at low $Q^2$ can be, however, better understood when expressed in terms of structure form factors, particularly near the pseudothreshold, when the magnitude of the photon three-momentum vanishes ($|{\bf q}| =0$). At the pseudothreshold the invariant four-momentum square became $q^2= (M_R-M_N)^2$, well in the timelike region $Q^2 =- q^2 < 0$ [$M_N$ and $M_R$ are the mass of the nucleon and of the resonance, respectively]. In the helicity amplitude representation, the amplitudes have well-defined dependences on $|{\bf q}|$, near the pseudothreshold, and there are correlations between different amplitudes. Those constraints are often ignored in the empirical parametrizations of the helicity amplitudes. In the present work, we show that the structure of the transition current near the pseudothreshold has an impact on the parametrizations of the data. We present a method which modifies analytic parametrizations of the data at low $Q^2$, in order to take into account the constraints of the transition amplitudes near the pseudothreshold. The model dependence of the parametrizations on the low-$Q^2$ data is studied in detail.


title: Sensor Switching Control Under Attacks Detectable by Finite Sample Dynamic Watermarking Tests

abstract: Control system security is enhanced by the ability to detect malicious attacks on sensor measurements. Dynamic watermarking can detect such attacks on linear time-invariant (LTI) systems. However, existing theory focuses on attack detection and not on the use of watermarking in conjunction with attack mitigation strategies. In this paper, we study the problem of switching between two sets of sensors: One set of sensors has high accuracy but is vulnerable to attack, while the second set of sensors has low accuracy but cannot be attacked. The problem is to design a sensor switching strategy based on attack detection by dynamic watermarking. This requires new theory because existing results are not adequate to control or bound the behavior of sensor switching strategies that use finite data. To overcome this, we develop new finite sample hypothesis tests for dynamic watermarking in the case of bounded disturbances, using the modern theory of concentration of measure for random matrices. Our resulting switching strategy is validated with a simulation analysis in an autonomous driving setting, which demonstrates the strong performance of our proposed policy.


title: Adaptively Sparse Transformers

abstract: Attention mechanisms have become ubiquitous in NLP. Recent architectures, notably the Transformer, learn powerful context-aware word representations through layered, multi-headed attention. The multiple heads learn diverse types of word relationships. However, with standard softmax attention, all attention heads are dense, assigning a non-zero weight to all context words. In this work, we introduce the adaptively sparse Transformer, wherein attention heads have flexible, context-dependent sparsity patterns. This sparsity is accomplished by replacing softmax with $\alpha$-entmax: a differentiable generalization of softmax that allows low-scoring words to receive precisely zero weight. Moreover, we derive a method to automatically learn the $\alpha$ parameter -- which controls the shape and sparsity of $\alpha$-entmax -- allowing attention heads to choose between focused or spread-out behavior. Our adaptively sparse Transformer improves interpretability and head diversity when compared to softmax Transformers on machine translation datasets. Findings of the quantitative and qualitative analysis of our approach include that heads in different layers learn different sparsity preferences and tend to be more diverse in their attention distributions than softmax Transformers. Furthermore, at no cost in accuracy, sparsity in attention heads helps to uncover different head specializations.


title: Numerical analysis of a semilinear fractional diffusion equation

abstract: This paper considers the numerical analysis of a semilinear fractional diffusion equation with nonsmooth initial data. A new Gr\"onwall's inequality and its discrete version are proposed. By the two inequalities, error estimates in three Sobolev norms are derived for a spatial semi-discretization and a full discretization, which are optimal with respect to the regularity of the solution. A sharp temporal error estimate on graded temporal grids is also rigorously established. In addition, the spatial accuracy $\scriptstyle O(h^2(t^{-\alpha} + \ln(1/h)\!)\!) $ in the pointwise $ \scriptstyle L^2(\Omega) $-norm is obtained for a spatial semi-discretization. Finally, several numerical results are provided to verify the theoretical results.


title: Inverse Radon transform and the transverse-momentum dependent functions

abstract: We revisit the standard representation of the (inverse) Radon transform which is well-known in the mathematical literature. We extend this representation to the case involving the parton distributions. We have found the new additional contribution which is essentially related to the generalized transverse-momentum dependent parton distribution and double-distribution functions. We discuss the possible relationship of this term with the Sivers function.


title: Singularity theorems and the inclusion of torsion in affine theories of gravity

abstract: We extend the scope of the Raychaudhuri-Komar singularity theorem of General Relativity to affine theories of gravity with and without torsion. We first generalize the existing focusing theorems using time-like and null congruences of curves which are hypersurface orthogonal, showing how the presence of torsion affects the formation of focal points in Lorentzian manifolds. Considering the energy conservation on a given affine gravity theory, we prove new singularity theorems for accelerated curves in the cases of Lorentzian manifolds containing perfect fluids or scalar field matter sources.


title: Enumeration and Extensions of Word-representants

abstract: Given a finite word $w$ over a finite alphabet $V$, consider the graph with vertex set $V$ and with an edge between two elements of $V$ if and only if the two elements alternate in the word $w$. Such a graph is said to be word-representable or 11-representable by the word $w$; this latter terminology arises from the phenomenon that the condition of two elements $x$ and $y$ alternating in a word $w$ is the same as the condition of the subword of $w$ induced by $x$ and $y$ avoiding the pattern 11. In this paper, we first study minimal length words which word-represent graphs, giving an explicit formula for both the length and the number of such words in the case of trees and cycles. We then extend the notion of word-representability (or 11-representability) of graphs to $t$-representability of graphs, for any pattern $t$ on two letters. We prove that every graph is $t$-representable for any pattern $t$ on two letters (except for possibly one class of $t$). Finally, we pose a few open problems for future consideration.


title: Lasing in $\mathrm{Si}_3\mathrm{N}_4$-organic hybrid (SiNOH) waveguides

abstract: Lasers that operate in the visible wavelength range and can be easily integrated on a chip together with photonic integrated circuits are of particular interest for biophotonic applications. We present a new class of laser sources that emit at 600 nm wavelength and can be efficiently integrated on the $\mathrm{Si}_3\mathrm{N}_4$ platform using low-cost fabrication techniques. The $\mathrm{Si}_3\mathrm{N}_4$-organic hybrid lasers are optically pumped from top without the need of a precise alignment. We investigate different laser geometries based on spiral shaped ring resonators and distributed feedback resonators. We discuss the laser properties with a focus on their application in the field of bioscience.


title: Approximating Stacked and Bidirectional Recurrent Architectures with the Delayed Recurrent Neural Network

abstract: Recent work has shown that topological enhancements to recurrent neural networks (RNNs) can increase their expressiveness and representational capacity. Two popular enhancements are stacked RNNs, which increases the capacity for learning non-linear functions, and bidirectional processing, which exploits acausal information in a sequence. In this work, we explore the delayed-RNN, which is a single-layer RNN that has a delay between the input and output. We prove that a weight-constrained version of the delayed-RNN is equivalent to a stacked-RNN. We also show that the delay gives rise to partial acausality, much like bidirectional networks. Synthetic experiments confirm that the delayed-RNN can mimic bidirectional networks, solving some acausal tasks similarly, and outperforming them in others. Moreover, we show similar performance to bidirectional networks in a real-world natural language processing task. These results suggest that delayed-RNNs can approximate topologies including stacked RNNs, bidirectional RNNs, and stacked bidirectional RNNs - but with equivalent or faster runtimes for the delayed-RNNs.


title: Photonic chip-based resonant supercontinuum

abstract: Supercontinuum generation in optical fibers is one of the most dramatic nonlinear effects discovered, allowing short pulses to be converted into multi-octave spanning coherent spectra. However, generating supercontinua that are both coherent and broadband requires pulses that are simultaneously ultrashort with high peak power. This results in a reducing efficiency with increasing pulse repetition rate, that has hindered supercontinua at microwave line spacing, i.e. 10s of GHz. Soliton microcombs by contrast, can generate octave-spanning spectra, but with good conversion efficiency only at vastly higher repetition rates in the 100s of GHz. Here, we bridge this efficiency gap with resonant supercontinuum, allowing supercontinuum generation using input pulses with an ultra-low 6 picojoule energy, and duration of 1 picosecond, 10-fold longer than what is typical. By applying synchronous pulse-driving to a dispersion-engineered, low-loss Si$_3$N$_4$ photonic chip microresonator, we generate dissipative Kerr solitons with a strong dispersive wave, both bound to the input pulse. This creates a smooth, flattened 2,200 line frequency comb, with an electronically detectable repetition rate of 28 GHz, constituting the largest bandwidth-line-count product for any microcomb generated to date. Strikingly, we observe that solitons exist in a weakly bound state with the input pulse, stabilizing their repetition rate, but simultaneously allowing noise transfer from one to the other to be suppressed even for offset frequencies 100 times lower than the linear cavity decay rate. We demonstrate that this nonlinear filtering can be enhanced by pulse-driving asynchronously, in order to preserve the coherence of the comb. Taken together, our work establishes resonant supercontinuum as a promising route to broadband and coherent spectra.


title: High-resolution 3D refractive index microscopy of multiple-scattering samples from intensity images

abstract: Optical diffraction tomography (ODT) reconstructs a samples volumetric refractive index (RI) to create high-contrast, quantitative 3D visualizations of biological samples. However, standard implementations of ODT use interferometric systems, and so are sensitive to phase instabilities, complex mechanical design, and coherent noise. Furthermore, their reconstruction framework is typically limited to weakly-scattering samples, and thus excludes a whole class of multiple-scattering samples. Here, we implement a new 3D RI microscopy technique that utilizes a computational multi-slice beam propagation method to invert the optical scattering process and reconstruct high-resolution (NA>1.0) 3D RI distributions of multiple-scattering samples. The method acquires intensity-only measurements from different illumination angles, and then solves a non-linear optimization problem to recover the sample 3D RI distribution. We experimentally demonstrate reconstruction of samples with varying amounts of multiple scattering: a 3T3 fibroblast cell, a cluster of C. elegans embryos, and a whole C. elegans worm, with lateral and axial resolutions of 250 nm and 900 nm, respectively.


title: Racial Disparities in Voting Wait Times: Evidence from Smartphone Data

abstract: Equal access to voting is a core feature of democratic government. Using data from millions of smartphone users, we quantify a racial disparity in voting wait times across a nationwide sample of polling places during the 2016 U.S. presidential election. Relative to entirely-white neighborhoods, residents of entirely-black neighborhoods waited 29% longer to vote and were 74% more likely to spend more than 30 minutes at their polling place. This disparity holds when comparing predominantly white and black polling places within the same states and counties, and survives numerous robustness and placebo tests. We shed light on the mechanism for these results and discuss how geospatial data can be an effective tool to both measure and monitor these disparities going forward.


title: Meta-Learning with Warped Gradient Descent

abstract: Learning an efficient update rule from data that promotes rapid learning of new tasks from the same distribution remains an open problem in meta-learning. Typically, previous works have approached this issue either by attempting to train a neural network that directly produces updates or by attempting to learn better initialisations or scaling factors for a gradient-based update rule. Both of these approaches pose challenges. On one hand, directly producing an update forgoes a useful inductive bias and can easily lead to non-converging behaviour. On the other hand, approaches that try to control a gradient-based update rule typically resort to computing gradients through the learning process to obtain their meta-gradients, leading to methods that can not scale beyond few-shot task adaptation. In this work, we propose Warped Gradient Descent (WarpGrad), a method that intersects these approaches to mitigate their limitations. WarpGrad meta-learns an efficiently parameterised preconditioning matrix that facilitates gradient descent across the task distribution. Preconditioning arises by interleaving non-linear layers, referred to as warp-layers, between the layers of a task-learner. Warp-layers are meta-learned without backpropagating through the task training process in a manner similar to methods that learn to directly produce updates. WarpGrad is computationally efficient, easy to implement, and can scale to arbitrarily large meta-learning problems. We provide a geometrical interpretation of the approach and evaluate its effectiveness in a variety of settings, including few-shot, standard supervised, continual and reinforcement learning.


title: On a property of harmonic measure on simply connected domains

abstract: Let $D \subset \mathbb{C}$ be a domain with $0 \in D$. For $R>0$, let ${{\hat \omega }_D}\left( {R} \right)$ denote the harmonic measure of $ D \cap \left\{ {\left| z \right| = R} \right\}$ at $0$ with respect to the domain $ D \cap \left\{ {\left| z \right| < R} \right\} $ and ${\omega_D}\left( {R} \right)$ denote the harmonic measure of $\partial D \cap \left\{ {\left| z \right| \ge R} \right\}$ at $0$ with respect to $D$. The behavior of the functions ${\omega_D}$ and ${{\hat \omega }_D}$ near $\infty$ determines (in some sense) how large $D$ is. However, it is not known whether the functions ${\omega_D}$ and ${{\hat \omega }_D}$ always have the same behavior when $R$ tends to $\infty$. Obviously, ${\omega_D}\left( {R} \right) \le {{\hat \omega }_D}\left( {R} \right)$ for every $R>0$. Thus, the arising question, first posed by Betsakos, is the following: Does there exist a positive constant $C$ such that for all simply connected domains $D$ with $0 \in D$ and all $R>0$, \[{\omega_D}\left( {R} \right) \ge C{{\hat \omega }_D}\left( {R} \right)? \] In general, we prove that the answer is negative by means of two different counter-examples. However, under additional assumptions involving the geometry of $D$, we prove that the answer is positive. We also find the value of the optimal constant for starlike domains.


title: The z=7.54 Quasar ULAS J1342+0928 Is Hosted by a Galaxy Merger

abstract: We present ALMA 0.28"x0.20" (1.4kpc x 1.0kpc) resolution observations of the [CII] 158um line and underlying dust continuum in the host galaxy of the most distant quasar currently known, ULAS J1342+0928 at z=7.5413. Both the [CII] and continuum are detected and spatially resolved. The [CII] line is ~1.5 times more extended than the continuum emission, showing an elongated and complex structure of approximately 3.2kpc x 6.4kpc. Two separate peaks are clearly seen (>6sigma each) in three 100 km/s width [CII] channel maps. The velocity field of the [CII] gas does not show evidence of a coherent rotation field but rather chaotic motion, reminiscent of an on-going merger.


title: Discontinuous Galerkin discretization for quantum simulation of chemistry

abstract: Methods for electronic structure based on Gaussian and molecular orbital discretizations offer a well established, compact representation that forms much of the foundation of correlated quantum chemistry calculations on both classical and quantum computers. Despite their ability to describe essential physics with relatively few basis functions, these representations can suffer from a quartic growth of the number of integrals. Recent results have shown that, for some quantum and classical algorithms, moving to representations with diagonal two-body operators can result in dramatically lower asymptotic costs, even if the number of functions required increases significantly. We introduce a way to interpolate between the two regimes in a systematic and controllable manner, such that the number of functions is minimized while maintaining a block diagonal structure of the two-body operator and desirable properties of an original, primitive basis. Techniques are analyzed for leveraging the structure of this new representation on quantum computers. Empirical results for hydrogen chains suggest a scaling improvement from $O(N^{4.5})$ in molecular orbital representations to $O(N^{2.6})$ in our representation for quantum evolution in a fault-tolerant setting, and exhibit a constant factor crossover at 15 to 20 atoms. Moreover, we test these methods using modern density matrix renormalization group methods classically, and achieve excellent accuracy with respect to the complete basis set limit with a speedup of 1-2 orders of magnitude with respect to using the primitive or Gaussian basis sets alone. These results suggest our representation provides significant cost reductions while maintaining accuracy relative to molecular orbital or strictly diagonal approaches for modest-sized systems in both classical and quantum computation for correlated systems.


title: Physical modelling of galaxy clusters and Bayesian inference in astrophysics

abstract: I compare the mass values obtained with data taken from the Arcminute Microkelvin Imager (AMI) radio interferometer system and from the Planck satellite. The former of these uses a Bayesian analysis pipeline that parameterises a cluster in terms of its physical quantities, and models the dark matter \& baryonic components of a cluster using Navarro-Frenk-White (NFW) and generalised-NFW profiles respectively. I also analyse simulated AMI data with input values based on PwS mass estimates. I then compare three cluster models using AMI data for the 54 cluster sample. The two observational models considered only model the gas content of the cluster. To compare the physical and observational models I consider their posterior parameter estimates, including the calculation of a metric defined between two probability distributions. The models' fit to the cluster data is evaluated by looking at the Bayesian evidence values. Improvements to the physical modelling of galaxy clusters are then considered, either by relaxing some of the assumptions underlying the physical model, or by introducing a new profile for the dark matter component of clusters. The final part of the cluster analysis work focuses on Bayesian analysis using a joint likelihood function of data from both AMI and the Planck satellite simultaneously. Finally, a new Bayesian inference algorithm based on nested sampling is presented. The algorithm, named the "geometric nested sampler", is an adaption of the Metropolis-Hastings nested sampler and makes use of the geometrical interpretation of sets of parameters to sample from their domains efficiently. The geometric nested sampler is tested on several toy models as well as a model representing the emission of gravitational waves from binary black hole mergers.


title: Ramsey Goodness of Paths in Random Graphs

abstract: We say that a graph $G$ is Ramsey for $H_1$ versus $H_2$, and write $G \to (H_1,H_2)$, if every red-blue colouring of the edges of $G$ contains either a red copy of $H_1$ or a blue copy of $H_2$. In this paper we study the threshold for the event that the Erd\H{o}s--R\'enyi random graph $G(N,p)$ is Ramsey for a clique versus a path. We show that $$G\big( (1 + \varepsilon) rn,p \big) \to (K_{r+1},P_n)$$ with high probability if $p \gg n^{-2 / (r + 1)}$, and $$G\big( rn + t, p \big) \to (K_{r+1},P_n)$$ with high probability if $p \gg n^{-2 / (r + 2)}$ and $t \gg 1/p$. Both of these results are sharp (in different ways), since with high probability $G(Cn,p) \not\to (K_{r+1}, P_n)$ for any constant $C > 0$ if $p \ll n^{-2/(r + 1)}$, and $G(rn + t, p) \not\to (K_{r+1}, P_n)$ if $t \ll 1/p$, for any $0 < p \le 1$.


title: Interactive Task and Concept Learning from Natural Language Instructions and GUI Demonstrations

abstract: Natural language programming is a promising approach to enable end users to instruct new tasks for intelligent agents. However, our formative study found that end users would often use unclear, ambiguous or vague concepts when naturally instructing tasks in natural language, especially when specifying conditionals. Existing systems have limited support for letting the user teach agents new concepts or explaining unclear concepts. In this paper, we describe a new multi-modal domain-independent approach that combines natural language programming and programming-by-demonstration to allow users to first naturally describe tasks and associated conditions at a high level, and then collaborate with the agent to recursively resolve any ambiguities or vagueness through conversations and demonstrations. Users can also define new procedures and concepts by demonstrating and referring to contents within GUIs of existing mobile apps. We demonstrate this approach in PUMICE, an end-user programmable agent that implements this approach. A lab study with 10 users showed its usability.


title: Animated Stickies: Fast Video Projection Mapping onto a Markerless Plane through a Direct Closed-Loop Alignment

abstract: This paper presents a fast projection mapping method for moving image content projected onto a markerless planar surface using a low-latency Digital Micromirror Device (DMD) projector. By adopting a closed-loop alignment approach, in which not only the surface texture but also the projected image is tracked by a camera, the proposed method is free from a calibration or position adjustment between the camera and projector. We designed fiducial patterns to be inserted into a fast flapping sequence of binary frames of the DMD projector, which allows the simultaneous tracking of the surface texture and a fiducial geometry separate from a single image captured by the camera. The proposed method implemented on a CPU runs at 400 fps and enables arbitrary video contents to be "stuck" onto a variety of textured surfaces.


title: Hyperbolic String Vertices

abstract: The string vertices of closed string field theory are subsets of the moduli spaces of punctured Riemann surfaces that satisfy a geometric version of the Batalin-Vilkovisky master equation. We present a homological proof of existence of string vertices and their uniqueness up to canonical transformations. Using hyperbolic metrics on surfaces with geodesic boundaries we give an exact construction of string vertices as sets of surfaces with systole greater than or equal to $L$ with $L\leq 2\, \hbox{arcsinh}\, 1$. Intrinsic hyperbolic collars prevent the appearance of short geodesics upon sewing. The surfaces generated by Feynman diagrams are naturally endowed with Thurston metrics: hyperbolic on the vertices and flat on the propagators. For the classical theory the length $L$ is arbitrary and, as $L\to \infty$ hyperbolic vertices become the minimal-area vertices of closed string theory.


title: On the Boundary Layer Equations with Phase Transition in the Kinetic Theory of Gases

abstract: Consider the steady Boltzmann equation with slab symmetry for a monatomic, hard sphere gas in a half space. At the boundary of the half space, it is assumed that the gas is in contact with its condensed phase. The present paper discusses the existence and uniqueness of a uniformly decaying boundary layer type solution of the Boltzmann equation in this situation, in the vicinity of the Maxwellian equilibrium with zero bulk velocity, with the same temperature as that of the condensed phase, and whose pressure is the saturating vapor pressure at the temperature of the interface. This problem has been extensively studied first by Y. Sone, K. Aoki and their collaborators, by means of careful numerical simulations. See section 2 of [C. Bardos, F. Golse, Y. Sone: J. Stat. Phys. 124 (2006), 275-300] for a very detailed presentation of these works. More recently T.-P. Liu and S.-H. Yu [Arch. Rational Mech. Anal. 209 (2013), 869-997] have proposed an extensive mathematical strategy to handle the problems studied numerically by Y. Sone, K. Aoki and their group. The present paper offers an alternative, possibly simpler proof of one of the results discussed in [T.P. Liu, S.-H. Yu, loc. cit.]


title: Efficient Pourbaix diagrams of many-element compounds

abstract: Pourbaix diagrams have long been an essential tool for determining the phase stability of solids and their associated ionic species under electrochemical conditions. In recent years, Pourbaix diagrams have been used for applications ranging from corrosion-resistance alloy design to electrocatalysis, and the data from which they are generated has been enhanced by the availability of materials data in various online databases. However, generation of multi-element Pourbaix diagrams has a critical bottleneck which makes 3-element systems difficult to analyze quickly, and 4 and 5 element systems intractable. In this work, we present a method for constructing Pourbaix diagrams which uses a pre-processing step to circumvent the most egregious computational bottleneck and make many-element Pourbaix diagrams computationally efficient.


title: Equivalence groupoid of a class of general Burgers-Korteweg-de Vries equations with space-dependent coefficients

abstract: We describe the equivalence groupoid of the class of general Burgers - Korteweg - de Vries equations with space-dependent coefficients. This class is shown to reduce by a family of equivalence transformations to a subclass whose usual equivalence group is four-dimensional. Classified are admissible transformations of this subclass and singled out are its subclasses admitting maximal nontrivial conditional equivalence groups. All of them turn out to have dimension higher than four. In particular, a few new examples of nontrivial cases of normalization in the generalized sense of classes of differential equations appeared this way.


title: Relativistic quantum oscillators in the global monopole spacetime

abstract: We investigated the effects of the global monopole spacetime on the Dirac and Klein-Gordon relativistic quantum oscillators. In order to do this, we solve the Dirac and Klein-Gordon equations analytically and discuss the influence of this background which is characterized by the curvature of the spacetime on the energy profiles of these oscillators. In addition, we introduce a hard-wall potential and, for a particular case, determine the energy spectrum for relativistic quantum oscillators in this background.


title: Limit theorems for generalized density-dependent Markov chains and bursty stochastic gene regulatory networks

abstract: Stochastic gene regulatory networks with bursting dynamics can be modeled mesocopically as a generalized density-dependent Markov chain (GDDMC) or macroscopically as a piecewise-deterministic Markov process (PDMP). Here we prove a limit theorem showing that each family of GDDMCs will converge to a PDMP as the system size tends to infinity. Moreover, under a simple dissipative condition, we prove the existence and uniqueness of the stationary distribution and the exponential ergodicity for the PDMP limit via the coupling method. Further extensions and applications to single-cell stochastic gene expression kinetics and bursty stochastic gene regulatory networks are also discussed and the convergence of the stationary distribution of the GDDMC model to that of the PDMP model is also proved.


title: The arithmetic basilica: a quadratic PCF arboreal Galois group

abstract: The arboreal Galois group of a polynomial $f$ over a field $K$ encodes the action of Galois on the iterated preimages of a root point $x_0\in K$, analogous to the action of Galois on the $\ell$-power torsion of an abelian variety. We compute the arboreal Galois group of the postcritically finite polynomial $f(z) = z^2 - 1$ when the field $K$ and root point $x_0$ satisfy a simple condition. We call the resulting group the arithmetic basilica group because of its relation to the basilica group associated with the complex dynamics of $f$. For $K=\mathbb{Q}$, our condition holds for infinitely many choices of $x_0$.


title: Handling Syntactic Divergence in Low-resource Machine Translation

abstract: Despite impressive empirical successes of neural machine translation (NMT) on standard benchmarks, limited parallel data impedes the application of NMT models to many language pairs. Data augmentation methods such as back-translation make it possible to use monolingual data to help alleviate these issues, but back-translation itself fails in extreme low-resource scenarios, especially for syntactically divergent languages. In this paper, we propose a simple yet effective solution, whereby target-language sentences are re-ordered to match the order of the source and used as an additional source of training-time supervision. Experiments with simulated low-resource Japanese-to-English, and real low-resource Uyghur-to-English scenarios find significant improvements over other semi-supervised alternatives.


title: Maximum Amplification of Enstrophy in 3D Navier-Stokes Flows

abstract: This investigation concerns a systematic search for potentially singular behavior in 3D Navier-Stokes flows. Enstrophy serves as a convenient indicator of the regularity of solutions to the Navier Stokes system --- as long as this quantity remains finite, the solutions are guaranteed to be smooth and satisfy the equations in the classical (pointwise) sense. However, there are no estimates available with finite a priori bounds on the growth of enstrophy and hence the regularity problem for the 3D Navier-Stokes system remains open. In order to quantify the maximum possible growth of enstrophy, we consider a family of PDE optimization problems in which initial conditions with prescribed enstrophy $\mathcal{E}_0$ are sought such that the enstrophy in the resulting Navier-Stokes flow is maximized at some time $T$. Such problems are solved computationally using a large-scale adjoint-based gradient approach derived in the continuous setting. By solving these problems for a broad range of values of $\mathcal{E}_0$ and $T$, we demonstrate that the maximum growth of enstrophy is in fact finite and scales in proportion to $\mathcal{E}_0^{3/2}$ as $\mathcal{E}_0$ becomes large. Thus, in such worst-case scenario the enstrophy still remains bounded for all times and there is no evidence for formation of singularity in finite time. We also analyze properties of the Navier-Stokes flows leading to the extreme enstrophy values and show that this behavior is realized by a series of vortex reconnection events.


title: Single-cell stochastic gene expression kinetics with coupled positive-plus-negative feedback

abstract: Here we investigate single-cell stochastic gene expression kinetics in a minimal coupled gene circuit with positive-plus-negative feedback. A triphasic stochastic bifurcation upon the increasing ratio of the positive and negative feedback strengths is observed, which reveals a strong synergistic interaction between positive and negative feedback loops. We discover that coupled positive-plus-negative feedback amplifies gene expression mean but reduces gene expression noise over a wide range of feedback strengths when promoter switching is relatively slow, stabilizing gene expression around a relatively high level. In addition, we study two types of macroscopic limits of the discrete chemical master equation model: the Kurtz limit applies to proteins with large burst frequencies and the L\'{e}vy limit applies to proteins with large burst sizes. We derive the analytic steady-state distributions of the protein abundance in a coupled gene circuit for both the discrete model and its two macroscopic limits, generalizing the results obtained in [Chaos 26:043108, 2016]. We also obtain the analytic time-dependent protein distribution for the classical Friedman-Cai-Xie random bursting model proposed in [Phys. Rev. Lett. 97:168302, 2006]. Our analytic results are further applied to study the structure of gene expression noise in a coupled gene circuit and a complete decomposition of noise in terms of five different biophysical origins is provided.

categories: q-bio.MN math.PR physics.bio-ph q-bio.QM


title: Declarative Programming for Microcontrollers -- Datalog on Arduino

abstract: In this paper we describe an approach to programming microcontrollers based on the Arduino platform using Datalog as a clear and concise description language for system behaviors. The application areas of cheap and easily programmable microcontrollers, like robotics, home automation, and IoT devices hold mainstream appeal and are often used as motivation in natural science and technology teaching. The choice of programming languages for microcontrollers is severely limited, especially with regard to rule-based declarative languages. We use an approach that is based on the Dedalus language augmented with operations that allow for side-effects and takes the limited resources of a microcontroller into account. Our compiler and runtime environment allow to run Datalog programs on Arduino-based systems.


title: Improving Performance of an Analog Electronic Device Using Quantum Error Correction

abstract: The use of analog classical systems for computation is generally thought to be a difficult proposition due to the susceptibility of these devices to noise and the lack of a clear framework for achieving fault-tolerance. We present experimental results for the application of quantum error correction (QEC) techniques to a prototype analog computational device called a quantum emulation device. It is shown that for the gates tested (transversal $Z$, $X$ and $SH$) there is a marked improvement in the performance characteristics of the gate operations following error correction using the 5-Qubit Perfect code. In the case of the $Z$ gate, the median fidelity improved from 0.995 to 0.99998, a reduction in the gate error by over two orders of magnitude. Other transverse gates similarly show strong improvements.


title: An active smartphone authentication method based on daily cyclical activity

abstract: Smartphones have become an important tool for people's daily lives, which brings higher security requirements in high-risk application areas, for example, mobile payment. Although the combination of physical password, fingerprint and facial recognition have improved the security to a certain extent, there still exists a high risk of being decrepted. This paper attempts an algorithm which is more suitable for studying human partial periodic activity, namely Prophet algorithm. This algorithm has strong robustness for missing data and trend change, and can deal with outliers well. The experimental results on the UniMiB SHAR DATA show that the user simply needs to do 5 cycles of specified actions to realize the prediction of the next time series. The Error analysis of cross validation was applied to 4 different indicators, and the Mean Squared Error of the optimal result "Jumping" behavior was only 8.20%. With these appealing features, The main contribution of this paper is to propose a smart phone user identification system based on behavioral activity cycle, which can be replicated in other behavioral studies. Another outstanding feature of such a system is the capability of fitting models using small data set by exploiting behavioral characteristics derived from periodicity and thus reducing dependence on sensor scanning frequency, therefore the system balances among energy consumption, data quantity and fitting accuracy.


title: Characterizing magnetic field morphologies in three Serpens protostellar cores with ALMA

abstract: With the aim of characterizing the dynamical processes involved in the formation of young protostars, we present high angular resolution ALMA dust polarization observations of the Class 0 protostellar cores Serpens SMM1, Emb 8(N), and Emb 8. With spatial resolutions ranging from 150 to 40 au at 870 {\mu}m, we find unexpectedly high values of the polarization fraction along the outflow cavity walls in Serpens Emb8(N). We use 3 mm and 1 mm molecular tracers to investigate outflow and dense gas properties and their correlation with the polarization. These observations allow us to investigate the physical processes involved in the Radiative Alignment Torques (RATs) acting on dust grains along the outflow cavity walls, which experience irradiation from accretion processes and outflow shocks. The inner core of SMM1-a presents a polarization pattern with a poloidal magnetic field at the bases of the two lobes of the bipolar outflow. To the south of SMM1-a we see two polarized filaments, one of which seems to trace the redshifted outflow cavity wall. The other may be an accretion streamer of material infalling onto the central protostar. We propose that the polarized emission we see at millimeter wavelengths along the irradiated cavity walls can be reconciled with the expectations of RAT theory if the aligned grains present at '<' 500 au scales in Class 0 envelopes have grown larger than the 0.1 {\mu}m size of ISM dust grains. Our observations allow us to constrain the star-forming sources magnetic field morphologies within the central cores, along the outflow cavity walls, and in possible accretion streamers.