title: Learning to solve TV regularized problems with unrolled algorithms

abstract: Total Variation (TV) is a popular regularization strategy that promotes piece-wise constant signals by constraining the $\ell_1$-norm of the first order derivative of the estimated signal. The resulting optimization problem is usually solved using iterative algorithms such as proximal gradient descent, primal-dual algorithms or ADMM. However, such methods can require a very large number of iterations to converge to a suitable solution. In this paper, we accelerate such iterative algorithms by unfolding proximal gradient descent solvers in order to learn their parameters for 1D TV regularized problems. While this could be done using the synthesis formulation, we demonstrate that this leads to slower performances. The main difficulty in applying such methods in the analysis formulation lies in proposing a way to compute the derivatives through the proximal operator. As our main contribution, we develop and characterize two approaches to do so, describe their benefits and limitations, and discuss the regime where they can actually improve over iterative procedures. We validate those findings with experiments on synthetic and real data.


title: Model-based Policy Optimization with Unsupervised Model Adaptation

abstract: Model-based reinforcement learning methods learn a dynamics model with real data sampled from the environment and leverage it to generate simulated data to derive an agent. However, due to the potential distribution mismatch between simulated data and real data, this could lead to degraded performance. Despite much effort being devoted to reducing this distribution mismatch, existing methods fail to solve it explicitly. In this paper, we investigate how to bridge the gap between real and simulated data due to inaccurate model estimation for better policy optimization. To begin with, we first derive a lower bound of the expected return, which naturally inspires a bound maximization algorithm by aligning the simulated and real data distributions. To this end, we propose a novel model-based reinforcement learning framework AMPO, which introduces unsupervised model adaptation to minimize the integral probability metric (IPM) between feature distributions from real and simulated data. Instantiating our framework with Wasserstein-1 distance gives a practical model-based approach. Empirically, our approach achieves state-of-the-art performance in terms of sample efficiency on a range of continuous control benchmark tasks.


title: Deep-learning interatomic potential for irradiation damage simulations in MoS2 with ab initial accuracy

abstract: Potentials that could accurately describe the irradiation damage processes are highly desired to figure out the atomic-level response of various newly-discovered materials under irradiation environments. In this work, we introduce a deep-learning interatomic potential for monolayer MoS2 by combining all-electron calculations, an active-learning sampling method and a hybrid deep-learning model. This potential could not only give an overall good performance on the predictions of near-equilibrium material properties including lattice constants, elastic coefficients, energy stress curves, phonon spectra, defect formation energy and displacement threshold, but also reproduce the ab initial irradiation damage processes with high quality. Further irradiation simulations indicate that one single highenergy ion could generate a large nanopore with a diameter of more than 2 nm, or a series of multiple nanopores, which is qualitatively verified by the subsequent 500 keV Au+ ion irradiation experiments. This work provides a promising and feasible approach to simulate irradiation effects in enormous newly-discovered materials with unprecedented accuracy.


title: RONELD: Robust Neural Network Output Enhancement for Active Lane Detection

abstract: Accurate lane detection is critical for navigation in autonomous vehicles, particularly the active lane which demarcates the single road space that the vehicle is currently traveling on. Recent state-of-the-art lane detection algorithms utilize convolutional neural networks (CNNs) to train deep learning models on popular benchmarks such as TuSimple and CULane. While each of these models works particularly well on train and test inputs obtained from the same dataset, the performance drops significantly on unseen datasets of different environments. In this paper, we present a real-time robust neural network output enhancement for active lane detection (RONELD) method to identify, track, and optimize active lanes from deep learning probability map outputs. We first adaptively extract lane points from the probability map outputs, followed by detecting curved and straight lanes before using weighted least squares linear regression on straight lanes to fix broken lane edges resulting from fragmentation of edge maps in real images. Lastly, we hypothesize true active lanes through tracking preceding frames. Experimental results demonstrate an up to two-fold increase in accuracy using RONELD on cross-dataset validation tests.


title: Use of Uncertain Additional Information in Newsvendor Models

abstract: The newsvendor problem is a popular inventory management problem in supply chain management and logistics. Solutions to the newsvendor problem determine optimal inventory levels. This model is typically fully determined by a purchase and sale prices and a distribution of random market demand. From a statistical point of view, this problem is often considered as a quantile estimation of a critical fractile which maximizes anticipated profit. The distribution of demand is a random variable and is often estimated on historic data. In an ideal situation, when the probability distribution of the demand is known, one can determine the quantile of a critical fractile minimizing a particular loss function. Since maximum likelihood estimation is asymptotically efficient, under certain regularity assumptions, the maximum likelihood estimators are used for the quantile estimation problem. Then, the Cramer-Rao lower bound determines the lowest possible asymptotic variance. Can one find a quantile estimate with a smaller variance then the Cramer-Rao lower bound? If a relevant additional information is available then the answer is yes. Additional information may be available in different forms. This manuscript considers minimum variance and minimum mean squared error estimation for incorporating additional information for estimating optimal inventory levels. By a more precise assessment of optimal inventory levels, we maximize expected profit


title: Freeform imaging systems: Fermat's principle unlocks 'first time right' design

abstract: For more than 150 years, scientists have advanced aberration theory to describe, analyse, and eliminate imperfections that disturb the imaging quality of optical components and systems. Simultaneously, they have developed optical design methods for and manufacturing techniques of imaging systems with ever-increasing complexity and performance up to the point where they are now including optical elements that are unrestricted in their surface shape. These so-called optical freeform elements offer degrees of freedom that can greatly extend the functionalities and further boost the specifications of state-of-the-art imaging systems. However, the drastically increased number of surface coefficients of these freeform surfaces poses severe challenges for the optical design process, such that the deployment of freeform optics remained limited until today. In this paper, we present a deterministic direct optical design method for freeform imaging systems based on differential equations derived from Fermat's principle and solved using power series. The method allows calculating the optical surface coefficients that ensure minimal image blurring for each individual order of aberrations. We demonstrate the systematic, deterministic, scalable and holistic character of our method with catoptric and catadioptric design examples. As such we introduce a disruptive methodology to design optical imaging systems 'first time right', we eliminate the 'trial and error' approach in present-day optical design, and we pave the way to a fast-track uptake of freeform elements to create the next-generation high-end optics.


title: A topological approach to undefinability in algebraic extensions of $\mathbb{Q}$

abstract: For any subset $Z \subseteq \mathbb{Q}$, consider the set $S_Z$ of subfields $L\subseteq \overline{\mathbb{Q}}$ which contain a co-infinite subset $C \subseteq L$ that is universally definable in $L$ such that $C \cap \mathbb{Q}=Z$. Placing a natural topology on the set $\text{Sub}(\overline{\mathbb{Q}})$ of subfields of $\overline{\mathbb{Q}}$, we show that if $Z$ is not thin in $\mathbb{Q}$, then $S_Z$ is meager in $\text{Sub}(\overline{\mathbb{Q}})$. Here, thin and meager both mean "small", in terms of arithmetic geometry and topology, respectively. For example, this implies that only a meager set of fields $L$ have the property that the ring of algebraic integers $\mathcal{O}_L$ is universally definable in $L$. The main tools are Hilbert's Irreducibility Theorem and a new normal form theorem for existential definitions. The normal form theorem, which may be of independent interest, says roughly that every $\exists$-definable subset of an algebraic extension of $\mathbb Q$ is a finite union of single points and projections of hypersurfaces defined by absolutely irreducible polynomials.


title: Optimizing Radiotherapy Plans for Cancer Treatment with Tensor Networks

abstract: We present a novel application of Tensor Network methods in cancer treatment as a potential tool to solve the dose optimization problem in radiotherapy. In particular, the Intensity-Modulated Radiation Therapy (IMRT) technique - that allows treating irregular and inhomogeneous tumors while reducing the radiation toxicity on healthy organs - is based on the optimization of the radiation beamlets intensities. The optimization aims to maximize the delivery of the therapy dose to cancer while avoiding the organs at risk to prevent their damage by the radiation. Here, we map the dose optimization problem into the search of the ground state of an Ising-like Hamiltonian, describing a system of long-range interacting qubits. Finally, we apply a Tree Tensor Network algorithm to find the ground-state of the Hamiltonian. In particular, we present an anatomical scenario exemplifying a prostate cancer treatment. A similar approach can be applied to future hybrid classical-quantum algorithms, paving the way for the use of quantum technologies in future medical treatments.

categories: physics.med-ph cond-mat.stat-mech physics.comp-ph quant-ph


title: Survey on Causal-based Machine Learning Fairness Notions

abstract: Addressing the problem of fairness is crucial to safely use machine learning algorithms to support decisions with a critical impact on people's lives such as job hiring, child maltreatment, disease diagnosis, loan granting, etc. Several notions of fairness have been defined and examined in the past decade, such as statistical parity and equalized odds. The most recent fairness notions, however, are causal-based and reflect the now widely accepted idea that using causality is necessary to appropriately address the problem of fairness. This paper examines an exhaustive list of causal-based fairness notions and study their applicability in real-world scenarios. As the majority of causal-based fairness notions are defined in terms of non-observable quantities (e.g., interventions and counterfactuals), their deployment in practice requires to compute or estimate those quantities using observational data. This paper offers a comprehensive report of the different approaches to infer causal quantities from observational data including identifiability (Pearl's SCM framework) and estimation (potential outcome framework). The main contributions of this survey paper are (1) a guideline to help selecting a suitable fairness notion given a specific real-world scenario, and (2) a ranking of the fairness notions according to Pearl's causation ladder indicating how difficult it is to deploy each notion in practice.


title: Strain-mediated high conductivity in ultrathin antiferromagnetic metallic nitrides

