title: Feasible IV Regression without Excluded Instruments

abstract: The relevance condition of Integrated Conditional Moment (ICM) estimators is significantly weaker than the conventional IV's in at least two respects: (1) consistent estimation without excluded instruments is possible, provided endogenous covariates are non-linearly mean-dependent on exogenous covariates, and (2) endogenous covariates may be uncorrelated with but mean-dependent on instruments. These remarkable properties notwithstanding, multiplicative-kernel ICM estimators suffer diminished identification strength, large bias, and severe size distortions even for a moderately sized instrument vector. This paper proposes a computationally fast linear ICM estimator that better preserves identification strength in the presence of multiple instruments and a test of the ICM relevance condition. Monte Carlo simulations demonstrate a considerably better size control in the presence of multiple instruments and a favourably competitive performance in general. An empirical example illustrates the practical usefulness of the estimator, where estimates remain plausible when no excluded instrument is used.


title: Dark Energy and the Time Dependence of Fundamental Particle Constants

abstract: The cosmic time dependencies of $G$, $\alpha$, $\hbar$ and of Standard Model parameters like the Higgs vev and elementary particle masses are studied in the framework of a new dark energy interpretation. Due to the associated time variation of rulers, many effects turn out to be invisible. However, a rather large time dependence is claimed to arise in association with dark energy measurements, and smaller ones in connection with the Standard Model.


title: Euclid After Computer Proof-checking

abstract: Euclid pioneered the concept of a mathematical theory developed from axioms by a series of justified proof steps. From the outset there were critics and improvers. In this century the use of computers to check proofs for correctness sets a new standard of rigor. How does Euclid stand up under such an examination? And what does the exercise have to teach us about geometry, mathematical foundations, and the relation of logic to truth?


title: Expansion of bundles of light rays in the Lema\^itre -- Tolman models

abstract: The locus of $\theta \equiv {k^{\mu}};_{\mu} = 0$ for bundles of light rays emitted at noncentral points is investigated for Lema\^{\i}tre -- Tolman (L--T) models. The three loci that coincide for a central emission point: (1) maxima of $R$ along the rays, (2) $\theta = 0$, (3) $R = 2M$ are all different for a noncentral emitter. If an extremum of $R$ along a nonradial ray exists, then it must lie in the region $R > 2M$. In $2M < R \leq 3M$ it can only be a maximum; in $R > 3M$ both minima and maxima can exist. The intersection of (1) with the equatorial hypersurface (EHS) $\vartheta = \pi/2$ is numerically determined for an exemplary toy model (ETM), for two typical emitter locations. The equation of (2) is derived for a general L--T model, and its intersection with the EHS in the ETM is numerically determined for the same two emitter locations. Typically, $\theta$ has no zeros or two zeros along a ray, and becomes $+ \infty$ at the Big Crunch (BC). The only rays on which $\theta \to - \infty$ at the BC are the radial ones. Along rays on the boundaries between the no-zeros and the two-zeros regions $\theta$ has one zero, but still tends to $+ \infty$ at the BC. When the emitter is sufficiently close to the center, $\theta$ has 4 or 6 zeros along some rays (resp. 3 or 5 on the boundary rays). For noncentral emitters in a collapsing L--T model, $R = 2M$ is still the ultimate barrier behind which events become invisible from outside; loci (1) and (2) are not such barriers.


title: Asymptotic and Finite-time Cluster Synchronization of Neural Networks via Two Different Controllers

abstract: In this paper, by using pinning impulse controller and hybrid controller respectively, the research difficulties of asymptotic synchronization and finite time cluster synchronization of time-varying delayed neural networks are studied. On the ground of Lyapunov stability theorem and Lyapunov-Razumikhin method, a novel sufficient criterion on asymptotic cluster synchronization of time-varying delayed neural networks is obtained. Utilizing Finite time stability theorem and hybrid control technology, a sufficient criterion on finite-time cluster synchronization is also obtained. In order to deal with time-varying delay and save control cost, pinning pulse control is introduced to promote the realization of asymptotic cluster synchronization. Following the idea of pinning control scheme, we design a progressive hybrid control to promote the realization of finite time cluster synchronization. Finally, an example is given to illustrate the theoretical results.


title: On quantum states with a finite-dimensional approximation property

abstract: We consider a class (convex set) of quantum states containing all finite rank states and infinite rank states with the sufficient rate of decreasing of eigenvalues (in particular, all Gaussian states). Quantum states from this class are characterized by the property (called the FA-property) that allows to obtain several results concerning finite-dimensional approximation of basic entropic and information characteristics of quantum systems and channels. We obtain a simple sufficient condition of the FA-property. We show that this property implies finiteness of the von Neumann entropy, but leave unsolved the question concerning the converse implication. We obtain uniform approximation results for characteristics depending on a pair (channel, input state) and for characteristics depending on a pair (channel, input ensemble). We establish the uniform continuity of the above characteristics as functions of a channel w.r.t. the strong convergence provided that the FA-property holds either for the input state or for the average state of input ensemble.


title: Evaluation of soccer team defense based on prediction models of ball recovery and being attacked: A pilot study

abstract: With the development of measurement technology, data on the movements of actual games in various sports can be obtained and used for planning and evaluating the tactics and strategy. Defense in team sports is generally difficult to be evaluated because of the lack of statistical data. Conventional evaluation methods based on predictions of scores are considered unreliable because they predict rare events throughout the game. Besides, it is difficult to evaluate various plays leading up to a score. In this study, we propose a method to evaluate team defense from a comprehensive perspective related to team performance by predicting ball recovery and being attacked, which occur more frequently than goals, using player actions and positional data of all players and the ball. Using data from 45 soccer matches, we examined the relationship between the proposed index and team performance in actual matches and throughout a season. Results show that the proposed classifiers predicted the true events (mean F1 score $>$ 0.483) better than the existing classifiers which were based on rare events or goals (mean F1 score $<$ 0.201). Also, the proposed index had a moderate correlation with the long-term outcomes of the season ($r =$ 0.397). These results suggest that the proposed index might be a more reliable indicator rather than winning or losing with the inclusion of accidental factors.


title: Regularized Covariance Estimation for Polarization Radar Detection in Compound Gaussian Sea Clutter

abstract: This paper investigates regularized estimation of Kronecker-structured covariance matrices (CM) for polarization radar in sea clutter scenarios where the data are assumed to follow the complex, elliptically symmetric (CES) distributions with a Kronecker-structured CM. To obtain a well-conditioned estimate of the CM, we add penalty terms of Kullback-Leibler divergence to the negative log-likelihood function of the associated complex angular Gaussian (CAG) distribution. This is shown to be equivalent to regularizing Tyler's fixed-point equations by shrinkage. A sufficient condition that the solution exists is discussed. An iterative algorithm is applied to solve the resulting fixed-point iterations and its convergence is proved. In order to solve the critical problem of tuning the shrinkage factors, we then introduce two methods by exploiting oracle approximating shrinkage (OAS) and cross-validation (CV). The proposed estimator, referred to as the robust shrinkage Kronecker estimator (RSKE), is shown to achieve better performance compared with several existing methods when the training samples are limited. Simulations are conducted for validating the RSKE and demonstrating its high performance by using the IPIX 1998 real sea data.


title: Detecting Anomalous Swarming Agents with Graph Signal Processing

abstract: Collective motion among biological organisms such as insects, fish, and birds has motivated considerable interest not only in biology but also in distributed robotic systems. In a robotic or biological swarm, anomalous agents (whether malfunctioning or nefarious) behave differently than the normal agents and attempt to hide in the "chaos" of the swarm. By defining a graph structure between agents in a swarm, we can treat the agents' properties as a graph signal and use tools from the field of graph signal processing to understand local and global swarm properties. Here, we leverage this idea to show that anomalous agents can be effectively detected using their impacts on the graph Fourier structure of the swarm.


title: Optimal periodic resource allocation in reactive dynamical systems: Application to microalgal production

abstract: In this article, we focus on a periodic resource allocation problem applied to a dynamical system which comes from a biological system. More precisely, we consider a system with $N$ resources and $N$ activities, each activity use the allocated resource to evolve up to a given time $T > 0$ where a control (represented by a given permutation) will be applied on the system to reallocate the resources. The goal is to find the optimal control strategies which optimize the cost or the benefit of the system. This problem can be illustrated by an industrial biological application, namely, the optimization of a mixing strategy to enhance the growth rate in a microalgal raceway system. A mixing device, such as a paddle wheel, is considered to control the rearrangement of the depth of the algae cultures, hence the light perceived at each lap. We prove that if the dynamics of the system is periodic, then the period corresponds to one reallocation whatever the order of the involved permutation matrix is. A nonlinear optimization problem for one reallocation process is then introduced. Since $N!$ permutations need to be tested in the general case, it can be numerically solved only for a limited number of $N$. To overcome this difficulty, we introduce a second optimization problem which provides a suboptimal solution of the initial problem, but whose solution can be determined explicitly. A sufficient condition to characterize cases where the two problems have the same solution is given. Some numerical experiments are performed to assess the benefit of optimal strategies in various settings.


title: The rational Sklyanin algebra and the Wilson and para-Racah polynomials

abstract: The relation between Wilson and para-Racah polynomials and representations of the degenerate rational Sklyanin algebra is established. Second order Heun operators on quadratic grids with no diagonal terms are determined. These special or S-Heun operators lead to the rational degeneration of the Sklyanin algebra; they also entail the contiguity and structure operators of the Wilson polynomials. The finite-dimensional restriction yields a representation that acts on the para-Racah polynomials.


title: Inferred vs traditional personality assessment: are we predicting the same thing?

