title: A fidelity measure for quantum states based on the matrix geometric mean

abstract: Uhlmann's fidelity function is one of the most widely used similarity measures in quantum theory. One definition of this function is that it is the minimum classical fidelity associated with a quantum-to-classical measurement procedure of two quantum states. In 2010, Matsumoto introduced another fidelity function which is dual to Uhlmann's in the sense that it is the maximimum classical fidelity associated with a classical-to-quantum preparation procedure for two quantum states. Matsumoto's fidelity can also be defined using the well-established notion definition of the matrix geometric mean. In this work, we examine Matsumoto's fidelity through the lens of semidefinite programming to give simple proofs that it possesses many desirable properties for a similarity measure, including monotonicity under quantum channels, joint concavity, and unitary invariance. Finally, we provide a geometric interpretation of this fidelity in terms of the Riemannian space of positive definite matrices, and show how this picture can be useful in understanding some of its peculiar properties.


title: MomentumRNN: Integrating Momentum into Recurrent Neural Networks

abstract: Designing deep neural networks is an art that often involves an expensive search over candidate architectures. To overcome this for recurrent neural nets (RNNs), we establish a connection between the hidden state dynamics in an RNN and gradient descent (GD). We then integrate momentum into this framework and propose a new family of RNNs, called {\em MomentumRNNs}. We theoretically prove and numerically demonstrate that MomentumRNNs alleviate the vanishing gradient issue in training RNNs. We study the momentum long-short term memory (MomentumLSTM) and verify its advantages in convergence speed and accuracy over its LSTM counterpart across a variety of benchmarks. We also demonstrate that MomentumRNN is applicable to many types of recurrent cells, including those in the state-of-the-art orthogonal RNNs. Finally, we show that other advanced momentum-based optimization methods, such as Adam and Nesterov accelerated gradients with a restart, can be easily incorporated into the MomentumRNN framework for designing new recurrent cells with even better performance. The code is available at https://github.com/minhtannguyen/MomentumRNN.


title: Guided Mid-IR and Near-IR Light within a Hybrid Hyperbolic-Material/Silicon Waveguide Heterostructure

abstract: Silicon waveguides have enabled large-scale manipulation and processing of near-infrared optical signals on chip. Yet, expanding the bandwidth of guided waves to other frequencies would further increase the functionality of silicon as a photonics platform. Frequency multiplexing by integrating additional architectures is one approach to the problem, but this is challenging to design and integrate within the existing form factor due to scaling with the free-space wavelength. Here, we demonstrate that a hexagonal boron nitride (hBN)/silicon hybrid waveguide can enable dual-band operation at both mid-infrared (6.5-7.0 um) and telecom (1.55 um) frequencies, respectively. Our device is realized via lithography-free transfer of hBN onto a silicon waveguide, maintaining near-infrared operation, while mid-infrared waveguiding of the hyperbolic phonon polaritons (HPhPs) in hBN is induced by the index contrast between the silicon waveguide and the surrounding air, thereby eliminating the need for deleterious etching of the hBN. We verify the behavior of HPhP waveguiding in both straight and curved trajectories, and validate their propagation characteristics within an analytical waveguide theoretical framework. This approach exemplifies a generalizable approach based on integrating hyperbolic media with silicon photonics for realizing frequency multiplexing in on-chip photonic systems.


title: Maximum $w$-cyclic holely group divisible packings with block size three and applications to optical orthogonal codes

abstract: In this paper we investigate combinatorial constructions for $w$-cyclic holely group divisible packings with block size three (briefly by $3$-HGDPs). For any positive integers $u,v,w$ with $u\equiv0,1~(\bmod~3)$, the exact number of base blocks of a maximum $w$-cyclic $3$-HGDP of type $(u,w^v)$ is determined. This result is used to determine the exact number of codewords in a maximum three-dimensional $(u\times v\times w,3,1)$ optical orthogonal code with at most one optical pulse per spatial plane and per wavelength plane.


title: Incorporating User Micro-behaviors and Item Knowledge into Multi-task Learning for Session-based Recommendation

abstract: Session-based recommendation (SR) has become an important and popular component of various e-commerce platforms, which aims to predict the next interacted item based on a given session. Most of existing SR models only focus on exploiting the consecutive items in a session interacted by a certain user, to capture the transition pattern among the items. Although some of them have been proven effective, the following two insights are often neglected. First, a user's micro-behaviors, such as the manner in which the user locates an item, the activities that the user commits on an item (e.g., reading comments, adding to cart), offer fine-grained and deep understanding of the user's preference. Second, the item attributes, also known as item knowledge, provide side information to model the transition pattern among interacted items and alleviate the data sparsity problem. These insights motivate us to propose a novel SR model MKM-SR in this paper, which incorporates user Micro-behaviors and item Knowledge into Multi-task learning for Session-based Recommendation. Specifically, a given session is modeled on micro-behavior level in MKM-SR, i.e., with a sequence of item-operation pairs rather than a sequence of items, to capture the transition pattern in the session sufficiently. Furthermore, we propose a multi-task learning paradigm to involve learning knowledge embeddings which plays a role as an auxiliary task to promote the major task of SR. It enables our model to obtain better session representations, resulting in more precise SR recommendation results. The extensive evaluations on two benchmark datasets demonstrate MKM-SR's superiority over the state-of-the-art SR models, justifying the strategy of incorporating knowledge learning.


title: Potential Field Guided Actor-Critic Reinforcement Learning

abstract: In this paper, we consider the problem of actor-critic reinforcement learning. Firstly, we extend the actor-critic architecture to actor-critic-N architecture by introducing more critics beyond rewards. Secondly, we combine the reward-based critic with a potential-field-based critic to formulate the proposed potential field guided actor-critic reinforcement learning approach (actor-critic-2). This can be seen as a combination of the model-based gradients and the model-free gradients in policy improvement. State with large potential field often contains a strong prior information, such as pointing to the target at a long distance or avoiding collision by the side of an obstacle. In this situation, we should trust potential-field-based critic more as policy evaluation to accelerate policy improvement, where action policy tends to be guided. For example, in practical application, learning to avoid obstacles should be guided rather than learned by trial and error. State with small potential filed is often lack of information, for example, at the local minimum point or around the moving target. At this time, we should trust reward-based critic as policy evaluation more to evaluate the long-term return. In this case, action policy tends to explore. In addition, potential field evaluation can be combined with planning to estimate a better state value function. In this way, reward design can focus more on the final stage of reward, rather than reward shaping or phased reward. Furthermore, potential field evaluation can make up for the lack of communication in multi-agent cooperation problem, i.e., multi-agent each has a reward-based critic and a relative unified potential-field-based critic with prior information. Thirdly, simplified experiments on predator-prey game demonstrate the effectiveness of the proposed approach.


title: Algebraic stability theorem for derived categories of zigzag persistence modules

abstract: We study distances on zigzag persistence modules from the viewpoint of derived categories and Auslander--Reiten quivers. The derived category of ordinary persistence modules is derived equivalent to that of arbitrary zigzag persistence modules, depending on a classical tilting module. Through this derived equivalence, we define and compute distances on the derived category of arbitrary zigzag persistence modules and prove an algebraic stability theorem. We also compare our distance with the distance for purely zigzag persistence modules introduced by Botnan--Lesnick and the sheaf-theoretic convolution distance due to Kashiwara--Schapira.


title: Finite-thickness effect and spin polarization of the even-denominator fractional quantum Hall states

abstract: The spin-polarized even-denominator fractional quantum Hall (FQH) states in the second Landau level (LL), i.e. 5/2 and 7/2, may possess novel quasi-particle excitations obeying non-Abelian statistics. However, the spin polarization of the 7/2 FQH state has not been investigated experimentally and the spin polarization of the 5/2 FQH state from tilted field experiments remains controversial. Using a piezo-driven sample rotator with the lowest electron temperature down to 25 mK, we studied the energy gap of the even-denominator FQH states in the second LL by precise control of the tilted angles with a resolution less than 0.1{\deg}. We observed two different energy gap dependences on the in-plane magnetic field for 5/2, 7/2, other FQH states (7/3 and 8/3) in the second LL and reentrant integer quantum Hall (RIQH) states in the third LL. Though the transition fields vary from states, their corresponding in-plane magnetic lengths are comparable to the quantum well thickness of the sample, which may result from the influence of the finite-thickness effect. At low in-plane magnetic fields, before the conjectured finite-thickness effect starts to dominate, the energy gaps of both 5/2 and 7/2 states show a non-decreasing behavior, supporting spin-polarized ground states. Our results also suggest that the 7/3, 8/3 FQH states, and the RIQH states in the third LL are spin-polarized or partially spin-polarized.


title: Learning Bayesian Networks with Annealing Machine

abstract: Recent studies have reported that annealing machines are capable of solving combinatorial optimization problems with high accuracy. Annealing machines can potentially be applied to score-based Bayesian network structure learning. However, the bit capacity of an annealing machine is currently limited. To utilize the annealing technology, converting score-based learning problems into quadratic unconstrained binary optimizations within the bit capacity is necessary. In this paper, we propose an efficient conversion method with the advanced identification of candidate parent sets and their decomposition. We also provide an integer programming problem to find the decomposition that minimizes the number of required bits. Experimental results on $7$ benchmark datasets with variables from $75$ to $223$ show that our approach requires less bits than the $100$K bit capacity of the fourth-generation Fujitsu Digital Annealer, a fully coupled annealing machine developed with semiconductor technology. Moreover, we demonstrate that the Digital Annealer with our conversion method outperforms existing algorithms on score maximization. These results highlight the utility of annealing processors in learning Bayesian networks.


title: Classical inequalities for pseudo-integral

abstract: In this paper, we have derived certain classical inequalities, namely, Young's, H\"older's, Minkowski's and Hermite-Hadamard inequalities for pseudo-integral (also known as $g$-integral). For Young's, H\"older's, Minkowski's inequalities, both the cases $p>1$ as well as $p<1,\,p\ne 0$ have been covered. Moreover, in the case of Hermite-Hadamard inequality, a refinement has also been proved and as a special case, $g$-analogue of geometric-logarithmic-arithmetic inequality has been deduced.