abstract: Strain engineering provides the ability to control the ground states and associated phase transition in the epitaxial films. However, the systematic study of intrinsic characters and their strain dependency in transition-metal nitrides remains challenging due to the difficulty in fabricating the stoichiometric and high-quality films. Here we report the observation of electronic state transition in highly crystalline antiferromagnetic CrN films with strain and reduced dimensionality. Shrinking the film thickness to a critical value of ~ 30 unit cells, a profound conductivity reduction accompanied by unexpected volume expansion is observed in CrN films. The electrical conductivity is observed surprisingly when the CrN layer as thin as single unit cell thick, which is far below the critical thickness of most metallic films. We found that the metallicity of an ultrathin CrN film recovers from an insulating behavior upon the removal of as-grown strain by fabrication of first-ever freestanding nitride films. Both first-principles calculations and linear dichroism measurements reveal that the strain-mediated orbital splitting effectively customizes the relatively small bandgap at the Fermi level, leading to exotic phase transition in CrN. The ability to achieve highly conductive nitride ultrathin films by harness strain-controlling over competing phases can be used for utilizing their exceptional characteristics.


title: Learning a Low-dimensional Representation of a Safe Region for Safe Reinforcement Learning on Dynamical Systems

abstract: For safely applying reinforcement learning algorithms on high-dimensional nonlinear dynamical systems, a simplified system model is used to formulate a safe reinforcement learning framework. Based on the simplified system model, a low-dimensional representation of the safe region is identified and is used to provide safety estimates for learning algorithms. However, finding a satisfying simplified system model for complex dynamical systems usually requires a considerable amount of effort. To overcome this limitation, we propose in this work a general data-driven approach that is able to efficiently learn a low-dimensional representation of the safe region. Through an online adaptation method, the low-dimensional representation is updated by using the feedback data such that more accurate safety estimates are obtained. The performance of the proposed approach for identifying the low-dimensional representation of the safe region is demonstrated with a quadcopter example. The results show that, compared to previous work, a more reliable and representative low-dimensional representation of the safe region is derived, which then extends the applicability of the safe reinforcement learning framework.


title: Thermal statistics of small magnets

abstract: While the canonical ensemble has been tremendously successful in capturing thermal statistics of macroscopic systems, deviations from canonical behavior exhibited by small systems are not well understood. Here, using a small two dimensional Ising magnet embedded inside a larger Ising magnet heat bath, we characterize the failures of the canonical ensemble when describing small systems. We find significant deviations from the canonical behavior for small systems near and below the critical point of the two dimensional Ising model. Notably, the agreement with the canonical ensemble is driven not by the system size but by the statistical decoupling between the system and its surrounding. A superstatistical framework wherein we allow the temperature of the small magnet to vary is able to capture its thermal statistics with significantly higher accuracy than the Gibbs-Boltzmann distribution. We discuss future directions.


title: A Versatile Crack Inspection Portable System based on Classifier Ensemble and Controlled Illumination

abstract: This paper presents a novel setup for automatic visual inspection of cracks in ceramic tile as well as studies the effect of various classifiers and height-varying illumination conditions for this task. The intuition behind this setup is that cracks can be better visualized under specific lighting conditions than others. Our setup, which is designed for field work with constraints in its maximum dimensions, can acquire images for crack detection with multiple lighting conditions using the illumination sources placed at multiple heights. Crack detection is then performed by classifying patches extracted from the acquired images in a sliding window fashion. We study the effect of lights placed at various heights by training classifiers both on customized as well as state-of-the-art architectures and evaluate their performance both at patch-level and image-level, demonstrating the effectiveness of our setup. More importantly, ours is the first study that demonstrates how height-varying illumination conditions can affect crack detection with the use of existing state-of-the-art classifiers. We provide an insight about the illumination conditions that can help in improving crack detection in a challenging real-world industrial environment.


title: Four-component Full Configuration Interaction Quantum Monte Carlo for Relativistic Correlated Electron Problems

abstract: An adaptation of the full configuration interaction quantum Monte Carlo (FCIQMC) method is presented, for correlated electron problems containing heavy elements and the presence of significant relativistic effects. The modified algorithm allows for the sampling of the four-component spinors of the Dirac--Coulomb(--Breit) Hamiltonian within the relativistic no-pair approximation. The loss of spin symmetry and the general requirement for complex-valued Hamiltonian matrix elements are the most immediate considerations in expanding the scope of FCIQMC into the relativistic domain, and the alternatives for their efficient implementation are motivated and demonstrated. For the canonical correlated four-component chemical benchmark application of Thallium Hydride, we show that the necessary modifications do not particularly adversely affect the convergence of the systematic (initiator) error to the exact correlation energy for FCIQMC calculations, which is primarily dictated by the sparsity of the wave function, allowing the computational effort to somewhat bypass the formal increases in Hilbert space dimension for these problems. We apply the method to the larger problem of the spectroscopic constants of Tin Oxide, correlating 28 electrons in 122 Kramers-paired spinors, finding good agreement with experimental and prior theoretical relativistic studies.


title: Multilayer Network Analysis for Improved Credit Risk Prediction

abstract: We present a multilayer network model for credit risk assessment. Our model accounts for multiple connections between borrowers (such as their geographic location and their economic activity) and allows for explicitly modelling the interaction between connected borrowers. We develop a multilayer personalized PageRank algorithm that allows quantifying the strength of the default exposure of any borrower in the network. We test our methodology in an agricultural lending framework, where it has been suspected for a long time default correlates between borrowers when they are subject to the same structural risks. Our results show there are significant predictive gains just by including centrality multilayer network information in the model, and these gains are increased by more complex information such as the multilayer PageRank variables. The results suggest default risk is highest when an individual is connected to many defaulters, but this risk is mitigated by the size of the neighbourhood of the individual, showing both default risk and financial stability propagate throughout the network.


title: A sequence of algebraic integer relation numbers which converges to 4

abstract: Let $\alpha \in \mathbb{R}$ and let $$A=\begin{bmatrix} 1 & 1 \\ 0 & 1\end{bmatrix} \ \text{and} \ B_{\alpha} = \begin{bmatrix} 1 & 0 \\ \alpha & 1\end{bmatrix}.$$ The subgroup $G_\alpha$ of $\mathrm{SL}_2(\mathbb{R})$ is a group generated by the matrices $A$ and $B_\alpha$. In this paper, we investigate the property of the group $G_\alpha.$ We construct a generalization of the Farey graph for the subgroup $G_\alpha.$ This graph determines whether the group $G_\alpha$ is a free group of rank $2$. More precisely, the group $G_\alpha$ is a free group of rank $2$ if and only if the graph is tree. In particular, we show that if $1/2$ is a vertex of the graph, then $G_\alpha$ is not a free group of rank $2$. Using this, we construct a sequence of real numbers so that the sequence converges to $4$ and each number has the corresponding group that is not a free group of rank $2$. It turns out that the real numbers are algebraic integers.


title: Domain Generalized Person Re-Identification via Cross-Domain Episodic Learning

abstract: Aiming at recognizing images of the same person across distinct camera views, person re-identification (re-ID) has been among active research topics in computer vision. Most existing re-ID works require collection of a large amount of labeled image data from the scenes of interest. When the data to be recognized are different from the source-domain training ones, a number of domain adaptation approaches have been proposed. Nevertheless, one still needs to collect labeled or unlabelled target-domain data during training. In this paper, we tackle an even more challenging and practical setting, domain generalized (DG) person re-ID. That is, while a number of labeled source-domain datasets are available, we do not have access to any target-domain training data. In order to learn domain-invariant features without knowing the target domain of interest, we present an episodic learning scheme which advances meta learning strategies to exploit the observed source-domain labeled data. The learned features would exhibit sufficient domain-invariant properties while not overfitting the source-domain data or ID labels. Our experiments on four benchmark datasets confirm the superiority of our method over the state-of-the-arts.


title: Adelic superrigidity and profinitely solitary lattices

abstract: By arithmeticity and superrigidity, a commensurability class of lattices in a higher rank Lie group is defined by a unique algebraic group over a unique number subfield of $\mathbb{R}$ or $\mathbb{C}$. We prove an adelic version of superrigidity which implies that two such commensurability classes define the same profinite commensurability class if and only if the algebraic groups are adelically isomorphic. We discuss noteworthy consequences on profinite rigidity questions.


title: A tutorial comparing different covariate balancing methods with an application evaluating the causal effects of substance use treatment programs for adolescents

abstract: Randomized controlled trials are the gold standard for measuring causal effects. However, they are often not always feasible, and causal treatment effects must be estimated from observational data. Observational studies do not allow robust conclusions about causal relationships unless statistical techniques account for the imbalance of pretreatment confounders across groups while key assumptions hold. Propensity score and balance weighting (PSBW) are useful techniques that aim to reduce the imbalances between treatment groups by weighting the groups to look alike on the observed confounders. There are many methods available to estimate PSBW. However, it is unclear a priori which will achieve the best trade-off between covariate balance and effective sample size. Moreover, it is critical to assess the validity of key assumptions required for robust estimation of the needed treatment effects, including the overlap and no unmeasured confounding assumptions. We present a step-by-step guide to covariate balancing strategies, including how to evaluate overlap, obtain estimates of PSBW, check for covariate balance, and assess sensitivity to unobserved confounding. We compare the performance of several estimation methods using a case study examining the relative effectiveness of substance use treatment programs and provide a user-friendly web application that can implement the proposed steps.


title: M Theory from F Theory

abstract: We write down a $GL(D+1)$ (D for the dimension of string theory) manifest fundamental brane worldvolume current algebra description of M theory, which consists of a pair of vector field $X^m$ and dual 2-form field $X_{mn}$, compositing together to parametrize the spacetime, with a selfduality condition for sectioning. The worldvolume of the brane itself is a (D+2) dimensional object, and the background spacetime after sectioning has dimension (D+1). We summarize the features of the algebra. The field contents of the corresponding background geometry, the usual vielbein $e_a{}^m$ and the 3-form $A_{mnp}$, could be identified as different blocks of the composite spacetime vielbein, by solving the orthogonality condition. Their behaviour under gauge transformation are also determined by the corresponding rules of the composite vielbein. Then by solving F theory $\mathcal{V}$ constraints, we reduce the number of worldvolume and show how to recover M theory from F theory.


title: On bodies floating in equilibrium in every orientation

