title: Exploiting RapidWright in the Automatic Generation of Application-Specific FPGA Overlays

abstract: Overlay architectures implemented on FPGA devices have been proposed as a means to increase FPGA adoption in general-purpose computing. They provide the benefits of software such as flexibility and programmability, thus making it easier to build dedicated compilers. However, existing overlays are generic, resource and power hungry with performance usually an order of magnitude lower than bare metal implementations. As a result, FPGA overlays have been confined to research and some niche applications. In this paper, we introduce Application-Specific FPGA Overlays (AS-Overlays), which can provide bare-metal performance to FPGA overlays, thus opening doors for broader adoption. Our approach is based on the automatic extraction of hardware kernels from data flow applications. Extracted kernels are then leveraged for application-specific generation of hardware accelerators. Reconfiguration of the overlay is done with RapidWright which allows to bypass the HDL design flow. Through prototyping, we demonstrated the viability and relevance of our approach. Experiments show a productivity improvement up to 20x compared to the state of the art FPGA overlays, while achieving over 1.33x higher Fmax than direct FPGA implementation and the possibility of lower resource and power consumption compared to bare metal.


title: Identifying Multiple Populations in M71 Using CN

abstract: We have observed the CN features at $\sim$3800 ~\AA~ and 4120~ \AA~ as well as the CH band at $\sim$4300~\AA~ for 145 evolved stars in the Galactic globular cluster M71 using the multi-object spectrograph, Hydra, on the WIYN-3.5m telescope. We use these measurements to create two $\delta$CN indices finding that both distributions are best fit by two populations, a CN-enhanced and CN-normal. We find that 42 $\pm$ 4\% of the RGB stars in our sample are CN-enhanced. The percentage of CN-enhanced is 40 $\pm$ 13\% for the AGB and 33 $\pm$ 9\% for the HB stars, which suggests there are no missing second generation stars at these stages of stellar evolution. The two generations also separate in magnitude and color on the HB, which allows us to find the difference in He abundance between the two populations by fitting appropriate ZAHBs. The broad range of distances from the cluster's center covered by our sample allows us to study the dependence of the ratio of the number of first to second population stars on the distance from the cluster's center, and we find that this ratio does not vary radially and that the two populations are spatially mixed. Finally, we compare our identification of multiple populations with the classification based on the Na-O anti-correlation and the HST UV photometry, and we find good agreement with both methods.


title: Viscous motion of spherical nanoparticles that scatter laser radiation in the Rayleigh regime

abstract: The mechanism of transverse radiation viscosity for nanospheres moving in laser field is analyzed. It is demonstrated that in the process of light scattering by these particles besides the force Fs accelerating them in the direction of radiation propagation and the gradient force Fg that is due to the spatial inhomogeneity of the light field, there are forces Fvisc that slow down the movement of particles in the transverse directions. These light viscosity forces are due to the Doppler shift in frequency of scattered radiation. The general expressions for these forces acting on particles that scatter radiation in the Rayleigh regime are derived and applied to estimate their effect on levitated nanospheres and also on slow electrons moving in the laser and magnetic fields. The possible experiments for observation the effects of light viscosity is discussed.


title: Frictional boundary layer effect on vortex condensation in rotating turbulent convection

abstract: We perform direct numerical simulations of rotating Rayleigh--B\'enard convection of fluids with low ($Pr=0.1$) and high ($Pr=5$) Prandtl numbers in a horizontally periodic layer with no-slip top and bottom boundaries. At both Prandtl numbers, we demonstrate the presence of an upscale transfer of kinetic energy that leads to the development of domain-filling vortical structures. Sufficiently strong buoyant forcing and rotation foster the quasi-two-dimensional turbulent state of the flow, despite the formation of plume-like vertical disturbances promoted by so-called Ekman pumping from the viscous boundary layer.


title: Phase transformation in boron under shock compression

abstract: Using first-principles molecular dynamics, we calculated the equation of state and shock Hugoniot of various boron phases. We find a large mismatch between Hugoniots based on existing knowledge of the equilibrium phase diagram and those measured by shock experiments, which could be reconciled if the $\alpha$-B$_{12}$/$\beta \rightarrow\gamma$-B$_{28}$ transition is significantly over-pressurized in boron under shock compression. Our results also indicate that there exists an anomaly and negative Clapeyron slope along the melting curve of boron at 100 GPa and 1500--3000 Kelvin. These results enable in-depth understanding of matter under shock compression, in particular the significance of compression-rate dependence of phase transitions and kinetic effects in experimental measurements.


title: Structured climate financing: valuation of CDOs on inhomogeneous asset pools

abstract: Recently, a number of structured funds have emerged as public-private partnerships with the intent of promoting investment in renewable energy in emerging markets. These funds seek to attract institutional investors by tranching the asset pool and issuing senior notes with a high credit quality. Financing of renewable energy (RE) projects is achieved via two channels: small RE projects are financed indirectly through local banks that draw loans from the fund's assets, whereas large RE projects are directly financed from the fund. In a bottom-up Gaussian copula framework, we examine the diversification properties and RE exposure of the senior tranche. To this end, we introduce the LH++ model, which combines a homogeneous infinitely granular loan portfolio with a finite number of large loans. Using expected tranche percentage notional (which takes a similar role as the default probability of a loan), tranche prices and tranche sensitivities in RE loans, we analyse the risk profile of the senior tranche. We show how the mix of indirect and direct RE investments in the asset pool affects the sensitivity of the senior tranche to RE investments and how to balance a desired sensitivity with a target credit quality and target tranche size.


title: Analysis of quasar magnitudes

abstract: Since their discovery the analysis of quasar magnitudes has generally required some form of evolution. Assuming that quasars do not have evolution this paper shows that they have a well-defined intrinsic magnitude distribution that is independent of cosmological models. However the average apparent magnitudes are essentially constant which means that the only cosmological information they contain is that the absolute flux density has a power law distribution. Thus quasar magnitudes, by themselves, are essentially useless for cosmological investigations.


title: Expression of the Holtsmark function in terms of hypergeometric $_2F_2$ and Airy $\mathrm{Bi}$ functions

abstract: The Holtsmark distribution has applications in plasma physics, for the electric-microfield distribution involved in spectral line shapes for instance, as well as in astrophysics for the distribution of gravitating bodies. It is one of the few examples of a stable distribution for which a closed-form expression of the probability density function is known. However, the latter is not expressible in terms of elementary functions. In the present work, we mention that the Holtsmark probability density function can be expressed in terms of hypergeometric function $_2F_2$ and of Airy function of the second kind $\mathrm{Bi}$ and its derivative. The new formula is simpler than the one proposed by Lee involving $_2F_3$ and $_3F_4$ hypergeometric functions.

categories: math-ph astro-ph.GA cs.GL math.MP physics.atom-ph


title: Graph Cepstrum: Spatial Feature Extracted from Partially Connected Microphones

abstract: In this paper, we propose an effective and robust method of spatial feature extraction for acoustic scene analysis utilizing partially synchronized and/or closely located distributed microphones. In the proposed method, a new cepstrum feature utilizing a graph-based basis transformation to extract spatial information from distributed microphones, while taking into account whether any pairs of microphones are synchronized and/or closely located, is introduced. Specifically, in the proposed graph-based cepstrum, the log-amplitude of a multichannel observation is converted to a feature vector utilizing the inverse graph Fourier transform, which is a method of basis transformation of a signal on a graph. Results of experiments using real environmental sounds show that the proposed graph-based cepstrum robustly extracts spatial information with consideration of the microphone connections. Moreover, the results indicate that the proposed method more robustly classifies acoustic scenes than conventional spatial features when the observed sounds have a large synchronization mismatch between partially synchronized microphone groups.


title: Relational Semigroups and Object-Free Categories

abstract: This note relates axioms for partial semigroups and monoids with those for small object-free categories, either with multiple monoidal units or with source and target maps. We discuss the adjunction of a zero element to both kinds of category and provide examples that separate the algebras considered.


title: Mixing convex-optimization bounds for maximum-entropy sampling

abstract: The maximum-entropy sampling problem is a fundamental and challenging combinatorial-optimization problem, with application in spatial statistics. It asks to find a maximum-determinant order-$s$ principal submatrix of an order-$n$ covariance matrix. Exact solution methods for this NP-hard problem are based on a branch-and-bound framework. Many of the known upper bounds for the optimal value are based on convex optimization. We present a methodology for "mixing" these bounds to achieve better bounds.


title: Learning Unitaries by Gradient Descent

abstract: We study the hardness of learning unitary transformations in $U(d)$ via gradient descent on time parameters of alternating operator sequences. We provide numerical evidence that, despite the non-convex nature of the loss landscape, gradient descent always converges to the target unitary when the sequence contains $d^2$ or more parameters. Rates of convergence indicate a "computational phase transition." With less than $d^2$ parameters, gradient descent converges to a sub-optimal solution, whereas with more than $d^2$ parameters, gradient descent converges exponentially to an optimal solution.


title: Recent Developments of Muon g-2 from Lattice QCD

abstract: One of the most promising quantities for the search of signatures of physics beyond the Standard Model is the anomalous magnetic moment $g-2$ of the muon, where a comparison of the experimental result with the Standard Model estimate yields a deviation of about $3.5~\sigma$. On the theory side, the largest uncertainty arises from the hadronic sector, namely the hadronic vacuum polarisation and the hadronic light-by-light scattering. I review recent progress in calculating the hadronic contributions to the muon $g-2$ from the lattice and discuss the prospects and challenges to match the precision of the upcoming experiments.


title: An efficient automated data analytics approach to large scale computational comparative linguistics

