title: Exploring Inherent Properties of the Monophonic Melody of Songs

abstract: Melody is one of the most important components in music. Unlike other components in music theory, such as harmony and counterpoint, computable features for melody is urgently in need. These features are highly demanded as data-driven methods dominating the fields such as musical information retrieval and automatic music composition. To boost the performance of deep-learning-related musical tasks, we propose a set of interpretable features on monophonic melody for computational purposes. These features are defined not only in mathematical form, but also with some considerations on composers 'intuition. For example, the Melodic Center of Gravity can reflect the sentence-wise contour of the melody, the local / global melody dynamics quantifies the dynamics of a melody that couples pitch and time in a sentence. We found that these features are considered by people universally in many genres of songs, even for atonal composition practices. Hopefully, these melodic features can provide nov el inspiration for future researchers as a tool in the field of MIR and automatic composition.


title: FedNER: Privacy-preserving Medical Named Entity Recognition with Federated Learning

abstract: Medical named entity recognition (NER) has wide applications in intelligent healthcare. Sufficient labeled data is critical for training accurate medical NER model. However, the labeled data in a single medical platform is usually limited. Although labeled datasets may exist in many different medical platforms, they cannot be directly shared since medical data is highly privacy-sensitive. In this paper, we propose a privacy-preserving medical NER method based on federated learning, which can leverage the labeled data in different platforms to boost the training of medical NER model and remove the need of exchanging raw data among different platforms. Since the labeled data in different platforms usually has some differences in entity type and annotation criteria, instead of constraining different platforms to share the same model, we decompose the medical NER model in each platform into a shared module and a private module. The private module is used to capture the characteristics of the local data in each platform, and is updated using local labeled data. The shared module is learned across different medical platform to capture the shared NER knowledge. Its local gradients from different platforms are aggregated to update the global shared module, which is further delivered to each platform to update their local shared modules. Experiments on three publicly available datasets validate the effectiveness of our method.


title: Sharp Liouville Theorems

abstract: Consider the equation div$(\varphi^2 \nabla \sigma)=0$ in $\mathbb{R}^N,$ where $\varphi>0$. Berestycki, Caffarelli and Nirenberg proved that if there exists $C>0$ such that $\int_{B_R}(\varphi \sigma)^2 \leq CR^2$ for every $R\geq 1$ then $\sigma$ is necessarily constant. In this paper we provide necessary and sufficient conditions on $0<\Psi\in C([1,\infty))$ for which this result remains true if we replace $R^2$ with $\Psi(R)$ in any dimension $N$. In the case of the convexity of $\Psi$ for large $R>1$ and $\Psi'>0$, this condition is equivalent to $\displaystyle{\int_1^\infty\frac{1}{\Psi'}=\infty}$.


title: Effective theory for self-interacting dark matter and massive spin-2 mediators

abstract: We consider the effective theory for self-interacting dark matter with arbitrary spin and go beyond the previous discussion in the literature by introducing a massive spin-2 particle as the mediator for the dark matter self-scattering. We present the effective self-interactions for dark matter in the leading order expansions with the momentum transfer and the dark matter velocity. We compare between the Born cross section and the non-perturbative cross section in the leading order approximation of the effective Yukawa interaction. As a result, we find that there is a wide range of dark matter and spin-2 particle masses for a velocity-dependent self-scattering to solve small-scale problems at galaxies as well as satisfy the bounds from galaxy clusters at the same time.


title: Improving Irregularly Sampled Time Series Learning with Dense Descriptors of Time

abstract: Supervised learning with irregularly sampled time series have been a challenge to Machine Learning methods due to the obstacle of dealing with irregular time intervals. Some papers introduced recently recurrent neural network models that deals with irregularity, but most of them rely on complex mechanisms to achieve a better performance. This work propose a novel method to represent timestamps (hours or dates) as dense vectors using sinusoidal functions, called Time Embeddings. As a data input method it and can be applied to most machine learning models. The method was evaluated with two predictive tasks from MIMIC III, a dataset of irregularly sampled time series of electronic health records. Our tests showed an improvement to LSTM-based and classical machine learning models, specially with very irregular data.


title: Sun-Heliosphere Observation-based Ionization Rates Model

abstract: The solar wind (SW) and the extreme ultraviolet (EUV) radiation modulate fluxes of interstellar and heliospheric particles inside the heliosphere both in time and in space. Understanding this modulation is necessary to correctly interpret measurements of particles of interstellar origin inside the heliosphere. We present a revision of heliospheric ionization rates and provide the Sun-Heliosphere Observation-based Ionization Rates (SHOIR) model based on the currently available data. We calculate the total ionization rates using revised SW and solar EUV data. We study the in-ecliptic variation of the SW parameters, the latitudinal structure of the SW speed and density, and the reconstruction of the photoionization rates. The revision most affects the SW out of the ecliptic plane during solar maximum and the estimation of the photoionization rates, the latter due to a change of the reference data. The revised polar SW is slower and denser during the solar maximum of solar cycle (SC) 24. The current estimated total ionization rates are higher than the previous ones for H, O, and Ne, and lower for He. The changes for the in-ecliptic total ionization rates are less than 10% for H and He, up to 20% for O, and up to 35% for Ne. Additionally, the changes are not constant in time and vary as a function of time and latitude.


title: U-Det: A Modified U-Net architecture with bidirectional feature network for lung nodule segmentation

abstract: Early diagnosis and analysis of lung cancer involve a precise and efficient lung nodule segmentation in computed tomography (CT) images. However, the anonymous shapes, visual features, and surroundings of the nodule in the CT image pose a challenging problem to the robust segmentation of the lung nodules. This article proposes U-Det, a resource-efficient model architecture, which is an end to end deep learning approach to solve the task at hand. It incorporates a Bi-FPN (bidirectional feature network) between the encoder and decoder. Furthermore, it uses Mish activation function and class weights of masks to enhance segmentation efficiency. The proposed model is extensively trained and evaluated on the publicly available LUNA-16 dataset consisting of 1186 lung nodules. The U-Det architecture outperforms the existing U-Net model with the Dice similarity coefficient (DSC) of 82.82% and achieves results comparable to human experts.


title: Self-Supervised Light Field Reconstruction Using Shearlet Transform and Cycle Consistency

abstract: The image-based rendering approach using Shearlet Transform (ST) is one of the state-of-the-art Densely-Sampled Light Field (DSLF) reconstruction methods. It reconstructs Epipolar-Plane Images (EPIs) in image domain via an iterative regularization algorithm restoring their coefficients in shearlet domain. Consequently, the ST method tends to be slow because of the time spent on domain transformations for dozens of iterations. To overcome this limitation, this letter proposes a novel self-supervised DSLF reconstruction method, CycleST, which applies ST and cycle consistency to DSLF reconstruction. Specifically, CycleST is composed of an encoder-decoder network and a residual learning strategy that restore the shearlet coefficients of densely-sampled EPIs using EPI reconstruction and cycle consistency losses. Besides, CycleST is a self-supervised approach that can be trained solely on Sparsely-Sampled Light Fields (SSLFs) with small disparity ranges ($\leqslant$ 8 pixels). Experimental results of DSLF reconstruction on SSLFs with large disparity ranges (16 - 32 pixels) from two challenging real-world light field datasets demonstrate the effectiveness and efficiency of the proposed CycleST method. Furthermore, CycleST achieves ~ 9x speedup over ST, at least.


title: The TESS light curve of AI Phoenicis

abstract: Accurate masses and radii for normal stars derived from observations of detached eclipsing binary stars are of fundamental importance for testing stellar models and may be useful for calibrating free parameters in these model if the masses and radii are sufficiently precise and accurate. We aim to measure precise masses and radii for the stars in the bright eclipsing binary AI Phe, and to quantify the level of systematic error in these estimates. We use several different methods to model the TESS light curve of AI Phe combined with spectroscopic orbits from multiple sources to estimate precisely the stellar masses and radii together with robust error estimates. We find that the agreement between different methods for the light curve analysis is very good but some methods underestimate the errors on the model parameters. The semi-amplitudes of the spectroscopic orbits derived from spectra obtained with modern echelle spectrographs are consistent to within 0.1%. The masses of the stars in AI Phe are $M_1 = 1.1938 \pm 0.0008 M_{\odot}$ and $M_2 = 1.2438 \pm 0.0008M_{\odot}$, and the radii are $R_1 = 1.8050 \pm 0.0022 R_{\odot}$ and $R_2 = 2.9332 \pm 0.0023 R_{\odot}$. We conclude that it is possible to measure accurate masses and radii for stars in bright eclipsing binary stars to a precision of 0.2% or better using photometry from TESS and spectroscopy obtained with modern echelle spectrographs. We provide recommendations for publishing masses and radii of eclipsing binary stars at this level of precision.


title: Chaotic dynamics of graphene and graphene nanoribbons

abstract: We study the chaotic dynamics of graphene structures, considering both a periodic, defect free, graphene sheet and graphene nanoribbons (GNRs) of various widths. By numerically calculating the maximum Lyapunov exponent, we quantify the chaoticity for a spectrum of energies in both systems. We find that for all cases, the chaotic strength increases with the energy density, and that the onset of chaos in graphene is slow, becoming evident after more than $10^4$ natural oscillations of the system. For the GNRs, we also investigate the impact of the width and chirality (armchair or zigzag edges) on their chaotic behavior. Our results suggest that due to the free edges the chaoticity of GNRs is stronger than the periodic graphene sheet, and decreases by increasing width, tending asymptotically to the bulk value. In addition, the chaotic strength of armchair GNRs is higher than a zigzag ribbon of the same width. Further, we show that the composition of ${}^{12}C$ and ${}^{13}C$ carbon isotopes in graphene has a minor impact on its chaotic strength.


title: Dynamic Averaging Load Balancing on Cycles

