title: Stokes, Gibbs and volume computation of semi-algebraic sets

abstract: We consider the problem of computing the Lebesgue volume of compact basic semi-algebraic sets. In full generality, it can be approximated as closely as desired by a converging hierarchy of upper bounds obtained by applying the Moment-SOS (sums of squares) methodology to a certain infinite-dimensional linear program (LP). At each step one solves a semidefinite relaxation of the LP which involves pseudo-moments up to a certain degree. Its dual computes a polynomial of same degree which approximates from above the discontinuous indicator function of the set, hence with a typical Gibbs phenomenon which results in a slow convergence of the associated numerical scheme. Drastic improvements have been observed by introducing in the initial LP additional linear moment constraints obtained from a certain application of Stokes' theorem for integration on the set. However and so far there was no rationale to explain this behavior. We provide a refined version of this extended LP formulation. When the set is the smooth super-level set of a single polynomial, we show that the dual of this refined LP has an optimal solution which is a continuous function.Therefore in this dual one now approximates a continuous function by a polynomial, hence with no Gibbs phenomenon, which explains and improves the already observed drastic acceleration of the convergence of the hierarchy. Interestingly, the technique of proof involves recent results on Poisson's partial differential equation (PDE).


title: A formal model of Algorand smart contracts

abstract: We develop a formal model of Algorand stateless smart contracts (stateless ASC1.) We exploit our model to prove fundamental properties of the Algorand blockchain, and to establish the security of some archetypal smart contracts. While doing this, we highlight various design patterns supported by Algorand. We perform experiments to validate the coherence of our formal model w.r.t. the actual implementation.


title: Stein Variational Gaussian Processes

abstract: We show how to use Stein variational gradient descent (SVGD) to carry out inference in Gaussian process (GP) models with non-Gaussian likelihoods and large data volumes. Markov chain Monte Carlo (MCMC) is extremely computationally intensive for these situations, but the parametric assumptions required for efficient variational inference (VI) result in incorrect inference when they encounter the multi-modal posterior distributions that are common for such models. SVGD provides a non-parametric alternative to variational inference which is substantially faster than MCMC. We prove that for GP models with Lipschitz gradients the SVGD algorithm monotonically decreases the Kullback-Leibler divergence from the sampling distribution to the true posterior. Our method is demonstrated on benchmark problems in both regression and classification, a multimodal posterior, and an air quality example with 550,134 spatiotemporal observations, showing substantial performance improvements over MCMC and VI.


title: Cosmological Insights into the Early Accretion of \textit{r}-Process-Enhanced stars. I. A Comprehensive Chemo-dynamical Analysis of LAMOST J1109+0754

abstract: This study presents a comprehensive chemo-dynamical analysis of LAMOST J1109+0754, a bright (V = 12.8), extremely metal-poor (\abund{Fe}{H} = $-3.17$) star, with a strong \textit{r}-process enhancement (\abund{Eu}{Fe} = +0.94 $\pm$ 0.12). Our results are based on the 7-D measurements supplied by $Gaia$ and the chemical composition derived from a high-resolution ($R\sim 110,000$), high signal-to-noise ratio ($S/N \sim 60)$ optical spectrum obtained by the 2.4\,m Automated Planet Finder Telescope at Lick Observatory. We obtain chemical abundances of 31 elements (from lithium to thorium). The abundance ratios (\abund{X}{Fe}) of the light-elements (Z $\leqslant 30$) suggest a massive Population\,III progenitor in the 13.4-29.5\,M$_\odot$ mass range. The heavy-element ($30 <$ Z $\leqslant 90$) abundance pattern of J1109+075 agrees extremely well with the scaled-Solar \textit{r}-process signature. We have developed a novel approach to trace the kinematic history and orbital evolution of J1109+0754 with a c\textbf{O}smologically de\textbf{RI}ved tim\textbf{E}-varyi\textbf{N}g Galactic po\textbf{T}ential (the ORIENT) constructed from snapshots of a simulated Milky-Way analog taken from the \texttt{Illustris-TNG} simulation. The orbital evolution within this Milky Way-like galaxy, along with the chemical-abundance pattern implies that J1109+0754 likely originated in a low-mass dwarf galaxy located $\sim$ 60\,kpc from the center of the Galaxy, which was accreted $\sim$ 6 - 7\,Gyr ago, and that the star now belongs to the outer-halo population.


title: On the Convergence of the Multipole Expansion Method

abstract: The multipole expansion method (MEM) is a spatial discretization technique that is widely used in applications that feature scattering of waves from circular cylinders. Moreover, it also serves as a key component in several other numerical methods in which scattering computations involving arbitrarily shaped objects are accelerated by enclosing the objects in artificial cylinders. A fundamental question is that of how fast the approximation error of the MEM converges to zero as the truncation number goes to infinity. Despite the fact that the MEM was introduced in 1913, and has been in widespread usage as a numerical technique since as far back as 1955, to the best of the authors' knowledge, a precise characterization of the asymptotic rate of convergence of the MEM has not been obtained. In this work, we provide a resolution to this issue. While our focus in this paper is on the Dirichlet scattering problem, this is merely for convenience and our results actually establish convergence rates that hold for all MEM formulations irrespective of the specific boundary conditions or boundary integral equation solution representation chosen.


title: On the Graphon Mean Field Game Equations: Individual Agent Affine Dynamics and Mean Field Dependent Performance Functions

abstract: This paper establishes unique solvability of a class of Graphon Mean Field Game equations. The special case of a constant graphon yields the result for the Mean Field Game equations.


title: Direct computation of nonlinear mapping via normal form for reduced-order models of finite element nonlinear structures

abstract: The direct computation of the third-order normal form for a geometrically nonlinear structure discretised with the finite element (FE) method, is detailed. The procedure allows to define a nonlinear mapping in order to derive accurate reduced-order models (ROM) relying on invariant manifold theory. The proposed reduction strategy is direct and simulation free, in the sense that it allows to pass from physical coordinates (FE nodes) to normal coordinates, describing the dynamics in an invariant-based span of the phase space. The number of master modes for the ROM is not a priori limited since a complete change of coordinate is proposed. The underlying theory ensures the quality of the predictions thanks to the invariance property of the reduced subspace, together with their curvatures in phase space that accounts for the nonresonant nonlinear couplings. The method is applied to a beam discretised with 3D elements and shows its ability in recovering internal resonance at high energy. Then a fan blade model is investigated and the correct prediction given by the ROMs are assessed and discussed. A method is proposed to approximate an aggregate value for the damping, that takes into account the damping coefficients of all the slave modes, and also using the Rayleigh damping model as input. Frequency-response curves for the beam and the blades are then exhibited, showing the accuracy of the proposed method.


title: GEFA: Early Fusion Approach in Drug-Target Affinity Prediction

abstract: Predicting the interaction between a compound and a target is crucial for rapid drug repurposing. Deep learning has been successfully applied in drug-target affinity (DTA) problem. However, previous deep learning-based methods ignore modeling the direct interactions between drug and protein residues. This would lead to inaccurate learning of target representation which may change due to the drug binding effects. In addition, previous DTA methods learn protein representation solely based on a small number of protein sequences in DTA datasets while neglecting the use of proteins outside of the DTA datasets. We propose GEFA (Graph Early Fusion Affinity), a novel graph-in-graph neural network with attention mechanism to address the changes in target representation because of the binding effects. Specifically, a drug is modeled as a graph of atoms, which then serves as a node in a larger graph of residues-drug complex. The resulting model is an expressive deep nested graph neural network. We also use pre-trained protein representation powered by the recent effort of learning contextualized protein representation. The experiments are conducted under different settings to evaluate scenarios such as novel drugs or targets. The results demonstrate the effectiveness of the pre-trained protein embedding and the advantages our GEFA in modeling the nested graph for drug-target interaction.


title: Many-body collisional dynamics of impurities injected into a double-well trapped Bose-Einstein condensate

abstract: We unravel the many-body dynamics of a harmonically trapped impurity colliding with a bosonic medium confined in a double-well upon quenching the initially displaced harmonic trap to the center of the double-well. We reveal that the emerging correlation dynamics crucially depends on the impurity-medium interaction strength allowing for a classification into different dynamical response regimes. For strong attractive impurity-medium couplings the impurity is bound to the bosonic bath, while for intermediate attractions it undergoes an effective tunneling. In the case of weak attractive or repulsive couplings the impurity penetrates the bosonic bath and performs a dissipative oscillatory motion. Further increasing the impurity-bath repulsion results in the pinning of the impurity between the density peaks of the bosonic medium, a phenomenon that is associated with a strong impurity-medium entanglement. For strong repulsions the impurity is totally reflected by the bosonic medium. To unravel the underlying microscopic excitation processes accompanying the dynamics we employ an effective potential picture. We extend our results to the case of two bosonic impurities and demonstrate the existence of a qualitatively similar impurity dynamics.


title: Adaptive Multi-modal Fusion Hashing via Hadamard Matrix

abstract: Hashing plays an important role in information retrieval, due to its low storage and high speed of processing. Among the techniques available in the literature, multi-modal hashing, which can encode heterogeneous multi-modal features into compact hash codes, has received particular attention. Most of the existing multi-modal hashing methods adopt the fixed weighting factors to fuse multiple modalities for any query data, which cannot capture the variation of different queries. Besides, many methods introduce hyper-parameters to balance many regularization terms that make the optimization harder. Meanwhile, it is time-consuming and labor-intensive to set proper parameter values. The limitations may significantly hinder their promotion in real applications. In this paper, we propose a simple, yet effective method that is inspired by the Hadamard matrix. The proposed method captures the multi-modal feature information in an adaptive manner and preserves the discriminative semantic information in the hash codes. Our framework is flexible and involves a very few hyper-parameters. Extensive experimental results show the method is effective and achieves superior performance compared to state-of-the-art algorithms.