abstract: Machine learning methods are widely used by researchers to predict psychological characteristics from digital records. To find out whether automatic personality estimates retain the properties of the original traits, we reviewed 220 recent articles. First, we put together the predictive quality estimates from a subset of the studies which declare separation of training, validation, and testing phases, which is critical for ensuring the correctness of quality estimates in machine learning. Only 20% of the reviewed papers met this criterion. To compare the reported quality estimates, we converted them to approximate Pearson correlations. The credible upper limits for correlations between predicted and self-reported personality traits vary in a range between 0.42 and 0.48, depending on the specific trait. The achieved values are substantially below the correlations between traits measured with distinct self-report questionnaires. This suggests that we cannot readily interpret personality predictions as estimates of the original traits or expect predicted personality traits to reproduce known relationships with life outcomes regularly. Next, we complement quality estimates evaluation with evidence on psychometric properties of predicted traits. The few existing results suggest that predicted traits are less stable with time and have lower effective dimensionality than self-reported personality. The predictive text-based models perform substantially worse outside their training domains but stay above a random baseline. The evidence on the relationships between predicted traits and external variables is mixed. Predictive features are difficult to use for validation, due to the lack of prior hypotheses. Thus, predicted personality traits fail to retain important properties of the original characteristics. This calls for the cautious use and targeted validation of the predictive models.


title: Physics potential for the H$\rightarrow$ZZ$^{*}$ decay at the CEPC

abstract: The precision of the yield measurement of the Higgs boson decaying into a pair of $Z$ bosons process at the Circular Electron Positron Collider (CEPC) is evaluated. Including the recoil $Z$ boson associated with the Higgs production (Higgsstrahlung) total three $Z$ bosons are involved for this channel, from which final states characterized by the presence of a pair of leptons, quarks, and neutrinos are chosen for the signal. Two analysis approaches are compared and the final precision of ${\sigma}_{\mathrm{ZH}}{\cdot}$BR($H \rightarrow ZZ^{*}$) is estimated to be 7.9% using a multivariate analysis technique, based on boosted decision trees. The relative precision of the Higgs boson width, using this $H \rightarrow ZZ^{*}$ decay topology, is estimated by combining the obtained result with the precision of the inclusive \textit{ZH} cross section measurement.


title: Adaptation to a heterogeneous patchy environment with nonlocal dispersion

abstract: In this work, we provide an asymptotic analysis of the solutions to an elliptic integro-differential equation. This equation describes the evolutionary equilibria of a phenotypically structured population, subject to selection, mutation, and both local and non-local dispersion in a spatially heterogeneous, possibly patchy, environment. Considering small effects of mutations, we provide an asymptotic description of the equilibria of the phenotypic density. This asymptotic description involves a Hamilton-Jacobi equation with constraint coupled with an eigenvalue problem. Based on such analysis, we characterize some qualitative properties of the phenotypic density at equilibrium depending on the heterogeneity of the environment. In particular, we show that when the heterogeneity of the environment is low, the population concentrates around a single phenotypic trait leading to a unimodal phenotypic distribution. On the contrary, a strong fragmentation of the environment leads to multi-modal phenotypic distributions.


title: SILT: Efficient transformer training for inter-lingual inference

abstract: The ability of transformers to perform precision tasks such as question answering, Natural Language Inference (NLI) or summarising, have enabled them to be ranked as one of the best paradigm to address Natural Language Processing (NLP) tasks. NLI is one of the best scenarios to test these architectures, due to the knowledge required to understand complex sentences and established relationships between a hypothesis and a premise. Nevertheless, these models suffer from incapacity to generalise to other domains or difficulties to face multilingual and interlingual scenarios. The leading pathway in the literature to address these issues involve designing and training extremely large architectures, which leads to unpredictable behaviours and to establish barriers which impede broad access and fine tuning. In this paper, we propose a new architecture called Siamese Inter-Lingual Transformer (SILT), to efficiently align multilingual embeddings for Natural Language Inference, allowing for unmatched language pairs to be processed. SILT leverages siamese pre-trained multi-lingual transformers with frozen weights where the two input sentences attend each other to later be combined through a matrix alignment method. The experimental results carried out in this paper evidence that SILT allows to reduce drastically the number of trainable parameters while allowing for inter-lingual NLI and achieving state-of-the-art performance on common benchmarks. We make our code and dataset available at https://github.com/jahuerta92/siamese-inter-lingual-transformer.


title: Accretive Computation of Global Transformations

abstract: Global transformations form a categorical framework adapting graph transformations to describe fully synchronous rule systems on a given data structure.In this work we focus on data structures that can be captured as presheaves and study the computational aspects of such synchronous rule systems.To obtain an online algorithm, a complete study of the sub-steps within each synchronous step is done at the semantic level.This leads to the definition of accretive rule systems and a local criterion to characterize these systems.Finally an online computation algorithm for theses systems is given.


title: X-ray emitting structures in the Vela SNR: ejecta anisotropies and progenitor stellar wind residuals

abstract: The Vela supernova remnant (SNR) shows several ejecta fragments protruding beyond the forward shock (shrapnel). Recent studies have revealed high Si abundance in two shrapnel (A and G), located in opposite directions with respect to the SNR center. This suggests the possible existence of a Si-rich jet-counterjet structure. We analyzed an XMM-Newton observation of a bright clump, behind shrapnel G, which lies along the direction connecting A and G. The aim is to study the physical and chemical properties of this clump to ascertain whether it is part of this putative jet-like structure. We produced background-corrected and adaptively-smoothed count-rate images and median photon energy maps, and performed a spatially resolved spectral analysis. We identified two structures with different physical properties. The first one is remarkably elongated along the direction connecting A and G. Its X-ray spectrum is much softer than that of the other two shrapnel, to the point of hindering the determination of the Si abundance, however its physical and chemical properties are consistent with those of shrapnel A and G. The second structure, running along the southeast-northwest direction, has a higher temperature and appears like a thin filament. By analyzing the ROSAT data, we have found that this filament is part of a very large and coherent structure that we identified in the western rim of the shell. We obtained a thorough description of the tail of Shrapnel G. In addition we discovered a coherent and very extended feature that we interpret as a signature of an earlier interaction of the remnant with the stellar wind of its progenitor star. The peculiar Ne/O ratio we found in the wind residual may be suggestive of a Wolf-Rayet progenitor for Vela SNR, though further analysis is required to address this point.


title: On $L^{2}$-harmonic forms of complete almost K\"{a}hler manifold

abstract: In this article, we study the $L^{2}$-harmonic forms on the complete $2n$-dimensional almost K\"{a}her manifold $X$. We observe that the $L^{2}$-harmonic forms can decomposition into Lefschetz powers of primitive forms. Therefore we can extend vanishing theorems of $d$(bounded) (resp. $d$(sublinear)) K\"{a}hler manifold proved by Gromov (resp. Cao-Xavier, Jost-Zuo) to almost K\"{a}hlerian case, that is, the spaces of all harmonic $(p,q)$-forms on $X$ vanishing unless $p+q=n$. We also give a lower bound on the spectra of the Laplace operator to sharpen the Lefschetz vanishing theorem on $d$(bounded) case.


title: Congestion control in high-speed networks using the probabilistic estimation approach

abstract: Nowadays, the bulk of Internet traffic uses TCP protocol for reliable transmission. But the standard TCP's performance is very poor in High Speed Networks (HSN) and hence the core gigabytes links are usually underutilization. This problem has roots in conservative nature of TCP, especially in its Additive Increase Multiplicative Decrease (AIMD) phase. In other words, since TCP can't figure out precisely the congestion status of the network, it follows a conservative strategy to keep the network from overwhelming. We believe that precisely congestion estimation in the network can solve this problem by avoiding unnecessary conservation. To this end, this paper proposes an algorithm which considers packet loss and delay information jointly and employs a probabilistic approach to accurately estimation of congestion status in the network. To examine the proposed scheme performance, extensive simulations have been performed in the NS-2 environment. Simulation results reveal that the proposed algorithm has better performance than existing algorithms in terms of bottleneck utilization, stability, throughput and fairness.


title: Constructive exact control of semilinear 1D heat equations

abstract: The exact distributed controllability of the semilinear heat equation $\partial_{t}y-\Delta y + g(y)=f \,1_{\omega}$ posed over multi-dimensional and bounded domains, assuming that $g\in C^1(\mathbb{R})$ satisfies the growth condition $\limsup_{r\to \infty} g(r)/(\vert r\vert \ln^{3/2}\vert r\vert)=0$ has been obtained by Fern\'andez-Cara and Zuazua in 2000. The proof based on a non constructive fixed point arguments makes use of precise estimates of the observability constant for a linearized heat equation. In the one dimensional setting, assuming that $g^\prime$ does not grow faster than $\beta \ln^{3/2}\vert r\vert$ at infinity for $\beta>0$ small enough and that $g^\prime$ is uniformly H\"older continuous on $\mathbb{R}$ with exponent $p\in [0,1]$, we design a constructive proof yielding an explicit sequence converging to a controlled solution for the semilinear equation, at least with order $1+p$ after a finite number of iterations.


title: A Robust Extrinsic Calibration Framework for Vehicles with Unscaled Sensors

abstract: Accurate extrinsic sensor calibration is essential for both autonomous vehicles and robots. Traditionally this is an involved process requiring calibration targets, known fiducial markers and is generally performed in a lab. Moreover, even a small change in the sensor layout requires recalibration. With the anticipated arrival of consumer autonomous vehicles, there is demand for a system which can do this automatically, after deployment and without specialist human expertise. To solve these limitations, we propose a flexible framework which can estimate extrinsic parameters without an explicit calibration stage, even for sensors with unknown scale. Our first contribution builds upon standard hand-eye calibration by jointly recovering scale. Our second contribution is that our system is made robust to imperfect and degenerate sensor data, by collecting independent sets of poses and automatically selecting those which are most ideal. We show that our approach's robustness is essential for the target scenario. Unlike previous approaches, ours runs in real time and constantly estimates the extrinsic transform. For both an ideal experimental setup and a real use case, comparison against these approaches shows that we outperform the state-of-the-art. Furthermore, we demonstrate that the recovered scale may be applied to the full trajectory, circumventing the need for scale estimation via sensor fusion.


