title: $R^*$: A robust MCMC convergence diagnostic with uncertainty using decision tree classifiers

abstract: Markov chain Monte Carlo (MCMC) has transformed Bayesian model inference over the past three decades: mainly because of this, Bayesian inference is now a workhorse of applied scientists. Under general conditions, MCMC sampling converges asymptotically to the posterior distribution, but this provides no guarantees about its performance in finite time. The predominant method for monitoring convergence is to run multiple chains and monitor individual chains' characteristics and compare these to the population as a whole: if within-chain and between-chain summaries are comparable, then this is taken to indicate that the chains have converged to a common stationary distribution. Here, we introduce a new method for diagnosing convergence based on how well a machine learning classifier model can successfully discriminate the individual chains. We call this convergence measure $R^*$. In contrast to the predominant $\widehat{R}$, $R^*$ is a single statistic across all parameters that indicates lack of mixing, although individual variables' importance for this metric can also be determined. Additionally, $R^*$ is not based on any single characteristic of the sampling distribution; instead it uses all the information in the chain, including that given by the joint sampling distribution, which is currently largely overlooked by existing approaches. We recommend calculating $R^*$ using two different machine learning classifiers - gradient-boosted regression trees and random forests - which each work well in models of different dimensions. Because each of these methods outputs a classification probability, as a byproduct, we obtain uncertainty in $R^*$. The method is straightforward to implement and could be a complementary additional check on MCMC convergence for applied analyses.


title: Duals of semisimple Poisson-Lie groups and cluster theory of moduli spaces of G-local systems

abstract: We study the dual ${\rm G}^\ast$ of a standard semisimple Poisson-Lie group ${\rm G}$ from a perspective of cluster theory. We show that the coordinate ring $\mathcal{O}({\rm G}^\ast)$ can be naturally embedded into a cluster Poisson algebra with a Weyl group action. We prove that $\mathcal{O}({\rm G}^\ast)$ admits a natural basis which has positive integer structure coefficients and satisfies an invariance property with respect to a braid group action. We continue the study of the moduli space $\mathscr{P}_{{\rm G},\mathbb{S}}$ of ${\rm G}$-local systems introduced in \cite{GS3}, and prove that the coordinate ring of $\mathscr{P}_{{\rm G}, \mathbb{S}}$ coincides with its underlying cluster Poisson algebra.


title: Rotational excitations in rare-earth nuclei: a comparative study within three cranking models with different mean fields and the treatments of pairing correlations

abstract: High-spin rotational bands in rare-earth Er ($Z=68$), Tm ($Z=69$) and Yb ($Z=70$) isotopes are investigated by three different nuclear models. These are (i) the cranked relativistic Hartree-Bogoliubov (CRHB) approach with approximate particle number projection by means of the Lipkin-Nogami (LN) method, (ii) the cranking covariant density functional theory (CDFT) with pairing correlations treated by a shell-model-like approach (SLAP) or the so called particle-number conserving (PNC) method, and (iii) cranked shell model (CSM) based on the Nilsson potential with pairing correlations treated by the PNC method. A detailed comparison between these three models in the description of the ground state rotational bands of even-even Er and Yb isotopes is performed. The similarities and differences between these models in the description of the moments of inertia, the features of band crossings, equilibrium deformations and pairing energies of even-even nuclei under study are discussed. These quantities are considered as a function of rotational frequency and proton and neutron numbers. The changes in the properties of the first band crossings with increasing neutron number in this mass region are investigated. On average, a comparable accuracy of the description of available experimental data is achieved in these models. However, the differences between model predictions become larger above the first band crossings. Because of time-consuming nature of numerical calculations in the CDFT-based models, a systematic study of the rotational properties of both ground state and excited state bands in odd-mass Tm nuclei is carried out only by the PNC-SCM. With few exceptions, the rotational properties of experimental 1-quasiparticle and 3-quasiparticle bands in $^{165,167,169,171}$Tm are reproduced reasonably well.


title: Hardness of Bounded Distance Decoding on Lattices in $\ell_p$ Norms

abstract: $ \newcommand{\Z}{\mathbb{Z}} \newcommand{\eps}{\varepsilon} \newcommand{\cc}[1]{\mathsf{#1}} \newcommand{\NP}{\cc{NP}} \newcommand{\problem}[1]{\mathrm{#1}} \newcommand{\BDD}{\problem{BDD}} $Bounded Distance Decoding $\BDD_{p,\alpha}$ is the problem of decoding a lattice when the target point is promised to be within an $\alpha$ factor of the minimum distance of the lattice, in the $\ell_{p}$ norm. We prove that $\BDD_{p, \alpha}$ is $\NP$-hard under randomized reductions where $\alpha \to 1/2$ as $p \to \infty$ (and for $\alpha=1/2$ when $p=\infty$), thereby showing the hardness of decoding for distances approaching the unique-decoding radius for large $p$. We also show fine-grained hardness for $\BDD_{p,\alpha}$. For example, we prove that for all $p \in [1,\infty) \setminus 2\Z$ and constants $C > 1, \eps > 0$, there is no $2^{(1-\eps)n/C}$-time algorithm for $\BDD_{p,\alpha}$ for some constant $\alpha$ (which approaches $1/2$ as $p \to \infty$), assuming the randomized Strong Exponential Time Hypothesis (SETH). Moreover, essentially all of our results also hold (under analogous non-uniform assumptions) for $\BDD$ with preprocessing, in which unbounded precomputation can be applied to the lattice before the target is available. Compared to prior work on the hardness of $\BDD_{p,\alpha}$ by Liu, Lyubashevsky, and Micciancio (APPROX-RANDOM 2008), our results improve the values of $\alpha$ for which the problem is known to be $\NP$-hard for all $p > p_1 \approx 4.2773$, and give the very first fine-grained hardness for $\BDD$ (in any norm). Our reductions rely on a special family of "locally dense" lattices in $\ell_{p}$ norms, which we construct by modifying the integer-lattice sparsification technique of Aggarwal and Stephens-Davidowitz (STOC 2018).


title: Generating Electronic Health Records with Multiple Data Types and Constraints

abstract: Sharing electronic health records (EHRs) on a large scale may lead to privacy intrusions. Recent research has shown that risks may be mitigated by simulating EHRs through generative adversarial network (GAN) frameworks. Yet the methods developed to date are limited because they 1) focus on generating data of a single type (e.g., diagnosis codes), neglecting other data types (e.g., demographics, procedures or vital signs) and 2) do not represent constraints between features. In this paper, we introduce a method to simulate EHRs composed of multiple data types by 1) refining the GAN model, 2) accounting for feature constraints, and 3) incorporating key utility measures for such generation tasks. Our analysis with over $770,000$ EHRs from Vanderbilt University Medical Center demonstrates that the new model achieves higher performance in terms of retaining basic statistics, cross-feature correlations, latent structural properties, feature constraints and associated patterns from real data, without sacrificing privacy.


title: Self-Supervised Log Parsing

abstract: Logs are extensively used during the development and maintenance of software systems. They collect runtime events and allow tracking of code execution, which enables a variety of critical tasks such as troubleshooting and fault detection. However, large-scale software systems generate massive volumes of semi-structured log records, posing a major challenge for automated analysis. Parsing semi-structured records with free-form text log messages into structured templates is the first and crucial step that enables further analysis. Existing approaches rely on log-specific heuristics or manual rule extraction. These are often specialized in parsing certain log types, and thus, limit performance scores and generalization. We propose a novel parsing technique called NuLog that utilizes a self-supervised learning model and formulates the parsing task as masked language modeling (MLM). In the process of parsing, the model extracts summarizations from the logs in the form of a vector embedding. This allows the coupling of the MLM as pre-training with a downstream anomaly detection task. We evaluate the parsing performance of NuLog on 10 real-world log datasets and compare the results with 12 parsing techniques. The results show that NuLog outperforms existing methods in parsing accuracy with an average of 99% and achieves the lowest edit distance to the ground truth templates. Additionally, two case studies are conducted to demonstrate the ability of the approach for log-based anomaly detection in both supervised and unsupervised scenario. The results show that NuLog can be successfully used to support troubleshooting tasks. The implementation is available at https://github.com/nulog/nulog.


title: Vector Valued Transference

abstract: Our principal result is the following. Let $X$ and $Y$ be Banach spaces, let $G$ be a locally compact abelian group, and let $K$ be an operator valued kernel defined on $G$ with values in the space of bounded linear operators from $X$ to $Y$. Suppose that $R$ and $\tilde{R}$ are representations of $G$ on $X$ and $Y$ respectively that intertwine the values of $K$. Then, under suitable boundedness conditions on $R, \tilde{R}$ and $K$, the formula $$T_Kx = \int_GK(u)R_{-u}xdu $$ defines a bounded linear operator $T_K$ from $X$ to $Y$ with norm controlled by norm of convolution by $K$ as a mapping from $L^p_X(G)$ into $L^p_Y(G)$, (for all values of $p$ in the range $1\le p < \infty$.) A number of applications to the geometry of Banach spaces are given. Several results are proved in the setting of abstract commutative harmonic analysis. We outline the proof of the affirmative resolution of a conjecture of Rubio de Francia. This technique of transference is used to obtain dimension free estimates for certain operators in an $\Bbb{R}^n$ setting.


title: On the Use of Quantum Entanglement in Secure Communications: A Survey

abstract: Quantum computing and quantum communications are exciting new frontiers in computing and communications. Indeed, the massive investments made by the governments of the US, China, and EU in these new technologies are not a secret and are based on the expected potential of these technologies to revolutionize communications, computing, and security. In addition to several field trials and hero experiments, a number of companies such as Google and IBM are actively working in these areas and some have already reported impressive demonstrations in the past few years. While there is some skepticism about whether quantum cryptography will eventually replace classical cryptography, the advent of quantum computing could necessitate the use of quantum cryptography as the ultimate frontier of secure communications. This is because, with the amazing speeds demonstrated with quantum computers, breaking cryptographic keys might no longer be a daunting task in the next decade or so. Hence, quantum cryptography as the ultimate frontier in secure communications might not be such a far-fetched idea. It is well known that Heisenberg's Uncertainty Principle is essentially a "negative result" in Physics and Quantum Mechanics. It turns out that Heisenberg's Uncertainty Principle, one of the most interesting results in Quantum Mechanics, could be the theoretical basis and the main scientific principle behind the ultimate frontier in quantum cryptography or secure communications in conjunction with Quantum Entanglement.


