title: Transforming Data Flow Diagrams for Privacy Compliance (Long Version)

abstract: Recent regulations, such as the European General Data Protection Regulation (GDPR), put stringent constraints on the handling of personal data. Privacy, like security, is a non-functional property, yet most software design tools are focused on functional aspects, using for instance Data Flow Diagrams (DFDs). In previous work, a conceptual model was introduced where DFDs could be extended into so-called Privacy-Aware Data Flow Diagrams (PA-DFDs) with the aim of adding specific privacy checks to existing DFDs. In this paper, we provide an explicit algorithm and a proof-of-concept implementation to transform DFDs into PA-DFDs. Our tool assists software engineers in the critical but error-prone task of systematically inserting privacy checks during design (they are automatically added by our tool) while still allowing them to inspect and edit the. PA-DFD if necessary. We have also identified and addressed ambiguities and inaccuracies in the high-level transformation proposed in previous work. We apply our approach to two realistic applications from the construction and online retail sectors.


title: Normal trace for vector fields of bounded mean oscillation

abstract: We introduce various spaces of vector fields of bounded mean oscillation ($BMO$) defined in a domain so that normal trace on the boundary is bounded when its divergence is well controlled. The behavior of "normal" component and "tangential" component may be different for our $BMO$ vector fields. As a result zero extension of the normal component stays in $BMO$ although such property may not hold for tangential components.


title: Quantum corrections to slow-roll inflation: scalar and tensor modes

abstract: Inflation is often described through the dynamics of a scalar field, slow-rolling in a suitable potential. Ultimately, this inflaton must be identified as the expectation value of a quantum field, evolving in a quantum effective potential. The shape of this potential is determined by the underlying tree-level potential, dressed by quantum corrections from the scalar field itself and the metric perturbations. Following [1], we compute the effective scalar field equations and the corrected Friedmann equations to quadratic order in both scalar field, scalar metric and tensor perturbations. We identify the quantum corrections from different sources at leading order in slow-roll, and estimate their magnitude in benchmark models of inflation. We comment on the implications of non-minimal coupling to gravity in this context.


title: Frequency theorem for parabolic equations and its relation to inertial manifolds theory

abstract: We obtain a version of the Frequency Theorem (a theorem on solvability of certain operator inequalities), which allows to construct quadratic Lyapunov functionals for semilinear parabolic equations. We show that the well-known Spectral Gap Condition, which was used in the theory of inertial manifolds by C. Foias, R. Temam and G. R. Sell, is a particular case of some frequency inequality, which arises within the Frequency Theorem. In particular, this allows to construct inertial manifolds for semilinear parabolic equations (including also some non-autonomous problems) in the context of a more general geometric theory developed in our adjacent works. This theory is based on quadratic Lyapunov functionals and generalizes the frequency-domain approach used by R. A. Smith. We also discuss the optimality of frequency inequalities and its relationship with known old and recent results in the field.


title: Revisiting Pixel-Wise Supervision for Face Anti-Spoofing

abstract: Face anti-spoofing (FAS) plays a vital role in securing face recognition systems from the presentation attacks (PAs). As more and more realistic PAs with novel types spring up, it is necessary to develop robust algorithms for detecting unknown attacks even in unseen scenarios. However, deep models supervised by traditional binary loss (e.g., `0' for bonafide vs. `1' for PAs) are weak in describing intrinsic and discriminative spoofing patterns. Recently, pixel-wise supervision has been proposed for the FAS task, intending to provide more fine-grained pixel/patch-level cues. In this paper, we firstly give a comprehensive review and analysis about the existing pixel-wise supervision methods for FAS. Then we propose a novel pyramid supervision, which guides deep models to learn both local details and global semantics from multi-scale spatial context. Extensive experiments are performed on five FAS benchmark datasets to show that, without bells and whistles, the proposed pyramid supervision could not only improve the performance beyond existing pixel-wise supervision frameworks, but also enhance the model's interpretability (i.e., locating the patch-level positions of PAs more reasonably). Furthermore, elaborate studies are conducted for exploring the efficacy of different architecture configurations with two kinds of pixel-wise supervisions (binary mask and depth map supervisions), which provides inspirable insights for future architecture/supervision design.


title: A fixed-point approach for decaying solutions of difference equations

abstract: A boundary value problem associated to the difference equation with advanced argument \begin{equation} \label{*}\Delta\bigl (a_{n}\Phi(\Delta x_{n})\bigr)+b_{n}\Phi(x_{n+p} )=0,\ \ n\geq1 \tag{$*$} \end{equation} is presented, where $\Phi(u)=|u|^{\alpha}$sgn $u,$ $\alpha>0,p$ is a positive integer and the sequences $a,b,$ are positive. We deal with a particular type of decaying solutions of (\ref{*}), that is the so-called intermediate solutions (see below for the definition) . In particular, we prove the existence of these type of solutions for (\ref{*}) by reducing it to a suitable boundary value problem associated to a difference equation without deviating argument. Our approach is based on a fixed point result for difference equations, which originates from existing ones stated in the continuous case. Some examples and suggestions for future researches complete the paper.


title: Weil--Petersson geodesics on the modular surface

abstract: We consider the Weil--Petersson (WP) metric on the modular surface. We lift WP geodesics to the universal cover of the modular surface and analyse geometric properties of the lifts as paths in the hyperbolic metric on the universal cover. For any pair of distinct points in the thick part of the universal cover, we prove that the WP and hyperbolic geodesic segments that connect the pair, fellow-travel in the thick part and all deviations between these segments arise during cusp excursions. Furthermore, we give a quantitative analysis of the deviation during an excursion. We leverage the fellow traveling to derive a correspondence between recurrent WP and hyperbolic geodesic rays from a base-point. We show that the correspondence can be promoted to a homeomorphism on the circle of directions. By analysing cuspidal winding of a typical WP geodesic ray, we show that the homeomorphism pushes forward a Lebesgue measure on the circle to a singular measure. In terms of continued fraction coefficients, the singularity boils down to a comparison that we prove, namely, the average coefficient is bounded along a typical WP ray but unbounded along a typical hyperbolic ray.


title: Low-Power IoT Communication Security: On the Performance of DTLS and TLS 1.3

abstract: Similarly to elsewhere on the Internet, practical security in the Internet of Things (IoT) is achieved by combining an array of mechanisms, at work at all layers of the protocol stack, in system software, and in hardware. Standard protocols such as Datagram Transport Layer Security (DTLS 1.2) and Transport Layer Security (TLS 1.2) are often recommended to secure communications to/from IoT devices. Recently, the TLS 1.3 standard was released and DTLS 1.3 is in the final stages of standardization. In this paper, we give an overview of version 1.3 of these protocols, and we provide the first experimental comparative performance analysis of different implementations and various configurations of these protocols, on real IoT devices based on low-power microcontrollers. We show how different implementations lead to different compromises. We measure and compare bytes-over-the-air, memory footprint, and energy consumption. We show that, when DTLS/TLS 1.3 requires more resources than DTLS/TLS 1.2, this additional overhead is quite reasonable. We also observe that, in some configurations, DTLS/TLS 1.3 actually decreases overhead and resource consumption. All in all, our study indicates that there is still room to optimize the existing implementations of these protocols.


title: Adaptive Smoothing Spline Estimator for the Function-on-Function Linear Regression Model

abstract: In this paper, we propose an adaptive smoothing spline (AdaSS) estimator for the function-on-function linear regression model where each value of the response, at any domain point, depends on the full trajectory of the predictor. The AdaSS estimator is obtained by the optimization of an objective function with two spatially adaptive penalties, based on initial estimates of the partial derivatives of the regression coefficient function. This allows the proposed estimator to adapt more easily to the true coefficient function over regions of large curvature and not to be undersmoothed over the remaining part of the domain. A novel evolutionary algorithm is developed ad hoc to obtain the optimization tuning parameters. Extensive Monte Carlo simulations have been carried out to compare the AdaSS estimator with competitors that have already appeared in the literature before. The results show that our proposal mostly outperforms the competitor in terms of estimation and prediction accuracy. Lastly, those advantages are illustrated also on two real-data benchmark examples.


title: Experimental measurement of the divergent quantum metric of an exceptional point

abstract: The geometry of Hamiltonian's eigenstates is encoded in the quantum geometric tensor (QGT). It contains both the Berry curvature, central to the description of topological matter and the quantum metric. So far the full QGT has been measured only in Hermitian systems, where the role of the quantum metric is mostly shown to determine corrections to physical effects. On the contrary, in non-Hermitian systems, and in particular near exceptional points, the quantum metric is expected to diverge and to often play a dominant role, for example on the enhanced sensing and on wave packet dynamics. In this work, we report the first experimental measurement of the quantum metric in a non-Hermitian system. The specific platform under study is an organic microcavity with exciton-polariton eigenstates, which demonstrate exceptional points. We measure the quantum metric's divergence and we determine the scaling exponent $n=-1.01\pm0.08$, which is in agreement with theoretical predictions for the second-order exceptional points.


title: On weak metric dimension of digraphs

abstract: Using the two way distance, we introduce the concepts of weak metric dimension of a strongly connected digraph $\Gamma$. We first establish lower and upper bounds for the number of arcs in $\Gamma$ by using the diameter and weak metric dimension of $\Gamma$, and characterize all digraphs attaining the lower or upper bound. Then we study a digraph with weak metric dimension $1$ and classify all vertex-transitive digraphs having weak metric dimension $1$. Finally, all digraphs of order $n$ with weak metric dimension $n-1$ or $n-2$ are determined.


title: Effect of inter-system crossing rates and optical illumination on the polarization of nuclear spins nearby nitrogen-vacancy centers

abstract: Several efforts have been made to polarize the nearby nuclear environment of nitrogen vacancy (NV) centers for quantum metrology and quantum information applications. Different methods showed different nuclear spin polarization efficiencies and rely on electronic spin polarization associated to the NV center, which in turn crucially depends on the inter-system crossing. Recently, the rates involved in the inter-system crossing have been measured leading to different transition rate models. Here, we consider the effect of these rates on several nuclear polarization methods based on the level anti-crossing, and precession of the nuclear population while the electronic spin is in the ms = 0 and ms = 1 spin states. We show that the nuclear polarization depends on the power of optical excitation used to polarize the electronic spin. The degree of nuclear spin polarization is different for each transition rate model. Therefore, the results presented here are relevant for validating these models and for polarizing nuclear spins. Furthermore, we analyze the performance of each method by considering the nuclear position relative to the symmetry axis of the NV center.


title: New method of verifying cryptographic protocols based on the process model

abstract: A cryptographic protocol (CP) is a distributed algorithm designed to provide a secure communication in an insecure environment. CPs are used, for example, in electronic payments, electronic voting procedures, database access systems, etc. Errors in the CPs can lead to great financial and social damage, therefore it is necessary to use mathematical methods to justify the correctness and safety of the CPs. In this paper, a new mathematical model of a CP is introduced, which allows one to describe both the CPs and their properties. It is shown how, on the base of this model, it is possible to solve the problems of verification of CPs.


title: Measurement of the CKM angle $\gamma$ and $B^0_s$-$\bar{B}^0_s$ mixing frequency with $B^0_s \rightarrow D_s^\mp h^\pm \pi^\pm \pi^\mp$ decays

abstract: The CKM angle $\gamma$ is measured for the first time from mixing-induced $CP$ violation between $B^0_s \rightarrow D_s^\mp K^\pm \pi^\pm \pi^\mp$ and $\bar{B}^0_s \rightarrow D_s^\pm K^\mp \pi^\mp \pi^\pm$ decays reconstructed in proton-proton collision data corresponding to an integrated luminosity of 9 ${\rm fb}^{-1}$ recorded with the LHCb detector. A time-dependent amplitude analysis is performed to extract the $CP$-violating weak phase $\gamma-2\beta_s$ and, subsequently, $\gamma$ by taking the $B^0_s$-$\bar{B}^0_s$ mixing phase $\beta_{s}$ as an external input. The measurement yields $\gamma = (44 \pm 12)^\circ$ modulo $180^\circ$, where statistical and systematic uncertainties are combined. An alternative model-independent measurement, integrating over the five-dimensional phase space of the decay, yields $\gamma = (44^{\,+\,20}_{\,-\,13})^\circ$ modulo $180^\circ$. Moreover, the $B^0_s$-$\bar{B}^0_s$ oscillation frequency is measured from the flavour-specific control channel $B^0_s \rightarrow D_s^- \pi^+ \pi^+ \pi^-$ to be $\Delta m_s = (17.757 \pm 0.007 \,({\rm stat.}) \pm 0.008 \,({\rm syst.})) \text{ps}^{-1}$, consistent with and more precise than the current world-average value.


title: Young stellar population gradients in central cluster galaxies from NUV and optical spectroscopy

abstract: Central cluster galaxies are the largest and most massive galaxies in the Universe. Although they host very old stellar populations, several studies found the existence of blue cores in some BCGs indicating ongoing star formation. We analyse VLT/X-Shooter stacked spectra of 6 nearby massive central galaxies with high central velocity dispersions ($\sigma$>300 km/s) at different galactocentric distances. We quantify the young stellar population out to 4 kpc by fitting near-UV and optical absorption line indices with predictions of composite stellar populations modelled by an old and a young stellar component. We also use IMF-sensitive indices since these galaxies have been found to host a bottom-heavy IMF in their central regions. We derive negative young stellar populations gradients, with mass fractions of stars younger than 1 Gyr decreasing with galactocentric distance, from 0.70% within 0.8 kpc to zero beyond 2 kpc. We also measure the mass fraction in young stars for individual galaxies in the highest S/N central regions. All the galaxies have young components of less than one percent. Our results clearly suggest that the star formation in massive central cluster galaxies takes place in their galaxy cores (<2 kpc), which, with deeper gravitational potential wells, are capable of retaining more gas. Among the possible sources for the gas required to form these young stars, our results are consistent with an in-situ origin via stellar evolution, which is sufficient to produce the observed young stellar populations.


title: Efficient Sampling for Predictor-Based Neural Architecture Search

abstract: Recently, predictor-based algorithms emerged as a promising approach for neural architecture search (NAS). For NAS, we typically have to calculate the validation accuracy of a large number of Deep Neural Networks (DNNs), what is computationally complex. Predictor-based NAS algorithms address this problem. They train a proxy model that can infer the validation accuracy of DNNs directly from their network structure. During optimization, the proxy can be used to narrow down the number of architectures for which the true validation accuracy must be computed, what makes predictor-based algorithms sample efficient. Usually, we compute the proxy for all DNNs in the network search space and pick those that maximize the proxy as candidates for optimization. However, that is intractable in practice, because the search spaces are often very large and contain billions of network architectures. The contributions of this paper are threefold: 1) We define a sample efficiency gain to compare different predictor-based NAS algorithms. 2) We conduct experiments on the NASBench-101 dataset and show that the sample efficiency of predictor-based algorithms decreases dramatically if the proxy is only computed for a subset of the search space. 3) We show that if we choose the subset of the search space on which the proxy is evaluated in a smart way, the sample efficiency of the original predictor-based algorithm that has access to the full search space can be regained. This is an important step to make predictor-based NAS algorithms useful, in practice.


