title: Enhanced cosmological perturbations and the merger rate of PBH binaries

abstract: The rate of merger events observed by LIGO/Virgo can be used in order to probe the fraction $f$ of dark mater in the form of Primordial Black Holes. Here, we consider the merger rate of PBH binaries, accounting for the effect of cosmological perturbations on their initial eccentricity $e$. The torque on the binaries receives significant contributions from a wide range of scales, that goes from the size of the horizon at the time when the binary forms, down to the co-moving size of the binary. In scenarios where PBH are formed from adiabatic perturbations, it is natural to expect an enhancement of the power spectrum $P_\Phi$ at small scales, where it is poorly constrained observationally. The effect can then be quite significant. For instance, a nearly flat spectrum with amplitude $P_\Phi \gtrsim 10^{-7}$ on scales smaller than $\sim 10 Mpc^{-1}$ gives a contribution $\langle j^2 \rangle \sim 10^3 P_\Phi$, where $j= (1-e^2)^{1/2}$ is the dimensionless angular momentum parameter of the binaries. This contribution can dominate over tidal torques from neighboring PBHs for any value of $f$. Current constraints allow for a power spectrum as large as $P_\Phi \sim 10^{-5}$ at the intermediate scales $10^3-10^5 Mpc^{-1}$, comparable to the co-moving size of the binaries at the time of formation. In particular, this can relax current bounds on the PBH abundance based on the observed LIGO/Virgo merger rate, allowing for a fraction $f\sim 10\%$ of dark matter in PBH of mass $\sim 30 M_\odot$. We investigate the differential merger rate $\Delta\Gamma(m_1,m_2)$, as a function of the masses of the binary components, and the corresponding ``universality" coefficient $\alpha = -(m_1+m_2)^2 \partial^2 \ln \Delta\Gamma/\partial m_1\partial m_2$, concluding that merger rates may provide valuable information on the spectrum of primordial cosmological perturbations at currently uncharted lengthscales.


title: Exponential collapse with variable time scale driven by a scalar field

abstract: We study the dynamic collapse driven by a scalar field, when a relativistic observer falls co-moving with the collapse and cross the horizon of a Schwarzschild black-hole (BH), at $t=t_0$. During the collapse the scale of time is considered as variable. Back-reaction effects and gravitational waves produced during the exponential collapse are studied. We demonstrate that back-reaction effects act as the source of gravitational waves emitted during the collapse, and wavelengths of gravitational waves (GW) are in the range: $\lambda \ll r_s\equiv {e^{-2h_0t_0}\over 2 h_0}$, that is, smaller than the Schwarzschild radius. We demonstrate that during all the collapse the global topology of the space-time remains hyperbolic when the observer cross the horizon.


title: Semantic Driven Fielded Entity Retrieval

abstract: A common approach for knowledge-base entity search is to consider an entity as a document with multiple fields. Models that focus on matching query terms in different fields are popular choices for searching such entity representations. An instance of such a model is FSDM (Fielded Sequential Dependence Model). We propose to integrate field-level semantic features into FSDM. We use FSDM to retrieve a pool of documents, and then to use semantic field-level features to re-rank those documents. We propose to represent queries as bags of terms as well as bags of entities, and eventually, use their dense vector representation to compute semantic features based on query document similarity. Our proposed re-ranking approach achieves significant improvement in entity retrieval on the DBpedia-Entity (v2) dataset over existing FSDM model. Specifically, for all queries we achieve 2.5% and 1.2% significant improvement in NDCG@10 and NDCG@100, respectively.


title: Multiple competition-based FDR control for peptide detection

abstract: Competition-based FDR control has been commonly used for over a decade in the computational mass spectrometry community (Elias and Gygi, 2007). Recently, the approach has gained significant popularity in other fields after Barber and Candes (2015) laid its theoretical foundation in a more general setting that included the feature selection problem. In both cases, the competition is based on a head-to-head comparison between an observed score and a corresponding decoy / knockoff. Keich and Noble (2017b) recently demonstrated some advantages of using multiple rather than a single decoy when addressing the problem of assigning peptide sequences to observed mass spectra. In this work, we consider a related problem -- detecting peptides based on a collection of mass spectra -- and we develop a new framework for competition-based FDR control using multiple null scores. Within this framework, we offer several methods, all of which are based on a novel procedure that rigorously controls the FDR in the finite sample setting. Using real data to study the peptide detection problem we show that, relative to existing single-decoy methods, our approach can increase the number of discovered peptides by up to 50% at small FDR thresholds.


title: Coexistence of 11.2Tb/s Carrier-Grade Classical Channels and a DV-QKD Channel over a 7-Core Multicore Fibre

abstract: We successfully demonstrate coexistence of record-high 11.2 Tb/s (56x200Gb/s) classical channels with a discrete-variable-QKD channel over a multicore fibre. Continuous secret key generation is confirmed together with classical channel performance below the SDFEC limit and a minimum quantum channel spacing of 17nm in the C-band.


title: Gravitational effects on geonium and free electron $\mathrm{g}_s$-factor measurements in a Penning trap

abstract: We present a theoretical analysis of an electron confined by a Penning trap, also known as geonium, that is affected by gravity. In particular, we investigate the gravitational influence on the electron dynamics and the electromagnetic field of the trap. We consider the special case of a homogeneous gravitational field, which is represented by Rindler spacetime. In this spacetime the Hamiltonian of an electron with anomalous magnetic moment is constructed. Based on this Hamiltonian and the exact solution to Maxwell equations for the field of a Penning trap in Rindler spacetime, we derived the transition energies of geonium up to the relativistic corrections of $1/\mathrm{c}^2$. These transition energies are used to obtain an extension of the well known $\mathrm{g}_s$-factor formula introduced by L. S. Brown and G. Gabrielse [Rev. Mod. Phys. 58, 233 1986].


title: Large-scale multiscale modeling of phase transformation in nanocrystalline materials: Atomistic and Phase-Field methods

abstract: In this research, atomistic molecular dynamics simulations are combined with mesoscopic phase-field computational methods in order to investigate phase-transformation in polycrystalline Aluminum microstructure. In fact, microstructural computational modeling of engineering materials could help to optimize their mechanical properties for industrial applications (e.g. directional solidification for turbine blades). As a result, a multiscale modeling approach is developed to find a relation between manufacturing variables (e.g. temperature) and microstructural properties of crystalline materials (e.g. grain size), which could be used to develop an advanced manufacturing process for sensitive applications. The results show that atomistic modeling of grain growth could be used as a first-principle approach in order to study phase transformation's kinetics, which could capture morphology of polycrystalline materials more accurately. On the other hand, phase-field mesoscopic approach needs less computational efforts, but still it relies on semi-empirical data to capture accurate phase transformation regimes, which makes this approach suitable for rapid examining of new manufacturing conditions as well as its effects on microstructural properties of polycrystalline materials.


title: Giant and Linear Magnetoresistance in Liquid Metals at Ambient Temperature

abstract: Disorder-induced magnetoresistance has been reported in a range of solid metals and semiconductors, however, the underlying physical mechanism is still under debate because it is difficult to experimentally control. Liquid metals, due to lack of long-range order, offers an ideal model system where many forms of disorder can be deactivated by freezing the liquid. Here we report non-saturating magnetoresistance discovered in the liquid state of three metals: Ga, Ga-In-Sn and Bi-Pb-Sn-In alloys. The giant magnetoresistance appears above the respective melting points and has a maximum of 2500% at 14 Tesla. The reduced diamagnetism in the liquid state implies that a short-mean free path of the electron, induced by the spatial distribution of the liquid structure, is a key factor. A potential technological merit of this liquidtronic magnetoresistance is that it naturally operates at higher temperatures.

categories: cond-mat.mtrl-sci cond-mat.mes-hall cond-mat.other cond-mat.str-el


title: Reproducibility in Machine Learning for Health

abstract: Machine learning algorithms designed to characterize, monitor, and intervene on human health (ML4H) are expected to perform safely and reliably when operating at scale, potentially outside strict human supervision. This requirement warrants a stricter attention to issues of reproducibility than other fields of machine learning. In this work, we conduct a systematic evaluation of over 100 recently published ML4H research papers along several dimensions related to reproducibility. We find that the field of ML4H compares poorly to more established machine learning fields, particularly concerning data and code accessibility. Finally, drawing from success in other fields of science, we propose recommendations to data providers, academic publishers, and the ML4H research community in order to promote reproducible research moving forward.


title: The carry propagation of the successor function

abstract: Given any numeration system, we call carry propagation at a number $N$ the number of digits that are changed when going from the representation of $N$ to the one of $N+1$, and amortized carry propagation the limit of the mean of the carry propagations at the first $N$ integers, when $N$ tends to infinity, if this limit exists. In the case of the usual base $p$ numeration system, it can be shown that the limit indeed exists and is equal to $p/(p-1)$. We recover a similar value for those numeration systems we consider and for which the limit exists. We address the problem of the existence of the amortized carry propagation in non-standard numeration systems of various kinds: abstract numeration systems, rational base numeration systems, greedy numeration systems and beta-numeration. We tackle the problem by three different types of techniques: combinatorial, algebraic, and ergodic. For each kind of numeration systems that we consider, the relevant method allows for establishing sufficient conditions for the existence of the carry propagation and examples show that these conditions are close to being necessary conditions.


title: Correct estimate of the probability of false detection of the matched filter in weak-signal detection problems. III (Peak distribution method versus the Gumbel distribution method)

abstract: The matched filter (MF) represents one of the main tools to detect signals from known sources embedded in the noise. In the Gaussian case the noise is assumed to be the realization of a Gaussian random field (GRF). The most important property of the MF, the maximization of the probability of detection subject to a constant probability of false detection or false alarm (PFA), makes it one of the most popular techniques. However, the MF technique relies upon the a priori knowledge of the number and the position of the searched signals in the GRF which usually are not available. A typical way out is to assume that the position of a signal coincides with one of the peaks in the matched filtered data. A detection is claimed when the probability that a given peak is due only to the noise (i.e. the PFA) is smaller than a prefixed threshold. In this case the probability density function (PDF) of the amplitudes has to be used for the computation of the PFA, which is different from the Gaussian. Moreover, the probability that a detection is false depends on the number of peaks present in the filtered GRF, the greater the number of peaks in a GRF, the higher the probability of peaks due to the noise that exceed the detection threshold. If not taken into account, the PFA can be severely underestimated. Many solutions proposed to this problem are non-parametric hence not able to exploit all the available information. This limitation has been overcome by means of two efficient parametric approaches, one based on the PDF of the peak amplitudes of a smooth and isotropic GRF whereas the other uses the Gumbel distribution (the asymptotic PDF of the corresponding extreme). Simulations and ALMA maps show that, although the two methods produce almost identical results, the first is more flexible and allows us to check the reliability of the detection procedure.