title: Deep connections between learning from limited labels & physical parameter estimation -- inspiration for regularization

abstract: Recently established equivalences between differential equations and the structure of neural networks enabled some interpretation of training of a neural network as partial-differential-equation (PDE) constrained optimization. We add to the previously established connections, explicit regularization that is particularly beneficial in the case of single large-scale examples with partial annotation. We show that explicit regularization of model parameters in PDE constrained optimization translates to regularization of the network output. Examination of the structure of the corresponding Lagrangian and backpropagation algorithm do not reveal additional computational challenges. A hyperspectral imaging example shows that minimum prior information together with cross-validation for optimal regularization parameters boosts the segmentation accuracy.


title: Degenerate Affine Flag Varieties and Quiver Grassmannians

abstract: We study finite dimensional approximations to degenerate versions of affine flag varieties using quiver Grassmannians for cyclic quivers. We prove that they admit cellular decompositions parametrized by affine Dellac configurations, and that their irreducible components are normal Cohen-Macaulay varieties with rational singularities.


title: Search for the high-spin members of the $\alpha$:2n:$\alpha$ band in $^{10}$Be

abstract: Method: Measuring excitation functions for $^{6}$He+$\alpha$ scattering, populating states in the excitation energy range from 4.5 MeV to 8 MeV in $^{10}$Be using a $^6$He rare-isotope beam and a thick helium gas target. Results: No new excited states in $^{10}$Be have been observed. Stringent limitation on the possible degree of $\alpha$-clustering of the hypothetical yrast 6$^+$ state has been obtained. Conclusions: The high-spin members of the $\alpha$:2n:$\alpha$ molecular-like rotational band configuration, that is considered to have a 0$^+$ bandhead at 6.18 MeV, either do not exist or have small overlap with the $^{6}$He(g.s.)+$\alpha$ channel.


title: Breast Cancer Detection Using Convolutional Neural Networks

abstract: Breast cancer is prevalent in Ethiopia that accounts 34% among women cancer patients. The diagnosis technique in Ethiopia is manual which was proven to be tedious, subjective, and challenging. Deep learning techniques are revolutionizing the field of medical image analysis and hence in this study, we proposed Convolutional Neural Networks (CNNs) for breast mass detection so as to minimize the overheads of manual analysis. CNN architecture is designed for the feature extraction stage and adapted both the Region Proposal Network (RPN) and Region of Interest (ROI) portion of the faster R-CNN for the automated breast mass abnormality detection. Our model detects mass region and classifies them into benign or malignant abnormality in mammogram(MG) images at once. For the proposed model, MG images were collected from different hospitals, locally.The images were passed through different preprocessing stages such as gaussian filter, median filter, bilateral filters and extracted the region of the breast from the background of the MG image. The performance of the model on test dataset is found to be: detection accuracy 91.86%, sensitivity of 94.67% and AUC-ROC of 92.2%.


title: Predicting the ultimate outcome of the COVID-19 outbreak in Italy

abstract: During the COVID-19 outbreak, it is essential to monitor the effectiveness of measures taken by governments on the course of the epidemic. Here we show that there is already a sufficient amount of data collected in Italy to predict the outcome of the process. We show that using the proper metric, the data from Hubei Province and Italy has striking similarity, which enables us to calculate the expected number of confirmed cases and the number of deaths by the end of the process. Our predictions will improve as new data points are generated day by day, which can help to make further public decisions. The method is based on the data analysis of logistic growth equations describing the process on the macroscopic level. At the time of writing of the first version, the number of fatalities in Italy was expected to be 6000, and the predicted end of the crisis was April 15, 2020. In this new version, we discuss what changed in the two weeks which passed since then. The trend changed drastically on March 17, 2020, when the Italian health system reached its capacity limit. Without this limit, probably 3500 more people would have died. Instead, due to the limitations, 17.000 people are expected to die now, which is a five-fold increase. The predicted end of the crisis now shifted to May 8, 2020.


title: Regularization of linear and nonlinear ill-posed problems by mollification

abstract: In this paper, we address the problem of approximating solutions of ill-posed problems using mollification. We quickly review existing mollification regularization methods and provide two new approximate solutions to a general ill-posed equation $T(f) =g$ where $T$ can be nonlinear. The regularized solutions we define extend the work of Bonnefond and Mar\'echal \cite{xapi}, and trace their origins in the variational formulation of mollification, which to the best of our knowledge, was first introduced by Lannes et al. \cite{lannes}. In addition to consistency results, for the first time, we provide some convergence rates for a mollification method defined through a variational formulation.


title: Big Code != Big Vocabulary: Open-Vocabulary Models for Source Code

abstract: Statistical language modeling techniques have successfully been applied to large source code corpora, yielding a variety of new software development tools, such as tools for code suggestion, improving readability, and API migration. A major issue with these techniques is that code introduces new vocabulary at a far higher rate than natural language, as new identifier names proliferate. Both large vocabularies and out-of-vocabulary issues severely affect Neural Language Models (NLMs) of source code, degrading their performance and rendering them unable to scale. In this paper, we address this issue by: 1) studying how various modelling choices impact the resulting vocabulary on a large-scale corpus of 13,362 projects; 2) presenting an open vocabulary source code NLM that can scale to such a corpus, 100 times larger than in previous work; and 3) showing that such models outperform the state of the art on three distinct code corpora (Java, C, Python). To our knowledge, these are the largest NLMs for code that have been reported. All datasets, code, and trained models used in this work are publicly available.


title: The value of randomized strategies in distributionally robust risk averse network interdiction games

abstract: Conditional Value at Risk (CVaR) is widely used to account for the preferences of a risk-averse agent in the extreme loss scenarios. To study the effectiveness of randomization in interdiction games with an interdictor that is both risk and ambiguity averse, we introduce a distributionally robust network interdiction game where the interdictor randomizes over the feasible interdiction plans in order to minimize the worst-case CVaR of the flow with respect to both the unknown distribution of the capacity of the arcs and his mixed strategy over interdicted arcs. The flow player, on the contrary, maximizes the total flow in the network. By using the budgeted uncertainty set, we control the degree of conservatism in the model and reformulate the interdictor's non-linear problem as a bi-convex optimization problem. For solving this problem to any given optimality level, we devise a spatial branch and bound algorithm that uses the McCormick inequalities and reduced reformulation linearization technique (RRLT) to obtain convex relaxation of the problem. We also develop a column generation algorithm to identify the optimal support of the convex relaxation which is then used in the coordinate descent algorithm to determine the upper bounds. The efficiency and convergence of the spatial branch and bound algorithm is established in the numerical experiments. Further, our numerical experiments show that randomized strategies can have significantly better in-sample and out-of-sample performance than optimal deterministic ones.


title: Combinatory Chemistry: Towards a Simple Model of Emergent Evolution

abstract: An explanatory model for the emergence of evolvable units must display emerging structures that (1) preserve themselves in time (2) self-reproduce and (3) tolerate a certain amount of variation when reproducing. To tackle this challenge, here we introduce Combinatory Chemistry, an Algorithmic Artificial Chemistry based on a minimalistic computational paradigm named Combinatory Logic. The dynamics of this system comprise very few rules, it is initialised with an elementary tabula rasa state, and features conservation laws replicating natural resource constraints. Our experiments show that a single run of this dynamical system with no external intervention discovers a wide range of emergent patterns. All these structures rely on acquiring basic constituents from the environment and decomposing them in a process that is remarkably similar to biological metabolisms. These patterns include autopoietic structures that maintain their organisation, recursive ones that grow in linear chains or binary-branching trees, and most notably, patterns able to reproduce themselves, duplicating their number at each generation.


title: Fluctuations of a membrane nanotube covered with an actin sleeve

abstract: Many biological functions rely on the reshaping of cell membranes, in particular into nanotubes, which are covered in vivo by dynamic actin networks. Nanotubes are subject to thermal fluctuations, but the effect of these on cell functions is unknown. Here, we form nanotubes from liposomes using an optically trapped bead adhering to the liposome membrane. From the power spectral density of this bead, we study the nanotube fluctuations in the range of membrane tensions measured in vivo. We show that an actin sleeve covering the nanotube damps its high frequency fluctuations because of the network viscoelasticity. Our work paves the way for further studies on the effect of nanotube fluctuations in cellular functions.


title: The physical limits of nanoLEDs and nanolasers for optical communications

abstract: Nanoscale light sources are being intensively investigated for their potential to enable low-energy, high-density optical communication and sensing systems. Both nano-light-emitting diodes (nanoLEDs) and nanolasers have been considered, based on advanced nanophotonic concepts such as photonic crystals and plasmonic structures, with dimensions well into the sub-micrometer domain. With decreasing dimensions, light-matter interaction becomes stronger, potentially leading to efficient and ultrafast radiative emission, both in the spontaneous and stimulated regime. These features have created wide expectations for the practical prospects of such nanoscale light sources, in particular for optical interconnects. In this article we examine the limits to the downscaling of LEDs and lasers, and ask ourselves which type of source is most suited to ultralow-power optical communications. Based on simple physical considerations on the scaling of spontaneous and stimulated emission rates for semiconductor active regions at room temperature, we analyze the speed and energy limits for nanoLEDs and nanolasers as a function of their size. The role of spontaneous emission enhancement (Purcell effect) in practical nanophotonic sources is also revisited. The main conclusion is that nanoLEDs reach a fundamental energy/speed limit for data rates exceeding a few Gb/s, whereas nanolasers with active dimensions in the range of few 100s nm may enable direct modulation rates larger than 40 Gb/s at power levels adequate for short-distance and low-energy optical interconnects.


title: Diffusion in a biased washboard potential revisited