abstract: We consider the following dynamic load-balancing process: given an underlying graph $G$ with $n$ nodes, in each step $t\geq 0$, one unit of load is created, and placed at a randomly chosen graph node. In the same step, the chosen node picks a random neighbor, and the two nodes balance their loads by averaging them. We are interested in the expected gap between the minimum and maximum loads at nodes as the process progresses, and its dependence on $n$ and on the graph structure. Similar variants of the above graphical balanced allocation process have been studied by Peres, Talwar, and Wieder, and by Sauerwald and Sun for regular graphs. These authors left as open the question of characterizing the gap in the case of \emph{cycle graphs} in the \emph{dynamic} case, where weights are created during the algorithm's execution. For this case, the only known upper bound is of $\mathcal{O}( n \log n )$, following from a majorization argument due to Peres, Talwar, and Wieder, which analyzes a related graphical allocation process. In this paper, we provide an upper bound of $\mathcal{O} ( \sqrt n \log n )$ on the expected gap of the above process for cycles of length $n$. We introduce a new potential analysis technique, which enables us to bound the difference in load between $k$-hop neighbors on the cycle, for any $k \leq n / 2$. We complement this with a "gap covering" argument, which bounds the maximum value of the gap by bounding its value across all possible subsets of a certain structure, and recursively bounding the gaps within each subset. We provide analytical and experimental evidence that our upper bound on the gap is tight up to a logarithmic factor.


title: Power Assisted Trend Following

abstract: 'The trend is your friend' is a common saying, the difficulty lies in determining if and when you are in a trend. Is the trend strong enough to trade? When does the trend reverse and how are you going to determine this? We will try and answer at least some of these questions here. We are deriving a novel indicator to measure the power of a trend using digital signal processing techniques, separating the Signal from the Noise. We apply these to examples as well as real data and evaluate the accuracy of these and the relation to PNL performance of the 'Volatility Index' trend following algorithm devised by J. Welles Wilder Jr. in 1978.


title: Neutron-hole strength in N = 81 nuclei

abstract: A systematic study of neutron-hole strength in the N = 81 nuclei 137Ba, 139Ce, 141Nd and 143Sm is reported. The single-neutron removal reactions (p,d) and (3He,4He) were measured at energies of 23 and 34 MeV, respectively. Spectroscopic factors were extracted from measured cross sections through a distorted-wave Born approximation analysis and centroids of single-particle strength have been established. The change in these centroid energies as a function of proton number have been compared to calculations of the monopole shift for the s1/2 and h11/2 orbitals, where the majority of the strength has been observed. Significant fragmentation of strength was observed for the d and g7/2 orbitals, particularly for the latter orbital which is deeply bound, with summed strengths that indicate a significant amount lies outside of the measured excitation energy range.


title: Graham's Formula for Valuing Growth Stocks

abstract: Benjamin Graham introduced a very simple formula for valuing a growth stock in 1962. How does it work and why? What is a sensible way to calculate this across many stocks and provide a scoring system to compare stocks amongst each other? We are presenting a methodology here which is put into practice.


title: Distributed and Democratized Learning: Philosophy and Research Challenges

abstract: Due to the availability of huge amounts of data and processing abilities, current artificial intelligence (AI) systems are effective in solving complex tasks. However, despite the success of AI in different areas, the problem of designing AI systems that can truly mimic human cognitive capabilities such as artificial general intelligence, remains largely open. Consequently, many emerging cross-device AI applications will require a transition from traditional centralized learning systems towards large-scale distributed AI systems that can collaboratively perform multiple complex learning tasks. In this paper, we propose a novel design philosophy called democratized learning (Dem-AI) whose goal is to build large-scale distributed learning systems that rely on the self-organization of distributed learning agents that are well-connected, but limited in learning capabilities. Correspondingly, inspired by the societal groups of humans, the specialized groups of learning agents in the proposed Dem-AI system are self-organized in a hierarchical structure to collectively perform learning tasks more efficiently. As such, the Dem-AI learning system can evolve and regulate itself based on the underlying duality of two processes which we call specialized and generalized processes. In this regard, we present a reference design as a guideline to realize future Dem-AI systems, inspired by various interdisciplinary fields. Accordingly, we introduce four underlying mechanisms in the design such as plasticity-stability transition mechanism, self-organizing hierarchical structuring, specialized learning, and generalization. Finally, we establish possible extensions and new challenges for the existing learning approaches to provide better scalable, flexible, and more powerful learning systems with the new setting of Dem-AI.


title: The HI mass function of group galaxies in the ALFALFA survey

abstract: We estimate the HI mass function (HIMF) of galaxies in groups based on thousands of ALFALFA (Arecibo Legacy Fast ALFA survey) HI detections within the galaxy groups of four widely used SDSS (Sloan Digital Sky Survey) groups catalogues. Although differences between the catalogues mean that there is no one definitive group galaxy HIMF, in general we find that the low-mass slope is flat, in agreement with studies based on small samples of individual groups, and that the 'knee' mass is slightly higher than that of the global HIMF of the full ALFALFA sample. We find that the observed fraction of ALFALFA galaxies in groups is approximately 22 per cent. These group galaxies were removed from the full ALFALFA source catalogue to calculate the field HIMF using the remaining galaxies. Comparison between the field and group HIMFs reveals that group galaxies make only a small contribution to the global HIMF as most ALFALFA galaxies are in the field, but beyond the HIMF 'knee' group galaxies dominate. Finally we attempt to separate the group galaxy HIMF into bins of group halo mass, but find that too few low-mass galaxies are detected in the most massive groups to tightly constrain the slope, owing to the rarity of such groups in the nearby Universe where low-mass galaxies are detectable with existing HI surveys.


title: Second Order Twist Contributions to the Balitsky-Kovchegov Equation at small-x: Deterministic and Stochastic Pictures

abstract: We study the second order twist corrections to a toy model of a dipole-dipole interaction in the context of a both deterministic and stochastic effects. This work is done in the high $N_C$ limit in the Bjorken picture. We show that the correction to the second twist terms of the stochastic picture suggest additional importance of the second twist correction in the stochastic model as compared to the deterministic model.


title: Gaussian null coordinates, near-horizon geometry and conserved charges on the horizon of extremal non-dilatonic black $p$-branes

abstract: In this paper, we examine the emergence of conserved charges on the horizon of a particular class of extremal non-dilatonic black $p$-branes (which reduce to extremal dilatonic black holes in $D=4$ dimensions upon toroidal compactification) in the presence of a probe massless scalar field in the bulk. This result is achieved by writing the black $p$-brane geometry in a Gaussian null coordinate system which allows us to get a non-singular near-horizon geometry description. We find that the near-horizon geometry is $AdS_{p+2}\times S^2$ and that the $AdS_{p+2}$ section has an internal structure which can be seen as a warped product of $AdS_{2}\times S^{p}$ in Gaussian null coordinates. We show that the bulk scalar field satisfying the field equations is expanded in terms of non-normalizable and normalizable modes, which for certain suitable quantization conditions are well-behaved at the boundary of $AdS_{p+2}$ space. Furthermore, we show that picking the normalizable modes results in the existence of conserved quantities on the horizon. We discuss the impact of these conserved quantities in the late time regime.


title: Expectations for the confirmation of Proxima c from a long-term radial velocity follow-up

abstract: Proxima c, a candidate second planet orbiting Proxima Centauri, was detected with the radial velocity method. The announced long orbital period (5.21$^{+0.26}_{-0.22}$ years), and small semi-amplitude of the induced Doppler signal (1.2$\pm$0.4 m/s), make this detection challenging and the target worthy of a follow-up in the next years. We intend to evaluate the impact of future data on the statistical significance of the detection through realistic simulated radial velocities to be added to the published dataset, spanning up to one orbital period of Proxima c in the time range 2019-2023. We find that the detection significance of Proxima c increases depending not only on the amount of data collected, but also on the number of instruments used, and especially on the timespan covered by the observational campaign. However, on average we do not get strong statistical evidence and we predict that, in the best-case scenario, in the next 5 years the detection of Proxima c can become significant at 4$\sigma$ level. If instead Proxima c does not exist, the detected signal may lower its significance down to 2$\sigma$.


title: Modeling the fringe field of permanent magnet multipoles using numerical simulations

abstract: Permanent Magnet multipoles (PMM) are widely used in accelerators to either focus particle beams or confine plasma in ion sources. The real magnetic field created by PMM is calculated by magnetic field simulation software and then used in particle tracking codes by means of 3 dimensional magnetic field map. A common alternative is to use the so-called 'hard edge' model, which gives an approximation of the magnetic field inside the PMM assuming a null fringe field. This work proposes an investigation of the PMM fringe field properties. An analytical model of PMM magnetic field is developed using the Fourier multipole expansion. A general axial potential function with a unique parameter ${\lambda}$, able to reproduce the actual PMM magnetic field (including its two fringe fields) with an explicit dependence on the PMM length is proposed. An analytical first order model including the axial fringe field is derived. This simple model complies with the Maxwell equations (curl(B)=0 and div(B)=0) and can replace advantageously the 'hard edge' model when fast analytical calculation are required. Higher order analytical multipole expansion model quality is assessed by means of ${\chi^2}$ estimators. The general dependence of the potential function parameter ${\lambda}$ is given as a function of the PMM geometry for quadrupole, hexapole and multipole, allowing to use the developed model in simulation programs where the multipole geometry is an input parameter.


title: The effect of diffusive nuclear burning in neutron star envelopes on cooling in accreting systems

abstract: Valuable information about the neutron star interior can be obtained by comparing observations of thermal radiation from a cooling neutron star crust with theoretical models. Nuclear burning of lighter elements that diffuse to deeper layers of the envelope can alter the relation between surface and interior temperatures and can change the chemical composition over time. We calculate new temperature relations and consider two effects of diffusive nuclear burning (DNB) for H-C envelopes. First, we consider the effect of a changing envelope composition and find that hydrogen is consumed on short timescales and our temperature evolution simulations correspond to those of a hydrogen-poor envelope within ~100 days. The transition from a hydrogen-rich to a hydrogen-poor envelope is potentially observable in accreting NS systems as an additional initial decline in surface temperature at early times after the outburst. Second, we find that DNB can produce a non-negligible heat flux, such that the total luminosity can be dominated by DNB in the envelope rather than heat from the deep interior. However, without continual accretion, heating by DNB in H-C envelopes is only relevant for <1-80 days after the end of an accretion outburst, as the amount of light elements is rapidly depleted. Comparison to crust cooling data shows that DNB does not remove the need for an additional shallow heating source. We conclude that solving the time-dependent equations of the burning region in the envelope self-consistently in thermal evolution models instead of using static temperature relations would be valuable in future cooling studies.