abstract: This research project aimed to overcome the challenge of analysing human language relationships, facilitate the grouping of languages and formation of genealogical relationship between them by developing automated comparison techniques. Techniques were based on the phonetic representation of certain key words and concept. Example word sets included numbers 1-10 (curated), large database of numbers 1-10 and sheep counting numbers 1-10 (other sources), colours (curated), basic words (curated). To enable comparison within the sets the measure of Edit distance was calculated based on Levenshtein distance metric. This metric between two strings is the minimum number of single-character edits, operations including: insertions, deletions or substitutions. To explore which words exhibit more or less variation, which words are more preserved and examine how languages could be grouped based on linguistic distances within sets, several data analytics techniques were involved. Those included density evaluation, hierarchical clustering, silhouette, mean, standard deviation and Bhattacharya coefficient calculations. These techniques lead to the development of a workflow which was later implemented by combining Unix shell scripts, a developed R package and SWI Prolog. This proved to be computationally efficient and permitted the fast exploration of large language sets and their analysis.


title: Espaces FC(g(F)) et endoscopie

abstract: Let F be a p-adic field and let G be a connected reductive group defined over F. We assume p is big. Denote g the Lie algebra of G. We normalize suitably a Fourier-transform on the space of smooth functions with compact support on g(F), which to f associates f^. In a preceding paper, we have defined the space FC(g(F)) of functions f such that the orbital integrals of f and of f^ are 0 for each element of g(F) which is not topologically nilpotent. These spaces are compatible with endoscopic transfer. We assume here that G is absolutely quasi-simple and simply connected. We define a decomposition of the space FC(g(F)) in a direct sum of subspaces such that the endoscopic transfer becomes (more or less) clear on each subspace. In particular, if G is quasi-split, we describe the subspace of "stable" elements in FC(g(F)).


title: Solvability for non-smooth Schr\"{o}dinger equations with singular potentials and square integrable data

abstract: We develop a holomorphic functional calculus for first-order operators $DB$ to solve boundary value problems for Schr\"{o}dinger equations $-\mathrm{div}\, A \nabla u + a V u = 0$ in the upper half-space $\mathbb{R}^{n+1}_+$ with $n\in\mathbb{N}$. This relies on quadratic estimates for $DB$, which are proved for coefficients $A,a,V$ that are independent of the transversal direction to the boundary, and comprised of a complex-elliptic pair $(A,a)$ that are bounded and measurable, and a singular potential $V$ in either $L^{n/2}(\mathbb{R}^n)$ or the reverse H\"{o}lder class $B^{q}(\mathbb{R}^n)$ with $q\geq\max\{\tfrac{n}{2},2\}$. In the latter case, square function bounds are also shown to be equivalent to non-tangential maximal function bounds. This allows us to prove that the (Dirichlet) Regularity and Neumann boundary value problems with $L^2(\mathbb{R}^n)$-data are well-posed if and only if certain boundary trace operators defined by the functional calculus are isomorphisms. We prove this property when the principal coefficient matrix $A$ has either a Hermitian or block structure. More generally, the set of all complex coefficients for which the boundary value problems are well-posed is shown to be open.


title: A spectral deferred correction method for incompressible flow with variable viscosity

abstract: This paper presents a semi-implicit spectral deferred correction (SDC) method for incompressible Navier-Stokes problems with variable viscosity and time-dependent boundary conditions. The proposed method integrates elements of velocity- and pressure-correction schemes, which yields a simpler pressure handling and a smaller splitting error than the SDPC method of Minion & Saye (J. Comput. Phys. 375: 797-822, 2018). Combined with the discontinuous Galerkin spectral-element method for spatial discretization it can in theory reach arbitrary order of accuracy in time and space. Numerical experiments in three space dimensions demonstrate up to order 12 in time and 17 in space for constant as well as varying, solution-dependent viscosity. Compared to SDPC the present method yields a substantial improvement of accuracy and robustness against order reduction caused by time-dependent boundary conditions.


title: Real-Time Deployment Aspects of C-Band and Millimeter-Wave 5G-NR Systems

abstract: Fifth-generation (5G) new radio (NR) deployments are being rolled out in both the C-band (3.3 - 5.0 GHz) and millimeter-wave (mmWave) band (24.5 - 29.5 GHz). For outdoor scenarios, the C-band is expected to provide wide area coverage and throughput uniformity, whereas the mmWave band is expected to provide ultra-high throughput to dedicated areas within the C-band coverage. Due to the differences in the frequency bands, both systems are expected to be designed with different transmit and receive parameters, naturally resulting in performance variations proportional to the chosen parameters. Unlike many previous works, this paper presents measurement evaluations in central Auckland, New Zealand, from a pre-commercial deployment of a single-user, single-cell 5G-NR system operating in both bands. The net throughput, coverage reliability, and channel rank are analyzed across the two bands with baseband and analog beamforming. Our results show that the C-band coverage is considerably better than mmWave, with a consistently higher channel rank. Furthermore, the spatial stationarity region (SSR) for the azimuth angles-of-departure (AODs) is characterized, and a model derived from the measured beam identities is presented. The SSR of azimuth AODs is seen to closely follow a gamma distribution.


title: Efficient computation of extreme excursion probabilities for dynamical systems

abstract: We develop a novel computational method for evaluating the extreme excursion probabilities arising for random initialization of nonlinear dynamical systems. The method uses a Markov chain Monte Carlo or a Laplace approximation approach to construct a biasing distribution that in turn is used in an importance sampling procedure to estimate the extreme excursion probabilities. The prior and likelihood of the biasing distribution are obtained by using Rice's formula from excursion probability theory. We use Gaussian mixture biasing distributions and approximate the non-Gaussian initial excitation by the method of moments to circumvent the linearity and Gaussianity assumptions needed by excursion probability theory. We demonstrate the effectiveness of this computational framework for nonlinear dynamical systems of up to 100 dimensions.


title: Verifying Tree Ensembles by Reasoning about Potential Instances

abstract: Imagine being able to ask questions to a black box model such as "Which adversarial examples exist?", "Does a specific attribute have a disproportionate effect on the model's prediction?" or "What kind of predictions could possibly be made for a partially described example?" This last question is particularly important if your partial description does not correspond to any observed example in your data, as it provides insight into how the model will extrapolate to unseen data. These capabilities would be extremely helpful as they would allow a user to better understand the model's behavior, particularly as it relates to issues such as robustness, fairness, and bias. In this paper, we propose such an approach for an ensemble of trees. Since, in general, this task is intractable we present a strategy that (1) can prune part of the input space given the question asked to simplify the problem; and (2) follows a divide and conquer approach that is incremental and can always return some answers and indicates which parts of the input domains are still uncertain. The usefulness of our approach is shown on a diverse set of use cases.


title: Zeta Functions and the (Linear) Logic of Markov Processes

abstract: The author introduced models of linear logic known as ''Interaction Graphs'' which generalise Girard's various geometry of interaction constructions. In this work, we establish how these models essentially rely on a deep connection between zeta functions and the execution of programs, expressed as a cocycle. This is first shown in the simple case of graphs, before begin lifted to dynamical systems. Focussing on probabilistic models, we then explain how the notion of graphings used in Interaction Graphs captures a natural class of sub-Markov processes. We then extend the realisability constructions and the notion of zeta function to provide a realisability model of second-order linear logic over the set of all (discrete-time) sub-Markov processes.


title: Rotational spectroscopy of singly $^{13}$C substituted isotopomers of propyne and determination of a semi-empirical equilibrium structure

abstract: Submillimeter spectra of three isotopomers of propyne containing one $^{13}$C atom were recorded in natural isotopic composition in the region of 426 GHz to 785 GHz. Additional measurements were carried out near 110 GHz. Combining these with earlier data resulted in greatly improved spectroscopic parameters which permit reliable extrapolations up to about 1.5 THz. Coupled cluster quantum-chemical calculations were carried out in order to assess the differences between equilibrium and ground state rotational parameters of these and many other isotopic species to evaluate semi-empirical equilibrium structural parameters. In addition, we estimated the main spectroscopic parameters of the isotopomers of propyne with two $^{13}$C atoms, which have not yet been studied in the laboratory, but which may be detectable in astronomical sources with a large amount of $^{13}$C compared to the dominant $^{12}$C.


title: Detecting inspiratory holds in bedside waveforms using a probabilistic modeling approach

abstract: In this paper we propose an algorithm for finding inspiratory holds in bedside waveforms using a probabilistic modeling approach. Inspiratory holds contain useful information about the properties of the patient. Finding inspiratory holds in bedside waveforms is a problem the author ran into when working on research with flow and pressure measurements from mechanical ventilators. In this paper, we propose a simple solution for detecting the inspiratory holds.


title: A Baker-Campbell-Hausdorff formula for the logarithm of permutations

abstract: The dynamics-from-permutations of classical Ising spins is studied for a chain of four spins. We obtain the Hamiltonian operator which is equivalent to the unitary permutation matrix that encodes assumed pairwise exchange interactions. It is shown how this can be summarized by an exact terminating Baker-Campbell-Hausdorff formula, which relates the Hamiltonian to a product of exponentiated two-spin exchange permutations. We briefly comment upon physical motivation and implications of this study.


title: Gravitating superconducting solitons in the (3+1)-dimensional Einstein gauged non-linear sigma-model

abstract: In this paper, we construct the first analytic examples of (3+1)-dimensional self-gravitating regular cosmic tube solutions which are superconducting, free of curvature singularities and with non-trivial topological charge in the Einstein-SU(2) non-linear sigma-model. These gravitating topological solitons at a large distance from the axis look like a (boosted) cosmic string with an angular defect given by the parameters of the theory, and near the axis, the parameters of the solutions can be chosen so that the metric is singularity free and without angular defect. The curvature is concentrated on a tube around the axis. These solutions are similar to the Cohen-Kaplan global string but regular everywhere, and the non-linear sigma-model regularizes the gravitating global string in a similar way as a non-Abelian field regularizes the Dirac monopole. Also, these solutions can be promoted to those of the fully coupled Einstein-Maxwell non-linear sigma-model in which the non-linear sigma-model is minimally coupled both to the U(1) gauge field and to General Relativity. The analysis shows that these solutions behave as superconductors as they carry a persistent current even when the U(1) field vanishes. Such persistent current cannot be continuously deformed to zero as it is tied to the topological charge of the solutions themselves. The peculiar features of the gravitational lensing of these gravitating solitons are shortly discussed.


title: A Private Quantum Bit String Commitment