abstract: The celebrated Sutherland-Einstein relation for systems at thermal equilibrium states that spread of trajectories of Brownian particles is an increasing function of temperature. Here, we scrutinize diffusion of underdamped Brownian motion in a biased periodic potential and analyse regimes in which a diffusion coefficient decreases with increasing temperature within finite temperature window. Comprehensive numerical simulations of the corresponding Langevin equation performed with unprecedented resolution allow us to construct phase diagram for the occurrence of the non-monotonic temperature dependence of the diffusion coefficient. We discuss the relation of the latter effect with the phenomenon of giant diffusion.


title: QCD phase transition in the chiral limit

abstract: We present a lattice QCD based determination of the chiral phase transition temperature in QCD with two massless (up and down) and one strange quark having its physical mass. We propose and calculate two novel estimators for the chiral transition temperature for several values of the light quark masses, corresponding to Goldstone pion masses in the range of $58~{\rm MeV}\lesssim m_\pi\lesssim 163~{\rm MeV}$. The chiral phase transition temperature is determined by extrapolating to vanishing pion mass using universal scaling analysis. After thermodynamic, continuum and chiral extrapolations we find the chiral phase transition temperature $T_c^0=132^{+3}_{-6}$ MeV. We also present some preliminary calculations on interplay of effective $U_A(1)$ restoration and chiral phase transition towards chiral limit.


title: Semi-supervised Contrastive Learning Using Partial Label Information

abstract: In semi-supervised learning, information from unlabeled examples is used to improve the model learned from labeled examples. In some learning problems, partial label information can be inferred from otherwise unlabeled examples and used to further improve the model. In particular, partial label information exists when subsets of training examples are known to have the same label, even though the label itself is missing. By encouraging the model to give the same label to all such examples through contrastive learning objectives, we can potentially improve its performance. We call this encouragement Nullspace Tuning because the difference vector between any pair of examples with the same label should lie in the nullspace of a linear model. In this paper, we investigate the benefit of using partial label information using a careful comparison framework over well-characterized public datasets. We show that the additional information provided by partial labels reduces test error over good semi-supervised methods usually by a factor of 2, up to a factor of 5.5 in the best case. We also show that adding Nullspace Tuning to the newer and state-of-the-art MixMatch method decreases its test error by up to a factor of 1.8.


title: Accretion Properties of PDS 70b with MUSE

abstract: We report a new evaluation of the accretion properties of PDS~70b obtained with VLT/MUSE. The main difference from previous studies in Haffert et al. (2019) and Aoyama & Ikoma (2019) is in the mass accretion rate. Simultaneous multiple line observations, such as H$\alpha$ and H$\beta$, can better constrain the physical properties of an accreting planet. While we clearly detected H$\alpha$ emissions from PDS~70b, no H$\beta$ emissions were detected. We estimate the line flux of H$\beta$ with a 3-$\sigma$ upper limit to be 2.3~$\times$~10$^{-16}$~erg~s$^{-1}$~cm$^{-2}$. The flux ratio $F_{\rm H\beta}$/$F_{\rm H\alpha}$ for PDS~70b is $<$~0.28. Numerical investigations by Aoyama et al. (2018) suggest that $F_{\rm H\beta}$/$F_{\rm H\alpha}$ should be close to unity if the extinction is negligible. We attribute the reduction of the flux ratio to the extinction, and estimate the extinction of H$\alpha$ ($A_{\rm H\alpha}$) for PDS~70b to be $>$~2.0~mag using the interstellar extinction value. %The expected $A_{\rm H\alpha}$ value in the gap of the protoplanetary disk at the PDS~70b location is 2.4~mag, which is consistent with the estimated extinction. By combining with the H$\alpha$ linewidth and the dereddening line luminosity of H$\alpha$, %we derive the PDS~70b dynamical mass and mass accretion rate to be \hashimotor{12~$\pm$~3~$M_{\rm Jup}$} and $\gtrsim$~5~$\times$~10$^{-7}$~$M_{\rm Jup}$~yr$^{-1}$, respectively. we derive the PDS~70b mass accretion rate to be $\gtrsim$~5~$\times$~10$^{-7}$~$M_{\rm Jup}$~yr$^{-1}$. The PDS~70b mass accretion rate is an order of magnitude larger than that of PDS~70. We found that the filling factor $f_{\rm f}$ (the fractional area of the planetary surface emitting H$\alpha$) is $\gtrsim$0.01, which is similar to the typical stellar value. The small value of $f_{\rm f}$ indicates that the H$\alpha$ emitting areas are localized at the surface of PDS~70b.


title: 3D medical image segmentation with labeled and unlabeled data using autoencoders at the example of liver segmentation in CT images

abstract: Automatic segmentation of anatomical structures with convolutional neural networks (CNNs) constitutes a large portion of research in medical image analysis. The majority of CNN-based methods rely on an abundance of labeled data for proper training. Labeled medical data is often scarce, but unlabeled data is more widely available. This necessitates approaches that go beyond traditional supervised learning and leverage unlabeled data for segmentation tasks. This work investigates the potential of autoencoder-extracted features to improve segmentation with a CNN. Two strategies were considered. First, transfer learning where pretrained autoencoder features were used as initialization for the convolutional layers in the segmentation network. Second, multi-task learning where the tasks of segmentation and feature extraction, by means of input reconstruction, were learned and optimized simultaneously. A convolutional autoencoder was used to extract features from unlabeled data and a multi-scale, fully convolutional CNN was used to perform the target task of 3D liver segmentation in CT images. For both strategies, experiments were conducted with varying amounts of labeled and unlabeled training data. The proposed learning strategies improved results in $75\%$ of the experiments compared to training from scratch and increased the dice score by up to $0.040$ and $0.024$ for a ratio of unlabeled to labeled training data of about $32 : 1$ and $12.5 : 1$, respectively. The results indicate that both training strategies are more effective with a large ratio of unlabeled to labeled training data.


title: Sensor Selection With Cost Constraints for Dynamically Relevant Bases

abstract: We consider cost-constrained sparse sensor selection for full-state reconstruction, applying a well-known greedy algorithm to dynamical systems for which the usual singular value decomposition (SVD) basis may not be available or preferred. We apply the cost-modified, column-pivoted QR decomposition to a physically relevant basis -- the pivots correspond to sensor locations, and these locations are penalized with a heterogeneous cost function. In considering different bases, we are able to account for the dynamics of the particular system, yielding sensor arrays that are nearly Pareto optimal in sensor cost and performance in the chosen basis. This flexibility extends our framework to include actuation and dynamic estimation, and to select sensors without training data. We provide three examples from the physical and engineering sciences and evaluate sensor selection in three dynamically relevant bases: truncated balanced modes for control systems, dynamic mode decomposition (DMD) modes, and a basis of analytic modes. We find that these bases all yield effective sensor arrays and reconstructions for their respective systems. When possible, we compare to results using an SVD basis and evaluate tradeoffs between methods.


title: Measurement of the heat flux normalised spin Seebeck coefficient of thin films as a function of temperature

abstract: The spin Seebeck effect (SSE) has generated interest in the thermoelectric and magnetic communities for potential high efficiency energy harvesting applications, and spintronic communities as a source of pure spin current. To understand the underlying mechanisms requires characterisation of potential materials across a range of temperatures, however, for thin films the default measurement of an applied temperature gradient (across the sample) has been shown to be compromised by the presence of thermal resistances. Here, we demonstrate a method to perform low temperature SSE measurements where instead of monitoring the temperature gradient, the heat flux passing through the sample is measured using two calibrated heat flux sensors. This has the advantage of measuring the heat loss through the sample as well as providing a reliable method to normalise the SSE response of thin film samples. We demonstrate this method with an $\text{SiO}_{2}/\text{Fe}_{3}O_{4}/\text{Pt}$ sample, where a semiconducting-insulating transition occurs at the Verwey transition, $T_{\text{V}}$, of $\text{Fe}_{3}\text{O}_{4}$ and quantify the thermomagnetic response above and below $T_{\text{V}}$.


title: Improving predictions by nonlinear regression models from outlying input data

abstract: When applying machine learning/statistical methods to the environmental sciences, nonlinear regression (NLR) models often perform only slightly better and occasionally worse than linear regression (LR). The proposed reason for this conundrum is that NLR models can give predictions much worse than LR when given input data which lie outside the domain used in model training. Continuous unbounded variables are widely used in environmental sciences, whence not uncommon for new input data to lie far outside the training domain. For six environmental datasets, inputs in the test data were classified as "outliers" and "non-outliers" based on the Mahalanobis distance from the training input data. The prediction scores (mean absolute error, Spearman correlation) showed NLR to outperform LR for the non-outliers, but often underperform LR for the outliers. An approach based on Occam's Razor (OR) was proposed, where linear extrapolation was used instead of nonlinear extrapolation for the outliers. The linear extrapolation to the outlier domain was based on the NLR model within the non-outlier domain. This NLR$_{\mathrm{OR}}$ approach reduced occurrences of very poor extrapolation by NLR, and it tended to outperform NLR and LR for the outliers. In conclusion, input test data should be screened for outliers. For outliers, the unreliable NLR predictions can be replaced by NLR$_{\mathrm{OR}}$ or LR predictions, or by issuing a "no reliable prediction" warning.


title: Amplified Spontaneous Emission and Random Lasing in MAPbBr$_3$ Halide Perovskite Single Crystals

abstract: Halide perovskites are a promising optical gain media with high tunability and simple solution synthesis. In this study, two gain regimes, namely amplified spontaneous emission and random lasing, are demonstrated in same MAPbBr$_3$ halide perovskite single crystal. For this, photoluminescence is measured at a temperature of 4 K with pulsed femtosecond pumping by UV light with a 80 MHz repetition rate. Random lasing is observed in areas of the sample where a random resonator was formed due to cracks and crystal imperfections. In more homogeneous regions of the sample, the dominant regime is amplified spontaneous emission. These two regimes are reliably distinguished by the line width, the mode structure, the growth of the intensity after the threshold, and the degree of polarization of the radiation. The spectral localization of the stimulated emission well below the bound exciton resonance raises a question concerning the origin of the emission in halide perovskite lasers.