title: Can neutron disappearance/reappearance experiments definitively rule out the existence of hidden braneworlds endowed with a copy of the Standard Model?

abstract: Many works, aiming to explain the origin of dark matter or dark energy, consider the existence of hidden (brane)worlds parallel to our own visible world - our usual universe - in a multidimensional bulk. Hidden braneworlds allow for hidden copies of the Standard Model. For instance, atoms hidden in a hidden brane could exist as dark matter candidates. As a way to constrain such hypotheses, the possibility for neutron-hidden neutron swapping can be tested thanks to disappearance-reappearance experiments also known as passing-through-walls neutron experiments. The neutron-hidden neutron coupling $g$ can be constrained from those experiments. While $g$ could be arbitrarily small, previous works involving a $M_4 \times R_1$ bulk, with DGP branes, show that $g$ then possesses a value which is reachable experimentally. It is of crucial interest to know if a reachable value for $g$ is universal or not and to estimate its magnitude. Indeed, it would allow, in a near future, to reject definitively - or not - the existence of hidden braneworlds from experiments. In the present paper, we explore this issue by calculating $g$ for DGP branes, for $M_4 \times S_1/Z_2$, $M_4 \times R_2$ and $M_4 \times T^2$ bulks. As a major result, no disappearance-reappearance experiment would definitively universally rules out the existence of hidden worlds endowed with their own copy of Standard Model particles, excepted for specific scenarios with conditions reachable in future experiments.


title: Investigation of the Effect of Fear and Stress on Password Choice (Extended Version)

abstract: Background. The current cognitive state, such as cognitive effort and depletion, incidental affect or stress may impact the strength of a chosen password unconsciously. Aim. We investigate the effect of incidental fear and stress on the measured strength of a chosen password. Method. We conducted two experiments with within-subject designs measuring the Zxcvbn \textsf{log10} number of guesses as strength of chosen passwords as dependent variable. In both experiments, participants were signed up to a site holding their personal data and, for the second run a day later, asked under a security incident pretext to change their password. (a) Fear. $N_\mathsf{F} = 34$ participants were exposed to standardized fear and happiness stimulus videos in random order. (b) \textbf{Stress.} $N_\mathsf{S} = 50$ participants were either exposed to a battery of standard stress tasks or left in a control condition in random order. The Zxcvbn password strength was compared across conditions. Results. We did not observe a statistically significant difference in mean Zxcvbn password strengths on fear (Hedges' $g_{\mathsf{av}} = -0.11$, 95\% CI $[-0.45, 0.23]$) or stress (and control group, Hedges' $g_{\mathsf{av}} = 0.01$, 95\% CI $[-0.31, 0.33]$). However, we found a statistically significant cross-over interaction of stress and TLX mental demand. Conclusions. While having observed negligible main effect size estimates for incidental fear and stress, we offer evidence towards the interaction between stress and cognitive effort that vouches for further investigation.


title: Two-dimensional simulated tempering for the isobaric-isothermal ensemble with fast on-the-fly weight determination

abstract: We propose a method to extend the fast on-the-fly weight determination scheme for simulated tempering to two-dimensional space including not only temperature but also pressure. During the simulated tempering simulation, weight parameters for temperature-update and pressure-update are self-updated independently according to the trapezoidal rule. In order to test the effectiveness of the algorithm, we applied our proposed method to a peptide, chignolin, in explicit water. After setting all weight parameters to zero, the weight parameters were quickly determined during the simulation. The simulation realised a uniform random walk in the entire temperature-pressure space.

categories: physics.comp-ph cond-mat.stat-mech physics.bio-ph physics.chem-ph


title: Ethical conceptual replication of visualization research considering sources of methodological bias and practical significance

abstract: General design principles for visualization have been relatively well-established based on a combination of cognitive and perceptual theory and empirical evaluations over the past 20 years. To determine how these principles hold up across use contexts and end-users, I argue that we should emphasize conceptual replication focused on determining practical significance and reducing methodological biases. This shift in thinking aims to determine how design principles interact with methodological approaches, laying the groundwork for visualization meta-science.


title: A Systematic Review on Model Watermarking for Neural Networks

abstract: Machine learning (ML) models are applied in an increasing variety of domains. The availability of large amounts of data and computational resources encourages the development of ever more complex and valuable models. These models are considered intellectual property of the legitimate parties who have trained them, which makes their protection against stealing, illegitimate redistribution, and unauthorized application an urgent need. Digital watermarking presents a strong mechanism for marking model ownership and, thereby, offers protection against those threats. This work presents a taxonomy identifying and analyzing different classes of watermarking schemes for ML models. It introduces a unified threat model to allow structured reasoning on and comparison of the effectiveness of watermarking methods in different scenarios. Furthermore, it systematizes desired security requirements and attacks against ML model watermarking. Based on that framework, representative literature from the field is surveyed to illustrate the taxonomy. Finally, shortcomings and general limitations of existing approaches are discussed, and an outlook on future research directions is given.


title: Integration of Formal Proof into Unified Assurance Cases with Isabelle/SACM

abstract: Assurance cases are often required to certify critical systems. The use of formal methods in assurance can improve automation, increase confidence, and overcome errant reasoning. However, assurance cases can never be fully formalised, as the use of formal methods is contingent on models that are validated by informal processes. Consequently, assurance techniques should support both formal and informal artifacts, with explicated inferential links between them. In this paper, we contribute a formal machine-checked interactive language, called Isabelle/SACM, supporting the computer-assisted construction of assurance cases compliant with the OMG Structured Assurance Case Meta-Model. The use of Isabelle/SACM guarantees well-formedness, consistency, and traceability of assurance cases, and allows a tight integration of formal and informal evidence of various provenance. In particular, Isabelle brings a diverse range of automated verification techniques that can provide evidence. To validate our approach, we present a substantial case study based on the Tokeneer secure entry system benchmark. We embed its functional specification into Isabelle, verify its security requirements, and form a modular security case in Isabelle/SACM that combines the heterogeneous artifacts. We thus show that Isabelle is a suitable platform for critical systems assurance.


title: A Decade of Evidence of Trend Following Investing in Cryptocurrencies

abstract: Cryptocurrency markets have many of the characteristics of 20th century commodities markets, making them an attractive candidate for trend following strategies. We present a decade of evidence from the infancy of bitcoin, showcasing the potential investor returns in cryptocurrency trend following, 255% walkforward annualised returns. We find that cryptocurrencies offer similar returns characteristics to commodities with similar risk-adjusted returns, and strong bear market diversification against traditional equities. Code available at https://github.com/Globe-Research/bittrends.


title: An Empirical Study on the Impact of Deep Parameters on Mobile App Energy Usage

abstract: Improving software performance through configuration parameter tuning is a common activity during software maintenance. Beyond traditional performance metrics like latency, mobile app developers are interested in reducing app energy usage. Some mobile apps have centralized locations for parameter tuning, similar to databases and operating systems, but it is common for mobile apps to have hundreds of parameters scattered around the source code. The correlation between these "deep" parameters and app energy usage is unclear. Researchers have studied the energy effects of deep parameters in specific modules, but we lack a systematic understanding of the energy impact of mobile deep parameters. In this paper we empirically investigate this topic, combining a developer survey with systematic energy measurements. Our motivational survey of 25 Android developers suggests that developers do not understand, and largely ignore, the energy impact of deep parameters. To assess the potential implications of this practice, we propose a deep parameter energy profiling framework that can analyze the energy impact of deep parameters in an app. Our framework identifies deep parameters, mutates them based on our parameter value selection scheme, and performs reliable energy impact analysis. Applying the framework to 16 popular Android apps, we discovered that deep parameter-induced energy inefficiency is rare. We found only 2 out of 1644 deep parameters for which a different value would significantly improve its app's energy efficiency. A detailed analysis found that most deep parameters have either no energy impact, limited energy impact, or an energy impact only under extreme values. Our study suggests that it is generally safe for developers to ignore the energy impact when choosing deep parameter values in mobile apps.


title: SOUP: Spatial-Temporal Demand Forecasting and Competitive Supply

abstract: We consider a setting with an evolving set of requests for transportation from an origin to a destination before a deadline and a set of agents capable of servicing the requests. In this setting, an assignment authority is to assign agents to requests such that the average idle time of the agents is minimized. An example is the scheduling of taxis (agents) to meet incoming requests for trips while ensuring that the taxis are empty as little as possible. In this paper, we study the problem of spatial-temporal demand forecasting and competitive supply (SOUP). We address the problem in two steps. First, we build a granular model that provides spatial-temporal predictions of requests. Specifically, we propose a Spatial-Temporal Graph Convolutional Sequential Learning (ST-GCSL) algorithm that predicts the service requests across locations and time slots. Second, we provide means of routing agents to request origins while avoiding competition among the agents. In particular, we develop a demand-aware route planning (DROP) algorithm that considers both the spatial-temporal predictions and the supplydemand state. We report on extensive experiments with realworld and synthetic data that offer insight into the performance of the solution and show that it is capable of outperforming the state-of-the-art proposals.


title: Hadronic X-ray Flares from Blazars