title: Understanding Flaky Tests: The Developer's Perspective

abstract: Flaky tests are software tests that exhibit a seemingly random outcome (pass or fail) when run against the same, identical code. Previous work has examined fixes to flaky tests and has proposed automated solutions to locate as well as fix flaky tests--we complement it by examining the perceptions of software developers about the nature, relevance, and challenges of this phenomenon. We asked 21 professional developers to classify 200 flaky tests they previously fixed, in terms of the nature of the flakiness, the origin of the flakiness, and the fixing effort. We complement this analysis with information about the fixing strategy. Subsequently, we conducted an online survey with 121 developers with a median industrial programming experience of five years. Our research shows that: The flakiness is due to several different causes, four of which have never been reported before, despite being the most costly to fix; flakiness is perceived as significant by the vast majority of developers, regardless of their team's size and project's domain, and it can have effects on resource allocation, scheduling, and the perceived reliability of the test suite; and the challenges developers report to face regard mostly the reproduction of the flaky behavior and the identification of the cause for the flakiness. Data and materials [https://doi.org/10.5281/zenodo.3265785].


title: Electron enrichment of zigzag edges of armchair-oriented graphene nano-ribbons increases their stability and induces pinning of Fermi level

abstract: Zigzag edges of neutral armchair-oriented Graphene Nano-Ribbons show states strongly localized at those edges. They behave as free radicals that can capture electrons during processing, increasing ribbon's stability. Thus, charging and its consequences should be investigated.Total energy calculations of finite ribbons using spin polarized Density Functional Theory (DFT) show that ribbon's charging is feasible. Energies for Pariser-Parr-Pople (PPP) model Hamiltonian are compatible with DFT allowing the study of larger systems. Results for neutral ribbons indicate: i) the fundamental gap of spin polarized (non polarized) solutions is larger (smaller) than experimental data, ii) the ground state is spin polarized, a characteristic still not observed experimentally. Total energy of GNRs decreases with the number of captured electrons reaching a minimum for a number that mainly depends on zigzag edges size. The following changes with respect to neutral GNRs are noted: i) the ground state is not spin polarized, ii) fundamental gap is in-between that of spin polarized and non polarized solutions of neutral ribbons, iii) while in neutral ribbons valence and conduction band onsets vs. the fundamental gap, linearly and symmetrically approach mid-gap with slope 0.5, charging induces Fermi level pinning, i.e., the slopes of the valence and conduction bands being about 0.1 and 0.9, in agreement with experiment.


title: How we do things with words: Analyzing text as social and cultural data

abstract: In this article we describe our experiences with computational text analysis. We hope to achieve three primary goals. First, we aim to shed light on thorny issues not always at the forefront of discussions about computational text analysis methods. Second, we hope to provide a set of best practices for working with thick social and cultural concepts. Our guidance is based on our own experiences and is therefore inherently imperfect. Still, given our diversity of disciplinary backgrounds and research practices, we hope to capture a range of ideas and identify commonalities that will resonate for many. And this leads to our final goal: to help promote interdisciplinary collaborations. Interdisciplinary insights and partnerships are essential for realizing the full potential of any computational text analysis that involves social and cultural concepts, and the more we are able to bridge these divides, the more fruitful we believe our work will be.


title: Multi-Frequency Atom-Photon Interactions

abstract: We present a formalism that enables the analytic calculation of the interaction of a spin-half particle with a polychromatic electromagnetic field. This powerful new approach provides a clear physical picture even for cases with highly degenerate energy levels, which are complicated to interpret in the standard dressed-atom picture. Typically semi-classical methods are used for such problems (leading to equations that are solved by Floquet theory). Our formalism is derived from quantum electrodynamics and thus is more widely applicable. In particular it makes accessible the intermediate regime between quantum and semi-classical dynamics. We give examples of the application to multi-frequency multi-photon processes in strong fields by deriving the Hamiltonians of such systems, and also to the dynamics of weak fields at long times for which semi-classical methods are insufficient.


title: Augmenting Self-attention with Persistent Memory

abstract: Transformer networks have lead to important progress in language modeling and machine translation. These models include two consecutive modules, a feed-forward layer and a self-attention layer. The latter allows the network to capture long term dependencies and are often regarded as the key ingredient in the success of Transformers. Building upon this intuition, we propose a new model that solely consists of attention layers. More precisely, we augment the self-attention layers with persistent memory vectors that play a similar role as the feed-forward layer. Thanks to these vectors, we can remove the feed-forward layer without degrading the performance of a transformer. Our evaluation shows the benefits brought by our model on standard character and word level language modeling benchmarks.


title: On Injectivity of Quantum Finite Automata

abstract: We consider notions of freeness and ambiguity for the acceptance probability of Moore-Crutchfield Measure Once Quantum Finite Automata (MO-QFA). We study the injectivity problem of determining if the acceptance probability function of a MO-QFA is injective over all input words, i.e., giving a distinct probability for each input word. We show that the injectivity problem is undecidable for 8 state MO-QFA, even when all unitary matrices and the projection matrix are rational and the initial state vector is real algebraic. We also show undecidability of this problem when the initial vector is rational, although with a huge increase in the number of states. We utilize properties of quaternions, free rotation groups, representations of tuples of rationals as linear sums of radicals and a reduction of the mixed modification of Post's correspondence problem, as well as a new result on rational polynomial packing functions which may be of independent interest.


title: Neutrino Masses from the Point of View of Economy and Simplicity

abstract: In the framework of such basic principles as local gauge invariance, unification of the weak and electromagnetic interactions and spontaneous symmetry breaking in the Standard Model the most economical and simplest possibilities are realized. We discuss the problem of neutrino masses from the point of view of economy and simplicity. It is unlikely that neutrino masses are of the same SM origin as masses of leptons and quarks. The Weinberg effective Lagrangian is the simplest and the most economical, beyond the Standard Model mechanism of the generation of small Majorana neutrino masses. The resolution of the sterile neutrino anomaly and observation of the neutrinoless double $\beta$-decay would be crucial tests of this mechanism.


title: Degeneracy Index and Poincar\'e-Hopf Theorem

abstract: A degenerate dynamical system is characterized by a state-dependent multiplier of the time derivative of the state in the time evolution equation. It can give rise to Hamiltonian systems whose symplectic structure possesses a non-constant rank throughout the phase space. Around points where the multiplier becomes singular, flow can experience abrupt and irreversible changes. We introduce a topological index for degenerate dynamical systems around these {\it degeneracy points} and show that it refines and extends the usual topological index in accordance with the Poincar\'e-Hopf Theorem.


title: Memory of Motion for Warm-starting Trajectory Optimization

abstract: Trajectory optimization for motion planning requires good initial guesses to obtain good performance. In our proposed approach, we build a memory of motion based on a database of robot paths to provide good initial guesses. The memory of motion relies on function approximators and dimensionality reduction techniques to learn the mapping between the tasks and the robot paths. Three function approximators are compared: $k$-Nearest Neighbor, Gaussian Process Regression, and Bayesian Gaussian Mixture Regression. In addition, we show that the memory can be used as a metric to choose between several possible goals, and using an ensemble method to combine different function approximators results in a significantly improved warm-starting performance. We demonstrate the proposed approach with motion planning examples on the dual-arm robot PR2 and the humanoid robot Atlas.


title: Generalizing from a few environments in safety-critical reinforcement learning

abstract: Before deploying autonomous agents in the real world, we need to be confident they will perform safely in novel situations. Ideally, we would expose agents to a very wide range of situations during training, allowing them to learn about every possible danger, but this is often impractical. This paper investigates safety and generalization from a limited number of training environments in deep reinforcement learning (RL). We find RL algorithms can fail dangerously on unseen test environments even when performing perfectly on training environments. Firstly, in a gridworld setting, we show that catastrophes can be significantly reduced with simple modifications, including ensemble model averaging and the use of a blocking classifier. In the more challenging CoinRun environment we find similar methods do not significantly reduce catastrophes. However, we do find that the uncertainty information from the ensemble is useful for predicting whether a catastrophe will occur within a few steps and hence whether human intervention should be requested.


title: Some decay properties of hidden-charm pentaquarks as baryon-meson molecules

abstract: It is pointed out that the previously suggested interpretation of the hidden-charm pentaquarks $P_c(4312)$, $P_c(4440)$ and $P_c(4457)$ as molecular bound states of $\Sigma_c \bar D$ and $\Sigma_c \bar D^*$ baryon-meson pairs gives rise to specific relations for their decays. In particular, the heavy quark spin symmetry predicts ratios of rates of decays of each of the molecules to $J/\psi p$ and to $\eta_c p$ as well as of the decays to $\Lambda_c \bar D$ and to $\Lambda_c \bar D^*$. Experimental studies of these relations would thus provide an indicative probe of the molecular structure.


title: Taking drift-diffusion analysis from the study of turbulent flows to the study of particulate matter smog and air pollutants dynamics

abstract: Drift-diffusion analysis has been introduced in physics as a method to study turbulent flows. In the current study, it is proposed to use the method to identify underlying dynamical models of particulate matter smog, ozone and nitrogen dioxide concentrations. Data from Chiangmai are considered, which is a major city in the northern part of Thailand that recently has witnessed a dramatic increase of hospitalization that are assumed to be related to extreme air pollution levels. Three variants of the drift-diffusion analysis method (kernel-density, binning, linear approximation) are considered. It is shown that all three variants explain the annual pollutant peaks during the first half of the year by assuming that the parameters of the physical-chemical evolution equations of the pollutants vary periodically throughout the year. Therefore, our analysis provides evidence that the underlying dynamical models of the three pollutants being considered are explicitly time-dependent.