title: Epilocal: a real-time tool for local epidemic monitoring

abstract: We describe Epilocal, a simple R program designed to automatically download the most recent data on reported infected SARS-CoV-2 cases for all Italian provinces and regions, and to provide a simple descriptive analysis. For each province the cumulative number of reported infected cases is available each day. In addition, the current numbers of hospitalized patients (separately for intensive care or not) and the cumulative number of deceased individuals are available at the region level. The data are analyzed through Poisson generalized linear models with logarithmic link function and polynomial regression on time. For cumulative data, we also consider a logistic parameterisation of the hazard function. Automatic model selection is performed to choose among the different model specifications, based on the statistical significance of the corresponding estimated parameters and on goodness-of-fit assessment. The chosen model is used to produce up-to-today estimates of the growth rate of the counts. Results are plotted on a map of the country to allow for a visual assessment of the geographic distribution of the areas with differential prevalence and rates of growth.


title: The limits of sustained self-excitation and stable periodic pulse trains in the Yamada model with delayed optical feedback

abstract: We consider the Yamada model for an excitable or self-pulsating laser with saturable absorber, and study the effects of delayed optical self-feedback in the excitable case. More specifically, we are concerned with the generation of stable periodic pulse trains via repeated self-excitation after passage through the delayed feedback loop, as well as their bifurcations. We show that onset and termination of such pulse trains correspond to the simultaneous bifurcation of countably many fold periodic orbits with infinite period in this delay differential equation. We employ numerical continuation and the concept of reappearance of periodic solutions to show that these bifurcations coincide with codimension-two points along families of connecting orbits and fold periodic orbits in a related advanced differential equation. These points include heteroclinic connections between steady states, as well as homoclinic bifurcations with non-hyperbolic equilibria. Tracking these codimension-two points in parameter space reveals the critical parameter values for the existence of periodic pulse trains. We use the recently developed theory of temporal dissipative solitons to infer necessary conditions for the stability of such pulse trains.


title: Virtually RFRS Mapping Tori and Coherence

abstract: Let $G$ be a finitely presented group that can be written as an extension \[ 1 \longrightarrow K \longrightarrow G \longrightarrow F_2 \longrightarrow 1 \] where $K$ is either the finitely generated free group $F_n$, $n > 2$ or the fundamental group of a closed surface of genus $g > 1$. We prove that if the image of the monodromy map $\rho \colon F_2 \to \operatorname{Out(K)}$ contains an element $\varphi \in \operatorname{Out(K)}$ such that the mapping torus $K \rtimes_{\varphi} \Bbb{Z}$ is virtually residually finite rationally solvable (for instance whenever the mapping torus is hyperbolic), then $G$ is not coherent. This applies, in particular, when the image is a purely pseudo--Anosov free subgroups of the mapping class group.


title: Medical Isotope Production with the IsoDAR Cyclotron

abstract: Authors describe technical advances that will allow the IsoDAR cyclotron -- being developed for neutrino physics research -- to produce medical isotopes more efficiently than existing cyclotrons can.


title: Getting to 99% Accuracy in Interactive Segmentation

abstract: Interactive object cutout tools are the cornerstone of the image editing workflow. Recent deep-learning based interactive segmentation algorithms have made significant progress in handling complex images and rough binary selections can typically be obtained with just a few clicks. Yet, deep learning techniques tend to plateau once this rough selection has been reached. In this work, we interpret this plateau as the inability of current algorithms to sufficiently leverage each user interaction and also as the limitations of current training/testing datasets. We propose a novel interactive architecture and a novel training scheme that are both tailored to better exploit the user workflow. We also show that significant improvements can be further gained by introducing a synthetic training dataset that is specifically designed for complex object boundaries. Comprehensive experiments support our approach, and our network achieves state of the art performance.


title: IBL-Infinity Model of String Topology from Perturbative Chern-Simons Theory

abstract: Following Cieliebak, Fukaya, Latschev and Volkov, we construct an IBL-infinity chain model for equivariant string topology on cyclic Hochschild cochains of de Rham cohomology. We study its properties and perform explicit computations.


title: Segmentation of brain tumor on magnetic resonance imaging using a convolutional architecture

abstract: The brain is a complex organ controlling cognitive process and physical functions. Tumors in the brain are accelerated cell growths affecting the normal function and processes in the brain. MRI scans provides detailed images of the body being one of the most common tests to diagnose brain tumors. The process of segmentation of brain tumors from magnetic resonance imaging can provide a valuable guide for diagnosis, treatment planning and prediction of results. Here we consider the problem brain tumor segmentation using a Deep learning architecture for use in tumor segmentation. Although the proposed architecture is simple and computationally easy to train, it is capable of reaching $IoU$ levels of 0.95.


title: Random bipartite posets and extremal problems

abstract: Previously, Erd\H{o}s, Kierstead and Trotter investigated the dimension of random height~$2$ partially ordered sets. Their research was motivated primarily by two goals: (1)~analyzing the relative tightness of the F\"{u}redi-Kahn upper bounds on dimension in terms of maximum degree; and (2)~developing machinery for estimating the expected dimension of a random labeled poset on $n$ points. For these reasons, most of their effort was focused on the case $0<p\le 1/2$. While bounds were given for the range $1/2\le p <1$, the relative accuracy of the results in the original paper deteriorated as $p$ approaches~$1$. Motivated by two extremal problems involving conditions that force a poset to contain a large standard example, we were compelled to revisit this subject, but now with primary emphasis on the range $1/2\le p<1$. Our sharpened analysis shows that as $p$ approaches~$1$, the expected value of dimension increases and then decreases, answering in the negative a question posed in the original paper. Along the way, we apply inequalities of Talagrand and Janson, establish connections with latin rectangles and the Euler product function, and make progress on both extremal problems.


title: Boosting Unconstrained Face Recognition with Auxiliary Unlabeled Data

abstract: In recent years, significant progress has been made in face recognition, which can be partially attributed to the availability of large-scale labeled face datasets. However, since the faces in these datasets usually contain limited degree and types of variation, the resulting trained models generalize poorly to more realistic unconstrained face datasets. While collecting labeled faces with larger variations could be helpful, it is practically infeasible due to privacy and labor cost. In comparison, it is easier to acquire a large number of unlabeled faces from different domains, which could be used to regularize the learning of face representations. We present an approach to use such unlabeled faces to learn generalizable face representations, where we assume neither the access to identity labels nor domain labels for unlabeled images. Experimental results on unconstrained datasets show that a small amount of unlabeled data with sufficient diversity can (i) lead to an appreciable gain in recognition performance and (ii) outperform the supervised baseline when combined with less than half of the labeled data. Compared with the state-of-the-art face recognition methods, our method further improves their performance on challenging benchmarks, such as IJB-B, IJB-C and IJB-S.


title: Finite-time Identification of Stable Linear Systems: Optimality of the Least-Squares Estimator

abstract: We present a new finite-time analysis of the estimation error of the Ordinary Least Squares (OLS) estimator for stable linear time-invariant systems. We characterize the number of observed samples (the length of the observed trajectory) sufficient for the OLS estimator to be $(\varepsilon,\delta)$-PAC, i.e., to yield an estimation error less than $\varepsilon$ with probability at least $1-\delta$. We show that this number matches existing sample complexity lower bounds [1,2] up to universal multiplicative factors (independent of ($\varepsilon,\delta)$ and of the system). This paper hence establishes the optimality of the OLS estimator for stable systems, a result conjectured in [1]. Our analysis of the performance of the OLS estimator is simpler, sharper, and easier to interpret than existing analyses. It relies on new concentration results for the covariates matrix.

categories: math.ST cs.LG cs.SY eess.SY stat.ML stat.TH


title: Magnetic Properties and Electronic Structure of Magnetic Topological Insulator MnBi$_2$Se$_4$

abstract: The intrinsic magnetic topological insulators MnBi$_2$X$_4$ (X = Se, Te) are promising candidates in realizing various novel topological states related to symmetry breaking by magnetic order. Although much progress had been made in MnBi$_2$Te$_4$, the study of MnBi$_2$Se$_4$ has been lacking due to the difficulty of material synthesis of the desired trigonal phase. Here, we report the synthesis of multilayer trigonal MnBi$_2$Se$_4$ with alternating-layer molecular beam epitaxy. Atomic-resolution scanning transmission electron microscopy (STEM) and scanning tunneling microscopy (STM) identify a well-ordered multilayer van der Waals (vdW) crystal with septuple-layer base units in agreement with the trigonal structure. Systematic thickness-dependent magnetometry studies illustrate the layered antiferromagnetic ordering as predicted by theory. Angle-resolved photoemission spectroscopy (ARPES) reveals the gapless Dirac-like surface state of MnBi$_2$Se$_4$, which demonstrates that MnBi$_2$Se$_4$ is a topological insulator above the magnetic ordering temperature. These systematic studies show that MnBi$_2$Se$_4$ is a promising candidate for exploring the rich topological phases of layered antiferromagnetic topological insulators.


title: Neural Networks for Encoding Dynamic Security-Constrained Optimal Power Flow

abstract: This paper introduces a framework to capture previously intractable optimization constraints and transform them to a mixed-integer linear program, through the use of neural networks. We encode the feasible space of optimization problems characterized by both tractable and intractable constraints, e.g. differential equations, to a neural network. Leveraging an exact mixed-integer reformulation of neural networks, we solve mixed-integer linear programs that accurately approximate solutions to the originally intractable non-linear optimization problem. We apply our methods to the AC optimal power flow problem (AC-OPF), where directly including dynamic security constraints renders the AC-OPF intractable. Our proposed approach has the potential to be significantly more scalable than traditional approaches. We demonstrate our approach for power system operation considering N-1 security and small-signal stability, showing how it can efficiently obtain cost-optimal solutions which at the same time satisfy both static and dynamic security constraints.


title: Towards a Peer-to-Peer Energy Market: an Overview