abstract: Ulam's problem 19 from the Scottish Book asks: {\it is a solid of uniform density which floats in water in every position necessarily a sphere?} We obtain several results related to this problem.


title: On the Generalisation Capabilities of Fingerprint Presentation Attack Detection Methods in the Short Wave Infrared Domain

abstract: Nowadays, fingerprint-based biometric recognition systems are becoming increasingly popular. However, in spite of their numerous advantages, biometric capture devices are usually exposed to the public and thus vulnerable to presentation attacks (PAs). Therefore, presentation attack detection (PAD) methods are of utmost importance in order to distinguish between bona fide and attack presentations. Due to the nearly unlimited possibilities to create new presentation attack instruments (PAIs), unknown attacks are a threat to existing PAD algorithms. This fact motivates research on generalisation capabilities in order to find PAD methods that are resilient to new attacks. In this context, we evaluate the generalisability of multiple PAD algorithms on a dataset of 19,711 bona fide and 4,339 PA samples, including 45 different PAI species. The PAD data is captured in the short wave infrared domain and the results discuss the advantages and drawbacks of this PAD technique regarding unknown attacks.


title: Solving relaxations of MAP-MRF problems: Combinatorial in-face Frank-Wolfe directions

abstract: We consider the problem of solving LP relaxations of MAP-MRF inference problems, and in particular the method proposed recently in (Swoboda, Kolmogorov 2019; Kolmogorov, Pock 2021). As a key computational subroutine, it uses a variant of the Frank-Wolfe (FW) method to minimize a smooth convex function over a combinatorial polytope. We propose an efficient implementation of this subproutine based on in-face Frank-Wolfe directions, introduced in (Freund et al. 2017) in a different context. More generally, we define an abstract data structure for a combinatorial subproblem that enables in-face FW directions, and describe its specialization for tree-structured MAP-MRF inference subproblems. Experimental results indicate that the resulting method is the current state-of-art LP solver for some classes of problems. Our code is available at https://pub.ist.ac.at/~vnk/papers/IN-FACE-FW.html.


title: Introducing and Applying Newtonian Blurring: An Augmented Dataset of 126,000 Human Connectomes at braingraph.org

abstract: Gaussian blurring is a well-established method for image data augmentation: it may generate a large set of images from a small set of pictures for training and testing purposes for Artificial Intelligence (AI) applications. When we apply AI for non-imagelike biological data, hardly any related method exists. Here we introduce the "Newtonian blurring" in human braingraph (or connectome) augmentation: Started from a dataset of 1053 subjects, we first repeat a probabilistic weighted braingraph construction algorithm 10 times for describing the connections of distinct cerebral areas, then take 7 repetitions in every possible way, delete the lower and upper extremes, and average the remaining 7-2=5 edge-weights for the data of each subject. This way we augment the 1053 graph-set to 120 x 1053 = 126,360 graphs. In augmentation techniques, it is an important requirement that no artificial additions should be introduced into the dataset. Gaussian blurring and also this Newtonian blurring satisfy this goal. The resulting dataset of 126,360 graphs, each in 5 resolutions (i.e., 631,800 graphs in total), is freely available at the site https://braingraph.org/cms/download-pit-group-connectomes/. Augmenting with Newtonian blurring may also be applicable in other non-image related fields, where probabilistic processing and data averaging are implemented.


title: Against All Odds: Winning the Defense Challenge in an Evasion Competition with Diversification

abstract: Machine learning-based systems for malware detection operate in a hostile environment. Consequently, adversaries will also target the learning system and use evasion attacks to bypass the detection of malware. In this paper, we outline our learning-based system PEberus that got the first place in the defender challenge of the Microsoft Evasion Competition, resisting a variety of attacks from independent attackers. Our system combines multiple, diverse defenses: we address the semantic gap, use various classification models, and apply a stateful defense. This competition gives us the unique opportunity to examine evasion attacks under a realistic scenario. It also highlights that existing machine learning methods can be hardened against attacks by thoroughly analyzing the attack surface and implementing concepts from adversarial learning. Our defense can serve as an additional baseline in the future to strengthen the research on secure learning.


title: Bayesian Neural Networks with Soft Evidence

abstract: Bayes's rule deals with hard evidence, that is, we can calculate the probability of event $A$ occuring given that event $B$ has occurred. Soft evidence, on the other hand, involves a degree of uncertainty about whether event $B$ has actually occurred or not. Jeffrey's rule of conditioning provides a way to update beliefs in the case of soft evidence. We provide a framework to learn a probability distribution on the weights of a neural network trained using soft evidence by way of two simple algorithms for approximating Jeffrey conditionalization. We propose an experimental protocol for benchmarking these algorithms on empirical datasets and find that Jeffrey based methods are competitive or better in terms of accuracy yet show improvements in calibration metrics upwards of 20% in some cases, even when the data contains mislabeled points.


title: Theory of the low-temperature longitudinal spin Seebeck effect

abstract: Using a simplified microscopic model of coupled spin and lattice excitations in a ferromagnetic insulator we evaluate the magnetic-field dependence of the spin Seebeck effect at low temperatures. The model includes Heisenberg exchange coupling, a harmonic lattice potential, and a pseudo-dipolar exchange interaction. Our approach goes beyond previous work [Phys. Rev. B 98, 134421 (2018)] in that it does not rely on the a priori assumption of a fast equilibration of the magnon and phonon distributions. Our theory shows that singular features in the magnetic-field dependence of the spin Seebeck effect at low temperatures observed by Kikkawa et al. [Phys. Rev. Lett. 117, 207203 (2016)] are independent of the relative strength of magnon-impurity and phonon-impurity scattering.


title: Teacher-Student Competition for Unsupervised Domain Adaptation

abstract: With the supervision from source domain only in class-level, existing unsupervised domain adaptation (UDA) methods mainly learn the domain-invariant representations from a shared feature extractor, which causes the source-bias problem. This paper proposes an unsupervised domain adaptation approach with Teacher-Student Competition (TSC). In particular, a student network is introduced to learn the target-specific feature space, and we design a novel competition mechanism to select more credible pseudo-labels for the training of student network. We introduce a teacher network with the structure of existing conventional UDA method, and both teacher and student networks compete to provide target pseudo-labels to constrain every target sample's training in student network. Extensive experiments demonstrate that our proposed TSC framework significantly outperforms the state-of-the-art domain adaptation methods on Office-31 and ImageCLEF-DA benchmarks.


title: Rayleigh scattering of a linearly polarized light: scenario of the complete experiment

abstract: The process of the elastic scattering of photons on atoms, known as the Rayleigh scattering, is investigated. Expressing the scattering observables in terms of the electric and magnetic complex scattering amplitudes, we work over the scenarios for their independent benchmarking by experiments. In particular, the combination of the cross section and polarization transfer of initially linearly polarized light provides an opportunity for a complete experiment for any but fixed scattering angle. It allows us to deduce the modulus of the electric and magnetic amplitudes together with the phase difference individually. The findings are illustrated by the theoretical calculations of the scattering of 150 - 250 keV photons on lead atoms.


title: Machine Learning Evaluation of the Echo-Chamber Effect in Medical Forums

abstract: We propose the Echo-Chamber Effect assessment of an online forum. Sentiments perceived by the forum readers are at the core of the analysis; a complete message is the unit of the study. We build 14 models and apply those to represent discussions gathered from an online medical forum. We use four multi-class sentiment classification applications and two Machine Learning algorithms to evaluate prowess of the assessment models.


title: Burton-Cabrera-Frank theory for surfaces with alternating step types

abstract: Burton-Cabrera-Frank (BCF) theory has proven to be a versatile framework to relate surface morphology and dynamics during crystal growth to the underlying mechanisms of adatom diffusion and attachment at steps. For an important class of crystal surfaces, including the basal planes of hexagonal close-packed and related systems, the steps in a sequence on a vicinal surface can exhibit properties that alternate from step to step. Here we develop BCF theory for such surfaces, relating observables such as alternating terrace widths as a function of growth conditions to the kinetic coefficients for adatom attachment at steps. We include the effects of step transparency and step-step repulsion. A general solution is obtained for the dynamics of the terrace widths, assuming quasi-steady-state adatom distributions on the terraces. An explicit simplified analytical solution is obtained under widely applicable approximations. From this we obtain expressions for the full-steady-state terrace fraction as a function of growth rate. Fits of the theoretical predictions to recent experimental determinations of the steady-state and dynamics of terrace fractions on GaN (0001) surfaces during organo-metallic vapor phase epitaxy give values of the kinetic coefficients for this system. In Appendices, we also connect a model for diffusion between kinks on steps to the model for diffusion between steps on terraces, which quantitatively relates step transparency to the kinetics of atom attachment at kinks, and consider limiting cases of diffusion-limited, attachment-limited, and mixed kinetics.


title: Robust High Dimensional Expectation Maximization Algorithm via Trimmed Hard Thresholding

abstract: In this paper, we study the problem of estimating latent variable models with arbitrarily corrupted samples in high dimensional space ({\em i.e.,} $d\gg n$) where the underlying parameter is assumed to be sparse. Specifically, we propose a method called Trimmed (Gradient) Expectation Maximization which adds a trimming gradients step and a hard thresholding step to the Expectation step (E-step) and the Maximization step (M-step), respectively. We show that under some mild assumptions and with an appropriate initialization, the algorithm is corruption-proofing and converges to the (near) optimal statistical rate geometrically when the fraction of the corrupted samples $\epsilon$ is bounded by $ \tilde{O}(\frac{1}{\sqrt{n}})$. Moreover, we apply our general framework to three canonical models: mixture of Gaussians, mixture of regressions and linear regression with missing covariates. Our theory is supported by thorough numerical results.


title: GANs for learning from very high class conditional noisy labels