title: Structure learning for zero-inflated counts, with an application to single-cell RNA sequencing data

abstract: The problem of estimating the structure of a graph from observed data is of growing interest in the context of high-throughput genomic data, and single-cell RNA sequencing in particular. These, however, are challenging applications, since the data consist of high-dimensional counts with high variance and over-abundance of zeros. Here, we present a general framework for learning the structure of a graph from single-cell RNA-seq data, based on the zero-inflated negative binomial distribution. We demonstrate with simulations that our approach is able to retrieve the structure of a graph in a variety of settings and we show the utility of the approach on real data.


title: On equivalence of two non-Riemannian curvatures in warped product Finsler metrics

abstract: In this paper we study warped product Finsler metrics and show that the notion of isotropic $E$-curvature and isotropic $S$-curvature are equivalent for this class of metrics.


title: Effective and Sparse Count-Sketch via k-means clustering

abstract: Count-sketch is a popular matrix sketching algorithm that can produce a sketch of an input data matrix X in O(nnz(X))time where nnz(X) denotes the number of non-zero entries in X. The sketched matrix will be much smaller than X while preserving most of its properties. Therefore, count-sketch is widely used for addressing high-dimensionality challenge in machine learning. However, there are two main limitations of count-sketch: (1) The sketching matrix used count-sketch is generated randomly which does not consider any intrinsic data properties of X. This data-oblivious matrix sketching method could produce a bad sketched matrix which will result in low accuracy for subsequent machine learning tasks (e.g.classification); (2) For highly sparse input data, count-sketch could produce a dense sketched data matrix. This dense sketch matrix could make the subsequent machine learning tasks more computationally expensive than on the original sparse data X. To address these two limitations, we first show an interesting connection between count-sketch and k-means clustering by analyzing the reconstruction error of the count-sketch method. Based on our analysis, we propose to reduce the reconstruction error of count-sketch by using k-means clustering algorithm to obtain the low-dimensional sketched matrix. In addition, we propose to solve k-mean clustering using gradient descent with -L1 ball projection to produce a sparse sketched matrix. Our experimental results based on six real-life classification datasets have demonstrated that our proposed method achieves higher accuracy than the original count-sketch and other popular matrix sketching algorithms. Our results also demonstrate that our method produces a sparser sketched data matrix than other methods and therefore the prediction cost of our method will be smaller than other matrix sketching methods.


title: Minimal Virtual Machines on IoT Microcontrollers: The Case of Berkeley Packet Filters with rBPF

abstract: Virtual machines (VM) are widely used to host and isolate software modules. However, extremely small memory and low-energy budgets have so far prevented wide use of VMs on typical microcontroller-based IoT devices. In this paper, we explore the potential of two minimal VM approaches on such low-power hardware. We design rBPF, a register-based VM based on extended Berkeley Packet Filters (eBPF). We compare it with a stack-based VM based on WebAssembly (Wasm) adapted for embedded systems. We implement prototypes of each VM, hosted in the IoT operating system RIOT. We perform measurements on commercial off-the-shelf IoT hardware. Unsurprisingly, we observe that both Wasm and rBPF virtual machines yield execution time and memory overhead, compared to not using a VM. We show however that this execution time overhead is tolerable for low-throughput, low-energy IoT devices. We further show that, while using a VM based on Wasm entails doubling the memory budget for a simple networked IoT application using a 6LoWPAN/CoAP stack, using a VM based on rBPF requires only negligible memory overhead (less than 10% more memory). rBPF is thus a promising approach to host small software modules, isolated from OS software, and updatable on-demand, over low-power networks.


title: Decaying positive global solutions of second order difference equations with mean curvature operator

abstract: A boundary value problem on an unbounded domain, associated to difference equations with the Euclidean mean curvature operator is considered. The existence of solutions which are positive on the whole domain and decaying at infinity is examined by proving new Sturm comparison theorems for linear difference equations and using a fixed point approach based on a linearization device. %The process from the continuous problem to discrete one is examined, too. The process of discretization of the boundary value problem on the unbounded domain is examined, and some discrepancies between the discrete and the continuous case are pointed out, too.


title: Non-Invertible-Element Constacyclic Codes over Finite PIRs