abstract: This work focuses on the electric power market, comparing the status quo with the recent trend towards the increase in distributed self-generation capabilities by prosumers. Starting from the existing tension between the intrinsically hierarchical current structure of the electricity distribution network and the substantially distributed and self-organising nature of the self-generation, we explore the limitations imposed by the current conditions. Initially, we introduce a potential multi-layered architecture for a Peer-to-Peer (P2P) energy market, discussing the fundamental aspects of local production and local consumption as part of a microgrid. Secondly, we analyse the consequent changes for the different users' roles, also in connection with some incentive models connected with the decentralisation of the power production. To give a full picture to the reader, we also scrutinise relevant elements of energy trading, such as Smart Contract and grid stability. Thirdly, we present an example of a typical P2P settlement, showcasing the role of all the previously analysed aspects. To conclude, we performed a review of relevant activities in this domain, to showcase where existing projects are going and what are the most important themes covered. Being this a work in progress, many open questions are still on the table and will be addressed in the next stages of the research. Eventually, by providing a reference model as base for further discussions and improvements, we would like to engage ourselves in a dialog with the different users and the broad community, oriented towards a more fair and ecological-friendly solution for the electricity market of the future.


title: Dynamics of tritrophic interaction with volatile compounds in plants

abstract: In this paper we will consider a mathematical model that describes, the tritrophic interaction between plants, herbivores and their natural enemies, where volatiles organic compounds (VOCs) released by plants play an important role. We show positivity and boundedness of the system solutions, existence of positive equilibrium and its local stability, we analyse global stability of positive equilibrium via the geometrical approach of Li and Muldowney. We pay attention to parameters in order to discuss different types of bifurcations. Finally, we present some numerical simulations to justify our analytical results.


title: Efficient Statistical Model for Predicting Electromagnetic Wave Distribution in Coupled Enclosures

abstract: The Random Coupling Model (RCM) has been successfully applied to predicting the statistics of currents and voltages at ports in complex electromagnetic (EM) enclosures operating in the short wavelength limit. Recent studies have extended the RCM to systems of multi-mode aperture-coupled enclosures. However, as the size (as measured in wavelengths) of a coupling aperture grows, the coupling matrix used in the RCM increases as well, and the computation becomes more complex and time-consuming. A simple Power Balance Model (PWB) can provide fast predictions for the \textit{averaged} power density of waves inside electrically-large systems for a wide range of cavity and coupling scenarios. However, the important interference induced fluctuations of the wavefield retained in the RCM is absent in PWB. Here we aim to combine the best aspects of each model to create a hybrid treatment and study the EM fields in coupled enclosure systems. The proposed hybrid approach provides both mean and fluctuation information of the EM fields without the full computational complexity of coupled-cavity RCM. We compare the hybrid model predictions with experiments on linear cascades of over-moded cavities. We find good agreement over a set of different loss parameters and for different coupling strengths between cavities. The range of validity and applicability of the hybrid method are tested and discussed.


title: Many cliques with few edges and bounded maximum degree

abstract: Generalized Tur\'an problems have been a central topic of study in extremal combinatorics throughout the last few decades. One such problem, maximizing the number of cliques of a fixed order in a graph with fixed number of vertices and bounded maximum degree, was recently completely resolved by Chase. Kirsch and Radcliffe raised a natural variant of this problem where the number of edges is fixed instead of the number of vertices. In this paper, we determine the maximum number of cliques of a fixed order in a graph with fixed number of edges and bounded maximum degree, resolving a conjecture by Kirsch and Radcliffe. We also give a complete characterization of the extremal graphs.


title: Chandra Resolves the Double FU Orionis System RNO 1B/1C in X-rays

abstract: We present new Chandra X-ray observations of the close pair of young stars RNO 1B and 1C (6$''$ separation) located in the L1287 cloud. RNO 1B erupted in 1978 - 1990 and is classified as an FU Orionis star (FUor). RNO 1C also shows most of the properties of an FUor but no eruption has yet been seen. Only a few dozen FUors are known and the presence of two such objects with a small angular separation is rare, suggesting a common origin. Both stars were faintly detected by Chandra and we summarize their X-ray properties within the framework of other previously detected FUors. We also report other X-ray detections in L1287 including the deeply-embedded young star RNO 1G, the jet-like radio source VLA 3, and an enigmatic hard flaring source with no 2MASS counterpart that was only detected in the second of two Chandra exposures.


title: Co-Optimizing Performance and Memory FootprintVia Integrated CPU/GPU Memory Management, anImplementation on Autonomous Driving Platform

abstract: Cutting-edge embedded system applications, such as self-driving cars and unmanned drone software, are reliant on integrated CPU/GPU platforms for their DNNs-driven workload, such as perception and other highly parallel components. In this work, we set out to explore the hidden performance implication of GPU memory management methods of integrated CPU/GPU architecture. Through a series of experiments on micro-benchmarks and real-world workloads, we find that the performance under different memory management methods may vary according to application characteristics. Based on this observation, we develop a performance model that can predict system overhead for each memory management method based on application characteristics. Guided by the performance model, we further propose a runtime scheduler. By conducting per-task memory management policy switching and kernel overlapping, the scheduler can significantly relieve the system memory pressure and reduce the multitasking co-run response time. We have implemented and extensively evaluated our system prototype on the NVIDIA Jetson TX2, Drive PX2, and Xavier AGX platforms, using both Rodinia benchmark suite and two real-world case studies of drone software and autonomous driving software.


title: AC Measurement of the Nernst effect of thin films at low temperatures

abstract: We describe an alternating current method to measure the Nernst effect in superconducting thin films at low temperatures. The Nernst effect is an important tool in the understanding superconducting fluctuations and, in particular, vortex motion near critical points. However, in most materials, the Nernst signal in a typical experimental setup rarely exceeds a few $\mu$V, in some cases being as low as a few nV. DC measurements of such small signals require extensive signal processing and protection against stray pickups and offsets, limiting the sensitivity of such measurements to $>$5nV. Here we describe a method utilizing a one-heater-two-thermometer setup with the heating element and thermometers fabricated on-chip with the sample, which helped to reduce thermal load and temperature lag between the substrate and thermometer. Using AC heating power and 2$\omega$ measurement, we are able to achieve sub-nanovolt sensitivity in 20-30nm MoGe thin films on glass substrate, compared to a sensitivity of $\sim$9nV using DC techniques on the same setup.


title: Spontaneous Symmetry Breakdown of Scale Symmetries

abstract: We discuss spontaneous symmetry breakdown (SSB) of both global and local scale symmetries in scalar-tensor gravity with two scalar fields, one of which couples nonminimally to scalar curvature while the other is a normal scalar field. In case of a global scale symmetry, by moving from the Jordan frame to the Einstein frame, a normal scalar field becomes massive while the dilaton remains massless after the SSB. In case of a local scale symmetry, we take a gauge fixing condition for the local scale invariance, $a R + b \phi^2 = k$, which was found in our previous study of a Weyl's quadratic gravity. Together with locally scale invariant potential terms in a classical action, this gauge condition generates a Higgs potential whose vacuum expectation value (VEV) produces the Einstein-Hilbert action with an $R^2$ term in the lowest level of approximation. One interesting aspect in this model is that a massless Nambu-Goldstone boson associated with a scale invariance is absorbed into the metric tensor field and consequently an $R^2$ term is induced in the action.


title: BrazilDAM: A Benchmark dataset for Tailings Dam Detection

abstract: In this work we present BrazilDAM, a novel public dataset based on Sentinel-2 and Landsat-8 satellite images covering all tailings dams cataloged by the Brazilian National Mining Agency (ANM). The dataset was built using georeferenced images from 769 dams, recorded between 2016 and 2019. The time series were processed in order to produce cloud free images. The dams contain mining waste from different ore categories and have highly varying shapes, areas and volumes, making BrazilDAM particularly interesting and challenging to be used in machine learning benchmarks. The original catalog contains, besides the dam coordinates, information about: the main ore, constructive method, risk category, and associated potential damage. To evaluate BrazilDAM's predictive potential we performed classification essays using state-of-the-art deep Convolutional Neural Network (CNNs). In the experiments, we achieved an average classification accuracy of 94.11% in tailing dam binary classification task. In addition, others four setups of experiments were made using the complementary information from the original catalog, exhaustively exploiting the capacity of the proposed dataset.


title: Data-Driven Attack Detection for Linear Systems

abstract: This paper studies the attack detection problem in a data-driven and model-free setting, for deterministic systems with linear and time-invariant dynamics. Differently from existing studies that leverage knowledge of the system dynamics to derive security bounds and monitoring schemes, we focus on the case where the system dynamics, as well as the attack strategy and attack location, are unknown. We derive fundamental security limitations as a function of only the observed data and without estimating the system dynamics (in fact, no assumption is made on the identifiability of the system). In particular, (i) we derive detection limitations as a function of the informativity and length of the observed data, (ii) provide a data-driven characterization of undetectable attacks, and (iii) construct a data-driven detection monitor. Surprisingly, and in accordance with recent studies on data-driven control, our results show that model-based and data-driven security techniques share the same fundamental limitations, provided that the collected data remains sufficiently informative.


title: Big mapping class groups: an overview

abstract: We survey recent developments on mapping class groups of surfaces of infinite topological type.


title: Laser Spectroscopic Technique for Direct Identification of a Single Virus I: FASTER CARS

abstract: From the famous 1918 H1N1 influenza to the present COVID-19 pandemic, the need for improved virial detection techniques is all too apparent. The aim of the present paper is to show that identification of individual virus particles in clinical sample materials quickly and reliably is near at hand. First of all, our team has developed techniques for identification of virions based on a modular atomic force microscopy (AFM). Furthermore, Femtosecond Adaptive Spectroscopic Techniques with Enhanced Resolution via Coherent Anti-Stokes Raman Scattering (FASTER CARS) [1] using tip-enhanced techniques markedly improves the sensitivity.

categories: physics.med-ph physics.app-ph physics.bio-ph physics.optics


title: ParKCa: Causal Inference with Partially Known Causes