title: Neutron Diffraction Evidence for Local Spin Canting, Weak Jahn-Teller Distortion, and Magnetic Compensation in Ti$_{1-x}$Mn$_{x}$Co$_2$O$_4$ Spinel

abstract: A systematic study using neutron diffraction and magnetic susceptibility are reported on Mn substituted ferrimagnetic inverse spinel Ti$_{1-x}$Mn$_{x}$Co$_2$O$_4$ in the temperature interval 1.6 K $\leq$ $T$ $\leq$ 300 K. Our neutron diffraction study reveals cooperative distortions of the $T$O$_6$ octahedral for all the Jahn-Teller active ions $T$ = Mn$^{3+}$, Ti$^{3+}$ and Co$^{3+}$, which are confirmed by the X-ray photoelectron spectroscopy. Two specific compositions ($x$ = 0.2 and 0.4) have been chosen because of their unique features: noncollinear Yafet-Kittel type ordering, and weak tetragonal distortion with ${c/a}$ $<$ 1, in which the apical bond length $d_c$($T_B$-O) is longer than the equatorial $d_{ab}$($T_B$-O) due to the splitting of the $e_g$ level of Mn$^{3+}$ ions into $d_{x^2-y^2}$ and $d_{z^2}$. For $x$ = 0.4, the distortion in the $T_B$O$_6$ octahedra is stronger as compared to $x$ = 0.2 because of the higher content of trivalent Mn. Ferrimagnetic ordering in $x$ = 0.4 and $x$ = 0.2 sets in at 110.3 and 78.2 K, respectively due to the unequal magnetic moments of cations, where Ti$^{3+}$, Mn$^{3+}$, and Co$^{3+}$ occupying the octahedral, whereas, Co$^{2+}$ sits in the tetrahedral site. In addition, weak antiferromagnetic component could be observed lying perpendicular to the ferrimagnetic component. The analysis of static and dynamic magnetic susceptibilities combined with the heat-capacity data reveals a magnetic compensation phenomenon at $T_{COMP}$ = 25.4 K in $x$ = 0.2 and a reentrant spin-glass behaviour in $x$ = 0.4 with a freezing temperature $\sim$110.1 K. The compensation phenomenon is characterized by sign reversal of magnetization and bipolar exchange bias effect below $T_{COMP}$ with its magnitude depending on the direction of external magnetic field and the cooling protocol.


title: Drying of porous media by concurrent drainage and evaporation: A pore network modeling study

abstract: Drainage and evaporation can occur simultaneously during the drying of porous media, but the interactions between these processes and their effects on drying are rarely studied. In this work, we develop a pore network model that considers drainage, evaporation, and rarefied multi-component gas transport in porous media with nanoscale pores. Using this model, we investigate the drying of a liquid solvent-saturated porous medium enabled by the flow of purge gas through it. Simulations show that drying progresses in three stages, and the solvent removal by drainage effects (evaporation effects) becomes increasingly weak (strong) as drying progresses through these stages. Interestingly, drainage can contribute considerably to solvent removal even after evaporation effects become very strong, especially when the applied pressure difference across the porous medium is low. We show that these phenomena are the results of the coupling between the drainage and evaporation effects and this coupling depends on the operating conditions and the stage of drying.


title: \.Insans{\i}z Ara\c{c}larla D\"uzlemsel Olmayan Ara\c{c}lar{\i}n Taranmas{\i}

abstract: The importance of area coverage with unmanned vehicles, in other words, traveling an area with an unmanned vehicle such as a robot or a UAV completely or partially with minimum cost, is increasing with the increase in usage of such vehicles today. Area coverage with unmanned vehicles is used today in the exploration of an area with UAVs, sweeping mines with robots, cleaning ground with robots in large shopping malls, mowing lawn in a large area etc. The problem has versions such as area coverage with a single unmanned vehicle, area coverage with multiple unmanned vehicles, on-line area coverage (The map of the area that will be covered is not known before starting the coverage) with unmanned vehicles etc. In addition, the area may have obstacles that the vehicles cannot move over. Naturally, many researches are working on the problem and a lot of researches have been done on the problem until today. Spanning tree coverage is one of the major approaches to the problem. In this approach, at the basic level, the planar area is divided into identical squares according to the range of sight of the vehicle, and centers of these squares are assumed to be vertexes of a graph. The vertexes of this graph are connected with the edges with unit costs and after finding the minimum spanning tree of the graph, the vehicle strolls around the spanning tree. The method we propose suggests a way to cover a non-planar area with unmanned vehicles. The method we propose also takes advantage of the spanning-tree coverage approach, but instead of assigning unit costs to the edges, we assigned a weight to each edge using slopes between vertexes those the edges connect. We have gotten noticeably better results than the results we got when we did not consider the slope between two squares and used the classical spanning tree approach.


title: Drift-Adjusted And Arbitrated Ensemble Framework For Time Series Forecasting

abstract: Time Series Forecasting is at the core of many practical applications such as sales forecasting for business, rainfall forecasting for agriculture and many others. Though this problem has been extensively studied for years, it is still considered a challenging problem due to complex and evolving nature of time series data. Typical methods proposed for time series forecasting modeled linear or non-linear dependencies between data observations. However it is a generally accepted notion that no one method is universally effective for all kinds of time series data. Attempts have been made to use dynamic and weighted combination of heterogeneous and independent forecasting models and it has been found to be a promising direction to tackle this problem. This method is based on the assumption that different forecasters have different specialization and varying performance for different distribution of data and weights are dynamically assigned to multiple forecasters accordingly. However in many practical time series data-set, the distribution of data slowly evolves with time. We propose to employ a re-weighting based method to adjust the assigned weights to various forecasters in order to account for such distribution-drift. An exhaustive testing was performed against both real-world and synthesized time-series. Experimental results show the competitiveness of the method in comparison to state-of-the-art approaches for combining forecasters and handling drift.


title: Health State Estimation

abstract: Life's most valuable asset is health. Continuously understanding the state of our health and modeling how it evolves is essential if we wish to improve it. Given the opportunity that people live with more data about their life today than any other time in history, the challenge rests in interweaving this data with the growing body of knowledge to compute and model the health state of an individual continually. This dissertation presents an approach to build a personal model and dynamically estimate the health state of an individual by fusing multi-modal data and domain knowledge. The system is stitched together from four essential abstraction elements: 1. the events in our life, 2. the layers of our biological systems (from molecular to an organism), 3. the functional utilities that arise from biological underpinnings, and 4. how we interact with these utilities in the reality of daily life. Connecting these four elements via graph network blocks forms the backbone by which we instantiate a digital twin of an individual. Edges and nodes in this graph structure are then regularly updated with learning techniques as data is continuously digested. Experiments demonstrate the use of dense and heterogeneous real-world data from a variety of personal and environmental sensors to monitor individual cardiovascular health state. State estimation and individual modeling is the fundamental basis to depart from disease-oriented approaches to a total health continuum paradigm. Precision in predicting health requires understanding state trajectory. By encasing this estimation within a navigational approach, a systematic guidance framework can plan actions to transition a current state towards a desired one. This work concludes by presenting this framework of combining the health state and personal graph model to perpetually plan and assist us in living life towards our goals.


title: Infinite populations of migrants as complex systems: self-regulation

abstract: A model is proposed and studied describing an infinite population of point migrants arriving in and departing from $X\subseteq \mathbf{R}^d$, $d\geq 1$. Both these acts occur at random with state-dependent rates. That is, depending on their geometry the existing migrants repel and attract the newcomers, which makes the population a complex system. Its states are probability measures on an appropriate configuration space, and their evolution $\mu_0 \to \mu_t$ is obtained by solving the corresponding Fokker-Planck equation. The main result is the conclusion that this evolution of states preserves their sub-Poissonicity, and hence a local self-regulation (suppression of clustering) takes place due to the inter-particle repulsion -- no matter of how small range. Further possibilities to study the proposed model with the help of this result are also discussed.


title: New heuristics for burning graphs

abstract: The concept of graph burning and burning number ($bn(G)$) of a graph G was introduced recently [1]. Graph burning models the spread of contagion (fire) in a graph in discrete time steps. $bn(G)$ is the minimum time needed to burn a graph $G$.The problem is NP-complete. In this paper, we develop first heuristics to solve the problem in general (connected) graphs. In order to test the performance of our algorithms, we applied them on some graph classes with known burning number such as theta graphs, we tested our algorithms on DIMACS and BHOSLIB that are known benchmarks for NP-hard problems in graph theory. We also improved the upper bound for burning number on general graphs in terms of their distance to cluster. Then we generated a data set of 2000 random graphs with known distance to cluster and tested our heuristics on them.


title: Predicting the number of coauthors for researchers: A learning model

abstract: Predicting the number of coauthors for researchers contributes to understanding the development of team science. However, it is an elusive task due to diversity in the collaboration patterns of researchers. This study provides a learning model for the dynamics of this variable; the parameters are learned from empirical data that consist of the number of publications and the number of coauthors at given time intervals. The model is based on relationship between the annual number of new coauthors and time given an annual number of publications, the relationship between the annual number of publications and time given a historical number of publications, and Lotka's law. The assumptions of the model are validated by applying it on the high-quality dblp dataset. The effectiveness of the model is tested on the dataset by satisfactory fittings on the evolutionary trend of the number of coauthors for researchers, the distribution of this variable, and the occurrence probability of collaboration events. Due to its regression nature, the model has the potential to be extended to assess the confidence level of the prediction results and thus has applicability to other empirical research.


title: Detection of Information Hiding at Anti-Copying 2D Barcodes