abstract: The detection of a high-energy neutrino from the flaring blazar TXS 0506+056 and the subsequent discovery of a neutrino excess from the same direction have strengthened the hypothesis that blazars are cosmic neutrino sources. The lack, however, of $\gamma$-ray flaring activity during the latter period challenges the standard scenario of correlated $\gamma$-ray and high-energy neutrino emission in blazars. We propose instead that TeV-PeV neutrinos are produced in coincidence with X-ray flares that are powered by proton synchrotron radiation. In this case, neutrinos are produced by photomeson interactions of protons with their own synchrotron radiation, while MeV to GeV $\gamma$-rays are the result of synchrotron-dominated electromagnetic cascades developed in the source. Using a time-dependent approach, we find that this "pure hadronic flaring" hypothesis has several interesting consequences. The X-ray flux is a good proxy for the all-flavor neutrino flux, while certain neutrino-rich X-ray flares may be dark in GeV-TeV $\gamma$-rays. Lastly, hadronic X-ray flares are accompanied by an equally bright MeV component that is detectable by proposed missions like e-ASTROGAM and AMEGO. We then applied this scenario to the extreme blazar 3HSP J095507.9+355101 that has been associated with IceCube-200107A while undergoing an X-ray flare. We showed that the number of muon and antimuon neutrinos above 100 TeV during hadronic flares can be up to $\sim3-10$ times higher than the expected number in standard leptohadronic models. Still, frequent hadronic flaring activity is necessary for explaining the detected neutrino event IceCube-200107A.


title: p-Determinants and monodromy of differential operators

abstract: We prove that $p$-determinants of a certain class of differential operators can be lifted to power series over $\mathbb{Q}$. We compute these power series in terms of monodromy of the corresponding differential operators.


title: Higher-order contact mechanics

abstract: We present a complete theory of higher-order autonomous contact mechanics, which allows us to describe higher-order dynamical systems with dissipation. The essential tools for the theory are the extended higher-order tangent bundles, ${\rm T}^kQ\times{\mathbb R}$, whose geometric structures are previously introduced in order to state the Lagrangian and Hamiltonian formalisms for these kinds of systems, including their variational formulation. The variational principle, the contact forms, and the geometric dynamical equations are obtained by using those structures and generalizing the standard formulation of contact Lagrangian and Hamiltonian systems. As an alternative approach, we develop a unified description that encompasses the Lagrangian and Hamiltonian equations as well as their relationship through the Legendre map; all of them are obtained from the contact dynamical equations and the constraint algorithm that is implemented because, in this formalism, the dynamical systems are always singular. Some interesting examples are finally analyzed using these geometric formulations.


title: Multimessenger Constraints on Intergalactic Magnetic Fields from the Flare of TXS 0506+056

abstract: The origin of magnetic fields in the universe is an open problem. Seed magnetic fields possibly produced in early times may have survived up to the present day close to their original form, providing an untapped window to the primeval universe. The recent observations of high-energy neutrinos from the blazar TXS 0506+056 in association with an electromagnetic counterpart in a broad range of wavelengths can be used to probe intergalactic magnetic fields via the time delay between the neutrinos and gamma rays as well as the time dependence of the gamma-ray fluxes. Using extensive three-dimensional Monte Carlo simulations, we present a novel method to constrain these fields. We apply it to TXS 0506+056 and, for the first time, derive constraints on both the magnetic-field strength and its coherence length, considering six orders of magnitude for each.


title: Active particles with polar alignment in ring-shaped confinement

abstract: We study steady-state properties of a suspension of active, nonchiral and chiral, Brownian particles with polar alignment and steric interactions confined within a ring-shaped (annulus) confinement in two dimensions. Exploring possible interplays between polar interparticle alignment, geometric confinement and the surface curvature, being incorporated here on minimal levels, we report a surface-population reversal effect, whereby active particles migrate from the outer concave boundary of the annulus to accumulate on its inner convex boundary. This contrasts the conventional picture, implying stronger accumulation of active particles on concave boundaries relative to the convex ones. The population reversal is caused by both particle alignment and surface curvature, disappearing when either of these factors is absent. We explore the ensuing consequences for the chirality-induced current and swim pressure of active particles and analyze possible roles of system parameters, such as the mean number density of particles and particle self-propulsion, chirality and alignment strengths.


title: First-principles modeling of plasmons in aluminum under ambient and extreme conditions

abstract: The theoretical understanding of plasmon behavior is crucial for an accurate interpretation of inelastic scattering diagnostics in many experiments. We highlight the utility of linear-response time-dependent density functional theory (LR-TDDFT) as a first-principles framework for consistently modeling plasmon properties. We provide a comprehensive analysis of plasmons in aluminum from ambient to warm dense matter conditions and assess typical properties such as the dynamical structure factor, the plasmon dispersion, and the plasmon lifetime. We compare our results with scattering measurements and with other TDDFT results as well as models such as the random phase approximation, the Mermin approach, and the dielectric function obtained using static local field corrections of the uniform electron gas parametrized from path integral Monte Carlo simulations. We conclude that results for the plasmon dispersion and lifetime are inconsistent between experiment and theories and that the common practice of extracting and studying plasmon dispersion relations is an insufficient procedure to capture the complicated physics contained in the dynamic structure factor in its full breadth.


title: A Comprehensive Survey of the Tactile Internet: State of the art and Research Directions

abstract: The Internet has made several giant leaps over the years, from a fixed to a mobile Internet, then to the Internet of Things, and now to a Tactile Internet. The Tactile Internet goes far beyond data, audio and video delivery over fixed and mobile networks, and even beyond allowing communication and collaboration among things. It is expected to enable haptic communication and allow skill set delivery over networks. Some examples of potential applications are tele-surgery, vehicle fleets, augmented reality and industrial process automation. Several papers already cover many of the Tactile Internet-related concepts and technologies, such as haptic codecs, applications, and supporting technologies. However, none of them offers a comprehensive survey of the Tactile Internet, including its architectures and algorithms. Furthermore, none of them provides a systematic and critical review of the existing solutions. To address these lacunae, we provide a comprehensive survey of the architectures and algorithms proposed to date for the Tactile Internet. In addition, we critically review them using a well-defined set of requirements and discuss some of the lessons learned as well as the most promising research directions.


title: Integration of Roadside Camera Images and Weather Data for Monitoring Winter Road Surface Conditions

abstract: During the winter season, real-time monitoring of road surface conditions is critical for the safety of drivers and road maintenance operations. Previous research has evaluated the potential of image classification methods for detecting road snow coverage by processing images from roadside cameras installed in RWIS (Road Weather Information System) stations. However, there are a limited number of RWIS stations across Ontario, Canada; therefore, the network has reduced spatial coverage. In this study, we suggest improving performance on this task through the integration of images and weather data collected from the RWIS stations with images from other MTO (Ministry of Transportation of Ontario) roadside cameras and weather data from Environment Canada stations. We use spatial statistics to quantify the benefits of integrating the three datasets across Southern Ontario, showing evidence of a six-fold increase in the number of available roadside cameras and therefore improving the spatial coverage in the most populous ecoregions in Ontario. Additionally, we evaluate three spatial interpolation methods for inferring weather variables in locations without weather measurement instruments and identify the one that offers the best tradeoff between accuracy and ease of implementation.


title: Lorentz Resonance in the Homogenization of Plasmonic Crystals

abstract: We explain the Lorentz resonances in plasmonic crystals that consist of 2D nano dielectric inclusions as the interaction between resonant material properties and geometric resonances of electrostatic nature. One example of such plasmonic crystals are graphene nanosheets that are periodically arranged within a non-magnetic bulk dielectric. We identify local geometric resonances on the length scale of the small scale period. From a materials perspective, the graphene surface exhibits a dispersive surface conductance captured by the Drude model. Together these phenomena conspire to generate Lorentz resonances at frequencies controlled by the surface geometry and the surface conductance. The Lorentz resonances found in the frequency response of the effective dielectric tensor of the bulk metamaterial is shown to be given by an explicit formula, in which material properties and geometric resonances are decoupled. This formula is rigorous and obtained directly from corrector fields describing local electrostatic fields inside the heterogeneous structure. Our analytical findings can serve as an efficient computational tool to describe the general frequency dependence of periodic optical devices. As a concrete example, we investigate two prototypical geometries composed of nanotubes and nanoribbons.


title: Vertical Power Flow Forecast with LSTMs Using Regular Training Update Strategies

abstract: The strong growth of renewable energy sources and the high volatility in power generation of these sources, as well as the increasing amount of volatile energy consumption is leading to major challenges in the electrical grid. In order to ensure safety and reliability in the electricity grid, the power flow in the grid needs to be observed to prevent overloading. Furthermore, the energy supply and consumption need to be continuously balanced to ensure the security of energy supply. Therefore a high quality of power flow forecasts for the next few hours within the grid are needed. In this paper we investigate forecasts of the vertical power flow at transformer between the medium and high voltage grid. Forecasting the vertical power flow is challenging due to constantly changing characteristics of the power flow at the transformer. This is mainly a result of dynamic grid topologies, changes in the installed assets, maintenance of the transformer itself as well as the volatile generation. In this paper we present a novel approach to deal with these challenges. For the multi step time series forecasts a Long-Short Term Memory (LSTM) is used. In our presented approach an update process where the model is retrained regularly is investigated and compared to baseline models. The model is retrained as soon as a sufficient amount of new measurements are available. This retraining should capture changes in the characteristic of the transformer that the model has not yet seen in the past and therefore cannot be predicted by the model. For the regular update process we investigate different strategies where especially the number of used epochs are considered, but also different learning rates are used. We show that our new approach significantly outperforms the investigated baseline approaches.


title: Transient Classification in low SNR Gravitational Wave data using Deep Learning

abstract: The recent advances in Gravitational-wave astronomy have greatly accelerated the study of Multimessenger astrophysics. There is a need for the development of fast and efficient algorithms to detect non-astrophysical transients and noises due to the rate and scale at which the data is being provided by LIGO and other gravitational wave observatories. These transients and noises can interfere with the study of gravitational waves and binary mergers and induce false positives. Here, we propose the use of deep learning algorithms to detect and classify these transient signals. Traditional statistical methods are not well designed for dealing with temporal signals but supervised deep learning techniques such as RNN-LSTM and deep CNN have proven to be effective for solving problems such as time-series forecasting and time-series classification. We also use unsupervised models such as Total variation, Principal Component Analysis, Support Vector Machine, Wavelet decomposition or Random Forests for feature extraction and noise reduction and then study the results obtained by RNN-LSTM and deep CNN for classifying the transients in low-SNR signals. We compare the results obtained by the combination of various unsupervised models and supervised models. This method can be extended to real-time detection of transients and merger signals using deep-learning optimized GPU's for early prediction and study of various astronomical events. We will also explore and compare other machine learning models such as MLP, Stacked Autoencoder, Random forests, extreme learning machine, Support Vector machine and logistic regression classifier.