abstract: We use Generative Adversarial Networks (GANs) to design a class conditional label noise (CCN) robust scheme for binary classification. It first generates a set of correctly labelled data points from noisy labelled data and 0.1% or 1% clean labels such that the generated and true (clean) labelled data distributions are close; generated labelled data is used to learn a good classifier. The mode collapse problem while generating correct feature-label pairs and the problem of skewed feature-label dimension ratio ($\sim$ 784:1) are avoided by using Wasserstein GAN and a simple data representation change. Another WGAN with information-theoretic flavour on top of the new representation is also proposed. The major advantage of both schemes is their significant improvement over the existing ones in presence of very high CCN rates, without either estimating or cross-validating over the noise rates. We proved that KL divergence between clean and noisy distribution increases w.r.t. noise rates in symmetric label noise model; can be extended to high CCN rates. This implies that our schemes perform well due to the adversarial nature of GANs. Further, use of generative approach (learning clean joint distribution) while handling noise enables our schemes to perform better than discriminative approaches like GLC, LDMI and GCE; even when the classes are highly imbalanced. Using Friedman F test and Nemenyi posthoc test, we showed that on high dimensional binary class synthetic, MNIST and Fashion MNIST datasets, our schemes outperform the existing methods and demonstrate consistent performance across noise rates.


title: Variograms for spatial functional data with phase variation

abstract: Spatial, amplitude and phase variations in spatial functional data are confounded. Conclusions from the popular functional trace variogram, which quantifies spatial variation, can be misleading when analysing misaligned functional data with phase variation. To remedy this, we describe a framework that extends amplitude-phase separation methods in functional data to the spatial setting, with a view towards performing clustering and spatial prediction. We propose a decomposition of the trace variogram into amplitude and phase components and quantify how spatial correlations between functional observations manifest in their respective amplitude and phase components. This enables us to generate separate amplitude and phase clustering methods for spatial functional data, and develop a novel spatial functional interpolant at unobserved locations based on combining separate amplitude and phase predictions. Through simulations and real data analyses, we found that the proposed methods result in more accurate predictions and more interpretable clustering results.


title: Towards cosmological constraints from the compressed modal bispectrum: a robust comparison of real-space bispectrum estimators

abstract: Higher-order clustering statistics, like the galaxy bispectrum, can add complementary cosmological information to what is accessible with two-point statistics, like the power spectrum. While the standard way of measuring the bispectrum involves estimating a bispectrum value in a large number of Fourier triangle bins, the compressed modal bispectrum approximates the bispectrum as a linear combination of basis functions and estimates the expansion coefficients on the chosen basis. In this work, we compare the two estimators by using parallel pipelines to analyze the real-space halo bispectrum measured in a suite of $N$-body simulations corresponding to a total volume of $\sim 1{,}000 \,h^{-3}\,{\rm Gpc}^3$, with covariance matrices estimated from 10,000 mock halo catalogs. We find that the modal bispectrum yields constraints that are consistent and competitive with the standard bispectrum analysis: for the halo bias and shot noise parameters within the tree-level halo bispectrum model up to $k_{\rm max} \approx 0.06 \, (0.10) \,h\,{\rm Mpc}^{-1}$, only 6 (10) modal expansion coefficients are necessary to obtain constraints equivalent to the standard bispectrum estimator using $\sim$ 20 to 1,600 triangle bins, depending on the bin width. For this work, we have implemented a modal estimator pipeline using Markov Chain Monte Carlo simulations for the first time, and we discuss in detail how the parameter posteriors and modal expansion are robust to, or sensitive to, several user settings within the modal bispectrum pipeline. The combination of the highly efficient compression that is achieved and the large number of mock catalogs available allows us to quantify how our modal bispectrum constraints depend on the number of mocks that are used to estimate covariance matrices and the functional form of the likelihood.


title: EPTAS for $k$-means Clustering of Affine Subspaces

abstract: We consider a generalization of the fundamental $k$-means clustering for data with incomplete or corrupted entries. When data objects are represented by points in $\mathbb{R}^d$, a data point is said to be incomplete when some of its entries are missing or unspecified. An incomplete data point with at most $\Delta$ unspecified entries corresponds to an axis-parallel affine subspace of dimension at most $\Delta$, called a $\Delta$-point. Thus we seek a partition of $n$ input $\Delta$-points into $k$ clusters minimizing the $k$-means objective. For $\Delta=0$, when all coordinates of each point are specified, this is the usual $k$-means clustering. We give an algorithm that finds an $(1+ \epsilon)$-approximate solution in time $f(k,\epsilon, \Delta) \cdot n^2 \cdot d$ for some function $f$ of $k,\epsilon$, and $\Delta$ only.


title: Characterizations of weighted (b,c) inverse

abstract: The notion of weighted $(b,c)$-inverse of an element in rings were introduced, very recently [Comm. Algebra, 48 (4) (2020): 1423-1438]. In this paper, we further elaborate on this theory by establishing a few characterizations of this inverse and their relationships with other $(v, w)$-weighted $(b,c)$-inverses. We introduce some necessary and sufficient conditions for the existence of the hybrid $(v, w)$-weighted $(b,c)$-inverse and annihilator $(v, w)$-weighted $(b,c)$-inverse of elements in rings. In addition to this, we explore a few sufficient conditions for the reverse-order law of the annihilator $(v, w)$-weighted $(b,c)$-inverses.


title: Learning to Reconstruct and Segment 3D Objects

abstract: To endow machines with the ability to perceive the real-world in a three dimensional representation as we do as humans is a fundamental and long-standing topic in Artificial Intelligence. Given different types of visual inputs such as images or point clouds acquired by 2D/3D sensors, one important goal is to understand the geometric structure and semantics of the 3D environment. Traditional approaches usually leverage hand-crafted features to estimate the shape and semantics of objects or scenes. However, they are difficult to generalize to novel objects and scenarios, and struggle to overcome critical issues caused by visual occlusions. By contrast, we aim to understand scenes and the objects within them by learning general and robust representations using deep neural networks, trained on large-scale real-world 3D data. To achieve these aims, this thesis makes three core contributions from object-level 3D shape estimation from single or multiple views to scene-level semantic understanding.


title: Warrior1: A Performance Sanitizer for C++

abstract: This paper presents Warrior1, a tool that detects performance anti-patterns in C++ libraries. Many programs are slowed down by many small inefficiencies. Large-scale C++ applications are large, complex, and developed by large groups of engineers over a long period of time, which makes the task of identifying inefficiencies difficult. Warrior1 was designed to detect the numerous small performance issues that are the result of inefficient use of C++ libraries. The tool detects performance anti-patterns such as map double-lookup, vector reallocation, short lived objects, and lambda object capture by value. Warrior1 is implemented as an instrumented C++ standard library and an off-line diagnostics tool. The tool is very effective in detecting issues. We demonstrate that the tool is able to find a wide range of performance anti-patterns in a number of popular performance sensitive open source projects.


title: Latency-aware and -predictable Communication with Open Protocol Stacks for Remote Drone Control

abstract: In order to create cooperating swarms of Unmanned Autonomous Vehicles (UAVs) that also interact with various other systems and devices, open and free communication systems are mandatory. This paper presents an implementation of such a communication system to incorporate the Crazyflie nano-drone as a UAV platform. The protocol stack leverages the open Predictably Reliable Real-time Transport (PRRT) protocol that adds latency-awareness and -predictability to stacks composed of standard Internet protocols. To enable the drone to receive and reply to control commands via Wi-Fi, it has been extended with a Raspberry Pi that runs two variants of the Crazybridge---a software to connect the control board to the network. To evaluate how practical this solution is for the use in control applications, the communication has been analysed with a focus on the latency properties. Our investigations show that despite using the open protocol stack---and hence opting out of specialised implementations---the resulting latencies are in the same order of magnitude (4 to 9 ms) as the latency of the proprietary link.


title: Parallel and Distributed algorithms for ML problems

abstract: In this paper we make a survey of modern parallel and distributed approaches to solve sum-type convex minimization problems come from ML applications.


title: Brain Atlas Guided Attention U-Net for White Matter Hyperintensity Segmentation

abstract: White Matter Hyperintensities (WMH) are the most common manifestation of cerebral small vessel disease (cSVD) on the brain MRI. Accurate WMH segmentation algorithms are important to determine cSVD burden and its clinical consequences. Most of existing WMH segmentation algorithms require both fluid attenuated inversion recovery (FLAIR) images and T1-weighted images as inputs. However, T1-weighted images are typically not part of standard clinicalscans which are acquired for patients with acute stroke. In this paper, we propose a novel brain atlas guided attention U-Net (BAGAU-Net) that leverages only FLAIR images with a spatially-registered white matter (WM) brain atlas to yield competitive WMH segmentation performance. Specifically, we designed a dual-path segmentation model with two novel connecting mechanisms, namely multi-input attention module (MAM) and attention fusion module (AFM) to fuse the information from two paths for accurate results. Experiments on two publicly available datasets show the effectiveness of the proposed BAGAU-Net. With only FLAIR images and WM brain atlas, BAGAU-Net outperforms the state-of-the-art method with T1-weighted images, paving the way for effective development of WMH segmentation. Availability:https://github.com/Ericzhang1/BAGAU-Net


title: Simulating Photodissociation Reactions in Bad Cavities with the Lindblad Equation

abstract: Optical cavities, e.g. as used in organic polariton experiments, often employ low finesse mirrors or plasmonic structures. The photon lifetime in these setups is comparable to the timescale of the nuclear dynamics governing the photochemistry. This highlights the need for including the effect of dissipation in the molecular simulations. In this study, we perform wave packet dynamics with the Lindblad master equation, to study the effect of a finite photon lifetime on the dissociation of the MgH$^+$ molecule model system. Photon lifetimes of several different orders of magnitude are considered to encompass an ample range of effects inherent to lossy cavities.


title: Continuous monitoring of cognitive load using advanced computerized analysis of brain signals during virtual simulator training for laparoscopic surgery, reflects laparoscopic dexterity. A comparative study using a novel wireless device