title: Obj-GloVe: Scene-Based Contextual Object Embedding

abstract: Recently, with the prevalence of large-scale image dataset, the co-occurrence information among classes becomes rich, calling for a new way to exploit it to facilitate inference. In this paper, we propose Obj-GloVe, a generic scene-based contextual embedding for common visual objects, where we adopt the word embedding method GloVe to exploit the co-occurrence between entities. We train the embedding on pre-processed Open Images V4 dataset and provide extensive visualization and analysis by dimensionality reduction and projecting the vectors along a specific semantic axis, and showcasing the nearest neighbors of the most common objects. Furthermore, we reveal the potential applications of Obj-GloVe on object detection and text-to-image synthesis, then verify its effectiveness on these two applications respectively.


title: Analytic and directional wavelet packets in the space of periodic signals

abstract: The paper presents a versatile library of analytic and quasi-analytic complex-valued wavelet packets (WPs) which originate from discrete splines of arbitrary orders. The real parts of the quasi-analytic WPs are the regular spline-based orthonormal WPs designed in [2]. The imaginary parts are the so-called complementary orthonormal WPs, which, unlike the symmetric regular WPs, they are antisymmetric. Tensor products of 1D quasi-analytic WPs provide a diversity of 2D WPs oriented in multiple directions. For example, a set of the fourth-level WPs comprises 62 different directions. The designed computational scheme in the paper enables us to get fast and easy implementation of the WP transforms. The presented WPs proved to be efficient in signal/image processing applications such as restoration of images degraded by either additive noise or missing of up to 90% of their pixels.


title: Exploring effective charge in electromigration using machine learning

abstract: The effective charge of an element is a parameter characterizing the electromgration effect, which can determine the reliability of interconnection in electronic technologies. In this work, machine learning approaches were employed to model the effective charge (z*) as a linear function of physically meaningful elemental properties. Average 5-fold (leave-out-alloy-group) cross-validation yielded root-mean-square-error divided by whole data set standard deviation (RMSE/$\sigma$) values of 0.37 $\pm$ 0.01 (0.22 $\pm$ 0.18), respectively, and $R^2$ values of 0.86. Extrapolation to z* of totally new alloys showed limited but potentially useful predictive ability. The model was used in predicting z* for technologically relevant host-impurity pairs.


title: HOnnotate: A method for 3D Annotation of Hand and Object Poses

abstract: We propose a method for annotating images of a hand manipulating an object with the 3D poses of both the hand and the object, together with a dataset created using this method. Our motivation is the current lack of annotated real images for this problem, as estimating the 3D poses is challenging, mostly because of the mutual occlusions between the hand and the object. To tackle this challenge, we capture sequences with one or several RGB-D cameras and jointly optimize the 3D hand and object poses over all the frames simultaneously. This method allows us to automatically annotate each frame with accurate estimates of the poses, despite large mutual occlusions. With this method, we created HO-3D, the first markerless dataset of color images with 3D annotations for both the hand and object. This dataset is currently made of 77,558 frames, 68 sequences, 10 persons, and 10 objects. Using our dataset, we develop a single RGB image-based method to predict the hand pose when interacting with objects under severe occlusions and show it generalizes to objects not seen in the dataset.


title: Hybrid static potential flux tubes from SU(2) and SU(3) lattice gauge theory

abstract: We compute chromoelectric and chromomagnetic flux densities for hybrid static potentials in SU(2) and SU(3) lattice gauge theory. In addition to the ordinary static potential with quantum numbers $\Lambda_\eta^\epsilon = \Sigma_g^+$, we present numerical results for seven hybrid static potentials corresponding to $\Lambda_\eta^{(\epsilon)} = \Sigma_u^+, \Sigma_g^-, \Sigma_u^-, \Pi_g, \Pi_u, \Delta_g, \Delta_u$, where the flux densities of five of them are studied for the first time in this work. We observe hybrid static potential flux tubes, which are significantly different from that of the ordinary static potential. They are reminiscent of vibrating strings, with localized peaks in the flux densities that can be interpreted as valence gluons.


title: Event-by-Event Multiplicity Correlations in $^{252}$Cf(sf)

abstract: Excited nuclear fragments are emitted during nuclear fission. The de-excitation of these fission fragments takes place as sequential emission of neutrons followed by photons. A correlation between neutron and photon multiplicities accompanying fission is thus expected. Fission event generators based on established statistical nuclear physics models predict a negative event-by-event correlation in neutron-photon multiplicity. A survey of published experimental results of an event-by-event covariance between the neutron and photon multiplicities emitted following the spontaneous fission of $^{252}$Cf is presented. Analytic unfolding expressions are developed in this work to determine the bias introduced by background sources, particle misclassification, pulse pileup, and inelastic photon production. The published experimental data are re-analyzed using these unfolding techniques and are found to be in qualitative agreement with the predictions of model-based calculations. In particular, we have concluded that there exists a significant event-by-event neutron-photon emission competition following the spontaneous fission of $^{252}$Cf.


title: Themis: Fair and Efficient GPU Cluster Scheduling

abstract: Modern distributed machine learning (ML) training workloads benefit significantly from leveraging GPUs. However, significant contention ensues when multiple such workloads are run atop a shared cluster of GPUs. A key question is how to fairly apportion GPUs across workloads. We find that established cluster scheduling disciplines are a poor fit because of ML workloads' unique attributes: ML jobs have long-running tasks that need to be gang-scheduled, and their performance is sensitive to tasks' relative placement. We propose Themis, a new scheduling framework for ML training workloads. It's GPU allocation policy enforces that ML workloads complete in a finish-time fair manner, a new notion we introduce. To capture placement sensitivity and ensure efficiency, Themis uses a two-level scheduling architecture where ML workloads bid on available resources that are offered in an auction run by a central arbiter. Our auction design allocates GPUs to winning bids by trading off efficiency for fairness in the short term but ensuring finish-time fairness in the long term. Our evaluation on a production trace shows that Themis can improve fairness by more than 2.25X and is ~5% to 250% more cluster efficient in comparison to state-of-the-art schedulers.


title: A transition from parabolic to conical shape as a common effect in nearby AGN jets

abstract: Observational studies of collimation in jets in active galactic nuclei (AGN) are a key to understanding their formation and acceleration processes. We have performed an automated search for jet shape transitions in a sample of 367 AGN using VLBA data at 15 GHz and 1.4 GHz. This search has found ten out of 29 nearby jets at redshifts z<0.07 with a transition from a parabolic to conical shape, while the full analyzed sample is dominated by distant AGN with a typical z about 1. The ten AGN are UGC00773, NGC1052, 3C111, 3C120, TXS0815-094, Mrk180, PKS1514+00, NGC6251, 3C371, and BL Lac. We conclude that the geometry transition may be a common effect in AGN jets. It can be observed only when sufficient linear resolution is obtained. Supplementing these results with previously reported shape breaks in the nearby AGN 1H0323+342 and M87, we estimate that the break occurs at 10^5-10^6 gravitational radii from the nucleus. We suggest that the jet shape transition happens when the bulk plasma kinetic energy flux becomes equal to the Poynting energy flux, while the ambient medium pressure is assumed to be governed by Bondi accretion. In general, the break point may not coincide with the Bondi radius. The observational data supports our model predictions on the jet acceleration and properties of the break point.


title: Optimal lower bounds for Donaldson's J-functional

abstract: In this paper we provide an explicit formula for the optimal lower bound of Donaldson's J-functional, in the sense of finding explicitly the optimal constant in the definition of coercivity, which always exists and takes negative values in general. This constant is positive precisely if the J-equation admits a solution, and the explicit formula has a number of applications. First, this leads to new existence criteria for constant scalar curvature K\"ahler (cscK) metrics in terms of Tian's alpha invariant. Moreover, we use the above formula to discuss Calabi dream manifolds and an analogous notion for the J-equation, and show that for surfaces the optimal bound is an explicitly computable rational function which typically tends to minus infinity as the underlying class approaches the boundary of the K\"ahler cone, even when the underlying K\"ahler classes admit cscK metrics. As a final application we show that if the Lejmi-Sz\'ekelyhidi conjecture holds, then the optimal bound coincides with its algebraic counterpart, the set of J-semistable classes equals the closure of the set of uniformly J-stable classes in the K\"ahler cone, and there exists an optimal degeneration for uniform J-stability.


title: Cosmology and dark energy from joint gravitational wave-GRB observations

abstract: Gravitational-wave (GW) detectors can contribute to the measurement of cosmological parameters and to testing the dark-energy sector of alternatives to $\Lambda$CDM, by using standard sirens. In this paper we focus on binary neutron stars with a counterpart detected through a gamma-ray burst (GRB), both at a second-generation network made by advanced LIGO+advanced Virgo+LIGO India+Kagra, and at third-generation (3G) detectors, discussing in particular the cases of a single Einstein Telescope (ET), and of a network of ET plus two Cosmic Explorer (CE). We construct mock catalogs of standard sirens, using different scenarios for the local merger rate and for the detection of the electromagnetic counterpart. For 3G detectors we estimate the coincidences with a GRB detector with the characteristics of the proposed THESEUS mission. We discuss how these standard sirens with a GRB counterpart can improve the determination of cosmological parameters (and particularly of $H_0$) in $\Lambda$CDM, and we then study how to extract information on dark energy, considering both a non-trivial dark energy equation of state and modified GW propagation. We find that a 2G detector network can already reach, over several years of data taking, an interesting sensitivity to modified GW propagation, while a single ET detector would have a remarkable potential for discovery. We also find that, to fully exploit the potential of a ET+CE+CE network, it is necessary a much stronger program of search for electromagnetic counterparts (or else to resort to statistical methods for standard sirens), and furthermore gravitational lensing can become a limiting factor.


title: Neutrino Masses and Mixing from Double Covering of Finite Modular Groups

