title: Inverse Structural Design of Graphene/Boron Nitride Hybrids by Regressional GAN

abstract: Inverse design of materials with desired properties is currently laborious and heavily relies on intuition of researchers through a trial-and-error process. The massive combinational spaces due to the constituent elements and their structural configurations are too overwhelming to be all searched even by high-throughput computations. Herein, we demonstrated a novel regressional generative adversarial network (RGAN) for inverse design of representative two-dimensional materials, graphene and boron-nitride (BN) hybrids. A significant novelty of the proposed RGAN is that it combines the supervised and regressional convolutional neural network (CNN) with the traditional unsupervised GAN, thus overcoming the common technical barrier in the traditional GANs, which cannot generate data associated with given continuous quantitative labels. The proposed RGAN enables to autonomously generate graphene/BN hybrids with any given bandgaps. Moreover, the generated structures exhibit high fidelity, yielding bandgaps within ~ 10% MAEF of the desired bandgaps as cross-validated by density functional theory (DFT) calculations. Further analysis by principle component analysis (PCA) and modified locally linear embedding (MLLE) methods on the latent features encoded by the regressor reveals that the generator has successfully generated structures that followed the statistical distribution of the real structures. It implies the possibility of the RGAN in recognizing physical rules hidden in the high-dimensional data. This new inverse design methodology would speed up the discovery and development of other 2D materials and beyond.


title: The effects of ecological selection on species diversity and trait distribution: predictions and an empirical test

abstract: Ecological selection is a major driver of community assembly. Selection is classified as stabilizing when species with intermediate trait values gain the highest reproductive success, whereas selection is considered directional when fitness is highest for species with extreme trait values. Previous studies have investigated the effects of different selection types on trait distribution, but the effects of selection on species diversity have remained unclear. Here, we propose a framework for inferring the type and strength of selection by studying species diversity and trait distribution together against null expectations. We use a simulation model to confirm our prediction that directional selection should lead to lower species diversity than stabilizing selection despite a similar effect on trait community-weighted variance. We apply the framework to a mesocosm system of annual plants to test whether differences in species diversity between two habitats that vary in productivity are related to differences in selection on seed mass. We show that, in both habitats, species diversity was lower than the null expectation, but that species diversity was lower in the more productive habitat. We attribute this difference to strong directional selection for large-seeded species in the productive habitat as indicated by trait community-weighted-mean being higher and community-weighted variance being lower than the null expectations. In the less productive habitat, we found that community-weighted variance was higher than expected by chance, suggesting that seed mass could be a driver of niche partitioning under such conditions. Altogether, our results suggest that viewing species diversity and trait distribution as interrelated patterns driven by the same process, ecological selection, is helpful in understanding community assembly.


title: An improved muon track reconstruction for IceCube

abstract: IceCube is a cubic-kilometer Cherenkov telescope operating at the South Pole. One of its main objectives is to detect astrophysical neutrinos and identify their sources. High-energy muon neutrinos are identified through the secondary muons produced via charge current interactions with the ice. The present best-performing directional reconstruction of the muon track is a maximum likelihood method which uses the arrival time distribution of Cherenkov photons registered by the experiment's photomultipliers. A known systematic shortcoming of the prevailing method is to assume a continuous energy loss along the muon track. This contribution discusses a generalized Ansatz where the expected arrival time distribution is parametrized by a stochastic muon energy loss pattern. This more realistic parametrization of the muon energy loss profile leads to an improvement of about 20% to the muon angular resolution of IceCube.


title: Estimation of perceptual scales using ordinal embedding

abstract: In this paper, we address the problem of measuring and analysing sensation, the subjective magnitude of one's experience. We do this in the context of the method of triads: the sensation of the stimulus is evaluated via relative judgments of the form: "Is stimulus S_i more similar to stimulus S_j or to stimulus S_k?". We propose to use ordinal embedding methods from machine learning to estimate the scaling function from the relative judgments. We review two relevant and well-known methods in psychophysics which are partially applicable in our setting: non-metric multi-dimensional scaling (NMDS) and the method of maximum likelihood difference scaling (MLDS). We perform an extensive set of simulations, considering various scaling functions, to demonstrate the performance of the ordinal embedding methods. We show that in contrast to existing approaches our ordinal embedding approach allows, first, to obtain reasonable scaling function from comparatively few relative judgments, second, the estimation of non-monotonous scaling functions, and, third, multi-dimensional perceptual scales. In addition to the simulations, we analyse data from two real psychophysics experiments using ordinal embedding methods. Our results show that in the one-dimensional, monotonically increasing perceptual scale our ordinal embedding approach works as well as MLDS, while in higher dimensions, only our ordinal embedding methods can produce a desirable scaling function. To make our methods widely accessible, we provide an R-implementation and general rules of thumb on how to use ordinal embedding in the context of psychophysics.


title: Clustering Longitudinal Life-Course Sequences Using Mixtures of Exponential-Distance Models

abstract: Sequence analysis is an increasingly popular approach for analysing life courses represented by ordered collections of activities experienced by subjects over time. Here, we analyse a survey data set containing information on the career trajectories of a cohort of Northern Irish youths tracked between the ages of 16 and 22. We propose a novel, model-based clustering approach suited to the analysis of such data from a holistic perspective, with the aims of estimating the number of typical career trajectories, identifying the relevant features of these patterns, and assessing the extent to which such patterns are shaped by background characteristics. Several criteria exist for measuring pairwise dissimilarities among categorical sequences. Typically, dissimilarity matrices are employed as input to heuristic clustering algorithms. The family of methods we develop instead clusters sequences directly using mixtures of exponential-distance models. Basing the models on weighted variants of the Hamming distance metric permits closed-form expressions for parameter estimation. Simultaneously allowing the component membership probabilities to depend on fixed covariates and accommodating sampling weights in the clustering process yields new insights on the Northern Irish data. In particular, we find that school examination performance is the single most important predictor of cluster membership.


title: Constrained Thompson Sampling for Real-Time Electricity Pricing with Grid Reliability Constraints

abstract: We consider the problem of an aggregator attempting to learn customers' load flexibility models while implementing a load shaping program by means of broadcasting daily dispatch signals. We adopt a multi-armed bandit formulation to account for the stochastic and unknown nature of customers' responses to dispatch signals. We propose a constrained Thompson sampling heuristic, Con-TS-RTP, that accounts for various possible aggregator objectives (e.g., to reduce demand at peak hours, integrate more intermittent renewable generation, track a desired daily load profile, etc) and takes into account the operational constraints of a distribution system to avoid potential grid failures as a result of uncertainty in the customers' response. We provide a discussion on the regret bounds for our algorithm as well as a discussion on the operational reliability of the distribution system's constraints being upheld throughout the learning process.


title: Case Study: Disclosure of Indirect Device Fingerprinting in Privacy Policies

abstract: Recent developments in online tracking make it harder for individuals to detect and block trackers. Some sites have deployed indirect tracking methods, which attempt to uniquely identify a device by asking the browser to perform a seemingly-unrelated task. One type of indirect tracking, Canvas fingerprinting, causes the browser to render a graphic recording rendering statistics as a unique identifier. In this work, we observe how indirect device fingerprinting methods are disclosed in privacy policies, and consider whether the disclosures are sufficient to enable website visitors to block the tracking methods. We compare these disclosures to the disclosure of direct fingerprinting methods on the same websites. Our case study analyzes one indirect fingerprinting technique, Canvas fingerprinting. We use an existing automated detector of this fingerprinting technique to conservatively detect its use on Alexa Top 500 websites that cater to United States consumers, and we examine the privacy policies of the resulting 28 websites. Disclosures of indirect fingerprinting vary in specificity. None described the specific methods with enough granularity to know the website used Canvas fingerprinting. Conversely, many sites did provide enough detail about usage of direct fingerprinting methods to allow a website visitor to reliably detect and block those techniques. We conclude that indirect fingerprinting methods are often difficult to detect and are not identified with specificity in privacy policies. This makes indirect fingerprinting more difficult to block, and therefore risks disturbing the tentative armistice between individuals and websites currently in place for direct fingerprinting. This paper illustrates differences in fingerprinting approaches, and explains why technologists, technology lawyers, and policymakers need to appreciate the challenges of indirect fingerprinting.


title: Enabling and Exploiting Partition-Level Parallelism (PALP) in Phase Change Memories

abstract: Phase-change memory (PCM) devices have multiple banks to serve memory requests in parallel. Unfortunately, if two requests go to the same bank, they have to be served one after another, leading to lower system performance. We observe that a modern PCM bank is implemented as a collection of partitions that operate mostly independently while sharing a few global peripheral structures, which include the sense amplifiers (to read) and the write drivers (to write). Based on this observation, we propose PALP, a new mechanism that enables partition-level parallelism within each PCM bank, and exploits such parallelism by using the memory controller's access scheduling decisions. PALP consists of three new contributions. First, we introduce new PCM commands to enable parallelism in a bank's partitions in order to resolve the read-write bank conflicts, with minimal changes needed to PCM logic and its interface. Second, we propose simple circuit modifications that introduce a new operating mode for the write drivers, in addition to their default mode of serving write requests. When configured in this new mode, the write drivers can resolve the read-read bank conflicts, working jointly with the sense amplifiers. Finally, we propose a new access scheduling mechanism in PCM that improves performance by prioritizing those requests that exploit partition-level parallelism over other requests, including the long outstanding ones. While doing so, the memory controller also guarantees starvation-freedom and the PCM's running-average-power-limit (RAPL). We evaluate PALP with workloads from the MiBench and SPEC CPU2017 Benchmark suites. Our results show that PALP reduces average PCM access latency by 23%, and improves average system performance by 28% compared to the state-of-the-art approaches.