abstract: Simulation-based training is an effective tool for acquiring practical skills, specifically to train new surgeons in a controlled and hazard-free environment, it is however important to measure participants cognitive load to decide whether they are ready to go into a real surgery. In the present study we measured performance on a surgery simulator of medical students and interns, while their brain activity was monitored by a mobile EEG device. 38 medical studentswere underwent 3 experiments undergoing a task with Simbionix simulator, while their brain activity was measured using a single-channel EEG device (Aurora by Neurosteer). On each experiment, participants performed 3 repeats of a simulator task using laparoscopic hands. The retention between tasks was different on each experiment, to examine changes in performance and cognitive load biomarkers that occur during the task or as a results of night sleep consolidation. The participants behavioral performance improved with trial repetition in all 3 experiments. In Exps. 1 & 2, the theta band activity significantly decreased with better individual performance, as exhibited by some of the behavioral measurements of the simulator. The novel VC9 biomarker (previously shown to correlate with cognitive load), exhibited a significant decrease with better individual performance shown by all behavioral measurements. In correspondence with previous research, theta decreased with lower cognitive load and higher performance and the novel biomarker, VC9, showed higher sensitivity to load changes. Together, these measurements might be for neuroimaging assessment of cognitive load while performing simulator laparoscopic tasks. This could potentially be expanded to evaluate efficacy of different medical simulations to provide more efficient training to medical staff and to measure cognitive and mental load in real laparoscopic surgeries.


title: A factor-adjusted multiple testing of general alternatives

abstract: Factor-adjusted multiple testing is used for handling strong correlated tests. Since most of previous works control the false discovery rate under sparse alternatives, we develop a two-step method, namely the AdaFAT, for any true false proportion. In this paper, the proposed procedure is adjusted by latent factor loadings. Under the existence of explanatory variables, a uniform convergence rate of the estimated factor loadings is given. We also show that the power of AdaFAT goes to one along with the controlled false discovery rate. The performance of the proposed procedure is examined through simulations calibrated by China A-share market.


title: Investigation on structural, electronic and magnetic properties of Co2FeGe Heusler alloy: experiment and theory

abstract: Experimental and computational studies were performed on Co2FeGe Heusler alloy. It was found that the alloy has very high experimental magnetic moment of 6.1 muB/f.u., curie temperature of 1073K and very high spin-wave stiffness constant of 10.4 nm2-meV, which indicates that the magnetic moment is very high and do not vary with change in temperature in the range 0-300K. The alloy strictly follows Slater-Pauling (SP) rule and the minor experimental deviation from its SP value is justified by doing full-potential density functional calculations which gives more accurate result when electron-electron correlation parameter (U) is taken into account with conventional GGA method. Effect of lattice strain and electron correlation on individual atomic moments, total magnetic moment and spin-polarization is studied in detail and can be concluded that they have a role in the deviation of the experimental results from the expected theoretical values.


title: Subdomain Separability in Global Optimization

abstract: We propose a generalization of separability in the context of global optimization. Our results apply to objective functions implemented as differentiable computer programs. They are presented in the context of a simple branch and bound method. The often significant search space reduction can be expected to yield an acceleration of any global optimization method. We show how to utilize interval derivatives calculated by adjoint algorithmic differentiation to examine the monotonicity of the objective with respect to so called structural separators and how to verify the latter automatically.


title: The scaling limit of the directed polymer with power-law tail disorder

abstract: In this paper, we study the so-called intermediate disorder regime for a directed polymer in a random environment with heavy-tail. Consider a simple symmetric random walk $(S_n)_{n\geq 0}$ on $\mathbb{Z}^d$, with $d\geq 1$, and modify its law using Gibbs weights in the product form $\prod_{n=1}^{N} (1+\beta\eta_{n,S_n})$, where $(\eta_{n,x})_{n\ge 0, x\in \mathbb{Z}^d}$ is a field of i.i.d. random variables whose distribution satisfies $\mathbb{P}(\eta>z) \sim z^{-\alpha}$ as $z\to\infty$, for some $\alpha\in(0,2)$. We prove that if $\alpha< \min(1+\frac{2}{d},2)$, when sending $N$ to infinity and rescaling the disorder intensity by taking $\beta=\beta_N \sim N^{-\gamma}$ with $\gamma =\frac{d}{2\alpha}(1+\frac{2}{d}-\alpha)$, the distribution of the trajectory under diffusive scaling converges in law towards a random limit, which is the continuum polymer with L\'evy $\alpha$-stable noise constructed in the companion paper arXiv:2007.06484.


title: Multiclass Wound Image Classification using an Ensemble Deep CNN-based Classifier

abstract: Acute and chronic wounds are a challenge to healthcare systems around the world and affect many people's lives annually. Wound classification is a key step in wound diagnosis that would help clinicians to identify an optimal treatment procedure. Hence, having a high-performance classifier assists the specialists in the field to classify the wounds with less financial and time costs. Different machine learning and deep learning-based wound classification methods have been proposed in the literature. In this study, we have developed an ensemble Deep Convolutional Neural Network-based classifier to classify wound images including surgical, diabetic, and venous ulcers, into multi-classes. The output classification scores of two classifiers (patch-wise and image-wise) are fed into a Multi-Layer Perceptron to provide a superior classification performance. A 5-fold cross-validation approach is used to evaluate the proposed method. We obtained maximum and average classification accuracy values of 96.4% and 94.28% for binary and 91.9\% and 87.7\% for 3-class classification problems. The results show that our proposed method can be used effectively as a decision support system in classification of wound images or other related clinical applications.


title: Multi-Modal Super Resolution for Dense Microscopic Particle Size Estimation

abstract: Particle Size Analysis (PSA) is an important process carried out in a number of industries, which can significantly influence the properties of the final product. A ubiquitous instrument for this purpose is the Optical Microscope (OM). However, OMs are often prone to drawbacks like low resolution, small focal depth, and edge features being masked due to diffraction. We propose a powerful application of a combination of two Conditional Generative Adversarial Networks (cGANs) that Super Resolve OM images to look like Scanning Electron Microscope (SEM) images. We further demonstrate the use of a custom object detection module that can perform efficient PSA of the super-resolved particles on both, densely and sparsely packed images. The PSA results obtained from the super-resolved images have been benchmarked against human annotators, and results obtained from the corresponding SEM images. The proposed models show a generalizable way of multi-modal image translation and super-resolution for accurate particle size estimation.


title: A Geometrical Model for the Evolution of Spherical Planetary Nebulae Based on Thin-Shell Formalism

abstract: A spherical planetary nebula is described as a geometric model. The nebula itself is considered as a thin-shell which visualized as a boundary of two spacetimes. The inner and outer curvature tensors of the thin-shell are found in order to get an expression of the energy-momentum tensor on the thin-shell. The energy density and pressure expressions are derived using the energy-momentum tensor. The time evolution of the radius of the thin-shell is obtained in terms of the energy density. The model is tested by using a simple power function for decreasing energy density and the evolution pattern of the planetary nebula is attained.


title: Stochastic recursions on directed random graphs

abstract: For a directed graph $G(V_n, E_n)$ on the vertices $V_n = \{1,2, \dots, n\}$, we study the distribution of a Markov chain $\{ {\bf R}^{(k)}: k \geq 0\}$ on $\mathbb{R}^n$ such that the $i$th component of ${\bf R}^{(k)}$, denoted $R_i^{(k)}$, corresponds to the value of the process on vertex $i$ at time $k$. We focus on processes $\{ {\bf R}^{(k)}: k \geq 0\}$ where the value of $R_i^{(k+1)}$ depends only on the values $\{ R_j^{(k)}: j \to i\}$ of its inbound neighbors, and possibly on vertex attributes. We then show that, provided $G(V_n, E_n)$ converges in the local weak sense to a marked Galton-Watson process, the dynamics of the process for a uniformly chosen vertex in $V_n$ can be coupled, for any fixed $k$, to a process $\{ \mathcal{R}_\emptyset^{(r)}: 0 \leq r \leq k\}$ constructed on the limiting marked Galton-Watson tree. Moreover, we derive sufficient conditions under which $\mathcal{R}^{(k)}_\emptyset$ converges, as $k \to \infty$, to a random variable $\mathcal{R}^*$ that can be characterized in terms of the attracting endogenous solution to a branching distributional fixed-point equation. Our framework can also be applied to processes $\{ {\bf R}^{(k)}: k \geq 0\}$ whose only source of randomness comes from the realization of the graph $G(V_n, E_n)$.


title: Faster Convergence of Stochastic Gradient Langevin Dynamics for Non-Log-Concave Sampling

abstract: We provide a new convergence analysis of stochastic gradient Langevin dynamics (SGLD) for sampling from a class of distributions that can be non-log-concave. At the core of our approach is a novel conductance analysis of SGLD using an auxiliary time-reversible Markov Chain. Under certain conditions on the target distribution, we prove that $\tilde O(d^4\epsilon^{-2})$ stochastic gradient evaluations suffice to guarantee $\epsilon$-sampling error in terms of the total variation distance, where $d$ is the problem dimension. This improves existing results on the convergence rate of SGLD (Raginsky et al., 2017; Xu et al., 2018). We further show that provided an additional Hessian Lipschitz condition on the log-density function, SGLD is guaranteed to achieve $\epsilon$-sampling error within $\tilde O(d^{15/4}\epsilon^{-3/2})$ stochastic gradient evaluations. Our proof technique provides a new way to study the convergence of Langevin-based algorithms and sheds some light on the design of fast stochastic gradient-based sampling algorithms.


title: Better Distractions: Transformer-based Distractor Generation and Multiple Choice Question Filtering

abstract: For the field of education, being able to generate semantically correct and educationally relevant multiple choice questions (MCQs) could have a large impact. While question generation itself is an active research topic, generating distractors (the incorrect multiple choice options) receives much less attention. A missed opportunity, since there is still a lot of room for improvement in this area. In this work, we train a GPT-2 language model to generate three distractors for a given question and text context, using the RACE dataset. Next, we train a BERT language model to answer MCQs, and use this model as a filter, to select only questions that can be answered and therefore presumably make sense. To evaluate our work, we start by using text generation metrics, which show that our model outperforms earlier work on distractor generation (DG) and achieves state-of-the-art performance. Also, by calculating the question answering ability, we show that larger base models lead to better performance. Moreover, we conducted a human evaluation study, which confirmed the quality of the generated questions, but showed no statistically significant effect of the QA filter.