abstract: This paper concerns the problem of detecting the use of information hiding at anti-copying 2D barcodes. Prior hidden information detection schemes are either heuristicbased or Machine Learning (ML) based. The key limitation of prior heuristics-based schemes is that they do not answer the fundamental question of why the information hidden at a 2D barcode can be detected. The key limitation of prior MLbased information schemes is that they lack robustness because a printed 2D barcode is very much environmentally dependent, and thus an information hiding detection scheme trained in one environment often does not work well in another environment. In this paper, we propose two hidden information detection schemes at the existing anti-copying 2D barcodes. The first scheme is to directly use the pixel distance to detect the use of an information hiding scheme in a 2D barcode, referred as to the Pixel Distance Based Detection (PDBD) scheme. The second scheme is first to calculate the variance of the raw signal and the covariance between the recovered signal and the raw signal, and then based on the variance results, detects the use of information hiding scheme in a 2D barcode, referred as to the Pixel Variance Based Detection (PVBD) scheme. Moreover, we design advanced IC attacks to evaluate the security of two existing anti-copying 2D barcodes. We implemented our schemes and conducted extensive performance comparison between our schemes and prior schemes under different capturing devices, such as a scanner and a camera phone. Our experimental results show that the PVBD scheme can correctly detect the existence of the hidden information at both the 2LQR code and the LCAC 2D barcode. Moreover, the probability of successfully attacking of our IC attacks achieves 0.6538 for the 2LQR code and 1 for the LCAC 2D barcode.


title: Data-Aided LS Channel Estimation in Massive MIMO Turbo-Receiver

abstract: In this paper, we propose a new algorithm of iterative least squared (LS) channel estimation for 64 antennas Massive Multiple Input, Multiple Output (MIMO) turbo-receiver. The algorithm employs log-likelihood ratios (LLR) of low-density parity-check (LDPC) decoder and minimum mean square error (MMSE) estimator to achieve soft data symbols. These soft data symbols are further MMSE-weighted again and combined with pilot symbols to achieve a modified LS channel estimate. The modified LS estimate is employed by the same channel estimation unit to enhance turbo-receiver performance via channel re-estimation, as a result, the proposed approach has low complexity and fits any channel estimation solution, which is quite valuable in practice. We analyze both hard and soft algorithm versions and present simulation results of 5G turbo-receiver in the 3D-UMa model of the QuaDRiGa 2.0 channel. Simulation results demonstrate up to 0.3dB performance gain compared to the unweighted hard data symbols utilization in the LS channel re-calculation.


title: Bayesian approach to inverse scattering with topological priors

abstract: We propose a Bayesian inference framework to estimate uncertainties in inverse scattering problems. Given the observed data, the forward model and their uncertainties, we find the posterior distribution over a finite parameter field representing the objects. To construct the prior distribution we use a topological sensitivity analysis. We demonstrate the approach on the Bayesian solution of 2D inverse problems in light and acoustic holography with synthetic data. Statistical information on objects such as their center location, diameter size, orientation, as well as material properties, are extracted by sampling the posterior distribution. Assuming the number of objects known, comparison of the results obtained by Markov Chain Monte Carlo sampling and by sampling a Gaussian distribution found by linearization about the maximum a posteriori estimate show reasonable agreement. The latter procedure has low computational cost, which makes it an interesting tool for uncertainty studies in 3D. However, MCMC sampling provides a more complete picture of the posterior distribution and yields multi-modal posterior distributions for problems with larger measurement noise. When the number of objects is unknown, we devise a stochastic model selection framework.

categories: math.NA cs.NA math.OC math.ST physics.data-an stat.TH


title: Compact Schur-Weyl duality: real Lie groups and the cyclotomic Brauer algebra

abstract: We show that the centraliser of the maximal compact subgroup of the real orthogonal or symplectic groups acting on tensors of their standard representation are isomorphic to cyclotomic Brauer algebras. We also show that for the symplectic group this cyclotomic Brauer algebra splits into summands of walled Brauer algebras.


title: The early phase of the COVID-19 outbreak in Lombardy, Italy

abstract: In the night of February 20, 2020, the first case of novel coronavirus disease (COVID-19) was confirmed in the Lombardy Region, Italy. In the week that followed, Lombardy experienced a very rapid increase in the number of cases. We analyzed the first 5,830 laboratory-confirmed cases to provide the first epidemiological characterization of a COVID-19 outbreak in a Western Country. Epidemiological data were collected through standardized interviews of confirmed cases and their close contacts. We collected demographic backgrounds, dates of symptom onset, clinical features, respiratory tract specimen results, hospitalization, contact tracing. We provide estimates of the reproduction number and serial interval. The epidemic in Italy started much earlier than February 20, 2020. At the time of detection of the first COVID-19 case, the epidemic had already spread in most municipalities of Southern-Lombardy. The median age for of cases is 69 years (range, 1 month to 101 years). 47% of positive subjects were hospitalized. Among these, 18% required intensive care. The mean serial interval is estimated to be 6.6 days (95% CI, 0.7 to 19). We estimate the basic reproduction number at 3.1 (95% CI, 2.9 to 3.2). We estimated a decreasing trend in the net reproduction number starting around February 20, 2020. We did not observe significantly different viral loads in nasal swabs between symptomatic and asymptomatic. The transmission potential of COVID-19 is very high and the number of critical cases may become largely unsustainable for the healthcare system in a very short-time horizon. We observed a slight decrease of the reproduction number, possibly connected with an increased population awareness and early effect of interventions. Aggressive containment strategies are required to control COVID-19 spread and catastrophic outcomes for the healthcare system.


title: Stability of the Inverse Problem for Dini Continuous Conductivities in the Plane

abstract: We show that the inverse problem of Calderon for conductivities in a two-dimensional Lipschitz domain is stable in a class of conductivities that are Dini continuous. This extends previous stability results when the conductivities are known to be Holder continuous.


title: Crime Prediction Using Spatio-Temporal Data

abstract: A crime is a punishable offence that is harmful for an individual and his society. It is obvious to comprehend the patterns of criminal activity to prevent them. Research can help society to prevent and solve crime activates. Study shows that only 10 percent offenders commits 50 percent of the total offences. The enforcement team can respond faster if they have early information and pre-knowledge about crime activities of the different points of a city. In this paper, supervised learning technique is used to predict crimes with better accuracy. The proposed system predicts crimes by analyzing data-set that contains records of previously committed crimes and their patterns. The system stands on two main algorithms - i) decision tree, and ii) k-nearest neighbor. Random Forest algorithm and Adaboost are used to increase the accuracy of the prediction. Finally, oversampling is used for better accuracy. The proposed system is feed with a criminal-activity data set of twelve years of San Francisco city.


title: Hadronic weak decays of $\Lambda_c$ in the quark model

abstract: The hadronic weak decays of $\Lambda_c$ are studied in the framework of a constituent quark model. With the combined analysis of the Cabbibo-favored processes, $\Lambda_c\to \Lambda\pi^+$, $\Sigma^0\pi^+$ and $\Sigma^+\pi^0$, we confirm that the non-factorizable transition mechanisms play a crucial role in the understanding of their compatible branching ratios. We emphasize that the SU(3) flavor symmetry breaking effects, which is generally at the order of $1\sim 2\%$, can be amplified by the destructive interferences among the pole terms in the diagrams with internal conversion. Some contributions are sensitive to the spatial distribution of the scalar-isoscalar light-quark sector in the $\Lambda_c$, and its overlap with the light quarks in the final state hyperon. Namely, a compact diquark configuration is disfavored.


title: A catalog of intermediate duration Type I X-ray bursts observed with the INTEGRAL satellite

abstract: We present a catalog of long duration bursts observed with the Joint European X-ray Monitor (JEM-X) and IBIS/ISGRI instruments onboard the INTEGRAL satellite. The fourteen bursts have e-folding times ranging from 55 s to ? 17 min, and are therefore classified as intermediate-duration bursts, caused by the ignition of an unusually thick helium layer. Though seven events have already been reported in literature, we have systematically re-analyzed the whole sample. We find three new photospheric radius expansion (PRE) bursts, which are not reported in the literature, allowing us to provide a new estimate of the distances to these sources. We apply the enhanced persistent emission method (also known as the fa method) on sources with detectable persistent emission prior to a burst, in order to follow the evolution of the accretion rate during the burst. Although we do not get significantly better fits, the evolution of the fa factor shows an indicative behavior, which we discuss.


title: The changing-type SN 2014C may come from an 11-$M_\odot$ star stripped by binary interaction and violent eruption

abstract: SN 2014C was an unprecedented supernova (SN) that displayed a metamorphosis from Type Ib to Type IIn over $\sim$200 days. This transformation is consistent with a helium star having exploded in a cavity surrounded by a dense shell of the progenitor's stripped hydrogen envelope. For at least 5 years post-explosion, the ejecta continued to interact with an outer, extended component of circumstellar medium (CSM) that was ejected even before the dense shell. It is still unclear, however, what kind of progenitor could have undergone such a complicated mass-loss history before it produced this peculiar SN. In this paper, we report a new analysis of SN 2014C's host star cluster based on data from the Hubble Space Telescope (HST). By carefully fitting its spectral energy distribution (SED), we derive a precise cluster age of 20.0$^{+3.5}_{-2.6}$ Myr, which corresponds to the progenitor's lifetime assuming coevolution. Combined with binary stellar evolution models, we find that SN 2014C's progenitor may have been an $\sim$11-$M_\odot$ star in a relatively wide binary system. The progenitor's envelope was partially stripped by Case C or Case BC mass transfer via binary interaction, followed by a violent eruption that ejected the last hydrogen layer before terminal explosion. Thus, SN 2014C, in common with SNe 2006jc and 2015G, may be a third example that violent eruptions, with mass-loss rates matching luminous blue variable (LBV) giant eruptions, can also occur in much lower-mass massive stars if their envelopes are partially or completely stripped in interacting binaries.


title: Towards Laser Ion Acceleration With Holed Targets

abstract: Although the interaction of a flat-foil with currently available laser intensities is now considered a routine process, during the last decade emphasis is given to targets with complex geometries aiming on increasing the ion energy. This work presents a target geometry where two symmetric side-holes and a central-hole are drilled into the foil. A study of the various side-holes and central-hole length combinations is performed with 2-dimensional particle-in-cell simulations for polyethylene targets and a laser intensity of 5.2x10^21 W cm^-2. The holed-targets show a remarkable increase of the conversion efficiency, which corresponds to a different target configuration for electrons, protons and carbon ions. Furthermore, diffraction of the laser pulse leads to a directional high energy electron beam, with a temperature of ~40 MeV or seven times higher than in the case of a flat-foil. The higher conversion efficiency consequently leads to a significant enhancement of the maximum proton energy from holed-targets.