abstract: We propose an entanglement-based quantum bit string commitment protocol whose composability is proven in the random oracle model. This protocol has the additional property of preserving the privacy of the committed message. Even though this property is not resilient against man-in-the-middle attacks, this threat can be circumvented by considering that the parties communicate through an authenticated channel. The protocol remains secure (but not private) if we realize the random oracles as physical unclonable functions in the so-called bad PUF model with access before the opening phase.


title: Comment on "Individual external dose monitoring of all citizens of Date City by passive dosimeter 5 to 51 months after the Fukushima NPP accident (series): 1."

abstract: We point out numerous inconsistencies and inappropriate statements in M. Miyazaki and H. Hayano, Journal of Radiological Protection 37, 1 (2016).


title: A Tool for Conducting User Studies on Mobile Devices

abstract: With the ever-growing interest in the area of mobile information retrieval and the ongoing fast development of mobile devices and, as a consequence, mobile apps, an active research area lies in studying users' behavior and search queries users submit on mobile devices. However, many researchers require to develop an app that collects useful information from users while they search on their phones or participate in a user study. In this paper, we aim to address this need by providing a comprehensive Android app, called Omicron, which can be used to collect mobile query logs and perform user studies on mobile devices. Omicron, at its current version, can collect users' mobile queries, relevant documents, sensor data as well as user activity and interaction data in various study settings. Furthermore, we designed Omicron in such a way that it is conveniently extendable to conduct more specific studies and collect other types of sensor data. Finally, we provide a tool to monitor the participants and their data both during and after the collection process.


title: Non-uniqueness for reflected rough differential equations

abstract: We give an example of a reflected diffferential equation which may have infinitely many solutions if the driving signal is rough enough (e.g. of infinite $p$-variation, for some $p>2$). For this equation, we identify a sharp condition on the modulus of continuity of the signal under which uniqueness holds. L\'evy's modulus for Brownian motion turns out to be a boundary case. We further show that in our example, non-uniqueness holds almost surely when the driving signal is a fractional Brownian motion with Hurst index $H < \frac{1}{2}$. The considered equation is driven by a two-dimensional signal with one component of bounded variation, so that rough path theory is not needed to make sense of the equation.


title: Plasma-environment effects on K lines of astrophysical interest III. IPs, K thresholds, radiative rates, and Auger widths in Fe ix - Fe xvi

abstract: Aims. In the context of black-hole accretion disks, we aim to compute the plasma-environment effects on the atomic parameters used to model the decay of K-vacancy states in moderately charged iron ions, namely Fe ix - Fe xvi. Methods. We used the fully relativistic multiconfiguration Dirac-Fock (MCDF) method approximating the plasma electron-nucleus and electron-electron screenings with a time-averaged Debye-Huckel potential. Results. We report modified ionization potentials, K-threshold energies, wavelengths, radiative emission rates, and Auger widths for plasmas characterized by electron temperatures and densities in the ranges $10^5$ - $10^7$ K and $10^{18}$ - $10^{22}$ cm$^{-3}$. Conclusions. This study confirms that the high-resolution X-ray spectrometers onboard the future XRISM and ATHENA space missions will be capable of detecting the lowering of the K edges of these ions due to the extreme plasma conditions occurring in accretion disks around compact objects.


title: Remarks on Sobolev-Morrey-Campanato spaces defined on $C^{0,\gamma}$ domains

abstract: We discuss a few old results concerning embedding theorems for Campanato and Sobolev-Morrey spaces adapting the formulations to the case of domains of class $C^{0,\gamma}$, and we present more recent results concerning the extension of functions from Sobolev-Morrey spaces defined on those domains. As a corollary of the extension theorem we obtain an embedding theorem for Sobolev-Morrey spaces on arbitrary $C^{0,\gamma}$ domains.


title: Assessment on Thermal Transport Properties of Group-III Nitrides: A Classical Molecular Dynamics Study with Transferable Tersoff Type Inter-atomic Potentials

abstract: In this study, by means of classical molecular dynamics simulations, we investigated the thermal transport properties of hexagonal single-layer, zinc-blend and wurtzite phases of BN, AlN, and GaN crystals, which are very promising for the application and design of high-quality electronic devices. With this in mind, we generated fully transferable Tersoff-type empirical inter-atomic potential parameter sets by utilizing an optimization procedure based on particle swarm optimization. The predicted thermal properties as well as the structural, mechanical and vibrational properties of all materials are in very good agreement with existing experimental and first-principles data. The impact of isotopes on thermal transport is also investigated and between $\sim$10 and 50\% reduction in phonon thermal transport with random isotope distribution is observed in BN and GaN crystals. Our investigation distinctly shows that the generated parameter sets are fully transferable and very useful in exploring the thermal properties of systems containing these nitrides.


title: Processive and Distributive Non-Equilibrium Networks Discriminate in Alternate Limits

abstract: We study biochemical reaction networks capable of product discrimination inspired by biological proofreading mechanisms. At equilibrium, product discrimination, the selective formation of a "correct" product with respect to an "incorrect product", is fundamentally limited by the free energy difference between the two products. However, biological systems often far exceed this limit, by using discriminatory networks that expend free energy to maintain non-equilibrium steady states. Non-equilibrium systems are notoriously difficult to analyze and no systematic methods exist for determining parameter regimes which maximize discrimination. Here we introduce a measure that can be computed directly from the biochemical rate constants and which provides a necessary condition for proofreading. Non-equilibrium systems can discriminate using binding energy (energetic) differences or using activation energy (kinetic) differences, but in both cases, proofreading is optimal when dissipation is maximized, so long as the necessary condition we introduce is satisfied. We show that these constraints can be in opposition and that this explains why the error rate is a non-monotonic function of the irreversible drive in the original proofreading scheme of Hopfield and Ninio. Finally, we introduce mixed networks, in which one product is favored energetically and the other kinetically. In such networks, sensitive product switching can be achieved simply by spending free energy. Biologically, this corresponds to the ability to select between products by driving a single reaction without network fine tuning. This may be used to explore alternate product spaces in challenging environments.


title: Interaction between low-level jets and wind farms in a stable atmospheric boundary layer

abstract: Low-level jets (LLJs) are the wind maxima in the lower regions of the atmosphere with a high wind energy potential. Here we use large-eddy simulations to study the effect of LLJ height on the flow dynamics in a wind farm with $10\times4$ turbines. We change the LLJ height and atmospheric thermal stratification by varying the surface cooling rate. We find that the first row power production is higher in the presence of a LLJ compared to a neutral reference case without LLJ. Besides, we show that the first row power production increases with decreasing LLJ height. Due to the higher turbulence intensity, the wind turbine wakes recover faster in a neutral boundary layer than in a stably stratified one. However, for strong thermal stratification with a low-height LLJ, the wake recovery can be faster than for the neutral reference case as energy can be entrained from the LLJ. Flow visualizations reveal that under stable stratification the growth of wind farm's internal boundary layer is restricted and the wind flows around the wind farm. Wind farms extract energy from LLJs through wake meandering and turbulent entrainment depending on the LLJ height. Both effects are advantageous for wake recovery, which is beneficial for the performance of downwind turbines. This finding is confirmed by an energy budget analysis, which reveals a significant increase in the kinetic energy flux in the presence of a LLJ. The jet strength reduces as it passes through consecutive turbine rows. For strong stratification, the combined effect of buoyancy destruction and turbulence dissipation is larger than the turbulent entrainment. Therefore, the power production of turbines in the back of the wind farm is relatively low for strong atmospheric stratifications.


title: On the Coverage Performance of Boolean-Poisson Cluster Models for Wireless Sensor Networks

abstract: In this paper, we consider wireless sensor networks (WSNs) with sensor nodes exhibiting clustering in their deployment. We model the coverage region of such WSNs by Boolean Poisson cluster models (BPCM) where sensors nodes' location is according to a Poisson cluster process (PCP) and each sensor has an independent sensing range around it. We consider two variants of PCP, in particular \matern and Thomas cluster process to form Boolean \matern and Thomas cluster models. We first derive the capacity functional of these models. Using the derived expressions, we compute the sensing probability of an event and compare it with sensing probability of a WSN modeled by a Boolean Poisson model where sensors are deployed according to a Poisson point process. We also derive the power required for each cluster to collect data from all of its sensors for the three considered WSNs. We show that a BPCM WSN has less power requirement in comparison to the Boolean Poisson WSN, but it suffers from lower coverage, leading to a trade-off between per-cluster power requirement and the sensing performance. A cluster process with desired clustering may provide better coverage while maintaining low power requirements.


title: Predicting Goal-directed Attention Control Using Inverse-Reinforcement Learning

abstract: Understanding how goal states control behavior is a question ripe for interrogation by new methods from machine learning. These methods require large and labeled datasets to train models. To annotate a large-scale image dataset with observed search fixations, we collected 16,184 fixations from people searching for either microwaves or clocks in a dataset of 4,366 images (MS-COCO). We then used this behaviorally-annotated dataset and the machine learning method of Inverse-Reinforcement Learning (IRL) to learn target-specific reward functions and policies for these two target goals. Finally, we used these learned policies to predict the fixations of 60 new behavioral searchers (clock = 30, microwave = 30) in a disjoint test dataset of kitchen scenes depicting both a microwave and a clock (thus controlling for differences in low-level image contrast). We found that the IRL model predicted behavioral search efficiency and fixation-density maps using multiple metrics. Moreover, reward maps from the IRL model revealed target-specific patterns that suggest, not just attention guidance by target features, but also guidance by scene context (e.g., fixations along walls in the search of clocks). Using machine learning and the psychologically-meaningful principle of reward, it is possible to learn the visual features used in goal-directed attention control.


title: A study of defect-based error estimates for the Krylov approximation of $\varphi$-functions