abstract: We extend the even weight modular forms of modular invariant approach to general integral weight modular forms. We find that the modular forms of integral weights and level $N$ can be arranged into irreducible representations of the homogeneous finite modular group $\Gamma'_N$ which is the double covering of $\Gamma_N$. The lowest weight 1 modular forms of level 3 are constructed in terms of Dedekind eta-function, and they transform as a doublet of $\Gamma'_3 \cong T'$. The modular forms of weights 2, 3, 4, 5 and 6 are presented. We build a model of lepton masses and mixing based on $T'$ modular symmetry.


title: Secure Computation in Decentralized Data Markets

abstract: Decentralized data markets gather data from many contributors to create a joint data cooperative governed by market stakeholders. The ability to perform secure computation on decentralized data markets would allow for useful insights to be gained while respecting the privacy of data contributors. In this paper, we design secure protocols for such computation by utilizing secure multi-party computation techniques including garbled circuit evaluation and homomorphic encryption. Our proposed solutions are efficient and capable of performing arbitrary computation, but we report performance on two specific applications in the healthcare domain to emphasize the applicability of our methods to sensitive datasets.


title: An innovative adaptive kriging approach for efficient binary classification of mechanical problems

abstract: Kriging is an efficient machine-learning tool, which allows to obtain an approximate response of an investigated phenomenon on the whole parametric space. Adaptive schemes provide a the ability to guide the experiment yielding new sample point positions to enrich the metamodel. Herein a novel adaptive scheme called Monte Carlo-intersite Voronoi (MiVor) is proposed to efficiently identify binary decision regions on the basis of a regression surrogate model. The performance of the innovative approach is tested for analytical functions as well as some mechanical problems and is furthermore compared to two regression-based adaptive schemes. For smooth problems, all three methods have comparable performances. For highly fluctuating response surface as encountered e.g. for dynamics or damage problems, the innovative MiVor algorithm performs very well and provides accurate binary classification with only a few observation points.


title: Long-time asymptotics for evolutionary crystal dislocation models

abstract: We consider a family of evolution equations that generalize the Peierls-Nabarro model for crystal dislocations. They can be seen as semilinear parabolic reaction-diffusion equations in which the diffusion is regulated by a fractional Laplace operator of order $2 s \in (0, 2)$ acting in one space dimension and the reaction is determined by a $1$-periodic multi-well potential. We construct solutions of these equations that represent the typical propagation of $N \ge 2$ equally oriented dislocations of size $1$. For large times, the dislocations occur around points that evolve according to a repulsive dynamical system. When $s \in (1/2, 1)$, these solutions are shown to be asymptotically stable with respect to odd perturbations.


title: Slow-light frequency combs and dissipative Kerr solitons in coupled-cavity waveguides

abstract: We study Kerr frequency combs and dissipative Kerr solitons in silicon photonic crystal coupled-cavity waveguides (CCW) with globally optimized dispersion at telecom wavelengths. The corresponding threshold for comb generation is found to explicitly depend on the main CCW figures of merit, namely, mode volume, normal mode quality factor and slow-light group index. Our analysis is carried out by solving the non-linear dynamics of the CCW Bloch modes in presence of Kerr non-linearity and two-photon absorption. Our results open the way to CCW comb generation via dispersion engineering and slow-light enhancement.


title: Symmetry Configuration Mapping for Representing Quantum Systems on Quantum Computers

abstract: Quantum computing has the potential to significantly speed up complex computational tasks, and arguably the most promising application area for near-term quantum computers is the simulation of quantum mechanics. To make the most of our limited quantum computing resources, we need new and more compact algorithms and mappings. Whereas previous work, including recent demonstrations, has focused primarily on general mappings for fermionic systems, we propose instead to construct customized mappings tailored to the considered quantum mechanical systems. Specifically, we take advantage of existing symmetry, which we build into the mappings a priori to obtain optimal compactness. To demonstrate this approach, we have performed quantum computing calculations of the fluorine molecule, in which we have mapped 16 active spin-orbitals to 4 qubits. This is a four-fold reduction in the qubit requirement, as compared to the standard general mappings. Moreover, our compact system-to-qubits mappings are robust against noise that breaks symmetry, thereby reducing non-statistical errors in the computations. Furthermore, many systems, including F2, are described by real Hamiltonians, allowing us to also reduce the number of single-qubit operations in the hardware-efficient ansatz for the quantum variational eigensolver by roughly a factor of three.


title: Some remarks on "Convergence of Picard's iteration using projection algorithm for noncyclic contractions" [Indag. Math. 30 (2019) 227--239]

abstract: In this note, at first we prove that the existence of best proximity points for cyclic relatively nonexpansive mappings is equivalent to the existence of best proximity pairs for noncyclic relatively nonexpansive mappings in the setting of strictly convex Banach spaces by using the projection operator. In this way, we conclude that a main result of the paper "Proximal normal structure and relatively nonexpansive mappings", Studia Math., (171~(2005) 283--293) immediately follows. We then discuss the convergence of best proximity pairs for noncyclic contractions by applying the convergence of iterative sequences for cyclic contractions and show that the convergence method of a recent paper published in Indag. Math., 30(1) (2019) 227--239 is obtained exactly from Picard's iteration sequence.


title: Efficient Isomorphism for $S_d$-graphs and $T$-graphs

abstract: An $H$-graph is one representable as the intersection graph of connected subgraphs of a suitable subdivision of a fixed graph $H$, introduced by Bir\'{o}, Hujter and Tuza (1992). An $H$-graph is proper if the representing subgraphs of $H$ can be chosen incomparable by the inclusion. In this paper, we focus on the isomorphism problem for $S_d$-graphs and $T$-graphs, where $S_d$ is the star with $d$ rays and $T$ is an arbitrary fixed tree. Answering an open problem of Chaplick, T\"{o}pfer, Voborn\'{\i}k and Zeman (2016), we provide an FPT-time algorithm for testing isomorphism and computing the automorphism group of $S_d$-graphs when parameterized by~$d$, which involves the classical group-computing machinery by Furst, Hopcroft, and Luks (1980). We also show that the isomorphism problem of $S_d$-graphs is at least as hard as the isomorphism problem of posets of bounded width, for which no efficient combinatorial-only algorithm is known to date. Then we extend our approach to an XP-time algorithm for isomorphism of $T$-graphs when parameterized by the size of $T$. Lastly, we contribute a simple FPT-time combinatorial algorithm for isomorphism testing in the special case of proper $S_d$- and $T$-graphs.


title: Constraining Dark Matter -- Dark Radiation interactions with CMB, BAO, and Lyman-$\alpha$

abstract: Several interesting Dark Matter (DM) models invoke a dark sector leading to two types of relic particles, possibly interacting with each other: non-relativistic DM, and relativistic Dark Radiation (DR). These models have interesting consequences for cosmological observables, and could in principle solve problems like the small-scale cold DM crisis, Hubble tension, and/or low $\sigma_8$ value. Their cosmological behaviour is captured by the ETHOS parametrisation, which includes a DR-DM scattering rate scaling like a power-law of the temperature, $T^n$. Scenarios with $n=0$, $2$, or $4$ can easily be realised in concrete dark sector set-ups. Here we update constraints on these three scenarios using recent CMB, BAO, and high-resolution Lyman-$\alpha$ data. We introduce a new Lyman-$\alpha$ likelihood that is applicable to a wide range of cosmological models with a suppression of the matter power spectrum on small scales. For $n=2$ and $4$, we find that Lyman-$\alpha$ data strengthen the CMB+BAO bounds on the DM-DR interaction rate by many orders of magnitude. However, models offering a possible solution to the missing satellite problem are still compatible with our new bounds. For $n=0$, high-resolution Lyman-$\alpha$ data bring no stronger constraints on the interaction rate than CMB+BAO data, except for extremely small values of the DR density. Using CMB+BAO data and a theory-motivated prior on the minimal density of DR, we find that the $n=0$ model can reduce the Hubble tension from $4.1\sigma$ to $2.7\sigma$, while simultaneously accommodating smaller values of the $\sigma_8$ and $S_8$ parameters hinted by cosmic shear data.


title: Seismic data denoising and deblending using deep learning

abstract: An important step of seismic data processing is removing noise, including interference due to simultaneous and blended sources, from the recorded data. Traditional methods are time-consuming to apply as they often require manual choosing of parameters to obtain good results. We use deep learning, with a U-net model incorporating a ResNet architecture pretrained on ImageNet and further trained on synthetic seismic data, to perform this task. The method is applied to common offset gathers, with adjacent offset gathers of the gather being denoised provided as additional input channels. Here we show that this approach leads to a method that removes noise from several datasets recorded in different parts of the world with moderate success. We find that providing three adjacent offset gathers on either side of the gather being denoised is most effective. As this method does not require parameters to be chosen, it is more automated than traditional methods.

categories: eess.IV cs.LG physics.comp-ph physics.geo-ph


title: Boundary properties of fractional objects: flexibility of linear equations and rigidity of minimal graphs

abstract: The main goal of this article is to understand the trace properties of nonlocal minimal graphs in~$\R^3$, i.e. nonlocal minimal surfaces with a graphical structure. We establish that at any boundary points at which the trace from inside happens to coincide with the exterior datum, also the tangent planes of the traces necessarily coincide with those of the exterior datum. This very rigid geometric constraint is in sharp contrast with the case of the solutions of the linear equations driven by the fractional Laplacian, since we also show that, in this case, the fractional normal derivative can be prescribed arbitrarily, up to a small error. We remark that, at a formal level, the linearization of the trace of a nonlocal minimal graph is given by the fractional normal derivative of a fractional Laplace problem, therefore the two problems are formally related. Nevertheless, the nonlinear equations of fractional mean curvature type present very specific properties which are strikingly different from those of other problems of fractional type which are apparently similar, but diverse in structure, and the nonlinear case given by the nonlocal minimal graphs turns out to be significantly more rigid than its linear counterpart.


title: Some preliminary observations on a defect Navier-Stokes system

abstract: Some implications of the simplest accounting of defects of compatibility in the velocity field on the structure of the classical Navier-Stokes equations are explored, leading to connections between classical elasticity, the elastic theory of defects, plasticity theory, and classical fluid mechanics.

categories: physics.flu-dyn math-ph math.MP physics.class-ph


title: Lepton flavour universality in charged-current $B$ decays