title: Stance Control Inspired by Cerebellum Stabilizes Reflex-Based Locomotion on HyQ Robot

abstract: Advances in legged robotics are strongly rooted in animal observations. A clear illustration of this claim is the generalization of Central Pattern Generators (CPG), first identified in the cat spinal cord, to generate cyclic motion in robotic locomotion. Despite a global endorsement of this model, physiological and functional experiments in mammals have also indicated the presence of descending signals from the cerebellum, and reflex feedback from the lower limb sensory cells, that closely interact with CPGs. To this day, these interactions are not fully understood. In some studies, it was demonstrated that pure reflex-based locomotion in the absence of oscillatory signals could be achieved in realistic musculoskeletal simulation models or small compliant quadruped robots. At the same time, biological evidence has attested the functional role of the cerebellum for predictive control of balance and stance within mammals. In this paper, we promote both approaches and successfully apply reflex-based dynamic locomotion, coupled with a balance and gravity compensation mechanism, on the state-of-art HyQ robot. We discuss the importance of this stability module to ensure a correct foot lift-off and maintain a reliable gait. The robotic platform is further used to test two different architectural hypotheses inspired by the cerebellum. An analysis of experimental results demonstrates that the most biologically plausible alternative also leads to better results for robust locomotion.


title: Flexible placements of graphs with rotational symmetry

abstract: We study the existence of an $n$-fold rotationally symmetric placement of a symmetric graph in the plane allowing a continuous deformation that preserves the symmetry and the distances between adjacent vertices. We show that such a flexible placement exists if and only if the graph has a NAC-colouring satisfying an additional property on the symmetry; a NAC-colouring is a surjective edge colouring by two colours such that every cycle is either monochromatic, or there are at least two edges of each colour.


title: Berger And Turaev's Proof of Herman's Conjecture of Positive Entropy

abstract: We present the proof of Berger and Turaev of Herman's positive entropy conjecture. In every neighbourhood of identity in the set of smooth symplectic diffeomorphisms of the 2-dimensional disc, there exists a diffeomorphism with positive metric entropy.


title: Understanding the Schrodinger equation as a kinematic statement: A probability-first approach to quantum

abstract: Quantum technology is seeing a remarkable explosion in interest due to a wave of successful commercial technology. As a wider array of engineers and scientists are needed, it is time we rethink quantum educational paradigms. Current approaches often start from classical physics, linear algebra, or differential equations. This chapter advocates for beginning with probability theory. In the approach outlined in this chapter, there is less in the way of explicit axioms of quantum mechanics. Instead the historically problematic measurement axiom is inherited from probability theory where many philosophical debates remain. Although not a typical route in introductory material, this route is nonetheless a standard vantage on quantum mechanics. This chapter outlines an elementary route to arrive at the Schr\"odinger equation by considering allowable transformations of quantum probability functions (density matrices). The central tenet of this chapter is that probability theory provides the best conceptual and mathematical foundations for introducing the quantum sciences.


title: Characterizing gonality for two-component stable curves

abstract: It is a well-known result that a stable curve of compact type over $\mathbb{C}$ having two components is hyperelliptic if and only if both components are hyperelliptic and the point of intersection is a Weierstrass point for each of them. With the use of admissible covers, we generalize this characterization in two ways: for stable curves of higher gonality having two smooth components and one node; and for hyperelliptic and trigonal stable curves having two smooth non rational components and any number of nodes.


title: A Generalized Euler Probability Distribution

abstract: From a new class of q-deformed coherent states we introduce a generalization of the Euler probability distribution for which the main statistical parameters are obtained explicitly. As application, we discuss the corresponding photon counting statistics with respect to all parameters labeling the coherent states under consideration.


title: Physiologically Driven Storytelling: Concept and Software Tool

abstract: We put forth Physiologically Driven Storytelling, a new approach to interactive storytelling where narratives adaptively unfold based on the reader's physiological state. We first describe a taxonomy framing how physiological signals can be used to drive interactive systems both as input and output. We then propose applications to interactive storytelling and describe the implementation of a software tool to create Physiological Interactive Fiction (PIF). The results of an online study (N=140) provided guidelines towards augmenting the reading experience. PIF was then evaluated in a lab study (N=14) to determine how physiological signals can be used to infer a reader's state. Our results show that breathing, electrodermal activity, and eye tracking can help differentiate positive from negative tones, and monotonous from exciting events. This work demonstrates how PIF can support storytelling in creating engaging content and experience tailored to the reader. Moreover, it opens the space to future physiologically driven systems within broader application areas.


title: Memory effect and BMS symmetries for extreme black holes

abstract: We study horizon shells and soldering freedom for extreme black holes and how supertranslation-like Bondi-Metzner-Sachs (BMS) symmetries appear as soldering transformations. Further, for a null shell placed infinitesimally close to the horizon of an extreme Reissner-Nordstr$\ddot{o}$m (RN) black hole, we show superrotation-like symmetries also arise as soldering freedom. Next, considering the interaction of impulsive gravitational waves supported at the horizon shell with test particles, we study how the "memory" (or the imprints) of BMS-like symmetries gets encoded in the geodesics (test particles) crossing the shell. Our study shows, timelike test particles get displaced from their initial plane when they cross the horizon shell. For a null geodesic congruence crossing the horizon shell, the optical tensors corresponding to the congruence suffer jumps. In both the cases, the changes are induced by BMS parameters that constitute the gravity wave and matter degrees of freedom of the shell.


title: Fast generalized Nash equilibrium seeking under partial-decision information

abstract: We address the generalized Nash equilibrium seeking problem in a partial-decision information scenario, where each agent can only exchange information with some neighbors, although its cost function possibly depends on the strategies of all agents. The few existing methods build on projected pseudo-gradient dynamics, and require either double-layer iterations or conservative conditions on the step sizes. To overcome both these flaws and improve efficiency, we design the first fully-distributed single-layer algorithms based on proximal best-response. Our schemes are fixed-step and allow for inexact updates, which is crucial for reducing the computational complexity. Under standard assumptions on the game primitives, we establish convergence to a variational equilibrium (with linear rate for games without coupling constraints) by recasting our algorithms as proximal-point methods, opportunely preconditioned to distribute the computation among the agents. Since our analysis hinges on a restricted monotonicity property, we also provide new general results that significantly extend the domain of applicability of proximal-point methods. Besides, the operator-theoretic approach favors the implementation of provably correct acceleration schemes that can further improve the convergence speed. Finally, the potential of our algorithms is demonstrated numerically, revealing much faster convergence with respect to projected pseudo-gradient methods and validating our theoretical findings.


title: Forecasting the Ambient Solar Wind with Numerical Models. II. An Adaptive Prediction System for Specifying Solar Wind Speed Near the Sun

abstract: The ambient solar wind flows and fields influence the complex propagation dynamics of coronal mass ejections in the interplanetary medium and play an essential role in shaping Earth's space weather environment. A critical scientific goal in the space weather research and prediction community is to develop, implement and optimize numerical models for specifying the large-scale properties of solar wind conditions at the inner boundary of the heliospheric model domain. Here we present an adaptive prediction system that fuses information from in situ measurements of the solar wind into numerical models to better match the global solar wind model solutions near the Sun with prevailing physical conditions in the vicinity of Earth. In this way, we attempt to advance the predictive capabilities of well-established solar wind models for specifying solar wind speed, including the Wang-Sheeley-Arge (WSA) model. In particular, we use the Heliospheric Upwind eXtrapolation (HUX) model for mapping the solar wind solutions from the near-Sun environment to the vicinity of Earth. In addition, we present the newly developed Tunable HUX (THUX) model which solves the viscous form of the underlying Burgers equation. We perform a statistical analysis of the resulting solar wind predictions for the time 2006-2015. The proposed prediction scheme improves all the investigated coronal/heliospheric model combinations and produces better estimates of the solar wind state at Earth than our reference baseline model. We discuss why this is the case, and conclude that our findings have important implications for future practice in applied space weather research and prediction.


title: Non-homogeneous initial boundary value problems for the biharmonic Schr\"odinger equation on an interval

abstract: In this paper we consider the initial boundary value problem (IBVP) for the nonlinear biharmonic Schr\"odinger equation posed on a bounded interval $(0,L)$ with non-homogeneous Navier or Dirichlet boundary conditions, respectively. For Navier boundary IBVP, we set up its local well-posedness if the initial data lies in $H^s(0, L)$ with $s\geq 0$ and $s\neq n+1/2, n\in \mathbb{N}$, and the boundary data are selected from the appropriate spaces with optimal regularities, i.e., the $j$-th order data are chosen in $H_{loc}^{(s+3-j)/4}(\mathbb {R}^+)$, for $j=0,2$. For Dirichlet boundary IBVP the corresponding local well-posedness is obtained when $s>10/7$ and $s\neq n+1/2, n\in \mathbb{N}$, and the boundary data are selected from the appropriate spaces with optimal regularities, i.e., the $j$-th order data are chosen in $H_{loc}^{(s+3-j)/4}(\mathbb {R}^+)$, for $j=0,1$.


title: Selecting Relevant Features from a Multi-domain Representation for Few-shot Classification

abstract: Popular approaches for few-shot classification consist of first learning a generic data representation based on a large annotated dataset, before adapting the representation to new classes given only a few labeled samples. In this work, we propose a new strategy based on feature selection, which is both simpler and more effective than previous feature adaptation approaches. First, we obtain a multi-domain representation by training a set of semantically different feature extractors. Then, given a few-shot learning task, we use our multi-domain feature bank to automatically select the most relevant representations. We show that a simple non-parametric classifier built on top of such features produces high accuracy and generalizes to domains never seen during training, which leads to state-of-the-art results on MetaDataset and improved accuracy on mini-ImageNet.


title: On a sharp lemma of Cassels and Montgomery on manifolds