abstract: Prior recent work, devoted to the study of polynomial Krylov techniques for the approximation of the action of the matrix exponential ${\rm e}^{tA}v$, is extended to the case of associated $\varphi$-functions (which occur within the class of exponential integrators). In particular, a~posteriori error bounds and estimates, based on the notion of the defect (residual) of the Krylov approximation are considered. Computable error bounds and estimates are discussed and analyzed. This includes a new error bound which favorably compares to existing error bounds in specific cases. The accuracy of various error bounds is characterized in relation to corresponding Ritz values of $A$. Ritz values yield properties of the spectrum of $A$ (specific properties are known a~priori, e.g. for Hermitian or skew-Hermitian matrices) in relation to the actual starting vector $v$ and can be computed. This gives theoretical results together with criteria to quantify the achieved accuracy on the run. For other existing error estimates the reliability and performance is studied by similar techniques. Effects of finite precision (floating point arithmetic) are also taken into account.


title: Hydrogen bonding and percolation in propan-2-ol -- water liquid mixtures: X-ray diffraction experiments and computer simulations

abstract: Synchrotron X-ray diffraction measurements have been conducted on aqueous mixtures of propan-2-ol (a.k.a. isopropanol, or 2-propanol), for alcohol contents between 10 and 90 molar %, from room temperature down to 230 K. Molecular dynamics simulations, by using an all-atom parametrization of the propan-2-ol molecule and the well-known TIP4P/2005 water model, were able to provide semi-quantitative descriptions of the measured total structure factors. Various quantities related to hydrogen bonding, like hydrogen bond numbers, size distribution of cyclic entities and cluster size distributions, have been determined from the particle co-ordinates obtained from the simulations. The percolation threshold at room temperature could be estimated to be between isopropanol concentrations of 62 and 74 molar %, whereas at very low temperature, calculations yielded a value above 90 molar %.


title: Experimental characterization of the energetics of quantum logic gates

abstract: We characterize the energetic footprint of a two-qubit quantum gate from the perspective of non-equilibrium quantum thermodynamics. We experimentally reconstruct the statistics of energy and entropy fluctuations following the implementation of a controlled-unitary gate, linking them to the performance of the gate itself and the phenomenology of Landauer principle at the single-quantum level. Our work thus addresses the energetic cost of operating quantum circuits, a problem that is crucial for the grounding of the upcoming quantum technologies.


title: Quarkonium measurements in nucleus-nucleus collisions with ALICE

abstract: Heavy quarks are produced in the early stages of nucleus-nucleus collisions and can therefore provide important insight into the Quark-Gluon Plasma (QGP). Quarkonia are proposed as crucial probes to study the QGP. The extent of the medium modification for heavy-quark quarkonium production in heavy-ion collisions is measured in terms of a nuclear modification factor $R_{\rm AA}$, defined as the quarkonium yield in heavy-ion collisions divided by the relative quarkonium cross section in pp collisions scaled by the nuclear overlap function. A possible path-length dependent quarkonium dissociation, as well as a contribution of (re-)generation of quarkonia from heavy quarks in the medium, would lead to an azimuthal anisotropy of quarkonium production relative to the reaction plane. In this contribution, the recent ALICE measurements of quarkonium in Pb-Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV will be discussed for both mid- and forward rapidity. The dependence of $R_{\rm AA}$ on centrality and $p_{\rm T}$ for J/$\psi$, $\Upsilon$(1S), $\Upsilon$(2S), as well as the J/$\psi$ elliptic flow $v_{2}$ will be shown. The experimental data and the current theoretical model calculations will be also discussed.


title: Quantum Theory of Radiative Decay Rate and Frequency Shift of Surface Plasmon Modes

abstract: In this paper we study, in the time domain, the interaction between localized surface plasmons and photons in arbitrarily shaped metal nanoparticles, by using the Hopfield approach to quantize the plasmon modes, where the electron oscillations are represented by a harmonic matter field linearly coupled to the electromagnetic radiation. The plasmon - photon coupling gives rise to dressed plasmon modes. We have found that the radiation does not induce a significant coupling among the different quasi-electrostatic plasmon modes for particles of size up to the plasma wavelength, but causes a frequency shift and an exponential decay in time of the modes. By solving the equations governing the expectation values of the plasmon creation and annihilation operators, we obtain a new closed-form full-wave expression for the decay rate and for the frequency shift of the plasmon modes. It is non-perturbative and it only depends on the surface charge distribution of the quasi-electrostatic plasmon modes. We validate the expression against the Mie theory for a nano-sphere of radius comparable to the plasma wavelength. Eventually, we investigate the decay rate and the frequency shift of the plasmon modes in isolated and interacting nanoparticle of non-canonical shape, as their size increases up to the plasma wavelength.


title: A Heteroscedastic Uncertainty Model for Decoupling Sources of MRI Image Quality

abstract: Quality control (QC) of medical images is essential to ensure that downstream analyses such as segmentation can be performed successfully. Currently, QC is predominantly performed visually at significant time and operator cost. We aim to automate the process by formulating a probabilistic network that estimates uncertainty through a heteroscedastic noise model, hence providing a proxy measure of task-specific image quality that is learnt directly from the data. By augmenting the training data with different types of simulated k-space artefacts, we propose a novel cascading CNN architecture based on a student-teacher framework to decouple sources of uncertainty related to different k-space augmentations in an entirely self-supervised manner. This enables us to predict separate uncertainty quantities for the different types of data degradation. While the uncertainty measures reflect the presence and severity of image artefacts, the network also provides the segmentation predictions given the quality of the data. We show models trained with simulated artefacts provide informative measures of uncertainty on real-world images and we validate our uncertainty predictions on problematic images identified by human-raters.


title: Hausdorff dimension of frequency sets of univoque sequences

abstract: For integer $m\ge3$, we study the dynamical system $(\Lambda_m,\sigma_m)$ where $\Lambda_m$ is the set $\{w\in\{0,1\}^\mathbb{N}: w$ does not contain $0^m$ or $1^m\}$ and $\sigma_m$ is the shift map on $\{0,1\}^\mathbb{N}$ restricted to $\Lambda_m$, study the Bernoulli-type measures on $\Lambda_m$ and find out the unique equivalent $\sigma_m$-invariant ergodic probability measure. As an application, we obtain the Hausdorff dimension of the set of univoque sequences, the Hausdorff dimension of the set of sequences in which the lengths of consecutive $0$'s and consecutive $1$'s are bounded, and the Hausdorff dimension of their frequency subsets.


title: Dynamo in weakly collisional nonmagnetized plasmas impeded by Landau damping of magnetic fields

abstract: We perform fully kinetic simulations of flows known to produce dynamo in magnetohydrodynamics (MHD), considering scenarios with low Reynolds number and high magnetic Prandtl number, relevant for galaxy cluster scale fluctuation dynamos. We find that Landau damping on the electrons leads to a rapid decay of magnetic perturbations, impeding the dynamo. This collisionless damping process operates on spatial scales where electrons are nonmagnetized, reducing the range of scales where the magnetic field grows in high magnetic Prandtl number fluctuation dynamos. When electrons are not magnetized down to the resistive scale, the magnetic energy spectrum is expected to be limited by the scale corresponding to magnetic Landau damping or, if smaller, the electron gyroradius scale, instead of the resistive scale. In simulations we thus observe decaying magnetic fields where resistive MHD would predict a dynamo.


title: Two-Sample Testing for Event Impacts in Time Series

abstract: In many application domains, time series are monitored to detect extreme events like technical faults, natural disasters, or disease outbreaks. Unfortunately, it is often non-trivial to select both a time series that is informative about events and a powerful detection algorithm: detection may fail because the detection algorithm is not suitable, or because there is no shared information between the time series and the events of interest. In this work, we thus propose a non-parametric statistical test for shared information between a time series and a series of observed events. Our test allows identifying time series that carry information on event occurrences without committing to a specific event detection methodology. In a nutshell, we test for divergences of the value distributions of the time series at increasing lags after event occurrences with a multiple two-sample testing approach. In contrast to related tests, our approach is applicable for time series over arbitrary domains, including multivariate numeric, strings or graphs. We perform a large-scale simulation study to show that it outperforms or is on par with related tests on our task for univariate time series. We also demonstrate the real-world applicability of our approach on datasets from social media and smart home environments.


title: Tunable pseudo-magnetic fields for polaritons in strained metasurfaces

abstract: Artificial magnetic fields are revolutionizing our ability to manipulate neutral particles, by enabling the emulation of exotic phenomena once thought to be exclusive to charged particles. In particular, pseudo-magnetic fields generated by nonuniform strain in artificial lattices have attracted considerable interest because of their simple geometrical origin. However, to date, these strain-induced pseudo-magnetic fields have failed to emulate the tunability of real magnetic fields because they are dictated solely by the strain configuration. Here, we overcome this apparent limitation for polaritons supported by strained metasurfaces, which can be realized with classical dipole antennas or quantum dipole emitters. Without altering the strain configuration, we unveil how one can tune the pseudo-magnetic field by modifying the electromagnetic environment via an enclosing photonic cavity which modifies the nature of the interactions between the dipoles. Remarkably, due to the competition between short-range Coulomb interactions and long-range photon-mediated interactions, we find that the pseudo-magnetic field can be entirely switched off at a critical cavity height for any strain configuration. Consequently, by varying only the cavity height, we demonstrate a tunable Lorentz-like force that can be switched on/off and an unprecedented collapse and revival of polariton Landau levels. Unlocking this tunable pseudo-magnetism for the first time poses new intriguing questions beyond the paradigm of conventional tight-binding physics.


title: Efficient Differentially Private $F_0$ Linear Sketching