abstract: In this paper we introduce the notion of $\lambda$-constacyclic codes over finite rings $R$ for arbitary element $\lambda$ of $R$. We study the non-invertible-element constacyclic codes (NIE-constacyclic codes) over finite principal ideal rings (PIRs). We determine the algebraic structures of all NIE-constacyclic codes over finite chain rings, give the unique form of the sets of the defining polynomials and obtain their minimum Hamming distances. A general form of the duals of NIE-constacyclic codes over finite chain rings is also provided. In particular, we give a necessary and sufficient condition for the dual of an NIE-constacyclic code to be an NIE-constacyclic code. Using the Chinese Remainder Theorem, we study the NIE-constacyclic codes over finite PIRs. Furthermore, we construct some optimal NIE-constacyclic codes over finite PIRs in the sense that they achieve the maximum possible minimum Hamming distances for some given lengths and cardinalities.


title: Role of inhibitory neurons in temporal correlations of critical and supercritical spontaneous activity

abstract: Experimental and numerical results suggest that the brain can be viewed as a system acting close to a critical point, as confirmed by scale-free distributions of relevant quantities in a variety of different systems and models. Less attention has received the investigation of the temporal correlation functions in brain activity in different, healthy and pathological, conditions. Here we perform this analysis by means of a model with short and long-term plasticity which implements the novel feature of different recovery rates for excitatory and inhibitory neurons, found experimentally. We evidence the important role played by inhibitory neurons in the supercritical state: We detect an unexpected oscillatory behaviour of the correlation decay, whose frequency depends on the fraction of inhibitory neurons and their connectivity degree. This behaviour can be rationalized by the observation that bursts in activity become more frequent and with a smaller amplitude as inhibition becomes more relevant.


title: Red Supergiants in M31 and M33 I. The Complete Sample

abstract: The aim of this paper is to establish a complete sample of red supergiants (RSGs) in M31 and M33. The member stars of the two galaxies are selected from the near-infrared (NIR) point sources after removing the foreground dwarfs from their obvious branch in the $J-H/H-K$ diagram with the archival photometric data taken by the UKIRT/WFCAM. This separation by NIR colors of dwarfs from giants is confirmed by the optical/infrared color-color diagrams ($r-z/z-H$ and $B-V/V-R$), and the Gaia measurement of parallax and proper motion. The RSGs are then identified by their outstanding location in the members' $J-K/K$ diagram due to high luminosity and low effective temperature. The resultant sample has 5,498 and 3,055 RSGs in M31 and M33 respectively, which should be complete because the lower limiting $K$ magnitude of RSGs in both cases is brighter than the complete magnitude of the UKIRT photometry. Analysis of the control fields finds that the pollution rate in the RSGs sample is less than 1\%. The by-product is the complete sample of oxygen-rich asymptotic giant branch stars (AGBs), carbon-rich AGBs, thermally pulsing AGBs and extreme AGBs. In addition, the tip-RGB is determined together with its implication on the distance modulus to M31 and M33.


title: A decentralized approach towards secure firmware updates and testing over commercial IoT Devices

abstract: Internet technologies have made a paradigm shift in the fields of computing and data science and one such paradigm defining change is the Internet of Things or IoT. Nowadays, thousands of household appliances use integrated smart devices which allow remote monitoring and control and also allow intensive computational work such as high end AI-integrated smart security systems with sustained alerts for the user. The update process of these IoT devices usually lack the ability of checking the security of centralized servers, which may be compromised and host malicious firmware files as it is presumed that the servers are secure during deployment. The solution for this problem can be solved using a decentralized database to hold the hashes and the firmware. This paper discusses the possible implications of insecure servers used to host the firmwares of commercial IoT products, and aims to provide a blockchain based decentralized solution to host firmware files with the property of immutability, and controlled access to the firmware upload functions so as to stop unauthorized use. The paper sheds light over possible hardware implementations and the use of cryptographically secure components in such secure architecture models.


title: Intervention fatigue is the primary cause of strong secondary waves in the COVID-19 pandemic

abstract: As of November 2020, the number of COVID-19 cases is increasing rapidly in many countries. In Europe, the virus spread slowed considerably in the late spring due to strict lockdown, but a second wave of the pandemic grew throughout the fall. In this study, we first reconstruct the time evolution of the effective reproduction numbers ${\cal R}(t)$ for each country by integrating the equations of the classic SIR model. We cluster countries based on the estimated ${\cal R}(t)$ through a suitable time series dissimilarity. The result suggests that simple dynamical mechanisms determine how countries respond to changes in COVID-19 case counts. Inspired by these results, we extend the SIR model to include a social response to explain the number $X(t)$ of new confirmed daily cases. As a first-order model, we assume that the social response is on the form $d_t{\cal R}=-\nu (X-X^*)$, where $X^*$ is a threshold for response. The response rate $\nu$ depends on whether $X^*$ is below or above this threshold, on three parameters $\nu_1,\;\nu_2,\,\nu_3,$, and on $t$. When $X<X^*$, $\nu= \nu_1$, describes the effect of relaxed intervention when the incidence rate is low. When $X>X^*$, $\nu=\nu_2\exp{(-\nu_3t)}$, models the impact of interventions when incidence rate is high. The parameter $\nu_3$ represents the fatigue, i.e., the reduced effect of intervention as time passes. The proposed model reproduces typical evolving patterns of COVID-19 epidemic waves observed in many countries. Estimating the parameters $\nu_1,\,\nu_2,\,\nu_3$ and initial conditions, such as ${\cal R}_0$, for different countries helps to identify important dynamics in their social responses. One conclusion is that the leading cause of the strong second wave in Europe in the fall of 2020 was not the relaxation of interventions during the summer, but rather the general fatigue to interventions developing in the fall.


title: Derived Langlands VI: Monomial resolutions and $2$-variable L-functions

abstract: In this brief essay a construction of the $2$-variable L-function of Langlands is sketched in terms of monomial resolutions of admissible representations of reductive locally $p$-adic Lie groups.


title: Dzyaloshinskii-Moriya interaction induced by an ultrashort electromagnetic pulse: Application to coherent (anti)ferromagnetic skyrmion nucleation

abstract: We show how a Dzyaloshinskii-Moriya interaction can be generated in an ultrathin metal film from the induced internal electric field created by an ultrashort electromagnetic pulse. This interaction does not require structural inversion-symmetry breaking, and its amplitude can be tuned depending on the amplitude of the field. We perform first-principles calculations to estimate the strength of the field-induced magnetoelectric coupling for ferromagnetic Fe, Co, and Ni, and antiferromagnetic Mn, as well as FePt alloys. Last, using atomistic simulations, we demonstrate how an isolated antiferromagnetic skyrmion can be coherently nucleated from the collinear background by an ultrashort pulse in electric field on a 100-fs timescale.


title: Uniaxial N\'eel Vector Control in Perovskite Oxide Thin Films by Anisotropic Strain Engineering

abstract: Antiferromagnetic thin films typically exhibit a multi-domain state, and control of the antiferromagnetic N\'eel vector is challenging as antiferromagnetic materials are robust to magnetic perturbations. By relying on anisotropic in-plane strain engineering of epitaxial thin films of the prototypical antiferromagnetic material LaFeO3, uniaxial N\'eel vector control is demonstrated. Orthorhombic (011)- and (101)-oriented DyScO3, GdScO3 and NdGaO3 substrates are used to engineer different anisotropic in-plane strain states. The anisotropic in-plane strain stabilises structurally monodomain monoclinic LaFeO3 thin films. The uniaxial N\'eel vector is found along the tensile strained b axis, contrary to bulk LaFeO3 having the N\'eel vector along the shorter a axis, and no magnetic domains are found. Hence, anisotropic strain engineering is a viable tool for designing unique functional responses, further enabling antiferromagnetic materials for mesoscopic device technology.


title: Using Machine Learning to Create an Early Warning System for Welfare Recipients

abstract: Using high-quality nation-wide social security data combined with machine learning tools, we develop predictive models of income support receipt intensities for any payment enrolee in the Australian social security system between 2014 and 2018. We show that off-the-shelf machine learning algorithms can significantly improve predictive accuracy compared to simpler heuristic models or early warning systems currently in use. Specifically, the former predicts the proportion of time individuals are on income support in the subsequent four years with greater accuracy, by a magnitude of at least 22% (14 percentage points increase in the R2), compared to the latter. This gain can be achieved at no extra cost to practitioners since the algorithms use administrative data currently available to caseworkers. Consequently, our machine learning algorithms can improve the detection of long-term income support recipients, which can potentially provide governments with large savings in accrued welfare costs.


title: Growth of Perturbations using Lambert$W$ Equation of State

abstract: Recently, a novel equation of state (EoS) parameter for dark energy has been introduced which deals with a special mathematical function, known as the Lambert$W$ function. In this paper, we study the effect on the growth of perturbations for the Lambert$W$ dark energy model. We perform the analysis for two different approaches. In the first case we consider the universe to be filled with two different fluid components, namely, the baryonic matter component and the Lambert$W$ dark energy component, while in the second case we consider that there is a single fluid component in the universe whose equation of state parameter is described by the Lambert$W$ function. We then compare the growth rates of Lambert$W$ model with that for a standard $\Lambda$CDM model as well as the CPL model. Our results indicate that the presence of Lambert$W$ dynamical dark energy sector changes the growth rate and affects the matter fluctuations in the universe to a great extent.


title: Infrared small target detection based on isotropic constraint under complex background

abstract: Infrared search and tracking (IRST) system has been widely concerned and applied in the area of national defence. Small target detection under complex background is a very challenging task in the development of system algorithm. Low signal-to-clutter ratio (SCR) of target and the interference caused by irregular background clutter make it difficult to get an accurate result. In this paper, small targets are considered to have two characteristics of high contrast and isotropy, and we propose a multilayer gray difference (MGD) method constrained by isotropy. Firstly, the suspected regions are obtained through MGD, and then the eigenvalues of the original image's Hessian matrix are calculated to obtain the isotropy parameter of each region. Finally, those regions do not meet the isotropic constraint condition are suppressed. Experiments show that the proposed method is effective and superior to several common methods in terms of signal-to-clutter ratio gain (SCRG) and receiver operating characteristic (ROC) curve.


title: The rotation rate of solar active and ephemeral regions -- I. Dependence on morphology and peak magnetic flux