title: On the Combinatorics of Placing Balls into Ordered Bins

abstract: In this paper, we use techniques of enumerative combinatorics to study the following problem: we count the number of ways to split $n$ balls into nonempty, ordered bins so that the most crowded bin has exactly $k$ balls. We find closed forms for three of the different cases that can arise: $k > \frac{n}{2}$, $k = \frac{n}{2}$, and when there exists $j < k$ such that $n = 2k + j$. As an immediate result of our proofs, we find a closed form for the number of positive integer solutions to $x_1 + x_2 + \dots + x_{\ell} = n$ with the attained maximum of $\{x_1, x_2, \dots, x_{\ell}\}$ being equal to $k$, when $n$ and $k$ have one of the aforementioned algebraic relationships to each other. The problem is generalized to find a formula that enumerates the total number of ways without specific conditions on $n, \ell, k$. Subsequently, various additional identities and estimates related to this enumeration are proven and interpreted.


title: Drug Repurposing for COVID-19 via Knowledge Graph Completion

abstract: Objective: To discover candidate drugs to repurpose for COVID-19 using literature-derived knowledge and knowledge graph completion methods. Methods: We propose a novel, integrative, and neural network-based literature-based discovery (LBD) approach to identify drug candidates from both PubMed and COVID-19-focused research literature. Our approach relies on semantic triples extracted using SemRep (via SemMedDB). We identified an informative subset of semantic triples using filtering rules and an accuracy classifier developed on a BERT variant, and used this subset to construct a knowledge graph. Five SOTA, neural knowledge graph completion algorithms were used to predict drug repurposing candidates. The models were trained and assessed using a time slicing approach and the predicted drugs were compared with a list of drugs reported in the literature and evaluated in clinical trials. These models were complemented by a discovery pattern-based approach. Results: Accuracy classifier based on PubMedBERT achieved the best performance (F1= 0.854) in classifying semantic predications. Among five knowledge graph completion models, TransE outperformed others (MR = 0.923, Hits@1=0.417). Some known drugs linked to COVID-19 in the literature were identified, as well as some candidate drugs that have not yet been studied. Discovery patterns enabled generation of plausible hypotheses regarding the relationships between the candidate drugs and COVID-19. Among them, five highly ranked and novel drugs (paclitaxel, SB 203580, alpha 2-antiplasmin, pyrrolidine dithiocarbamate, and butylated hydroxytoluene) with their mechanistic explanations were further discussed. Conclusion: We show that an LBD approach can be feasible for discovering drug candidates for COVID-19, and for generating mechanistic explanations. Our approach can be generalized to other diseases as well as to other clinical questions.


title: Neuronal graphs: a graph theory primer for microscopic, functional networks of neurons recorded by Calcium imaging

abstract: Connected networks are a fundamental structure of neurobiology. Understanding these networks will help us elucidate the neural mechanisms of computation. Mathematically speaking these networks are `graphs' - structures containing objects that are connected. In neuroscience, the objects could be regions of the brain, e.g. fMRI data, or be individual neurons, e.g. calcium imaging with fluorescence microscopy. The formal study of graphs, graph theory, can provide neuroscientists with a large bank of algorithms for exploring networks. Graph theory has already been applied in a variety of ways to fMRI data but, more recently, has begun to be applied at the scales of neurons, e.g. from functional calcium imaging. In this primer we explain the basics of graph theory and relate them to features of microscopic functional networks of neurons from calcium imaging - neuronal graphs. We explore recent examples of graph theory applied to calcium imaging and we highlight some areas where researchers new to the field could go awry.


title: End-to-End Text-to-Speech using Latent Duration based on VQ-VAE

abstract: Explicit duration modeling is a key to achieving robust and efficient alignment in text-to-speech synthesis (TTS). We propose a new TTS framework using explicit duration modeling that incorporates duration as a discrete latent variable to TTS and enables joint optimization of whole modules from scratch. We formulate our method based on conditional VQ-VAE to handle discrete duration in a variational autoencoder and provide a theoretical explanation to justify our method. In our framework, a connectionist temporal classification (CTC) -based force aligner acts as the approximate posterior, and text-to-duration works as the prior in the variational autoencoder. We evaluated our proposed method with a listening test and compared it with other TTS methods based on soft-attention or explicit duration modeling. The results showed that our systems rated between soft-attention-based methods (Transformer-TTS, Tacotron2) and explicit duration modeling-based methods (Fastspeech).


title: Limits on the Existence of sub-MeV Sterile Neutrinos from the Decay of $^7$Be in Superconducting Quantum Sensors

abstract: Sterile neutrinos are natural extensions to the standard model of particle physics and provide a possible portal to the dark sector. We report a new search for the existence of sub-MeV sterile neutrinos using the decay-momentum reconstruction technique in the decay of $^7$Be. The experiment measures the total energy of the $^7$Li daughter atom from the electron capture decay of $^7$Be implanted into sensitive superconducting tunnel junction (STJ) quantum sensors. This first experiment presents data from a single STJ operated at a low count rate for a net total of 28 days, and provides exclusion limits on sterile neutrinos in the mass range from 100 to 850 keV that improve upon previous work by up to an order of magnitude.


title: Nonreciprocal transition between two indirectly coupled energy levels

abstract: We propose a theoretical scheme to realize nonreciprocal transition between two energy levels that can not coupled directly. Suppose they are coupled indirectly by two auxiliary levels with a cyclic four-level configuration, and the four transitions in the cyclic configuration are controlled by external fields. The indirectly transition become nonreciprocal when the time reversal symmetry of the system is broken by the synthetic magnetic flux, i.e., the total phase of the external driving fields through the cyclic four-level configuration. The nonreciprocal transition can be identified by the elimination of a spectral line in the spontaneous emission spectrum. Our work introduces a feasible way to observe nonreciprocal transition in a wide range of multi-level systems, including natural atoms or ions with parity symmetry.


title: Inspection-on-the-fly using Hybrid Physical Interaction Control for Aerial Manipulators

abstract: Inspection for structural properties (surface stiffness and coefficient of restitution) is crucial for understanding and performing aerial manipulations in unknown environments, with little to no prior knowledge on their state. Inspection-on-the-fly is the uncanny ability of humans to infer states during manipulation, reducing the necessity to perform inspection and manipulation separately. This paper presents an infrastructure for inspection-on-the-fly method for aerial manipulators using hybrid physical interaction control. With the proposed method, structural properties (surface stiffness and coefficient of restitution) can be estimated during physical interactions. A three-stage hybrid physical interaction control paradigm is presented to robustly approach, acquire and impart a desired force signature onto a surface. This is achieved by combining a hybrid force/motion controller with a model-based feed-forward impact control as intermediate phase. The proposed controller ensures a steady transition from unconstrained motion control to constrained force control, while reducing the lag associated with the force control phase. And an underlying Operational Space dynamic configuration manager permits complex, redundant vehicle/arm combinations. Experiments were carried out in a mock-up of a Dept. of Energy exhaust shaft, to show the effectiveness of the inspection-on-the-fly method to determine the structural properties of the target surface and the performance of the hybrid physical interaction controller in reducing the lag associated with force control phase.


title: Wavefronts and modal structure of long surface and internal ring waves on a parallel shear current

abstract: We study long surface and internal ring waves propagating in a stratified fluid over a parallel shear flow. The far-field modal and amplitude equations for the ring waves are presented in dimensional form. We re-derive them from the formulation for plane waves tangent to the ring wave, which opens a way to obtaining important characteristics of the ring waves (e.g. group speed) and to constructing more general `hybrid solutions' consisting of a part of a ring wave and two tangent plane waves. The modal equations constitute a new spectral problem, and are analysed for a number of examples of surface ring waves in a homogeneous fluid and internal ring waves in a stratified fluid. The detailed analysis is developed for the case of a two-layered fluid with a linear shear current where we study their wavefronts and two-dimensional modal structure. Comparisons are made between the modal functions of the surface waves in a homogeneous and two-layered fluids, as well as the interfacial waves described exactly and in the rigid-lid approximation. We also analyse the wavefronts of surface and interfacial waves for a family of power-law upper-layer currents, which can be used to model wind generated currents, river inflows and exchange flows in straits. A global and local measure of the deformation of wavefronts are introduced and evaluated.


title: Evaluation of Energy Resolution by Changing Angle and Position of Incident Photon in a LYSO Calorimeter

abstract: In this paper, we investigate the effect on energy resolution from changing the angle and the position of incidence photon for a 5 $\times$ 5 crystal matrix built with (25 $\times$ 25 $\times$ 200) mm$^{3}$ LYSO scintillating crystals. Those crystals have been proposed for the electromagnetic calorimeter of the Turkish Accelerator Center-Particle Factory (TAC-PF) detector. The energy resolution was obtained as $\sigma_{E}/E = 0.42 \% / \sqrt{E/GeV} \oplus 1.60 \%$ at the center of the matrix in the energy range of 50 MeV to 2 GeV. When we examined the dependence of the energy resolution on the incidence angle of the photon, resolution began to deteriorate at angles greater than $2^{\circ}$ on the 5 $\times$ 5 crystal matrix. Moreover, energy resolution at the corners of the central crystal was worse than at the center of the central crystal by a factor of 1.3 at 50 MeV and 1.1 at 2 GeV.


title: Incorporating Terminology Constraints in Automatic Post-Editing

abstract: Users of machine translation (MT) may want to ensure the use of specific lexical terminologies. While there exist techniques for incorporating terminology constraints during inference for MT, current APE approaches cannot ensure that they will appear in the final translation. In this paper, we present both autoregressive and non-autoregressive models for lexically constrained APE, demonstrating that our approach enables preservation of 95% of the terminologies and also improves translation quality on English-German benchmarks. Even when applied to lexically constrained MT output, our approach is able to improve preservation of the terminologies. However, we show that our models do not learn to copy constraints systematically and suggest a simple data augmentation technique that leads to improved performance and robustness.