abstract: Tests of lepton flavour universality in charged-current $B$ decays offer an excellent opportunity to test the Standard Model, and show hints of new physics in analyses performed by the LHCb, Belle and BaBar experiments. These proceedings present the results from the LHCb collaboration on measurements of $\mathcal{R}(D^*)$ and $\mathcal{R}(J/\psi)$. It also presents the latest semileptonic tag measurement of $\mathcal{R}(D)$ and $\mathcal{R}(D^*)$ by the Belle collaboration. The latest HFLAV average shows a discrepancy of 3.1$\sigma$ between the Standard Model predictions and combined measurements of $\mathcal{R}(D)$ and $\mathcal{R}(D^*)$.


title: Fermion-fermion duality in 3+1 dimensions

abstract: Dualities play a central role in both quantum field theories and condensed matter systems. Recently, a web of dualities has been discovered in 2+1 dimensions. Here, we propose in particular a generalization of the Son's fermion-fermion duality to 3+1 dimensions. We show that the action of charged Dirac fermions coupled to an external electromagnetic field is dual to an action of neutral fermions minimally coupled to an emergent vector gauge field. This dual action contains also a further tensor (Kalb-Ramond) gauge field coupled to the emergent and electromagnetic vector potentials. We firstly demonstrate the duality in the massive case. We then show the duality in the case of massless fermions starting from a lattice model and employing the slave-rotor approach already used in the 2+1-dimensional duality [Burkov, Phys. Rev. B 99, 035124 (2019)]. We finally apply this result to 3D Dirac semimetals in the low-energy regime. Besides the implications in topological phases of matter, our results shed light on the possible existence of a novel web of dualities in 3+1-dimensional (non-supersymmetric) quantum field theories.


title: Discrete effect on the anti-bounce-back boundary condition of lattice Bhatnagar-Gross-Krook model for convection-diffusion equations

abstract: The discrete effect on the boundary condition has been a fundamental topic for the lattice Boltzmann method in simulating heat and mass transfer problems. In previous works based on the halfway anti-bounce-back (ABB) boundary condition for convection-diffusion equations (CDEs), it is reported that the discrete effect cannot be commonly removed in the Bhatnagar-Gross-Krook (BGK) model except for a special value of relaxation time. Targeting this point in the present paper, we still proceed within the framework of BGK model for two-dimensional CDEs, and analyze the discrete effect on a non-halfway ABB boundary condition which incorporates the effect of the distance ratio. By analyzing an unidirectional diffusion problem with a parabolic distribution, the theoretical derivations with three different discrete velocity models show that the numerical slip is a combined function of the relaxation time and the distance ratio. Different from previous works, we definitely find that the relaxation time can be freely adjusted by the distance ratio in a proper range to eliminate the numerical slip. Some numerical simulations are carried out to validate the theoretical derivations, and the numerical results for the cases of straight and curved boundaries confirm our theoretical analysis. Finally, it should be noted that the present analysis can be extended from the BGK model to other lattice Boltzmann (LB) collision models for CDEs, which can broaden the parameter range of the relaxation time to approach 0.5.


title: Optimistic Bull or Pessimistic Bear: Adaptive Deep Reinforcement Learning for Stock Portfolio Allocation

abstract: Portfolio allocation is crucial for investment companies. However, getting the best strategy in a complex and dynamic stock market is challenging. In this paper, we propose a novel Adaptive Deep Deterministic Reinforcement Learning scheme (Adaptive DDPG) for the portfolio allocation task, which incorporates optimistic or pessimistic deep reinforcement learning that is reflected in the influence from prediction errors. Dow Jones 30 component stocks are selected as our trading stocks and their daily prices are used as the training and testing data. We train the Adaptive DDPG agent and obtain a trading strategy. The Adaptive DDPG's performance is compared with the vanilla DDPG, Dow Jones Industrial Average index and the traditional min-variance and mean-variance portfolio allocation strategies. Adaptive DDPG outperforms the baselines in terms of the investment return and the Sharpe ratio.


title: FEM-based Scalp-to-Cortex data mapping via the solution of the Cauchy problem

abstract: We propose an approach and the numerical algorithm for pre-processing of the electroencephalography (EEG) data, enabling to generate an accurate mapping of the potential from the measurement area - scalp - to the brain surface. The algorithm based on the solution of ill-posed Cauchy problem for the Laplace's equation using tetrahedral finite elements linear approximation. Application of the proposed algorithm sufficiently increases the spatial resolution of the EEG technique, making it comparable with much more complicated electrocorticography (ECoG) method.


title: Adaptive Approximate Bayesian Computation Tolerance Selection

abstract: Approximate Bayesian Computation (ABC) methods are increasingly used for inference in situations in which the likelihood function is either computationally costly or intractable to evaluate. Extensions of the basic ABC rejection algorithm have improved the computational efficiency of the procedure and broadened its applicability. The ABC-Population Monte Carlo (ABC-PMC) approach of Beaumont et al. (2009) has become a popular choice for approximate sampling from the posterior. ABC-PMC is a sequential sampler with an iteratively decreasing value of the tolerance, which specifies how close the simulated data need to be to the real data for acceptance. We propose a method for adaptively selecting a sequence of tolerances that improves the computational efficiency of the algorithm over other common techniques. In addition we define a stopping rule as a by-product of the adaptation procedure, which assists in automating termination of sampling. The proposed automatic ABC-PMC algorithm can be easily implemented and we present several examples demonstrating its benefits in terms of computational efficiency.


title: Semiclassical vibrational spectroscopy with Hessian databases

abstract: We report on a new approach to ease the computational overhead of ab initio on-the-fly semiclassical dynamics simulations for vibrational spectroscopy. The well known bottleneck of such computations lies in the necessity to estimate the Hessian matrix for propagating the semiclassical pre-exponential factor at each step along the dynamics. The procedure proposed here is based on the creation of a dynamical database of Hessians and associated molecular geometries able to speed up calculations while preserving the accuracy of results at a satisfactory level. This new approach can be interfaced to both analytical potential energy surfaces and on-the-fly dynamics, allowing one to study even large systems previously not achievable. We present results obtained for semiclassical vibrational power spectra of methane, glycine, and N-acetyl-L-phenylalaninyl-L-methionine-amide, a molecule of biological interest made of 46 atoms.


title: Best k-layer neural network approximations

abstract: We show that the empirical risk minimization (ERM) problem for neural networks has no solution in general. Given a training set $s_1, \dots, s_n \in \mathbb{R}^p$ with corresponding responses $t_1,\dots,t_n \in \mathbb{R}^q$, fitting a $k$-layer neural network $\nu_\theta : \mathbb{R}^p \to \mathbb{R}^q$ involves estimation of the weights $\theta \in \mathbb{R}^m$ via an ERM: \[ \inf_{\theta \in \mathbb{R}^m} \; \sum_{i=1}^n \lVert t_i - \nu_\theta(s_i) \rVert_2^2. \] We show that even for $k = 2$, this infimum is not attainable in general for common activations like ReLU, hyperbolic tangent, and sigmoid functions. A high-level explanation is like that for the nonexistence of best rank-$r$ approximations of higher-order tensors --- the set of parameters is not a closed set --- but the geometry involved for best $k$-layer neural networks approximations is more subtle. In addition, we show that for smooth activations $\sigma(x)= 1/\bigl(1 + \exp(-x)\bigr)$ and $\sigma(x)=\tanh(x)$, such failure to attain an infimum can happen on a positive-measured subset of responses. For the ReLU activation $\sigma(x)=\max(0,x)$, we completely classifying cases where the ERM for a best two-layer neural network approximation attains its infimum. As an aside, we obtain a precise description of the geometry of the space of two-layer neural networks with $d$ neurons in the hidden layer: it is the join locus of a line and the $d$-secant locus of a cone.


title: A New Action for Nonrelativistic Particle in Curved Background

abstract: We obtain a new form for the action of a nonrelativistic particle coupled to Newtonian gravity. The result is different from that existing in the literature which, as shown here, is riddled with problems and inconsistencies. The present derivation is based on the formalism of galilean gauge theory, introduced by us as an alternative method of analysing nonrelativistic symmetries in gravitational background.


title: Pion and kaon parton distribution functions from basis light front quantization and QCD evolution

abstract: We investigate the parton distribution functions (PDFs) of the pion and the kaon by combining quantum chromodynamics (QCD) evolution with the basis light front quantization. The initial PDFs result from the light front wave functions obtained by diagonalizing the effective Hamiltonian consisting of the holographic QCD confinement potential, a complementary longitudinal confinement potential, and the color-singlet Nambu--Jona-Lasinio interactions. The valence-quark PDF of the pion, after QCD evolution, is consistent with the result from the E-0615 experiment at Fermilab. Meanwhile, the pion structure function calculated from the PDFs agrees with the ZEUS and the H1 experiments at DESY-HERA for large $x$. Additionally, the ratio of the up quark PDF of the kaon to that of the pion is in agreement with the NA-003 experiment at CERN. We also present the cross section for the pion-nucleus induced Drell-Yan process with the obtained pion PDFs supplemented by the PDFs of the target nuclei.


title: Untangling direct and domain-mediated interactions between nicotinic acetylcholine receptors in DHA-rich membranes

abstract: At the neuromuscular junction (NMJ), the nicotinic acetylcholine receptor (nAChR) self-associates to give rise to rapid muscle movement. While lipid domains have maintained nAChR aggregates in-vitro, their specific roles in nAChR clustering are currently unknown. In the present study, we carried out coarse-grained molecular dynamics simulations (CG-MD) of 1-4 nAChR molecules in two membrane environments: One mixture containing domain-forming, homoacidic lipids, and a second mixture consisting of heteroacidic lipids. Spontaneous dimerization of nAChRs was up to ten times more likely in domain-forming membranes; however, the effect was not significant in four-protein systems, suggesting that lipid domains are less critical to nAChR oligomerization when protein concentration is higher. With regard to lipid preferences, nAChRs consistently partitioned into liquid-disordered domains occupied by the omega-3 ($\omega$-3) fatty acid, Docosahexaenoic acid (DHA); enrichment of DHA boundary lipids increased with protein concentration, particularly in homoacidic membranes. This result suggests dimer formation blocks access of saturated chains and cholesterol, but not polyunsaturated chains, to boundary lipid sites.