title: Thermodynamically consistent nonlinear viscoplastic formulation with well-conditioned recovery of the inviscid solution: Theory and implicit integration algorithm with exact solution for the linear case

abstract: In this work, a consistent viscoplasticity formulation is derived from thermodynamical principles and employing the concept of continuum elastic corrector rate. The proposed model is developed based on the principle of maximum viscoplastic dissipation for determining the flow direction. The model uses both the equivalent viscoplastic strain and its rate as state variables. Power balance and energy balance give, respectively, separate evolution equations for the equivalent viscoplastic strain rate and the viscoplastic strain, the former written in terms of inviscid rates. Several key points distinguish our formulation from other proposals. First, the viscoplastic strain rate (instead of a yield function) consistently distinguishes conservative from dissipative behaviours during reverse loading; and the discrete implicit integration algorithm is an immediate implementation of the continuum theory based on the mentioned principles. Second, the inviscid solution is recovered in a well-conditioned manner by simply setting the viscosity to zero. Indeed, inviscid plasticity, viscoelasticity and viscoplasticity are particular cases of our formulation and integration algorithm, and are recovered just by setting the corresponding parameters to zero (viscosity or yield stress). Third, the linear viscoplasticity solution is obtained in an exact manner for proportional loading cases, independently of the time step employed. Four, general nonlinear models (Perzyna, Norton, etc) may be immediately incorporated as particular cases both in the theory and the computational implementation.


title: Delay Characterization of Mobile Edge Computing for 6G Time-Sensitive Services

abstract: Time-sensitive services (TSSs) have been widely envisioned for future sixth generation (6G) wireless communication networks. Due to its inherent low-latency advantage, mobile edge computing (MEC) will be an indispensable enabler for TSSs. The random characteristics of the delay experienced by users are key metrics reflecting the quality of service (QoS) of TSSs. Most existing studies on MEC have focused on the average delay. Only a few research efforts have been devoted to other random delay characteristics, such as the delay bound violation probability and the probability distribution of the delay, by decoupling the transmission and computation processes of MEC. However, if these two processes could not be decoupled, the coupling will bring new challenges to analyzing the random delay characteristics. In this paper, an MEC system with a limited computation buffer at the edge server is considered. In this system, the transmission process and computation process form a feedback loop and could not be decoupled. We formulate a discrete-time two-stage tandem queueing system. Then, by using the matrix-geometric method, we obtain the estimation methods for the random delay characteristics, including the probability distribution of the delay, the delay bound violation probability, the average delay and the delay standard deviation. The estimation methods are verified by simulations. The random delay characteristics are analyzed by numerical experiments, which unveil the coupling relationship between the transmission process and computation process for MEC. These results will largely facilitate elaborate allocation of communication and computation resources to improve the QoS of TSSs.


title: COVID-19: Forecasting mortality given mobility trend data and non-pharmaceutical interventions