abstract: Using magnetic field maps acquired by the Helioseismic and Magnetic Imager on board the Solar Dynamics Observatory we measured rotation rates of 864 active and 322 ephemeral regions observed between 2010 and 2016. We found smaller magnetic tracers to show a tendency to rotate faster as compared to larger ones. Thus, ephemeral regions exhibit on average the fastest rotation rate. We further divided active regions into three classes. Class A comprised magnetic bipoles obeying Hale's polarity law, Joy's law, and exhibiting more coherent leading polarity in comparison with the following one. The second class B included active regions violating at least one of the aforementioned empirical laws. The third class U comprised unipolar active regions. We found no significant difference between the rotation rates of active regions of classes A and B. In contrast, unipolar active regions exhibited on average lower rotation rate and narrower distribution of the rotation rate differences. Assuming the rotation rate to indicate the anchoring depth of the magnetic structure within the convection zone, we supposed that active regions of classes A and B might be anchored throughout the entire convective envelope while unipolar active regions a rooted within a thin layer located either near the base of the convection zone or at a shallow near-surface depth.


title: Constant-round Multi-party Quantum Computation for Constant Parties

abstract: One of the central themes in classical cryptography is multi-party computation, which performs joint computation on multiple participants' data while maintaining data privacy. The extension to the quantum regime was proposed in 2002, but despite two decades of research, the current state-of-the-art multi-party quantum computation protocol for a fixed number of parties (even 2 parties) requires unbounded rounds of communication, which greatly limit its practical usage. In this work, we propose the first constant-round multi-party quantum computation protocol for a fixed number of parties against specious adversaries, thereby significantly reducing the required number of rounds. Our work constitutes a key step towards practical implementation of secure multi-party quantum computation, and opens the door for practitioners to be involved in this exciting field. The result of our work has wide implications to quantum fidelity evaluation, quantum machine learning, quantum voting, and distributed quantum information processing.


title: On the possibility to probe the flare productivity of an active region on the early stage of emergence

abstract: A long-standing problem is to predict the future flare productivity of an active region (AR) when it is on the stage of early emergence. The aim of this study is to probe two parameters of the photospheric magnetic field, both derived during the emergence phase of an AR, and to compare them with the flare productivity of the well-developed AR. The parameters are: (i) the index of the magnetic power spectrum (the slope of the spectrum) on the stage of emergence, and (ii) the flux emergence rate. Analysis of 243 emerging ARs showed that the magnetic power index increases from values typical for quiet-Sun regions to the values typical for mature ARs within a day while the emergence proceeds for several days; frequently, after the increase, the value of the power index undulates around some mean value with fluctuations being several time less than the growth of the power index during the emergence onset. For a subset of 34 flare-productive ARs we found no correlation between the power spectrum index at the stage of emergence and flare index derived from the entire interval of the AR's presence on the disc. At the same time, the flux emergence rate correlates well with the flare index (the Pearson's correlation coefficient is 0.74). We conclude that high flux emergence rate is a necessary condition for an AR to produce strong flares in the future, and so, the flux emergence rate can be used to probe the future flare productivity of the AR.


title: How Far Are We from Robust Voice Conversion: A Survey

abstract: Voice conversion technologies have been greatly improved in recent years with the help of deep learning, but their capabilities of producing natural sounding utterances in different conditions remain unclear. In this paper, we gave a thorough study of the robustness of known VC models. We also modified these models, such as the replacement of speaker embeddings, to further improve their performances. We found that the sampling rate and audio duration greatly influence voice conversion. All the VC models suffer from unseen data, but AdaIN-VC is relatively more robust. Also, the speaker embedding jointly trained is more suitable for voice conversion than those trained on speaker identification.


title: A unified diagrammatic approach to topological fixed point models

abstract: We introduce a systematic mathematical language for describing fixed point models and apply it to the study to topological phases of matter. The framework is reminiscent of state-sum models and lattice topological quantum field theories, but is formalised and unified in terms of tensor networks. In contrast to existing tensor network ansatzes for the study of ground states of topologically ordered phases, the tensor networks in our formalism represent discrete path integrals in Euclidean space-time. This language is more directly related to the Hamiltonian defining the model than other approaches, via a Trotterization of the respective imaginary time evolution. We introduce our formalism by simple examples, and demonstrate its full power by expressing known families of models in 2+1 dimensions in their most general form, namely string-net models and Kitaev quantum doubles based on weak Hopf algebras. To elucidate the versatility of our formalism, we also show how fermionic phases of matter can be described and provide a framework for topological fixed point models in 3+1 dimensions.


title: Two types of size Ramsey numbers for matchings of small order

abstract: For simple graphs $G$ and $H$, their size Ramsey number $\hat{r}(G,H)$ is the smallest possible size of $F$ such that for any red-blue coloring of its edges, $F$ contains either a red $G$ or a blue $H$. Similarly, we can define the connected size Ramsey number ${\hat{r}}_c(G,H)$ by adding the prerequisite that $F$ must be connected. In this paper, we explore the relationships between these size Ramsey numbers and give some results on their values for certain classes of graphs. We are mainly interested in the cases where $G$ is either a $2K_2$ or a $3K_2$, and where $H$ is either a cycle $C_n$ or a union of paths $nP_m$. Additionally, we improve an upper bound regarding the values of $\hat{r}(tK_2,P_m)$ and ${\hat{r}}_c(tK_2,P_m)$ for certain $t$ and $m$.


title: The images of special generic maps of an elementary class

abstract: The class of special generic maps contains Morse functions with exactly two singular points, characterizing spheres topologically which are not $4$-dimensional and the $4$-dimensional unit sphere. This class is for higher dimensional versions of such functions. Canonical projections of unit spheres are simplest examples and suitable manifolds diffeomorphic to ones represented as connected sums of products of spheres admit such maps. They have been found to restrict the topologies and the differentiable structures of the manifolds strongly. The present paper focuses on images of special generic maps on closed manifolds. They are smoothly immersed compact manifolds whose dimensions are same as those of the targets. Some studies imply that they have much information on homology groups and cohomology rings. We present new construction and explicit examples of special generic maps by investigating the images and the present paper is essentially on construction and explicit algebraic topological and differential topological studies of immersed compact manifolds.


title: Dendritic trafficking: synaptic scaling and structural plasticity

abstract: Neuronal circuits internally regulate electrical signaling via a host of homeostatic mechanisms. Two prominent mechanisms, synaptic scaling and structural plasticity, are believed to maintain average activity within an operating range by modifying the strength and spatial extent of network connectivity using negative feedback. However, both mechanisms operate on relatively slow timescales and thus face fundamental limits due to delays. We show that these mechanisms fulfill complementary roles in maintaining stability in a large network. In particular, even relatively, slow growth dynamics improves performance significantly beyond synaptic scaling alone.


title: A phase field formulation for dissolution-driven stress corrosion cracking

abstract: We present a new theoretical and numerical framework for modelling mechanically-assisted corrosion in elastic-plastic solids. Both pitting and stress corrosion cracking (SCC) can be captured, as well as the pit-to-crack transition. Localised corrosion is assumed to be dissolution-driven and a formulation grounded upon the film rupture-dissolution-repassivation mechanism is presented to incorporate the influence of film passivation. The model incorporates, for the first time, the role of mechanical straining as the electrochemical driving force, accelerating corrosion kinetics. The computational complexities associated with tracking the evolving metal-electrolyte interface are resolved by making use of a phase field paradigm, enabling an accurate approximation of complex SCC morphologies. The coupled electro-chemo-mechanical formulation is numerically implemented using the finite element method and an implicit time integration scheme; displacements, phase field order parameter and concentration are the primary variables. Five case studies of particular interest are addressed to showcase the predictive capabilities of the model, revealing an excellent agreement with analytical solutions and experimental measurements. By modelling these paradigmatic 2D and 3D boundary value problems we show that our formulation can capture: (i) the transition from activation-controlled corrosion to diffusion-controlled corrosion, (ii) the sensitivity of interface kinetics to mechanical stresses and strains, (iii) the role of film passivation in reducing corrosion rates, and (iv) the dependence of the stability of the passive film to local strain rates. The influence of these factors in driving the shape change of SCC defects, including the pit-to-crack transition, is a natural outcome of the model, laying the foundations for a mechanistic assessment of engineering materials and structures.


title: Towards SAR Tomographic Inversion via Sparse Bayesian Learning

abstract: Existing SAR tomography (TomoSAR) algorithms are mostly based on an inversion of the SAR imaging model, which are often computationally expensive. Previous study showed perspective of using data-driven methods like KPCA to decompose the signal and reduce the computational complexity. This paper gives a preliminary demonstration of a new data-driven method based on sparse Bayesian learning. Experiments on simulated data show that the proposed method significantly outperforms KPCA methods in estimating the steering vectors of the scatterers. This gives a perspective of data-drive approach or combining it with model-driven approach for high precision tomographic inversion of large areas.


title: Astrophysical bounds on mirror dark matter derived from binary pulsars timing data

abstract: Mirror Dark Matter (MDM) has been considered as an elegant framework for a particle theory of Dark Matter (DM). It is supposed that there exists a dark sector which is mirror of the ordinary matter. Some MDM models allow particle interactions mirror and ordinary matter, in addition to the gravitational interaction. The possibility of neutron to mirror neutron transition has recently been discussed both from theoretical and experimental perspectives. This paper is based on a previous work in which we obtained stringent upper limits on the possibility of converting neutrons to mirror neutrons in the interiors of neutron stars, by using timing data of binary pulsars. Such a transition would imply mass loss in neutron stars leading to a significant change of orbital period of neutron star binary systems. The observational bounds on the period changes of such binaries, therefore put strong limits on the above transition rate and hence on the neutron -- mirror-neutron mixing parameter $\epsilon'$. Our limits are much stronger than the values required to explain the neutron decay anomaly via $n-n'$ mixing.


title: FOWD: A Free Ocean Wave Dataset for Data Mining and Machine Learning