title: Penalized Variable Selection in Multi-Parameter Regression Survival Modelling

abstract: Multi-parameter regression (MPR) modelling refers to the approach whereby covariates are allowed to enter the model through multiple distributional parameters simultaneously. This is in contrast to the standard approaches where covariates enter through a single parameter (e.g., a location parameter). Penalized variable selection has received a significant amount of attention in recent years: methods such as the least absolute shrinkage and selection operator (LASSO), smoothly clipped absolute deviation (SCAD), and adaptive LASSO are used to simultaneously select variables and estimate their regression coefficients. Therefore, in this paper, we develop penalized multi-parameter regression methods and investigate their associated performance through simulation studies and real data; as an example, we consider the Weibull MPR model.


title: "Machine LLRning": Learning to Softly Demodulate

abstract: Soft demodulation, or demapping, of received symbols back into their conveyed soft bits, or bit log-likelihood ratios (LLRs), is at the very heart of any modern receiver. In this paper, a trainable universal neural network-based demodulator architecture, dubbed "LLRnet", is introduced. LLRnet facilitates an improved performance with significantly reduced overall computational complexity. For instance for the commonly used quadrature amplitude modulation (QAM), LLRnet demonstrates LLR estimates approaching the optimal log maximum a-posteriori inference with an order of magnitude less operations than that of the straightforward exact implementation. Link-level simulation examples for the application of LLRnet to 5G-NR and DVB-S.2 are provided. LLRnet is a (yet another) powerful example for the usefulness of applying machine learning to physical layer design.


title: Classification of Cardiac Arrhythmias from Single Lead ECG with a Convolutional Recurrent Neural Network

abstract: While most heart arrhythmias are not immediately harmful, they can lead to severe complications. In particular, atrial fibrillation, the most common arrhythmia, is characterized by fast and irregular heart beats and increases the risk of suffering a stroke. To detect such abnormal heart conditions, we propose a system composed of two main parts: a smart vest with two cooperative sensors to collect ECG data and a neural network architecture to classify heart rhythms. The smart vest uses two dry bi-electrodes to record a single lead ECG signal. The biopotential signal is then streamed via a gateway to the cloud where a neural network detects and classifies the heart arrhythmias. We selected an architecture that combines convolutional and recurrent layers. The convolutional layers extract relevant features from sliding windows of ECG and the recurrent layer aggregates them for a final softmax layer that performs the classification. Our neural network achieves an accuracy of 87.50% on the dataset of the challenge of Computing in Cardiology 2017.


title: Method of diagnosing heart disease based on deep learning ECG signal

abstract: The traditional method of diagnosing heart disease on ECG signal is artificial observation. Some have tried to combine expertise and signal processing to classify ECG signal by heart disease type. However, the currency is not so sufficient that it can be used in medical applications. We develop an algorithm that combines signal processing and deep learning to classify ECG signals into Normal AF other rhythm and noise, which help us solve this problem. It is demonstrated that we can obtain the time-frequency diagram of ECG signal by wavelet transform, and use DNN to classify the time-frequency diagram to find out the heart disease that the signal collector may have. Overall, an accuracy of 94 percent is achieved on the validation set. According to the evaluation criteria of PhysioNet/Computing in Cardiology (CinC) in 2017, the F1 score of this method is 0.957, which is higher than the first place in the competition in 2017.


title: Electroencephalogram (EEG) for Delineating Objective Measure of Autism Spectrum Disorder (ASD) (Extended Version)

abstract: Autism Spectrum Disorder (ASD) is a developmental disorder that often impairs a child's normal development of the brain. According to CDC, it is estimated that 1 in 6 children in the US suffer from development disorders, and 1 in 68 children in the US suffer from ASD. This condition has a negative impact on a person's ability to hear, socialize and communicate. Overall, ASD has a broad range of symptoms and severity; hence the term spectrum is used. One of the main contributors to ASD is known to be genetics. Up to date, no suitable cure for ASD has been found. Early diagnosis is crucial for the long-term treatment of ASD, but this is challenging due to the lack of a proper objective measures. Subjective measures often take more time, resources, and have false positives or false negatives. There is a need for efficient objective measures that can help in diagnosing this disease early as possible with less effort. EEG measures the electric signals of the brain via electrodes placed on various places on the scalp. These signals can be used to study complex neuropsychiatric issues. Studies have shown that EEG has the potential to be used as a biomarker for various neurological conditions including ASD. This chapter will outline the usage of EEG measurement for the classification of ASD using machine learning algorithms.


title: Learning for Detection: MIMO-OFDM Symbol Detection through Downlink Pilots

abstract: Reservoir computing (RC) is a special recurrent neural network which consists of a fixed high dimensional feature mapping and trained readout weights. In this paper, we introduce a new RC structure for multiple-input, multiple-output orthogonal frequency-division multiplexing (MIMO-OFDM) symbol detection, namely windowed echo state network (WESN). The theoretical analysis shows that adding buffers in input layers can bring an enhanced short-term memory (STM) to the underlying neural network. Furthermore, a unified training framework is developed for the WESN MIMO-OFDM symbol detector using both comb and scattered pilot patterns that are compatible with the structure adopted in 3GPP LTE/LTE-Advanced systems. Complexity analysis suggests the advantages of WESN based symbol detector over state-of-the-art symbol detectors such as the linear minimum mean square error (LMMSE) detection and the sphere decoder, when the system is employed with a large number of OFDM sub-carriers. Numerical evaluations illustrate the advantage of the introduced WESN-based symbol detector and demonstrate that the improvement of STM can significantly improve symbol detection performance as well as effectively mitigate model mismatch effects compared to existing methods.


title: Enabling NB-IoT on Unlicensed Spectrum

abstract: The deployment of Internet of Things (IoT) technologies in unlicensed spectrum is a candidate feature for 5G to support massive connections of IoT devices. Current IoT unlicensed band technologies, such as Sigfox and LoRa, are all at an early stage of deployment without a significant market share. In this context, the MulteFire (MF) Alliance has envisioned to adapt the cellular NB-IoT design to operate within the Sub-1 GHz unlicensed spectrum. However, the diverse regulatory requirements within this unlicensed band put a hurdle to the world-wide deployment of unlicensed band IoT technologies. To settle this challenge, MF has designed a specific framework for narrow-band (NB)-IoT systems operating on unlicensed spectrum (NB-IoT-U), which can be utilized under both the Federal Communications Commission (FCC) and European Telecommunication Standards Institute (ETSI) regulatory requirements. In this paper, enhanced synchronization and physical broadcasting signals are proposed based upon the framework designed by MF with the aim to allow a more robust detection, and to fulfil the coverage targets set for this technology. Furthermore, in order to allow the system to operate as a frequency hopping spread spectrum (FHSS) system, a novel frequency hopping pattern generator compliant with the regulatory requirements is designed, and its performance is evaluated.


title: Analytical Modeling of UAV-to-Vehicle Propagation Channels in Built-Up Areas

abstract: This letter presents an analytical path loss model for air-ground (AG) propagation between unmanned aerial vehicles (UAVs) and ground-based vehicles. We consider built-up areas, such as the ones defined by ITU-R. The three-dimensional (3D) path loss model is based on propagation conditions and essential parameters are derived by using geometric methods. Owing to the generality, the analytical model is capable of arbitrary deployments of buildings, such as suburban, urban and dense urban. The analytical model is evaluated numerically, and validations conducted by ray-tracing simulations show the high accuracy of the proposed model. The closed-form analytical formulas provide a useful tool for quick and accurate prediction of UAV-to-vehicle propagation channels.


title: Maximum-Entropy States for Magnetized Ion Transport

abstract: For a plasma with fixed total energy, number of particles, and momentum, the distribution function that maximizes entropy is a Boltzmann distribution. If, in addition, the rearrangement of charge is constrained, as happens on ion-ion collisional timescales for cross-field multiple-species transport, the maximum-entropy state is instead given by the classic impurity pinch relation. The maximum-entropy derivation, unlike previous approaches, does not rely on the details of the collision operator, only on the presence of certain conservation properties.


title: Metal Object Detection Based on Load Impedance and Input Power Characteristics for High-dimensional WPT System

abstract: High-dimensional wireless power transfer (WPT) systems have received increasing attention for charging mobile devices. With the receivers have higher spatial freedom, the systems are more susceptible to the metal objects in the surroundings. However, conventional methods for metal object detection (MOD) can't satisfy the requirements of safety and stability of new systems. This paper proposed a metal detection method which is more suitable for high-dimensional WPT systems. The key feature of the proposed method is combining load impedance and input power characteristics curves to determine whether the receiver is a metal object or a coil. And the influence of the metal is discussed in the circuit and magnetic model. The method is theoretically proven by the simulation calculation. And the experiments verified that the method can accurately detect the metal objects by setting the threshold curves.


title: Artificial Intelligence Enhances the Performance of Chaos-based Wireless Communication

abstract: Some new findings for chaos-based wireless communication systems have been identified recently. First, chaos has proven to be the optimal communication waveform because chaotic signals can achieve the maximum signal to noise ratio at receiver with the simplest matched filter. Second, the information transmitted in chaotic signals is not modified by the multipath wireless channel. Third, chaos properties can be used to relief inter-symbol interference (ISI) caused by multipath propagation. Although recent work has reported the method of obtaining the optimal threshold to eliminate the ISI in chaos-based wireless communication, its practical implementation is still a challenge. By knowing the channel parameters and all symbols, especially the future symbol to be transmitted in advance, it is almost an impossible task in the practical communication systems. Owning to Artificial intelligence (AI) recent developments, Convolutional Neural Network (CNN) with deep learning structure is being proposed to predict future symbols based on the received signal, so as to further reduce ISI and obtain better bit error rate (BER) performance as compared to that used the existing sub-optimal threshold. The feature of the method involves predicting the future symbol and obtaining a better threshold suitable for time variant channel. Numerical simulation and experimental results validate our theory and the superiority of the proposed method.


title: Tucker Tensor Decomposition on FPGA