abstract: We develop a novel hybrid epidemiological model and a specific methodology for its calibration to distinguish and assess the impact of mobility restrictions (given by Apple's mobility trends data) from other complementary non-pharmaceutical interventions (NPIs) used to control the spread of COVID-19. Using the calibrated model, we estimate that mobility restrictions contribute to 47 % (US States) and 47 % (worldwide) of the overall suppression of the disease transmission rate using data up to 13/08/2020. The forecast capacity of our model was evaluated doing four-weeks ahead predictions. Using data up to 30/06/20 for calibration, the mean absolute percentage error (MAPE) of the prediction of cumulative deceased individuals was 5.0 % for the United States (51 states) and 6.7 % worldwide (49 countries). This MAPE was reduced to 3.5% for the US and 3.8% worldwide using data up to 13/08/2020. We find that the MAPE was higher for the total confirmed cases at 11.5% worldwide and 10.2% for the US States using data up to 13/08/2020. Our calibrated model achieves an average R-Squared value for cumulative confirmed and deceased cases of 0.992 using data up to 30/06/20 and 0.98 using data up to 13/08/20.


title: Realisability for Infinitary Intuitionistic Set Theory

abstract: We introduce a realisability semantics for infinitary intuitionistic set theory that is based on Ordinal Turing Machines (OTMs). We show that our notion of OTM-realisability is sound with respect to certain systems of infinitary intuitionistic logic, and that all axioms of infinitary Kripke-Platek set theory are realised. Finally, we use a variant of our notion of realisability to show that the propositional admissible rules of (finitary) intuitionistic Kripke-Platek set theory are exactly the admissible rules of intuitionistic propositional logic.


title: Sharp Sobolev estimates for concentration of solutions to an aggregation-diffusion equation

abstract: We consider the drift-diffusion equation $u_t-\epsilon\Delta u + \nabla \cdot(u\nabla K^*u)=0$ in the whole space with global-in-time solutions bounded in all Sobolev spaces; for simplicity, we restrict ourselves to the model case $K(x)=-|x|$. We quantify the mass concentration phenomenon, a genuinely nonlinear effect, for radially symmetric solutions of this equation for small diffusivity $\epsilon$ studied in our previous paper [3], obtaining optimal sharp upper and lower bounds for Sobolev norms.


title: Goal-Directed Occupancy Prediction for Lane-Following Actors

abstract: Predicting the possible future behaviors of vehicles that drive on shared roads is a crucial task for safe autonomous driving. Many existing approaches to this problem strive to distill all possible vehicle behaviors into a simplified set of high-level actions. However, these action categories do not suffice to describe the full range of maneuvers possible in the complex road networks we encounter in the real world. To combat this deficiency, we propose a new method that leverages the mapped road topology to reason over possible goals and predict the future spatial occupancy of dynamic road actors. We show that our approach is able to accurately predict future occupancy that remains consistent with the mapped lane geometry and naturally captures multi-modality based on the local scene context while also not suffering from the mode collapse problem observed in prior work.


title: Intelligent Risk Alarm for Asthma Patients using Artificial Neural Networks

abstract: Asthma is a chronic disease of the airways of the lungs. It results in inflammation and narrowing of the respiratory passages, which prevents air flow into the airways and leads to frequent bouts of shortness of breath with wheezing accompanied by coughing and phlegm after exposure to inhalation of substances that provoke allergic reactions or irritation of the respiratory system. Data mining in healthcare system is very important in diagnosing and understanding data, so data mining aims to solve basic problems in diagnosing diseases due to the complexity of diagnosing asthma. Predicting chemicals in the atmosphere is very important and one of the most difficult problems since the last century. In this paper, the impact of chemicals on asthma patient will be presented and discussed. Sensor system called MQ5 will be used to examine the smoke and nitrogen content in the atmosphere. MQ5 will be inserted in a wristwatch that checks the smoke and nitrogen content in the patients place, the system shall issue a warning alarm if this gas affects the person with asthma. It will be based on the Artificial Neural Networks (ANN) algorithm that has been built using data that containing a set of chemicals such as carbon monoxide, NMHC (GT) acid gas, C6H6 (GT) Gasoline, NOx (GT) Nitrogen Oxide, and NO2 (GT) Nitrogen Dioxide. The temperature and humidity will be also used as they can negatively affect asthma patient. Finally, the rating model was evaluated and achieved 99.58% classification accuracy.


title: COVID-19 Pandemic Prediction using Time Series Forecasting Models

abstract: Millions of people have been infected and lakhs of people have lost their lives due to the worldwide ongoing novel Coronavirus (COVID-19) pandemic. It is of utmost importance to identify the future infected cases and the virus spread rate for advance preparation in the healthcare services to avoid deaths. Accurately forecasting the spread of COVID-19 is an analytical and challenging real-world problem to the research community. Therefore, we use day level information of COVID-19 spread for cumulative cases from whole world and 10 mostly affected countries; US, Spain, Italy, France, Germany, Russia, Iran, United Kingdom, Turkey, and India. We utilize the temporal data of coronavirus spread from January 22, 2020 to May 20, 2020. We model the evolution of the COVID-19 outbreak, and perform prediction using ARIMA and Prophet time series forecasting models. Effectiveness of the models are evaluated based on the mean absolute error, root mean square error, root relative squared error, and mean absolute percentage error. Our analysis can help in understanding the trends of the disease outbreak, and provide epidemiological stage information of adopted countries. Our investigations show that ARIMA model is more effective for forecasting COVID-19 prevalence. The forecasting results have potential to assist governments to plan policies to contain the spread of the virus.


title: Tarsier: Evolving Noise Injection in Super-Resolution GANs

abstract: Super-resolution aims at increasing the resolution and level of detail within an image. The current state of the art in general single-image super-resolution is held by NESRGAN+, which injects a Gaussian noise after each residual layer at training time. In this paper, we harness evolutionary methods to improve NESRGAN+ by optimizing the noise injection at inference time. More precisely, we use Diagonal CMA to optimize the injected noise according to a novel criterion combining quality assessment and realism. Our results are validated by the PIRM perceptual score and a human study. Our method outperforms NESRGAN+ on several standard super-resolution datasets. More generally, our approach can be used to optimize any method based on noise injection.


title: Do We Really Sample Right In Model-Based Diagnosis?

abstract: Statistical samples, in order to be representative, have to be drawn from a population in a random and unbiased way. Nevertheless, it is common practice in the field of model-based diagnosis to make estimations from (biased) best-first samples. One example is the computation of a few most probable possible fault explanations for a defective system and the use of these to assess which aspect of the system, if measured, would bring the highest information gain. In this work, we scrutinize whether these statistically not well-founded conventions, that both diagnosis researchers and practitioners have adhered to for decades, are indeed reasonable. To this end, we empirically analyze various sampling methods that generate fault explanations. We study the representativeness of the produced samples in terms of their estimations about fault explanations and how well they guide diagnostic decisions, and we investigate the impact of sample size, the optimal trade-off between sampling efficiency and effectivity, and how approximate sampling techniques compare to exact ones.


title: Improved Dimensionality Reduction of various Datasets using Novel Multiplicative Factoring Principal Component Analysis (MPCA)

abstract: Principal Component Analysis (PCA) is known to be the most widely applied dimensionality reduction approach. A lot of improvements have been done on the traditional PCA, in order to obtain optimal results in the dimensionality reduction of various datasets. In this paper, we present an improvement to the traditional PCA approach called Multiplicative factoring Principal Component Analysis (MPCA). The advantage of MPCA over the traditional PCA is that a penalty is imposed on the occurrence space through a multiplier to make negligible the effect of outliers in seeking out projections. Here we apply two multiplier approaches, total distance and cosine similarity metrics. These two approaches can learn the relationship that exists between each of the data points and the principal projections in the feature space. As a result of this, improved low-rank projections are gotten through multiplying the data iteratively to make negligible the effect of corrupt data in the training set. Experiments were carried out on YaleB, MNIST, AR, and Isolet datasets and the results were compared to results gotten from some popular dimensionality reduction methods such as traditional PCA, RPCA-OM, and also some recently published methods such as IFPCA-1 and IFPCA-2.


title: Fast computation of hyperelliptic curve isogenies in odd characteristic

abstract: Let p be an odd prime number and g $\ge$ 2 be an integer. We present an algorithm for computing explicit rational representations of isogenies between Jacobians of hyperelliptic curves of genus g over an extension K of the field of p-adic numbers Qp. It relies on an efficient resolution, with a logarithmic loss of p-adic precision, of a first order system of differential equations.


title: Spectral Fundamentals and Characterizations of Signed Directed Graphs

abstract: The spectral properties of signed directed graphs, which may be naturally obtained by assigning a sign to each edge of a directed graph, have received substantially less attention than those of their undirected and/or unsigned counterparts. To represent such signed directed graphs, we use a striking equivalence to $\mathbb{T}_6$-gain graphs to formulate a Hermitian adjacency matrix, whose entries are the unit Eisenstein integers $\exp(k\pi i/3),$ $k\in \mathbb{Z}_6.$ Many well-known results, such as (gain) switching and eigenvalue interlacing, naturally carry over to this paradigm. We show that non-empty signed directed graphs whose spectra occur uniquely, up to isomorphism, do not exist, but we provide several infinite families whose spectra occur uniquely up to switching equivalence. Intermediate results include a classification of all signed digraphs with rank $2,3$, and a deep discussion of signed digraphs with extremely few (1 or 2) non-negative (eq. non-positive) eigenvalues.


title: Influence of lateral confinement on granular flows: comparison between shear-driven and gravity-driven flows

abstract: The properties of confined granular flows are studied through discrete numerical simulations. Two types of flows with different boundaries are compared: (i) gravity-driven flows topped with a free surface and over a base where erosion balances accretion (ii) shear-driven flows with a constant pressure applied at their top and a bumpy bottom moving at constant velocity. In both cases we observe shear localization over or/and under a creep zone. We show that, although the different boundaries induce different flow properties (e.g. shear localization of transverse velocity profiles), the two types of flow share common properties like (i) a power law relation between the granular temperature and the shear rate (whose exponent varies from 1 for dense flows to 2 for dilute flows) and (ii) a weakening of friction at the sidewalls which gradually decreases with the depth within the flow.


title: Broadband Frequency and Spatial On-Demand Tailoring of Topological Wave Propagation Harnessing Piezoelectric Metamaterials

abstract: Recent studies have developed tunable topological elastic metamaterials to maximize performance in the presence of varying external conditions, adapt to changing operating requirements, and enable new functionalities such as a programmable wave path. However, a challenge remains to achieve a tunable topological metamaterial that is comprehensively adaptable in both the frequency and spatial domains and is effective over a broad frequency bandwidth that includes a subwavelength regime. To advance the state of the art, this research presents a piezoelectric metamaterial with the capability to concurrently tailor the frequency, path, and mode shape of topological waves using resonant circuitry. In the research presented in this manuscript, the plane wave expansion method is used to detect a frequency tunable subwavelength Dirac point in the band structure of the periodic unit cell and discover an operating region over which topological wave propagation can exist. Dispersion analyses for a finite strip illuminate how circuit parameters can be utilized to adjust mode shapes corresponding to topological edge states. A further evaluation provides insight into how increased electromechanical coupling and lattice reconfiguration can be exploited to enhance the frequency range for topological wave propagation, increase achievable mode localization, and attain additional edge states. Topological guided wave propagation that is subwavelength in nature and adaptive in path, localization, and frequency is illustrated in numerical simulations of thin plate structures. Outcomes from the presented work indicate that the easily integrable and comprehensively tunable proposed metamaterial could be employed in applications requiring a multitude of functions over a broad frequency bandwidth.


title: Finding a Maximum Minimal Separator: Graph Classes and Fixed-Parameter Tractability

abstract: We study the problem of finding a maximum cardinality minimal separator of a graph. This problem is known to be NP-hard even for bipartite graphs. In this paper, we strengthen this hardness by showing that for planar bipartite graphs, the problem remains NP-hard. Moreover, for co-bipartite graphs and for line graphs, the problem also remains NP-hard. On the positive side, we give an algorithm deciding whether an input graph has a minimal separator of size at least $k$ that runs in time $2^{O(k)}n^{O(1)}$. We further show that a subexponential parameterized algorithm does not exist unless the Exponential Time Hypothesis (ETH) fails. Finally, we discuss a lower bound for polynomial kernelizations of this problem.


title: Double Oracle Algorithm for Computing Equilibria in Continuous Games

abstract: Many efficient algorithms have been designed to recover Nash equilibria of various classes of finite games. Special classes of continuous games with infinite strategy spaces, such as polynomial games, can be solved by semidefinite programming. In general, however, continuous games are not directly amenable to computational procedures. In this contribution, we develop an iterative strategy generation technique for finding a Nash equilibrium in a whole class of continuous two-person zero-sum games with compact strategy sets. The procedure, which is called the double oracle algorithm, has been successfully applied to large finite games in the past. We prove the convergence of the double oracle algorithm to a Nash equilibrium. Moreover, the algorithm is guaranteed to recover an approximate equilibrium in finitely-many steps. Our numerical experiments show that it outperforms fictitious play on several examples of games appearing in the literature. In particular, we provide a detailed analysis of experiments with a version of the continuous Colonel Blotto game.


title: Randomized Progressive Hedging methods for Multi-stage Stochastic Programming

abstract: Progressive Hedging is a popular decomposition algorithm for solving multi-stage stochastic optimization problems. A computational bottleneck of this algorithm is that all scenario subproblems have to be solved at each iteration. In this paper, we introduce randomized versions of the Progressive Hedging algorithm able to produce new iterates as soon as a single scenario subproblem is solved. Building on the relation between Progressive Hedging and monotone operators, we leverage recent results on randomized fixed point methods to derive and analyze the proposed methods. Finally, we release the corresponding code as an easy-to-use Julia toolbox and report computational experiments showing the practical interest of randomized algorithms, notably in a parallel context. Throughout the paper, we pay a special attention to presentation, stressing main ideas, avoiding extra-technicalities, in order to make the randomized methods accessible to a broad audience in the Operations Research community.


title: Gravitational wave lensing beyond general relativity: birefringence, echoes and shadows

abstract: Gravitational waves (GW), as light, are gravitationally lensed by intervening matter, deflecting their trajectories, delaying their arrival and occasionally producing multiple images. In theories beyond general relativity (GR), new gravitational degrees of freedom add an extra layer of complexity and richness to GW lensing. We develop a formalism to compute GW propagation beyond GR over general space-times, including kinetic interactions with new fields. Our framework relies on identifying the dynamical propagation eigenstates (linear combinations of the metric and additional fields) at leading order in a short-wave expansion. We determine these eigenstates and the conditions under which they acquire a different propagation speed around a lens. Differences in speed between eigenstates cause birefringence phenomena, including time delays between the metric polarizations (orthogonal superpositions of $h_+,h_\times$) observable without an electromagnetic counterpart. In particular, GW echoes are produced when the accumulated delay is larger than the signal's duration, while shorter time delays produce a scrambling of the wave-form. We also describe the formation of GW shadows as non-propagating metric components are sourced by the background of the additional fields around the lens. As an example, we apply our methodology to quartic Horndeski theories with Vainshtein screening and show that birefringence effects probe a region of the parameter space complementary to the constraints from the multi-messenger event GW170817. In the future, identified strongly lensed GWs and binary black holes merging near dense environments, such as active galactic nuclei, will fulfill the potential of these novel tests of gravity.


title: Brain Tumor Segmentation using 3D-CNNs with Uncertainty Estimation

abstract: Automation of brain tumors in 3D magnetic resonance images (MRIs) is key to assess the diagnostic and treatment of the disease. In recent years, convolutional neural networks (CNNs) have shown improved results in the task. However, high memory consumption is still a problem in 3D-CNNs. Moreover, most methods do not include uncertainty information, which is specially critical in medical diagnosis. This work proposes a 3D encoder-decoder architecture, based on V-Net \cite{vnet} which is trained with patching techniques to reduce memory consumption and decrease the effect of unbalanced data. We also introduce voxel-wise uncertainty, both epistemic and aleatoric using test-time dropout and data-augmentation respectively. Uncertainty maps can provide extra information to expert neurologists, useful for detecting when the model is not confident on the provided segmentation.


title: Fractional vertex-arboricity of planar graphs

abstract: We initiate a systematic study of the fractional vertex-arboricity of planar graphs and demonstrate connections to open problems concerning both fractional coloring and the size of the largest induced forest in planar graphs. In particular, the following three long-standing conjectures concern the size of a largest induced forest in a planar graph, and we conjecture that each of these can be generalized to the setting of fractional vertex-arboricity. In 1979, Albertson and Berman conjectured that every planar graph has an induced forest on at least half of its vertices, in 1987, Akiyama and Watanabe conjectured that every bipartite planar graph has an induced forest on at least five-eighths of its vertices, and in 2010, Kowalik, Lu\v{z}ar, and \v{S}krekovski conjectured that every planar graph of girth at least five has an induced forest on at least seven-tenths of its vertices. We make progress toward the fractional generalization of the latter of these, by proving that every planar graph of girth at least five has fractional vertex-arboricity at most $2 - 1/324$.


title: Sound, Complete, Linear-Space, Best-First Diagnosis Search

abstract: Various model-based diagnosis scenarios require the computation of the most preferred fault explanations. Existing algorithms that are sound (i.e., output only actual fault explanations) and complete (i.e., can return all explanations), however, require exponential space to achieve this task. As a remedy, to enable successful diagnosis on memory-restricted devices and for memory-intensive problem cases, we propose RBF-HS, a diagnostic search method based on Korf's well-known RBFS algorithm. RBF-HS can enumerate an arbitrary fixed number of fault explanations in best-first order within linear space bounds, without sacrificing the desirable soundness or completeness properties. Evaluations using real-world diagnosis cases show that RBF-HS, when used to compute minimum-cardinality fault explanations, in most cases saves substantial space (up to 98 %) while requiring only reasonably more or even less time than Reiter's HS-Tree, a commonly used and as generally applicable sound, complete and best-first diagnosis search.


title: A fresh look at the RR Lyrae population in the Draco dwarf spheroidal galaxy with Gaia

abstract: We present a catalogue of 285 RR Lyrae stars (RRLs) in the Draco dwarf spheroidal galaxy (dSph), obtained by combining data from a number of different surveys including the second data release (DR2) of the European Space Agency (ESA) cornerstone mission Gaia. We have determined individual distances to the RRLs in our sample using for the first time a Gaia G-band luminosity-metallicity relation (MG - [Fe/H]) and study the structure of the Draco dSph as traced by its RRL population. We find that the RRLs located in the western/south-western region of Draco appear to be closer to us, which may be a clue of interaction between Draco and the Milky Way (MW). The average distance modulus of Draco measured with the RRLs is 19.53 +/- 0.07 mag, corresponding to a distance of 80.5 +/- 2.6 kpc, in good agreement with previous determinations in the literature. Based on the pulsation properties of the RRLs we confirm the Oosterhoff-intermediate nature of Draco. We present an additional sample of 41 candidate RRLs in Draco, which we selected from the Gaia DR2 catalogue based on the uncertainty of their G-band magnitudes. Additional epoch data that will become available in the Gaia third data release (DR3) will help to confirm whether these candidates are bona-fide Draco RRLs.


title: Tuning Word2vec for Large Scale Recommendation Systems

abstract: Word2vec is a powerful machine learning tool that emerged from Natural Lan-guage Processing (NLP) and is now applied in multiple domains, including recom-mender systems, forecasting, and network analysis. As Word2vec is often used offthe shelf, we address the question of whether the default hyperparameters are suit-able for recommender systems. The answer is emphatically no. In this paper, wefirst elucidate the importance of hyperparameter optimization and show that un-constrained optimization yields an average 221% improvement in hit rate over thedefault parameters. However, unconstrained optimization leads to hyperparametersettings that are very expensive and not feasible for large scale recommendationtasks. To this end, we demonstrate 138% average improvement in hit rate with aruntime budget-constrained hyperparameter optimization. Furthermore, to makehyperparameter optimization applicable for large scale recommendation problemswhere the target dataset is too large to search over, we investigate generalizinghyperparameters settings from samples. We show that applying constrained hy-perparameter optimization using only a 10% sample of the data still yields a 91%average improvement in hit rate over the default parameters when applied to thefull datasets. Finally, we apply hyperparameters learned using our method of con-strained optimization on a sample to the Who To Follow recommendation serviceat Twitter and are able to increase follow rates by 15%.


title: Style-invariant Cardiac Image Segmentation with Test-time Augmentation

abstract: Deep models often suffer from severe performance drop due to the appearance shift in the real clinical setting. Most of the existing learning-based methods rely on images from multiple sites/vendors or even corresponding labels. However, collecting enough unknown data to robustly model segmentation cannot always hold since the complex appearance shift caused by imaging factors in daily application. In this paper, we propose a novel style-invariant method for cardiac image segmentation. Based on the zero-shot style transfer to remove appearance shift and test-time augmentation to explore diverse underlying anatomy, our proposed method is effective in combating the appearance shift. Our contribution is three-fold. First, inspired by the spirit of universal style transfer, we develop a zero-shot stylization for content images to generate stylized images that appearance similarity to the style images. Second, we build up a robust cardiac segmentation model based on the U-Net structure. Our framework mainly consists of two networks during testing: the ST network for removing appearance shift and the segmentation network. Third, we investigate test-time augmentation to explore transformed versions of the stylized image for prediction and the results are merged. Notably, our proposed framework is fully test-time adaptation. Experiment results demonstrate that our methods are promising and generic for generalizing deep segmentation models.


title: Exciton-polariton mediated interaction between two nitrogen-vacancy color centers in diamond using two-dimensional transition metal dichalcogenides

abstract: In this paper, starting from a quantum master equation, we discuss the interaction between two negatively charged Nitrogen-vacancy color centers in diamond via exciton-polaritons propagating in a two-dimensional transition metal dichalcogenide layer in close proximity to a diamond crystal. We focus on the optical 1.945 eV transition and model the Nitrogen-vacancy color centers as two-level (artificial) atoms. We find that the interaction parameters and the energy levels renormalization constants are extremely sensitive to the distance of the Nitrogen-vacancy centers to the transition metal dichalcogenide layer. Analytical expressions are obtained for the spectrum of the exciton-polaritons and for the damping constants entering the Lindblad equation. The conditions for occurrence of exciton mediated superradiance are discussed.


title: Senior Thesis for Haverford College Convex Optimization, Newton's Method and Interior Point Method

abstract: This paper consists of four general parts: convex sets; convex functions; convex optimization; and the interior-point algorithm. I will start by introducing the definition of convex sets and give three common convex set examples which will be used later in this paper, then prove the significant separating and supporting hyperplane theorems. Stepping into convex functions, in addition to offering definitions, I will also prove the first order and second-order conditions for convexity of a function, and then introduce a couple of examples that will be used in a convex optimization problem later. Next, I will provide the official definition of convex optimization problems and prove some characteristics they have, including the existence (through optimality criterion) and the uniqueness of a solution. I will also generate two convex optimization problems, one of which cannot be simply solved and requires additional skills. Afterward, I will introduce duality, for the sake of constructing the interior-point method. In the last section, I will first present the descent method and Newton's method, which serve as the foundation of the interior-point method. Then, I will show how to use the logarithmic barrier function and the central path to build up the interior-point method.


title: Isolation Distributional Kernel: A New Tool for Point & Group Anomaly Detection

abstract: We introduce Isolation Distributional Kernel as a new way to measure the similarity between two distributions. Existing approaches based on kernel mean embedding, which convert a point kernel to a distributional kernel, have two key issues: the point kernel employed has a feature map with intractable dimensionality; and it is {\em data independent}. This paper shows that Isolation Distributional Kernel (IDK), which is based on a {\em data dependent} point kernel, addresses both key issues. We demonstrate IDK's efficacy and efficiency as a new tool for kernel based anomaly detection for both point and group anomalies. Without explicit learning, using IDK alone outperforms existing kernel based point anomaly detector OCSVM and other kernel mean embedding methods that rely on Gaussian kernel. For group anomaly detection,we introduce an IDK based detector called IDK$^2$. It reformulates the problem of group anomaly detection in input space into the problem of point anomaly detection in Hilbert space, without the need for learning. IDK$^2$ runs orders of magnitude faster than group anomaly detector OCSMM.We reveal for the first time that an effective kernel based anomaly detector based on kernel mean embedding must employ a characteristic kernel which is data dependent.


title: End-to-End Prediction of Parcel Delivery Time with Deep Learning for Smart-City Applications

abstract: The acquisition of massive data on parcel delivery motivates postal operators to foster the development of predictive systems to improve customer service. Predicting delivery times successive to being shipped out of the final depot, referred to as last-mile prediction, deals with complicating factors such as traffic, drivers' behaviors, and weather. This work studies the use of deep learning for solving a real-world case of last-mile parcel delivery time prediction. We present our solution under the IoT paradigm and discuss its feasibility on a cloud-based architecture as a smart city application. We focus on a large-scale parcel dataset provided by Canada Post, covering the Greater Toronto Area (GTA). We utilize an origin-destination (OD) formulation, in which routes are not available, but only the start and end delivery points. We investigate three categories of convolutional-based neural networks and assess their performances on the task. We further demonstrate how our modeling outperforms several baselines, from classical machine learning models to referenced OD solutions. Specifically, we show that a ResNet architecture with 8 residual blocks displays the best trade-off between performance and complexity. We perform a thorough error analysis across the data and visualize the deep features learned to better understand the model behavior, making interesting remarks on data predictability. Our work provides an end-to-end neural pipeline that leverages parcel OD data as well as weather to accurately predict delivery durations. We believe that our system has the potential not only to improve user experience by better modeling their anticipation but also to aid last-mile postal logistics as a whole.


title: An Algebraic Approach to Fourier Transformation

abstract: The notion of Fourier transformation is described from an algebraic perspective that lends itself to applications in Symbolic Computation. We build the algebraic structures on the basis of a given Heisenberg group (in the general sense of nilquadratic groups enjoying a splitting property); this includes in particular the whole gamut of Pontryagin duality. The free objects in the corresponding categories are determined, and various examples are given. As a first step towards Symbolic Computation, we study two constructive examples in some detail -- the Gaussians (with and without polynomial factors) and the hyperbolic secant algebra.


title: Explaining Chemical Toxicity using Missing Features

abstract: Chemical toxicity prediction using machine learning is important in drug development to reduce repeated animal and human testing, thus saving cost and time. It is highly recommended that the predictions of computational toxicology models are mechanistically explainable. Current state of the art machine learning classifiers are based on deep neural networks, which tend to be complex and harder to interpret. In this paper, we apply a recently developed method named contrastive explanations method (CEM) to explain why a chemical or molecule is predicted to be toxic or not. In contrast to popular methods that provide explanations based on what features are present in the molecule, the CEM provides additional explanation on what features are missing from the molecule that is crucial for the prediction, known as the pertinent negative. The CEM does this by optimizing for the minimum perturbation to the model using a projected fast iterative shrinkage-thresholding algorithm (FISTA). We verified that the explanation from CEM matches known toxicophores and findings from other work.


title: Grain Surface Classification via Machine Learning Methods

abstract: In this study, radar signals were analyzed to classify grain surface types by using machine learning methods. Radar backscatter signals were recorded using a vector network analyzer between 18-40 GHz. A total of 5681 measurements of A scan signals were collected. The proposed method framework consists of two parts. First Order Statistical features are obtained by applying Fast Fourier Transform (FFT), Discrete Cosine Transform (DCT), Discrete Wavelet Transform (DWT) on backscatter signals in the first part of the framework. Classification process of these features was carried out with Support Vector Machine (SVM). In the second part of the proposed framework, two dimensional matrices in complex form were obtained by applying Short Time Fourier Transform (STFT) on the signals. Gray-Level Co-Occurrence Matrix (GLCM) and Gray-Level Run-Length Matrix (GLRLM) were obtained and feature extraction process was completed. Classification process was carried out with DVM. 10-k cross validation was applied. The highest performance was achieved with STFT+GLCM+SVM.


title: Integrating Battery Aging in the Optimization for Bidirectional Charging of Electric Vehicles

abstract: Smart charging of Electric Vehicles (EVs) reduces operating costs, allows more sustainable battery usage, and promotes the rise of electric mobility. In addition, bidirectional charging and improved connectivity enables efficient power grid support. Today, however, uncoordinated charging, e.g. governed by users' habits, is still the norm. Thus, the impact of upcoming smart charging applications is mostly unexplored. We aim to estimate the expenses inherent with smart charging, e.g. battery aging costs, and give suggestions for further research. Using typical on-board sensor data we concisely model and validate an EV battery. We then integrate the battery model into a realistic smart charging use case and compare it with measurements of real EV charging. The results show that i) the temperature dependence of battery aging requires precise thermal models for charging power greater than 7 kW, ii) disregarding battery aging underestimates EVs' operating costs by approx. 30%, and iii) the profitability of Vehicle-to-Grid (V2G) services based on bidirectional power flow, e.g. energy arbitrage, depends on battery aging costs and the electricity price spread.