abstract: The occurrence of extreme (rogue) waves in the ocean is for the most part still shrouded in mystery, as the rare nature of these events makes them difficult to analyze with traditional methods. Modern data mining and machine learning methods provide a promising way out, but they typically rely on the availability of massive amounts of well-cleaned data. To facilitate the application of such data-hungry methods to surface ocean waves, we developed FOWD, a freely available wave dataset and processing framework. FOWD describes the conversion of raw observations into a catalogue that maps characteristic sea state parameters to observed wave quantities. Specifically, we employ a running window approach that respects the non-stationary nature of the oceans, and extensive quality control to reduce bias in the resulting dataset. We also supply a reference Python implementation of the FOWD processing toolkit, which we use to process the entire CDIP buoy data catalogue containing over 4 billion waves. In a first experiment, we find that, when the full elevation time series is available, surface elevation kurtosis and maximum wave height are the strongest univariate predictors for rogue wave activity. When just a spectrum is given, crest-trough correlation, spectral bandwidth, and mean period fill this role.


title: Flux-density measurements of the high-redshift blazar PSO J047.4478+27.2992 at 4.7 and 8.2 GHz with RATAN-600

abstract: We report the first detection at 4.7 and 8.2 GHz with the RATAN-600 radio telescope of the source PSO J047.4478+27.2992, which has been announced as a blazar at z=6.1. The average flux density is $12\pm3$ and $8\pm4$ mJy at 4.7 and 8.2 GHz respectively, and an upper limit is estimated as 3 mJy at 11.2 GHz. The light curve at 4.7 GHz, obtained on a time-scale of four months, exhibits moderate variability of the source (fractional variability $F_{var}=0.28\pm0.02$). The new RATAN-600 measurements extend previous literature data to higher frequencies, indicating a single power-law radio spectrum with $\alpha_{0.147-8.2}$=$-0.51\pm0.1$. The monochromatic radio luminosity at 4.7 GHz $P_{\nu} \sim 2 \times 10^{27}$ W Hz$^{-1}$ is consistent with values for high-redshift quasars at $z\geq 3$.


title: Picking BERT's Brain: Probing for Linguistic Dependencies in Contextualized Embeddings Using Representational Similarity Analysis

abstract: As the name implies, contextualized representations of language are typically motivated by their ability to encode context. Which aspects of context are captured by such representations? We introduce an approach to address this question using Representational Similarity Analysis (RSA). As case studies, we investigate the degree to which a verb embedding encodes the verb's subject, a pronoun embedding encodes the pronoun's antecedent, and a full-sentence representation encodes the sentence's head word (as determined by a dependency parse). In all cases, we show that BERT's contextualized embeddings reflect the linguistic dependency being studied, and that BERT encodes these dependencies to a greater degree than it encodes less linguistically-salient controls. These results demonstrate the ability of our approach to adjudicate between hypotheses about which aspects of context are encoded in representations of language.


title: An Optimal Triangle Projector with Prescribed Area and Orientation, Application to Position-Based Dynamics

abstract: The vast majority of mesh-based modelling applications iteratively transform the mesh vertices under prescribed geometric conditions. This occurs in particular in methods cycling through the constraint set such as Position-Based Dynamics (PBD). A common case is the approximate local area preservation of triangular 2D meshes under external editing constraints. At the constraint level, this yields the nonconvex optimal triangle projection under prescribed area problem, for which there does not currently exist a direct solution method. In current PBD implementations, the area preservation constraint is linearised. The solution comes out through the iterations, without a guarantee of optimality, and the process may fail for degenerate inputs where the vertices are colinear or colocated. We propose a closed-form solution method and its numerically robust algebraic implementation. Our method handles degenerate inputs through a two-case analysis of the problem's generic ambiguities. We show in a series of experiments in area-based 2D mesh editing that using optimal projection in place of area constraint linearisation in PBD speeds up and stabilises convergence.


title: Fuzzy Stochastic Timed Petri Nets for Causal properties representation

abstract: Imagery is frequently used to model, represent and communicate knowledge. In particular, graphs are one of the most powerful tools, being able to represent relations between objects. Causal relations are frequently represented by directed graphs, with nodes denoting causes and links denoting causal influence. A causal graph is a skeletal picture, showing causal associations and impact between entities. Common methods used for graphically representing causal scenarios are neurons, truth tables, causal Bayesian networks, cognitive maps and Petri Nets. Causality is often defined in terms of precedence (the cause precedes the effect), concurrency (often, an effect is provoked simultaneously by two or more causes), circularity (a cause provokes the effect and the effect reinforces the cause) and imprecision (the presence of the cause favors the effect, but not necessarily causes it). We will show that, even though the traditional graphical models are able to represent separately some of the properties aforementioned, they fail trying to illustrate indistinctly all of them. To approach that gap, we will introduce Fuzzy Stochastic Timed Petri Nets as a graphical tool able to represent time, co-occurrence, looping and imprecision in causal flow.


title: Sharp time decay estimates for the discrete Klein-Gordon equation

abstract: We establish sharp time decay estimates for the the Klein-Gordon equation on the cubic lattice in dimensions $d=2,3,4$. The $\ell^1\to\ell^{\infty}$ dispersive decay rate is $|t|^{-3/4}$ for $d=2$, $|t|^{-7/6}$ for $d=3$ and $|t|^{-3/2}\log|t|$ for $d=4$. These decay rates are faster than conjectured by Kevrekidis and Stefanov (2005). The proof relies on oscillatory integral estimates and proceeds by a detailed analysis of the the singularities of the associated phase function. We also prove new Strichartz estimates and discuss applications to nonlinear PDEs and spectral theory.


title: CLAWS: Clustering Assisted Weakly Supervised Learning with Normalcy Suppression for Anomalous Event Detection

abstract: Learning to detect real-world anomalous events through video-level labels is a challenging task due to the rare occurrence of anomalies as well as noise in the labels. In this work, we propose a weakly supervised anomaly detection method which has manifold contributions including1) a random batch based training procedure to reduce inter-batch correlation, 2) a normalcy suppression mechanism to minimize anomaly scores of the normal regions of a video by taking into account the overall information available in one training batch, and 3) a clustering distance based loss to contribute towards mitigating the label noise and to produce better anomaly representations by encouraging our model to generate distinct normal and anomalous clusters. The proposed method obtains83.03% and 89.67% frame-level AUC performance on the UCF Crime and ShanghaiTech datasets respectively, demonstrating its superiority over the existing state-of-the-art algorithms.


title: Continuous functions on Berkovich spaces

abstract: Let $k$ be a perfect complete valued field with a nontrivial non-archimedean norm $|\cdot|$ and $\omega\in k$ with $0<|\omega|<1.$ Let $X$ be a reduced and normal $k$-analytic space. Then $O^{\circ}\simeq \lim\limits_{\leftarrow}O^{\circ}/\omega^n.$


abstract: To eliminate the problems of large dimensional differences, big semantic gap, and mutual interference caused by hybrid features, in this paper, we propose a novel Multi-Features Guidance Network for partial-to-partial point cloud registration(MFG). The proposed network mainly includes four parts: keypoints' feature extraction, correspondences searching, correspondences credibility computation, and SVD, among which correspondences searching and correspondence credibility computation are the cores of the network. Unlike the previous work, we utilize the shape features and the spatial coordinates to guide correspondences search independently and fusing the matching results to obtain the final matching matrix. In the correspondences credibility computation module, based on the conflicted relationship between the features matching matrix and the coordinates matching matrix, we score the reliability for each correspondence, which can reduce the impact of mismatched or non-matched points. Experimental results show that our network outperforms the current state-of-the-art while maintaining computational efficiency.


title: The 2-Adic Complexity of Two Classes of Binary Sequences with Interleaved Structure

abstract: The autocorrelation values of two classes of binary sequences are shown to be good in [6]. We study the 2-adic complexity of these sequences. Our results show that the 2-adic complexity of such sequences is large enough to resist the attack of the rational approximation algorithm.


title: Tackling Domain-Specific Winograd Schemas with Knowledge-Based Reasoning and Machine Learning

abstract: The Winograd Schema Challenge (WSC) is a common-sense reasoning task that requires background knowledge. In this paper, we contribute to tackling WSC in four ways. Firstly, we suggest a keyword method to define a restricted domain where distinctive high-level semantic patterns can be found. A thanking domain was defined by key-words, and the data set in this domain is used in our experiments. Secondly, we develop a high-level knowledge-based reasoning method using semantic roles which is based on the method of Sharma [2019]. Thirdly, we propose an ensemble method to combine knowledge-based reasoning and machine learning which shows the best performance in our experiments. As a machine learning method, we used Bidirectional Encoder Representations from Transformers (BERT) [Kocijan et al., 2019]. Lastly, in terms of evaluation, we suggest a "robust" accuracy measurement by modifying that of Trichelair et al. [2018]. As with their switching method, we evaluate a model by considering its performance on trivial variants of each sentence in the test set.


title: Computational efficient deep neural network with difference attention maps for facial action unit detection

abstract: In this paper, we propose a computational efficient end-to-end training deep neural network (CEDNN) model and spatial attention maps based on difference images. Firstly, the difference image is generated by image processing. Then five binary images of difference images are obtained using different thresholds, which are used as spatial attention maps. We use group convolution to reduce model complexity. Skip connection and $\text{1}\times \text{1}$ convolution are used to ensure good performance even if the network model is not deep. As an input, spatial attention map can be selectively fed into the input of each block. The feature maps tend to focus on the parts that are related to the target task better. In addition, we only need to adjust the parameters of classifier to train different numbers of AU. It can be easily extended to varying datasets without increasing too much computation. A large number of experimental results show that the proposed CEDNN is obviously better than the traditional deep learning method on DISFA+ and CK+ datasets. After adding spatial attention maps, the result is better than the most advanced AU detection method. At the same time, the scale of the network is small, the running speed is fast, and the requirement for experimental equipment is low.