title: Fermi Liquid Theory Sheds Light on "Hot" EHL in 1L-MoS$_2$

abstract: 2D transition metal dichalcogenides (TMDCs) exhibit an electron-hole liquid phase transition at unusually high temperatures. Because these materials are atomically thin, optical excitation leads to material expansion. As a result, during the EHL phase transition the electronic band structure evolves due to both material thermal expansion and renormalization of the bands under high excitation densities. Specifically, these effects lead to indirect gap electronic band structure with a valence band maximum located at the $\Gamma$ valley. In this work we developed a methodology for analyzing the spectral evolution of the photoluminescence of suspended 1L-MoS$_2$ during the EHL phase transition by using Fermi liquid theory. The resulting analysis reveals valley-specific carrier densities, radiative recombination efficiencies, and intraband carrier relaxation kinetics in 1L-MoS$_2$. More broadly, the results outline a methodology for predicting critical EHL parameters, shedding light onto the EHL phase transition in 2D TDMCs.


title: Increasing Depth Leads to U-Shaped Test Risk in Over-parameterized Convolutional Networks

abstract: Recent works have demonstrated that increasing model capacity through width in over-parameterized neural networks leads to a decrease in test risk. For neural networks, however, model capacity can also be increased through depth, yet understanding the impact of increasing depth on test risk remains an open question. In this work, we demonstrate that the test risk of over-parameterized convolutional networks is a U-shaped curve (i.e. monotonically decreasing, then increasing) with increasing depth. We first provide empirical evidence for this phenomenon via image classification experiments using both ResNets and the convolutional neural tangent kernel (CNTK). We then present a novel linear regression framework for characterizing the impact of depth on test risk, and show that increasing depth leads to a U-shaped test risk for the linear CNTK. In particular, we prove that the linear CNTK corresponds to a depth-dependent linear transformation on the original space and characterize properties of this transformation. We then analyze over-parameterized linear regression under arbitrary linear transformations and, in simplified settings, provably identify the depths which minimize each of the bias and variance terms of the test risk.


title: Noise-driven Topological Changes in Chaotic Dynamics

abstract: Noise modifies the behavior of chaotic systems in both quantitative and qualitative ways. To study these modifications, the present work compares the topological structure of the deterministic Lorenz (1963) attractor with its stochastically perturbed version. The deterministic attractor is well known to be "strange" but it is frozen in time. When driven by multiplicative noise, the Lorenz model's random attractor (LORA) evolves in time. Algebraic topology sheds light on the most striking effects involved in such an evolution. In order to examine the topological structure of the snapshots that approximate LORA, we use Branched Manifold Analysis through Homologies (BraMAH) -- a technique originally introduced to characterize the topological structure of deterministically chaotic flows -- which is being extended herein to nonlinear noise-driven systems. The analysis is performed for a fixed realization of the driving noise at different time instants in time. The results suggest that LORA's evolution includes sharp transitions that appear as topological tipping points.


title: Correlation Functions for a Chain of Short Range Oscillators

abstract: We consider a system of harmonic oscillators with short range interactions and we study their correlation functions when the initial data is sampled with respect to the Gibbs measure. Such correlation functions display rapid oscillations that travel through the chain. We show that the correlation functions always have two fastest peaks which move in opposite directions and decay at rate $t^{-\frac{1}{3}}$ for position and momentum correlations and as $t^{-\frac{2}{3}}$ for energy correlations. The shape of these peaks is asymptotically described by the Airy function. Furthermore, the correlation functions have some non generic peaks with lower decay rates. In particular, there are peaks which decay at rate $t^{-\frac{1}{4}}$ for position and momentum correlators and with rate $t^{-\frac{1}{2}}$ for energy correlators. The shape of these peaks is described by the Pearcey integral. Crucial for our analysis is an appropriate generalisation of spacings, i.e. differences of the positions of neighbouring particles, that are used as spatial variables in the case of nearest neighbour interactions. Using the theory of circulant matrices we are able to introduce a quantity that retains both localisation and analytic viability. This also allows us to define and analyse some additional quantities used for nearest neighbour chains. Finally, we study numerically the evolution of the correlation functions after adding nonlinear perturbations to our model. Within the time range of our numerical simulations the asymptotic description of the linear case seems to persist for small nonlinear perturbations while stronger nonlinearities change shape and decay rates of the peaks significantly.

categories: math-ph cond-mat.stat-mech math.MP math.PR math.SP


title: Detection and Prediction of Equilibrium States in Kinetic Plasma Simulations via Mode Tracking using Reduced-Order Dynamic Mode Decomposition

abstract: A dynamic mode decomposition (DMD) based reduced-order model (ROM) is developed for tracking, detection, and prediction of kinetic plasma behavior. DMD is applied to the high-fidelity kinetic plasma model based on the electromagnetic particle-in-cell (EMPIC) algorithm to extract the underlying dynamics and key features of the model. In particular, the ability of DMD to reconstruct the spatial pattern of the self electric field from high-fidelity data and the effect of DMD extrapolated self-fields on charged particle dynamics are investigated. An in-line sliding-window DMD method is presented for identifying the transition from transient to equilibrium state based on the loci of DMD eigenvalues in the complex plane. The in-line detection of equilibrium state combined with time extrapolation ability of DMD has the potential to effectively expedite the simulation. Case studies involving electron beams and plasma ball are presented to assess the strengths and limitations of the proposed method.


title: Influence maximization on temporal networks

abstract: We consider the optimization problem of seeding a spreading process on a temporal network so that the expected size of the resulting outbreak is maximized. We frame the problem for a spreading process following the rules of the susceptible-infected-recovered model with temporal scale equal to the one characterizing the evolution of the network topology. We perform a systematic analysis based on a corpus of 12 real-world temporal networks and quantify the performance of solutions to the influence maximization problem obtained using different level of information about network topology and dynamics. We find that having perfect knowledge of the network topology but in a static and/or aggregated form is not helpful in solving the influence maximization problem effectively. Knowledge, even if partial, of the early stages of the network dynamics appears instead essential for the identification of quasioptimal sets of influential spreaders.


title: An upper bound on the topological complexity of discriminantal varieties

abstract: We give an upper bound on the topological complexity of varieties $\mathcal{V}$ obtained as complements in $\mathbb{C}^m$ of the zero locus of a polynomial. As an application, we determine the topological complexity of unordered configuration spaces of the plane.


title: Cooperative Multi-Sensor Detection under Variable-Length Coding

abstract: We investigate the testing-against-independence problem \mw{over a cooperative MAC} with two sensors and a single detector under an average rate constraint on the sensors-detector links. For this setup, we design a variable-length coding scheme that maximizes the achievable type-II error exponent when the type-I error probability is limited to $\epsilon$. Similarly to the single-link result, we show here that the optimal error exponent depends on $\epsilon$ and that variable-length coding allows to \mw{increase the rates over the optimal} fixed-length coding \mw{scheme} by the factor $(1-\epsilon)^{-1}$.


title: Brightest Cluster Galaxies: the centre can(not?) hold

abstract: We explore the persistence of the alignment of brightest cluster galaxies (BCGs) with their local environment. We find that a significant fraction of BCGs do not coincide with the centroid of the X-ray gas distribution and/or show peculiar velocities (they are not at rest with respect to the cluster mean). Despite this, we find that BCGs are generally aligned with the cluster mass distribution even when they have significant offsets from the X-ray centre and significant peculiar velocities. The large offsets are not consistent with simple theoretical models. To account for these observations BCGs must undergo mergers preferentially along their major axis, the main infall direction. Such BCGs may be oscillating within the cluster potential after having been displaced by mergers or collisions, or the dark matter halo itself may not yet be relaxed.


title: Aerial Mobile Manipulator System to Enable Dexterous Manipulations with Increased Precision

abstract: Problems associated with physical interactions using aerial mobile manipulators (AMM) are being independently addressed with respect to mobility and manipulability. Multirotor unmanned aerial vehicles (UAV) are a common choice for mobility while on-board manipulators are increasingly be used for manipulability. However, the dynamic coordination between the UAV and on-board manipulator remains a significant obstacle to enable dexterous manipulation with high precision. This paper presents an AMM system configuration to addresses both the mobility and manipulability issues together. A fully-actuated UAV is chosen to achieve dexterous aerial mobile manipulation, but is limited by the actuation range of the UAV. An on-board manipulator is employed to enhance the performance in terms of dexterity and precision at the end-effector. Experimental results on position keeping of the dexterous hexrotor by withstanding the disturbances caused by the motions of the on-board manipulator and external wind disturbances are presented. Preliminary simulation results on end-point tracking in a simple planar on-board manipulator case is presented.


title: Human psychophysical discrimination of spatially dependant Pancharatnam-Berry phases in optical spin-orbit states

abstract: We tested the ability of human observers to discriminate distinct profiles of spatially dependant geometric phases when directly viewing stationary structured light beams. Participants viewed polarization coupled orbital angular momentum (OAM) states, or ``spin-orbit'' states, in which the OAM was induced through Pancharatnam-Berry phases. The coupling between polarization and OAM in these beams manifests as spatially dependant polarization. Regions of uniform polarization are perceived as specifically oriented Haidinger's brushes, and study participants discriminated between two spin-orbit states based on the rotational symmetry in the spatial orientations of these brushes. Participants used self-generated eye movements to prevent adaptation to the visual stimuli. After initial training, the participants were able to correctly discriminate between two spin-orbit states, differentiated by OAM $=\pm1$, with an average success probability of $69\%$ ($S.D. = 22\%$, $p = 0.013$). These results support our previous observation that human observers can directly perceive spin-orbit states, and extend this finding to non-rotating beams, OAM modes induced via Pancharatnam-Berry phases, and the discrimination of states that are differentiated by OAM.


title: Spontaneous and explicit parity-time-symmetry breaking in drift wave instabilities