abstract: A powerful feature of linear sketches is that from sketches of two data vectors, one can compute the sketch of the difference between the vectors. This allows us to answer fine-grained questions about the difference between two data sets. In this work, we consider how to construct sketches for weighted $F_0$, i.e., the summed weights of the elements in the data set, that are small, differentially private, and computationally efficient. Let a weight vector $w\in(0,1]^u$ be given. For $x\in\{0,1\}^u$ we are interested in estimating $\Vert x\circ w\Vert_1$ where $\circ$ is the Hadamard product (entrywise product). Building on a technique of Kushilevitz et al.~(STOC 1998), we introduce a sketch (depending on $w$) that is linear over GF(2), mapping a vector $x\in \{0,1\}^u$ to $Hx\in\{0,1\}^\tau$ for a matrix $H$ sampled from a suitable distribution $\mathcal{H}$. Differential privacy is achieved by using randomized response, flipping each bit of $Hx$ with probability $p<1/2$. We show that for every choice of $0<\beta < 1$ and $\varepsilon=O(1)$ there exists $p<1/2$ and a distribution $\mathcal{H}$ of linear sketches of size $\tau = O(\log^2(u)\varepsilon^{-2}\beta^{-2})$ such that: 1) For random $H\sim\mathcal{H}$ and noise vector $\varphi$, given $Hx + \varphi$ we can compute an estimate of $\Vert x\circ w\Vert_1$ that is accurate within a factor $1\pm\beta$, plus additive error $O(\log(u)\varepsilon^{-2}\beta^{-2})$, with probability $1-1/u$, and 2) For every $H\sim\mathcal{H}$, $Hx + \varphi$ is $\varepsilon$-differentially private over the randomness in $\varphi$. The special case $w=(1,\dots,1)$ is unweighted $F_0$. Our results both improve the efficiency of existing methods for unweighted $F_0$ estimating and extend to a weighted generalization. We also give a distributed streaming implementation for estimating the size of the union between two input streams.


title: Global Hilbert expansion for the relativistic Vlasov-Maxwell-Boltzmann system

abstract: Consider the relativistic Vlasov-Maxwell-Boltzmann system describing the dynamics of an electron gas in the presence of a fixed ion background. Thanks to recent works \cite{Germain-Masmoudi-ASENS-2014, Guo-Ionescu-Pausader-JMP-2014} and \cite{Deng-Ionescu-Pausader-ARMA-2017}, we establish the global-in-time validity of its Hilbert expansion and derive the limiting relativistic Euler-Maxwell system as the mean free path goes to zero. Our method is based on the $L^2-L^{\infty}$ framework and the Glassey-Strauss Representation of the electromagnetic field, with auxiliary $H^1$ estimate and $W^{1,\infty}$ estimates to control the characteristic curves and corresponding $L^{\infty}$ norm.


title: Magnon valley Hall effect in CrI3-based vdW heterostructures

abstract: Magnonic excitations in the two-dimensional (2D) van der Waals (vdW) ferromagnet CrI3 are studied. We find that bulk magnons exhibit a non-trivial topological band structure without the need for Dzyaloshinskii-Moriya (DM) interaction. This is shown in vdW heterostructures, consisting of single-layer CrI3 on top of different 2D materials as MoTe2, HfS2 and WSe2. We find numerically that the proposed substrates modify substantially the out-of-plane magnetic anisotropy on each sublattice of the CrI3 subsystem. The induced staggered anisotropy, combined with a proper band inversion, leads to the opening of a topological gap of the magnon spectrum. Since the gap is opened non-symmetrically at the K+ and K- points of the Brillouin zone, an imbalance in the magnon population between these two valleys can be created under a driving force. This phenomenon is in close analogy to the so-called valley Hall effect (VHE), and thus termed as magnon valley Hall effect (MVHE). In linear response to a temperature gradient we quantify this effect by the evaluation of the temperature-dependence of the magnon thermal Hall effect. These findings open a different avenue by adding the valley degrees of freedom besides the spin, in the study of magnons.


title: A framework for large-scale mapping of human settlement extent from Sentinel-2 images via fully convolutional neural networks

abstract: Human settlement extent (HSE) information is a valuable indicator of world-wide urbanization as well as the resulting human pressure on the natural environment. Therefore, mapping HSE is critical for various environmental issues at local, regional, and even global scales. This paper presents a deep-learning-based framework to automatically map HSE from multi-spectral Sentinel-2 data using regionally available geo-products as training labels. A straightforward, simple, yet effective fully convolutional network-based architecture, Sen2HSE, is implemented as an example for semantic segmentation within the framework. The framework is validated against both manually labelled checking points distributed evenly over the test areas, and the OpenStreetMap building layer. The HSE mapping results were extensively compared to several baseline products in order to thoroughly evaluate the effectiveness of the proposed HSE mapping framework. The HSE mapping power is consistently demonstrated over 10 representative areas across the world. We also present one regional-scale and one country-wide HSE mapping example from our framework to show the potential for upscaling. The results of this study contribute to the generalization of the applicability of CNN-based approaches for large-scale urban mapping to cases where no up-to-date and accurate ground truth is available, as well as the subsequent monitor of global urbanization.


title: An Autonomous Intrusion Detection System Using an Ensemble of Advanced Learners

abstract: An intrusion detection system (IDS) is a vital security component of modern computer networks. With the increasing amount of sensitive services that use computer network-based infrastructures, IDSs need to be more intelligent and autonomous. Aside from autonomy, another important feature for an IDS is its ability to detect zero-day attacks. To address these issues, in this paper, we propose an IDS which reduces the amount of manual interaction and needed expert knowledge and is able to yield acceptable performance under zero-day attacks. Our approach is to use three learning techniques in parallel: gated recurrent unit (GRU), convolutional neural network as deep techniques and random forest as an ensemble technique. These systems are trained in parallel and the results are combined under two logics: majority vote and "OR" logic. We use the NSL-KDD dataset to verify the proficiency of our proposed system. Simulation results show that the system has the potential to operate with a very low technician interaction under the zero-day attacks. We achieved 87:28% accuracy on the NSL-KDD's "KDDTest+" dataset and 76:61% accuracy on the challenging "KDDTest-21" with lower training time and lower needed computational resources.


title: Singularity formation for the Serre-Green-Naghdi equations and applications to abcd-Boussinesq systems

abstract: In this work we prove that the solution of the Serre-Green-Naghdi equation cannot be globally defined when the interface reaches the impervious bottom tangentially. As a consequence, our result complements the paper \emph{Camassa, R., Falqui, G., Ortenzi, G., Pedroni, M., \& Thomson, C. Hydrodynamic models and confinement effects by horizontal boundaries. Journal of Nonlinear Science, 29(4), 1445-1498, 2019.} Furthermore, we also prove that the solution to the $abcd-$Boussinesq system can change sign in finite time. Finally, we provide with a proof of a scenario of finite time singularity for the $abcd-$Boussinesq system. These latter mathematical results are related to the numerics in \emph{Bona, \& Chen, Singular solutions of a Boussinesq system for water waves. J. Math. Study, 49(3), 205-220, 2016}.


title: TQFT, Symmetry Breaking, and Finite Gauge Theory in 3+1D

abstract: We derive a canonical form for 2-group gauge theory in 3+1D which shows they are either equivalent to Dijkgraaf-Witten theory or to the so-called "EF1" topological order of Lan-Wen. According to that classification, recently argued from a different point of view by Johnson-Freyd, this amounts to a very large class of all 3+1D TQFTs. We use this canonical form to compute all possible anomalies of 2-group gauge theory which may occur without spontaneous symmetry breaking, providing a converse of the recent symmetry-enforced-gaplessness constraints of C\'ordova-Ohmori and also uncovering some possible new examples. On the other hand, in cases where the anomaly is matched by a TQFT, we try to provide the simplest possible such TQFT. For example, with anomalies involving time reversal, $\mathbb{Z}_2$ gauge theory almost always works.


title: Centralized and distributed online learning for sparse time-varying optimization

abstract: The development of online algorithms to track time-varying systems has drawn a lot of attention in the last years, in particular in the framework of online convex optimization. Meanwhile, sparse time-varying optimization has emerged as a powerful tool to deal with widespread applications, ranging from dynamic compressed sensing to parsimonious system identification. In most of the literature on sparse time-varying problems, some prior information on the system's evolution is assumed to be available. In contrast, in this paper, we propose an online learning approach, which does not employ a given model and is suitable for adversarial frameworks. Specifically, we develop centralized and distributed algorithms, and we theoretically analyze them in terms of dynamic regret, in an online learning perspective. Further, we propose numerical experiments that illustrate their practical effectiveness.


title: Causal Structure Discovery from Distributions Arising from Mixtures of DAGs

abstract: We consider distributions arising from a mixture of causal models, where each model is represented by a directed acyclic graph (DAG). We provide a graphical representation of such mixture distributions and prove that this representation encodes the conditional independence relations of the mixture distribution. We then consider the problem of structure learning based on samples from such distributions. Since the mixing variable is latent, we consider causal structure discovery algorithms such as FCI that can deal with latent variables. We show that such algorithms recover a "union" of the component DAGs and can identify variables whose conditional distribution across the component DAGs vary. We demonstrate our results on synthetic and real data showing that the inferred graph identifies nodes that vary between the different mixture components. As an immediate application, we demonstrate how retrieval of this causal information can be used to cluster samples according to each mixture component.


title: A few topics in BFKL phenomenology at hadron colliders

abstract: This is a personal recollection of several results involving the phenomenological study of the multi-Regge limit of scattering amplitudes. None of them would have been possible without the encouragement and constant support from Lev Nikolaevich Lipatov.


title: Bayes-optimal Methods for Finding the Source of a Cascade

abstract: We study the problem of estimating the source of a network cascade. The cascade starts from a single vertex at time 0 and spreads over time, but only a noisy version of the propagation is observable. The goal is then to design a stopping time and estimator that will estimate the source well while ensuring the cost of the cascade to the system is not too large. We rigorously formulate a Bayesian approach to the problem. If vertices can be labelled by vectors in Euclidean space (which is natural in geo-spatial networks), the optimal estimator is the conditional mean estimator, and we derive an explicit form for the optimal stopping time under minimal assumptions on the cascade dynamics. We study the performance of the optimal stopping time on lattices, and show that a computationally efficient but suboptimal stopping time which compares the posterior variance to a threshold has near-optimal performance.


title: Extremal parameters and their duals for Fuchsian boundary maps