title: The antiferromagnetic $S=1/2$ Heisenberg model on the ${\text{C}}_{60}$ fullerene geometry

abstract: We solve the quantum-mechanical antiferromagnetic Heisenberg model with spins positioned on vertices of the truncated icosahedron using the density-matrix renormalization group (DMRG). This describes magnetic properties of the undoped C$_{60}$ fullerene at half filling in the limit of strong on-site interaction $U$. We calculate the ground state and correlation functions for all possible distances, the lowest singlet and triplet excited states, as well as thermodynamic properties, namely the specific heat and spin susceptibility. We find that unlike the exactly solvable C$_{20}$ to C$_{32}$, the lowest excited state is a triplet rather than a singlet, indicating a reduced frustration due to the presence of many hexagon faces and the separation of the pentagon faces. This implies that frustration may be tuneable within the fullerenes by changing their size. The spin-spin correlations are much stronger along the hexagon bonds and rapidly decrease with distance, so that the molecule is large enough not to be correlated across its whole extent. The specific heat shows a high-temperature peak and a low-temperature shoulder reminiscent of the Kagom\'e lattice, while the spin susceptibility shows a single broad peak and is very close to the one of C$_{20}$.


title: ($q$-)Supercongruences hit again

abstract: Using an intrinsic $q$-hypergeometric strategy, we generalise Dwork-type congruences $H(p^{s+1})/H(p^s)\equiv H(p^s)/H(p^{s-1})\pmod{p^3}$ for $s=1,2,\dots$ and $p$ a prime, when $H(N)$ are truncated hypergeometric sums corresponding to the periods of rigid Calabi--Yau threefolds.

categories: math.NT math-ph math.CA math.CO math.MP math.QA


title: On stability of nonzero set-point for non linear impulsive control systems

abstract: The interest in non-linear impulsive systems (NIS) has been growing due to its impact in application problems such as disease treatments (diabetes, HIV, influenza, among many others), where the control action (drug administration) is given by short-duration pulses followed by time periods of null values. Within this framework the concept of equilibrium needs to be extended (redefined) to allows the system to keep orbiting (between two consecutive pulses) in some state space regions out of the origin, according to usual objectives of most real applications. Although such regions can be characterized by means of a discrete-time system obtained by sampling the NIS at the impulsive times, no agreements have reached about their asymptotic stability (AS). This paper studies the asymptotic stability of control equilibrium orbits for NSI, based on the underlying discrete time system, in order to establish the conditions under which the AS for the latter leads to the AS for the former. Furthermore, based on the latter AS characterization, an impulsive Model Predictive Control (i-MPC) that feasibly stabilizes the non-linear impulsive system is presented. Finally, the proposed stable MPC is applied to two control problems of interest: the intravenous bolus administration of Lithium and the administration of antiretrovirals for HIV treatments.


title: Unequal Representations: Analyzing Intersectional Biases in Word Embeddings Using Representational Similarity Analysis

abstract: We present a new approach for detecting human-like social biases in word embeddings using representational similarity analysis. Specifically, we probe contextualized and non-contextualized embeddings for evidence of intersectional biases against Black women. We show that these embeddings represent Black women as simultaneously less feminine than White women, and less Black than Black men. This finding aligns with intersectionality theory, which argues that multiple identity categories (such as race or sex) layer on top of each other in order to create unique modes of discrimination that are not shared by any individual category.


title: A Convenient Infinite Dimensional Framework for Generative Adversarial Learning

abstract: In recent years, generative adversarial networks (GANs) have demonstrated impressive experimental results while there are only a few works that foster statistical learning theory for GANs. In this work, we propose an infinite dimensional theoretical framework for generative adversarial learning. We assume that the probability density functions of the underlying measure are uniformly bounded, $k$-times $\alpha$-H\"{o}lder differentiable ($C^{k,\alpha}$) and uniformly bounded away from zero. Under these assumptions, we show that the Rosenblatt transformation induces an optimal generator, which is realizable in the hypothesis space of $C^{k,\alpha}$-generators. With a consistent definition of the hypothesis space of discriminators, we further show that the Jensen-Shannon divergence between the distribution induced by the generator from the adversarial learning procedure and the data generating distribution converges to zero. Under certain regularity assumptions on the density of the data generating process, we also provide rates of convergence based on chaining and concentration.


title: M\"obius Strip Microlasers: a Testbed for Non-Euclidean Photonics

abstract: We report on experiments with M\"obius strip microlasers which were fabricated with high optical quality by direct laser writing. A M\"obius strip, i.e., a band with a half twist, exhibits the fascinating property that it has a single nonorientable surface and a single boundary. We provide evidence that, in contrast to conventional ring or disk resonators, a M\"obius strip cavity cannot sustain whispering gallery modes (WGM). Comparison between experiments and 3D finite difference time domain (FDTD) simulations reveals that the resonances are localized on periodic geodesics.


title: Normal CM-fields with class number one

abstract: We show that assuming the generalized Riemann hypothesis there are no normal CM-fields with class number one of degree 64 and 96. This is done by constructing complete tables of normal CM-fields using discriminant bounds of Lee--Kwon. This solves the class number one problem for normal CM-fields assuming GRH. Using the same technique to solve the relative class number one problem in degrees 16, 32, 56 and 82, also the corresponding relative class number one problem is solved assuming GRH.


title: AI Discovering a Coordinate System of Chemical Elements: Dual Representation by Variational Autoencoders

abstract: The periodic table is a fundamental representation of chemical elements that plays essential theoretical and practical roles. The research article discusses the experiences of unsupervised training of neural networks to represent elements on the 2D latent space based on their electron configurations. To emphasize chemical properties of the elements, the original data of electron configurations has been realigned towards valence orbitals. Recognizing seven shells and four subshells, the input data has been arranged as 7x4 images. Latent space representation has been performed using a convolutional beta variational autoencoder (beta-VAE). Despite discrete and sparse input data, the beta-VAE disentangles elements of different periods, blocks, groups, and types. The unsupervised representation of elements on the latent space reveals pairwise symmetries of periods and elements related to the invariance of quantum numbers of corresponding elements. In addition, it isolates outliers that turned out to be known cases of Madelung's rule violations for lanthanide and actinide elements. Considering the generative capabilities of beta-VAE, the supervised machine learning has been set to find out if there are insightful patterns distinguishing electron configurations between real elements and decoded artificial ones. Also, the article addresses the capability of dual representation by autoencoders. Conventionally, autoencoders represent observations of input data on the latent space. By transposing and duplicating original input data, it is possible to represent variables on the latent space which can lead to the discovery of meaningful patterns among input variables. Applying that unsupervised learning for transposed data of electron configurations, the order of input variables that has been arranged by the encoder on the latent space has turned out to exactly match the sequence of Madelung's rule.


title: SEA: Sentence Encoder Assembly for Video Retrieval by Textual Queries

abstract: Retrieving unlabeled videos by textual queries, known as Ad-hoc Video Search (AVS), is a core theme in multimedia data management and retrieval. The success of AVS counts on cross-modal representation learning that encodes both query sentences and videos into common spaces for semantic similarity computation. Inspired by the initial success of previously few works in combining multiple sentence encoders, this paper takes a step forward by developing a new and general method for effectively exploiting diverse sentence encoders. The novelty of the proposed method, which we term Sentence Encoder Assembly (SEA), is two-fold. First, different from prior art that use only a single common space, SEA supports text-video matching in multiple encoder-specific common spaces. Such a property prevents the matching from being dominated by a specific encoder that produces an encoding vector much longer than other encoders. Second, in order to explore complementarities among the individual common spaces, we propose multi-space multi-loss learning. As extensive experiments on four benchmarks (MSR-VTT, TRECVID AVS 2016-2019, TGIF and MSVD) show, SEA surpasses the state-of-the-art. In addition, SEA is extremely ease to implement. All this makes SEA an appealing solution for AVS and promising for continuously advancing the task by harvesting new sentence encoders.


title: Leveraging Architectural Support of Three Page Sizes with Trident

abstract: Large pages are commonly deployed to reduce address translation overheads for big-memory workloads. Modern x86-64 processors from Intel and AMD support two large page sizes -- 1GB and 2MB. However, previous works on large pages have primarily focused on 2MB pages, partly due to lack of substantial evidence on the profitability of 1GB pages to real-world applications. We argue that in fact, inadequate system software support is responsible for a decade of underutilized hardware support for 1GB pages. Through extensive experimentation on a real system, we demonstrate that 1GB pages can improve performance over 2MB pages, and when used in tandem with 2MB pages for an important set of applications; the support for the latter is crucial but missing in current systems. Our design and implementation of \trident{} in Linux fully exploit hardware supported large pages by dynamically and transparently allocating 1GB, 2MB, and 4KB pages as deemed suitable. \trident{} speeds up eight memory-intensive applications by {$18\%$}, on average, over Linux's use of 2MB pages. We also propose \tridentpv{}, an extension to \trident{} that effectively virtualizes 1GB pages via copy-less promotion and compaction in the guest OS. Overall, this paper shows that even GB-sized pages have considerable practical significance with adequate software enablement, in turn motivating architects to continue investing/innovating in large pages.


title: Smoothing does not give a selection principle for transport equations with bounded autonomous fields

abstract: We give an example of a bounded divergence free autonomous vector field in $\mathbb R^3$ (and of a nonautonomous bounded divergence free vector field in $\mathbb R^2$) and of a bounded initial data for which the Cauchy problem for the corresponding transport equation has $2$ distinct solutions. We then show that both solutions are limits of classical solutions of transport equations for appropriate smoothings of the vector fields and of the initial data.


title: Stochastic effects of waves on currents in the ocean mixed layer