title: eduPIC: an introductory particle based code for radio-frequency plasma simulation

abstract: For the self-consistent description of various plasma sources operated in the low-pressure (nonlocal, kinetic) regime, the Particle-In-Cell simulation approach, combined with the Monte Carlo treatment of collision processes (PIC/MCC), has become an important tool during the past decades. PIC/MCC simulation codes have been developed and maintained by many research groups, some of these codes are available to the community as freeware resources. While this computational approach has already been present for a number of decades, the rapid evolution of the computing infrastructure makes it increasingly more popular and accessible, as simulations of simple systems can be executed now on personal computers or laptops. During the past few years we have experienced an increasing interest in lectures and courses dealing with the basics of particle simulations, including the PIC/MCC technique. In a response to this, this paper (i) provides a tutorial on the physical basis and the algorithms of the PIC/MCC technique and (ii) presents a basic (spatially one-dimensional) electrostatic PIC/MCC simulation code for Capacitively Coupled Plasmas, whose source is made freely available in various programming languages. We share the code in C/C++ versions, as well as in a version written in Rust, which is a rapidly emerging computational language. Our code intends to be a "starting tool" for those who are interested in learning the details of the PIC/MCC technique and would like to develop the "skeleton" code further, for their research purposes.


title: HADES RV Programme with HARPS-N at TNG XIII. A sub-Neptune around the M dwarf GJ 720 A

abstract: Context. The high number of super-Earth and Earth-like planets in the habitable zone (HZ) detected around M-dwarf stars in the last years has revealed these stellar objects to be the key for planetary radial velocity (RV) searches. Aims. Using the HARPS-N spectrograph within The HArps-n red Dwarf Exoplanet Survey (HADES) we reach the precision needed to detect small planets with a few Earth masses using the RV technique. Methods. We obtained 138 HARPS-N RV measurements between 2013 May and 2020 September of GJ 720 A, classified as an M0.5V star located at a distance of 15.56 pc. To characterize the stellar variability and to discern the periodic variation due to the Keplerian signals from those related to stellar activity, the HARPS-N spectroscopic activity indicators and the simultaneous photometric observations were analyzed. The combined analysis of HARPS-N RVs and activity indicators let us to address the nature of the periodic signals. The final model and the orbital planetary parameters were obtained by fitting simultaneously the stellar variability and the Keplerian signal using a Gaussian process regression and following a Bayesian criterion. Results. The HARPS-N RV periodic signals around 40 d and 100 d have counterparts at the same frequencies in HARPS-N activity indicators and photometric light curves. Then we attribute these periodicities to stellar activity the former period being likely associated with the stellar rotation. GJ 720 A shows the most significant signal at 19.466$\pm$0.005 d with no counterparts in any stellar activity indices. We hence ascribe this RV signal, having a semiamplitude of 4.72$\pm$0.27 m/s , to the presence of a sub-Neptune mass planet. The planet GJ 720 Ab has a minimum mass of 13.64$\pm$0.79 M$_{\oplus}$, it is in circular orbit at 0.119$\pm$0.002 AU from its parent star, and lies inside the inner boundary of the HZ around its parent star.


title: Extending Representation Formulae for Boundary Voltage Perturbations of Low Volume Fraction to Very Contrasted Conductivity Inhomogeneities

abstract: Imposing either Dirichlet or Neumann boundary conditions on the boundary of a smooth bounded domain $\Omega$, we study the perturbation incurred by the voltage potential when the conductivity is modified in a set of small measure. We consider $\left(\gamma_{n}\right)_{n\in\mathbb{N}}$, a sequence of perturbed conductivity matrices differing from a smooth $\gamma_{0}$ background conductivity matrix on a measurable set well within the domain, and we assume $\left(\gamma_{n}-\gamma_{0}\right)\gamma_{n}^{-1}\left(\gamma_{n}-\gamma_{0}\right)\to0$ in $L^{1}(\Omega)$. Adapting the limit measure, we show that the general representation formula introduced for bounded contrasts in \citep{capdeboscq-vogelius-03a} can be extended to unbounded sequencesof matrix valued conductivities.


title: UniParma at SemEval-2021 Task 5: Toxic Spans Detection Using CharacterBERT and Bag-of-Words Model

abstract: With the ever-increasing availability of digital information, toxic content is also on the rise. Therefore, the detection of this type of language is of paramount importance. We tackle this problem utilizing a combination of a state-of-the-art pre-trained language model (CharacterBERT) and a traditional bag-of-words technique. Since the content is full of toxic words that have not been written according to their dictionary spelling, attendance to individual characters is crucial. Therefore, we use CharacterBERT to extract features based on the word characters. It consists of a CharacterCNN module that learns character embeddings from the context. These are, then, fed into the well-known BERT architecture. The bag-of-words method, on the other hand, further improves upon that by making sure that some frequently used toxic words get labeled accordingly. With a 4 percent difference from the first team, our system ranked 36th in the competition. The code is available for further re-search and reproduction of the results.


title: Quantitative De Giorgi methods in kinetic theory