abstract: Let $\left( \mathcal{M},g\right) $ be a $d$-dimensional compact connected Riemannian manifold and let $\left\{ \varphi_{m}\right\}_{m=0}^{+\infty}$ be a complete sequence of orthonormal eigenfunctions of the Laplace-Beltrami operator on $\mathcal{M}$. We show that there exists a positive constant $C$ such that for all integers $N$ and $X$ and for all finite sequences of $N$ points in $\mathcal{M}$, $\left\{ x\left( j\right) \right\}_{j=1}^{N}$, and positive weights $\left\{ a_{j}\right\}_{j=1}^{N}$ we have \[ \sum_{m=0}^{X} | \sum_{j=1}^{N} a_{j} \varphi_{m} ( x( j) ) | ^{2}\geq \max \{ CX\sum_{j=1}^{N}a_{j}^{2},( \sum_{j=1}^{N}a_{j}) ^{2}\}.\]


title: Ordered Functional Decision Diagrams: A Functional Semantics For Binary Decision Diagrams

abstract: We introduce a novel framework, termed $\lambda$DD, that revisits Binary Decision Diagrams from a purely functional point of view. The framework allows to classify the already existing variants, including the most recent ones like Chain-DD and ESRBDD, as implementations of a special class of ordered models. We enumerate, in a principled way, all the models of this class and isolate its most expressive model. This new model, termed $\lambda$DD-O-NUCX, is suitable for both dense and sparse Boolean functions, and is moreover invariant by negation. The canonicity of $\lambda$DD-O-NUCX is formally verified using the Coq proof assistant. We furthermore give bounds on the size of the different diagrams: the potential gain achieved by more expressive models can be at most linear in the number of variables n.


title: High-redshift cosmography: auxiliary variables versus Pad\'e polynomials

abstract: Cosmography becomes non-predictive when cosmic data span beyond the red shift limit $z\simeq1 $. This leads to a \emph{strong convergence issue} that jeopardizes its viability. In this work, we critically compare the two main solutions of the convergence problem, i.e. the $y$-parametrizations of the redshift and the alternatives to Taylor expansions based on Pad\'e series. In particular, among several possibilities, we consider two widely adopted parametrizations, namely $y_1=1-a$ and $y_2=\arctan(a^{-1}-1)$, being $a$ the scale factor of the Universe. We find that the $y_2$-parametrization performs relatively better than the $y_1$-parametrization over the whole redshift domain. Even though $y_2$ overcomes the issues of $y_1$, we get that the most viable approximations of the luminosity distance $d_L(z)$ are given in terms of Pad\'e approximations. In order to check this result by means of cosmic data, we analyze the Pad\'e approximations up to the fifth order, and compare these series with the corresponding $y$-variables of the same orders. We investigate two distinct domains involving Monte Carlo analysis on the Pantheon Superovae Ia data, $H(z)$ and shift parameter measurements. We conclude that the (2,1) Pad\'e approximation is statistically the optimal approach to explain low and high-redshift data, together with the fifth-order $y_2$-parametrization. At high redshifts, the (3,2) Pad\'e approximation cannot be fully excluded, while the (2,2) Pad\'e one is essentially ruled out.


title: Modeling tumor growth: a simple individual-based model and its analysis

abstract: Initiation and development of a malignant tumor is a complex phenomenon that has critical stages determining its long time behavior. This phenomenon is mathematically described by means of various models: from simple heuristic models to those employing stochastic processes. In this chapter, we discuss some aspects of such modeling by analyzing a simple individual-based model, in which tumor cells are presented as point particles drifting in $\mathbf{R}_{+}:=[0,+\infty)$ towards the origin with unit speed. At the origin, each of them splits into two new particles that instantly appear in $\mathbf{R}_{+}$ at random positions. During their drift the particles are subject to a random death before splitting. In this model, trait $x\in \mathbf{R}_{+}$ of a given cell corresponds to time to its division and the death is caused by therapeutic factors. On its base we demonstrate how to derive a condition -- involving the therapy related death rate and cell cycle distribution parameters -- under which the tumor size remains bounded in time, which practically means combating the disease.


title: Ising model on random triangulations of the disk: phase transition

abstract: In [arXiv:1806.06668], we have studied the Boltzmann random triangulation of the disk coupled to an Ising model on its faces with Dobrushin boundary condition at its critical temperature. In this paper, we investigate the phase transition of this model by extending our previous results to arbitrary temperature: We compute the partition function of the model at all temperatures, and derive several critical exponents associated with the infinite perimeter limit. We show that the model has a local limit at any temperature, whose properties depend drastically on the temperature. At high temperatures, the local limit is reminiscent of the uniform infinite half-planar triangulation (UIHPT) decorated with a subcritical percolation. At low temperatures, the local limit develops a bottleneck of finite width due to the energy cost of the main Ising interface between the two spin clusters imposed by the Dobrushin boundary condition. This change can be summarized by a novel order parameter with a nice geometric meaning. In addition to the phase transition, we also generalize our construction of the local limit from the two-step asymptotic regime used in [arXiv:1806.06668] to a more natural diagonal asymptotic regime. We obtain in this regime a scaling limit related to the length of the main Ising interface, which coincides with predictions from the continuum theory of quantum surfaces (a.k.a.\ Liouville quantum gravity).


title: M\"obius Functions of Some Annular Noncrossing Objects

abstract: We present the M\"{o}bius functions of several posets of annular noncrossing objects, namely a self-dual extension of the annular noncrossing permutations, minimal length annular partitioned permutations, and annular noncrossing partitons.


title: Entropy rigidity for 3D conservative Anosov flows and dispersing billiards

abstract: Given an integer $k \geq 5$, and a $C^k$ Anosov flow $\Phi$ on some compact connected $3$-manifold preserving a smooth volume, we show that the measure of maximal entropy (MME) is the volume measure if and only if $\Phi$ is $C^{k-\varepsilon}$-conjugate to an algebraic flow, for $\varepsilon>0$ arbitrarily small. Besides the rigidity, we also study the entropy flexibility, and show that the metric entropy with respect to the volume measure and the topological entropy of suspension flows over Anosov diffeomorphisms on the $2$-torus achieve all possible values subject to natural normalizations. Moreover, in the case of dispersing billiards, we show that if the measure of maximal entropy is the volume measure, then the Birkhoff Normal Form of regular periodic orbits with a homoclinic intersection is linear.


title: Dependence of Solar Wind Proton Temperature on the Polarisation Properties of Alfv\'enic Fluctuations at Ion-kinetic Scales

abstract: We use fluctuating magnetic helicity to investigate the polarisation properties of Alfv\'enic fluctuations at ion-kinetic scales in the solar wind as a function of $\beta_p$, the ratio of proton thermal pressure to magnetic pressure, and $\theta_{vB}$, the angle between the proton flow and local mean magnetic field, $\mathbf{B}_0$. Using almost 15 years of \textit{Wind} observations, we separate the contributions to helicity from fluctuations with wave-vectors, $\textbf{k}$, quasi-parallel and oblique to $\mathbf{B}_0$, finding that the helicity of Alfv\'enic fluctuations is consistent with predictions from linear Vlasov theory. This result suggests that the non-linear turbulent fluctuations at these scales share at least some polarisation properties with Alfv\'en waves. We also investigate the dependence of proton temperature in the $\beta_p$-$\theta_{vB}$ plane to probe for possible signatures of turbulent dissipation, finding that it correlates with $\theta_{vB}$. The proton temperature parallel to $\mathbf{B}_0$ is higher in the parameter space where we measure the helicity of right-handed Alfv\'enic fluctuations, and the temperature perpendicular to $\mathbf{B}_0$ is higher where we measure left-handed fluctuations. This finding is inconsistent with the general assumption that by sampling different $\theta_{vB}$ in the solar wind we can analyse the dependence of the turbulence distribution on $\theta_{kB}$, the angle between $\textbf{k}$ and $\mathbf{B}_0$. After ruling out both instrumental and expansion effects, we conclude that our results provide new evidence for the importance of local kinetic processes that depend on $\theta_{vB}$ in determining proton temperature in the solar wind.


title: SAT: Improving Adversarial Training via Curriculum-Based Loss Smoothing

abstract: Adversarial training (AT) has become a popular choice for training robust networks. However, it tends to sacrifice clean accuracy heavily in favor of robustness and suffers from a large generalization error. To address these concerns, we propose Smooth Adversarial Training (SAT), guided by our analysis on the eigenspectrum of the loss Hessian. We find that curriculum learning, a scheme that emphasizes on starting "easy" and gradually ramping up on the "difficulty" of training, smooths the adversarial loss landscape for a suitably chosen difficulty metric. We present a general formulation for curriculum learning in the adversarial setting and propose two difficulty metrics based on the maximal Hessian eigenvalue (H-SAT) and the softmax probability (P-SA). We demonstrate that SAT stabilizes network training even for a large perturbation norm and allows the network to operate at a better clean accuracy versus robustness trade-off curve compared to AT. This leads to a significant improvement in both clean accuracy and robustness compared to AT, TRADES, and other baselines. To highlight a few results, our best model improves normal and robust accuracy by 6% and 1% on CIFAR-100 compared to AT, respectively. On Imagenette, a ten-class subset of ImageNet, our model outperforms AT by 23% and 3% on normal and robust accuracy respectively.


title: Predictive modelling of disease propagation in a mobile, connected community

abstract: We present numerical results obtained from the modelling of a stochastic, highly connected and mobile community. The spread of attributes like health, disease among the community members is simulated using cellular automata on a planar 2 dimensional bounded space. With remarkably few assumptions, we are able to predict the future course of propagation of such disease as a function of time and the fine-tuning of parameters related to the interaction among the automata.


title: Studies on the Spectral Schwartz Distribution

abstract: The resolvent function of an operator in a Banach space is defined on an open subset of the complex plane and is holomorphic. It obeys the resolvent equation. A generalization of this equation to Schwartz distributions is defined and a Schwartz distribution, which satisfies that equation is called a resolvent distribution. In important cases the resolvent distribution is the whole plane. Its restriction to the subset, where it is continuous, is the usual resolvent function. Its complex conjugate derivative is ,but a factor, the spectral Schwartz distribution, which is carried by a subset of the spectral set of the operator. The spectral distribution yields a spectral decomposition. We have a generalized orthogonality relation. Completeness is defined in a natural way and is the case e.g. if the operator is bounded. The spectral distribution of a matrix and a unitary operator are given. If the the operator is a self-adjoint operator on a Hilbert space, the spectral distribution is the derivative of the spectral family. We calculate the spectral distribution and the eigen value problem of the multiplication operator and of the multiplication operator with a non symmetric rank one perturbation. The operator is not normal and may have discrete real or imaginary eigenvalues or a nontrivial Jordan decomposition.