abstract: This paper introduces an energy-preserving stochastic model for studying wave effects on currents in the ocean mixing layer. The model is called stochastic forcing by Lie transport (SFLT). The SFLT model is derived here from a stochastic constrained variational principle, so it has a Kelvin circulation theorem. The examples of SFLT given here treat 3D Euler fluid flow, rotating shallow water dynamics and the Euler-Boussinesq equations. In each example, one sees the effect of stochastic Stokes drift and material entrainment in the generation of fluid circulation. We also present an Eulerian-averaged SFLT model (EA SFLT), based on decomposing the Eulerian solutions of the energy-conserving SFLT model into sums of their expectations and fluctuations.


title: CCIC-WSN: An Architecture for Single Channel Cluster-based Information-Centric Wireless Sensor Networks

abstract: The promising vision of Information-Centric Networking (ICN) and of its realization, Named Data Networking (NDN), has attracted extensive attention in recent years in the context of the Internet of Things (IoT) and Wireless Sensor Networks (WSNs). However, a comprehensive NDN/ICN-based architectural design for WSNs, including specially tailored naming schemes and forwarding mechanisms, has yet to be explored. In this paper, we present single-Channel Cluster-based Information-Centric WSN (CCIC-WSN), an NDN/ICN-based framework to fulfill the requirements of cluster-based WSNs, such as communication between child nodes and cluster heads, association of new child nodes with cluster heads, discovery of the namespace of newly associated nodes, and child node mobility. Through an extensive simulation study, we demonstrate that CCIC-WSN achieves 71-90% lower energy consumption and 74-96% lower data retrieval delays than recently proposed frameworks for NDN/ICN-based WSNs under various evaluation settings.


title: Gender bias in magazines oriented to men and women: a computational approach

abstract: Cultural products are a source to acquire individual values and behaviours. Therefore, the differences in the content of the magazines aimed specifically at women or men are a means to create and reproduce gender stereotypes. In this study, we compare the content of a women-oriented magazine with that of a men-oriented one, both produced by the same editorial group, over a decade (2008-2018). With Topic Modelling techniques we identify the main themes discussed in the magazines and quantify how much the presence of these topics differs between magazines over time. Then, we performed a word-frequency analysis to validate this methodology and extend the analysis to other subjects that did not emerge automatically. Our results show that the frequency of appearance of the topics Family, Business and Women as sex objects, present an initial bias that tends to disappear over time. Conversely, in Fashion and Science topics, the initial differences between both magazines are maintained. Besides, we show that in 2012, the content associated with horoscope increased in the women-oriented magazine, generating a new gap that remained open over time. Also, we show a strong increase in the use of words associated with feminism since 2015 and specifically the word abortion in 2018. Overall, these computational tools allowed us to analyse more than 24,000 articles. Up to our knowledge, this is the first study to compare magazines in such a large dataset, a task that would have been prohibitive using manual content analysis methodologies.


title: Learning to Sample the Most Useful Training Patches from Images

abstract: Some image restoration tasks like demosaicing require difficult training samples to learn effective models. Existing methods attempt to address this data training problem by manually collecting a new training dataset that contains adequate hard samples, however, there are still hard and simple areas even within one single image. In this paper, we present a data-driven approach called PatchNet that learns to select the most useful patches from an image to construct a new training set instead of manual or random selection. We show that our simple idea automatically selects informative samples out from a large-scale dataset, leading to a surprising 2.35dB generalisation gain in terms of PSNR. In addition to its remarkable effectiveness, PatchNet is also resource-friendly as it is applied only during training and therefore does not require any additional computational cost during inference.


title: A robust DPG method for large domains

abstract: We observe a dramatic lack of robustness of the DPG method when solving problems on large domains and where stability is based on a Poincar\'e-type inequality. We show how robustness can be re-established by using appropriately scaled test norms. As model cases we study the Poisson problem and the Kirchhoff--Love plate bending model, and also include fully discrete variants where optimal test functions are approximated. Numerical experiments for both model problems, including an-isotropic domains and mixed boundary conditions, confirm our findings.


abstract: To counter the volatile nature of renewable energy sources, gas networks take a vital role. But, to ensure fulfillment of contracts under these circumstances, a vast number of possible scenarios, incorporating uncertain supply and demand, has to be simulated ahead of time. This many-query gas network simulation task can be accelerated by model reduction, yet, large-scale, nonlinear, parametric, hyperbolic partial differential(-algebraic) equation systems, modeling natural gas transport, are a challenging application for model order reduction algorithms. For this industrial application, we bring together the scientific computing topics of: mathematical modeling of gas transport networks, numerical simulation of hyperbolic partial differential equation, and parametric model reduction for nonlinear systems. This research resulted in the "morgen" (Model Order Reduction for Gas and Energy Networks) software platform, which enables modular testing of various combinations of models, solvers, and model reduction methods. In this work we present the theoretical background on systemic modeling and structured, data-driven, system-theoretic model reduction for gas networks, as well as the implementation of "morgen" and associated numerical experiments testing model reduction adapted to gas network models.

categories: math.OC cs.CE cs.NA cs.SY eess.SY math.NA


title: GIRAFFE: Representing Scenes as Compositional Generative Neural Feature Fields

abstract: Deep generative models allow for photorealistic image synthesis at high resolutions. But for many applications, this is not enough: content creation also needs to be controllable. While several recent works investigate how to disentangle underlying factors of variation in the data, most of them operate in 2D and hence ignore that our world is three-dimensional. Further, only few works consider the compositional nature of scenes. Our key hypothesis is that incorporating a compositional 3D scene representation into the generative model leads to more controllable image synthesis. Representing scenes as compositional generative neural feature fields allows us to disentangle one or multiple objects from the background as well as individual objects' shapes and appearances while learning from unstructured and unposed image collections without any additional supervision. Combining this scene representation with a neural rendering pipeline yields a fast and realistic image synthesis model. As evidenced by our experiments, our model is able to disentangle individual objects and allows for translating and rotating them in the scene as well as changing the camera pose.


title: Space-time POD-Galerkin approach for parametric flow control

abstract: In this contribution we propose reduced order methods to fast and reliably solve parametrized optimal control problems governed by time dependent nonlinear partial differential equations. Our goal is to provide a tool to deal with the time evolution of several nonlinear optimality systems in many-query context, where a system must be analysed for various physical and geometrical features. Optimal control can be used in order to fill the gap between collected data and mathematical model and it is usually related to very time consuming activities: inverse problems, statistics, etc. Standard discretization techniques may lead to unbearable simulations for real applications. We aim at showing how reduced order modelling can solve this issue. We rely on a space-time POD-Galerkin reduction in order to solve the optimal control problem in a low dimensional reduced space in a fast way for several parametric instances. The proposed algorithm is validated with a numerical test based on environmental sciences: a reduced optimal control problem governed by viscous Shallow Waters Equations parametrized not only in the physics features, but also in the geometrical ones. We will show how the reduced model can be useful in order to recover desired velocity and height profiles more rapidly with respect to the standard simulation, not losing accuracy.


title: Do You Live a Healthy Life? Analyzing Lifestyle by Visual Life Logging

abstract: A healthy lifestyle is the key to better health and happiness and has a considerable effect on quality of life and disease prevention. Current lifelogging/egocentric datasets are not suitable for lifestyle analysis; consequently, there is no research on lifestyle analysis in the field of computer vision. In this work, we investigate the problem of lifestyle analysis and build a visual lifelogging dataset for lifestyle analysis (VLDLA). The VLDLA contains images captured by a wearable camera every 3 seconds from 8:00 am to 6:00 pm for seven days. In contrast to current lifelogging/egocentric datasets, our dataset is suitable for lifestyle analysis as images are taken with short intervals to capture activities of short duration; moreover, images are taken continuously from morning to evening to record all the activities performed by a user. Based on our dataset, we classify the user activities in each frame and use three latent fluents of the user, which change over time and are associated with activities, to measure the healthy degree of the user's lifestyle. The scores for the three latent fluents are computed based on recognized activities, and the healthy degree of the lifestyle for the day is determined based on the scores for the latent fluents. Experimental results show that our method can be used to analyze the healthiness of users' lifestyles.


title: Chemical abundances of Seyfert 2 AGNs IV. Composite models calculated by photoionization + shocks

abstract: We build detailed composite models of photoionization and shock ionization based on the SUMA code to reproduce emission lines emitted from the Narrow Line Regions (NLR) of Seyfert 2 nuclei. The aim of this work is to investigate diagram AGN positions according to shock parameters, shock effects on the gas temperature and ionization structures and derive a semi-empirical abundance calibration based on emission-line ratios little sensitive to the shock presence. The models were used to reproduce optical (3000 < A < 7000) emission line intensities of 244 local (z < 0.4) Seyfert 2s, whose observational data were selected from Sloan Digital Sky Survey DR7. Our models suggest that shocks in Seyfert 2 nuclei have velocities in the range of 50-300 km/s and imply a narrower metallicity range (0.6 < (Z/Z) < 1.6) than those derived using pure photoionization models. Our results indicate that shock velocity in AGNs can not be estimated using standard optical line ratio diagrams, based on integrated spectra. Our models predict a different temperature structure and O+/O and O2+/O fractional abundances throughout the NLR clouds than those derived from pure photoionization models, mainly in shock-dominated objects. This suggests that, in order to minimize the shock effects, the combination of emission-lines emitted by ions with similar intermediate ionization potential could be good metallicity indicators. Finally, we derive two calibrations between the N/O abundance ratio and the N2O2=log([N II]6584/[O II]3727) and N2=log([N II]6584/H{\alpha}) indexes which agree with that derived from pure photoionization models.


title: Recurrent Multi-view Alignment Network for Unsupervised Surface Registration

abstract: Learning non-rigid registration in an end-to-end manner is challenging due to the inherent high degrees of freedom and the lack of labeled training data. In this paper, we resolve these two challenges simultaneously. First, we propose to represent the non-rigid transformation with a point-wise combination of several rigid transformations. This representation not only makes the solution space well-constrained but also enables our method to be solved iteratively with a recurrent framework, which greatly reduces the difficulty of learning. Second, we introduce a differentiable loss function that measures the 3D shape similarity on the projected multi-view 2D depth images so that our full framework can be trained end-to-end without ground truth supervision. Extensive experiments on several different datasets demonstrate that our proposed method outperforms the previous state-of-the-art by a large margin. The source codes are available at https://github.com/WanquanF/RMA-Net.