title: How Much Does It Hurt: A Deep Learning Framework for Chronic Pain Score Assessment

abstract: Chronic pain is defined as pain that lasts or recurs for more than 3 to 6 months, often long after the injury or illness that initially caused the pain has healed. The "gold standard" for chronic pain assessment remains self report and clinical assessment via a biopsychosocial interview, since there has been no device that can measure it. A device to measure pain would be useful not only for clinical assessment, but potentially also as a biofeedback device leading to pain reduction. In this paper we propose an end-to-end deep learning framework for chronic pain score assessment. Our deep learning framework splits the long time-course data samples into shorter sequences, and uses Consensus Prediction to classify the results. We evaluate the performance of our framework on two chronic pain score datasets collected from two iterations of prototype Pain Meters that we have developed to help chronic pain subjects better understand their health condition.


title: A Low Rank Quaternion Decomposition Algorithm and Its Application in Color Image Inpainting

abstract: In this paper, we propose a lower rank quaternion decomposition algorithm and apply it to color image inpainting. We introduce a concise form for the gradient of a real function in quaternion matrix variables. The optimality conditions of our quaternion least squares problem have a simple expression with this form. The convergence and convergence rate of our algorithm are established with this tool.


title: Evasive Windows Malware: Impact on Antiviruses and Possible Countermeasures