abstract: Methods for causal inference from observational data are an alternative for scenarios where collecting counterfactual data or realizing a randomized experiment is not possible. Adopting a stacking approach, our proposed method ParKCA combines the results of several causal inference methods to learn new causes in applications with some known causes and many potential causes. We validate ParKCA in two Genome-wide association studies, one real-world and one simulated dataset. Our results show that ParKCA can infer more causes than existing methods.


title: Nearest Neighbor Dirichlet Mixtures

abstract: There is a rich literature on Bayesian methods for density estimation, which characterize the unknown density as a mixture of kernels. Such methods have advantages in terms of providing uncertainty quantification in estimation, while being adaptive to a rich variety of densities. However, relative to frequentist locally adaptive kernel methods, Bayesian approaches can be slow and unstable to implement in relying on Markov chain Monte Carlo algorithms. To maintain most of the strengths of Bayesian approaches without the computational disadvantages, we propose a class of nearest neighbor-Dirichlet mixtures. The approach starts by grouping the data into neighborhoods based on standard algorithms. Within each neighborhood, the density is characterized via a Bayesian parametric model, such as a Gaussian with unknown parameters. Assigning a Dirichlet prior to the weights on these local kernels, we obtain a pseudo-posterior for the weights and kernel parameters. A simple and embarrassingly parallel Monte Carlo algorithm is proposed to sample from the resulting pseudo-posterior for the unknown density. Desirable asymptotic properties are shown, and the methods are evaluated in simulation studies and applied to a motivating data set in the context of classification.


title: Model reduction of linear hybrid systems

abstract: The paper proposes a model reduction algorithm for linear hybrid systems, i.e., hybrid systems with externally induced discrete events, with linear continuous subsystems, and linear reset maps. The model reduction algorithm is based on balanced truncation. Moreover, the paper also proves an analytical error bound for the difference between the input-output behaviors of the original and the reduced order model. This error bound is formulated in terms of singular values of the Gramians used for model reduction.


title: An End-to-end Framework For Low-Resolution Remote Sensing Semantic Segmentation

abstract: High-resolution images for remote sensing applications are often not affordable or accessible, especially when in need of a wide temporal span of recordings. Given the easy access to low-resolution (LR) images from satellites, many remote sensing works rely on this type of data. The problem is that LR images are not appropriate for semantic segmentation, due to the need for high-quality data for accurate pixel prediction for this task. In this paper, we propose an end-to-end framework that unites a super-resolution and a semantic segmentation module in order to produce accurate thematic maps from LR inputs. It allows the semantic segmentation network to conduct the reconstruction process, modifying the input image with helpful textures. We evaluate the framework with three remote sensing datasets. The results show that the framework is capable of achieving a semantic segmentation performance close to native high-resolution data, while also surpassing the performance of a network trained with LR inputs.


title: Limits on $w$CDM from the EFTofLSS with the PyBird code

abstract: We apply the Effective Field Theory of Large-Scale Structure to analyze the $w$CDM cosmological model. By using the full shape of the power spectrum and the BAO post-reconstruction measurements from BOSS, the Supernovae from Pantheon, and a prior from BBN, we set the competitive CMB-independent limit $w=-1.046_{-0.052}^{+0.055}$ at $68\%$ C.L.. After adding the Planck CMB data, we find $w=-1.023_{-0.030}^{+0.033}$ at $68\%$ C.L.. Our results are obtained using PyBird, a new, fast Python-based code which we make publicly available.


title: An Algorithm for Computing a Minimal Comprehensive Gr\"obner\, Basis of a Parametric Polynomial System

abstract: An algorithm to generate a minimal comprehensive Gr\"obner\, basis of a parametric polynomial system from an arbitrary faithful comprehensive Gr\"obner\, system is presented. A basis of a parametric polynomial ideal is a comprehensive Gr\"obner\, basis if and only if for every specialization of parameters in a given field, the specialization of the basis is a Gr\"obner\, basis of the associated specialized polynomial ideal. The key idea used in ensuring minimality is that of a polynomial being essential with respect to a comprehensive Gr\"obner\, basis. The essentiality check is performed by determining whether a polynomial can be covered for various specializations by other polynomials in the associated branches in a comprehensive Gr\"obner\, system. The algorithm has been implemented and successfully tried on many examples from the literature.


title: Minimal Kinematics: An All $k$ and $n$ Peek into ${\rm Trop}^+{\rm G}(k,n)$

abstract: In this note we present a formula for the Cachazo-Early-Guevara-Mizera (CEGM) generalized biadjoint amplitudes for all $k$ and $n$ on what we call the minimal kinematics. We prove that on the minimal kinematics, the scattering equations on the configuration space of $n$ points on $\mathbb{CP}^{k-1}$ has a unique solution, and that this solution is in the image of a Veronese embedding. The minimal kinematics is an all $k$ generalization of the one recently introduced by Early for $k=2$ and uses a choice of cyclic ordering. We conjecture an explicit formula for $m_n^{(k)}(\mathbb{I},\mathbb{I})$ which we have checked analytically through $n=10$ for all $k$. The answer is a simple rational function which has only simple poles; the poles have the combinatorial structure of the circulant graph ${\rm C}_n^{(1,2,\dots, k-2)}$. Generalized biadjoint amplitudes can also be evaluated using the positive tropical Grassmannian ${\rm Tr}^+{\rm G}(k,n)$ in terms of generalized planar Feynman diagrams. We find perfect agreement between both definitions for all cases where the latter is known in the literature. In particular, this gives the first strong consistency check on the $90\,608$ planar arrays for ${\rm Tr}^+{\rm G}(4,8)$ recently computed by Cachazo, Guevara, Umbert and Zhang. We also introduce another class of special kinematics called planar-basis kinematics which generalizes the one introduced by Cachazo, He and Yuan for $k=2$ and uses the planar basis recently introduced by Early for all $k$. Based on numerical computations through $n=8$ for all $k$, we conjecture that on the planar-basis kinematics $m_n^{(k)}(\mathbb{I},\mathbb{I})$ evaluates to the multidimensional Catalan numbers, suggesting the possibility of novel combinatorial interpretations. For $k=2$ these are the standard Catalan numbers.


title: Learning Nonlinear Loop Invariants with Gated Continuous Logic Networks (Extended Version)

abstract: Verifying real-world programs often requires inferring loop invariants with nonlinear constraints. This is especially true in programs that perform many numerical operations, such as control systems for avionics or industrial plants. Recently, data-driven methods for loop invariant inference have shown promise, especially on linear invariants. However, applying data-driven inference to nonlinear loop invariants is challenging due to the large numbers of and magnitudes of high-order terms, the potential for overfitting on a small number of samples, and the large space of possible inequality bounds. In this paper, we introduce a new neural architecture for general SMT learning, the Gated Continuous Logic Network (G-CLN), and apply it to nonlinear loop invariant learning. G-CLNs extend the Continuous Logic Network (CLN) architecture with gating units and dropout, which allow the model to robustly learn general invariants over large numbers of terms. To address overfitting that arises from finite program sampling, we introduce fractional sampling---a sound relaxation of loop semantics to continuous functions that facilitates unbounded sampling on real domain. We additionally design a new CLN activation function, the Piecewise Biased Quadratic Unit (PBQU), for naturally learning tight inequality bounds. We incorporate these methods into a nonlinear loop invariant inference system that can learn general nonlinear loop invariants. We evaluate our system on a benchmark of nonlinear loop invariants and show it solves 26 out of 27 problems, 3 more than prior work, with an average runtime of 53.3 seconds. We further demonstrate the generic learning ability of G-CLNs by solving all 124 problems in the linear Code2Inv benchmark. We also perform a quantitative stability evaluation and show G-CLNs have a convergence rate of $97.5\%$ on quadratic problems, a $39.2\%$ improvement over CLN models.


title: Remarks on the range and multiple range of random walk up to the time of exit