title: Zeno Dynamics of Quantum Chirps

abstract: We investigate the Zeno dynamics of quantum chirps. More specifically, we analyze the Zeno dynamics of a chirped solitary wave solution of the non-dimensional nonlinear Schr\"odinger equation (NLSE) with a gain/loss term. We show that the motion of the quantum chirp can be inhibited by frequent measurements in the observation domain. In order to assess the effect of Zeno dynamics on the quantum chirp imaging resolution, we analyze and compare its point spread functions (PSF) under free evolution and under Zeno dynamics. Additionally, we investigate the effect of observation frequency on the quantum chirp profile as well as on its probability of lingering in the observation domain. Our results can be immediately used to investigate the signal properties of quantum chirps and their possible jamming properties by Zeno dynamics in quantum radar technology. Additionally, many potential applications in atomic physics and optics including nonlinear phenomena and chirps can be investigated using the approach presented in our paper.


title: Chaotic diffusion for particles moving in a time dependent potential well

abstract: The chaotic diffusion for particles moving in a time dependent potential well is described by using two different procedures: (i) via direct evolution of the mapping describing the dynamics and ; (ii) by the solution of the diffusion equation. The dynamic of the diffusing particles is made by the use of a two dimensional, nonlinear area preserving map for the variables energy and time. The phase space of the system is mixed containing both chaos, periodic regions and invariant spanning curves limiting the diffusion of the chaotic particles. The chaotic evolution for an ensemble of particles is treated as random particles motion and hence described by the diffusion equation. The boundary conditions impose that the particles can not cross the invariant spanning curves, serving as upper boundary for the diffusion, nor the lowest energy domain that is the energy the particles escape from the time moving potential well. The diffusion coefficient is determined via the equation of the mapping while the analytical solution of the diffusion equation gives the probability to find a given particle with a certain energy at a specific time. The momenta of the probability describe qualitatively the behavior of the average energy obtained by numerical simulation, which is investigated either as a function of the time as well as some of the control parameters of the problem.


title: From identification of random contributions to determination of the optimum forecast of a soccer match

abstract: The forecasting of sports events is of broad interest from the applied but also from the theoretical perspective. In this work the question is addressed for the example of the German soccer Bundesliga how a theoretically optimum forecast of the goal difference of a match can be characterized. This involves a careful analysis of the random contributions in a match and its disentanglement from the informative contributions, resulting from the individual team strengths. An important aspect is the consideration of the time dependence of the team strength which turns out to mainly fluctuate around a team-specific value during the course of a season. Two types of time-dependent properties have to be distinguished, one being uncorrelated between different match-days, the other being correlated and thus accessible by an appropriate correlation analysis. For some performance indicators, which may be used to estimate the team strength, the quality of the respective forecast is compared to the theoretical optimum. Knowledge of the informative contribution allows one to conclude that the offensive team strength is more important than the defensive team strength for the final success.


title: Scattering-induced intensity reduction: large mass content with small grains in the inner region of the TW Hya disk

abstract: Dust continuum observation is one of the best methods to constrain the properties of protoplanetary disks. Recent theoretical studies have suggested that the dust scattering at the millimeter wavelength potentially reduces the observed intensity, which results in an underestimate in the dust mass. We investigate whether the dust scattering indeed reduces the observed continuum intensity by comparing the ALMA archival data of the TW Hya disk at Band 3, 4, 6, 7 and 9 to models obtained by radiative transfer simulations. We find that the model with scattering by 300 ${\rm \mu m}$-sized grains well reproduces the observed SED of the central part of the TW Hya disk while the model without scattering is also consistent within the errors of the absolute fluxes. To explain the intensity at Band 3, the dust surface density needs to be $\sim$ 10 ${\rm g\,cm^{-2}}$ at 10 au in the model with scattering, which is 26 times more massive than previously predicted. The model without scattering needs 2.3 times higher dust mass than the model with scattering because it needs lower temperature. At Band 7, scattering reduces the intensity by $\sim$ 35% which makes the disk looks optically thin even though it is optically thick. Our study suggests the TW Hya disk is still capable of forming cores of giant planets at where the current solar system planets exist.


title: Visual Navigation Among Humans with Optimal Control as a Supervisor

abstract: Real world visual navigation requires robots to operate in unfamiliar, human-occupied dynamic environments. Navigation around humans is especially difficult because it requires anticipating their future motion, which can be quite challenging. We propose an approach that combines learning-based perception with model-based optimal control to navigate among humans based only on monocular, first-person RGB images. Our approach is enabled by our novel data-generation tool, HumANav that allows for photorealistic renderings of indoor environment scenes with humans in them, which are then used to train the perception module entirely in simulation. Through simulations and experiments on a mobile robot, we demonstrate that the learned navigation policies can anticipate and react to humans without explicitly predicting future human motion, generalize to previously unseen environments and human behaviors, and transfer directly from simulation to reality. Videos describing our approach and experiments, as well as a demo of HumANav are available on the project website.


title: Storing and retrieving wavefronts with resistive temporal memory

abstract: We extend the reach of temporal computing schemes by developing a memory for multi-channel temporal patterns or "wavefronts." This temporal memory re-purposes conventional one-transistor-one-resistor (1T1R) memristor crossbars for use in an arrival-time coded, single-event-per-wire temporal computing environment. The memristor resistances and the associated circuit capacitances provide the necessary time constants, enabling the memory array to store and retrieve wavefronts. The retrieval operation of such a memory is naturally in the temporal domain and the resulting wavefronts can be used to trigger time-domain computations. While recording the wavefronts can be done using standard digital techniques, that approach has substantial translation costs between temporal and digital domains. To avoid these costs, we propose a spike timing dependent plasticity (STDP) inspired wavefront recording scheme to capture incoming wavefronts. We simulate these designs with experimentally validated memristor models and analyze the effects of memristor non-idealities on the operation of such a memory.


title: On the affinization of a nilpotent orbit cover

abstract: Let $\mathfrak{g}$ be a simple classical Lie algebra over $\mathbb{C}$ and $G$ be the adjoint group. Consider a nilpotent element $e\in \mathfrak{g}$, and the adjoint orbit $\mathbb{O}=Ge$. The formal slices to the codimension $2$ orbits in the closure $\overline{\mathbb{O}}\subset \mathfrak{g}$ are well-known due to the work of Kraft and Procesi. In this paper, we prove a similar result for the universal $G$-equivariant cover $\widetilde{\mathbb{O}}$ of $\mathbb{O}$. Namely, we describe the codimension $2$ singularities for its affinization $Spec(\mathbb{C}[\widetilde{\mathbb{O}}])$.


title: Intermediate-mass Early-type Disk Galaxies in the Virgo Cluster. II. Near-Infrared Spectra and Evidence for Differences in Evolution

abstract: We discuss near-infrared (NIR) slit spectra of six early-type disk galaxies in the Virgo Cluster. While the depth of the CO(2,0) band is consistent with the centers of these galaxies having solar or slightly sub-solar luminosity-weighted metallicities, galaxy-to-galaxy variations in the depth of the NaI 2.21um doublet are found, with NaI2.21um lines in three galaxies (NGC 4491, NGC 4584, and NGC 4620) that are deeper than those predicted for a solar chemical mixture and a solar-neighborhood mass function. Considering the depth of the NaD lines found by Davidge (2018a), the deep NaI2.21um lines are tentatively attributed to a bottom-heavy mass function. This is counter to observed trends between mass function slope and velocity dispersion. Two of the three galaxies with deep NaI2.21um (NGC 4584 and NGC 4620) have Sersic exponents that are consistent with a classical bulge. As for NGC 4491, its central NIR spectrum contains prominent emission lines. The relative strengths of Brgamma and H_2S(1), the presence of [FeII] emission, and the mid-infrared spectral-energy distribution are all consistent with intense recent star formation near the center of that galaxy. Wavelengths that contain the Ballick-Ramsey C_2 band at 1.76um are matched by moderately metal-poor E-MILES model spectra with an age of 2 Gyr. It is concluded that these galaxies have evolved in a cluster environment for at least several Gyr.


title: On the asymptotic stability of the sine-Gordon kink in the energy space

abstract: We consider the sine-Gordon (SG) equation in 1+1 dimensions. The kink is a static, non-symmetric exact solution to SG, stable in the energy space $H^1\times L^2$. It is well-known that the linearized operator around the kink has a simple kernel and no internal modes. However, it possesses an odd resonance at the bottom of the continuum spectrum, deeply related to the existence of the (in)famous wobbling kink, an explicit periodic-in-time solution of SG around the kink that contradicts the asymptotic stability of the kink in the energy space. In this paper, we further investigate the influence of resonances in the asymptotic stability question. We also discuss the relationship between breathers, wobbling kinks and resonances in the SG setting. By gathering B\"acklund transformations (BT) and Virial estimates around odd perturbations of the vacuum solution, we first identify the manifold of initial data around zero under which BTs are related to the wobbling kink solution. It turns out that (even) small breathers are deeply related to odd perturbations around the kink, including the wobbling kink itself. As a consequence of this result, using BTs we can construct a smooth manifold of initial data close to the kink, for which there is asymptotic stability in the energy space. The initial data has spatial symmetry of the form (kink + odd, even), non-resonant in principle, and not preserved by the flow. This asymptotic stability property holds despite the existence of wobbling kinks in SG. We also show that wobbling kinks are orbitally stable under odd data, and clarify some interesting connections between SG and $\phi^4$ at the level of linear B\"acklund transformations.


title: ASAS J174406+2446.8 is identified as a marginal-contact binary with a possible cool third body