abstract: The perpetual opposition between antiviruses and malware leads both parties to evolve continuously. On the one hand, antiviruses put in place solutions that are more and more sophisticated and propose more complex detection techniques in addition to the classic signature analysis. This sophistication leads antiviruses to leave more traces of their presence on the machine they protect. To remain undetected as long as possible, malware can avoid executing within such environments by hunting down the modifications left by the antiviruses. This paper aims at determining the possibilities for malware to detect the antiviruses and then evaluating the efficiency of these techniques on a panel of antiviruses that are the most used nowadays. We then collect samples showing this kind of behavior and propose to evaluate a countermeasure that creates false artifacts, thus forcing malware to evade.


title: A Note on Toroidal Maxwell-Cremona Correspondences

abstract: We explore toroidal analogues of the Maxwell-Cremona correspondence. Erickson and Lin [arXiv:2003.10057] showed the following correspondence for geodesic torus graphs $G$: a positive equilibrium stress for $G$, an orthogonal embedding of its dual graph $G^*$, and vertex weights such that $G$ is the intrinsic weighted Delaunay graph of its vertices. We extend their results to equilibrium stresses that are not necessarily positive, which correspond to orthogonal drawings of $G^*$ that are not necessarily embeddings. We also give a correspondence between equilibrium stresses and parallel drawings of the dual.


title: Mixed labyrinth fractals

abstract: Labyrinth fractals are self-similar fractals that were introduced and studied in recent work by Cristea and Steinsky. In the present paper we define and study more general objects, called mixed labyrinth fractals, that are in general not self-similar and are constructed by using sequences of labyrinth patterns. We show that mixed labyrinth fractals are dendrites and study properties of the paths in the graphs associated to prefractals, and of arcs in the fractal, e.g., the path length and the box counting dimension and length of arcs. We also consider more general objects related to mixed labyrinth fractals, formulate two conjectures about arc lengths, and establish connections to recent results on generalised Sierpinski carpets.


title: Internal Structure of the Jets from Young Stars Simulated at Plasma Focus Facilities

abstract: The laboratory simulations of jets from young stars that have been carried out for many years at plasma focus facilities allow the internal structure of the active regions emerging during the interaction of the jet with the surrounding plasma to be studied in detail. We have found a new wide class of solutions for the equations of ideal magnetohydrodynamics describing the closed axisymmetric stationary flows that are apparently realized in the active regions. Such flows are shown to well reproduce the internal structure of the plasma structures observed in laboratory simulations of astrophysical jets.


title: Search for integrable two-component versions of the lattice equations in the ABS-list

abstract: We search and classify two-component versions of the quad equations in the ABS list, under certain assumptions. The independent variables will be called $y,z$ and in addition to multilinearity and irreducibility the equation pair is required to have the following specific properties: (1) The two equations forming the pair are related by $y\leftrightarrow z$ exchange. (2) When $z=y$ both equations reduce to one of the equations in the ABS list. (3) Evolution in any corner direction is by a multilinear equation pair. One straightforward way to construct such two-component pairs is by taking some particular equation in the ABS list (in terms of $y$), using replacement $y \leftrightarrow z$ for some particular shifts, after which the other equation of the pair is obtained by property (1). This way we can get 8 pairs for each starting equation. One of our main results is that due to condition (3) this is in fact complete for H1, H3, Q1, Q3. (For H2 we have a further case, Q2, Q4 we did not check.) As for the CAC integrability test, for each choice of the bottom equations we could in principle have $8^2$ possible side-equations. However, we find that only equations constructed with an even number of $y \leftrightarrow z$ replacements are possible, and for each such equation there are two sets of "side" equation pairs that produce (the same) genuine B\"acklund transformation and Lax pair.


title: Restrained Italian domination in graphs

abstract: For a graph $G=(V(G),E(G))$, an Italian dominating function (ID function) $f:V(G)\rightarrow\{0,1,2\}$ has the property that for every vertex $v\in V(G)$ with $f(v)=0$, either $v$ is adjacent to a vertex assigned $2$ under $f$ or $v$ is adjacent to least two vertices assigned $1$ under $f$. The weight of an ID function is $\sum_{v\in V(G)}f(v)$. The Italian domination number is the minimum weight taken over all ID functions of $G$. In this paper, we initiate the study of a variant of ID functions. A restrained Italian dominating function (RID function) $f$ of $G$ is an ID function of $G$ for which the subgraph induced by $\{v\in V(G)\mid f(v)=0\}$ has no isolated vertices, and the restrained Italian domination number $\gamma_{rI}(G)$ is the minimum weight taken over all RID functions of $G$. We first prove that the problem of computing this parameter is NP-hard, even when restricted to bipartite graphs and chordal graphs as well as planar graphs with maximum degree five. We prove that $\gamma_{rI}(T)$ for a tree $T$ of order $n\geq3$ different from the double star $S_{2,2}$ can be bounded from below by $(n+3)/2$. Moreover, all extremal trees for this lower bound are characterized in this paper. We also give some sharp bounds on this parameter for general graphs and give the characterizations of graphs $G$ with small or large $\gamma_{rI}(G)$.


title: Emergence of complex data from simple local rules in a network game