title: Learning of Long-Horizon Sparse-Reward Robotic Manipulator Tasks with Base Controllers

abstract: Deep Reinforcement Learning (DRL) enables robots to perform some intelligent tasks end-to-end. However, there are still many challenges for long-horizon sparse-reward robotic manipulator tasks. On the one hand, a sparse-reward setting causes exploration inefficient. On the other hand, exploration using physical robots is of high cost and unsafe. In this paper, we propose a method of learning long-horizon sparse-reward tasks utilizing one or more existing traditional controllers named base controllers in this paper. Built upon Deep Deterministic Policy Gradients (DDPG), our algorithm incorporates the existing base controllers into stages of exploration, value learning, and policy update. Furthermore, we present a straightforward way of synthesizing different base controllers to integrate their strengths. Through experiments ranging from stacking blocks to cups, it is demonstrated that the learned state-based or image-based policies steadily outperform base controllers. Compared to previous works of learning from demonstrations, our method improves sample efficiency by orders of magnitude and improves the performance. Overall, our method bears the potential of leveraging existing industrial robot manipulation systems to build more flexible and intelligent controllers.


title: Flat replacements of homology theories

abstract: To a homology theory one can associate an additive site and a new homological functor with values in the category of additive sheaves on that site. If this category of sheaves can be shown to be equivalent to a category of comodules of a Hopf algebroid, then we obtain a new homology theory by composing with the underlying module functor. This new homology theory is always flat and we call it a flat replacement of the original theory. For example, Pstr\k{a}gowski has shown that complex cobordism is a flat replacement of singular homology. In this article we study the basic properties of the sites associated to homology theories and we prove an existence theorem for flat replacements.


title: EEG-GCNN: Augmenting Electroencephalogram-based Neurological Disease Diagnosis using a Domain-guided Graph Convolutional Neural Network

abstract: This paper presents a novel graph convolutional neural network (GCNN)-based approach for improving the diagnosis of neurological diseases using scalp-electroencephalograms (EEGs). Although EEG is one of the main tests used for neurological-disease diagnosis, the sensitivity of EEG-based expert visual diagnosis remains at $\sim$50\%. This indicates a clear need for advanced methodology to reduce the false negative rate in detecting abnormal scalp-EEGs. In that context, we focus on the problem of distinguishing the abnormal scalp EEGs of patients with neurological diseases, which were originally classified as 'normal' by experts, from the scalp EEGs of healthy individuals. The contributions of this paper are three-fold: 1) we present EEG-GCNN, a novel GCNN model for EEG data that captures both the spatial and functional connectivity between the scalp electrodes, 2) using EEG-GCNN, we perform the first large-scale evaluation of the aforementioned hypothesis, and 3) using two large scalp-EEG databases, we demonstrate that EEG-GCNN significantly outperforms the human baseline and classical machine learning (ML) baselines, with an AUC of 0.90.


title: Wide-angle Image Rectification: A Survey

abstract: Wide field-of-view (FOV) cameras, which capture a larger scene area than narrow FOV cameras, are used in many applications including 3D reconstruction, autonomous driving, and video surveillance. However, wide-angle images contain distortions that violate the assumptions underlying pinhole camera models, resulting in object distortion, difficulties in estimating scene distance, area, and direction, and preventing the use of off-the-shelf deep models trained on undistorted images for downstream computer vision tasks. Image rectification, which aims to correct these distortions, can solve these problems. In this paper, we comprehensively survey progress in wide-angle image rectification from transformation models to rectification methods. Specifically, we first present a detailed description and discussion of the camera models used in different approaches. Then, we summarize several distortion models including radial distortion and projection distortion. Next, we review both traditional geometry-based image rectification methods and deep learning-based methods, where the former formulate distortion parameter estimation as an optimization problem and the latter treat it as a regression problem by leveraging the power of deep neural networks. We evaluate the performance of state-of-the-art methods on public datasets and show that although both kinds of methods can achieve good results, these methods only work well for specific camera models and distortion types. We also provide a strong baseline model and carry out an empirical study of different distortion models on synthetic datasets and real-world wide-angle images. Finally, we discuss several potential research directions that are expected to further advance this area in the future.


title: Predicting of shear wave velocity using Artificial Neural Networks

abstract: Shear wave velocity is an important parameter for determining lithology, porosity and the dynamic properties in geo-mechanical studies. However, due to time and cost limitations, shear wave velocity is not available at all intervals and in all wells. In this paper, well logs with strong correlation to shear wave velocity were determined and used to predict the shear wave velocity for the Tano North Field. Four different methods were used to estimate the shear wave velocity under three different conditions. Then, based on obtained coefficient of determination and average absolute percent relative error between real and predicted values of shear wave velocity, the final results were compared. The results of this work demonstrated that the neural network based on multiple variables can estimate the shear wave velocity better than other methods used.


title: Reconfigurable Intelligent Surfaces in Challenging Environments: Underwater, Underground, Industrial and Disaster

abstract: Challenging environments comprise a range of scenarios, which share the fact that it is extremely difficult to establish a communication link using conventional technology due to many impairments typically associated with the propagation medium and increased signal scattering. Specifically, underwater and underground media are known to absorb electromagnetic radiation, which heavily affects the overall path loss. Industrial and disaster environments can be viewed as rich scattering environments with corresponding substantial multipath propagation leading to intersymbol interference and deterioration of signal quality. Although the challenges for the design of communication networks, and specifically the Internet of Things (IoT), in such environments are known, there is no common enabler or solution for all these applications. Reconfigurable intelligent surfaces (RISs) have been introduced to improve the signal propagation characteristics by focusing the signal power in the preferred direction, thus making the communication environment 'smart'. While the usual application of RIS is related to blockage avoidance, the very same technique can be used to reduce the effect of multipath and even partially compensate the signal absorption via passive beamforming. Due to the beneficial properties of RIS, its use in challenging environments can become the aforementioned enabler and a game changing technology. In this paper, we discuss potential use cases, deployment strategies and design aspects for RIS devices in underwater IoT, underground IoT as well as Industry 4.0 and emergency networks. Furthermore, we provide a potential hardware architecture and derive the expected signal quality. The numerical results reveal substantial performance gains of up to 20 dB per decade. In addition, novel research challenges to be addressed in this context are described.


title: Wavelet-based clustering for time-series trend detection

abstract: In this paper, we introduce a method performing clustering of time-series on the basis of their trend (increasing, stagnating/decreasing, and seasonal behavior). The clustering is performed using $k$-means method on a selection of coefficients obtained by discrete wavelet transform, reducing drastically the dimensionality. The method is applied on an use case for the clustering of a 864 daily sales revenue time-series for 61 retail shops. The results are presented for different mother wavelets. The importance of each wavelet coefficient and its level is discussed thanks to a principal component analysis along with a reconstruction of the signal from the selected wavelet coefficients.


title: Point vortices dynamics on a rotating sphere and modeling of global atmospheric vortices interaction

abstract: It is shown that the hydrodynamics equations for a thin spherical liquid layer are satisfied by the stream function of a pair of antipodal vortices-APV, in contrast to the stream function of a single point vortex on a sphere with a background of a uniform opposite sign vorticity. A simple zero solution of the equation of the absolute vorticity conservation is used for bypassing well-known nonlinear problem of a point vortices interaction with regular vorticity field and an exact solution for APVs dynamics problem on a rotating sphere is obtained. Due to this a new stable stationary solution for the dynamics of APV is obtained, which can model the dynamics of the global vortex structures such as atmospheric centers of action.


title: Automatic artifact removal of resting-state fMRI with Deep Neural Networks

abstract: Functional Magnetic Resonance Imaging (fMRI) is a non-invasive technique for studying brain activity. During an fMRI session, the subject executes a set of tasks (task-related fMRI study) or no tasks (resting-state fMRI), and a sequence of 3-D brain images is obtained for further analysis. In the course of fMRI, some sources of activation are caused by noise and artifacts. The removal of these sources is essential before the analysis of the brain activations. Deep Neural Network (DNN) architectures can be used for denoising and artifact removal. The main advantage of DNN models is the automatic learning of abstract and meaningful features, given the raw data. This work presents advanced DNN architectures for noise and artifact classification, using both spatial and temporal information in resting-state fMRI sessions. The highest performance is achieved by a voting schema using information from all the domains, with an average accuracy of over 98% and a very good balance between the metrics of sensitivity and specificity (98.5% and 97.5% respectively).


title: On Solar Photovoltaic Parameter Estimation: Global Optimality Analysis and a Simple Efficient Differential Evolution Method

abstract: A large variety of sophisticated metaheuristic methods have been proposed for photovoltaic parameter extraction. Our aim is not to develop another metaheuristic method but to investigate two practically important yet rarely studied issues: (i) whether existing results are already globally optimal; (ii) whether a significantly simpler metaheuristic can achieve equally good performance. We take the two widely used I-V curve datasets for case studies. The first issue is addressed using a branch and bound algorithm, which certifies the global minimum rigorously or locates a fairly tight upper bound, despite its intolerable slowness. These values are useful references for fair evaluation and further development of metaheuristics. Next, extensive examination and comparison reveal that, perhaps surprisingly, an elementary differential evolution (DE) algorithm can either attain the global minimum certified above or obtain the best-known result. More attractively, the simple DE algorithm takes only a fraction of the runtime of state-of-the-art metaheuristic methods and is particularly preferable in time-sensitive applications. This novel, unusual, and notable finding also indicates that the employment of increasingly complicated metaheuristics might be somewhat overkilling for regular PV parameter estimation. Finally, we discuss the implications of these results for future research and suggest the simple DE method as the first choice for industrial applications.