abstract: Tensor computation has emerged as a powerful mathematical tool for solving high-dimensional and/or extreme-scale problems in science and engineering. The last decade has witnessed tremendous advancement of tensor computation and its applications in machine learning and big data. However, its hardware optimization on resource-constrained devices remains an (almost) unexplored field. This paper presents an hardware accelerator for a classical tensor computation framework, Tucker decomposition. We study three modules of this architecture: tensor-times-matrix (TTM), matrix singular value decomposition (SVD), and tensor permutation, and implemented them on Xilinx FPGA for prototyping. In order to further reduce the computing time, a warm-start algorithm for the Jacobi iterations in SVD is proposed. A fixed-point simulator is used to evaluate the performance of our design. Some synthetic data sets and a real MRI data set are used to validate the design and evaluate its performance. We compare our work with state-of-the-art software toolboxes running on both CPU and GPU, and our work shows 2.16 - 30.2x speedup on the cardiac MRI data set.


title: Deep Learning for Hybrid 5G Services in Mobile Edge Computing Systems: Learn from a Digital Twin

abstract: In this work, we consider a mobile edge computing system with both ultra-reliable and low-latency communications services and delay tolerant services. We aim to minimize the normalized energy consumption, defined as the energy consumption per bit, by optimizing user association, resource allocation, and offloading probabilities subject to the quality-of-service requirements. The user association is managed by the mobility management entity (MME), while resource allocation and offloading probabilities are determined by each access point (AP). We propose a deep learning (DL) architecture, where a digital twin of the real network environment is used to train the DL algorithm off-line at a central server. From the pre-trained deep neural network (DNN), the MME can obtain user association scheme in a real-time manner. Considering that real networks are not static, the digital twin monitors the variation of real networks and updates the DNN accordingly. For a given user association scheme, we propose an optimization algorithm to find the optimal resource allocation and offloading probabilities at each AP. Simulation results show that our method can achieve lower normalized energy consumption with less computation complexity compared with an existing method and approach to the performance of the global optimal solution.


title: Dedekind sums arising from newform Eisenstein series

abstract: For primitive non-trivial Dirichlet characters $\chi_1$ and $\chi_2$, we study the weight zero newform Eisenstein series $E_{\chi_1,\chi_2}(z,s)$ at $s=1$. The holomorphic part of this function has a transformation rule that we express in finite terms as a generalized Dedekind sum. This gives rise to the explicit construction (in finite terms) of elements of $H^1(\Gamma_0(N), \mathbb{C})$. We also give a short proof of the reciprocity formula for this Dedekind sum.


title: Digital Electronics and Analog Photonics for Convolutional Neural Networks (DEAP-CNNs)

abstract: Convolutional Neural Networks (CNNs) are powerful and highly ubiquitous tools for extracting features from large datasets for applications such as computer vision and natural language processing. However, a convolution is a computationally expensive operation in digital electronics. In contrast, neuromorphic photonic systems, which have experienced a recent surge of interest over the last few years, propose higher bandwidth and energy efficiencies for neural network training and inference. Neuromorphic photonics exploits the advantages of optical electronics, including the ease of analog processing, and busing multiple signals on a single waveguide at the speed of light. Here, we propose a Digital Electronic and Analog Photonic (DEAP) CNN hardware architecture that has potential to be 2.8 to 14 times faster while maintaining the same power usage of current state-of-the-art GPUs.

categories: eess.SP cs.NE physics.app-ph physics.optics


title: iVAMS 2.0: Machine-Learning-Metamodel-Integrated Intelligent Verilog-AMS for Fast and Accurate Mixed-Signal Design Optimization

abstract: The gap between abstraction levels in analog design is a major obstacle for advancing analog and mixed-signal (AMS) design automation and computer-aided design (CAD). Intelligent models for low-level analog building blocks are needed to bridge the accuracy gap between behavioral and transistor-level simulations. The models should be able to accurately estimate the characteristics of the analog block over a large design space. Machine learning (ML) models based on actual silicon have the capabilities of capturing detailed characteristics of complex designs. In this paper, a ML model called Artificial Neural Network Metamodels (ANNM) have been explored to capture the highly nonlinear nature of analog blocks. The application of these intelligent models to multi-objective AMS block optimization is demonstrated. Parameterized behavioral models in Verilog-AMS based on the ANN metamodels are constructed for efficient AMS design exploration. To the best of the authors' knowledge this is the first paper to integrate ANN models in Verilog-AMS, which is called iVAMS 2.0. To demonstrate the application of iVAMS 2.0, this paper presents two case studies: an operational amplifier (OP-AMP) and a phase-locked loop (PLL). A biologically-inspired "firefly optimization algorithm" is applied to an OP-AMP design in the iVAMS 2.0 framework. The optimization process is sped up by 5580X due to the use of iVAMS with negligible loss in accuracy. Similarly, for a PLL design, the physical design aware ANNs are trained and used as metamodels to predict its frequency, locking time, and power. Thorough experimental results demonstrate that only 100 sample points are sufficient for ANNs to predict the output of circuits with 21 design parameters within 3% accuracy. A proposed artificial bee colony (ABC) based algorithm performs optimization over the ANN metamodels of the PLL.


title: MuFA (Multi-type Fourier Analyzer): A tool for batch generation of MuMax3 input scripts and multi-type Fourier analysis from micromagnetic simulation output data

abstract: We present a tool for batch generation of input scripts and multi-type Fourier analysis from simulation results for the micromagnetic software MuMax3. The introduction of graphical user interface and parameter-sweeping functionality strongly speed up the input scripts creation and accelerate model optimization processes consequently. Three types of important Fourier analysis methods are provided for the acquisition of the quantitative frequency compositions, the spin-wave dispersion curve and the spatial distribution of spin-wave powers at different frequencies, respectively. Since the Fourier analysis is accelerated by parallel computations, the time cost is reduced to an acceptable level even in the presentation of tens of gigabytes data. With the MuMax3 and our proposal, a complete micromagnetic simulating tool chain from scripts generation to post analysis has been developed.


title: Deep Learned Optical Multiplexing for Multi-Focal Plane Microscopy

abstract: To obtain microscope images at multiple focal planes, the distance between the objective and sample can be mechanically adjusted. Images are acquired sequentially at each axial distance. Digital refocusing with a light-emitting diode (LED) array microscope allows elimination of this mechanical movement. In an LED array microscope, the light source of a conventional widefield microscope is replaced with a 2-dimensional LED matrix. A stack of images is acquired from the LED array microscope by sequentially illuminating each LED and capturing an image. Previous work has shown that we can achieve digital refocusing by post-processing this LED image stack. Though mechanical scanning is eliminated, digital refocusing with an LED array microscope has low temporal resolution due to the acquisition of multiple images. In this work, we propose a new paradigm for multi-focal plane microscopy for live imaging, utilizing an LED array microscope and deep learning. In our deep learning approach, we look for a single LED illumination pattern that allows the information from multiple focal planes to be multiplexed into a single image. We jointly optimize this LED illumination pattern with the parameters of a post-processing deep neural network, using a training set of LED image stacks from fixed, not live, Dugesia japonica planarians. Once training is complete, we obtain multiple focal planes by inputting a single multiplexed LED image into the trained post-processing deep neural network. We demonstrate live imaging of a D. japonica planarian at 5 focal planes with our method.


title: Experimental Demonstration of Eight-dimensional Modulation Formats for Long-haul Optical Transmission

abstract: Two novel 5.5 bit/4D modulation formats are experimentally demonstrated to transmit over 9680 km of SSMF. A reach increase of 29.4% over the 6 bit/4D PM-8QAM is shown. Furthermore, increased nonlinear tolerance with respect to PM-8QAM is shown.


title: 2D Anisotropic KPZ at stationarity: scaling, tightness and non triviality