abstract: As one of the main subjects of investigation in data science, network science has been demonstrated a wide range of applications to real-world networks analysis and modeling. For example, the pervasive presence of structural or topological characteristics, such as the small-world phenomenon, small-diameter, scale-free properties, or fat-tailed degree distribution were one of the underlying pillars fostering the study of complex networks. Relating these phenomena with other emergent properties in complex systems became a subject of central importance. By introducing new implications on the interface between data science and complex systems science with the purpose of tackling some of these issues, in this article we present a model for a network game played by complex networks in which nodes are computable systems. In particular, we present and discuss how some network topological properties and simple local communication rules are able to generate a phase transition with respect to the emergence of incompressible data.

categories: cs.LO cs.GT cs.MA cs.SI cs.SY eess.SY physics.soc-ph


title: Safe Coverage of Moving Domains for Vehicles with Second Order Dynamics

abstract: Autonomous coverage of a specified area by robots operating in close proximity with each other has many potential applications such as real-time monitoring of rapidly changing environments, and search and rescue; however, coordination and safety are two fundamental challenges. For coordination, we propose a distributed controller for covering moving, compact domains for two types of vehicles with second order dynamics (double integrator and fixed-wing aircraft) with bounded input forces. This control policy is based on artificial potentials and alignment forces designed to promote desired vehicle-domain and inter-vehicle separations and relative velocities. We prove that certain coverage configurations are locally asymptotically stable. For safety, we establish energy conditions for collision free motion and utilize Hamilton-Jacobi (HJ) reachability theory for last-resort pairwise collision avoidance. We derive an analytical solution to the associated HJ partial differential equation corresponding to the collision avoidance problem between two double integrator vehicles. We demonstrate our approach in several numerical simulations involving the two types of vehicles covering convex and non-convex moving domains.


title: Axion-Radiation Conversion by Super and Normal Conductors

abstract: We have proposed a method for the detection of dark matter axion. It uses superconductor under strong magnetic field. As is well known, the dark matter axion induces oscillating electric field under magnetic field. The electric field is proportional to the magnetic field and makes charged particles oscillate in conductors. Then, radiations of electromagnetic fields are produced. Radiation flux depends on how large the electric field is induced and how large the number of charged particles is present in the conductors. We show that the electric field in superconductor is essentially identical to the one induced in vacuum. It is proportional to the magnetic field. It is only present in the surface because of Meissner effect. On the other hand, although the magnetic field can penetrates the normal conductor, the oscillating electric field is only present in the surface of the conductor because of the skin effect. The strength of the electric field induced in the surface is equal to the one in vacuum. We obtain the electric field in the superconductor by solving equations of electromagnetic fields coupled with axion and Cooper pair described by Ginzburg-Landau model. The electric field in the normal conductor is obtained by solving equations of electromagnetic fields in the conductor coupled with axion. We compare radiation flux from the cylindrical superconductor with that from the normal conductor with same size. We find that the radiation flux from the superconductor is a hundred times larger than the flux from the normal conductor. We also show that when we use superconducting resonant cavity, we obtain radiation energy generated in the cavity two times of the order of the magnitude larger than that in normal conducting resonant cavity.


title: Towards a Systematic Computational Framework for Modeling Multi-Agent Decision-Making at Micro Level for Smart Vehicles in a Smart World

abstract: We propose a multi-agent based computational framework for modeling decision-making and strategic interaction at micro level for smart vehicles in a smart world. The concepts of Markov game and best response dynamics are heavily leveraged. Our aim is to make the framework conceptually sound and computationally practical for a range of realistic applications, including micro path planning for autonomous vehicles. To this end, we first convert the would-be stochastic game problem into a closely related deterministic one by introducing risk premium in the utility function for each individual agent. We show how the sub-game perfect Nash equilibrium of the simplified deterministic game can be solved by an algorithm based on best response dynamics. In order to better model human driving behaviors with bounded rationality, we seek to further simplify the solution concept by replacing the Nash equilibrium condition with a heuristic and adaptive optimization with finite look-ahead anticipation. In addition, the algorithm corresponding to the new solution concept drastically improves the computational efficiency. To demonstrate how our approach can be applied to realistic traffic settings, we conduct a simulation experiment: to derive merging and yielding behaviors on a double-lane highway with an unexpected barrier. Despite assumption differences involved in the two solution concepts, the derived numerical solutions show that the endogenized driving behaviors are very similar. We also briefly comment on how the proposed framework can be further extended in a number of directions in our forthcoming work, such as behavioral calibration using real traffic video data, computational mechanism design for traffic policy optimization, and so on.


title: Anomalous Nonlinear Dynamics Behavior of Fractional Viscoelastic Structures

abstract: Fractional models and their parameters are sensitive to changes in the intrinsic micro-structures of anomalous materials. We investigate how such physics-informed models propagate the evolving anomalous rheology to the nonlinear dynamics of mechanical systems. In particular, we analyze the vibration of a fractional, geometrically nonlinear viscoelastic cantilever beam, under base excitation and free vibration, where the viscoelastic response is general through a distributed-order fractional model. We employ Hamilton's principle to obtain the corresponding equation of motion with the choice of specific material distribution functions that recover a fractional Kelvin-Voigt viscoelastic model of order $\alpha$. Through spectral decomposition in space, the resulting time-fractional partial differential equation reduces to a nonlinear time-fractional ordinary differential equation, in which the linear counterpart is numerically integrated by employing a direct L1-difference scheme. We further develop a semi-analytical scheme to solve the nonlinear system through a method of multiple scales, which yields a cubic algebraic equation in terms of the frequency. Our numerical results suggest a set of $\alpha$-dependent anomalous dynamic qualities, such as far-from-equilibrium power-law amplitude decay rates, super-sensitivity of amplitude response at free vibration, and bifurcation in steady-state amplitude at primary resonance.


title: Matrix-Monotonic Optimization Part II: Multi-Variable Optimization

abstract: In contrast to Part I of this treatise [1] that focuses on the optimization problems associated with single matrix variables, in this paper, we investigate the application of the matrix-monotonic optimization framework in the optimization problems associated with multiple matrix variables. It is revealed that matrix-monotonic optimization still works even for multiple matrix-variate based optimization problems, provided that certain conditions are satisfied. Using this framework, the optimal structures of the matrix variables can be derived and the associated multiple matrix-variate optimization problems can be substantially simplified. In this paper, several specific examples are given, which are essentially open problems. Firstly, we investigate multi-user multiple-input multiple-output (MU- MIMO) uplink communications under various power constraints. Using the proposed framework, the optimal structures of the precoding matrices at each user under various power constraints can be derived. Secondly, we considered the optimization of the signal compression matrices at each sensor under various power constraints in distributed sensor networks. Finally, we investigate the transceiver optimization for multi-hop amplify-and-forward (AF) MIMO relaying networks with imperfect channel state information (CSI) under various power constraints. At the end of this paper, several simulation results are given to demonstrate the accuracy of the proposed theoretical results.


title: Deep Learning of Individual Aesthetics

abstract: Accurate evaluation of human aesthetic preferences represents a major challenge for creative evolutionary and generative systems research. Prior work has tended to focus on feature measures of the artefact, such as symmetry, complexity and coherence. However, research models from Psychology suggest that human aesthetic experiences encapsulate factors beyond the artefact, making accurate computational models very difficult to design. The interactive genetic algorithm (IGA) circumvents the problem through human-in-the-loop, subjective evaluation of aesthetics, but is limited due to user fatigue and small population sizes. In this paper we look at how recent advances in deep learning can assist in automating personal aesthetic judgement. Using a leading artist's computer art dataset, we investigate the relationship between image measures, such as complexity, and human aesthetic evaluation. We use dimension reduction methods to visualise both genotype and phenotype space in order to support the exploration of new territory in a generative system. Convolutional Neural Networks trained on the artist's prior aesthetic evaluations are used to suggest new possibilities similar or between known high quality genotype-phenotype mappings. We integrate this classification and discovery system into a software tool for evolving complex generative art and design.


title: Latent Causal Socioeconomic Health Index

abstract: This research develops a model-based LAtent Causal Socioeconomic Health (LACSH) index at the national level. Motivated by the need for a holistic national well-being index, we build upon the latent health factor index (LHFI) approach that has been used to assess the unobservable ecological/ecosystem health. LHFI integratively models the relationship between metrics, latent health, and covariates that drive the notion of health. In this paper, the LHFI structure is integrated with spatial modeling and statistical causal modeling. Our efforts are focused on developing the integrated framework to facilitate the understanding of how an observational continuous variable might have causally affected a latent trait that exhibits spatial correlation. A novel visualization technique to evaluate covariate balance is also introduced for the case of a continuous policy (treatment) variable. Our resulting LACSH framework and visualization tool are illustrated through two global case studies on national socioeconomic health (latent trait), each with various metrics and covariates pertaining to different aspects of societal health, and the treatment variable being mandatory maternity leave days and government expenditure on healthcare, respectively. We validate our model by two simulation studies. All approaches are structured in a Bayesian hierarchical framework and results are obtained by Markov chain Monte Carlo techniques.


title: Gribov-Zwanziger confinement, high energy evolution and large impact parameter behaviour of the scattering amplitude

abstract: In this paper we derive the high energy evolution equation in the Gribov-Zwanziger approach, for the confinement of quarks and gluons. We demonstrate that the new equation generates an exponential decrease of the scattering amplitude at large impact parameter, and resolves the main difficulties of CGC (Colour Glass Condensate) high energy effective theory. Such behaviour occurs if the gluon propagator in Gribov-Zwanziger approach, does not vanish at small momenta. Solving the non-linear equation for deep inelastic scattering, we show that the suggested equation leads to a Froissart disc with radius ($R_F$), which increases as $ R_F \,\propto Y = \ln (1/x)$, and with a finite width for the distribution over $| b - R_F|$.