title: Multi-Antenna Relaying and Reconfigurable Intelligent Surfaces: End-to-End SNR and Achievable Rate

abstract: In this report, we summarize the end-to-end signal-to-noise ratio and the rate of half-duplex, full-duplex, amplify-and-forward, and decode-and-forward relay-aided communications, and well as the signal-to-noise ratio and the rate of the emerging technology known as reconfigurable intelligent surfaces.


title: Assessing the Impact of a User-Item Collaborative Attack on Class of Users

abstract: Collaborative Filtering (CF) models lie at the core of most recommendation systems due to their state-of-the-art accuracy. They are commonly adopted in e-commerce and online services for their impact on sales volume and/or diversity, and their impact on companies' outcome. However, CF models are only as good as the interaction data they work with. As these models rely on outside sources of information, counterfeit data such as user ratings or reviews can be injected by attackers to manipulate the underlying data and alter the impact of resulting recommendations, thus implementing a so-called shilling attack. While previous works have focused on evaluating shilling attack strategies from a global perspective paying particular attention to the effect of the size of attacks and attacker's knowledge, in this work we explore the effectiveness of shilling attacks under novel aspects. First, we investigate the effect of attack strategies crafted on a target user in order to push the recommendation of a low-ranking item to a higher position, referred to as user-item attack. Second, we evaluate the effectiveness of attacks in altering the impact of different CF models by contemplating the class of the target user, from the perspective of the richness of her profile (i.e., cold v.s. warm user). Finally, similar to previous work we contemplate the size of attack (i.e., the amount of fake profiles injected) in examining their success. The results of experiments on two widely used datasets in business and movie domains, namely Yelp and MovieLens, suggest that warm and cold users exhibit contrasting behaviors in datasets with different characteristics.


title: Magnetic signatures of domain walls in $s+is$ and $s+id$ superconductors: observability and what that can tell us about the superconducting order parameter

abstract: One of the defining features of spontaneously broken time-reversal symmetry (BTRS) is the existence of domain walls, the detection of which would be strong evidence for such systems. There is keen interest in BTRS currently, in part, due to recent muon spin rotation experiments, which have pointed towards $\textrm{Ba}_{1-x}\textrm{K}_x\textrm{Fe}_2\textrm{As}_2$ exhibiting a remarkable case of $s$-wave superconductivity with spontaneously broken time-reversal symmetry. A key question, however, is how to differentiate between the different theoretical models which describe such a state. Two particularly popular choices of model are $s+is$ and $s+id$ superconducting states. In this paper, we obtain solutions for domain walls in $s+is$ and $s+id$ systems, including the effects of lattice anisotropies. We show that, in general, both models exhibit spontaneous magnetic field, that extend along the entire length of the domain wall. We demonstrate the qualitative difference between the magnetic signatures of $s+is$ and $s+id$ domain walls and propose a procedure to extract the superconducting pairing symmetry from the magnetic-field response of domain walls.


title: Static and Dynamic Charged Black Holes

abstract: We consider a class of Einstein-Maxwell-dilaton theories in general dimensions and construct both static and dynamic charged black holes. We adopt the reverse engineering procedure and make a specific ansatz for the scalar field and then derive the necessary scalar potential and the non-minimal coupling function between the scalar and the Maxwell field. The resulting static black holes contain mass and electric charge as integration constants. We find that some of the static solutions can be promoted to become dynamical ones in the Eddington-Finkelstein-like coordinates. The collapse solutions describe the evolution from a smaller charged black hole to a larger black hole state, driven by the scalar field.


title: Equilibrium and stability of two-dimensional pinned drops

abstract: Superhydrophobicity relies on the stability of drops's interfaces pinned on sharp edges to sustain non-wetting (Cassie-Baxter) equilibrium states. Gibbs already pointed out that equilibrium is possible as long as the pinning angle at the edge falls between the equilibrium contact angles corresponding to the flanks of the edge. However, the lack of stability can restrict further the realizable equilibrium configurations. To find these limits we analyze here the equilibrium and stability of two-dimensional drops bounded by interfaces pinned on mathematically sharp edges. We are specifically interested on how the drop's stability depends on its size, which is measured with the Bond number $Bo = (\mathcal{W}_d/\ell_c)^2$, defined as the ratio of the drop's characteristic length scale $\mathcal{W}_d$ to the capillary length $\ell_c = \sqrt{\sigma/\rho g}$. Drops with a fixed volume become more stable as they shrink in size. On the contrary, open drops, i.e. capable of exchanging mass with a reservoir, are less stable as their associated Bond number decreases.


title: Main trends of the quasar main sequence -- effect of viewing angle