abstract: We consider hypoelliptic equations of kinetic Fokker-Planck type, also known as Kolmogorov or ultraparabolic equations, with rough coefficients in the drift-diffusion operator. We give novel short quantitative proofs of the De Giorgi intermediate-value Lemma as well as weak Harnack and Harnack inequalities. This implies H{\"o}lder continuity with quantitative estimates. The paper is self-contained.


title: Best Practices for Collecting Gender and Sex Data

abstract: The measurement and analysis of human sex and gender is a nuanced problem with many overlapping considerations including statistical bias, data privacy, and the ethical treatment of study subjects. Traditionally, human gender and sex have been categorized and measured with respect to an artificial binary system. The continuation of this tradition persists mainly because it is easy to replication and not, as we argue, because it produces the most valuable scientific information. Sex and gender identity data is crucial for many applications of statistical analysis and many modern scientists acknowledge the limitations of the current system. However, discrimination against sex and gender minorities poses very real privacy concerns when collecting and distributing gender and sex data. As such, extra thoughtfulness and care is essential to design safe and informative scientific studies. In this paper, we present statistically informed recommendations for the data collection and analysis of human subjects that not only respect each individual's identity and protect their privacy, but also establish standards for collecting higher quality data.


title: Hellinger and total variation distance in approximating L{\'e}vy driven SDEs

abstract: In this paper, we get some convergence rates in total variation distance in approximating discretized paths of L{\'e}vy driven stochastic differential equations, assuming that the driving process is locally stable. The particular case of the Euler approximation is studied. Our results are based on sharp local estimates in Hellinger distance obtained using Malliavin calculus for jump processes.


title: Textile D forms and $D_{4d}$

abstract: D-forms have in the past been created from inflexible materials, or considered as abstract mathematical objects. This paper describes a number of realisations of D-forms, and the related pita-forms, in textiles. Examples are given in which the created 3D space is purposefully employed and others in which ornamentation of the constituent surfaces is the highlighted feature. In particular, a set of biscornu has been fashioned to provide a 3D sampler of the axial point group $D_{4d}$ and its subgroups, using hitomezashi.


title: Numerical Simulations on Nonlinear Quantum Graphs with the GraFiDi Library

abstract: Nonlinear quantum graphs are metric graphs equipped with a nonlinear Schr{\"o}dinger equation. Whereas in the last ten years they have known considerable developments on the theoretical side, their study from the numerical point of view remains in its early stages. The goal of this paper is to present the Grafidi library, a Python library which has been developed with the numerical simulation of nonlinear Schr{\"o}dinger equations on graphs in mind. We will show how, with the help of the Grafidi library, one can implement the popular normalized gradient flow and nonlinear conjugate gradient flow methods to compute ground states of a nonlinear quantum graph. We will also simulate the dynamics of the nonlinear Schr{\"o}dinger equation with a Crank-Nicolson relaxation scheme and a Strang splitting scheme. Finally, in a series of numerical experiments on various types of graphs, we will compare the outcome of our numerical calculations for ground states with the existing theoretical results, thereby illustrating the versatility and efficiency of our implementations in the framework of the Grafidi library.


title: Machine learning initialization to accelerate Stokes profile inversions

abstract: In this work, we discuss the application of convolutional neural networks (CNNs) as a tool to advantageously initialize Stokes profile inversions. To demonstrate the usefulness of CNNs, we concentrate in this paper on the inversion of LTE Stokes profiles. We use observations taken with the spectropolarimeter onboard the Hinode spacecraft as a test benchmark. First, we carefully analyze the data with the SIR inversion code using a given initial atmospheric model. The code provides a set of atmospheric models that reproduce the observations. These models are then used to train a CNN. Afterwards, the same data are again inverted with SIR but using the trained CNN to provide the initial guess atmospheric models for SIR. The CNNs allow us to significantly reduce the number of inversion cycles when used to compute initial guess model atmospheres, decreasing the computational time for LTE inversions by a factor of two to four. CNN's alone are much faster than assisted inversions, but the latter are more robust and accurate. The advantages and limitations of machine learning techniques for estimating optimum initial atmospheric models for spectral line inversions are discussed. Finally, we describe a python wrapper for the SIR and DeSIRe codes that allows for the easy setup of parallel inversions. The assisted inversions can speed up the inversion process, but the efficiency and accuracy of the inversion results depend strongly on the solar scene and the data used for the CNN training. This method (assisted inversions) will not obviate the need for analyzing individual events with the utmost care but will provide solar scientists with a much better opportunity to sample large amounts of inverted data, which will undoubtedly broaden the physical discovery space.


title: Identification of materials with strong magneto-structural coupling using computational high-throughput screening

abstract: Important phenomena such as magnetostriction, magnetocaloric, and magnetoelectric effects arise from, or could be enhanced by, the coupling of magnetic and structural degrees of freedom. The coupling of spin and lattice also influence transport and structural properties in magnetic materials in particular around phase transitions. In this paper we propose a method for screening materials for a strong magneto-structural coupling by assessing the effect of the local magnetic configuration on the atomic forces using density functional theory (DFT). We have employed the disordered local moment approach in a supercell formulation to probe different magnetic local configurations and their forces and performed a high-throughput search on binary and ternary compounds available in the Crystallographic Open Database. We identify a list of materials with a strong spin-lattice coupling out of which several are already known to display magneto-lattice coupling-phenomena like Fe3O4 and CrN. Others, such as Mn2CrO4 and CaFe7O11 have been less studied and are yet to reveal their potentials in experiments and applications.


title: Representations of integers as sums of four polygonal numbers and partial theta functions

abstract: In this paper, we consider representations of integers as sums of at most four distinct $m$-gonal numbers (allowing a fixed number of repeats of each polygonal number occurring in the sum). We show that the number of such representations with non-negative parameters (hence counting the number of points in a regular $m$-gon) is asymptotically the same as $\frac{1}{16}$ times the number of such representations with arbitrary integer parameters (often called generalized polygonal numbers).


title: Ramanujan in Computing Technology

abstract: This paper is a tribute to the genius of the legendary Indian mathematician Srinivasa Ramanujan (22 December 1887 - 26 April 1920) in the centenary year of his death. The life story of Ramanujan is so well known that it needs no elaboration not even a summarisation. In his short life period he made substantial contributions to mathematical analysis, number theory, infinite series, and continued fractions, including solutions to mathematical problems then considered unsolvable. Ramanujan independently compiled nearly 3,900 results in the form of identities and equations. Many were completely novel; his original and highly unconventional results, such as the Ramanujan prime, the Ramanujan theta function, partition formulae and mock theta functions, have opened entire new areas of work and inspired a vast amount of further research. Nearly all his claims have now been proven correct. The focus of the paper is the increasing influence of the ideas propounded by Ramanujan in the development of computing technology. We shall discuss the application of certain infinite series discovered by Ramanujan in computing the value of the mathematical constant $\pi$. We shall also consider certain special graphs known as Ramanujan graphs and the reason for designating them as such. We shall examine how certain researchers are attempting to create an abstract machine which they call Ramanujan machine which is thought of as simulating the hypothesised thought process of Ramnujan. We shall also have a brief look at the applications of Ramanujan's discoveries in signal processing.


title: The Old and the New: Can Physics-Informed Deep-Learning Replace Traditional Linear Solvers?

abstract: Physics-Informed Neural Networks (PINN) are neural networks encoding the problem governing equations, such as Partial Differential Equations (PDE), as a part of the neural network. PINNs have emerged as a new essential tool to solve various challenging problems, including computing linear systems arising from PDEs, a task for which several traditional methods exist. In this work, we focus first on evaluating the potential of PINNs as linear solvers in the case of the Poisson equation, an omnipresent equation in scientific computing. We characterize PINN linear solvers in terms of accuracy and performance under different network configurations (depth, activation functions, input data set distribution). We highlight the critical role of transfer learning. Our results show that low-frequency components of the solution converge quickly as an effect of the F-principle. In contrast, an accurate solution of the high frequencies requires an exceedingly long time. To address this limitation, we propose integrating PINNs into traditional linear solvers. We show that this integration leads to the development of new solvers whose performance is on par with other high-performance solvers, such as PETSc conjugate gradient linear solvers, in terms of performance and accuracy. Overall, while the accuracy and computational performance are still a limiting factor for the direct use of PINN linear solvers, hybrid strategies combining old traditional linear solver approaches with new emerging deep-learning techniques are among the most promising methods for developing a new class of linear solvers.


title: Set-to-Sequence Methods in Machine Learning: a Review

abstract: Machine learning on sets towards sequential output is an important and ubiquitous task, with applications ranging from language modeling and meta-learning to multi-agent strategy games and power grid optimization. Combining elements of representation learning and structured prediction, its two primary challenges include obtaining a meaningful, permutation invariant set representation and subsequently utilizing this representation to output a complex target permutation. This paper provides a comprehensive introduction to the field as well as an overview of important machine learning methods tackling both of these key challenges, with a detailed qualitative comparison of selected model architectures.


title: Entanglement dualities in supersymmetry

abstract: We derive a general relation between the bosonic and fermionic entanglement in the ground states of supersymmetric quadratic Hamiltonians. For this, we construct canonical identifications between bosonic and fermionic subsystems. Our derivation relies on a unified framework to describe both, bosonic and fermionic Gaussian states in terms of so-called linear complex structures $J$. The resulting dualities apply to the full entanglement spectrum between the bosonic and the fermionic systems, such that the von Neumann entropy and arbitrary Renyi entropies can be related. We illustrate our findings in one and two-dimensional systems, including the paradigmatic Kitaev honeycomb model. While typically SUSY preserves features like area law scaling of the entanglement entropies on either side, we find a peculiar phenomenon, namely, an amplified scaling of the entanglement entropy ("super area law") in bosonic subsystems when the dual fermionic subsystems develop almost maximally entangled modes.


title: Modeling and simulation of nuclear architecture reorganization process using a phase field approach

abstract: We develop a special phase field/diffusive interface method to model the nuclear architecture reorganization process. In particular, we use a Lagrange multiplier approach in the phase field model to preserve the specific physical and geometrical constraints for the biological events. We develop several efficient and robust linear and weakly nonlinear schemes for this new model. To validate the model and numerical methods, we present ample numerical simulations which in particular reproduce several processes of nuclear architecture reorganization from the experiment literature.


title: Pseudo Supervised Solar Panel Mapping based on Deep Convolutional Networks with Label Correction Strategy in Aerial Images

abstract: Solar panel mapping has gained a rising interest in renewable energy field with the aid of remote sensing imagery. Significant previous work is based on fully supervised learning with classical classifiers or convolutional neural networks (CNNs), which often require manual annotations of pixel-wise ground-truth to provide accurate supervision. Weakly supervised methods can accept image-wise annotations which can help reduce the cost for pixel-level labelling. Inevitable performance gap, however, exists between weakly and fully supervised methods in mapping accuracy. To address this problem, we propose a pseudo supervised deep convolutional network with label correction strategy (PS-CNNLC) for solar panels mapping. It combines the benefits of both weak and strong supervision to provide accurate solar panel extraction. First, a convolutional neural network is trained with positive and negative samples with image-level labels. It is then used to automatically identify more positive samples from randomly selected unlabeled images. The feature maps of the positive samples are further processed by gradient-weighted class activation mapping to generate initial mapping results, which are taken as initial pseudo labels as they are generally coarse and incomplete. A progressive label correction strategy is designed to refine the initial pseudo labels and train an end-to-end target mapping network iteratively, thereby improving the model reliability. Comprehensive evaluations and ablation study conducted validate the superiority of the proposed PS-CNNLC.


title: Scalable Hypergraph Embedding System

abstract: Many problems such as node classification and link prediction in network data can be solved using graph embeddings. However, it is difficult to use graphs to capture non-binary relations such as communities of nodes. These kinds of complex relations are expressed more naturally as hypergraphs. While hypergraphs are a generalization of graphs, state-of-the-art graph embedding techniques are not adequate for solving prediction and classification tasks on large hypergraphs accurately in reasonable time. In this paper, we introduce HyperNetVec, a novel hierarchical framework for scalable unsupervised hypergraph embedding. HyperNetVec exploits shared-memory parallelism and is capable of generating high quality embeddings for real-world hypergraphs with millions of nodes and hyperedges in only a couple of minutes while existing hypergraph systems either fail for such large hypergraphs or may take days to produce the embeddings.


title: Stable Sheaves on K3 Surfaces via Wall-Crossing

abstract: We give a new proof of the following theorem: moduli spaces of stable complexes on a complex projective K3 surface, with primitive Mukai vector and with respect to a generic Bridgeland stability condition, are hyperk\"{a}hler varieties of $\mathrm{K3}^{[n]}$-type of expected dimension. We use derived equivalences, deformations and wall-crossing for Bridgeland stability to reduce to the case of the Hilbert scheme of points.


title: Physics-Informed Neural Network Method for Solving One-Dimensional Advection Equation Using PyTorch

abstract: Numerical solutions to the equation for advection are determined using different finite-difference approximations and physics-informed neural networks (PINNs) under conditions that allow an analytical solution. Their accuracy is examined by comparing them to the analytical solution. We used a machine learning framework like PyTorch to implement PINNs. PINNs approach allows training neural networks while respecting the PDEs as a strong constraint in the optimization as apposed to making them part of the loss function. In standard small-scale circulation simulations, it is shown that the conventional approach incorporates a pseudo diffusive effect that is almost as large as the effect of the turbulent diffusion model; hence the numerical solution is rendered inconsistent with the PDEs. This oscillation causes inaccuracy and computational uncertainty. Of all the schemes tested, only the PINNs approximation accurately predicted the outcome. We assume that the PINNs approach can transform the physics simulation area by allowing real-time physics simulation and geometry optimization without costly and time-consuming simulations on large supercomputers.


title: An open-source ABAQUS implementation of the scaled boundary finite element method to study interfacial problems using polyhedral meshes

abstract: The scaled boundary finite element method (SBFEM) is capable of generating polyhedral elements with an arbitrary number of surfaces. This salient feature significantly alleviates the meshing burden being a bottleneck in the analysis pipeline in the standard finite element method (FEM). In this paper, we implement polyhedral elements based on the SBFEM into the commercial finite element software ABAQUS. To this end, user elements are provided through the user subroutine UEL. Detailed explanations regarding the data structures and implementational aspects of the procedures are given. The focus of the current implementation is on interfacial problems and therefore, element-based surfaces are created on polyhedral user elements to establish interactions. This is achieved by an overlay of standard finite elements with negligible stiffness, provided in the ABAQUS element library, with polyhedral user elements. By means of several numerical examples, the advantages of polyhedral elements regarding the treatment of non-matching interfaces and automatic mesh generation are clearly demonstrated. Thus, the performance of ABAQUS for problems involving interfaces is augmented based on the availability of polyhedral meshes. Due to the implementation of polyhedral user elements, ABAQUS can directly handle complex geometries given in the form of digital images or stereolithography (STL) files. In order to facilitate the use of the proposed approach, the code of the UEL is published open-source and can be downloaded from https://github.com/ShukaiYa/SBFEM-UEL.


title: Numerical simulation and entropy dissipative cure of the carbuncle instability for the shallow water circular hydraulic jump

abstract: We investigate the numerical artifact known as a carbuncle, in the solution of the shallow water equations. We propose a new Riemann solver that is based on a local measure of the entropy residual and aims to avoid carbuncles while maintaining high accuracy. We propose a new challenging test problem for shallow water codes, consisting of a steady circular hydraulic jump that can be physically unstable. We show that numerical methods are prone to either suppress the instability completely or form carbuncles. We test existing cures for the carbuncle. In our experiments, only the proposed method is able to avoid unphysical carbuncles without suppressing the physical instability.

categories: physics.flu-dyn cs.NA math.NA nlin.CD physics.comp-ph


title: Hadronic Rescattering in pA and AA Collisions

abstract: In a recent article we presented a model for hadronic rescattering, and some results were shown for pp collisions at LHC energies. In order to extend the studies to pA and AA collisions, the Angantyr model for heavy-ion collisions is taken as the starting point. Both these models are implemented within the general-purpose Monte Carlo event generator Pythia, which makes the matching reasonably straightforward, and allows for detailed studies of the full space--time evolution. The rescattering rate is significantly higher than in pp, especially for central AA collisions, where the typical primary hadron rescatters several times. We study the impact of rescattering on a number of distributions, such as pT and eta spectra, and the space--time evolution of the whole collision process. Notably rescattering is shown to give a significant contribution to elliptic flow in XeXe and PbPb, and to give a nontrivial impact on charm production.


title: Multimodal End-to-End Sparse Model for Emotion Recognition

abstract: Existing works on multimodal affective computing tasks, such as emotion recognition, generally adopt a two-phase pipeline, first extracting feature representations for each single modality with hand-crafted algorithms and then performing end-to-end learning with the extracted features. However, the extracted features are fixed and cannot be further fine-tuned on different target tasks, and manually finding feature extraction algorithms does not generalize or scale well to different tasks, which can lead to sub-optimal performance. In this paper, we develop a fully end-to-end model that connects the two phases and optimizes them jointly. In addition, we restructure the current datasets to enable the fully end-to-end training. Furthermore, to reduce the computational overhead brought by the end-to-end model, we introduce a sparse cross-modal attention mechanism for the feature extraction. Experimental results show that our fully end-to-end model significantly surpasses the current state-of-the-art models based on the two-phase pipeline. Moreover, by adding the sparse cross-modal attention, our model can maintain performance with around half the computation in the feature extraction part.


title: Stickelberger series and Main Conjecture for function fields

abstract: Let F be a global function field of characteristic p with ring of integers A and let \Phi be a Hayes module on the Hilbert class field H(A) of F. We prove an Iwasawa Main Conjecture for the Z_p^\infty-extension F/F generated by the \mathfrak{p}-power torsion of \Phi (\mathfrak{p} a prime of A). The main tool is a Stickelberger series whose specialization provides a generator for the Fitting ideal of the class group of F. Moreover we prove that the same series, evaluated at complex or p-adic characters, interpolates the Goss Zeta-function or some p-adic L-function, thus providing the link between the algebraic structure (class groups) and the analytic functions, which is the crucial part of Iwasawa Main Conjecture.


title: Secure Hypersphere Range Query on Encrypted Data

abstract: Spatial queries like range queries, nearest neighbor, circular range queries etc. are the most widely used queries in the location-based applications. Building secure and efficient solutions for these queries in the cloud computing framework is critical and has been an area of active research. This paper focuses on the problem of Secure Circular Range Queries (SCRQ), where client submits an encrypted query (consisting of a center point and radius of the circle) and the cloud (storing encrypted data points) has to return the points lying inside the circle. The existing solutions for this problem suffer from various disadvantages such as high processing time which is proportional to square of the query radius, query generation phase which is directly proportional to the number of points covered by the query etc. This paper presents solution for the above problem which is much more efficient than the existing solutions. Three protocols are proposed with varying characteristics. It is shown that all the three protocols are secure. The proposed protocols can be extended to multiple dimensions and thus are able to handle Secure Hypersphere Range Queries (SHRQ) as well. Internally the proposed protocols use pairing-based cryptography and a concept of lookup table. To enable the efficient use of limited size lookup table, a new storage scheme is presented. The proposed storage scheme enables the protocols to handle query with much larger radius values. Using the SHRQ protocols, we also propose a mechanism to answer the Secure range Queries. Extensive performance evaluation has been done to evaluate the efficiency of the proposed protocols


title: Large-Scale Zero-Shot Image Classification from Rich and Diverse Textual Descriptions

abstract: We study the impact of using rich and diverse textual descriptions of classes for zero-shot learning (ZSL) on ImageNet. We create a new dataset ImageNet-Wiki that matches each ImageNet class to its corresponding Wikipedia article. We show that merely employing these Wikipedia articles as class descriptions yields much higher ZSL performance than prior works. Even a simple model using this type of auxiliary data outperforms state-of-the-art models that rely on standard features of word embedding encodings of class names. These results highlight the usefulness and importance of textual descriptions for ZSL, as well as the relative importance of auxiliary data type compared to algorithmic progress. Our experimental results also show that standard zero-shot learning approaches generalize poorly across categories of classes.


title: NGC 6611 601: A hot pre-main sequence spectroscopic binary containing a centrifugal magnetosphere host star

abstract: W 601 (NGC 6611 601) is one of the handful of known magnetic Herbig Ae/Be stars. We report the analysis of a large dataset of high-resolution spectropolarimetry. The star is a previously unreported spectroscopic binary, consisting of 2 B2 stars with a mass ratio of 1.8, masses of 12 M$_\odot$ and 6.2 $M_\odot$, in an eccentric 110-day orbit. The magnetic field belongs to the secondary, W 601 B. The H$\alpha$ emission is consistent with an origin in W 601 B's centrifugal magnetosphere; the star is therefore not a classical Herbig Be star in the sense that its emission is not formed in an accretion disk. However, the low value of $\log{g} = 3.8$ determined via spectroscopic analysis, and the star's membership in the young NGC 6611 cluster, are most consistent with it being on the pre-main sequence. The rotational period inferred from the variability of the H$\alpha$ line and the longitudinal magnetic field $\langle B_z \rangle$ is 1.13 d. Modelling of Stokes $V$ and $\langle B_z \rangle$ indicates a surface dipolar magnetic field $B_{\rm d}$ between 6 and $11$ kG. With its strong emission, rapid rotation, and strong surface magnetic field, W 601 B is likely a precursor to H$\alpha$-bright magnetic B-type stars such as $\sigma$ Ori E. By contrast, the primary is an apparently non-magnetic ($B_{\rm d} < 300$ G) pre-main sequence early B-type star. In accordance with expectations from magnetic braking, the non-magnetic primary is apparently more rapidly rotating than the magnetic star.


title: Fourier Transform of Percoll Gradients Boosts CNN Classification of Hereditary Hemolytic Anemias

abstract: Hereditary hemolytic anemias are genetic disorders that affect the shape and density of red blood cells. Genetic tests currently used to diagnose such anemias are expensive and unavailable in the majority of clinical labs. Here, we propose a method for identifying hereditary hemolytic anemias based on a standard biochemistry method, called Percoll gradient, obtained by centrifuging a patient's blood. Our hybrid approach consists on using spatial data-driven features, extracted with a convolutional neural network and spectral handcrafted features obtained from fast Fourier transform. We compare late and early feature fusion with AlexNet and VGG16 architectures. AlexNet with late fusion of spectral features performs better compared to other approaches. We achieved an average F1-score of 88% on different classes suggesting the possibility of diagnosing of hereditary hemolytic anemias from Percoll gradients. Finally, we utilize Grad-CAM to explore the spatial features used for classification.


title: DUETS: A Dataset of Reproducible Pairs ofJava Library-Clients

abstract: Software engineering researchers look for software artifacts to study their characteristics or to evaluate new techniques. In this paper, we introduce DUETS, a new dataset of software libraries and their clients. This dataset can be exploited to gain many different insights, such as API usage, usage inputs, or novel observations about the test suites of clients and libraries. DUETS is meant to support both static and dynamic analysis. This means that the libraries and the clients compile correctly, they are executable and their test suites pass. The dataset is composed of open-source projects that have more than five stars on GitHub. The final dataset contains 395 libraries and 2,874 clients. Additionally, we provide the raw data that we use to create this dataset, such as 34,560 pom.xml files or the complete file list from 34,560 projects. This dataset can be used to study how libraries are used by their clients or as a list of software projects that successfully build. The client's test suite can be used as an additional verification step for code transformation techniques that modify the libraries.


title: Automatic detection of impact craters on Al foils from the Stardust interstellar dust collector using convolutional neural networks

abstract: NASA's Stardust mission utilized a sample collector composed of aerogel and aluminum foil to return cometary and interstellar particles to Earth. Analysis of the aluminum foil begins with locating craters produced by hypervelocity impacts of cometary and interstellar dust. Interstellar dust craters are typically less than one micrometer in size and are sparsely distributed, making them difficult to find. In this paper, we describe a convolutional neural network based on the VGG16 architecture that achieves high specificity and sensitivity in locating impact craters in the Stardust interstellar collector foils. We evaluate its implications for current and future analyses of Stardust samples.


title: Factorization of the Loewner matrix pencil and its consequences

abstract: This paper starts by deriving a factorization of the Loewner matrix pencil that appears in the data-driven modeling approach known as the Loewner framework and explores its consequences. The first is that the associated quadruple constructed from the data yields a model without requiring further processing. The second consequence is related to how sensitive the eigenvalues of the Loewner pencil are to perturbations. Based on an explicit generalized eigenvalue decomposition of this pencil and by making use of perturbation theory of matrix pencils, we explore two types of eigenvalue sensitivities. The first one is defined with respect to unstructured perturbations of the Loewner pencil, while the second one is defined for structured perturbations. We also discuss how the choice of data affects the two sensitivities.


title: A variable corona for GRS 1915+105

abstract: Most models of the low frequency quasi periodic oscillations (QPOs) in low-mass X-ray binaries (LMXBs) explain the dynamical properties of those QPOs. On the other hand, in recent years reverberation models that assume a lamp-post geometry have been successfull in explaining the energy-dependent time lags of the broad-band noise component in stellar mass black-holes and active galactic nuclei. We have recently shown that Comptonisation can explain the spectral-timing properties of the kilo-hertz (kHz) QPOs observed in neutron star (NS) LMXBs. It is therefore worth exploring whether the same family of models would be as successful in explaining the low-frequency QPOs. In this work, we use a Comptonisation model to study the frequency dependence of the phase lags of the type-C QPO in the BH LMXB GRS 1915+105. The phase lags of the QPO in GRS 1915+105 make a transition from hard to soft at a QPO frequency of around 1.8 Hz. Our model shows that at high QPO frequencies a large corona of ~ 100-150 R_g covers most of the accretion disc and makes it 100% feedback dominated, thus producing soft lags. As the observed QPO frequency decreases, the corona gradually shrinks down to around 3-17 R_g, and at 1.8 Hz feedback onto the disc becomes inefficient leading to hard lags. We discuss how changes in the accretion geometry affect the timing properties of the type-C QPO.


title: A New Parameterized Family of Stochastic Particle Flow Filters

abstract: In this paper, we are interested in obtaining answers to the following questions for particle flow filters: Can we provide a theoretical guarantee that particle flow filters give correct results such as unbiased estimates? Are particle flows stable and under what conditions? Can we have one particle flow filter, rather than multiple seemingly different ones? To answer these questions, we first derive a parameterized family of stochastic particle flow filters, in which particle flows are driven by a linear combination of prior knowledge and measurement likelihood information. We then show that several particle flows existing in the literature are special cases of this family. We prove that the particle flows are unbiased under the assumption of linear measurement and Gaussian distributions, and that estimates constructed from the stochastic flows are consistent. We further establish several finite time stability concepts for this new family of stochastic particle flows. The results reported in this paper represent a significant development toward establishing a theoretical foundation for particle flow filters.


title: Low-frequency radio spectra of submillimetre galaxies in the Lockman Hole

abstract: We investigate the radio properties of a sample of 53 sources selected at 850 $\mu$m from the SCUBA-2 Cosmology Legacy Survey using new deep, low-frequency radio imaging of the Lockman Hole field from the Low Frequency Array. Combining these data with additional radio observations from the GMRT and the JVLA, we find a variety of radio spectral shapes and luminosities within our sample despite their similarly bright submillimetre flux densities. We characterise their spectral shapes in terms of multi-band radio spectral indices. Finding strong spectral flattening at low frequencies in ~20% of sources, we investigate the differences between sources with extremely flat low-frequency spectra and those with `normal' radio spectral indices. As there are no other statistically significant differences between the two subgroups of our sample as split by the radio spectral index, we suggest that any differences are undetectable in galaxy-averaged properties that we can observe with our unresolved images, and likely relate to galaxy properties that we cannot resolve, on scales $\lesssim$ 1 kpc. We attribute the observed spectral flattening in the radio to free-free absorption, proposing that those sources with significant low-frequency spectral flattening have a clumpy distribution of star-forming gas. We estimate an average spatial extent of absorbing material of at most several hundred parsecs to produce the levels of absorption observed in the radio spectra. This estimate is consistent with the highest-resolution observations of submillimetre galaxies in the literature, which find examples of non-uniform dust distributions on scales of ~100 pc, with evidence for clumps and knots in the interstellar medium. Additionally, we find two bright (> 6 mJy) submm sources undetected at all other wavelengths. We speculate that these objects may be very high redshift sources, likely residing at z > 4.


title: A note on damped wave equations with a nonlinear dissipation in non-cylindrical domains

abstract: In this paper, we study the large time behavior of a class of wave equation with a nonlinear dissipation in non-cylindrical domains. The result we obtained here relaxes the conditions for the nonlinear term coefficients (in precise, that is $\beta(t)|u|^\rho u$) in \cite{alb} and \cite{ha} (which require $\beta(t)$ to be a constant or $\beta(t)$ to be decreasing with time $t$) and has less restriction for the defined regions.


title: Pulse shape simulation of p-type Point Contact Germanium Detector for The MAJORONA DEMONSTRATOR Experiment

abstract: In the rare physics events search, liquid and solid detector plays an important role. Its analysis is based on the detection of excess events over the expected background or on the detection of an annual event rate modulation. Germanium detectors sensitivities have been demonstrated as efficient means to probe Weakly Interacting Massive Particles (WIMPs). Germanium detectors having sensitivity less than 100 eV opens new platform for the investigation of neutrinos and dark matter physics. There are various working configuration available for germanium detectors such as semi-coaxial, true coaxial and point contact. The point contact germanium detectors have been widely used for dark matter and low energy neutrino searches because of their sub KeV energy threshold and low electronic noise. Pulse shape of P-type point contact (P-PC) germanium detector has been simulated in the present work and compared to the pulse shape of Majorana Demonstrator(MJD).


title: Scalar isoscalar mesons and the scalar glueball from radiative $J/\psi$ decays

abstract: A coupled-channel analysis of BESIII data on radiative $J/\psi$ decays into $\pi\pi$, $K\bar K$, $\eta\eta$ and $\omega\phi$ has been performed. The partial-wave amplitude is constrained by a large number of further data. The analysis finds ten isoscalar scalar mesons. Their masses, widths and decay modes are determined. The scalar mesons are interpreted as mainly SU(3)-singlet and mainly octet states. Octet isoscalar scalar states are observed with significant yields only in the 1500-2100\,MeV mass region. Singlet scalar mesons are produced over a wide mass range but their yield peaks in the same mass region. The peak is interpreted as scalar glueball. Its mass and width are determined to $M=1865$\er25$^{+10}_{-30}$ {\rm MeV} and $\Gamma= 370$\er$50^{+30}_{-20}$ {\rm MeV}, its yield in radiative $J/\psi$ decays to ($5.8\pm 1.0)\,10^{-3}$.


title: Quantization of Calogero-Painlev\'e System and Multi-Particle Quantum Painlev\'e Equations II-VI

abstract: We consider the isomonodromic formulation of the Calogero-Painlev\'e multi-particle systems and proceed to their canonical quantization. We then proceed to the quantum Hamiltonian reduction on a special representation to radial variables, in analogy with the classical case and also with the theory of quantum Calogero equations. This quantized version is compared to the generalization of a result of Nagoya on integral representations of certain solutions of the quantum Painlev\'e equations. We also provide multi-particle generalizations of these integral representations.


title: Towards a Systematic Engineering of Industrial Domain-Specific Language

abstract: Domain-Specific Languages (DSLs) help practitioners in contributing solutions to challenges of specific domains. The efficient development of user-friendly DSLs suitable for industrial practitioners with little expertise in modelling still is challenging. For such practitioners, who often do not model on a daily basis, there is a need to foster reduction of repetitive modelling tasks and providing simplified visual representations of DSL parts. For industrial language engineers, there is no methodical support for providing such guidelines or documentation as part of reusable language modules. Previous research either addresses the reuse of languages or guidelines for modelling. For the efficient industrial deployment of DSLs, their combination is essential: the efficient engineering of DSLs from reusable modules that feature integrated documentation and guidelines for industrial practitioners. To solve these challenges, we propose a systematic approach for the industrial engineering of DSLs based on the concept of reusable DSL Building Blocks, which rests on several years of experience in the industrial engineering of DSLs and their deployment to various organizations. We investigated our approach via focus group methods consisting of five participants from industry and research qualitatively. Ultimately, DSL Building Blocks support industrial language engineers in developing better usable DSLs and industrial practitioners in more efficiently achieving their modelling.


title: Accelerating Radiation Therapy Dose Calculation with Nvidia GPUs

abstract: Radiation Treatment Planning (RTP) is the process of planning the appropriate external beam radiotherapy to combat cancer in human patients. RTP is a complex and compute-intensive task, which often takes a long time (several hours) to compute. Reducing this time allows for higher productivity at clinics and more sophisticated treatment planning, which can materialize in better treatments. The state-of-the-art in medical facilities uses general-purpose processors (CPUs) to perform many steps in the RTP process. In this paper, we explore the use of accelerators to reduce RTP calculating time. We focus on the step that calculates the dose using the Graphics Processing Unit (GPU), which we believe is an excellent candidate for this computation type. Next, we create a highly optimized implementation for a custom Sparse Matrix-Vector Multiplication (SpMV) that operates on numerical formats unavailable in state-of-the-art SpMV libraries (e.g., Ginkgo and cuSPARSE). We show that our implementation is several times faster than the baseline (up-to 4x) and has a higher operational intensity than similar (but different) versions such as Ginkgo and cuSPARSE.


title: Sublinear Average-Case Shortest Paths in Weighted Unit-Disk Graphs

abstract: We consider the problem of computing shortest paths in weighted unit-disk graphs in constant dimension $d$. Although the single-source and all-pairs variants of this problem are well-studied in the plane case, no non-trivial exact distance oracles for unit-disk graphs have been known to date, even for $d=2$. The classical result of Sedgewick and Vitter [Algorithmica '86] shows that for weighted unit-disk graphs in the plane the $A^*$ search has average-case performance superior to that of a standard shortest path algorithm, e.g., Dijkstra's algorithm. Specifically, if the $n$ corresponding points of a weighted unit-disk graph $G$ are picked from a unit square uniformly at random, and the connectivity radius is $r\in (0,1)$, $A^*$ finds a shortest path in $G$ in $O(n)$ expected time when $r=\Omega(\sqrt{\log n/n})$, even though $G$ has $\Theta((nr)^2)$ edges in expectation. In other words, the work done by the algorithm is in expectation proportional to the number of vertices and not the number of edges. In this paper, we break this natural barrier and show even stronger sublinear time results. We propose a new heuristic approach to computing point-to-point exact shortest paths in unit-disk graphs. We analyze the average-case behavior of our heuristic using the same random graph model as used by Sedgewick and Vitter and prove it superior to $A^*$. Specifically, we show that, if we are able to report the set of all $k$ points of $G$ from an arbitrary rectangular region of the plane in $O(k + t(n))$ time, then a shortest path between arbitrary two points of such a random graph on the plane can be found in $O(1/r^2 + t(n))$ expected time. In particular, the state-of-the-art range reporting data structures imply a sublinear expected bound for all $r=\Omega(\sqrt{\log n/n})$ and $O(\sqrt{n})$ expected bound for $r=\Omega(n^{-1/4})$ after only near-linear preprocessing of the point set.


title: Nudging Students Toward Better Software Engineering Behaviors

abstract: Student experiences in large undergraduate Computer Science courses are increasingly impacted by automated systems. Bots, or agents of software automation, are useful for efficiently grading and generating feedback. Current efforts at automation in CS education focus on supporting instructional tasks, but do not address student struggles due to poor behaviors, such as procrastination. In this paper, we explore using bots to improve the software engineering behaviors of students using developer recommendation choice architectures, a framework incorporating behavioral science concepts in recommendations to improve the actions of programmers. We implemented this framework in class-bot, a novel system designed to nudge students to make better choices while working on programming assignments. This work presents a preliminary evaluation integrating this tool in an introductory programming course. Our results show that class-bot is beneficial for improving student development behaviors increasing code quality and productivity.


title: Reconstructing random heterogeneous media through differentiable optimization

abstract: Microstructure reconstruction is a key enabler of process-structure-property linkages, a central topic in materials engineering. Revisiting classical optimization-based reconstruction techniques,they are recognized as a powerful framework to reconstruct random heterogeneous media, especially due to their generality and controllability. The stochasticity of the available approaches is, however, identified as a performance bottleneck. In this work, reconstruction is approached as a differentiable optimization problem, where the error of a generic prescribed descriptor is minimized under consideration of its derivative. As an exemplary descriptor, a suitable differentiable version of spatial correlations is formulated, along with a multigrid scheme to ensure scalability. The applicability of differentiable optimization realized through this descriptor is demonstrated using a wide variety of heterogeneous media, achieving exact statistical equivalence with errors as low as 0 in a short time. We conclude that, while still in an early stage of development, this approach has the potential to significantly alleviate the computational effort currently associated with reconstructing general random heterogeneous media.


title: Moroccan Dialect -Darija- Open Dataset

abstract: Darija Open Dataset (DODa) is an open-source project for the Moroccan dialect. With more than 10,000 entries DODa is arguably the largest open-source collaborative project for Darija-English translation built for Natural Language Processing purposes. In fact, besides semantic categorization, DODa also adopts a syntactic one, presents words under different spellings, offers verb-to-noun and masculine-to-feminine correspondences, contains the conjugation of hundreds of verbs in different tenses, and many other subsets to help researchers better understand and study Moroccan dialect. This data paper presents a description of DODa, its features, how it was collected, as well as a first application in Image Classification using ImageNet labels translated to Darija. This collaborative project is hosted on GitHub platform under MIT's Open-Source license and aims to be a standard resource for researchers, students, and anyone who is interested in Moroccan Dialect


title: Disconnection-Mediated Migration of Interfaces in Microstructures: I. continuum model

abstract: A long-standing goal of materials science is to understand, predict and control the evolution of microstructures in crystalline materials. Most microstructure evolution is controlled by interface motion; hence, the establishment of rigorous interface equations of motion is a universal goal of materials science. We present a new model for the motion of arbitrarily curved interfaces that respects the underlying crystallography of the two phases/domains meeting at the interface and is consistent with microscopic mechanisms of interface motion; i.e., disconnection migration (line defects in the interface with step and dislocation character). We derive the equation of motion for interface migration under the influence of a wide range of driving forces. In Part II of this paper [Salvalaglio, Han and Srolovitz, 2021], we implement the interface model and the equation of motion proposed in this paper in a diffuse interface simulation approach for complex morphology and microstructure evolution.


title: Disconnection-Mediated Migration of Interfaces in Microstructures: II. diffuse interface simulations

abstract: The motion of interfaces is an essential feature of microstructure evolution in crystalline materials. While atomic-scale descriptions provide mechanistic clarity, continuum descriptions are important for understanding microstructural evolution and upon which microscopic features it depends. We develop a microstructure evolution simulation approach that is linked to the underlying microscopic mechanisms of interface migration. We extend the continuum approach describing the disconnection-mediated motion of interfaces introduced in Part I [Han, Srolovitz and Salvalaglio, 2021] to a diffuse interface, phase-field model suitable for large-scale microstructure evolution. A broad range of numerical simulations showcases the capability of the method and the influence of microscopic interface migration mechanisms on microstructure evolution. These include, in particular, the effects of stress and its coupling to interface migration which arises from disconnections, showing how this leads to important differences from classical microstructure evolution represented by mean curvature flow.


title: Examining norms and social expectations surrounding exclusive breastfeeding: Evidence from Mali

abstract: Why do people engage in certain behavior. What are the effects of social expectations and perceptions of community behavior and beliefs on own behavior. Given that proper infant feeding practices are observable and have significant health impacts, we explore the relevance of these questions in the context of exclusive infant breastfeeding behavior using social norms theory. We make use of a primary survey of mothers of children below the age of two years in the Kayes and Sikasso region of Mali, which have a historically lower prevalence of exclusive breastfeeding. The findings from regression estimations, controlling for a host of potential confounding factors, indicate that expectations about the behavior of other community members can strongly predict individual exclusive breastfeeding. Beliefs about approval of the infant feeding behavior of the community though are found to be only modestly associated with it. In addition, mothers who hold false but positive beliefs about the community are found to exclusively breastfeed their kids. Further, using responses from randomly assigned vignettes where we experimentally manipulated the levels of social expectations, our data reveal a strong relationship between perceived prevalence of community level exclusive breastfeeding and individual behavior. This result indicates the existence of a potential causal relationship. We argue that our findings represent an important foundation for the design of policy interventions aimed at altering social expectations, and thus effecting a measurable change in individual behaviors. This type of intervention, by using social norm messaging to end negative behavior, avoids the use of coercive measures to effect behavior change in a cost-effective and efficient way.


title: Infection percolation: A dynamic network model of disease spreading

abstract: Models of disease spreading are critical for predicting infection growth in a population and evaluating public health policies. However, standard models typically represent the dynamics of disease transmission between individuals using macroscopic parameters that do not accurately represent person-to-person variability. To address this issue, we present a dynamic network model that provides a straightforward way to incorporate both disease transmission dynamics at the individual scale as well as the full spatiotemporal history of infection at the population scale. We find that disease spreads through a social network as a traveling wave of infection, followed by a traveling wave of recovery, with the onset and dynamics of spreading determined by the interplay between disease transmission and recovery. We use these insights to develop a scaling theory that predicts the dynamics of infection for diverse diseases and populations. Furthermore, we show how spatial heterogeneities in susceptibility to infection can either exacerbate or quell the spread of disease, depending on its infectivity. Ultimately, our dynamic network approach provides a simple way to model disease spreading that unifies previous findings and can be generalized to diverse diseases, containment strategies, seasonal conditions, and community structures.

categories: physics.soc-ph cond-mat.stat-mech cs.SI nlin.AO q-bio.PE


title: Radical Complexity

abstract: This is an informal and sketchy review of six topical, somewhat unrelated subjects in quantitative finance: rough volatility models; random covariance matrix theory; copulas; crowded trades; high-frequency trading & market stability; and "radical complexity" & scenario based (macro)economics. Some open questions and research directions are briefly discussed.


title: A Robust Tube-Based Smooth-MPC for Robot Manipulator Planning

abstract: Model Predictive Control (MPC) has shown the great performance of target optimization and constraint satisfaction. However, the heavy computation of the Optimal Control Problem (OCP) at each triggering instant brings the serious delay from state sampling to the control signals, which limits the applications of MPC in resource-limited robot manipulator systems over complicated tasks. In this paper, we propose a novel robust tube-based smooth-MPC strategy for nonlinear robot manipulator planning systems with disturbances and constraints. Based on piecewise linearization and state prediction, our control strategy improves the smoothness and optimizes the delay of the control process. By deducing the deviation of the real system states and the nominal system states, we can predict the next real state set at the current instant. And by using this state set as the initial condition, we can solve the next OCP ahead and store the optimal controls based on the nominal system states, which eliminates the delay. Furthermore, we linearize the nonlinear system with a given upper bound of error, reducing the complexity of the OCP and improving the response speed. Based on the theoretical framework of tube MPC, we prove that the control strategy is recursively feasible and closed-loop stable with the constraints and disturbances. Numerical simulations have verified the efficacy of the designed approach compared with the conventional MPC.


title: Deep Contextual Bandits for Fast Neighbor-Aided Initial Access in mmWave Cell-Free Networks

abstract: Access points (APs) in millimeter-wave (mmWave) and sub-THz-based user-centric (UC) networks will have sleep mode functionality. As a result of this, it becomes challenging to solve the initial access (IA) problem when the sleeping APs are activated to start serving users. In this paper, a novel deep contextual bandit (DCB) learning method is proposed to provide instant IA using information from the neighboring active APs. In the proposed approach, beam selection information from the neighboring active APs is used as an input to neural networks that act as a function approximator for the bandit algorithm. Simulations are carried out with realistic channel models generated using the Wireless Insight ray-tracing tool. The results show that the system can respond to dynamic throughput demands with negligible latency compared to the standard baseline 5G IA scheme. The proposed fast beam selection scheme can enable the network to use energy-saving sleep modes without compromising the quality of service due to inefficient IA


title: Weak transport equation on a bounded domain: stability theory apr\'es DiPerna-Lions

abstract: The article studies the transport equation that governes the motion of a fluid in a bounded domain, under the hypothesis of zero velocity at the boundary and supposing the incompressible nature of the fluid. Together with existence and uniqueness results, we study the DiPerna-Lions stability problem, in the case of a bounded domain.


title: Generating Annotated Training Data for 6D Object Pose Estimation in Operational Environments with Minimal User Interaction

abstract: Recently developed deep neural networks achieved state-of-the-art results in the subject of 6D object pose estimation for robot manipulation. However, those supervised deep learning methods require expensive annotated training data. Current methods for reducing those costs frequently use synthetic data from simulations, but rely on expert knowledge and suffer from the "domain gap" when shifting to the real world. Here, we present a proof of concept for a novel approach of autonomously generating annotated training data for 6D object pose estimation. This approach is designed for learning new objects in operational environments while requiring little interaction and no expertise on the part of the user. We evaluate our autonomous data generation approach in two grasping experiments, where we archive a similar grasping success rate as related work on a non autonomously generated data set.


title: Single Underwater Image Restoration by Contrastive Learning

abstract: Underwater image restoration attracts significant attention due to its importance in unveiling the underwater world. This paper elaborates on a novel method that achieves state-of-the-art results for underwater image restoration based on the unsupervised image-to-image translation framework. We design our method by leveraging from contrastive learning and generative adversarial networks to maximize mutual information between raw and restored images. Additionally, we release a large-scale real underwater image dataset to support both paired and unpaired training modules. Extensive experiments with comparisons to recent approaches further demonstrate the superiority of our proposed method.


title: On the orthogonality of generalized eigenspaces for the Ornstein--Uhlenbeck operator

abstract: We study the orthogonality of the generalized eigenspaces of an Ornstein--Uhlenbeck operator $\mathcal L$ in $\mathbb{R}^N$, with drift given by a real matrix $B$ whose eigenvalues have negative real parts. If $B$ has only one eigenvalue, we prove that any two distinct generalized eigenspaces of $\mathcal L$ are orthogonal with respect to the invariant Gaussian measure. Then we show by means of two examples that if $B$ admits distinct eigenvalues, the generalized eigenspaces of $\mathcal L$ may or may not be orthogonal.


title: ShipSRDet: An End-to-End Remote Sensing Ship Detector Using Super-Resolved Feature Representation

abstract: High-resolution remote sensing images can provide abundant appearance information for ship detection. Although several existing methods use image super-resolution (SR) approaches to improve the detection performance, they consider image SR and ship detection as two separate processes and overlook the internal coherence between these two correlated tasks. In this paper, we explore the potential benefits introduced by image SR to ship detection, and propose an end-to-end network named ShipSRDet. In our method, we not only feed the super-resolved images to the detector but also integrate the intermediate features of the SR network with those of the detection network. In this way, the informative feature representation extracted by the SR network can be fully used for ship detection. Experimental results on the HRSC dataset validate the effectiveness of our method. Our ShipSRDet can recover the missing details from the input image and achieves promising ship detection performance.


title: Comprehensive study of out-of-plane transport properties in BaFe$_{2}$As$_{2}$: Three-dimensional electronic state and effect of chemical substitution

abstract: We investigated the out-of-plane transport properties of parent and chemically substituted BaFe$_{2}$As$_{2}$ for various types of substitution. Based on the studies of Hall coefficient and chemical-substitution effect, we have clarified the origin for the unusual temperature dependence of out-of-plane resistivity $\rho_c(T)$ in the high-temperature paramagnetic-tetragonal phase. Electron (hole) carriers have an incoherent (coherent) character, which is responsible for non-metallic (metallic) $\rho_c(T)$. Although both of electron and hole contributions are almost comparable, a slightly larger contribution comes from electrons at high temperatures, while from holes at low temperatures, resulting in a maximum in $\rho_c(T)$. In the low-temperature antiferromagnetic-orthorhombic phase, the major effect of substitution is to increase the residual-resistivity component, as in the case for the in-plane transport. In particular, Co atoms substituted for Fe give rise to strong scattering with large $\mathit{ac}$ anisotropy. We found that K substitution induces a non-metallic behavior in $\rho_c(T)$ at low temperatures, which is likely due to a weakly localized nature along the $c$-axis direction.


title: Control of band structure of FeSe single crystals via biaxial strain

abstract: We performed systematic transport measurements on FeSe single crystals with applying in-plane biaxial strain $\varepsilon$ ranging from -0.96% to 0.23%. Biaxial strain was introduced by firmly gluing samples to various substrate materials with different thermal expansion. With increasing $\varepsilon$, structural and superconducting transition temperatures monotonically increased and decreased, respectively. We analyzed magneto-transport results using a compensated three-carrier model. The evaluated densities of hole and electron carriers systematically changed with strain. This indicates that we succeeded in controlling the band structure of single-crystalline FeSe.


title: Imaging non-radiative point defects buried in quantum wells using cathodoluminescence

abstract: Crystallographic point defects (PDs) can dramatically decrease the efficiency of optoelectronic semiconductor devices, many of which are based on quantum well (QW) heterostructures. However, spatially resolving individual non-radiative PDs buried in such QWs has so far not been demonstrated. Here, using high-resolution cathodoluminescence (CL) and a specific sample design, we spatially resolve, image, and analyse non-radiative PDs in InGaN/GaN QWs. We identify two different types of PD by their contrasting behaviour with temperature, and measure their densities from $10^{14}$ cm$^{-3}$ to as high as $10^{16}$ cm$^{-3}$. Our CL images clearly illustrate the interplay between PDs and carrier dynamics in the well: increasing PD concentration severely limits carrier diffusion lengths, while a higher carrier density suppresses the non-radiative behaviour of PDs. The results in this study are readily interpreted directly from CL images, and represent a significant advancement in nanoscale PD analysis.


title: Gorenstein weak global dimension is symmetric

abstract: We study the Gorenstein weak global dimension of associative rings and its relation to the Gorenstein global dimension. In particular, we prove the conjecture that the Gorenstein weak global dimension is a left-right symmetric invariant -- just like the (absolute) weak global dimension.


title: Reachable Distance Function for KNN Classification

abstract: Distance function is a main metrics of measuring the affinity between two data points in machine learning. Extant distance functions often provide unreachable distance values in real applications. This can lead to incorrect measure of the affinity between data points. This paper proposes a reachable distance function for KNN classification. The reachable distance function is not a geometric direct-line distance between two data points. It gives a consideration to the class attribute of a training dataset when measuring the affinity between data points. Concretely speaking, the reachable distance between data points includes their class center distance and real distance. Its shape looks like "Z", and we also call it a Z distance function. In this way, the affinity between data points in the same class is always stronger than that in different classes. Or, the intraclass data points are always closer than those interclass data points. We evaluated the reachable distance with experiments, and demonstrated that the proposed distance function achieved better performance in KNN classification.


title: Accuracy Gains from Privacy Amplification Through Sampling for Differential Privacy

abstract: Recent research in differential privacy demonstrated that (sub)sampling can amplify the level of protection. For example, for $\epsilon$-differential privacy and simple random sampling with sampling rate $r$, the actual privacy guarantee is approximately $r\epsilon$, if a value of $\epsilon$ is used to protect the output from the sample. In this paper, we study whether this amplification effect can be exploited systematically to improve the accuracy of the privatized estimate. Specifically, assuming the agency has information for the full population, we ask under which circumstances accuracy gains could be expected, if the privatized estimate would be computed on a random sample instead of the full population. We find that accuracy gains can be achieved for certain regimes. However, gains can typically only be expected, if the sensitivity of the output with respect to small changes in the database does not depend too strongly on the size of the database. We only focus on algorithms that achieve differential privacy by adding noise to the final output and illustrate the accuracy implications for two commonly used statistics: the mean and the median. We see our research as a first step towards understanding the conditions required for accuracy gains in practice and we hope that these findings will stimulate further research broadening the scope of differential privacy algorithms and outputs considered.


title: Molecular spin qudits for quantum simulation of light-matter interactions

abstract: We show that molecular spin qudits provide an ideal platform to simulate the quantum dynamics of photon fields strongly interacting with matter. The basic unit of the proposed molecular quantum simulator can be realized by a simple dimer of a spin 1/2 and a spin $S$ transition metal ion, solely controlled by microwave pulses. The spin $S$ ion is exploited to encode the photon field in a flexible architecture, which enables the digital simulation of a wide range of spin-boson models much more efficiently than by using a multi-qubit register. The effectiveness of our proposal is demonstrated by numerical simulations using realistic molecular parameters, whose prerequisites delineating possible chemical approaches are also discussed.