abstract: We describe arithmetic cross-sections for geodesic flow on compact surfaces of constant negative curvature using generalized Bowen-Series boundary maps and their natural extensions associated to cocompact torsion-free Fuchsian groups. If the boundary map parameters are extremal, that is, each is an endpoint of a geodesic that extends a side of the fundamental polygon, then the natural extension map has a domain with finite rectangular structure, and the associated arithmetic cross-section is parametrized by this set. This construction allows us to represent the geodesic flow as a special flow over a symbolic system of coding sequences. Moreover, each extremal parameter has a corresponding dual parameter choice such that the "past" of the arithmetic code of a geodesic is the "future" for the code using the dual parameter. This duality was observed for two classical parameter choices by Adler and Flatto; here we show constructively that every extremal parameter set has a dual.


title: Calderon-Zygmund type estimates for nonlocal PDE with H\"older continuous kernel

abstract: We study interior $L^p$-regularity theory, also known as Calderon-Zygmund theory, of the equation \[ \int_{\mathbb{R}^n} \int_{\mathbb{R}^n} \frac{K(x,y)\ (u(x)-u(y))\, (\varphi(x)-\varphi(y))}{|x-y|^{n+2s}}\, dx\, dy = \langle f, \varphi \rangle \quad \varphi \in C_c^\infty(\mathbb{R}^n). \] For $s \in (0,1)$, $t \in [s,2s]$, $p \in [2,\infty)$, $K$ an elliptic, symmetric, H\"older continuous kernel, if $f \in \left (H^{t,p'}_{00}(\Omega)\right )^\ast$, then the solution $u$ belongs to $H^{2s-t,p}_{loc}(\Omega)$ as long as $2s-t < 1$. The increase in differentiability is independent of the H\"older coefficient of $K$. For example, our result shows that if $f\in L^{p}_{loc}$ then $u\in H^{2s-\delta,p}_{loc}$ for any $\delta\in (0, s]$ as long as $2s-\delta < 1$. This is different than the classical analogue of divergence-form equations ${\rm div}(\bar{K} \nabla u) = f$ (i.e. $s=1$) where a $C^\gamma$-H\"older continuous coefficient $\bar{K}$ only allows for estimates of order $H^{1+\gamma}$. In fact, it is another appearance of the differential stability effect observed in many forms by many authors for this kind of nonlocal equations -- only that in our case we do not get a "small" differentiability improvement, but all the way up to $\min\{2s-t,1\}$. The proof argues by comparison with the (much simpler) equation \[ \int_{\mathbb{R}^n} K(z,z) (-\Delta)^{\frac{t}{2}} u(z) \, (-\Delta)^{\frac{2s-t}{2}} \varphi(z)\, dz = \langle g,\varphi\rangle \quad \varphi \in C_c^\infty(\mathbb{R}^n). \] and showing that as long as $K$ is H\"older continuous and $s,t, 2s-t \in (0,1)$ then the "commutator" \[ \int_{\mathbb{R}^n} K(z,z) (-\Delta)^{\frac{t}{2}} u(z) \, (-\Delta)^{\frac{2s-t}{2}} \varphi(z)\, dz - c\int_{\mathbb{R}^n} \int_{\mathbb{R}^n} \frac{K(x,y)\ (u(x)-u(y))\, (\varphi(x)-\varphi(y))}{|x-y|^{n+2s}}\, dx\, dy \] behaves like a lower order operator.


title: p-Value as the Strength of Evidence Measured by Confidence Distribution

abstract: The notion of p-value is a fundamental concept in statistical inference and has been widely used for reporting outcomes of hypothesis tests. However, p-value is often misinterpreted, misused or miscommunicated in practice. Part of the issue is that existing definitions of p-value are often derived from constructions under specific settings, and a general definition that directly reflects the evidence of the null hypothesis is not yet available. In this article, we first propose a general and rigorous definition of p-value that fulfills two performance-based characteristics. The performance-based definition subsumes all existing construction-based definitions of the p-value, and justifies their interpretations. The paper further presents a specific approach based on confidence distribution to formulate and calculate p-values. This specific way of computing p values has two main advantages. First, it is applicable for a wide range of hypothesis testing problems, including the standard one- and two-sided tests, tests with interval-type null, intersection-union tests, multivariate tests and so on. Second, it can naturally lead to a coherent interpretation of p-value as evidence in support of the null hypothesis, as well as a meaningful measure of degree of such support. In particular, it places a meaning of a large p-value, e.g. p-value of 0.8 has more support than 0.5. Numerical examples are used to illustrate the wide applicability and computational feasibility of our approach. We show that our proposal is effective and can be applied broadly, without further consideration of the form/size of the null space. As for existing testing methods, the solutions have not been available or cannot be easily obtained.


title: A Hybrid Quantum enabled RBM Advantage: Convolutional Autoencoders For Quantum Image Compression and Generative Learning

abstract: Understanding how the D-Wave quantum computer could be used for machine learning problems is of growing interest. Our work evaluates the feasibility of using the D-Wave as a sampler for machine learning. We describe a hybrid system that combines a classical deep neural network autoencoder with a quantum annealing Restricted Boltzmann Machine (RBM) using the D-Wave. We evaluate our hybrid autoencoder algorithm using two datasets, the MNIST dataset and MNIST Fashion dataset. We evaluate the quality of this method by using a downstream classification method where the training is based on quantum RBM-generated samples. Our method overcomes two key limitations in the current 2000-qubit D-Wave processor, namely the limited number of qubits available to accommodate typical problem sizes for fully connected quantum objective functions and samples that are binary pixel representations. As a consequence of these limitations we are able to show how we achieved nearly a 22-fold compression factor of grayscale 28 x 28 sized images to binary 6 x 6 sized images with a lossy recovery of the original 28 x 28 grayscale images. We further show how generating samples from the D-Wave after training the RBM, resulted in 28 x 28 images that were variations of the original input data distribution, as opposed to recreating the training samples. We formulated an MNIST classification problem using a deep convolutional neural network that used samples from a quantum RBM to train the MNIST classifier and compared the results with an MNIST classifier trained with the original MNIST training data set, as well as an MNIST classifier trained using classical RBM samples. Our hybrid autoencoder approach indicates advantage for RBM results relative to the use of a current RBM classical computer implementation for image-based machine learning and even more promising results for the next generation D-Wave quantum system.


title: Stability of Synchronized Steady State Solution of Diffusive Lotka-Volterra Predator-Prey Model

abstract: We study the reaction-diffusion Lotka-Volterra predator-prey model with Dirichlet boundary condition. In the case of equal diffusion rates and equal growth rates, the synchronized steady state solution is proved to be locally asymptotically stable.


title: The interplay between local and non-local master equations: exact and approximated dynamics

abstract: Master equations are a useful tool to describe the evolution of open quantum systems. In order to characterize the mathematical features and the physical origin of the dynamics, it is often useful to consider different kinds of master equations for the same system. Here, we derive an exact connection between the time-local and the integro-differential descriptions, focusing on the class of commutative dynamics. The use of the damping-basis formalism allows us to devise a general procedure to go from one master equation to the other and vice-versa, by working with functions of time and their Laplace transforms only. We further analyze the Lindbladian form of the time-local and the integro-differential master equations, where we account for the appearance of different sets of Lindbladian operators. In addition, we investigate a Redfield-like approximation, that transforms the exact integro-differential equation into a time-local one by means of a coarse graining in time. Besides relating the structure of the resulting master equation to those associated with the exact dynamics, we study the effects of the approximation on Markovianity. In particular, we show that, against expectation, the coarse graining in time can possibly introduce memory effects, leading to a violation of a divisibility property of the dynamics.


title: Rigid toric matrix Schubert varieties

abstract: For a given permutation $\pi \in S_N$, Fulton proves that the matrix Schubert variety $\overline{X_{\pi}} \cong Y_{\pi} \times \mathbb{C}^q$ can be defined via certain rank conditions encoded in the Rothe diagram of $\pi$. In the case where $Y_{\pi}:=\text{TV}(\sigma_{\pi})$ is toric (with respect to a $(\mathbb{C}^*)^{2N-1}$ action), we show that it can be described as an edge ideal of a bipartite graph $G^{\pi}$. We characterize the lower dimensional faces of the associated so-called edge cone $\sigma_{\pi}$ explicitly in terms of subgraphs of $G^{\pi}$ and present a combinatorial study for the first order deformations of $Y_{\pi}$. We prove that $Y_{\pi}$ is rigid if and only if the three-dimensional faces of $\sigma_{\pi}$ are all simplicial. Moreover, we reformulate this result in terms of Rothe diagram of $\pi$.


title: Non-reversibly updating a uniform [0,1] value for Metropolis accept/reject decisions

abstract: I show how it can be beneficial to express Metropolis accept/reject decisions in terms of comparison with a uniform [0,1] value, u, and to then update u non-reversibly, as part of the Markov chain state, rather than sampling it independently each iteration. This provides a small improvement for random walk Metropolis and Langevin updates in high dimensions. It produces a larger improvement when using Langevin updates with persistent momentum, giving performance comparable to that of Hamiltonian Monte Carlo (HMC) with long trajectories. This is of significance when some variables are updated by other methods, since if HMC is used, these updates can be done only between trajectories, whereas they can be done more often with Langevin updates. I demonstrate that for a problem with some continuous variables, updated by HMC or Langevin updates, and also discrete variables, updated by Gibbs sampling between updates of the continuous variables, Langevin with persistent momentum and non-reversible updates to u samples nearly a factor of two more efficiently than HMC. Benefits are also seen for a Bayesian neural network model in which hyperparameters are updated by Gibbs sampling.


title: AutoFCL: Automatically Tuning Fully Connected Layers for Handling Small Dataset