abstract: A method of Parity-Time (PT)-symmetry analysis is introduced to study the high dimensional, complicated parameter space of drift wave instabilities. We show that spontaneous PT-symmetry breaking leads to the Ion Temperature Gradient (ITG) instability of drift waves, and the collisional instability is the result of explicit PT-symmetry breaking. A new unstable drift wave induced by finite collisionality is identified. It is also found that gradients of ion temperature and density can destabilize the ion cyclotron waves when PT symmetry is explicitly broken by a finite collisionality.


title: Importance Reweighting for Biquality Learning

abstract: The field of Weakly Supervised Learning (WSL) has recently seen a surge of popularity, with numerous papers addressing different types of "supervision deficiencies", namely: poor quality, non adaptability, and insufficient quantity of labels. Regarding quality, label noise can be of different types, including completely-at-random, at-random or even not-at-random. All these kinds of label noise are addressed separately in the literature, leading to highly specialized approaches. This paper proposes an original, encompassing, view of Weakly Supervised Learning, which results in the design of generic approaches capable of dealing with any kind of label noise. For this purpose, an alternative setting called "Biquality data" is used. It assumes that a small trusted dataset of correctly labeled examples is available, in addition to an untrusted dataset of noisy examples. In this paper, we propose a new reweigthing scheme capable of identifying noncorrupted examples in the untrusted dataset. This allows one to learn classifiers using both datasets. Extensive experiments that simulate several types of label noise and that vary the quality and quantity of untrusted examples, demonstrate that the proposed approach outperforms baselines and state-of-the-art approaches.


title: Inferring respiratory and circulatory parameters from electrical impedance tomography with deep recurrent models

abstract: Electrical impedance tomography (EIT) is a noninvasive imaging modality that allows a continuous assessment of changes in regional bioimpedance of different organs. One of its most common biomedical applications is monitoring regional ventilation distribution in critically ill patients treated in intensive care units. In this work, we put forward a proof-of-principle study that demonstrates how one can reconstruct synchronously measured respiratory or circulatory parameters from the EIT image sequence using a deep learning model trained in an end-to-end fashion. We demonstrate that one can accurately infer absolute volume, absolute flow, normalized airway pressure and within certain limitations even the normalized arterial blood pressure from the EIT signal alone, in a way that generalizes to unseen patients without prior calibration. As an outlook with direct clinical relevance, we furthermore demonstrate the feasibility of reconstructing the absolute transpulmonary pressure from a combination of EIT and absolute airway pressure, as a way to potentially replace the invasive measurement of esophageal pressure. With these results, we hope to stimulate further studies building on the framework put forward in this work.


title: An Empirical Study for Vietnamese Constituency Parsing with Pre-training

abstract: In this work, we use a span-based approach for Vietnamese constituency parsing. Our method follows the self-attention encoder architecture and a chart decoder using a CKY-style inference algorithm. We present analyses of the experiment results of the comparison of our empirical method using pre-training models XLM-Roberta and PhoBERT on both Vietnamese datasets VietTreebank and NIIVTB1. The results show that our model with XLM-Roberta archived the significantly F1-score better than other pre-training models, VietTreebank at 81.19% and NIIVTB1 at 85.70%.


title: Optimism in the Face of Adversity: Understanding and Improving Deep Learning through Adversarial Robustness

abstract: Driven by massive amounts of data and important advances in computational resources, new deep learning systems have achieved outstanding results in a large spectrum of applications. Nevertheless, our current theoretical understanding on the mathematical foundations of deep learning lags far behind its empirical success. Towards solving the vulnerability of neural networks, however, the field of adversarial robustness has recently become one of the main sources of explanations of our deep models. In this article, we provide an in-depth review of the field of adversarial robustness in deep learning, and give a self-contained introduction to its main notions. But, in contrast to the mainstream pessimistic perspective of adversarial robustness, we focus on the main positive aspects that it entails. We highlight the intuitive connection between adversarial examples and the geometry of deep neural networks, and eventually explore how the geometric study of adversarial examples can serve as a powerful tool to understand deep learning. Furthermore, we demonstrate the broad applicability of adversarial robustness, providing an overview of the main emerging applications of adversarial robustness beyond security. The goal of this article is to provide readers with a set of new perspectives to understand deep learning, and to supply them with intuitive tools and insights on how to use adversarial robustness to improve it.


title: LoRa Performance Analysis with Superposed Signal Decoding

abstract: This paper considers the use of successive interference cancellation (SIC) to decode superposed signals in Long Range (LoRa) networks. We build over a known stochastic geometry model for LoRa networks and include the effect of recovering colliding packets through SIC. We derive closed-form expressions for the successful decoding of packets using SIC taking path loss, fading, noise, and interference into account, while we validate the model by means of Monte Carlo simulations. Results show that SIC-enabled LoRa networks improve worst-case reliability by up to 34%. We show that, for at least one test scenario, SIC increases by 159% the number of served users with the same worst-case reliability level.


title: Subsystem codes with high thresholds by gauge fixing and reduced qubit overhead

abstract: We introduce a technique that uses gauge fixing to significantly improve the quantum error correcting performance of subsystem codes. By changing the order in which check operators are measured, valuable additional information can be gained, and we introduce a new method for decoding which uses this information to improve performance. Applied to the subsystem toric code with three-qubit check operators, we increase the threshold under circuit-level depolarising noise from $0.67\%$ to $0.81\%$. The threshold increases further under a circuit-level noise model with small finite bias, up to $2.22\%$ for infinite bias. Furthermore, we construct families of finite-rate subsystem LDPC codes with three-qubit check operators and optimal-depth parity-check measurement schedules. To the best of our knowledge, these finite-rate subsystem codes outperform all known codes at circuit-level depolarising error rates as high as $0.2\%$, where they have a qubit overhead that is $4.3\times$ lower than the most efficient version of the surface code and $5.1\times$ lower than the subsystem toric code. Their threshold and pseudo-threshold exceeds $0.42\%$ for circuit-level depolarising noise, increasing to $2.4\%$ under infinite bias using gauge fixing.


title: Probabilistic Stirling numbers of the second kind and applications

abstract: Associated to each complex-valued random variable satisfying appropriate integrability conditions, we introduce a different generalization of the Stirling numbers of the second kind. Various equivalent definitions are provided. Attention, however, is focused on applications. Indeed, such numbers describe the moments of sums of i.i.d. random variables, determining their precise asymptotic behavior without making use of the central limit theorem. Such numbers also allow us to obtain explicit and simple Edgeworth expansions. Applications to L\'{e}vy processes and cumulants are discussed, as well.


title: Stability of 2-Parameter Persistent Homology

abstract: The \v{C}ech and Rips constructions of persistent homology are stable with respect to perturbations of the input data. However, neither is robust to outliers, and both can be insensitive to topological structure of high-density regions of the data. A natural solution is to consider 2-parameter persistence. This paper studies the stability of 2-parameter persistent homology: We show that several related density-sensitive constructions of bifiltrations from data satisfy stability properties accommodating the addition and removal of outliers. Specifically, we consider the multicover bifiltration, Sheehy's subdivision bifiltrations, and the degree bifiltrations. For the multicover and subdivision bifiltrations, we get 1-Lipschitz stability results closely analogous to the standard stability results for 1-parameter persistent homology. Our results for the degree bifiltrations are weaker, but they are tight, in a sense. As an application of our theory, we prove a law of large numbers for subdivision bifiltrations of random data.


title: PAC$^m$-Bayes: Narrowing the Empirical Risk Gap in the Misspecified Bayesian Regime

abstract: The Bayesian posterior minimizes the "inferential risk" which itself bounds the "predictive risk". This bound is tight when the likelihood and prior are well-specified. However since misspecification induces a gap, the Bayesian posterior predictive distribution may have poor generalization performance. This work develops a multi-sample loss (PAC$^m$) which can close the gap by spanning a trade-off between the two risks. The loss is computationally favorable and offers PAC generalization guarantees. Empirical study demonstrates improvement to the predictive distribution.


title: Super Cartan geometry and the super Ashtekar connection

abstract: This work is devoted to the geometric approach to supergravity. More precisely, we interpret $\mathcal{N}=1$, $D=4$ supergravity as a super Cartan geometry which provides a link between supergravity and Yang-Mills theory. To this end, we first review important aspects of the theory of supermanifolds and we establish a link between various different approaches. We then introduce super Cartan geometries using the concept of so-called enriched categories. This, among other things, will enable us to implement anticommutative fermionic fields. We will then also show that non-extended $D=4$ supergravity naturally arises in this framework. Finally, using this gauge-theoretic interpretation as well as the chiral structure of the underlying supersymmetry algebra, we will derive graded analoga of Ashtekar's self-dual variables and interpret them in terms of generalized super Cartan connections. This gives canonical chiral supergravity the structure of a Yang-Mills theory with gauge supergroup similar to the self-dual variables in ordinary first-order Einstein gravity which was first observed by F\"ul\"op \cite{Fulop:1993wi}. We then construct the parallel transport map corresponding to the super connection in mathematical rigorous way using again enriched categories. This provides the possibility of quantizing gravity and matter degrees of freedom in loop quantum gravity in a unified way.


title: Topological phase transition induced by band structure modulation in a Chern insulator

abstract: Here we study the systematic evolution of the topological properties of a Chern insulator in presence of an electronic dispersion that can be tuned smoothly from being Dirac-like till a semi-Dirac one and beyond. The band structure under such controlled deformation shows that the two Dirac points approach each other, merge at an intermediate point (the ${\mathbf{M}}$ point), where the low energy spectrum turns gapless, shows anisotropic Dirac features in the $k$-space and is denoted as the semi-Dirac limit, however a gap eventually opens up again in the spectrum. The Chern number phase diagram obtained via integrating the Berry curvature over the Brillouin zone (BZ) shows a gradual shrinking of the 'topological' lobes, and vanishes just beyond the semi-Dirac limit of the electronic dispersion. Thus there is a phase transition from a topological phase to a trivial phase across the semi-Dirac point. The vanishing of the anomalous Hall conductivity plateau and the merger of the chiral edge states with the bulk bands near the ${\mathbf{M}}$ point provide robust support of the observed phase transition.