abstract: ASAS J174406+2446.8 was originally found as a $\delta$ Scuti-type pulsating star with the period P=0.189068 $days$ by ASAS survey. However, the LAMOST stellar parameters reveal that it is far beyond the red edge of pulsational instability strip on the $\log g-T$ diagram of $\delta$ Scuti pulsating stars. To understand the physical properties of the variable star, we observed it by the 1.0-m Cassegrain reflecting telescope at Yunnan Observatories. Multi-color light curves in B, V, R$_{c}$ and I$_{c}$ bands were obtained and are analyzed by using the W-D program. It is found that this variable star is a shallow-contact binary with an EB-type light curve and an orbital period of 0.3781\,days rather than a $\delta$ Scuti star. It is a W-subtype contact binary with a mass ratio of $1.135(\pm0.019)$ and a fill-out factor of $10.4(\pm5.6)\,\%$. The situation of ASAS J174406+2446.8 resembles those of other EB-type marginal-contact binaries such as UU Lyn, II Per and GW Tau. All of them are at a key evolutionary phase from a semi-detached configuration to a contact system predicted by the thermal relaxation oscillation theory. The linear ephemeris was corrected by using 303 new determined times of light minimum. It is detected that the O - C curve shows a sinusoidal variationthat could be explained by the light-travel-time effect via the presence of a cool red dwarf. The present investigation reveals that some of the $\delta$ Scuti-type stars beyond the red edge of pulsating instability strip on the $\log g-T$ diagram are misclassified eclipsing binaries. To understand their structures and evolutionary states, more studies are required in the future.


title: Efficient Modelling of Ion Structure and Dynamics in Inorganic Metal Halide Perovskites

abstract: Metal halide perovskites (MHPs) are nowadays one of the most studied semiconductors due to their exceptional performance as active layers in solar cells. Although MHPs are excellent solid-state semiconductors, they are also ionic compounds, where ion migration plays a decisive role in their formation, their photovoltaic performance and their long-term stability. Given the above-mentioned complexity, molecular dynamics simulations based on classical force fields are especially suited to study MHP properties, such as lattice dynamics and ion migration. In particular, the possibility to model mixed compositions is important since they are the most relevant to optimize the optical band gap and the stability. With this intention, we employ DFT calculations and a genetic algorithm to develop a fully transferable classical force field valid for the benchmark inorganic perovskite compositional set CsPb(Br_xI_(1-x))_3 (x = 0,1/3,2/3,1). The resulting force field reproduces correctly, with a common set of parameters valid for all compositions, the experimental lattice parameter as a function of bromide/iodide ratio, the ion-ion distances and the XRD spectra of the pure and mixed structures. The simulated thermal conductivities and ion migration activation energies of the pure compounds are also in good agreement with experimental trends. Our molecular dynamics simulations make it possible to predict the compositional dependence of the ionic diffusion coefficient on bromide/iodide ratio and vacancy concentration. For vacancy concentrations of around 9 10^21 cm^-3, we obtained ionic diffusion coefficients at ambient temperature of 10^-11 and 10^-13 cm2/s for CsPbBr3 and CsPbI3, respectively. Interestingly, in comparison with the pure compounds, we predict a significantly lower activation energy for vacancy migration and faster diffusion for the mixed perovskites.


title: Traffic Abstractions of Nonlinear Homogeneous Event-Triggered Control Systems

abstract: In previous work, linear time-invariant event-triggered control (ETC) systems were abstracted to finite-state systems that capture the original systems' sampling behaviour. It was shown that these abstractions can be employed for scheduling of communication traffic in networks of ETC loops. In this paper, we extend this framework to the class of nonlinear homogeneous systems, however adopting a different approach in a number of steps. Finally, we discuss how the proposed methodology could be extended to general nonlinear systems.


title: Uniform Error Estimates for the Lanczos Method

abstract: The Lanczos method is one of the most powerful and fundamental techniques for solving an extremal symmetric eigenvalue problem. Convergence-based error estimates depend heavily on the eigenvalue gap. In practice, this gap is often relatively small, resulting in significant overestimates of error. One way to avoid this issue is through the use of uniform error estimates, namely, bounds that depend only on the dimension of the matrix and the number of iterations. In this work, we prove explicit upper and lower uniform error estimates for the Lanczos method. These lower bounds, paired with numerical results, imply that the maximum error of $m$ iterations of the Lanczos method over all $n \times n$ symmetric matrices does indeed depend on the dimension $n$ in practice. The improved bounds for extremal eigenvalues translate immediately to error estimates for the condition number of a symmetric positive definite matrix. In addition, we prove more specific results for matrices that possess some level of eigenvalue regularity or whose eigenvalues converge to some limiting empirical spectral distribution. Through numerical experiments, we show that the theoretical estimates of this paper do apply to practical computations for reasonably sized matrices.


title: Efficiency Guarantees for Parallel Incremental Algorithms under Relaxed Schedulers

abstract: Several classic problems in graph processing and computational geometry are solved via incremental algorithms, which split computation into a series of small tasks acting on shared state, which gets updated progressively. While the sequential variant of such algorithms usually specifies a fixed (but sometimes random) order in which the tasks should be performed, a standard approach to parallelizing such algorithms is to relax this constraint to allow for out-of-order parallel execution. This is the case for parallel implementations of Dijkstra's single-source shortest-paths algorithm (SSSP), and for parallel Delaunay mesh triangulation. While many software frameworks parallelize incremental computation in this way, it is still not well understood whether this relaxed ordering approach can still provide any complexity guarantees. In this paper, we address this problem, and analyze the efficiency guarantees provided by a range of incremental algorithms when parallelized via relaxed schedulers. We show that, for algorithms such as Delaunay mesh triangulation and sorting by insertion, schedulers with a maximum relaxation factor of $k$ in terms of the maximum priority inversion allowed will introduce a maximum amount of wasted work of $O(log(n) poly (k) ), $ where $n$ is the number of tasks to be executed. For SSSP, we show that the additional work is $O(poly (k) d_{max} / w_{min}), $ where $d_{\max}$ is the maximum distance between two nodes, and $w_{min}$ is the minimum such distance. In practical settings where $n \gg k$, this suggests that the overheads of relaxation will be outweighed by the improved scalability of the relaxed scheduler. On the negative side, we provide lower bounds showing that certain algorithms will inherently incur a non-trivial amount of wasted work due to scheduler relaxation, even for relatively benign relaxed schedulers.


title: Composition and Weight Pushing of Monotonic Subsequential Failure Transducers Representing Probabilistic Models

abstract: We present a construction for the composition of subsequential transducers (representing conditional probabilistic models) with subsequential failure transducers (representing probabilistic models). Under certain conditions, satisfied by the corresponding transduction devices, a more efficient construction is applicable that avoids the creation of unnecessary states. Furthermore, the weights of the resulting failure transducers can be efficiently redistributed via weight pushing in the $\langle \mathbb{R}_+, +, \times, 0, 1 \rangle$ and $\langle \mathbb{R}_+, \max, \times, 0, 1 \rangle$ semirings.


title: Out-of-Distribution Detection for Skin Lesion Images with Deep Isolation Forest

abstract: In this paper, we study the problem of out-of-distribution detection in skin disease images. Publicly available medical datasets normally have a limited number of lesion classes (e.g. HAM10000 has 8 lesion classes). However, there exists a few thousands of clinically identified diseases. Hence, it is important if lesions not in the training data can be differentiated. Toward this goal, we propose DeepIF, a non-parametric Isolation Forest based approach combined with deep convolutional networks. We conduct comprehensive experiments to compare our DeepIF with three baseline models. Results demonstrate state-of-the-art performance of our proposed approach on the task of detecting abnormal skin lesions.


title: Quintessence from higher curvature supergravity

abstract: In this contribution we revisit higher curvature N=1 supergravity and discuss the quintessence phase that can appear due to the $R^4$ terms. In particular we focus on the bosonic supersymmetric completion within the old-minimal and the new-minimal formulations.


title: A Correlated Random Coefficient Panel Model with Time-Varying Endogeneity

abstract: This paper studies a class of linear panel models with random coefficients. We do not restrict the joint distribution of the time-invariant unobserved heterogeneity and the covariates. We investigate identification of the average partial effect (APE) when fixed-effect techniques cannot be used to control for the correlation between the regressors and the time-varying disturbances. Relying on control variables, we develop a constructive two-step identification argument. The first step identifies nonparametrically the conditional expectation of the disturbances given the regressors and the control variables, and the second step uses ``between-group'' variations, correcting for endogeneity, to identify the APE. We propose a natural semiparametric estimator of the APE, show its $\sqrt{n}$ asymptotic normality and compute its asymptotic variance. The estimator is computationally easy to implement, and Monte Carlo simulations show favorable finite sample properties. Control variables arise in various economic and econometric models, and we propose applications of our argument in several models. As an empirical illustration, we estimate the average elasticity of intertemporal substitution in a labor supply model with random coefficients.


title: Lower bounds for discrete negative moments of the Riemann zeta function

abstract: We prove lower bounds for the discrete negative $2k$th moment of the derivative of the Riemann zeta function for all fractional $k\geqslant 0$. The bounds are in line with a conjecture of Gonek and Hejhal. Along the way, we prove a general formula for the discrete twisted second moment of the Riemann zeta function. This agrees with a conjecture of Conrey and Snaith.


title: Linear instability of viscoelastic pipe flow

abstract: A modal stability analysis shows that pressure-driven pipe flow of an Oldroyd-B fluid is linearly unstable to axisymmetric perturbations, in stark contrast to its Newtonian counterpart which is linearly stable at all Reynolds numbers. The dimensionless groups that govern stability are the Reynolds number, the elasticity number, and the ratio of solvent to solution viscosity. The unstable mode has a phase speed close to the base-state maximum over the entire unstable region in the relevant parameter space, implying that the unstable mode belongs to a class of viscoelastic center modes. Unlike the Newtonian transition which is dominated by nonlinear processes, the linear instability discussed here could be very relevant to the onset of turbulence in viscoelastic pipe flows. The prediction of an instability is, in fact, consistent with several experimental studies on pipe flow of polymer solutions, ranging from previous reports of early turbulence to the more recent discovery of elasto-inertial turbulence. The instability identified in this study comprehensively dispels the prevailing notion of pipe flow of viscoelastic fluids being linearly stable in the Reynolds-Weissenberg plane, marking a possible paradigm shift in our understanding of transition in rectilinear viscoelastic shearing flows.