abstract: We consider the scaling behavior of the range and $p$-multiple range, that is the number of points visited and the number of points visited exactly $p\geq 1$ times, of simple random walk on ${\mathbb Z}^d$, for dimensions $d\geq 2$, up to time of exit from a domain $D_N$ of the form $D_N = ND$ where $D\subset {\mathbb R}^d$, as $N\uparrow\infty$. Recent papers have discussed connections of the range and related statistics with the Gaussian free field, identifying in particular that the distributional scaling limit for the range, in the case $D$ is a cube in $d\geq 3$, is proportional to the exit time of Brownian motion. The purpose of this note is to give a concise, different argument that the scaled range and multiple range, in a general setting in $d\geq 2$, both weakly converge to proportional exit times of Brownian motion from $D$, and that the corresponding limit moments are `polyharmonic', solving a hierarchy of Poisson equations.


title: An Exploratory Study of Bot Commits

abstract: Background: Bots help automate many of the tasks performed by software developers and are widely used to commit code in various social coding platforms. At present, it is not clear what types of activities these bots perform and understanding it may help design better bots, and find application areas which might benefit from bot adoption. Aim: We aim to categorize the Bot Commits by the type of change (files added, deleted, or modified), find the more commonly changed file types, and identify the groups of file types that tend to get updated together. Method: 12,326,137 commits made by 461 popular bots (that made at least 1000 commits) were examined to identify the frequency and the type of files added/ deleted/ modified by the commits, and association rule mining was used to identify the types of files modified together. Result: Majority of the bot commits modify an existing file, a few of them add new files, while deletion of a file is very rare. Commits involving more than one type of operation are even rarer. Files containing data, configuration, and documentation are most frequently updated, while HTML is the most common type in terms of the number of files added, deleted, and modified. Files of the type "Markdown", "Ignore List", "YAML", "JSON" were the types that are updated together with other types of files most frequently. Conclusion: We observe that majority of bot commits involve single file modifications, and bots primarily work with data, configuration, and documentation files. A better understanding if this is a limitation of the bots and, if overcome, would lead to different kinds of bots remains an open question.


title: Deliberation Model Based Two-Pass End-to-End Speech Recognition

abstract: End-to-end (E2E) models have made rapid progress in automatic speech recognition (ASR) and perform competitively relative to conventional models. To further improve the quality, a two-pass model has been proposed to rescore streamed hypotheses using the non-streaming Listen, Attend and Spell (LAS) model while maintaining a reasonable latency. The model attends to acoustics to rescore hypotheses, as opposed to a class of neural correction models that use only first-pass text hypotheses. In this work, we propose to attend to both acoustics and first-pass hypotheses using a deliberation network. A bidirectional encoder is used to extract context information from first-pass hypotheses. The proposed deliberation model achieves 12% relative WER reduction compared to LAS rescoring in Google Voice Search (VS) tasks, and 23% reduction on a proper noun test set. Compared to a large conventional model, our best model performs 21% relatively better for VS. In terms of computational complexity, the deliberation decoder has a larger size than the LAS decoder, and hence requires more computations in second-pass decoding.


title: A proof of Fatou's interpolation theorem

abstract: For Fatou's interpolation theorem of 1906 we suggest a new elementary proof.


title: Applications of Alpha Space

abstract: We extend the definition of alpha space as introduced in [1] to two spacetime dimensions. We discuss how this can be used to find conformal block decompositions of known functions and how to easily recover several lightcone bootstrap results. In the second part of the paper we establish a connection between alpha space and the Lorentzian inversion formula of [2].


title: Dynamic Information Design: A Simple Problem on Optimal Sequential Information Disclosure

abstract: We study a dynamic information design problem in a finite-horizon setting consisting of two strategic and long-term optimizing agents, namely a principal (he) and a detector (she). The principal observes the evolution of a Markov chain that has two states, one "good" and one "bad" absorbing state, and has to decide how to sequentially disclose information to the detector. The detector's only information consists of the messages she receives from the principal. The detector's objective is to detect as accurately as possible the time of the jump from the good to the bad state. The principal's objective is to delay the detector as much as possible from detective the jump to the bad state. For this setting, we determine the optimal strategies of the principal and the detector. The detector's optimal strategy is described by time-varying thresholds on her posterior belief of the good state. We prove that it is optimal for the principal to give no information to the detector before a time threshold, run a mixed strategy to confuse the detector at the threshold time, and reveal the true state afterwards. We present an algorithm that determines both the optimal time threshold and the optimal mixed strategy that could be employed by the principal. We show, through numerical experiments, that this optimal sequential mechanism significantly outperforms any other information disclosure strategy presented in literature.


title: Maximizing Influence-based Group Shapley Centrality

abstract: One key problem in network analysis is the so-called influence maximization problem, which consists in finding a set $S$ of at most $k$ seed users, in a social network, maximizing the spread of information from $S$. This paper studies a related but slightly different problem: We want to find a set $S$ of at most $k$ seed users that maximizes the spread of information, when $S$ is added to an already pre-existing - but unknown - set of seed users $T$. We consider such scenario to be very realistic. Assume a central entity wants to spread a piece of news, while having a budget to influence $k$ users. This central authority may know that some users are already aware of the information and are going to spread it anyhow. The identity of these users being however completely unknown. We model this optimization problem using the Group Shapley value, a well-founded concept from cooperative game theory. While the standard influence maximization problem is easy to approximate within a factor $1-1/e-\epsilon$ for any $\epsilon>0$, assuming common computational complexity conjectures, we obtain strong hardness of approximation results for the problem at hand in this paper. Maybe most prominently, we show that it cannot be approximated within $1/n^{o(1)}$ under the Gap Exponential Time Hypothesis. Hence, it is unlikely to achieve anything better than a polynomial factor approximation. Nevertheless, we show that a greedy algorithm can achieve a factor of $\frac{1-1/e}{k}-\epsilon$ for any $\epsilon>0$, showing that not all is lost in settings where $k$ is bounded.


title: Pricing with Variance Gamma Information

abstract: In the information-based pricing framework of Brody, Hughston and Macrina, the market filtration $\{ \mathcal F_t\}_{t\geq 0}$ is generated by an information process $\{ \xi_t\}_{t\geq0}$ defined in such a way that at some fixed time $T$ an $\mathcal F_T$-measurable random variable $X_T$ is "revealed". A cash flow $H_T$ is taken to depend on the market factor $X_T$, and one considers the valuation of a financial asset that delivers $H_T$ at $T$. The value $S_t$ of the asset at any time $t\in[0,T)$ is the discounted conditional expectation of $H_T$ with respect to $\mathcal F_t$, where the expectation is under the risk neutral measure and the interest rate is constant. Then $S_{T^-} = H_T$, and $S_t = 0$ for $t\geq T$. In the general situation one has a countable number of cash flows, and each cash flow can depend on a vector of market factors, each associated with an information process. In the present work, we construct a new class of models for the market filtration based on the variance-gamma process. The information process is obtained by subordinating a particular type of Brownian random bridge with a gamma process. The filtration is taken to be generated by the information process together with the gamma bridge associated with the gamma subordinator. We show that the resulting extended information process has the Markov property and hence can be used to price a variety of different financial assets, several examples of which are discussed in detail.


title: Force approach for the pseudopotential lattice Boltzmann method

abstract: The pseudopotential method is one of the most popular extensions of the lattice Boltzmann method (LBM) for phase change and multiphase flow simulation. One attractive feature of the original proposed method consists on its simplicity of adding a force dependent on a nearest-neighbor potential function, which became known as the Shan-Chen interaction force. Some of the well known drawbacks implied by this method involves lack of thermodynamic consistency and impossibility to control the surface tension independently. In order to correct these deficiencies, different approaches were developed in the literature, such as multirange interactions potential, which involves larger stencils than nearest-neighbor approach, and modified forcing schemes. In this work, a strategy is developed to control the liquid-gas density ratio and the surface tension by means of an appropriate interaction force field using only nearest-neighbor interactions. The proposed procedure is devised starting from the desired pressure tensor, which allow for the control of the equilibrium multiphase properties such as liquid-gas coexistence curve and surface tension. Then, it is shown how to derive an external force field able to replicate the effects of this pressure tensor in the macroscopic conservation equations. The final step of our procedure is implementing this external force in the LBE by using the classical Guo forcing scheme. Numerical tests regarding static and dynamic flow conditions were performed. Results obtained from simulations showed good agreement with expected analytical values. Most divergent solution observed was the droplet oscillation period under certain flow conditions, which deviated 9% from expected analytical result. The observed results corroborate that the proposed method is able to replicate the desired macroscopic multiphase behaviour.


title: Ford Multi-AV Seasonal Dataset

abstract: This paper presents a challenging multi-agent seasonal dataset collected by a fleet of Ford autonomous vehicles at different days and times during 2017-18. The vehicles traversed an average route of 66 km in Michigan that included a mix of driving scenarios such as the Detroit Airport, freeways, city-centers, university campus and suburban neighbourhoods, etc. Each vehicle used in this data collection is a Ford Fusion outfitted with an Applanix POS-LV GNSS system, four HDL-32E Velodyne 3D-lidar scanners, 6 Point Grey 1.3 MP Cameras arranged on the rooftop for 360-degree coverage and 1 Pointgrey 5 MP camera mounted behind the windshield for the forward field of view. We present the seasonal variation in weather, lighting, construction and traffic conditions experienced in dynamic urban environments. This dataset can help design robust algorithms for autonomous vehicles and multi-agent systems. Each log in the dataset is time-stamped and contains raw data from all the sensors, calibration values, pose trajectory, ground truth pose, and 3D maps. All data is available in Rosbag format that can be visualized, modified and applied using the open-source Robot Operating System (ROS). We also provide the output of state-of-the-art reflectivity-based localization for bench-marking purposes. The dataset can be freely downloaded at our website.


title: An Exploratory Study to Find Motives Behind Cross-platform Forks from Software Heritage Dataset

abstract: The fork-based development mechanism provides the flexibility and the unified processes for software teams to collaborate easily in a distributed setting without too much coordination overhead.Currently, multiple social coding platforms support fork-based development, such as GitHub, GitLab, and Bitbucket. Although these different platforms virtually share the same features, they have different emphasis. As GitHub is the most popular platform and the corresponding data is publicly available, most of the current studies are focusing on GitHub hosted projects. However, we observed anecdote evidences that people are confused about choosing among these platforms, and some projects are migrating from one platform to another, and the reasons behind these activities remain unknown.With the advances of Software Heritage Graph Dataset (SWHGD),we have the opportunity to investigate the forking activities across platforms. In this paper, we conduct an exploratory study on 10popular open-source projects to identify cross-platform forks and investigate the motivation behind. Preliminary result shows that cross-platform forks do exist. For the 10 subject systems in this study, we found 81,357 forks in total among which 179 forks are on GitLab. Based on our qualitative analysis, we found that most of the cross-platform forks that we identified are mirrors of the repositories on another platform, but we still find cases that were created due to preference of using certain functionalities (e.g. Continuous Integration (CI)) supported by different platforms. This study lays the foundation of future research directions, such as understanding the differences between platforms and supporting cross-platform collaboration.


title: The Iterative Transformation Method

abstract: In a transformation method, the numerical solution of a given boundary value problem is obtained by solving one or more related initial value problems. Therefore, a transformation method, like a shooting method, is an initial value method. The main difference between a transformation and a shooting method is that the former is conceived and derive its formulation from the scaling invariance theory. This paper is concerned with the application of the iterative transformation method to several problems in the boundary layer theory. The iterative method is an extension of the T{\"o}pfer's non-iterative algorithm developed as a simple way to solve the celebrated Blasius problem. This iterative method provides a simple numerical test for the existence and uniqueness of solutions. Here we show how the method can be applied to problems with a homogeneous boundary conditions at infinity and in particular we solve the Sakiadis problem of boundary layer theory. Moreover, we show how to couple our method with Newton's root-finder. The obtained numerical results compare well with those available in the literature. The main aim here is that any method developed for the Blasius, or the Sakiadis, problem might be extended to more challenging or interesting problems. In this context, the iterative transformation method has been recently applied to compute the normal and reverse flow solutions of Stewartson for the Falkner-Skan model [Comput. \& Fluids, {\bf 73} (2013) pp. 202-209].


title: State of Charge Estimation of Parallel Connected Battery Cells via Descriptor System Theory

abstract: This manuscript presents an algorithm for individual Lithium-ion (Li-ion) battery cell state of charge (SOC) estimation when multiple cells are connected in parallel, using only terminal voltage and total current measurements. For battery packs consisting of thousands of cells, it is desirable to estimate individual SOCs by only monitoring the total current in order to reduce sensing cost. Mathematically, series connected cells yield dynamics given by ordinary differential equations under classical full voltage sensing. In contrast, parallel connected cells are evidently more challenging because the dynamics are governed by a nonlinear descriptor system, including differential equations and algebraic equations arising from voltage and current balance across cells. An observer with linear output error injection is formulated, where the individual cell SOCs and local currents are locally observable from the total current and voltage measurements. The asymptotic convergence of differential and algebraic states is established by considering local Lipschitz continuity property of system nonlinearities. Simulation results on LiNiMnCoO$_2$/Graphite (NMC) cells illustrate convergence for SOCs, local currents, and terminal voltage.


title: Optimal isotropic, reusable truss lattice material with near-zero Poisson's ratio

abstract: Cork is a natural amorphous material with near-zero Poisson's ratio that is ubiquitously used for sealing glass bottles. It is an anisotropic, transversally isotropic, composite that can hardly be scaled down. Here, we propose a new class of isotropic and reusable cork-like metamaterial that is designed from an hybrid truss-lattice material to show an isotropic Poisson's ratio close to zero. Optimization is conducted using a multi-objective genetic algorithm, assisted by an elliptical basis function neural network, and coupled with finite element simulations. The optimal micro-structured metamaterial, fabricated by two-photon lithography with a lattice constant of 300 \micro\meter, has an almost isotropic Poisson's ratio smaller than 0.08 in all directions. It can recover $96.6\%$ of its original shape after a compressional test exceeding $20\%$ strain.


title: Witnessing non-classicality beyond quantum theory

abstract: We propose a general argument to show that if a physical system can mediate locally the generation of entanglement between two quantum systems, then it itself must be non-classical. Remarkably, we do not assume any classical or quantum formalism to describe the mediating physical system: our result follows from general information-theoretic principles, drawn from the recently proposed constructor theory of information. This argument provides the indispensable theoretical basis for recently proposed tests of non-classicality in gravity, based on witnessing gravitationally-induced entanglement in quantum probes.


title: Tight Decomposition Functions for Continuous-Time Mixed-Monotone Systems with Disturbances

abstract: The vector field of a mixed-monotone system is decomposable via a decomposition function into increasing (cooperative) and decreasing (competitive) components, and this decomposition allows for, e.g., efficient computation of reachable sets and forward invariant sets. A main challenge in this approach, however, is identifying an appropriate decomposition function. In this work, we show that any continuous-time dynamical system with a Lipschitz continuous vector field is mixed-monotone, and we provide a construction for the decomposition function that yields the tightest approximation of reachable sets when used with the standard tools for mixed-monotone systems. Our construction is similar to that recently proposed by Yang and Ozay for computing decomposition functions of discrete-time systems [1] where we make appropriate modifications for the continuous-time setting and also extend to the case with unknown disturbance inputs. As in [1], our decomposition function construction requires solving an optimization problem for each point in the state-space; however, we demonstrate through example how tight decomposition functions can sometimes be calculated in closed form. As a second contribution, we show how under-approximations of reachable sets can be efficiently computed via the mixed-monotonicity property by considering the backward-time dynamics.


title: Quantum mechanics is *-algebras and tensor networks

abstract: We provide a systematic approach to quantum mechanics from an information-theoretic perspective using the language of tensor networks. Our formulation needs only a single kind of object, so-called positive *-tensors. Physical models translate experimental setups into networks of these *-tensors, and the evaluation of the resulting networks yields the probability distributions describing measurement outcomes. The idea behind our approach is similar to categorical formulations of quantum mechanics. However, our formulation is mathematically simpler and less abstract. Our presentation of the core formalism is completely self-contained and relies on minimal mathematical prerequesites. Therefore, we hope it is in principle also understandable to people without an extensive mathematical background. Additionally, we show how various types of models, like real-time evolutions or thermal systems can be translated into *-tensor networks.


title: Assessing Robustness to Noise: Low-Cost Head CT Triage

abstract: Automated medical image classification with convolutional neural networks (CNNs) has great potential to impact healthcare, particularly in resource-constrained healthcare systems where fewer trained radiologists are available. However, little is known about how well a trained CNN can perform on images with the increased noise levels, different acquisition protocols, or additional artifacts that may arise when using low-cost scanners, which can be underrepresented in datasets collected from well-funded hospitals. In this work, we investigate how a model trained to triage head computed tomography (CT) scans performs on images acquired with reduced x-ray tube current, fewer projections per gantry rotation, and limited angle scans. These changes can reduce the cost of the scanner and demands on electrical power but come at the expense of increased image noise and artifacts. We first develop a model to triage head CTs and report an area under the receiver operating characteristic curve (AUROC) of 0.77. We then show that the trained model is robust to reduced tube current and fewer projections, with the AUROC dropping only 0.65% for images acquired with a 16x reduction in tube current and 0.22% for images acquired with 8x fewer projections. Finally, for significantly degraded images acquired by a limited angle scan, we show that a model trained specifically to classify such images can overcome the technological limitations to reconstruction and maintain an AUROC within 0.09% of the original model.


title: Localization Efficiency in Massive MIMO Systems

abstract: In the next generation of wireless systems, Massive MIMO offers high angular resolution for localization. By virtue of large number of antennas, the Angle of Arrival (AoA) of User Terminals (UTs) can be estimated with high accuracy. According to Dense Multipath Component (DMC) channel model, local scatters around UTs can create different multipath signals for each antenna at the Base Station (BS). We obtain a deterministic form for the Cramer-Rao Lower Bound (CRLB) in a multi-user scenario when the contribution of the multipath signals is considered. We do this when the multipath signals are independent and identically distributed (i.i.d) with arbitrary distribution. Then, we redefine a localization efficiency function for a multi-user scenario and numerically optimize it with respect to (w.r.t) the number of antennas. We prove when only a subset of the available antennas is used, CRLB can be minimized w.r.t which set of antennas is used. Then, an antenna selection strategy that minimizes CRLB is proposed. As a benchmark, we apply the proposed antenna selection scheme to the MUltiple SIgnal Classification (MUSIC) algorithm and study its efficiency. Numerical results validate the accuracy of our analysis and show significant improvement in efficiency when the proposed antenna selection strategy is employed.


title: Distributed Small-Signal Stability Conditions for Inverter-Based Unbalanced Microgrids

abstract: The proliferation of inverter-based generation and advanced sensing, controls, and communication infrastructure have facilitated the accelerated deployment of microgrids. A coordinated network of microgrids can maintain reliable power delivery to critical facilities during extreme events. Low inertia offered by the power electronics interfaced energy resources however, can present significant challenges to ensuring stable operation of the microgrids. In this work, distributed small-signal stability conditions for inverter-based microgrids are developed that involve the droop controller parameters and the network parameters such as line impedances, loads, etc. The distributed closed-form parametric stability conditions derived in this paper can be verified in a computationally efficient manner, facilitating the reliable design and operations of networks of microgrids. Dynamic phasor models have been used to capture the effects of electromagnetic transients. Numerical results are presented, along with PSCAD simulations, to validate the analytical stability conditions. Effects of design choices, such as the conductor types, and inverter sizes, on the small-signal stability of inverter-based microgrids are investigated to identify interpretable stable or unstable region estimates.


title: Mixing Rates for Hamiltonian Monte Carlo Algorithms in Finite and Infinite Dimensions

abstract: We establish the geometric ergodicity of the preconditioned Hamiltonian Monte Carlo (HMC) algorithm defined on an infinite-dimensional Hilbert space, as developed in [Beskos et al., Stochastic Process. Appl., 2011]. This algorithm can be used as a basis to sample from certain classes of target measures which are absolutely continuous with respect to a Gaussian measure. Our work addresses an open question posed in [Beskos et al., Stochastic Process. Appl., 2011], and provides an alternative to a recent proof based on exact coupling techniques given in arXiv:1909.07962. The approach here establishes convergence in a suitable Wasserstein distance by using the weak Harris theorem together with a generalized coupling argument. We also show that a law of large numbers and central limit theorem can be derived as a consequence of our main convergence result. Moreover, our approach yields a novel proof of mixing rates for the classical finite-dimensional HMC algorithm. As such, the methodology we develop provides a flexible framework to tackle the rigorous convergence of other Markov Chain Monte Carlo algorithms. Additionally, we show that the scope of our result includes certain measures that arise in the Bayesian approach to inverse PDE problems, cf. [Stuart, Acta Numer., 2010]. Particularly, we verify all of the required assumptions for a certain class of inverse problems involving the recovery of a divergence free vector field from a passive scalar, arXiv:1808.01084v3.


title: Segmentation and Optimal Region Selection of Physiological Signals using Deep Neural Networks and Combinatorial Optimization

abstract: Physiological signals, such as the electrocardiogram and the phonocardiogram are very often corrupted by noisy sources. Usually, artificial intelligent algorithms analyze the signal regardless of its quality. On the other hand, physicians use a completely orthogonal strategy. They do not assess the entire recording, instead they search for a segment where the fundamental and abnormal waves are easily detected, and only then a prognostic is attempted. Inspired by this fact, a new algorithm that automatically selects an optimal segment for a post-processing stage, according to a criteria defined by the user is proposed. In the process, a Neural Network is used to compute the output state probability distribution for each sample. Using the aforementioned quantities, a graph is designed, whereas state transition constraints are physically imposed into the graph and a set of constraints are used to retrieve a subset of the recording that maximizes the likelihood function, proposed by the user. The developed framework is tested and validated in two applications. In both cases, the system performance is boosted significantly, e.g in heart sound segmentation, sensitivity increases 2.4% when compared to the standard approaches in the literature.