abstract: We address the effect of the viewing angle of the accretion disk plane and the geometry of the broad-line region (BLR) with the goal of interpreting the distribution of quasars along the main sequence (MS). We utilize photoionization code CLOUDY to model the BLR FeII emission, incorporating the grossly underestimated role of the form factor (f). We recover the dependence of the strength of the FeII emission in the optical (R$_{\rm{FeII}}$) on L$_{\rm{bol}}$/L$_{\rm{Edd}}$ ratio and related observational trends - as a function of the spectral energy distribution (SED) shape, cloud density, composition and intra-cloud dynamics, assumed following prior observational constraints. With this approach, we are now able to explain the diversity of quasars and the change of the quasar properties along the Main Sequence (MS). Our approach also explains the rarity of the highest FeII emitters known as the extreme xA sources and can be used as a predictive tool in future reverberation mapping studies of Type-1 AGNs. This approach further justifies the use of quasars as `cosmological probes'.


title: Rotational Doppler cooling and heating

abstract: Doppler cooling is a widely used technique to laser cool atoms and nanoparticles exploiting the Doppler shift involved in translational transformations. The rotational Doppler effect arising from rotational coordinate transformations should similarly enable optical manipulations of the rotational degrees of freedom in rotating nanosystems. Here, we show that rotational Doppler cooling and heating (RDC and RDH) effects embody rich and unexplored physics, such as a strong dependence on particle morphology. For geometrically confined particles, such as a nanorod that can represent diatomic molecules, RDC and RDH follow similar rules as their translational Doppler counterpart, where cooling and heating are always observed at red- or blue-detuned laser frequencies, respectively. Surprisingly, nanosystems that can be modeled as a solid particle shows a strikingly different response, where RDH appears in a frequency regime close to their resonances, while a detuned frequency produces cooling of rotation. We also predict that the RDH effect can lead to unprecedented spontaneous chiral symmetry breaking, whereby an achiral particle under linearly polarized illumination starts spontaneously rotating, rendering it nontrivial compared to the translational Doppler effect. Our results open up new exciting possibilities to control the rotational motion of molecules and nanoparticles.


title: Deterministic Epidemic Models For Ebola Infection With Time-dependent Controls

abstract: In this paper, we have studied epidemiological models for Ebola infection using nonlinear ordinary differential equations and optimal control theory. We considered optimal control analysis of SIR and SEIR models for the deadly Ebola infection using vaccination, treatment and educational campaign as time-dependent controls functions. We have applied indirect methods to study existing deterministic optimal control epidemic models for Ebola virus disease. These methods in optimal control are based on Hamiltonian function and the Pontryagin's maximum principle to construct adjoint equations and optimality systems. The forward-backward sweep numerical scheme with fourth-order Runge-Kutta method is used to solve the optimality system for the various control strategies. From our numerical illustrations, we can conclude that, effective educational campaigns and vaccination of susceptible individuals as were as effective treatments of infected individuals can help reduce the disease transmission.


title: Thermodynamical extension of a symplectic numerical scheme with half space and time shift demonstrated on rheological waves in solids

abstract: On the example of the Poynting-Thomson-Zener rheological model for solids, which exhibits both dissipation and wave propagation - with nonlinear dispersion relation -, we introduce and investigate a finite difference numerical scheme. Our goal is to demonstrate its properties and to ease the computations in later applications for continuum thermodynamical problems. The key element is the positioning of the discretized quantities with shifts by half space and time steps with respect to each other. The arrangement is chosen according to the spacetime properties of the quantities and of the equations governing them. Numerical stability, dissipative error and dispersive error are analysed in detail. With the best settings found, the scheme is capable of making precise and fast predictions. Finally, the proposed scheme is compared to a commercial finite element software, COMSOL, which demonstrates essential differences even on the simplest - elastic - level of modelling.


title: A Multi-level Clustering Approach for Anonymizing Large-Scale Physical Activity Data

abstract: Publishing physical activity data can facilitate reproducible health-care research in several areas such as population health management, behavioral health research, and management of chronic health problems. However, publishing such data also brings high privacy risks related to re-identification which makes anonymization necessary. One of the challenges in anonymizing physical activity data collected periodically is its sequential nature. The existing anonymization techniques work sufficiently for cross-sectional data but have high computational costs when applied directly to sequential data. This paper presents an effective anonymization approach, Multi-level Clustering based anonymization to anonymize physical activity data. Compared with the conventional methods, the proposed approach improves time complexity by reducing the clustering time drastically. While doing so, it preserves the utility as much as the conventional approaches.


title: Bloch wave approach to almost periodic homogenization and approximations of effective coefficients

abstract: Bloch wave homogenization is a spectral method for obtaining effective coefficients for periodically heterogeneous media. This method hinges on the direct integral decomposition of periodic operators, which is not available in a suitable form for almost periodic operators. In particular, the notion of Bloch eigenvalues and eigenvectors does not exist for almost periodic operators. However, we are able to recover the homogenization result in this case, by employing a sequence of periodic approximations to almost periodic operators. We also establish a rate of convergence for approximations of homogenized tensors for a class of almost periodic media. The results are supported by a numerical study.


title: Quantile Convolutional Neural Networks for Value at Risk Forecasting

abstract: This article presents a new method for forecasting Value at Risk. Convolutional neural networks can do time series forecasting, since they can learn local patterns in time. A simple modification enables them to forecast not the mean, but arbitrary quantiles of the distribution, and thus allows them to be applied to VaR-forecasting. The proposed model can learn from the price history of different assets, and it seems to produce fairly accurate forecasts.


title: Efficient and powerful equivalency test on combined mean and variance with application to diagnostic device comparison studies

abstract: In medical device comparison studies, equivalency test is commonly used to demonstrate two measurement methods agree up to a pre-specified performance goal based on the paired repeated measures. Such equivalency test often involves controlling the absolute differences that depend on both the mean and variance parameters, and poses some challenges for statistical analysis. For example, for the oximetry comparison study that motivates our research, FDA has clear guidelines approving an investigational pulse oximeter in comparison to a standard oximeter via testing the root mean squares (RMS), a composite measure of both mean and variance parameters. For the hypothesis testing of this composite measure, existing methods have been either exploratory or relying on the large-sample normal approximation with conservative and unsatisfactory performance. We develop a novel generalized pivotal test to rigorously and accurately test the system equivalency based on RMS. The proposed method has well-controlled type I error and favorable performance in our extensive numerical studies. When analyzing data from an oximetry comparison study, aiming to demonstrate performance equivalency between an FDA-cleared oximetry system and an investigational system, our proposed method resulted in a highly significant test result strongly supporting the system equivalency. We also provide efficient R programs for the proposed method in a publicly available R package. Considering that many practical equivalency studies of diagnostic devices are of small to medium sizes, our proposed method and software timely bridge an existing gap in the field.


title: A tree-based radial basis function method for noisy parallel surrogate optimization

abstract: Parallel surrogate optimization algorithms have proven to be efficient methods for solving expensive noisy optimization problems. In this work we develop a new parallel surrogate optimization algorithm (ProSRS), using a novel tree-based "zoom strategy" to improve the efficiency of the algorithm. We prove that if ProSRS is run for sufficiently long, with probability converging to one there will be at least one point among all the evaluations that will be arbitrarily close to the global minimum. We compare our algorithm to several state-of-the-art Bayesian optimization algorithms on a suite of standard benchmark functions and two real machine learning hyperparameter-tuning problems. We find that our algorithm not only achieves significantly faster optimization convergence, but is also 1-4 orders of magnitude cheaper in computational cost.


title: Counter-Rotation and High Velocity Outflow in the Parsec-Scale Molecular Torus of NGC 1068

abstract: We present 1.4 pc resolution observations of 256 GHz nuclear radio continuum and HCN ($J=3 \to 2$) in the molecular torus of NGC 1068. The integrated radio continuum emission has a flat spectrum consistent with free-free emission and resolves into an X-shaped structure resembling an edge-brightened bicone. HCN is detected in absorption against the continuum, and the absorption spectrum shows a pronounced blue wing that suggests a high-velocity molecular outflow with speeds reaching 450 km/s. Analysis of the off-nucleus emission line kinematics and morphology reveals two nested, rotating disk components. The inner disk, inside $r\sim 1.2$ pc, has kinematics consistent with the nearly edge-on, geometrically thin water megamaser disk in Keplerian rotation around a central mass of $1.66\times 10^7\,\mbox{M}_\odot$. The outer disk, which extends to $\sim 7$~pc radius, counter-rotates relative to the inner disk. The rotation curve of the outer disk is consistent with rotation around the same central mass as the megamaser disk but in the opposite sense. The morphology of the molecular gas is asymmetric around the nuclear continuum source. We speculate that the outer disk formed from more recently introduced molecular gas falling out of the host galaxy or from a captured dwarf satellite galaxy. In NGC 1068, we find direct evidence that the molecular torus consists of counter-rotating and misaligned disks on parsec scales.


title: Nonaffine lattice dynamics with the Ewald method reveals strongly nonaffine elasticity of {\alpha}-quartz

abstract: A lattice dynamical formalism based on nonaffine response theory is derived for non-centrosymmetric crystals, accounting for long-range interatomic interactions using the Ewald method. The framework takes equilibrated static configurations as input to compute the elastic constants in excellent agreement with both experimental data and calculations under strain. Besides this methodological improvement, which enables faster evaluation of elastic constants without the need of explicitly simulating the deformation process, the framework provides insights into the nonaffine contribution to the elastic constants of {\alpha}-quartz. It turns out that, due to the non-centrosymmetric lattice structure, the nonaffine (softening) correction to the elastic constants is very large, such that the overall elastic constants are at least 3-4 times smaller than the affine Born-Huang estimat


title: A way forward towards improvement of tensor force in \textit{pf}-shell

abstract: In many shell model interactions, the tensor force monopole matrix elements often retain systematic trends originating in the bare tensor force. In this work, however, we note for GX-interactions of \textit{pf}-shell that the seven out of ten T = 1 tensor force monopole matrix elements do not share these systematic. We ameliorate this disparity making use of Yukawa-type tensor force and spin-tensor decomposition. Furthermore, we modify the single-particle energy of $1p_{3/2}$ orbit and two TBMEs of $0f$-orbit,and test the revised interaction from Ca to Ge isotopes with various physics viewpoints. The results are found to be satisfactory with respect to the experimental data.


title: Minimal residual multistep methods for large stiff non-autonomous linear problems

abstract: The purpose of this work is to introduce a new idea of how to avoid the factorization of large matrices during the solution of stiff systems of ODEs. Starting from the general form of an explicit linear multistep method we suggest to adaptively choose its coefficients on each integration step in order to minimize the norm of the residual of an implicit BDF formula. Thereby we reduce the number of unknowns on each step from $n$ to $O(1)$, where $n$ is the dimension of the ODE system. We call this type of methods Minimal Residual Multistep (MRMS) methods. In the case of linear non-autonomous problem, besides the evaluations of the right-hand side of ODE, the resulting numerical scheme additionally requires one solution of a linear least-squares problem with a thin matrix per step. We show that the order of the method and its zero-stability properties coincide with those of the used underlying BDF formula. For the simplest analog of the implicit Euler method the properties of linear stability are investigated. Though the classical absolute stability analysis is not fully relevant to the MRMS methods, it is shown that this one-step method is applicable in stiff case. In the numerical experiment section we consider the fixed-step integration of a two-dimensional non-autonomous heat equation using the MRMS methods and their classical BDF counterparts. The starting values are taken from a preset slowly-varying exact solution. The comparison showed that both methods give similar numerical solutions, but in the case of large systems the MRMS methods are faster, and their advantage considerably increases with the growth of dimension. Python code with the experimantal code can be downloaded from the GitHub repository https://github.com/bfaleichik/mrms.


title: MobiSR: Efficient On-Device Super-Resolution through Heterogeneous Mobile Processors

abstract: In recent years, convolutional networks have demonstrated unprecedented performance in the image restoration task of super-resolution (SR). SR entails the upscaling of a single low-resolution image in order to meet application-specific image quality demands and plays a key role in mobile devices. To comply with privacy regulations and reduce the overhead of cloud computing, executing SR models locally on-device constitutes a key alternative approach. Nevertheless, the excessive compute and memory requirements of SR workloads pose a challenge in mapping SR networks on resource-constrained mobile platforms. This work presents MobiSR, a novel framework for performing efficient super-resolution on-device. Given a target mobile platform, the proposed framework considers popular model compression techniques and traverses the design space to reach the highest performing trade-off between image quality and processing speed. At run time, a novel scheduler dispatches incoming image patches to the appropriate model-engine pair based on the patch's estimated upscaling difficulty in order to meet the required image quality with minimum processing latency. Quantitative evaluation shows that the proposed framework yields on-device SR designs that achieve an average speedup of 2.13x over highly-optimized parallel difficulty-unaware mappings and 4.79x over highly-optimized single compute engine implementations.


title: Glitching pulsars: unraveling the interactions of general relativity with quantum fields in the strong field regimes

abstract: We present a modification of our previous model for the mechanisms underlying the glitch phenomena in pulsars and young neutron stars. Accordingly, pulsars are born with embryonic cores comprising of purely incompressible superconducting gluon-quark superfluid (henceforth SuSu-cores). As the ambient medium cools and spins down due to emission of magnetic dipole radiation, the mass and size of SuSu-cores are set to grow discretely with time, in accordance with the Onsager-Feynmann analysis of superfluidity. Presently, we propose that the spacetime embedding glitching pulsars is dynamical and of bimetric nature: inside SuSu-cores the spacetime must be flat, whereas the surrounding region, where the matter is compressible and dissipative, the spacetime is Schwarzschild. It is further proposed that the topological change of spacetime is derived by the strong nuclear force, whose operating length scales is found to increase with time to reach O(1) cm at the end of the luminous lifetimes of pulsars. The model presented here model is in line with the recent radio and GW observations of pulsars and NSs.


title: Variability of Ha emission in young stellar objects in the cluster IC 348

abstract: The aim of this work is to identify and study the stellar objects with variable EW(Ha) in the young stellar cluster IC 348. We performed photometric and slit-less observations at several epochs in order to reveal the variable objects. Significant variability of EW(Ha) was found in 90 out of 127 examined stars. From all epochs of observations, 32 objects were classified as CTT and 69 as WTT objects. The fraction of the variables in these samples is about 60%. We also identified 20 stellar objects, which showed not only a significant variability of the equivalent width, but which also change their apparent evolutionary stage. The analysis of data obtained over a wide wavelength range (from X-ray to mid-infrared) has shown that Ha activity and the measure of its variability are in good agreement with the activity of stellar objects measured with its other parameters, such as X-ray radiation and the mass accretion rate. The EW(Ha) differs not only between objects at different evolutionary stages, but also between variable and non-variable objects. The variables in the CTT and WTT samples are more active than non-variables although they have almost the same evolutionary age. Another distinct difference between these variables and non-variables is their average masses. The variables from both CTT and WTT samples are noticeably more massive than non-variables. Our data confirm the assumption made for other star formation regions, that the decay of accretion activity occurs more slowly for more massive CTT objects. Apparently, a similar trend is also present in WTT objects, which are at a later stage of evolution. The variability of the stellar objects, which change their evolutionary classes, at least in a fraction of them, is due to the fact that they are close binaries, which affects and modulates their Ha emission activity.


title: Large-Angle, Broadband and Multifunctional Gratings Based on Directively Radiating Waveguide Scatterers

abstract: Conventional surface-relief gratings are inefficient at deflecting normally-incident light by large angles. This constrains their use in many applications and limits the overall efficiency of any optical instrument integrating gratings. Here, we demonstrate a simple approach for the design of diffraction gratings that can be highly efficient for large deflection angles, while also offering additional functionality. The gratings are composed of a unit cell comprising a vertically-oriented asymmetric slot-waveguide. The unit cell shows oscillating unidirectional scattering behavior that can be precisely tuned as a function of the waveguide length. This occurs due to interference between multiple modes excited by the incident light. In contrast to metasurface-based gratings with multiple resonant sub-elements, a periodic arrangement of such non-resonant diffracting elements allows for broadband operation and a strong tolerance for variations in angle of incidence. Full-wave simulations show that our grating designs can exhibit diffraction efficiencies ranging from 94% for a deflection angle of 47$^\circ$ to 80% for deflection angle of 80$^\circ$. To demonstrate the multifunctionality of our grating design technique, we have also proposed a flat polarization beamsplitter, which allows for the separation of the two orthogonal polarizations by 80$^\circ$, with an efficiency of 80%.


title: Fractional Fermion Number and Hall Conductivity of Domain Walls

abstract: In this letter the fractional fermion number of thick domain walls is computed. The analysis is achieved by developing the heat kernel expansion of the spectral eta functon of the Dirac Hamiltonian governing the fermionic fluctuations around the domain wall. A formula is derived showing that a non null fermion number is always accompanied by a Hall conductivity induced on the wall. In the limit of thin and impenetrable walls the chiral bag boundary conditions arise, and the Hall conductivity is computed for this case as well.


title: Dynamics and Topology of Conformally Anosov Contact 3-Manifolds

abstract: We provide obstructions to the existence of conformally Anosov Reeb flows on a 3-manifold that partially generalize similar obstructions to Anosov Reeb flows. In particular, we show $\mathbb{S}^3$ does not admit conformally Anosov Reeb flows. We also give a Riemannian geometric condition on a metric compatible with a contact structure implying that a Reeb field is Anosov. From this we can give curvature conditions on a metric compatible with a contact structure that implies universal tightness of the contact structure among other things.


title: Prospects of designing gold-nanoparticles-based soft terahertz radiation sources and terahertz-to-infrared converters for concealed object detection technology

abstract: The two-phonon scheme of generation of terahertz (THz) photons by gold nanobars (GNBs) is considered. It is shown that in GNBs, by choosing their sizes, it is possible to provide conditions for converting the energy of longitudinal phonons with THz frequencies into the energy of THz photons. The prospects of designing GNBs-based soft THz radiation sources (frequencies: 0.14; 0.24; 0.41 and 0.70 THz) with a large flow cross-section (diameter ~40 cm) intended for detection of hidden objects under clothing to ensure security in public places (airports, railway stations, stadiums, etc.) are assessed. The choice of the above frequencies is a compromise between the requirements of low absorption of THz radiation by water vapor in air, good penetration through the fabric of clothing, favoring a sufficient resolution of the imaging system, and an abundance of corresponding longitudinal phonons, capable of exciting Fermi electrons in GNBs. Estimates of the characteristics of the terahertz-to-infrared converter based on gold nanospheres (GNSs), which could work in tandem with these sources of THz radiation -- as a means of visualization of hidden objects -- are also given.


title: Conductivity of the holographic p-wave superconductors with higher order corrections

abstract: We investigate the holographic $p$-wave superconductors in the presence of the higher order corrections on the gravity as well as on the gauge field side. On the gravity side, we add the Gauss-Bonnet curvature correction terms, while on the gauge field side we take the nonlinear Lagrangian in the form $F + b F^2$, where F is the Maxwell Lagrangian and b indicates the strength of nonlinearity. We employ the shooting method for the numerical calculations in order to obtain the ratio of the critical temperature $T_c$ over $\rho^ {1/(d-2)}$. We observe that by increasing the values of the mass and the nonlinear parameters the critical temperature decreases and thus the condensation becomes harder to form. In addition, the stronger Gauss-Bonnet parameter $\alpha$ hinders the superconducting phase in Gauss-Bonnet gravity. Furthermore, we calculate the electrical conductivity based on the holographic setup. The real and imaginary parts are related to each other by Kramers-Kronig relation which indicates a delta function and pole in low frequency regime, respectively. However, at enough large frequencies the trend of real part can be interpreted by $Re[\sigma] = \omega^{(d-4)}$. Moreover, in holographic model the ratio $\omega_g/T_c$ is always much larger than the BCS value $3.5$ due to the strong coupling of holographic superconductors. In both gravity kinds, decreasing the temperature or increasing the effect of nonlinearity shifts the gap frequency toward larger values. Besides, the gap frequency is occurred at larger values by enlarging the Gauss-Bonnet parameter. In general, the behavior of conductivity depends on the choice of the mass, the nonlinear and the Gauss-Bonnet parameters.


title: Multivaluedness Aspects in Self-organization, Complexity and Computations Investigations by Strong Anticipation

abstract: Since the introduction of strong anticipation by D.~Dubois the numerous investigations of concrete systems have been proposed. In proposed paper the new examples of discrete dynamical systems with anticipation are considered. The mathematical formulation of problems, possible analytical formulas for solutions and numerical examples of presumable solutions are proposed. One of the most interesting properties in such systems is presumable multivaluedness of the solutions. It can be considered from the point of view of dynamical chaos and complex behavior. We represent examples of periodic and complex solutions, attractor's properties and presumable applications in self-organization. The main peculiarity is the strong anticipation property. General new possibilities are the presumable multivaluedness of the dynamics of automata. Possible interpretations of such behavior of cellular automata are discussed. Further prospects for development of automata theory and hyper computation are proposed.


title: Chaos in balanced and unbalanced holographic s+p superconductors

abstract: In this work, we propose a toy model for a mixture of superconductors with competitive s and p modes using gauge/gravity duality. We demonstrate that the model undergoes phase transitions with the proper choice of different values for chemical potentials. We consider both balanced and unbalanced cases. We propose that the condensate field in the bulk toy model for a mixing s+p phase of high temperature superconductors undergoes a chaotic phase space scenario. Using a suitable measure function for chaos, we demonstrate the existence of chaotic dynamics via condensate fields. As an alternative method, we also investigate the phase portrait in the extended phase space for the condensates.


title: On the Magnetism of C14 Nb(Fe89.4Al10.6)2 Laves Phase Intermetallic Compound

abstract: C14 Nb(Fe89.4Al10.6)2 Laves phase intermetallic compound was investigated by DC magnetization (M) measurements performed in the temperature (T) interval of 20 to 175 K,under an applied magnetic field (H) ranging between 50 and 1250 Oe. Magnetization curves were recorded in the field-cooled (FC) and in the zero-field-cooled (ZFC) modes. They clearly showed an irreversible character that vanished at H=1250 Oe. Both magnetization curves exhibited well-defined peaks around T_N =72.3 K whose positions were H-independent, so they were identified as the compound's N\'eel temperature. The existence of irreversibility which decreases with H testify to a re-entrant character of magnetism in the studied compound. An increase of both MFC and MZFC observed below TN likely indicates a mixed i.e. ferromagnetic and antiferromagnetic ground magnetic state of the studied system.


title: Time Delay in the Swing Equation: A Variety of Bifurcations

abstract: The present paper addresses the swing equation with additional delayed damping as an example for pendulum-like systems. In this context, it is proved that recurring sub- and supercritical Hopf bifurcations occur if time delay is increased. To this end, a general formula for the first Lyapunov coefficient in second order systems with additional delayed damping and delay-free nonlinearity is given. In so far the paper extends results about stability switching of equilibria in linear time delay systems from Cooke and Grossman. In addition to the analytical results, periodic solutions are numerically dealt with. The numerical results demonstrate how a variety of qualitative behaviors is generated in the simple swing equation by only introducing time delay in a damping term.


title: Heterogeneities in electricity grids strongly enhance non-Gaussian features of frequency fluctuations under stochastic power input

abstract: Stochastic feed-in of fluctuating renewable energies is steadily increasing in modern electricity grids and this becomes an important risk factor for maintaining power grid stability. Here we study the impact of wind power feed-in on the short-term frequency fluctuations in power grids based on an IEEE test grid structure, the swing equation for the dynamics of voltage phase angles, and a series of measured wind speed data. External control measures are accounted for by adjusting the grid state to the average power feed-in on a time scale of one minute. The wind power is injected at a single node by replacing one of the conventional generator nodes in the test grid by a wind farm. We determine histograms of local frequencies for a large number of one-minute wind speed sequences taken from the measured data and for different injection nodes. These histograms exhibit a common type of shape, which can be described by a Gaussian distribution for small frequencies and a nearly exponentially decaying tail part. Non-Gaussian features become particularly pronounced for wind power injection at locations, which are weakly connected to the main grid structure. This effect is only present when taking into account the heterogeneities in transmission line and node properties of the grid, while it disappears upon homogenizing of these features. The standard deviation of the frequency fluctuations increases linearly with the average injected wind power.


title: Decision-facilitating information in hidden-action setups: An agent-based approach

abstract: The hidden-action model captures a fundamental problem of principal-agent theory and provides an optimal sharing rule when only the outcome but not the effort can be observed. However, the hidden-action model builds on various explicit and also implicit assumptions about the information of the contracting parties. This paper relaxes key assumptions regarding the availability of information included in the hidden-action model in order to study whether and, if so, how fast the optimal sharing rule is achieved and how this is affected by the various types of information employed in the principal-agent relation. Our analysis particularly focuses on information about the environment and about feasible actions for the agent. We follow an approach to transfer closed-form mathematical models into agent-based computational models and show that the extent of information about feasible options to carry out a task only has an impact on performance if decision makers are well informed about the environment, and that the decision whether to perform exploration or exploitation when searching for new feasible options only affects performance in specific situations. Having good information about the environment, on the contrary, appears to be crucial in almost all situations.


title: HATS: A Hierarchical Graph Attention Network for Stock Movement Prediction

abstract: Many researchers both in academia and industry have long been interested in the stock market. Numerous approaches were developed to accurately predict future trends in stock prices. Recently, there has been a growing interest in utilizing graph-structured data in computer science research communities. Methods that use relational data for stock market prediction have been recently proposed, but they are still in their infancy. First, the quality of collected information from different types of relations can vary considerably. No existing work has focused on the effect of using different types of relations on stock market prediction or finding an effective way to selectively aggregate information on different relation types. Furthermore, existing works have focused on only individual stock prediction which is similar to the node classification task. To address this, we propose a hierarchical attention network for stock prediction (HATS) which uses relational data for stock market prediction. Our HATS method selectively aggregates information on different relation types and adds the information to the representations of each company. Specifically, node representations are initialized with features extracted from a feature extraction module. HATS is used as a relational modeling module with initialized node representations. Then, node representations with the added information are fed into a task-specific layer. Our method is used for predicting not only individual stock prices but also market index movements, which is similar to the graph classification task. The experimental results show that performance can change depending on the relational data used. HATS which can automatically select information outperformed all the existing methods.


title: Bistability and Synchronization in Coupled Maps Models

abstract: Several coupled maps models are sketched and reviewed in this short communication. First, a discrete logistic type model that was proposed for the symbiotic interaction of two species. Second, a model of many of these symbiotic species evolving together under the coupling action of a global mean field. Third, a similar model with many of these species interacting in a symbiotic way with the nearest neighboring species in different kind of networks. The phenomena of bistability and synchronization is discussed in these models. Finally, some new results concerning the last model embedded in different topologies will be presented.


title: Route to hyperbolic hyperchaos in a nonautonomous time-delay system

abstract: We consider a self-oscillator whose excitation parameter is varied. Frequency of the variation is much smaller then the natural frequency of the oscillator so that oscillations in the system are periodically excited and decay. Also a time delay as added such that when the oscillations start to grow at a new excitation stage they are influenced via the delay line by the oscillations at the penultimate excitation stage. Due to a nonlinearity the seeding from the past arrives with a doubled phase so that oscillation phase changes from stage to stage according to chaotic Bernoulli-type map. As a result, the system operates as two coupled hyperbolic chaotic subsystems. Varying the relation between the delay time and the excitation period we affect a coupling strength between these subsystems as well as intensity of the phase doubling mechanism responsible for the hyperbolicity. Due to this, a transition from non-hyperbolic to hyperbolic hyperchaos occurs. The following steps of the transition scenario are revealed and analyzed: (a) an intermittency as alternation of long staying near a fixed point at the origin and short chaotic bursts; (b) chaotic oscillations with often visits to the fixed point; (c) plain hyperchaos without hyperbolicity after termination visiting the fixed point; (d) transformation of hyperchaos to hyperbolic form.


title: Branching random graph model of rough surfaces describes thermal properties of the effective molecular potential

abstract: Fluid properties near rough surfaces are crucial in describing fundamental surface phenomena and modern industrial material design implementations. One of the most powerful approaches to model real rough materials is based on the surface representation in terms of random geometry. Understanding the influence of random solid geometry on the low-temperature fluid thermodynamics is a cutting edge problem. Therefore this work extends recent studies bypassing high-temperature expansion and small heterogeneity scale. We introduce random branching trees whose topology reflects the hierarchical properties of a random solid geometry. This mathematical representation allows us to obtain averaged free energy using a statistical model of virtual clusters interacting through random ultrametric pairwise potentials. Our results demonstrate that a significant impact to fluid-solid interface energy is induced by the hierarchical structure of random geometry at low temperature. These calculations coincide with direct Monte Carlo simulations. Due to the study's interdisciplinary nature, the developed approach can be applied to a wide range of quenched disorder systems on random graphs.

categories: cond-mat.dis-nn cond-mat.mes-hall cond-mat.stat-mech math-ph math.MP


title: Enhancing quantum control by improving shape pulse generation

abstract: Most quantum processors requires pulse sequences for controlling quantum states. Here, we present an alternative algorithm for computing an optimal pulse sequence in order to perform a specific task, being an implementation of a quantum gate or a quantum state preparation. In our method, we reduced drastically the number of parameters to be fitted, by using a limited number of functions as the modulations for the amplitude and phase of the radio-frequency pulses, and employed approximations to make the algorithm fast and scalable. We demonstrate the success of the proposed algorithm, by performing several real experiments for 4, 7 and 12 quantum bits systems using NMR. In addition, we have also shown the efficiency of the algorithm, finding pulses for controlling with good fidelity the quantum states of spins in a fictional square bi-dimensional lattices containing 16, 36 and 100 qubits.


title: Pixel-wise Segmentation of Right Ventricle of Heart

abstract: One of the first steps in the diagnosis of most cardiac diseases, such as pulmonary hypertension, coronary heart disease is the segmentation of ventricles from cardiac magnetic resonance (MRI) images. Manual segmentation of the right ventricle requires diligence and time, while its automated segmentation is challenging due to shape variations and illdefined borders. We propose a deep learning based method for the accurate segmentation of right ventricle, which does not require post-processing and yet it achieves the state-of-the-art performance of 0.86 Dice coefficient and 6.73 mm Hausdorff distance on RVSC-MICCAI 2012 dataset. We use a novel adaptive cost function to counter extreme class-imbalance in the dataset. We present a comprehensive comparative study of loss functions, architectures, and ensembling techniques to build a principled approach for biomedical segmentation tasks.


title: Consistent Feature Construction with Constrained Genetic Programming for Experimental Physics

abstract: A good feature representation is a determinant factor to achieve high performance for many machine learning algorithms in terms of classification. This is especially true for techniques that do not build complex internal representations of data (e.g. decision trees, in contrast to deep neural networks). To transform the feature space, feature construction techniques build new high-level features from the original ones. Among these techniques, Genetic Programming is a good candidate to provide interpretable features required for data analysis in high energy physics. Classically, original features or higher-level features based on physics first principles are used as inputs for training. However, physicists would benefit from an automatic and interpretable feature construction for the classification of particle collision events. Our main contribution consists in combining different aspects of Genetic Programming and applying them to feature construction for experimental physics. In particular, to be applicable to physics, dimensional consistency is enforced using grammars. Results of experiments on three physics datasets show that the constructed features can bring a significant gain to the classification accuracy. To the best of our knowledge, it is the first time a method is proposed for interpretable feature construction with units of measurement, and that experts in high-energy physics validate the overall approach as well as the interpretability of the built features.


title: Evolutionary Computation, Optimization and Learning Algorithms for Data Science

abstract: A large number of engineering, science and computational problems have yet to be solved in a computationally efficient way. One of the emerging challenges is how evolving technologies grow towards autonomy and intelligent decision making. This leads to collection of large amounts of data from various sensing and measurement technologies, e.g., cameras, smart phones, health sensors, smart electricity meters, and environment sensors. Hence, it is imperative to develop efficient algorithms for generation, analysis, classification, and illustration of data. Meanwhile, data is structured purposefully through different representations, such as large-scale networks and graphs. We focus on data science as a crucial area, specifically focusing on a curse of dimensionality (CoD) which is due to the large amount of generated/sensed/collected data. This motivates researchers to think about optimization and to apply nature-inspired algorithms, such as evolutionary algorithms (EAs) to solve optimization problems. Although these algorithms look un-deterministic, they are robust enough to reach an optimal solution. Researchers do not adopt evolutionary algorithms unless they face a problem which is suffering from placement in local optimal solution, rather than global optimal solution. In this chapter, we first develop a clear and formal definition of the CoD problem, next we focus on feature extraction techniques and categories, then we provide a general overview of meta-heuristic algorithms, its terminology, and desirable properties of evolutionary algorithms.


title: Evolution of Ant Colony Optimization Algorithm -- A Brief Literature Review

abstract: Ant Colony Optimization (ACO) is a metaheuristic proposed by Marco Dorigo in 1991 based on behavior of biological ants. Pheromone laying and selection of shortest route with the help of pheromone inspired development of first ACO algorithm. Since, presentation of first such algorithm, many researchers have worked and published their research in this field. Though initial results were not so promising but recent developments have made this metaheuristic a significant algorithm in Swarm Intelligence. This research presents a brief overview of recent developments carried out in ACO algorithms in terms of both applications and algorithmic developments. For application developments, multi-objective optimization, continuous optimization and time-varying NP-hard problems have been presented. While to review articles based on algorithmic development, hybridization and parallel architectures have been investigated.


title: Small polarons and the Janus nature of $\text{TiO}_\text{2}(110)$

abstract: Polarons are ubiquitous in many semiconductors and have been linked with conductivity and optical response of materials for photovoltaics and heterogeneous catalysis, yet how surface polarons influence adsorption remains unclear. Here, by modelling the surface of rutile titania using density functional theory, we reveal the effect of small surface polarons on water adsorption, dissociation, and hydrogen bonding. On the one hand the presence of such polarons significantly suppresses dissociation of water molecules that are bonded directly to polaronic sites. On the other hand, polarons facilitate water dissociation at certain non-polaronic sites. Furthermore, polarons strengthen hydrogen bonds, which in turn affects water dissociation in hydrogen bonded overlayer structures. This study reveals that polarons at the rutile surface have complex, multi-faceted, effects on water adsorption, dissociation and hydrogen bonding, highlighting the importance of polarons on water structure and dynamics on such surfaces. We expect that many of the physical properties of surface polarons identified here will apply more generally to surfaces and interfaces that can host small polarons, beyond titania.


title: Network Time Synchronization of the Readout Electronics for a New Radioactive Gas Detection System

abstract: In systems with multiple radiation detectors, time synchronization of the data collected from different detectors is essential to reconstruct multi-detector events such as scattering and coincidences. In cases where the number of detectors exceeds the readout channels in a single data acquisition electronics module, multiple modules have to be synchronized, which is traditionally accomplished by distributing clocks and triggers via dedicated connections. To eliminate this added cabling complexity in the case of a new radioactive gas detection system prototype under development at the French Atomic Energy Commission, we implemented time synchronization between multiple XIA Pixie-Net detector readout modules through the existing Ethernet network, based on the IEEE 1588 precision time protocol. The detector system is dedicated to the measurement of radioactive gases at low activity and consists of eight large silicon pixels and two NaI(Tl) detectors, instrumented with a total of three 4-channel Pixie-Net modules. Detecting NaI(Tl)/silicon coincidences will make it possible to identify each radioisotope present in the sample. To allow these identifications at low activities, the Pixie-Net modules must be synchronized to a precision well below the targeted coincidence window of 500-1000 ns. Being equipped with an Ethernet PHY compatible with IEEE 1588 and synchronous Ethernet that outputs a locally generated but system-wide synchronized clock, the Pixie-Net can operate its analog to digital converters and digital processing circuitry with that clock and match time stamps for captured data across the three modules. Depending on the network configuration and synchronization method, the implementation is capable to achieve timing precisions between 300 ns and 200 ps.


title: Improving the Results of De novo Peptide Identification via Tandem Mass Spectrometry Using a Genetic Programming-based Scoring Function for Re-ranking Peptide-Spectrum Matches

abstract: De novo peptide sequencing algorithms have been widely used in proteomics to analyse tandem mass spectra (MS/MS) and assign them to peptides, but quality-control methods to evaluate the confidence of de novo peptide sequencing are lagging behind. A fundamental part of a quality-control method is the scoring function used to evaluate the quality of peptide-spectrum matches (PSMs). Here, we propose a genetic programming (GP) based method, called GP-PSM, to learn a PSM scoring function for improving the rate of confident peptide identification from MS/MS data. The GP method learns from thousands of MS/MS spectra. Important characteristics about goodness of the matches are extracted from the learning set and incorporated into the GP scoring functions. We compare GP-PSM with two methods including Support Vector Regression (SVR) and Random Forest (RF). The GP method along with RF and SVR, each is used for post-processing the results of peptide identification by PEAKS, a commonly used de novo sequencing method. The results show that GP-PSM outperforms RF and SVR and discriminates accurately between correct and incorrect PSMs. It correctly assigns peptides to 10% more spectra on an evaluation dataset containing 120 MS/MS spectra and decreases the false positive rate (FPR) of peptide identification.


title: A Fast and Efficient Stochastic Opposition-Based Learning for Differential Evolution in Numerical Optimization

abstract: A fast and efficient stochastic opposition-based learning (OBL) variant is proposed in this paper. OBL is a machine learning concept to accelerate the convergence of soft computing algorithms, which consists of simultaneously calculating an original solution and its opposite. Recently, a stochastic OBL variant called BetaCOBL was proposed, which is capable of controlling the degree of opposite solutions, preserving useful information held by original solutions, and preventing the waste of fitness evaluations. While it has shown outstanding performance compared to several state-of-the-art OBL variants, the high computational cost of BetaCOBL may hinder it from cost-sensitive optimization problems. Also, as it assumes that the decision variables of a given problem are independent, BetaCOBL may be ineffective for optimizing inseparable problems. In this paper, we propose an improved BetaCOBL that mitigates all the limitations. The proposed algorithm called iBetaCOBL reduces the computational cost from $O(NP^{2} \cdot D)$ to $O(NP \cdot D)$ ($NP$ and $D$ stand for population size and a dimension, respectively) using a linear time diversity measure. Also, the proposed algorithm preserves strongly dependent variables that are adjacent to each other using multiple exponential crossover. We used differential evolution (DE) variants to evaluate the performance of the proposed algorithm. The results of the performance evaluations on a set of 58 test functions show the excellent performance of iBetaCOBL compared to ten state-of-the-art OBL variants, including BetaCOBL.


title: Efficient training and design of photonic neural network through neuroevolution

abstract: Recently, optical neural networks (ONNs) integrated in photonic chips has received extensive attention because they are expected to implement the same pattern recognition tasks in the electronic platforms with high efficiency and low power consumption. However, the current lack of various learning algorithms to train the ONNs obstructs their further development. In this article, we propose a novel learning strategy based on neuroevolution to design and train the ONNs. Two typical neuroevolution algorithms are used to determine the hyper-parameters of the ONNs and to optimize the weights (phase shifters) in the connections. In order to demonstrate the effectiveness of the training algorithms, the trained ONNs are applied in the classification tasks for iris plants dataset, wine recognition dataset and modulation formats recognition. The calculated results exhibit that the training algorithms based on neuroevolution are competitive with other traditional learning algorithms on both accuracy and stability. Compared with previous works, we introduce an efficient training method for the ONNs and demonstrate their broad application prospects in pattern recognition, reinforcement learning and so on.


title: Adaptive Morley FEM for the von K\'{a}rm\'{a}n equations with optimal convergence rates

abstract: The adaptive nonconforming Morley finite element method (FEM) approximates a regular solution to the von K\'{a}rm\'{a}n equations with optimal convergence rates for sufficiently fine triangulations and small bulk parameter in the D\"orfler marking. This follows from the general axiomatic framework with the key arguments of stability, reduction, discrete reliability, and quasiorthogonality of an explicit residual-based error estimator. Particular attention is on the nonlinearity and the piecewise Sobolev embeddings required in the resulting trilinear form in the weak formulation of the nonconforming discretisation. The discrete reliability follows with a conforming companion for the discrete Morley functions from the medius analysis. The quasiorthogonality also relies on a novel piecewise $H^1$ a~priori error estimate and a careful analysis of the nonlinearity.


title: Graph based adaptive evolutionary algorithm for continuous optimization

abstract: he greatest weakness of evolutionary algorithms, widely used today, is the premature convergence due to the loss of population diversity over generations. To overcome this problem, several algorithms have been proposed, such as the Graph-based Evolutionary Algorithm (GEA) \cite{1} which uses graphs to model the structure of the population, but also memetic or differential evolution algorithms \cite{2,3}, or diversity-based ones \cite{4,5} have been designed. These algorithms are based on multi-populations, or often rather focus on the self-tuning parameters, however, they become complex to tune because of their high number of parameters. In this paper, our approach consists of an evolutionary algorithm that allows a dynamic adaptation of the search operators based on a graph in order to limit the loss of diversity and reduce the design complexity.


title: BGADAM: Boosting based Genetic-Evolutionary ADAM for Neural Network Optimization

abstract: For various optimization methods, gradient descent-based algorithms can achieve outstanding performance and have been widely used in various tasks. Among those commonly used algorithms, ADAM owns many advantages such as fast convergence with both the momentum term and the adaptive learning rate. However, since the loss functions of most deep neural networks are non-convex, ADAM also shares the drawback of getting stuck in local optima easily. To resolve such a problem, the idea of combining genetic algorithm with base learners is introduced to rediscover the best solutions. Nonetheless, from our analysis, the idea of combining genetic algorithm with a batch of base learners still has its shortcomings. The effectiveness of genetic algorithm can hardly be guaranteed if the unit models converge to close or the same solutions. To resolve this problem and further maximize the advantages of genetic algorithm with base learners, we propose to implement the boosting strategy for input model training, which can subsequently improve the effectiveness of genetic algorithm. In this paper, we introduce a novel optimization algorithm, namely Boosting based Genetic ADAM (BGADAM). With both theoretic analysis and empirical experiments, we will show that adding the boosting strategy into the BGADAM model can help models jump out the local optima and converge to better solutions.


title: Testing Robustness Against Unforeseen Adversaries

abstract: Adversarial robustness research primarily focuses on L_p perturbations, and most defenses are developed with identical training-time and test-time adversaries. However, in real-world applications developers are unlikely to have access to the full range of attacks or corruptions their system will face. Furthermore, worst-case inputs are likely to be diverse and need not be constrained to the L_p ball. To narrow in on this discrepancy between research and reality we introduce ImageNet-UA, a framework for evaluating model robustness against a range of unforeseen adversaries, including eighteen new non-L_p attacks. To perform well on ImageNet-UA, defenses must overcome a generalization gap and be robust to a diverse attacks not encountered during training. In extensive experiments, we find that existing robustness measures do not capture unforeseen robustness, that standard robustness techniques are beat by alternative training strategies, and that novel methods can improve unforeseen robustness. We present ImageNet-UA as a useful tool for the community for improving the worst-case behavior of machine learning systems.


title: Knowledge transfer in deep block-modular neural networks

abstract: Although deep neural networks (DNNs) have demonstrated impressive results during the last decade, they remain highly specialized tools, which are trained -- often from scratch -- to solve each particular task. The human brain, in contrast, significantly re-uses existing capacities when learning to solve new tasks. In the current study we explore a block-modular architecture for DNNs, which allows parts of the existing network to be re-used to solve a new task without a decrease in performance when solving the original task. We show that networks with such architectures can outperform networks trained from scratch, or perform comparably, while having to learn nearly 10 times fewer weights than the networks trained from scratch.


title: Exploiting a Stimuli Encoding Scheme of Spiking Neural Networks for Stream Learning

abstract: Stream data processing has gained progressive momentum with the arriving of new stream applications and big data scenarios. One of the most promising techniques in stream learning is the Spiking Neural Network, and some of them use an interesting population encoding scheme to transform the incoming stimuli into spikes. This study sheds lights on the key issue of this encoding scheme, the Gaussian receptive fields, and focuses on applying them as a pre-processing technique to any dataset in order to gain representativeness, and to boost the predictive performance of the stream learning methods. Experiments with synthetic and real data sets are presented, and lead to confirm that our approach can be applied successfully as a general pre-processing technique in many real cases.


title: Spiking Neural Networks and Online Learning: An Overview and Perspectives

abstract: Applications that generate huge amounts of data in the form of fast streams are becoming increasingly prevalent, being therefore necessary to learn in an online manner. These conditions usually impose memory and processing time restrictions, and they often turn into evolving environments where a change may affect the input data distribution. Such a change causes that predictive models trained over these stream data become obsolete and do not adapt suitably to new distributions. Specially in these non-stationary scenarios, there is a pressing need for new algorithms that adapt to these changes as fast as possible, while maintaining good performance scores. Unfortunately, most off-the-shelf classification models need to be retrained if they are used in changing environments, and fail to scale properly. Spiking Neural Networks have revealed themselves as one of the most successful approaches to model the behavior and learning potential of the brain, and exploit them to undertake practical online learning tasks. Besides, some specific flavors of Spiking Neural Networks can overcome the necessity of retraining after a drift occurs. This work intends to merge both fields by serving as a comprehensive overview, motivating further developments that embrace Spiking Neural Networks for online learning scenarios, and being a friendly entry point for non-experts.


title: Comparing reliability of grid-based Quality-Diversity algorithms using artificial landscapes

abstract: Quality-Diversity (QD) algorithms are a recent type of optimisation methods that search for a collection of both diverse and high performing solutions. They can be used to effectively explore a target problem according to features defined by the user. However, the field of QD still does not possess extensive methodologies and reference benchmarks to compare these algorithms. We propose a simple benchmark to compare the reliability of QD algorithms by optimising the Rastrigin function, an artificial landscape function often used to test global optimisation methods.


title: Reservoir-size dependent learning in analogue neural networks

abstract: The implementation of artificial neural networks in hardware substrates is a major interdisciplinary enterprise. Well suited candidates for physical implementations must combine nonlinear neurons with dedicated and efficient hardware solutions for both connectivity and training. Reservoir computing addresses the problems related with the network connectivity and training in an elegant and efficient way. However, important questions regarding impact of reservoir size and learning routines on the convergence-speed during learning remain unaddressed. Here, we study in detail the learning process of a recently demonstrated photonic neural network based on a reservoir. We use a greedy algorithm to train our neural network for the task of chaotic signals prediction and analyze the learning-error landscape. Our results unveil fundamental properties of the system's optimization hyperspace. Particularly, we determine the convergence speed of learning as a function of reservoir size and find exceptional, close to linear scaling. This linear dependence, together with our parallel diffractive coupling, represent optimal scaling conditions for our photonic neural network scheme.


title: Genetic Algorithm for the 0/1 Multidimensional Knapsack Problem

abstract: The 0/1 multidimensional knapsack problem is the 0/1 knapsack problem with m constraints which makes it difficult to solve using traditional methods like dynamic programming or branch and bound algorithms. We present a genetic algorithm for the multidimensional knapsack problem with Java and C++ code that is able to solve publicly available instances in a very short computational duration. Our algorithm uses iteratively computed Lagrangian multipliers as constraint weights to augment the greedy algorithm for the multidimensional knapsack problem and uses that information in a greedy crossover in a genetic algorithm. The algorithm uses several other hyperparameters which can be set in the code to control convergence. Our algorithm improves upon the algorithm by Chu and Beasley in that it converges to optimum or near optimum solutions much faster.


title: A NICER Look at Strong X-ray Obscuration in the Seyfert-2 Galaxy NGC 4388

abstract: We present an analysis of the time-averaged spectrum of the Seyfert-2 active galaxy NGC 4388, obtained by NICER. The intrinsic strength of the reflection spectrum in NGC 4388, the large collecting area and favorable pass band of NICER, and a net exposure of 105.6 ks yielded an exceptionally sensitive spectrum. Using two independent families of models, the intrinsic spectrum from the central engine is found to be highly obscured but not Compton-thick. Enforcing physical self-consistency within each model, the independent treatments give formally consistent results: N_H = 2.67 (-0.03,+0.02) E+23 cm^-2 or N_H = 2.64 (-0.03, +0.03) E+23 cm^-2. Past measurements made with Suzaku and XMM-Newton are in broad agreement with these column density values. A more recent measurement with NuSTAR (in late 2013) recorded a column density about twice as large; the robustness of this variability is reinforced by the use of consistent models and procedures. The neutral Fe K-alpha line in the NICER spectrum is nominally resolved and consistent with an origin in the optical broad line region (BLR). The data also require ionized absorption in the Fe K band, similar to the "warm absorbers" detected in Seyfert-1 active galactic nuclei (AGN). The low-energy spectrum is consistent with a set of ionized plasma components. We discuss these findings and note that the geometric inferences that derive from this analysis can be tested with XRISM and Athena.


title: Mapping of Local and Global Synapses on Spiking Neuromorphic Hardware

abstract: Spiking Neural Networks (SNNs) are widely deployed to solve complex pattern recognition, function approximation and image classification tasks. With the growing size and complexity of these networks, hardware implementation becomes challenging because scaling up the size of a single array (crossbar) of fully connected neurons is no longer feasible due to strict energy budget. Modern neromorphic hardware integrates small-sized crossbars with time-multiplexed interconnects. Partitioning SNNs becomes essential in order to map them on neuromorphic hardware with the major aim to reduce the global communication latency and energy overhead. To achieve this goal, we propose our instantiation of particle swarm optimization, which partitions SNNs into local synapses (mapped on crossbars) and global synapses (mapped on time-multiplexed interconnects), with the objective of reducing spike communication on the interconnect. This improves latency, power consumption as well as application performance by reducing inter-spike interval distortion and spike disorders. Our framework is implemented in Python, interfacing CARLsim, a GPU-accelerated application-level spiking neural network simulator with an extended version of Noxim, for simulating time-multiplexed interconnects. Experiments are conducted with realistic and synthetic SNN-based applications with different computation models, topologies and spike coding schemes. Using power numbers from in-house neuromorphic chips, we demonstrate significant reductions in energy consumption and spike latency over PACMAN, the widely-used partitioning technique for SNNs on SpiNNaker.


title: WikiCREM: A Large Unsupervised Corpus for Coreference Resolution

abstract: Pronoun resolution is a major area of natural language understanding. However, large-scale training sets are still scarce, since manually labelling data is costly. In this work, we introduce WikiCREM (Wikipedia CoREferences Masked) a large-scale, yet accurate dataset of pronoun disambiguation instances. We use a language-model-based approach for pronoun resolution in combination with our WikiCREM dataset. We compare a series of models on a collection of diverse and challenging coreference resolution problems, where we match or outperform previous state-of-the-art approaches on 6 out of 7 datasets, such as GAP, DPR, WNLI, PDP, WinoBias, and WinoGender. We release our model to be used off-the-shelf for solving pronoun disambiguation.


title: Refactoring Neural Networks for Verification

abstract: Deep neural networks (DNN) are growing in capability and applicability. Their effectiveness has led to their use in safety critical and autonomous systems, yet there is a dearth of cost-effective methods available for reasoning about the behavior of a DNN. In this paper, we seek to expand the applicability and scalability of existing DNN verification techniques through DNN refactoring. A DNN refactoring defines (a) the transformation of the DNN's architecture, i.e., the number and size of its layers, and (b) the distillation of the learned relationships between the input features and function outputs of the original to train the transformed network. Unlike with traditional code refactoring, DNN refactoring does not guarantee functional equivalence of the two networks, but rather it aims to preserve the accuracy of the original network while producing a simpler network that is amenable to more efficient property verification. We present an automated framework for DNN refactoring, and demonstrate its potential effectiveness through three case studies on networks used in autonomous systems.


title: IIB flux non-commutativity and the global structure of field theories

abstract: We discuss the origin of the choice of global structure for six dimensional $(2,0)$ theories and their compactifications in terms of their realization from IIB string theory on ALE spaces. We find that the ambiguity in the choice of global structure on the field theory side can be traced back to a subtle effect that needs to be taken into account when specifying boundary conditions at infinity in the IIB orbifold, namely the known non-commutativity of RR fluxes in spaces with torsion. As an example, we show how the classification of $\mathcal{N}=4$ theories by Aharony, Seiberg and Tachikawa can be understood in terms of choices of boundary conditions for RR fields in IIB. Along the way we encounter a formula for the fractional instanton number of $\mathcal{N}=4$ ADE theories in terms of the torsional linking pairing for rational homology spheres. We also consider six-dimensional $(1,0)$ theories, clarifying the rules for determining commutators of flux operators for discrete 2-form symmetries. Finally, we analyze the issue of global structure for four dimensional theories in the presence of duality defects.


title: Generating photon-added states without adding a photon

abstract: We show that a continuous range of nonclassical states of light can be generated using conditional measurements on the idler mode of an optical parametric amplifier. The output state is prepared by introducing a coherent state in the signal mode of the amplifier with a single photon in the idler mode, followed by a conditional measurement of a single photon in the output idler mode. By varying the gain of the amplifier, this approach can produce a coherent state, a photon-added state, a displaced number state, or a continuous range of other nonclassical states with intermediate properties. We note that this approach can generate a photon-added state even though the post-selected amplifier does not add any photons to the signal or idler modes. The ability to generate a continuous range of nonclassical states may have practical applications in quantum information processing.


title: Analog for the Wiener Lemma for Wolff-Denjoy Series

abstract: Let a function f with real poles be expanded in a Wolff-Denjoy series with positive coefficients. The main result of the note states that if we subtract its linear part from the function 1/f, then the remaining fractional part of this function will also expand into Wolff-Denjoy series (its poles are also real, and the coefficients of the series are negative). In other words, for Wolff-Denjoy series of the indicated form, an analogue of the well-known Wiener lemma in the theory of Fourier series is true up to a linear term. Applications of the result to operator theory are given.


title: Integrability and cycles of deformed ${\cal N}=2$ gauge theory

abstract: To analyse pure ${\cal N}=2$ $SU(2)$ gauge theory in the Nekrasov-Shatashvili (NS) limit (or deformed Seiberg-Witten (SW)), we use the Ordinary Differential Equation/Integrable Model (ODE/IM) correspondence, and in particular its (broken) discrete symmetry in its extended version with {\it two} singular irregular points. Actually, this symmetry appears to be 'manifestation' of the spontaneously broken $\mathbb{Z}_2$ R-symmetry of the original gauge problem and the two deformed SW cycles are simply connected to the Baxter's $T$ and $Q$ functions, respectively, of the Liouville conformal field theory at the self-dual point. The liaison is realised via a second order differential operator which is essentially the 'quantum' version of the square of the SW differential. Moreover, the constraints imposed by the broken $\mathbb{Z}_2$ R-symmetry acting on the moduli space (Bilal-Ferrari equations) seem to have their quantum counterpart in the $TQ$ and the $T$ periodicity relations, and integrability yields also a useful Thermodynamic Bethe Ansatz (TBA) for the cycles ($Y(\theta,\pm u)$ or their square roots, $Q(\theta,\pm u)$). A latere, two efficient asymptotic expansion techniques are presented. Clearly, the whole construction is extendable to gauge theories with matter and/or higher rank groups.


title: MuSHR: A Low-Cost, Open-Source Robotic Racecar for Education and Research

abstract: We present MuSHR, the Multi-agent System for non-Holonomic Racing. MuSHR is a low-cost, open-source robotic racecar platform for education and research, developed by the Personal Robotics Lab in the Paul G. Allen School of Computer Science & Engineering at the University of Washington. MuSHR aspires to contribute towards democratizing the field of robotics as a low-cost platform that can be built and deployed by following detailed, open documentation and do-it-yourself tutorials. A set of demos and lab assignments developed for the Mobile Robots course at the University of Washington provide guided, hands-on experience with the platform, and milestones for further development. MuSHR is a valuable asset for academic research labs, robotics instructors, and robotics enthusiasts.


title: Kinetic turbulence in shining pair plasma: intermittent beaming and thermalization by radiative cooling

abstract: High-energy astrophysical systems frequently contain collisionless relativistic plasmas that are heated by turbulent cascades and cooled by emission of radiation. Understanding the nature of this radiative turbulence is a frontier of extreme plasma astrophysics. In this paper, we use particle-in-cell simulations to study the effects of external inverse Compton radiation on turbulence driven in an optically thin, relativistic pair plasma. We focus on the statistical steady state (where injected energy is balanced by radiated energy) and perform a parameter scan spanning from low magnetization to high magnetization ($0.04 \lesssim \sigma \lesssim 11$). We demonstrate that the global particle energy distributions are quasi-thermal in all simulations, with only a modest population of nonthermal energetic particles (extending the tail by a factor of $\sim 2$). This indicates that nonthermal particle acceleration (observed in similar non-radiative simulations) is quenched by strong radiative cooling. The quasi-thermal energy distributions are well fit by analytic models in which stochastic particle acceleration (due to, e.g., second-order Fermi mechanism or gyroresonant interactions) is balanced by the radiation reaction force. Despite the efficient thermalization of the plasma, nonthermal energetic particles do make a conspicuous appearance in the anisotropy of the global momentum distribution as highly variable, intermittent beams (for high magnetization cases). The beamed high-energy particles are spatially coincident with intermittent current sheets, suggesting that localized magnetic reconnection may be a mechanism for kinetic beaming. This beaming phenomenon may explain rapid flares observed in various astrophysical systems (such as blazar jets, the Crab nebula, and Sagittarius A*).


title: UV considerations on scattering amplitudes in a web of theories

abstract: The scattering predictions of a web of theories including Yang-Mills (YM), gravity, bi-adjoint scalar, the non-linear sigma model (NLSM), Dirac-Born-Infeld-Volkov-Akulov (DBI-VA) and the special Galileon (sGal) form a class of special objects with two fascinating properties: they are related by the double-copy procedure, and they can be defined purely by on-shell constraints. We expand on both of these properties. First we show that NLSM tree-level amplitudes are fully determined by imposing color-dual structure together with cyclic invariance and locality. We then consider how hard-scaling can be used to constrain the predictions of these theories, as opposed to the usual soft-scaling. We probe the UV by generalizing the familiar BCFW shift off-shell to a novel single hard limit. We show that UV scalings are sufficient to fully constrain: 1. Bi-adjoint doubly-ordered amplitudes, assuming locality; 2. NLSM and BI, assuming locality and unitarity; 3. Special Galileon, assuming locality, unitarity, and a UV bound for the general Galileon vertex. We see how potentially distinct aspects of this UV behavior can be understood and unified via double-copy relations. Surprisingly, we find evidence that assuming unitarity for these theories may not be necessary, and can emerge via UV considerations and locality alone. These results complete the observations that, like IR considerations, UV scaling is sufficient to fully constrain a wide range of tree-level amplitudes, for both gauge, gravity, and effective field theories.


title: Good Fibrations through the Modal Prism

abstract: Homotopy type theory is a formal language for doing abstract homotopy theory -- the study of identifications. But in unmodified homotopy type theory, there is no way to say that these identifications come from identifying the path-connected points of a space. In other words, we can do abstract homotopy theory, but not algebraic topology. Shulman's Real Cohesive HoTT remedies this issue by introducing a system of modalities that relate the spatial structure of types to their homotopical structure. In this paper, we develop a theory of modal fibrations for a general modality, and apply it in particular to the shape modality of Real Cohesion. We then give examples of modal fibrations in Real Cohesive HoTT, and develop the theory of covering spaces.


title: More unlabelled data or label more data? A study on semi-supervised laparoscopic image segmentation

abstract: Improving a semi-supervised image segmentation task has the option of adding more unlabelled images, labelling the unlabelled images or combining both, as neither image acquisition nor expert labelling can be considered trivial in most clinical applications. With a laparoscopic liver image segmentation application, we investigate the performance impact by altering the quantities of labelled and unlabelled training data, using a semi-supervised segmentation algorithm based on the mean teacher learning paradigm. We first report a significantly higher segmentation accuracy, compared with supervised learning. Interestingly, this comparison reveals that the training strategy adopted in the semi-supervised algorithm is also responsible for this observed improvement, in addition to the added unlabelled data. We then compare different combinations of labelled and unlabelled data set sizes for training semi-supervised segmentation networks, to provide a quantitative example of the practically useful trade-off between the two data planning strategies in this surgical guidance application.


title: Deep Reinforcement Learning for Foreign Exchange Trading

abstract: Reinforcement learning can interact with the environment and is suitable for applications in decision control systems. Therefore, we used the reinforcement learning method to establish a foreign exchange transaction, avoiding the long-standing problem of unstable trends in deep learning predictions. In the system design, we optimized the Sure-Fire statistical arbitrage policy, set three different actions, encoded the continuous price over a period of time into a heat-map view of the Gramian Angular Field (GAF) and compared the Deep Q Learning (DQN) and Proximal Policy Optimization (PPO) algorithms. To test feasibility, we analyzed three currency pairs, namely EUR/USD, GBP/USD, and AUD/USD. We trained the data in units of four hours from 1 August 2018 to 30 November 2018 and tested model performance using data between 1 December 2018 and 31 December 2018. The test results of the various models indicated that favorable investment performance was achieved as long as the model was able to handle complex and random processes and the state was able to describe the environment, validating the feasibility of reinforcement learning in the development of trading strategies.


title: Hebbian Graph Embeddings

abstract: Representation learning has recently been successfully used to create vector representations of entities in language learning, recommender systems and in similarity learning. Graph embeddings exploit the locality structure of a graph and generate embeddings for nodes which could be words in a language, products of a retail website; and the nodes are connected based on a context window. In this paper, we consider graph embeddings with an error-free associative learning update rule, which models the embedding vector of node as a non-convex Gaussian mixture of the embeddings of the nodes in its immediate vicinity with some constant variance that is reduced as iterations progress. It is very easy to parallelize our algorithm without any form of shared memory, which makes it possible to use it on very large graphs with a much higher dimensionality of the embeddings. We study the efficacy of proposed method on several benchmark data sets and favorably compare with state of the art methods. Further, proposed method is applied to generate relevant recommendations for a large retailer.