title: Characterising authors on the extent of their paper acceptance: A case study of the Journal of High Energy Physics

abstract: New researchers are usually very curious about the recipe that could accelerate the chances of their paper getting accepted in a reputed forum (journal/conference). In search of such a recipe, we investigate the profile and peer review text of authors whose papers almost always get accepted at a venue (Journal of High Energy Physics in our current work). We find authors with high acceptance rate are likely to have a high number of citations, high $h$-index, higher number of collaborators etc. We notice that they receive relatively lengthy and positive reviews for their papers. In addition, we also construct three networks -- co-reviewer, co-citation and collaboration network and study the network-centric features and intra- and inter-category edge interactions. We find that the authors with high acceptance rate are more `central' in these networks; the volume of intra- and inter-category interactions are also drastically different for the authors with high acceptance rate compared to the other authors. Finally, using the above set of features, we train standard machine learning models (random forest, XGBoost) and obtain very high class wise precision and recall. In a followup discussion we also narrate how apart from the author characteristics, the peer-review system might itself have a role in propelling the distinction among the different categories which could lead to potential discrimination and unfairness and calls for further investigation by the system admins.


title: Exploring open cluster properties with Gaia and LAMOST

abstract: In Gaia DR2, the unprecedented high-precision level reached in sub-mas for astrometry and mmag for photometry. Using cluster members identified with these astrometry and photometry in Gaia DR2, we can obtain a reliable determination of cluster properties. However, because of the shortcoming of Gaia spectroscopic observation in dealing with densely crowded cluster region, the number of radial velocity and metallicity for cluster member stars from Gaia DR2 is still lacking. In this study, we aim to improve the cluster properties by combining the LAMOST spectra. In particular, we provide the list of cluster members with spectroscopic parameters as an add-value catalog in LAMOST DR5, which can be used to perform detailed study for a better understanding on the stellar properties, by using their spectra and fundamental properties from the host cluster. We cross-matched the spectroscopic catalog in LAMOST DR5 with the identified cluster members in Cantat-Gaudin et al.2018 and then used members with spectroscopic parameters to derive statistical properties of open clusters. We obtained a list of 8811 members with spectroscopic parameters and a catalog of 295 cluster properties. In addition, we study the radial and vertical metallicity gradient and age-metallicity relation with the compiled open clusters as tracers, finding slopes of -0.053$\pm$0.004 dex kpc$^{-1}$, -0.252$\pm$0.039 dex kpc$^{-1}$ and 0.022$\pm$0.008 dex Gyr$^{-1}$, respectively. Both slopes of metallicity distribution relation for young clusters (0.1 Gyr < Age < 2 Gyr) and the age-metallicity relation for clusters within 6 Gyr are consistent with literature results. In order to fully study the chemical evolution history in the disk, more spectroscopic observations for old and distant open clusters are needed for further investigation.


title: Longitudinal Self-Supervised Learning

abstract: Machine learning analysis of longitudinal neuroimaging data is typically based on supervised learning, which requires a large number of ground-truth labels to be informative. As ground-truth labels are often missing or expensive to obtain in neuroscience, we avoid them in our analysis by combing factor disentanglement with self-supervised learning to identify changes and consistencies across the multiple MRIs acquired of each individual over time. Specifically, we propose a new definition of disentanglement by formulating a multivariate mapping between factors (e.g., brain age) associated with an MRI and a latent image representation. Then, factors that evolve across acquisitions of longitudinal sequences are disentangled from that mapping by self-supervised learning in such a way that changes in a single factor induce change along one direction in the representation space. We implement this model, named Longitudinal Self-Supervised Learning (LSSL), via a standard autoencoding structure with a cosine loss to disentangle brain age from the image representation. We apply LSSL to two longitudinal neuroimaging studies to highlight its strength in extracting the brain-age information from MRI and revealing informative characteristics associated with neurodegenerative and neuropsychological disorders. Moreover, the representations learned by LSSL facilitate supervised classification by recording faster convergence and higher (or similar) prediction accuracy compared to several other representation learning techniques.


title: Quantum Gravity Witness via Entanglement of Masses: Casimir Screening

abstract: A recently proposed experimental protocol for Quantum Gravity induced Entanglement of Masses (QGEM) requires in principle realizable, but still very ambitious, set of parameters in matter-wave interferometry. Motivated by easing the experimental realization, in this paper, we consider the parameter space allowed by a slightly modified experimental design, which mitigates the Casimir potential between two spherical neutral test-masses by separating the two macroscopic interferometers by a thin conducting plate. Although this set-up will reintroduce a Casimir potential between the conducting plate and the masses, there are several advantages of this design. First, the quantum gravity induced entanglement between the two superposed masses will have no Casimir background. Secondly, the matter-wave interferometry itself will be greatly facilitated by allowing both the mass $10^{-16}-10^{-15}$kg and the superposition size $\Delta x \sim 20 \mu$m to be a one-two order of magnitude smaller than those proposed earlier, and thereby also two orders of magnitude smaller magnetic field gradient of $10^4$Tm$^{-1}$ to create that superposition through the Stern-Gerlach effect. In this context, we will further investigate the collisional decoherences and decoherence due to vibrational modes of the conducting plate.


title: Personalized Demand Response via Shape-Constrained Online Learning

abstract: This paper formalizes a demand response task as an optimization problem featuring a known time-varying engineering cost and an unknown (dis)comfort function. Based on this model, this paper develops a feedback-based projected gradient method to solve the demand response problem in an online fashion, where: i) feedback from the user is leveraged to learn the (dis)comfort function concurrently with the execution of the algorithm; and, ii) measurements of electrical quantities are used to estimate the gradient of the known engineering cost. To learn the unknown function, a shape-constrained Gaussian Process is leveraged; this approach allows one to obtain an estimated function that is strongly convex and smooth. The performance of the online algorithm is analyzed by using metrics such as the tracking error and the dynamic regret. A numerical example is illustrated to corroborate the technical findings.


title: Formal Verification of Access Control Model for My Health Record System

abstract: My Health Record system is the Australian Government's digital health record system that holds My Health Record. My Health Record is a secure online health record containing consumers' health information. The system aims to provide health care professionals with access to key health information, e.g. listing medicines, allergies and key diagnoses; radiology and pathology test results. The system (previously named Personally Controlled Electronic Health Record) enables consumers to decide how to share information with any of their health care providers who are registered and connected to the system. The My Health Record system operates under the Australian legislative framework My Health Records Act 2012. The Act establishes, inter alia, a privacy framework specifying which entities can collect, use and disclose certain information in the system and the penalties that can be imposed on improper collection, use and disclosure of this information. This paper presents the formal specification (from the legislation) and verification of the My Health Record regarding how consumers can control who access the information, and how the system adheres to such access. We rely on the correct-by-construction Event-B method to prove control and access properties of the system.


title: Gradient Projection for Solving Quadratic Programs with Standard Simplex Constraints

abstract: An important method to optimize a function on standard simplex is the active set algorithm, which requires the gradient of the function to be projected onto a hyperplane, with sign constraints on the variables that lie in the boundary of the simplex. We propose a new algorithm to efficiently project the gradient for this purpose. Furthermore, we apply the proposed gradient projection method to quadratic programs (QP) with standard simplex constraints, where gradient projection is used to explore the feasible region and, when we believe the optimal active set is identified, we switch to constrained conjugate gradient to accelerate convergence. Specifically, two different directions of gradient projection are used to explore the simplex, namely, the projected gradient and the reduced gradient. We choose one of the two directions according to the angle between the directions. Moreover, we propose two conditions for guessing the optimal active set heuristically. The first condition is that the working set remains unchanged for many iterations, and the second condition is that the angle between the projected gradient and the reduced gradient is small enough. Based on these strategies, a new active set algorithm for solving quadratic programs on standard simplex is proposed.


title: Effects of heterogeneous self-protection awareness on resource-epidemic coevolution dynamics

abstract: Recent studies have demonstrated that the allocation of individual resources has a significant influence on the dynamics of epidemic spreading. In the real scenario, individuals have a different level of awareness for self-protection when facing the outbreak of an epidemic. To investigate the effects of the heterogeneous self-awareness distribution on the epidemic dynamics, we propose a resource-epidemic coevolution model in this paper. We first study the effects of the heterogeneous distributions of node degree and self-awareness on the epidemic dynamics on artificial networks. Through extensive simulations, we find that the heterogeneity of self-awareness distribution suppresses the outbreak of an epidemic, and the heterogeneity of degree distribution enhances the epidemic spreading. Next, we study how the correlation between node degree and self-awareness affects the epidemic dynamics. The results reveal that when the correlation is positive, the heterogeneity of self-awareness restrains the epidemic spreading. While, when there is a significant negative correlation, strong heterogeneous or strong homogeneous distribution of the self-awareness is not conducive for disease suppression. We find an optimal heterogeneity of self-awareness, at which the disease can be suppressed to the most extent. Further research shows that the epidemic threshold increases monotonously when the correlation changes from most negative to most positive, and a critical value of the correlation coefficient is found. When the coefficient is below the critical value, an optimal heterogeneity of self-awareness exists; otherwise, the epidemic threshold decreases monotonously with the decline of the self-awareness heterogeneity. At last, we verify the results on four typical real-world networks and find that the results on the real-world networks are consistent with those on the artificial network.


title: Does Unsupervised Architecture Representation Learning Help Neural Architecture Search?

abstract: Existing Neural Architecture Search (NAS) methods either encode neural architectures using discrete encodings that do not scale well, or adopt supervised learning-based methods to jointly learn architecture representations and optimize architecture search on such representations which incurs search bias. Despite the widespread use, architecture representations learned in NAS are still poorly understood. We observe that the structural properties of neural architectures are hard to preserve in the latent space if architecture representation learning and search are coupled, resulting in less effective search performance. In this work, we find empirically that pre-training architecture representations using only neural architectures without their accuracies as labels considerably improve the downstream architecture search efficiency. To explain these observations, we visualize how unsupervised architecture representation learning better encourages neural architectures with similar connections and operators to cluster together. This helps to map neural architectures with similar performance to the same regions in the latent space and makes the transition of architectures in the latent space relatively smooth, which considerably benefits diverse downstream search strategies.


title: Non-parallel voice conversion based on source-to-target direct mapping

abstract: Recent works of utilizing phonetic posteriograms (PPGs) for non-parallel voice conversion have significantly increased the usability of voice conversion since the source and target DBs are no longer required for matching contents. In this approach, the PPGs are used as the linguistic bridge between source and target speaker features. However, this PPG-based non-parallel voice conversion has some limitation that it needs two cascading networks at conversion time, making it less suitable for real-time applications and vulnerable to source speaker intelligibility at conversion stage. To address this limitation, we propose a new non-parallel voice conversion technique that employs a single neural network for direct source-to-target voice parameter mapping. With this single network structure, the proposed approach can reduce both conversion time and number of network parameters, which can be especially important factors in embedded or real-time environments. Additionally, it improves the quality of voice conversion by skipping the phone recognizer at conversion stage. It can effectively prevent possible loss of phonetic information the PPG-based indirect method suffers. Experiments show that our approach reduces number of network parameters and conversion time by 41.9% and 44.5%, respectively, with improved voice similarity over the original PPG-based method.


title: Enhanced Stability of Antiferromagnetic Skyrmion during Its Motion by Anisotropic Dzyaloshinskii Moriya Interaction

abstract: Searching for new methods to enhance the stability of antiferromagnetic (AFM) skyrmion during its motion is an important issue for AFM spintronic devices. Herein, the spin polarized current-induced dynamics of a distorted AFM skyrmion is numerically studied, based on the Landau Lifshitz Gilbert simulations of the model with an anisotropic Dzyaloshinskii Moriya (DM) interaction. It is demonstrated that the DM interaction anisotropy induces the skyrmion deformation, which suppresses the distortion during the motion and enhances the stability of the skyrmion. Moreover, the effect of the DM interaction anisotropy on the skyrmion velocity is investigated in detail, and the simulated results are further explained by Thiele theory. This work unveils a promising strategy to enhance the stability and the maximum velocity of AFM skyrmion, benefiting future spintronic applications.


title: Targeted Pandemic Containment Through Identifying Local Contact Network Bottlenecks

abstract: Decision-making about pandemic mitigation often relies upon simulation modelling. Models of disease transmission through networks of contacts--between individuals or between population centres--are increasingly used for these purposes. Real-world contact networks are rich in structural features that influence infection transmission, such as tightly-knit local communities that are weakly connected to one another. In this paper, we propose a new flow-based edge-betweenness centrality method for detecting bottleneck edges that connect nodes in contact networks. In particular, we utilize convex optimization formulations based on the idea of diffusion with p-norm network flow. Using simulation models of COVID-19 transmission through real network data at both individual and county levels, we demonstrate that targeting bottleneck edges identified by the proposed method reduces the number of infected cases by up to 10% more than state-of-the-art edge-betweenness methods. Furthermore, the proposed method is orders of magnitude faster than existing methods.


title: Neural voice cloning with a few low-quality samples

abstract: In this paper, we explore the possibility of speech synthesis from low quality found data using only limited number of samples of target speaker. We try to extract only the speaker embedding from found data of target speaker unlike previous works which tries to train the entire text-to-speech system on found data. Also, the two speaker mimicking approaches which are adaptation and speaker-encoder-based are applied on newly released LibriTTS dataset and previously released VCTK corpus to examine the impact of speaker variety on clarity and target-speaker-similarity .


title: Vulnerable Road User Detection Using Smartphone Sensors and Recurrence Quantification Analysis

abstract: With the fast advancements of the Autonomous Vehicle (AV) industry, detection of Vulnerable Road Users (VRUs) using smartphones is critical for safety applications of Cooperative Intelligent Transportation Systems (C-ITSs). This study explores the use of low-power smartphone sensors and the Recurrence Quantification Analysis (RQA) features for this task. These features are computed over a thresholded similarity matrix extracted from nine channels: accelerometer, gyroscope, and rotation vector in each direction (x, y, and z). Given the high-power consumption of GPS, GPS data is excluded. RQA features are added to traditional time domain features to investigate the classification accuracy when using binary, four-class, and five-class Random Forest classifiers. Experimental results show a promising performance when only using RQA features with a resulted accuracy of 98. 34% and a 98. 79% by adding time domain features. Results outperform previous reported accuracy, demonstrating that RQA features have high classifying capability with respect to VRU detection.


title: Domain-adversarial training of multi-speaker TTS

abstract: Multi-speaker TTS has to learn both linguistic embedding and text embedding to generate speech of desired linguistic content in desired voice. However, it is unclear which characteristic of speech results from speaker and which part from linguistic content. In this paper, text embedding is forced to unlearn speaker dependent characteristic using gradient reversal layer to auxiliary speaker classifier that we introduce. We train a speaker classifier using angular margin softmax loss. In subjective evaluation, it is shown that the adversarial training of text embedding for unilingual multi-speaker TTS results in 39.9% improvement on similarity MOS and 40.1% improvement on naturalness MOS.


title: A Drone-based Networked System and Methods for Combating Coronavirus Disease (COVID-19) Pandemic

abstract: Coronavirus disease (COVID-19) is an infectious disease caused by a newly discovered coronavirus. It is similar to influenza viruses and raises concerns through alarming levels of spread and severity resulting in an ongoing pandemic worldwide. Within eight months (by August 2020), it infected 24.0 million persons worldwide and over 824 thousand have died. Drones or Unmanned Aerial Vehicles (UAVs) are very helpful in handling the COVID-19 pandemic. This work investigates the drone-based systems, COVID-19 pandemic situations, and proposes an architecture for handling pandemic situations in different scenarios using real-time and simulation-based scenarios. The proposed architecture uses wearable sensors to record the observations in Body Area Networks (BANs) in a push-pull data fetching mechanism. The proposed architecture is found to be useful in remote and highly congested pandemic areas where either the wireless or Internet connectivity is a major issue or chances of COVID-19 spreading are high. It collects and stores the substantial amount of data in a stipulated period and helps to take appropriate action as and when required. In real-time drone-based healthcare system implementation for COVID-19 operations, it is observed that a large area can be covered for sanitization, thermal image collection, and patient identification within a short period (2 KMs within 10 minutes approx.) through aerial route. In the simulation, the same statistics are observed with an addition of collision-resistant strategies working successfully for indoor and outdoor healthcare operations. Further, open challenges are identified and promising research directions are highlighted.


title: Time-Dependent Generator Coordinate Method for Many-Particle Tunneling

abstract: It has been known that the time-dependent Hartree-Fock (TDHF) method, or the time-dependent density functional theory (TDDFT), fails to describe many-body quantum tunneling. We overcome this problem by superposing a few time-dependent Slater determinants with the time-dependent generator coordinate method (TDGCM). We apply this method to scattering of two $\alpha$ particles in one dimension, and demonstrate that the TDGCM method yields a finite tunneling probability even at energies below the Coulomb barrier, at which the tunneling probability is exactly zero in the TDHF. This is the first case in which a many-particle tunneling is simulated with a microscopic real-time approach.


title: Smartphone Transportation Mode Recognition Using a Hierarchical Machine Learning Classifier and Pooled Features From Time and Frequency Domains

abstract: This paper develops a novel two-layer hierarchical classifier that increases the accuracy of traditional transportation mode classification algorithms. This paper also enhances classification accuracy by extracting new frequency domain features. Many researchers have obtained these features from global positioning system data; however, this data was excluded in this paper, as the system use might deplete the smartphone's battery and signals may be lost in some areas. Our proposed two-layer framework differs from previous classification attempts in three distinct ways: 1) the outputs of the two layers are combined using Bayes' rule to choose the transportation mode with the largest posterior probability; 2) the proposed framework combines the new extracted features with traditionally used time domain features to create a pool of features; and 3) a different subset of extracted features is used in each layer based on the classified modes. Several machine learning techniques were used, including k-nearest neighbor, classification and regression tree, support vector machine, random forest, and a heterogeneous framework of random forest and support vector machine. Results show that the classification accuracy of the proposed framework outperforms traditional approaches. Transforming the time domain features to the frequency domain also adds new features in a new space and provides more control on the loss of information. Consequently, combining the time domain and the frequency domain features in a large pool and then choosing the best subset results in higher accuracy than using either domain alone. The proposed two-layer classifier obtained a maximum classification accuracy of 97.02%.


title: SGD with shuffling: optimal rates without component convexity and large epoch requirements

abstract: We study without-replacement SGD for solving finite-sum optimization problems. Specifically, depending on how the indices of the finite-sum are shuffled, we consider the RandomShuffle (shuffle at the beginning of each epoch) and SingleShuffle (shuffle only once) algorithms. First, we establish minimax optimal convergence rates of these algorithms up to poly-log factors. Notably, our analysis is general enough to cover gradient dominated nonconvex costs, and does not rely on the convexity of individual component functions unlike existing optimal convergence results. Secondly, assuming convexity of the individual components, we further sharpen the tight convergence results for RandomShuffle by removing the drawbacks common to all prior arts: large number of epochs required for the results to hold, and extra poly-log factor gaps to the lower bound.


title: Constraining gluon density of pions at large $x$ by pion-induced $J/\psi$ production

abstract: The gluon distributions of the pion obtained from various global fits exhibit large variations among them. Within the framework of the color evaporation model, we show that the existing pion-induced $J/\psi$ production data, usually not included in the global fits, can impose useful additional constraints on the pion parton distribution functions (PDFs). In particular, these data can probe the pion's gluon densities at large $x$. Existing pion-induced $J/\psi$ data covering a broad range of beam momenta are compared with next-to-leading-order QCD calculations using various sets of pion PDFs. It is found that $J/\psi$ data measured at forward rapidity and at sufficiently high beam momentum are sensitive to the large-$x$ gluon distribution of pions. The current $J/\psi$ data favor the Sutton-Martin-Roberts-Stirling and Gluck-Reya-Vogt pion PDFs, containing significant gluon content at large $x$.


title: Two sites coherence and visibility

abstract: Wave-particle duality and the superposition of quantum mechanical states furnish quantum mechanics with unique features which distinguishes it from classical mechanics and give it the apparently counter-intuition interpretation. The two principles are responsible for the observation of the interference effects of quantum particles such as electrons, atoms and molecules. Visibility is a measure of the wave nature and can be though of as a "normalized" coherence quantifier. Reduction in the visibility arises from dephasing (decoherence): a process in which the relative phases get partially or totally destroyed leading to domination of the particle nature. We calculate the coherence and visibility of an ensemble of a single electron and two electrons on two sites using the density matrix formulation. For such a system of electrons, visibility and predictability (a particle nature measure) follow an oscillatory-complementary temporal evolution and neither of them depends on the single particle energy.


title: Patterns in Shi tableaux and Dyck paths

abstract: Shi tableaux are special binary fillings of certain Young diagrams which arise in the study of Shi hyperplane arrangements related to classical root systems. For type $A$, the set $\mathcal T$ of Shi tableaux naturally coincides with the set of Dyck paths, for which various notions of patterns have been introduced and studied over the years. In this paper we define a notion of pattern occurrence on $\mathcal T$ which, although it can be regarded as a pattern on Dyck paths, it is motivated by the underlying geometric structure of the tableaux. Our main goal in this work is to study the poset of Shi tableaux defined by pattern-containment. More precisely, we determine explicit formulas for upper and lower covers for each $T\in\mathcal T$, we consider pattern avoidance for the smallest non-trivial tableaux (size 2) and generalize these results to certain tableau of larger size. We conclude with open problems and possible future directions.


title: GA-guided mD-VcMD: A genetic-algorithm-based method for multi-dimensional virtual-system coupled molecular dynamics

abstract: We previously introduced a conformational sampling method, a multi-dimensional virtual-system coupled molecular dynamics (mD-VcMD), to enhance conformational sampling of a biomolecular system by computer simulations. Here, we present a new sampling method, subzone-based mD-VcMD, as an extension of mD-VcMD. Then, we further extend the subzone-based method using genetic algorithm (GA), and named it the GA-based mD-VcMD. Because the conformational space of the biomolecular system is vast, a single simulation cannot sample the conformational space throughout. Then, iterative simulations are performed to increase the sampled region gradually. The new methods have the following advantages: (1) The methods are free from a production run: I.e., all snapshots from all iterations can be used for analyses. (2) The methods are free from fine tuning of a weight function (probability distribution function or potential of mean force). (3) A simple procedure is available to assign a thermodynamic weight to snapshots sampled in spite that the weight function is not used to proceed the iterative simulations. Thus, a canonical ensemble (i.e., a thermally equilibrated ensemble) is generated from the resultant snapshots. (4) If a poorly-sampled region emerges in sampling, selective sampling can be performed focusing on the poorly-sampled region without breaking the proportion of the canonical ensemble. A free-energy landscape of the biomolecular system is obtainable from the canonical ensemble.


title: Planar Rectilinear Drawings of Outerplanar Graphs in Linear Time

abstract: We show how to test in linear time whether an outerplanar graph admits a planar rectilinear drawing, both if the graph has a prescribed plane embedding that the drawing has to respect and if it does not. Our algorithm returns a planar rectilinear drawing if the graph admits one.


title: Topological superconductivity in quasicrystals

abstract: We propose realization of non-Abelian topological superconductivity in two-dimensional quasicrystals by the same mechanism as in crystalline counterparts. Specifically, we study a two-dimensional electron gas in Penrose and Ammann-Beenker quasicrystals with Rashba spin-orbit coupling, perpendicular Zeeman magnetic field, and conventional $s$-wave superconductivity. We find that topological superconductivity with broken time-reversal symmetry is realized in both Penrose and Ammann-Beenker quasicrystals at low filling, where the Bott index is unity. The topological nature of this phase is confirmed by the existence of a zero-energy surface bound state and the chiral propagation of a wave packet projected onto the midgap bound state along the surfaces. Furthermore, we confirm the existence of a single Majorana zero mode each in a vortex at the center of the system and along the surfaces, signifying the non-Abelian character of the system when the Bott index is unity.


title: Enumerating Teams in First-Order Team Logics

abstract: We start the study of the enumeration complexity of different satisfiability problems in first-order team logics. Since many of our problems go beyond DelP, we use a framework for hard enumeration analogous to the polynomial hierarchy, which was recently introduced by Creignou et al. (Discret. Appl. Math. 2019). We show that the problem to enumerate all satisfying teams of a fixed formula in a given first-order structure is DelNP-complete for certain formulas of dependence logic and independence logic. For inclusion logic formulas, this problem is even in DelP. Furthermore, we study the variants of this problems where only maximal, minimal, maximum and minimum solutions, respectively, are considered. For the most part these share the same complexity as the original problem. An exception is the minimum-variant for inclusion logic, which is DelNP-complete.


title: Towards Flexible Device Participation in Federated Learning

abstract: Traditional federated learning algorithms impose strict requirements on the participation rates of devices, which limit the potential reach of federated learning. This paper extends the current learning paradigm to include devices that may become inactive, compute incomplete updates, and depart or arrive in the middle of training. We derive analytical results to illustrate how allowing more flexible device participation can affect the learning convergence when data is not independently and identically distributed (non-IID). We then propose a new federated aggregation scheme that converges even when devices may be inactive or return incomplete updates. We also study how the learning process can adapt to early departures or late arrivals, and analyze their impacts on the convergence.


title: AstroCatR: a Mechanism and Tool for Efficient Time Series Reconstruction of Large-Scale Astronomical Catalogues

abstract: Time series data of celestial objects are commonly used to study valuable and unexpected objects such as extrasolar planets and supernova in time domain astronomy. Due to the rapid growth of data volume, traditional manual methods are becoming extremely hard and infeasible for continuously analyzing accumulated observation data. To meet such demands, we designed and implemented a special tool named AstroCatR that can efficiently and flexibly reconstruct time series data from large-scale astronomical catalogues. AstroCatR can load original catalogue data from Flexible Image Transport System (FITS) files or databases, match each item to determine which object it belongs to, and finally produce time series datasets. To support the high-performance parallel processing of large-scale datasets, AstroCatR uses the extract-transform-load (ETL) preprocessing module to create sky zone files and balance the workload. The matching module uses the overlapped indexing method and an in-memory reference table to improve accuracy and performance. The output of AstroCatR can be stored in CSV files or be transformed other into formats as needed. Simultaneously, the module-based software architecture ensures the flexibility and scalability of AstroCatR. We evaluated AstroCatR with actual observation data from The three Antarctic Survey Telescopes (AST3). The experiments demonstrate that AstroCatR can efficiently and flexibly reconstruct all time series data by setting relevant parameters and configuration files. Furthermore, the tool is approximately 3X faster than methods using relational database management systems at matching massive catalogues.


title: Decorrelated Double Q-learning

abstract: Q-learning with value function approximation may have the poor performance because of overestimation bias and imprecise estimate. Specifically, overestimation bias is from the maximum operator over noise estimate, which is exaggerated using the estimate of a subsequent state. Inspired by the recent advance of deep reinforcement learning and Double Q-learning, we introduce the decorrelated double Q-learning (D2Q). Specifically, we introduce the decorrelated regularization item to reduce the correlation between value function approximators, which can lead to less biased estimation and low variance. The experimental results on a suite of MuJoCo continuous control tasks demonstrate that our decorrelated double Q-learning can effectively improve the performance.


title: Fractional Decomposition Tree Algorithm: A tool for studying the integrality gap of Integer Programs

abstract: We present a new algorithm, Fractional Decomposition Tree (FDT) for finding a feasible solution for an integer program (IP) where all variables are binary. FDT runs in polynomial time and is guaranteed to find a feasible integer solution provided the integrality gap is bounded. The algorithm gives a construction for Carr and Vempala's theorem that any feasible solution to the IP's linear-programming relaxation, when scaled by the instance integrality gap, dominates a convex combination of feasible solutions. FDT is also a tool for studying the integrality gap of IP formulations. We demonstrate that with experiments studying the integrality gap of two problems: optimally augmenting a tree to a 2-edge-connected graph and finding a minimum-cost 2-edge-connected multi-subgraph (2EC). We also give a simplified algorithm, Dom2IP, that more quickly determines if an instance has an unbounded integrality gap. We show that FDT's speed and approximation quality compare well to that of feasibility pump on moderate-sized instances of the vertex cover problem. For a particular set of hard-to-decompose fractional 2EC solutions, FDT always gave a better integer solution than the best previous approximation algorithm (Christofides).


title: Understanding the Role of Training Regimes in Continual Learning

abstract: Catastrophic forgetting affects the training of neural networks, limiting their ability to learn multiple tasks sequentially. From the perspective of the well established plasticity-stability dilemma, neural networks tend to be overly plastic, lacking the stability necessary to prevent the forgetting of previous knowledge, which means that as learning progresses, networks tend to forget previously seen tasks. This phenomenon coined in the continual learning literature, has attracted much attention lately, and several families of approaches have been proposed with different degrees of success. However, there has been limited prior work extensively analyzing the impact that different training regimes -- learning rate, batch size, regularization method-- can have on forgetting. In this work, we depart from the typical approach of altering the learning algorithm to improve stability. Instead, we hypothesize that the geometrical properties of the local minima found for each task play an important role in the overall degree of forgetting. In particular, we study the effect of dropout, learning rate decay, and batch size, on forming training regimes that widen the tasks' local minima and consequently, on helping it not to forget catastrophically. Our study provides practical insights to improve stability via simple yet effective techniques that outperform alternative baselines.


title: Three-dimensional hydrodynamic simulations of the upper atmosphere of $\pi$ Men c: comparison with Ly$\alpha$ transit observations

abstract: Aims: We aim at constraining the conditions of the wind and high-energy emission of the host star reproducing the non-detection of Ly$\alpha$ planetary absorption. Methods: We model the escaping planetary atmosphere, the stellar wind, and their interaction employing a multi-fluid, three-dimensional hydrodynamic code. We assume a planetary atmosphere composed of hydrogen and helium. We run models varying the stellar high-energy emission and stellar mass-loss rate, further computing for each case the Ly$\alpha$ synthetic planetary atmospheric absorption and comparing it with the observations. Results: We find that a non-detection of Ly$\alpha$ in absorption employing the stellar high-energy emission estimated from far-ultraviolet and X-ray data requires a stellar wind with a stellar mass-loss rate about six times lower than solar. This result is a consequence of the fact that, for $\pi$ Men c, detectable Ly$\alpha$ absorption can be caused exclusively by energetic neutral atoms, which become more abundant with increasing the velocity and/or the density of the stellar wind. By considering, instead, that the star has a solar-like wind, the non-detection requires a stellar ionising radiation about four times higher than estimated. This is because, despite the fact that a stronger stellar high-energy emission ionises hydrogen more rapidly, it also increases the upper atmosphere heating and expansion, pushing the interaction region with the stellar wind farther away from the planet, where the planet atmospheric density that remains neutral becomes smaller and the production of energetic neutral atoms less efficient. Conclusions: Comparing the results of our grid of models with what is expected and estimated for the stellar wind and high-energy emission, respectively, we support the idea that the atmosphere of $\pi$ Men c is likely not hydrogen-dominated.


title: Sur la r\'{e}partition jointe de la repr\'{e}sentation d'Ostrowski dans les classes de r\'{e}sidue

abstract: For two distinct integers $m_1,m_2\ge2$, we set $\alpha_1=[0;\overline{1,m_1}]$ and $\alpha_2=[0;\overline{1,m_2}]$ and we denote by $S_{\alpha_1}(n)$ and $S_{\alpha_2}(n)$ respectively the sum of digits functions in the Ostrowski $\alpha_1$ and $\alpha_2-$representations of $n$. Let $b_1,b_2 $ be positive integers satisfying $(b_1,m_1)=1$ and $(b_2,m_2)=1$, we obtain an estimation with an error term $O(N^{1-\delta})$ for the cardinal of the following set $$\Big\{ 0\leq n<N;\ S_{\alpha_1}(n)\equiv a_1\pmod{b_1},\ S_{\alpha_2}(n)\equiv a_2\pmod{b_2}\Big\},$$ for all integers $a_1$ and $a_2.$ Our result should be compared to that of B\'{e}sineau and Kim who treated the case of the $q-$representations in different bases (that are coprimes).


title: The eyes know it: FakeET -- An Eye-tracking Database to Understand Deepfake Perception

abstract: We present \textbf{FakeET}-- an eye-tracking database to understand human visual perception of \emph{deepfake} videos. Given that the principal purpose of deepfakes is to deceive human observers, FakeET is designed to understand and evaluate the ease with which viewers can detect synthetic video artifacts. FakeET contains viewing patterns compiled from 40 users via the \emph{Tobii} desktop eye-tracker for 811 videos from the \textit{Google Deepfake} dataset, with a minimum of two viewings per video. Additionally, EEG responses acquired via the \emph{Emotiv} sensor are also available. The compiled data confirms (a) distinct eye movement characteristics for \emph{real} vs \emph{fake} videos; (b) utility of the eye-track saliency maps for spatial forgery localization and detection, and (c) Error Related Negativity (ERN) triggers in the EEG responses, and the ability of the \emph{raw} EEG signal to distinguish between \emph{real} and \emph{fake} videos.


title: Synchronization Stability Analysis and Enhancement of Grid-Tied Multi-Converter Systems

abstract: Synchronization instability of grid-tied converters during grid faults is a serious concern. For phase-locked loop (PLL) synchronized converters, synchronization stability during grid faults requires that PLL can resynchronize with the power grid. It is of great significance for meeting grid code specifications such as reactive current support. However, the synchronization in the context of high-impedance weak grid is not always readily achieved. On one hand, the terminal voltage of converters is changeable due to being susceptible to the output current. One the other hand, there is interaction between multiple converters interconnected to the grid. Therefore, loss of synchronism (LOS) of converters during grid faults likely occurs. Previous research has investigated the LOS issue of single-converter infinite-bus (SCIB) system, but mostly overlooked the interaction between multiple converters. Previous efforts to address the issue of LOS are also in doubt for multi-converter systems. This paper develops a multi-converter infinite-bus (MCIB) system model and analyzes the root cause of LOS with considering multi-converter interaction during grid faults. Moreover, a feedforward compensated PLL (FFC-PLL) method applicable for multi-converter systems is developed to address the LOS issue. Simulations results verify the effectiveness of the method.


title: Needle in a Haystack: Label-Efficient Evaluation under Extreme Class Imbalance

abstract: Important tasks like record linkage and extreme classification demonstrate extreme class imbalance, with 1 minority instance to every 1 million or more majority instances. Obtaining a sufficient sample of all classes, even just to achieve statistically-significant evaluation, is so challenging that most current approaches yield poor estimates or incur impractical cost. Where importance sampling has been levied against this challenge, restrictive constraints are placed on performance metrics, estimates do not come with appropriate guarantees, or evaluations cannot adapt to incoming labels. This paper develops a framework for online evaluation based on adaptive importance sampling. Given a target performance metric and model for $p(y|x)$, the framework adapts a distribution over items to label in order to maximize statistical precision. We establish strong consistency and a central limit theorem for the resulting performance estimates, and instantiate our framework with worked examples that leverage Dirichlet-tree models. Experiments demonstrate an average MSE superior to state-of-the-art on fixed label budgets.


title: Maximal inequalities for stochastic convolutions and pathwise uniform convergence of time discretisation schemes

abstract: We prove a new Burkholder-Rosenthal type inequality for discrete-time processes taking values in a 2-smooth Banach space. As a first application we prove that if $(S(t,s))_{0\leq s\leq T}$ is a $C_0$-evolution family of contractions on a $2$-smooth Banach space $X$ and $(W_t)_{t\in [0,T]}$ is a cylindrical Brownian motion on a probability space $(\Omega,P)$, then for every $0<p<\infty$ there exists a constant $C_{p,X}$ such that for all progressively measurable processes $g: [0,T]\times \Omega\to X$ the process $(\int_0^t S(t,s)g_sdW_s)_{t\in [0,T]}$ has a continuous modification and $$E\sup_{t\in [0,T]}\Big\| \int_0^t S(t,s)g_sdW_s \Big\|^p\leq C_{p,X}^p \mathbb{E} \Bigl(\int_0^T \| g_t\|^2_{\gamma(H,X)}dt\Bigr)^{p/2}.$$ Moreover, for $2\leq p<\infty$ one may take $C_{p,X} = 10 D \sqrt{p},$ where $D$ is the constant in the definition of $2$-smoothness for $X$. Our result improves and unifies several existing maximal estimates and is even new in case $X$ is a Hilbert space. Similar results are obtained if the driving martingale $g_tdW_t$ is replaced by more general $X$-valued martingales $dM_t$. Moreover, our methods allow for random evolution systems, a setting which appears to be completely new as far as maximal inequalities are concerned. As a second application, for a large class of time discretisation schemes we obtain stability and pathwise uniform convergence of time discretisation schemes for solutions of linear SPDEs $$ du_t = A(t)u_tdt + g_tdW_t, \quad u_0 = 0,$$ Under spatial smoothness assumptions on the inhomogeneity $g$, contractivity is not needed and explicit decay rates are obtained. In the parabolic setting this sharpens several know estimates in the literature; beyond the parabolic setting this seems to provide the first systematic approach to pathwise uniform convergence to time discretisation schemes.


title: Quantum-over-classical Advantage in Solving Multiplayer Games

abstract: We study the applicability of quantum algorithms in computational game theory and generalize some results related to Subtraction games, which are sometimes referred to as one-heap Nim games. In quantum game theory, a subset of Subtraction games became the first explicitly defined class of zero-sum combinatorial games with provable separation between quantum and classical complexity of solving them. For a narrower subset of Subtraction games, an exact quantum sublinear algorithm is known that surpasses all deterministic algorithms for finding solutions with probability $1$. Typically, both Nim and Subtraction games are defined for only two players. We extend some known results to games for three or more players, while maintaining the same classical and quantum complexities: $\Theta\left(n^2\right)$ and $\tilde{O}\left(n^{1.5}\right)$ respectively.


title: RISCuer: A Reliable Multi-UAV Search and Rescue Testbed

abstract: We present the Robotics Intelligent Systems & Control (RISC) Lab multiagent testbed for reliable search and rescue and aerial transport in outdoor environments. The system consists of a team of three multirotor unmanned aerial vehicles (UAVs), which are capable of autonomously searching, picking up, and transporting randomly distributed objects in an outdoor field. The method involves vision based object detection and localization, passive aerial grasping with our novel design, GPS based UAV navigation, and safe release of the objects at the drop zone. Our cooperative strategy ensures safe spatial separation between UAVs at all times and we prevent any conflicts at the drop zone using communication enabled consensus. All computation is performed onboard each UAV. We describe the complete software and hardware architecture for the system and demonstrate its reliable performance using comprehensive outdoor experiments, and by comparing our results with some recent, similar works.


title: Precise Morrey regularity of the weak solutions to a kind of quasilinear systems with discontinuous data

abstract: We consider the Dirichlet problem for a class of quasilinear elliptic systems in domain with irregular boundary. The principal part satisfies componentwise coercivity condition and the nonlinear terms are Carath\'eodory maps having Morrey regularity in $x$ and verifying controlled growth conditions with respect to the other variables. We have obtained boundedness of the weak solution to the problem that permits to apply an iteration procedure in order to find optimal Morrey regularity of its gradient.


title: Early Detection of Retinopathy of Prematurity (ROP) in Retinal Fundus Images Via Convolutional Neural Networks

abstract: Retinopathy of prematurity (ROP) is an abnormal blood vessel development in the retina of a prematurely-born infant or an infant with low birth weight. ROP is one of the leading causes for infant blindness globally. Early detection of ROP is critical to slow down and avert the progression to vision impairment caused by ROP. Yet there is limited awareness of ROP even among medical professionals. Consequently, dataset for ROP is limited if ever available, and is in general extremely imbalanced in terms of the ratio between negative images and positive ones. In this study, we formulate the problem of detecting ROP in retinal fundus images in an optimization framework, and apply state-of-art convolutional neural network techniques to solve this problem. Experimental results based on our models achieve 100 percent sensitivity, 96 percent specificity, 98 percent accuracy, and 96 percent precision. In addition, our study shows that as the network gets deeper, more significant features can be extracted for better understanding of ROP.


title: Multi Layer Neural Networks as Replacement for Pooling Operations

abstract: Pooling operations, which can be calculated at low cost and serve as a linear or nonlinear transfer function for data reduction, are found in almost every modern neural network. Countless modern approaches have already tackled replacing the common maximum value selection and mean value operations, not to mention providing a function that allows different functions to be selected through changing parameters. Additional neural networks are used to estimate the parameters of these pooling functions.Consequently, pooling layers may require supplementary parameters to increase the complexity of the whole model. In this work, we show that one perceptron can already be used effectively as a pooling operation without increasing the complexity of the model. This kind of pooling allows for the integration of multi-layer neural networks directly into a model as a pooling operation by restructuring the data and, as a result, learnin complex pooling operations. We compare our approach to tensor convolution with strides as a pooling operation and show that our approach is both effective and reduces complexity. The restructuring of the data in combination with multiple perceptrons allows for our approach to be used for upscaling, which can then be utilized for transposed convolutions in semantic segmentation.


title: The structure of Zeckendorf expansions

abstract: In this paper we classify the Zeckendorf expansions according to their digit blocks. It turns out that if we consider these digit blocks as labels on the Fibonacci tree, then the numbers ending with a given digit block in their Zeckendorf expansion appear as compound Wythoff sequences in a natural way on this tree. Here the digit blocks consisting of only $0$'s are an exception. We also give a second description of these occurrence sequences as generalized Beatty sequences. Finally, we characterize the numbers with a fixed digit block occurring at an arbitrary fixed position in their Zeckendorf expansions, and determine their densities.


title: Generic Indic Text-to-speech Synthesisers with Rapid Adaptation in an End-to-end Framework

abstract: Building text-to-speech (TTS) synthesisers for Indian languages is a difficult task owing to a large number of active languages. Indian languages can be classified into a finite set of families, prominent among them, Indo-Aryan and Dravidian. The proposed work exploits this property to build a generic TTS system using multiple languages from the same family in an end-to-end framework. Generic systems are quite robust as they are capable of capturing a variety of phonotactics across languages. These systems are then adapted to a new language in the same family using small amounts of adaptation data. Experiments indicate that good quality TTS systems can be built using only 7 minutes of adaptation data. An average degradation mean opinion score of 3.98 is obtained for the adapted TTSes. Extensive analysis of systematic interactions between languages in the generic TTSes is carried out. x-vectors are included as speaker embedding to synthesise text in a particular speaker's voice. An interesting observation is that the prosody of the target speaker's voice is preserved. These results are quite promising as they indicate the capability of generic TTSes to handle speaker and language switching seamlessly, along with the ease of adaptation to a new language.


title: Towards Deeper Graph Neural Networks with Differentiable Group Normalization

abstract: Graph neural networks (GNNs), which learn the representation of a node by aggregating its neighbors, have become an effective computational tool in downstream applications. Over-smoothing is one of the key issues which limit the performance of GNNs as the number of layers increases. It is because the stacked aggregators would make node representations converge to indistinguishable vectors. Several attempts have been made to tackle the issue by bringing linked node pairs close and unlinked pairs distinct. However, they often ignore the intrinsic community structures and would result in sub-optimal performance. The representations of nodes within the same community/class need be similar to facilitate the classification, while different classes are expected to be separated in embedding space. To bridge the gap, we introduce two over-smoothing metrics and a novel technique, i.e., differentiable group normalization (DGN). It normalizes nodes within the same group independently to increase their smoothness, and separates node distributions among different groups to significantly alleviate the over-smoothing issue. Experiments on real-world datasets demonstrate that DGN makes GNN models more robust to over-smoothing and achieves better performance with deeper GNNs.


title: Effective boundary conditions for magnetohydrodynamic flows with thin Hartmann layers

abstract: Here we build some effective boundary conditions to be used in numerical calculations in order to avoid the thin meshing usually required in problems involving Hartmann layers near a locally plane wall. Wall model are provided for both tangential and normal electric current density and velocity. In particular, a condition on the normal derivative of the tangential velocity is derived. A wide variety of problems is covered as the only restriction is that the magnetic Reynolds number has to be large at the scale of the Hartmann layer. The cases of perfectly conducting or insulating wall are examined, as well as the case of a thin conducting wall. The newest result is a condition on the normal velocity accounting for inertial effects in the Hartmann layer.


title: Generation of polarized particle beams at relativistic laser intensities

abstract: The acceleration of polarized electrons, positrons, protons and ions in strong laser and plasma fields is a very attractive option to obtain polarized beams in the multi-MeV range. Recently, there has been substantial progress in the understanding of the dominant mechanisms leading to high degrees of polarization, in the numerical modelling of these processes and in their experimental implementation. This review article presents an overview on the current status of the field, on the concepts of polarized laser-plasma accelerators and of the beam polarimetry.


title: Measure-valued solutions and weak-strong uniqueness for the incompressible inviscid fluid-rigid body interaction

abstract: We consider a coupled system of partial and ordinary differential equations describing the interaction between an isentropic inviscid fluid and a rigid body moving freely inside the fluid. We prove the existence of measure-valued solutions which is generated by the vanishing viscosity limit of incompressible fluid-rigid body interaction system under some physically constitutive relations. Moreover, we show that the measure-value solution coincides with strong solution on the interval of its existence. This relies on the weak-strong uniqueness analysis.


title: Towards Robust Pattern Recognition: A Review

abstract: The accuracies for many pattern recognition tasks have increased rapidly year by year, achieving or even outperforming human performance. From the perspective of accuracy, pattern recognition seems to be a nearly-solved problem. However, once launched in real applications, the high-accuracy pattern recognition systems may become unstable and unreliable, due to the lack of robustness in open and changing environments. In this paper, we present a comprehensive review of research towards robust pattern recognition from the perspective of breaking three basic and implicit assumptions: closed-world assumption, independent and identically distributed assumption, and clean and big data assumption, which form the foundation of most pattern recognition models. Actually, our brain is robust at learning concepts continually and incrementally, in complex, open and changing environments, with different contexts, modalities and tasks, by showing only a few examples, under weak or noisy supervision. These are the major differences between human intelligence and machine intelligence, which are closely related to the above three assumptions. After witnessing the significant progress in accuracy improvement nowadays, this review paper will enable us to analyze the shortcomings and limitations of current methods and identify future research directions for robust pattern recognition.


title: Convolutional neural network based hierarchical autoencoder for nonlinear mode decomposition of fluid field data

abstract: We propose a customized convolutional neural network based autoencoder called a hierarchical autoencoder, which allows us to extract nonlinear autoencoder modes of flow fields while preserving the contribution order of the latent vectors. As preliminary tests, the proposed method is first applied to a cylinder wake at $Re_D$ = 100 and its transient process. It is found that the proposed method can extract the features of these laminar flow fields as the latent vectors while keeping the order of their energy content. The present hierarchical autoencoder is further assessed with a two-dimensional $y-z$ cross-sectional velocity field of turbulent channel flow at $Re_{\tau}$ = 180 in order to examine its applicability to turbulent flows. It is demonstrated that the turbulent flow field can be efficiently mapped into the latent space by utilizing the hierarchical model with a concept of ordered autoencoder mode family. The present results suggest that the proposed concept can be extended to meet various demands in fluid dynamics including reduced order modeling and its combination with linear theory-based methods by using its ability to arrange the order of the extracted nonlinear modes.


title: Generalized Weighted Survival and Failure Entropies and their Dynamic Versions

abstract: The weighted forms of generalized survival and failure entropies of order ($\alpha,\beta$) are proposed and some properties are obtained. We further propose the dynamic versions of weighted generalized survival and failures entropies and obtained some properties and bounds. Characterization for Rayleigh and power distributions are done by dynamic weighted generalized entropies. We further consider the empirical versions of generalized weighted survival and failure entropies and using the difference between theoretical and empirical survival entropies a test for exponentiality is considered.


title: Non-Negative Bregman Divergence Minimization for Deep Direct Density Ratio Estimation

abstract: Density ratio estimation (DRE) is at the core of various machine learning tasks such as anomaly detection and domain adaptation. In existing studies on DRE, methods based on Bregman divergence (BD) minimization have been extensively studied. However, BD minimization when applied with highly flexible models, such as deep neural networks, tends to suffer from what we call train-loss hacking, which is a source of overfitting caused by a typical characteristic of empirical BD estimators. In this paper, to mitigate train-loss hacking, we propose a non-negative correction for empirical BD estimators. Theoretically, we confirm the soundness of the proposed method through a generalization error bound. Through our experiments, the proposed methods show a favorable performance in inlier-based outlier detection.


title: Robust Sub-Gaussian Principal Component Analysis and Width-Independent Schatten Packing

abstract: We develop two methods for the following fundamental statistical task: given an $\epsilon$-corrupted set of $n$ samples from a $d$-dimensional sub-Gaussian distribution, return an approximate top eigenvector of the covariance matrix. Our first robust PCA algorithm runs in polynomial time, returns a $1 - O(\epsilon\log\epsilon^{-1})$-approximate top eigenvector, and is based on a simple iterative filtering approach. Our second, which attains a slightly worse approximation factor, runs in nearly-linear time and sample complexity under a mild spectral gap assumption. These are the first polynomial-time algorithms yielding non-trivial information about the covariance of a corrupted sub-Gaussian distribution without requiring additional algebraic structure of moments. As a key technical tool, we develop the first width-independent solvers for Schatten-$p$ norm packing semidefinite programs, giving a $(1 + \epsilon)$-approximate solution in $O(p\log(\tfrac{nd}{\epsilon})\epsilon^{-1})$ input-sparsity time iterations (where $n$, $d$ are problem dimensions).


title: Kernel Distributionally Robust Optimization

abstract: We propose kernel distributionally robust optimization (Kernel DRO) using insights from the robust optimization theory and functional analysis. Our method uses reproducing kernel Hilbert spaces (RKHS) to construct a wide range of convex ambiguity sets, which can be generalized to sets based on integral probability metrics and finite-order moment bounds. This perspective unifies multiple existing robust and stochastic optimization methods. We prove a theorem that generalizes the classical duality in the mathematical problem of moments. Enabled by this theorem, we reformulate the maximization with respect to measures in DRO into the dual program that searches for RKHS functions. Using universal RKHSs, the theorem applies to a broad class of loss functions, lifting common limitations such as polynomial losses and knowledge of the Lipschitz constant. We then establish a connection between DRO and stochastic optimization with expectation constraints. Finally, we propose practical algorithms based on both batch convex solvers and stochastic functional gradient, which apply to general optimization and machine learning tasks.


title: Confidence Interval for Off-Policy Evaluation from Dependent Samples via Bandit Algorithm: Approach from Standardized Martingales

abstract: This study addresses the problem of off-policy evaluation (OPE) from dependent samples obtained via the bandit algorithm. The goal of OPE is to evaluate a new policy using historical data obtained from behavior policies generated by the bandit algorithm. Because the bandit algorithm updates the policy based on past observations, the samples are not independent and identically distributed (i.i.d.). However, several existing methods for OPE do not take this issue into account and are based on the assumption that samples are i.i.d. In this study, we address this problem by constructing an estimator from a standardized martingale difference sequence. To standardize the sequence, we consider using evaluation data or sample splitting with a two-step estimation. This technique produces an estimator with asymptotic normality without restricting a class of behavior policies. In an experiment, the proposed estimator performs better than existing methods, which assume that the behavior policy converges to a time-invariant policy.


title: Characterizing Impacts of Heterogeneity in Federated Learning upon Large-Scale Smartphone Data

abstract: Federated learning (FL) is an emerging, privacy-preserving machine learning paradigm, drawing tremendous attention in both academia and industry. A unique characteristic of FL is heterogeneity, which resides in the various hardware specifications and dynamic states across the participating devices. Theoretically, heterogeneity can exert a huge influence on the FL training process, e.g., causing a device unavailable for training or unable to upload its model updates. Unfortunately, these impacts have never been systematically studied and quantified in existing FL literature. In this paper, we carry out the first empirical study to characterize the impacts of heterogeneity in FL. We collect large-scale data from 136k smartphones that can faithfully reflect heterogeneity in real-world settings. We also build a heterogeneity-aware FL platform that complies with the standard FL protocol but with heterogeneity in consideration. Based on the data and the platform, we conduct extensive experiments to compare the performance of state-of-the-art FL algorithms under heterogeneity-aware and heterogeneity-unaware settings. Results show that heterogeneity causes non-trivial performance degradation in FL, including up to 9.2% accuracy drop, 2.32x lengthened training time, and undermined fairness. Furthermore, we analyze potential impact factors and find that device failure and participant bias are two potential factors for performance degradation. Our study provides insightful implications for FL practitioners. On the one hand, our findings suggest that FL algorithm designers consider necessary heterogeneity during the evaluation. On the other hand, our findings urge system providers to design specific mechanisms to mitigate the impacts of heterogeneity.


title: Robust Design for Intelligent Reflecting Surfaces Assisted MISO Systems

abstract: In this work, we study the statistically robust beamforming design for an intelligent reflecting surfaces (IRS) assisted multiple-input single-output (MISO) wireless system under imperfect channel state information (CSI), where the channel estimation errors are assumed to be additive Gaussian. We aim at jointly optimizing the transmit/receive beamformers and IRS phase shifts to minimize the average mean squared error (MSE) at the user. In particular, to tackle the non-convex optimization problem, an efficient algorithm is developed by capitalizing on alternating optimization and majorization-minimization techniques. Simulation results show that the proposed scheme achieves robust MSE performance in the presence of CSI error, and substantially outperforms conventional non-robust methods.


title: The competition between the intrinsic and Rashba spin-orbit coupling and effects of correlations on Rashba SOC-driven transitions in the Kane-Mele model

abstract: We investigate, firstly, the effects of the Rashba SOC on the band structrue of the Kane-Mele model. The competition between the Rashba SOC and the intrinsic SOC can lead to the rich phenomenology. The Rashba SOC can drive the indirect and direct energy gap to close successively, but maintain the band touching between the valence band and the conduction band when the Rashba SOC is large enough to dominant the competition. We find that these touching points are located at $K$ and $K^{\prime}$ or/and some $2\pi/3$ rotationally symmetric points around $K$ and $K^{\prime}$ in the Brillouin zone. The indirect and direct energy gap closings correspond to the topologically trivial and non-trivial phase transitions respectively. For the small intrinsic SOC, the topologically non-trivial transition occurs when the ratio of the Rashba SOC to the intrinsic SOC is equal to the classical result, i.e. $2\sqrt{3}$. For the large intrinsic SOC, however, we find that the ratio decreases with the increasing intrinsic SOC. Secondly, using the slave-rotor mean field method we investigate the influences of the correlation on the Rashba SOC-driven topologically trivial and non-trivial transition in both the charge condensate and Mott regions. The topological Mott insulator with gapped or gapless spin excitations can arise from the interplay of the Rashba SOC and correlations.


title: Quantum Robust Fitting

abstract: Many computer vision applications need to recover structure from imperfect measurements of the real world. The task is often solved by robustly fitting a geometric model onto noisy and outlier-contaminated data. However, recent theoretical analyses indicate that many commonly used formulations of robust fitting in computer vision are not amenable to tractable solution and approximation. In this paper, we explore the usage of quantum computers for robust fitting. To do so, we examine and establish the practical usefulness of a robust fitting formulation inspired by Fourier analysis of Boolean functions. We then investigate a quantum algorithm to solve the formulation and analyse the computational speed-up possible over the classical algorithm. Our work thus proposes one of the first quantum treatments of robust fitting for computer vision.


title: Improved estimations of stochastic chemical kinetics by finite state expansion

abstract: Stochastic reaction networks are a fundamental model to describe interactions between species where random fluctuations are relevant. The master equation provides the evolution of the probability distribution across the discrete state space consisting of vectors of population counts for each species. However, since its exact solution is often elusive, several analytical approximations have been proposed. The deterministic rate equation (DRE) gives a macroscopic approximation as a compact system of differential equations that estimate the average populations for each species, but it may be inaccurate in the case of nonlinear interaction dynamics. Here we propose finite state expansion (FSE), an analytical method mediating between the microscopic and the macroscopic interpretations of a stochastic reaction network by coupling the master equation dynamics of a chosen subset of the discrete state space with the mean population dynamics of the DRE. An algorithm translates a network into an expanded one where each discrete state is represented as a further distinct species. This translation exactly preserves the stochastic dynamics, but the DRE of the expanded network can be interpreted as a correction to the original one. The effectiveness of FSE is demonstrated in models that challenge state-of-the-art techniques due to intrinsic noise, multi-scale populations, and multi-stability.


title: Backward phase-matched second-harmonic generation from stacked metasurfaces

abstract: We demonstrate phase-matched second-harmonic generation (SHG) from three-dimensional metamaterials consisting of stacked metasurfaces. To achieve phase matching, we utilize a novel mechanism based on phase engineering of the metasurfaces at the interacting wavelengths, facilitating phase-matched SHG in the unconventional backward direction. By stacking up to five metasurfaces, we obtain the expected factor of 25 enhancement in SHG efficiency. Our results motivate further investigations to achieve higher conversion efficiencies also with more complex wavefronts.


title: Universal Properties of a Run-and-Tumble Particle in Arbitrary Dimension

abstract: We consider an active run-and-tumble particle (RTP) in $d$ dimensions, starting from the origin and evolving over a time interval $[0,t]$. We examine three different models for the dynamics of the RTP: the standard RTP model with instantaneous tumblings, a variant with instantaneous runs and a general model in which both the tumblings and the runs are non-instantaneous. For each of these models, we use the Sparre Andersen theorem for discrete-time random walks to compute exactly the probability that the $x$ component does not change sign up to time $t$, showing that it does not depend on $d$. As a consequence of this result, we compute exactly other $x$-component properties, namely the distribution of the time of the maximum and the record statistics, showing that they are universal, i.e. they do not depend on $d$. Moreover, we show that these universal results hold also if the speed $v$ of the particle after each tumbling is random, drawn from a generic probability distribution. Our findings are confirmed by numerical simulations. Some of these results have been announced in a recent Letter [Phys. Rev. Lett. 124, 090603 (2020)].


title: Revisit of Semi-Implicit Schemes for Phase-Field Equations

abstract: It is a very common practice to use semi-implicit schemes in various computations, which treat selected linear terms implicitly and the nonlinear terms explicitly. For phase-field equations, the principal elliptic operator is treated implicitly to reduce the associated stability constraints while the nonlinear terms are still treated explicitly to avoid the expensive process of solving nonlinear equations at each time step. However, very few recent numerical analysis is relevant to semi-implicit schemes, while "stabilized" schemes have become very popular. In this work, we will consider semi-implicit schemes for the Allen-Cahn equation with {\em general potential} function. It will be demonstrated that the maximum principle is valid and the energy stability also holds for the numerical solutions. This paper extends the result of Tang \& Yang (J. Comput. Math., 34(5):471--481, 2016) which studies the semi-implicit scheme for the Allen-Cahn equation with {\em polynomial potentials}.


title: Intelligent Reflecting Surface Operation under Predictable Receiver Mobility: A Continuous Time Propagation Model

abstract: The operation of an intelligent reflecting surface (IRS) under predictable receiver mobility is investigated. We develop a continuous time system model for multipath channels and discuss the optimal IRS configuration with respect to received power, Doppler spread, and delay spread. It is shown that the received power can be maximized without adding Doppler spread to the system. In a numerical case study, we show that an IRS having the size of just two large billboards can improve the link budget of ground to Low Earth Orbit (LEO) satellite links by up to 6 dB. It also adds a second, almost equivalently strong, communication path that improves the link reliability.


title: New dynamical observer for a batch crystallization process based on solute concentration

abstract: In this paper a new observer is introduced to estimate the Crystal Size Distribution (CSD) only from the measurements of the solute concentration, temperature and a model of the growth rate. No model of the nucleation rate is needed. This approach is based on the use of a Kazantzis-Kravaris/Luenberger observer which exponentially estimates functionals of the CSD. Then, the full state is estimated by means of a Tikhonov regularization procedure. Numerical simulations are provided. Our approach relies on an infinite-dimensional observer, contrarily to the usual moment based observers.


title: CANOA: CAN Origin Authentication Through Power Side-Channel Monitoring

abstract: The lack of any sender authentication mechanism in place makes CAN (Controller Area Network) vulnerable to security threats. For instance, an attacker can impersonate an ECU (Electronic Control Unit) on the bus and send spoofed messages unobtrusively with the identifier of the impersonated ECU. To address this problem, we propose a novel sender authentication technique that uses power consumption measurements of the ECU to authenticate the sender of a message. When an ECU is transmitting, its power requirement is affected, and a characteristic pattern appears in its power consumption. Our technique exploits the power consumption of each ECU during the transmission of a message to determine whether the message actually originated from the purported sender. We evaluate our approach in both a lab setup and a real vehicle. We also evaluate our approach against factors that can impact the power consumption measurement of the ECU. The results of the evaluation show that the proposed technique is applicable in a broad range of operating conditions with reasonable computational power requirements and attaining good accuracy.


title: Sparse approximation of triangular transports. Part I: the finite dimensional case

abstract: For two probability measures $\rho$ and $\pi$ with analytic densities on the $d$-dimensional cube $[-1,1]^d$, we investigate the approximation of the unique triangular monotone Knothe-Rosenblatt transport $T:[-1,1]^d\to [-1,1]^d$, such that the pushforward $T_\sharp\rho$ equals $\pi$. It is shown that for $d\in\mathbb{N}$ there exist approximations $\tilde T$ of $T$, based on either sparse polynomial expansions or deep ReLU neural networks, such that the distance between $\tilde T_\sharp\rho$ and $\pi$ decreases exponentially. More precisely, we prove error bounds of the type $\exp(-\beta N^{1/d})$ (or $\exp(-\beta N^{1/(d+1)})$ for neural networks), where $N$ refers to the dimension of the ansatz space (or the size of the network) containing $\tilde T$; the notion of distance comprises the Hellinger distance, the total variation distance, the Wasserstein distance and the Kullback-Leibler divergence. Our construction guarantees $\tilde T$ to be a monotone triangular bijective transport on the hypercube $[-1,1]^d$. Analogous results hold for the inverse transport $S=T^{-1}$. The proofs are constructive, and we give an explicit a priori description of the ansatz space, which can be used for numerical implementations.


title: Projecting onto intersections of halfspaces and hyperplanes

abstract: It is well-known that the sequence of iterations of the composition of projections onto closed affine subspaces converges linearly to the projection onto the intersection of the affine subspaces when the sum of the corresponding linear subspaces is closed. Inspired by this, in this work, we systematically study the relation between the projection onto intersection of halfspaces and hyperplanes, and the composition of projections onto halfspaces and hyperplanes. In addition, as by-products, we provide the Karush-Kuhn-Tucker conditions for characterizing the optimal solution of convex optimization with finitely many equality and inequality constraints in Hilbert spaces and construct an explicit formula for the projection onto the intersection of hyperplane and halfspace.


title: Kidney segmentation in neck-to-knee body MRI of 40,000 UK Biobank participants

abstract: The UK Biobank is collecting extensive data on health-related characteristics of over half a million volunteers. The biological samples of blood and urine can provide valuable insight on kidney function, with important links to cardiovascular and metabolic health. Further information on kidney anatomy could be obtained by medical imaging. In contrast to the brain, heart, liver, and pancreas, no dedicated Magnetic Resonance Imaging (MRI) is planned for the kidneys. An image-based assessment is nonetheless feasible in the neck-to-knee body MRI intended for abdominal body composition analysis, which also covers the kidneys. In this work, a pipeline for automated segmentation of parenchymal kidney volume in UK Biobank neck-to-knee body MRI is proposed. The underlying neural network reaches a relative error of 3.8%, with Dice score 0.956 in validation on 64 subjects, close to the 2.6% and Dice score 0.962 for repeated segmentation by one human operator. The released MRI of about 40,000 subjects can be processed within two days, yielding volume measurements of left and right kidney. Algorithmic quality ratings enabled the exclusion of outliers and potential failure cases. The resulting measurements can be studied and shared for large-scale investigation of associations and longitudinal changes in parenchymal kidney volume.


title: Complex Dynamics in Simple Neural Networks: Understanding Gradient Flow in Phase Retrieval

abstract: Despite the widespread use of gradient-based algorithms for optimizing high-dimensional non-convex functions, understanding their ability of finding good minima instead of being trapped in spurious ones remains to a large extent an open problem. Here we focus on gradient flow dynamics for phase retrieval from random measurements. When the ratio of the number of measurements over the input dimension is small the dynamics remains trapped in spurious minima with large basins of attraction. We find analytically that above a critical ratio those critical points become unstable developing a negative direction toward the signal. By numerical experiments we show that in this regime the gradient flow algorithm is not trapped; it drifts away from the spurious critical points along the unstable direction and succeeds in finding the global minimum. Using tools from statistical physics we characterize this phenomenon, which is related to a BBP-type transition in the Hessian of the spurious minima.

categories: cs.LG cond-mat.dis-nn math.ST stat.ML stat.TH


title: Random forest estimation of conditional distribution functions and conditional quantiles

abstract: We propose a theoretical study of two realistic estimators of conditional distribution functions and conditional quantiles using random forests. The estimation process uses the bootstrap samples generated from the original dataset when constructing the forest. Bootstrap samples are reused to define the first estimator, while the second requires only the original sample, once the forest has been built. We prove that both proposed estimators of the conditional distribution functions are consistent uniformly a.s. To the best of our knowledge, it is the first proof of consistency including the bootstrap part. We also illustrate the estimation procedures on a numerical example.


title: Four quark operators for kaon bag parameter with gradient flow

abstract: To study the CP-violation using the $K_0-\bar{K}_0$ oscillation, we need the kaon bag parameter which represents QCD corrections in the leading Feynman diagrams. The lattice QCD provides us with the only way to evaluate the kaon bag parameter directly from the first principles of QCD. However, a calculation of relevant four quark operators with theoretically sound Wilson-type lattice quarks had to carry a numerically big burden of extra renormalizations and resolution of extra mixings due to the explicit chiral violation. Recently, the Small Flow-time eXpansion (SFtX) method was proposed as a general method based on the gradient flow to correctly calculate any renormalized observables on the lattice, irrespective of the explicit violations of related symmetries on the lattice. To apply the SFtX method, we need matching coefficients, which relate finite operators at small flow-times in the gradient flow scheme to renormalized observables in conventional renormalization schemes. In this paper, we calculate the matching coefficients for four quark operators and quark bi-linear operators, relevant to the kaon bag parameter.


title: Doubly random polytopes

abstract: A two-step model for generating random polytopes is considered. For parameters $d$, $m$, and $p$, the first step is to generate a simple polytope $P$ whose facets are given by $m$ uniform random hyperplanes tangent to the unit sphere in $\mathbb{R}^d$, and the second step is to sample each vertex of $P$ independently with probability $p$ and let $Q$ be the convex hull of the sampled vertices. We establish results on how well $Q$ approximates the unit sphere in terms of $m$ and $p$ as well as asymptotics on the combinatorial complexity of $Q$ for certain regimes of $p$.


title: Markov Random Geometric Graph (MRGG): A Growth Model for Temporal Dynamic Networks

abstract: We introduce Markov Random Geometric Graphs (MRGGs), a growth model for temporal dynamic networks. It is based on a Markovian latent space dynamic: consecutive latent points are sampled on the Euclidean Sphere using an unknown Markov kernel; and two nodes are connected with a probability depending on a unknown function of their latent geodesic distance. More precisely, at each stamp-time $k$ we add a latent point $X_k$ sampled by jumping from the previous one $X_{k-1}$ in a direction chosen uniformly $Y_k$ and with a length $r_k$ drawn from an unknown distribution called the latitude function. The connection probabilities between each pair of nodes are equal to the envelope function of the distance between these two latent points. We provide theoretical guarantees for the non-parametric estimation of the latitude and the envelope functions.We propose an efficient algorithm that achieves those non-parametric estimation tasks based on an ad-hoc Hierarchical Agglomerative Clustering approach. As a by product, we show how MRGGs can be used to detect dependence structure in growing graphs and to solve link prediction problems.