abstract: In this work we focus on the two-dimensional anisotropic KPZ (aKPZ) equation, which is formally given by \begin{equation*}\partial_t h =\frac{\nu}{2}\Delta h + \lambda((\partial_1 h)^2 - (\partial_2 h)^2) + \nu^\frac{1}{2}\xi,\end{equation*} where $\xi$ denotes a noise which is white in both space and time, and $\lambda$ and $\nu$ are positive constants. Due to the wild oscillations of the noise and the quadratic nonlinearity, the previous equation is classically ill-posed. It is not possible to linearise it via the Cole-Hopf transformation and the pathwise techniques for singular SPDEs (the theory of Regularity Structures by M. Hairer or the paracontrolled distributions approach of M. Gubinelli, P. Imkeller, N. Perkowski) are not applicable. In the present work, we consider a regularised version of aKPZ which preserves its invariant measure. We show that in order to have subsequential limits once the regularisation is removed, it is necessary to suitably renormalise $\lambda$ and $\nu$. Moreover, we prove that, in the regime suggested by the (non-rigorous) renormalisation group computations of [D.E. Wolf, "Kinetic roughening of vicinal surfaces'', Phys. Rev. Lett., 1991], i.e. $\nu$ constant and the coupling constant $\lambda$ converging to $0$ as the inverse of the square root logarithm, any limit differs from the solution to the linear equation obtained by simply dropping the nonlinearity in aKPZ.


title: Robust and Efficient Quadrotor Trajectory Generation for Fast Autonomous Flight

abstract: In this paper, we propose a robust and efficient quadrotor motion planning system for fast flight in 3-D complex environments. We adopt a kinodynamic path searching method to find a safe, kinodynamic feasible and minimum-time initial trajectory in the discretized control space. We improve the smoothness and clearance of the trajectory by a B-spline optimization, which incorporates gradient information from a Euclidean distance field (EDF) and dynamic constraints efficiently utilizing the convex hull property of B-spline. Finally, by representing the final trajectory as a non-uniform B-spline, an iterative time adjustment method is adopted to guarantee dynamically feasible and non-conservative trajectories. We validate our proposed method in various complex simulational environments. The competence of the method is also validated in challenging real-world tasks. We release our code as an open-source package.


title: Do wind-generated waves under steady forcing propagate primarily in the downwind direction?

abstract: Measuring the directional distribution of ocean waves is a difficult task in ocean wave research. For many decades, wind-generated waves are assumed to propagate primarily in the wind direction. The concept is applied either implicitly or explicitly in the formulation of the wind input function and the design of the wave directional distribution function. Photographs of the ocean surface wave fields under steady wind forcing, however, display very different directional properties. In the near shore region, the surface undulations are dominantly crosswind. As the fetch increases, the geometry of surface waves is crosshatched, indicating two dominant wave trains crossing each other. This pattern of wave propagation in off-wind directions is repeatedly photographed in different locations under different offshore (toward) wind events. The wavelengths of these airborne photographs are in the deca-meter length scale. Shorter waves with wavelengths ranging from a few centimeters to a few decimeters also show crosshatched patterns, based on pictures taken from land and ship. In the crosshatched cases, most pictures show that the wind vector bisecting the diamond shape of the surface wave geometry. The observed wave geometries and their fetch evolution expose a major flaw in our fundamental assumption (implicitly or explicitly) that "wind-generated waves propagate primarily in the wind direction." The assumption has been taken for granted for many decades in applications including basic research of wave dynamics and the design of directional wave spectral model for calculations of the ocean surface roughness and the force, energy, and momentum of a wave field.


title: Strong coupling Bose polarons in a two-dimensional gas

abstract: We study the properties of Bose polarons in two dimensions using quantum Monte Carlo techniques. Results for the binding energy, the effective mass, and the quasiparticle residue are reported for a typical strength of interactions in the gas and for a wide range of impurity-gas coupling strengths. A lower and an upper branch of the quasiparticle exist. The lower branch corresponds to an attractive polaron and spans from the regime of weak coupling where the impurity acts as a small density perturbation of the surrounding medium to deep bound states which involve many particles from the bath and extend as far as the healing length. The upper branch corresponds to an excited state where due to repulsion a low-density bubble forms around the impurity but might be unstable against decay into many-body bound states. Interaction effects strongly affect the quasiparticle properties of the polaron. In particular, in the strongly correlated regime, the impurity features a vanishing quasiparticle residue, signaling the transition from an almost free quasiparticle to a bound state involving many atoms from the bath.


title: Three-point functions at strong coupling in the BMN limit

abstract: We consider structure constants of single-trace operators at strong coupling in planar $\mathcal{N}=4$ SYM theory using the hexagon formalism. We concentrate on heavy-heavy-light correlators where the heavy operators are BMN operators, with large R-charges and finite anomalous dimensions, and the light one is a finite-charge chiral primary operator. They describe the couplings between two highly boosted strings and a supergravity mode in the bulk dual. In the hexagon framework, two sums over virtual magnons are needed to bind the hexagons together around the light operator. We evaluate these sums explicitly at strong coupling, for a certain choice of BMN operators, and show that they factorise into a ratio of Gamma functions and a simple stringy prefactor. The former originates from giant mirror magnons scanning the AdS geometry while the latter stems from small fluctuations around the BMN vacuum. The resulting structure constants have poles at positions where an enhanced mixing with double-trace operators is expected and zeros whenever the process is forbidden by supersymmetry. We also discuss the transition to the classical regime, when the length of the light operator scales like the string tension, where we observe similitudes with the Neumann coefficients of the pp-wave String Field Theory vertex.


title: $G$-fixed Hilbert schemes on $K3$ surfaces, modular forms, and eta products

abstract: Let $X$ be a complex $K3$ surface with an effective action of a group $G$ which preserves the holomorphic symplectic form. Let $$ Z_{X,G}(q) = \sum_{n=0}^{\infty} e\left(\operatorname{Hilb}^{n}(X)^{G} \right)\, q^{n-1} $$ be the generating function for the Euler characteristics of the Hilbert schemes of $G$-invariant length $n$ subschemes. We show that its reciprocal, $Z_{X,G}(q)^{-1}$ is the Fourier expansion of a modular cusp form of weight $\frac{1}{2} e(X/G)$ for the congruence subgroup $\Gamma_{0}(|G|)$. We give an explicit formula for $Z_{X,G}$ in terms of the Dedekind eta function for all 82 possible $(X,G)$. The key intermediate result we prove is of independent interest: it establishes an eta product identity for a certain shifted theta function of the root lattice of a simply laced root system. We extend our results to various refinements of the Euler characteristic, namely the Elliptic genus, the Chi-$y$ genus, and the motivic class.


title: What, When and Where of petitions submitted to the UK Government during a time of chaos

abstract: In times marked by political turbulence and uncertainty, as well as increasing divisiveness and hyperpartisanship, Governments need to use every tool at their disposal to understand and respond to the concerns of their citizens. We study issues raised by the UK public to the Government during 2015-2017 (surrounding the UK EU-membership referendum), mining public opinion from a dataset of 10,950 petitions (representing 30.5 million signatures). We extract the main issues with a ground-up natural language processing (NLP) method, latent Dirichlet allocation (LDA). We then investigate their temporal dynamics and geographic features. We show that whilst the popularity of some issues is stable across the two years, others are highly influenced by external events, such as the referendum in June 2016. We also study the relationship between petitions' issues and where their signatories are geographically located. We show that some issues receive support from across the whole country but others are far more local. We then identify six distinct clusters of constituencies based on the issues which constituents sign. Finally, we validate our approach by comparing the petitions' issues with the top issues reported in Ipsos MORI survey data. These results show the huge power of computationally analyzing petitions to understand not only what issues citizens are concerned about but also when and from where.

categories: cs.CY cs.SI physics.data-an physics.soc-ph


title: Lorentz Violation and Radiative Corrections in Gauge Theories

abstract: Various studies have already considered radiative corrections in Lorentz-violating models unveiling many instances where a minimal or nonminimal operator generates, via loop corrections, a contribution to the photon sector of the Standard-Model Extension. However, an important fraction of this literature does not follow the widely accepted conventions and notations of the Standard-Model Extension, and this obscures the comparison between different calculations as well as possible phenomenological consequences. After reviewing some of these works, we uncover one example where a well defined loop correction to the $k_{F}$ coefficient, already presented in the literature, allows us to improve the bounds on one specific coefficient of the fermion sector of the Lorentz-violating QED extension.


title: Taint analysis of the Bitcoin network

abstract: Determining the trust of an individual Bitcoin wallet is a difficult problem. There are no ratings, that offer vendors or exchanges meaningful information about the level of the taint of Bitcoins they are receiving. Lack of such information places exchanges liable in an event when the received Bitcoins are stolen or ill-gotten. In this paper, we try to solve this problem by introducing a Bitcoin address taint score called TaintRank. It provides insight into a specific wallet by taking the addresses it interacted with throughout history into consideration. This ranking method provides such Bitcoin exchange companies insight with whom they are trading.


title: The Bredon-Landweber region in $C_2$-equivariant stable homotopy groups

abstract: We use the $C_2$-equivariant Adams spectral sequence to compute part of the $C_2$-equivariant stable homotopy groups $\pi^{C_2}_{n,n}$. This allows us to recover results of Bredon and Landweber on the image of the geometric fixed-points map from the equivariant homotopy group $\pi^{C_2}_{n,n}$ to the classical $\pi_0$. We also recover results of Mahowald and Ravenel on the Mahowald root invariants of the elements $2^k$.


title: The use of the fractal Brouers-Sotolongo formalism to analyze the kinetics of drug release

abstract: We have applied the Brouers-Sotolongo fractal kinetic equation (BSf(t,n,{\alpha})), improving notably the precision, to nine cases reported recently in the literature on drug release. The reason of using this equation is that it contains as approximations some of the mostly used empirical formula used in that field. Moreover, this equation is now successfully employed for the investigation of sorption of contaminants in aqueous media. An important extension of the BSf(t,n,{\alpha}) has been the introduction of variation of the fractal time coefficient ({\alpha}(t^{\nu} )). This improvement can lead to a greater precision of the fits and deduce some hint on the nature of the drug release process which can give precious information to propose microscopic molecular ad hoc models. We, therefore, suggest the use of the BSf(t,n,{\alpha}(t^{\nu})) formula, as a first step, in any detailed investigation and practical application of drug release data both in vitro and in vivo studies starting with the Weibull and Hill approximations to follow properly the physical solution.


title: A Column Generation Approach to the Discrete Barycenter Problem

abstract: The discrete Wasserstein barycenter problem is a minimum-cost mass transport problem for a set of discrete probability measures. Although an exact barycenter is computable through linear programming, the underlying linear program can be extremely large. For worst-case input, a best known linear programming formulation is exponential in the number of variables, but has a low number of constraints, making it an interesting candidate for column generation. In this paper, we devise and study two column generation strategies: a natural one based on a simplified computation of reduced costs, and one through a Dantzig-Wolfe decomposition. For the latter, we produce efficiently solvable subproblems, namely, a pricing problem in the form of a classical transportation problem. The two strategies begin with an efficient computation of an initial feasible solution. While the structure of the constraints leads to the computation of the reduced costs of all remaining variables for setup, both approaches may outperform a computation using the full program in speed, and dramatically so in memory requirement. In our computational experiments, we exhibit that, depending on the input, either strategy can become a best choice.


title: A fast radio burst localised to a massive galaxy

abstract: Intense, millisecond-duration bursts of radio waves have been detected from beyond the Milky Way [1]. Their extragalactic origins are evidenced by their large dispersion measures, which are greater than expected for propagation through the Milky Way interstellar medium alone, and imply contributions from the intergalactic medium and potentially host galaxies [2]. Although several theories exist for the sources of these fast radio bursts, their intensities, durations and temporal structures suggest coherent emission from highly magnetised plasma [3,4]. Two sources have been observed to repeat [5,6], and one repeater (FRB 121102) has been localised to the largest star-forming region of a dwarf galaxy at a cosmological redshift of 0.19 [7, 8]. However, the host galaxies and distances of the so far non-repeating fast radio bursts are yet to be identified. Unlike repeating sources, these events must be observed with an interferometer with sufficient spatial resolution for arcsecond localisation at the time of discovery. Here we report the localisation of a fast radio burst (FRB 190523) to a few-arcsecond region containing a single massive galaxy at a redshift of 0.66. This galaxy is in stark contrast to the host of FRB 121102, being a thousand times more massive, with a greater than hundred times lower specific star-formation rate. The properties of this galaxy highlight the possibility of a channel for FRB production associated with older stellar populations.