abstract: Deep Convolutional Neural Networks (CNN) have evolved as popular machine learning models for image classification during the past few years, due to their ability to learn the problem-specific features directly from the input images. The success of deep learning models solicits architecture engineering rather than hand-engineering the features. However, designing state-of-the-art CNN for a given task remains a non-trivial and challenging task, especially when training data size is less. To address this phenomena, transfer learning has been used as a popularly adopted technique. While transferring the learned knowledge from one task to another, fine-tuning with the target-dependent Fully Connected (FC) layers generally produces better results over the target task. In this paper, the proposed AutoFCL model attempts to learn the structure of FC layers of a CNN automatically using Bayesian optimization. To evaluate the performance of the proposed AutoFCL, we utilize five pre-trained CNN models such as VGG-16, ResNet, DenseNet, MobileNet, and NASNetMobile. The experiments are conducted on three benchmark datasets, namely CalTech-101, Oxford-102 Flowers, and UC Merced Land Use datasets. Fine-tuning the newly learned (target-dependent) FC layers leads to state-of-the-art performance, according to the experiments carried out in this research. The proposed AutoFCL method outperforms the existing methods over CalTech-101 and Oxford-102 Flowers datasets by achieving the accuracy of 94.38% and 98.89%, respectively. However, our method achieves comparable performance on the UC Merced Land Use dataset with 96.83% accuracy. The source codes of this research are available at https://github.com/shabbeersh/AutoFCL.


title: Existence and stability of steady state solutions of reaction-diffusion equations with nonlocal delay effect

abstract: A general reaction-diffusion equation with spatiotemporal delay and homogeneous Dirichlet boundary condition is considered. The existence and stability of positive steady state solutions are proved via studying an equivalent reaction-diffusion system without nonlocal and delay structure and applying local and global bifurcation theory. The global structure of the set of steady states is characterized according to type of nonlinearities and diffusion coefficient. Our general results are applied to diffusive logistic growth models and Nicholson's blowflies type models.


title: Orbital Angular Momentum Multiplexing in Highly Reverberant Environments

abstract: Previous studies on orbital angular momentum (OAM) communication mainly considered line-of-sight environments. In this letter, however, it is found that OAM communication with high-order modulation can be achieved in highly reverberant environments by combining the OAM multiplexing with a spatial equalizer. The OAM multiplexing exhibits comparable performance of conventional multiple-input multiple-output (MIMO) system.


title: MindReading: An Ultra-Low-Power Photonic Accelerator for EEG-based Human Intention Recognition

abstract: A scalp-recording electroencephalography (EEG)-based brain-computer interface (BCI) system can greatly improve the quality of life for people who suffer from motor disabilities. Deep neural networks consisting of multiple convolutional, LSTM and fully-connected layers are created to decode EEG signals to maximize the human intention recognition accuracy. However, prior FPGA, ASIC, ReRAM and photonic accelerators cannot maintain sufficient battery lifetime when processing real-time intention recognition. In this paper, we propose an ultra-low-power photonic accelerator, MindReading, for human intention recognition by only low bit-width addition and shift operations. Compared to prior neural network accelerators, to maintain the real-time processing throughput, MindReading reduces the power consumption by 62.7\% and improves the throughput per Watt by 168\%.


title: Applications of Lambert-Tsallis and Lambert-Kaniadakis Functions in Differential and Difference Equations with Deformed Exponential Decay

abstract: The analysis of a dynamical system modelled by differential (continuum case) or difference equation (discrete case) with deformed exponential decay, here we consider Tsallis and Kaniadakis exponentials, may require the use of the recently proposed deformed Lambert functions: the Lambert-Tsallis and Lambert-Kaniadakis functions. In this direction, the present work studies the logistic map with deformed exponential decay, using the Lambert-Tsallis and the Lambert-Kaniadakis functions to determine the stable behaviour and the dynamic of the disentropy in the weak chaotic regime. Furthermore, we investigate the motion of projectile when the vertical motion is governed by a non-linear differential equation with Tsallis exponential in the coefficient of the second order derivative. In this case, we calculated the range of the projectile using the Lambert-Tsallis function.


title: Documentation of Machine Learning Software

abstract: Machine Learning software documentation is different from most of the documentations that were studied in software engineering research. Often, the users of these documentations are not software experts. The increasing interest in using data science and in particular, machine learning in different fields attracted scientists and engineers with various levels of knowledge about programming and software engineering. Our ultimate goal is automated generation and adaptation of machine learning software documents for users with different levels of expertise. We are interested in understanding the nature and triggers of the problems and the impact of the users' levels of expertise in the process of documentation evolution. We will investigate the Stack Overflow Q/As and classify the documentation related Q/As within the machine learning domain to understand the types and triggers of the problems as well as the potential change requests to the documentation. We intend to use the results for building on top of the state of the art techniques for automatic documentation generation and extending on the adoption, summarization, and explanation of software functionalities.


title: Deep-AIR: A Hybrid CNN-LSTM Framework forFine-Grained Air Pollution Forecast

abstract: Poor air quality has become an increasingly critical challenge for many metropolitan cities, which carries many catastrophicphysical and mental consequences on human health and quality of life. However, accurately monitoring and forecasting air qualityremains a highly challenging endeavour. Limited by geographically sparse data, traditional statistical models and newly emergingdata-driven methods of air quality forecasting mainly focused on the temporal correlation between the historical temporal datasets of airpollutants. However, in reality, both distribution and dispersion of air pollutants are highly location-dependant. In this paper, we proposea novel hybrid deep learning model that combines Convolutional Neural Networks (CNN) and Long Short Term Memory (LSTM)together to forecast air quality at high-resolution. Our model can utilize the spatial correlation characteristic of our air pollutant datasetsto achieve higher forecasting accuracy than existing deep learning models of air pollution forecast.


title: Artificial Intelligence Aided Next-Generation Networks Relying on UAVs

abstract: Artificial intelligence (AI) assisted unmanned aerial vehicle (UAV) aided next-generation networking is proposed for dynamic environments. In the AI-enabled UAV-aided wireless networks (UAWN), multiple UAVs are employed as aerial base stations, which are capable of rapidly adapting to the dynamic environment by collecting information about the users' position and tele-traffic demands, learning from the environment and acting upon the feedback received from the users. Moreover, AI enables the interaction amongst a swarm of UAVs for cooperative optimization of the system. As a benefit of the AI framework, several challenges of conventional UAWN may be circumvented, leading to enhanced network performance, improved reliability and agile adaptivity. As a further benefit, dynamic trajectory design and resource allocation are demonstrated. Finally, potential research challenges and opportunities are discussed.


title: Essentially Optimal Sparse Polynomial Multiplication

abstract: We present a probabilistic algorithm to compute the product of two univariate sparse polynomials over a field with a number of bit operations that is quasi-linear in the size of the input and the output. Our algorithm works for any field of characteristic zero or larger than the degree. We mainly rely on sparse interpolation and on a new algorithm for verifying a sparse product that has also a quasi-linear time complexity. Using Kronecker substitution techniques we extend our result to the multivariate case.


title: Effect of spatial average on the spatiotemporal pattern formation of reaction-diffusion systems

abstract: Some quantities in the reaction-diffusion models from cellular biology or ecology depend on the spatial average of density functions instead of local density functions. We show that such nonlocal spatial average can induce instability of constant steady state, which is different from classical Turing instability. For a general scalar equation with spatial average, the occurrence of the steady state bifurcation is rigorously proved, and the formula to determine the bifurcation direction and the stability of the bifurcating steady state is given. For the two-species model, spatially non-homogeneous time-periodic orbits could arise due to spatially non-homogeneous Hopf bifurcation from the constant equilibrium. Examples from a nonlocal cooperative Lotka-Volterra model and a nonlocal Rosenzweig-MacArthur predator-prey model are used to demonstrate the bifurcation of spatially non-homogeneous patterns.


title: Middle-mile Network Optimization in Rural Wireless Meshes

abstract: The status quo of limited broadband connectivity in rural areas motivates the need for fielding alternatives such as long-distance wireless mesh networks. A key aspect of fielding wireless meshes cost-effectively is planning how to connect the last-mile networks to the core network service providers (i.e., the network between the edge access terminals and the landline / optical fiber terminals) with minimal infrastructure cost and throughput constraints. This so-called middle-mile network optimization, which includes topology construction, tower height assignment, antenna and orientation selection, as well as transmit power assignment, is known to be a computationally hard problem. In this paper, we provide the first polynomial time approximation solution for a generalized version of the middle-mile network optimization problem, wherein point-to-point (i.e., WiFi p2p) links are deployed to bridge last-mile networks. Our solution has a cost performance ratio of $O(\ln{|A|}+\frac{|B|}{|A|}+\frac{|A|+|B|}{\gamma})$, where A and B respectively denote the number of terminals and non-terminals and $\gamma$ is the ratio of $\frac{link\ capacity}{terminal\ demand}$. Furthermore, our solution extends to hybrid networks, i.e., point-to-multipoint (i.e., WiFi p2mp) or omnidirectional (i.e., TV White Space) can serve as hyperlinks in addition to point-to-point links, to further reduce the cost of wireless links. We provide a complementary heuristic for our middle-mile network optimization solution that adds hyperlinks if and only if they reduce the cost.


title: Modeling Events and Events of Events in Software Engineering

abstract: A model is a simplified representation of portion of reality that hides a system s nonessential characteristics. It provides a means for reducing complexity as well as visualization and communication and a basis for building it. Most models involve graphic languages during many of the software lifecycle phases. A new model, called thinging machine (TM), has recently been developed as an extension of the input-process-output framework. The paper focuses on events in a TM, offering a new perspective that captures a system s dynamic behaviors and a means of diagrammatically modeling events. The event notion is an important factor in giving semantics to specifications and providing a natural way to specify the interfaces and observable behavior of system components. Specifically, five generic TM event processes are analyzed: create, process, receive, release, and transfer. All events can be mapped (or reduced) to the events of these five event processes


title: Fast Monte Carlo Dropout and Error Correction for Radio Transmitter Classification

abstract: Monte Carlo dropout may effectively capture model uncertainty in deep learning, where a measure of uncertainty is obtained by using multiple instances of dropout at test time. However, Monte Carlo dropout is applied across the whole network and thus significantly increases the computational complexity, proportional to the number of instances. To reduce the computational complexity, at test time we enable dropout layers only near the output of the neural network and reuse the computation from prior layers while keeping, if any, other dropout layers disabled. Additionally, we leverage the side information about the ideal distributions for various input samples to do `error correction' on the predictions. We apply these techniques to the radio frequency (RF) transmitter classification problem and show that the proposed algorithm is able to provide better prediction uncertainty than the simple ensemble average algorithm and can be used to effectively identify transmitters that are not in the training data set while correctly classifying transmitters it has been trained on.


title: Cosmographic analysis of redshift drift

abstract: Redshift drift is the phenomenon whereby the observed redshift between an emitter and observer comoving with the Hubble flow in an expanding FLRW universe will slowly evolve -- on a timescale comparable to the Hubble time. There are nevertheless serious astrometric proposals for actually observing this effect. We shall however pursue a more abstract theoretical goal, and perform a general cosmographic analysis of this effect, eschewing (for now) dynamical considerations in favour of purely kinematic symmetry considerations and Taylor series expansions based on FLRW spacetimes. We shall develop various exact results and series expansions for the redshift drift (and its derivative) in terms of the present day Hubble, deceleration, jerk, snap, crackle, and pop parameters, as well as the present day redshift of the source. In particular, potential observation of this redshift drift effect is intimately related to the universe exhibiting a nonzero deceleration parameter.


title: Tenderbake -- A Solution to Dynamic Repeated Consensus for Blockchains

abstract: First-generation blockchains provide probabilistic finality: a block can be revoked, albeit the probability decreases as the block sinks deeper into the chain. Recent proposals revisited committee-based BFT consensus to provide deterministic finality: as soon as a block is validated, it is never revoked. A distinguishing characteristic of these second-generation blockchains over classical BFT protocols is that committees change over time as the participation and the blockchain state evolve. In this paper, we push forward in this direction by proposing a formalization of the Dynamic Repeated Consensus problem and by providing generic procedures to solve it in the context of blockchains. Our approach is modular in that one can plug in different synchronizers and single-shot consensus instances. To offer a complete solution, we provide a concrete instantiation, called Tenderbake, and present a blockchain synchronizer and a single-shot consensus algorithm, working in a Byzantine and partially synchronous system model with eventually synchronous clocks. In contrast to recent proposals, our methodology is driven by the need to bound the message buffers. This is essential in preventing spamming and run-time memory errors. Moreover, Tenderbake processes can synchronize with each other without exchanging messages, leveraging instead the information stored in the blockchain.


title: Measurement of the permanent electric dipole moment of the neutron

abstract: We present the result of an experiment to measure the electric dipole moment (EDM) of the neutron at the Paul Scherrer Institute using Ramsey's method of separated oscillating magnetic fields with ultracold neutrons (UCN). Our measurement stands in the long history of EDM experiments probing physics violating time reversal invariance. The salient features of this experiment were the use of a Hg-199 co-magnetometer and an array of optically pumped cesium vapor magnetometers to cancel and correct for magnetic field changes. The statistical analysis was performed on blinded datasets by two separate groups while the estimation of systematic effects profited from an unprecedented knowledge of the magnetic field. The measured value of the neutron EDM is $d_{\rm n} = (0.0\pm1.1_{\rm stat}\pm0.2_{\rm sys})\times10^{-26}e\,{\rm cm}$.


title: Euclidean correlation functions of the topological charge density

abstract: We present first results of our study on the Euclidean topological charge density correlation function. In order to get a well defined topological charge density and to improve the signal of the correlation function at large separations we make use of the gradient flow. We investigate the flow-time dependence on fine quenched lattices. The final goal of this study is to perform a continuum extrapolation for the pure SU(3) plasma and to extract the related transport coefficient, the sphaleron rate.


title: Electro-optical ion trap for experiments with atom-ion quantum hybrid systems

abstract: In the development of atomic, molecular and optical (AMO) physics, atom-ion hybrid systems are characterized by the presence of a new tool in the experimental AMO toolbox: atom-ion interactions. One of the main limitations in state-of-the-art atom-ion experiments is represented by the micromotion component of the ions' dynamics in a Paul trap, as the presence of micromotion in atom-ion collisions results in a heating mechanism that prevents atom-ion mixtures from undergoing a coherent evolution. Here we report the design and the simulation of a novel ion trapping setup especially conceived for the integration with an ultracold atoms experiment. The ion confinement is realized by using an electro-optical trap based on the combination of an optical and an electrostatic field, so that no micromotion component will be present in the ions' dynamics. The confining optical field is generated by a deep optical lattice created at the crossing of a bow-tie cavity, while a static electric quadrupole ensures the ions' confinement in the plane orthogonal to the optical lattice. The setup is also equipped with a Paul trap for cooling the ions produced by photoionization of a hot atomic beam, and the design of the two ion traps facilitates the swapping of the ions from the Paul trap to the electro-optical trap.


title: Long route to consensus: Two stage coarsening in a binary choice voting model

abstract: Formation of consensus, in binary yes/no type of voting, is a well defined process. However, even in presence of clear incentives, the dynamics involved can be incredibly complex. Specifically, formations of large groups of similarly opinionated individuals could create a condition of `support-bubbles' or spontaneous polarization that renders consensus virtually unattainable (e.g., the question of the UK exiting the EU). There have been earlier attempts in capturing the dynamics of consensus formation in societies through simple $Z_2$-symmetric models hoping to capture the essential dynamics of average behaviorof a large number of individuals in a statistical sense. However, in absence of external noise, they tend to reach a frozen state with fragmented and polarized states i.e., two or more groups of similarly opinionated groups with frozen dynamics. Here we show in a kinetic exchange opinion model (KEM) considered on $L \times L$ square lattices, that while such frozen states could be avoided, an exponentially slow approach to consensus is manifested. Specifically, the system could either reach consensus in a time that scales as $L^2$ or a long lived metastable state (termed a domain-wall state) for which formation of consensus takes a time scaling as $L^{3.6}$. The latter behavior is comparable to some voter-like models with intermediate states studied previously. The late-time anomaly in the time scale is reflected in the persistence probability of the model. Finally, the interval of zero-crossing of the average opinion i.e., the time interval over which the average opinion does not change sign is shown to follow a scale free distribution, which is compared with that seen in the opinion surveys regarding Brexit and associated issues in the last 40 years. The issue of minority spreading is also addressed by calculating the exit probability.


title: On the problem of maximal $L^q$-regularity for viscous Hamilton-Jacobi equations

abstract: For $q>2, \gamma > 1$, we prove that maximal regularity of $L^q$ type holds for periodic solutions to $-\Delta u + |Du|^\gamma = f$ in $\mathbb{R}^d$, under the (sharp) assumption $q > d \frac{\gamma-1}\gamma$.


title: Optimal Controller Synthesis and Dynamic Quantizer Switching for Linear-Quadratic-Gaussian Systems

abstract: In networked control systems, often the sensory signals are quantized before being transmitted to the controller. Consequently, performance is affected by the coarseness of this quantization process. Modern communication technologies allow users to obtain resolution-varying quantized measurements based on the prices paid. In this paper, we consider optimal controller synthesis of a Quantized-Feedback Linear-Quadratic-Gaussian (QF-LQG) system where the measurements are to be quantized before being transmitted to the controller. The system is presented with several choices of quantizers, along with the cost of operating each quantizer. The objective is to jointly select the quantizers and the controller that would maintain an optimal balance between control performance and quantization cost. Under certain assumptions, this problem can be decoupled into two optimization problems: one for optimal controller synthesis and the other for optimal quantizer selection. We show that, similarly to the classical LQG problem, the optimal controller synthesis subproblem is characterized by Riccati equations. On the other hand, the optimal quantizer selection policy is found by solving a certain Markov-Decision-Process (MDP).


title: Data-driven biological network alignment that uses topological, sequence, and functional information

abstract: Many proteins remain functionally unannotated. Sequence alignment (SA) uncovers missing annotations by transferring functional knowledge between species' sequence-conserved regions. Because SA is imperfect, network alignment (NA) complements SA by transferring functional knowledge between conserved biological network, rather than just sequence, regions of different species. Existing NA assumes that it is topological similarity (isomorphic-like matching) between network regions that corresponds to the regions' functional relatedness. However, we recently found that functionally unrelated proteins are almost as topologically similar as functionally related proteins. So, we redefined NA as a data-driven framework, TARA, which learns from network and protein functional data what kind of topological relatedness (rather than similarity) between proteins corresponds to the proteins' functional relatedness. TARA used topological information (within each network) but not sequence information (between proteins across networks). Yet, its alignments yielded higher protein functional prediction accuracy than alignments of existing NA methods, even those that used both topological and sequence information. Here, we propose TARA++ that is also data-driven, like TARA and unlike other existing methods, but that uses across-network sequence information on top of within-network topological information, unlike TARA. To deal with the within-and-across-network analysis, we adapt social network embedding to the problem of biological NA. TARA++ outperforms protein functional prediction accuracy of existing methods.


title: Unsatisfiability Proofs for Weight 16 Codewords in Lam's Problem

abstract: In the 1970s and 1980s, searches performed by L. Carter, C. Lam, L. Thiel, and S. Swiercz showed that projective planes of order ten with weight 16 codewords do not exist. These searches required highly specialized and optimized computer programs and required about 2,000 hours of computing time on mainframe and supermini computers. In 2011, these searches were verified by D. Roy using an optimized C program and 16,000 hours on a cluster of desktop machines. We performed a verification of these searches by reducing the problem to the Boolean satisfiability problem (SAT). Our verification uses the cube-and-conquer SAT solving paradigm, symmetry breaking techniques using the computer algebra system Maple, and a result of Carter that there are ten nonisomorphic cases to check. Our searches completed in about 30 hours on a desktop machine and produced nonexistence proofs of about 1 terabyte in the DRAT (deletion resolution asymmetric tautology) format.


title: Nonexistence Certificates for Ovals in a Projective Plane of Order Ten

abstract: In 1983, a computer search was performed for ovals in a projective plane of order ten. The search was exhaustive and negative, implying that such ovals do not exist. However, no nonexistence certificates were produced by this search, and to the best of our knowledge the search has never been independently verified. In this paper, we rerun the search for ovals in a projective plane of order ten and produce a collection of nonexistence certificates that, when taken together, imply that such ovals do not exist. Our search program uses the cube-and-conquer paradigm from the field of satisfiability (SAT) checking, coupled with a programmatic SAT solver and the nauty symbolic computation library for removing symmetries from the search.