title: Incentive Alignment of Business Processes: a game theoretic approach

abstract: Many definitions of business processes refer to business goals, value creation, or profits/gains of sorts. Nevertheless, the focus of formal methods research on business processes, like the well-known soundness property, lies on correctness with regards to execution semantics of modeling languages. Among others, soundness requires proper completion of process instances. However, the question of whether participants have any interest in working towards completion (or in participating in the process) has not been addressed as of yet. In this work, we investigate whether inter-organizational business processes give participants incentives for achieving the common business goals---in short, whether incentives are aligned with the process. In particular, fair behavior should pay off and efficient completion of tasks should be rewarded. We propose a game-theoretic approach that relies on algorithms for solving stochastic games from the machine learning community. We describe a method for checking incentive alignment of process models with utility annotations for tasks, which can be used for a priori analysis of inter-organizational business processes. Last but not least, we show that the soundness property corresponds to a special case of incentive alignment.


title: The Spectral Norm of Random Lifts of Matrices

abstract: We study the spectral norm of random lifts of matrices. Given an $n\times n$ symmetric matrix $A$, and a centered distribution $\pi$ on $k\times k\ (k\ge 2)$ symmetric matrices with spectral norm at most $1$, let the matrix random lift $A^{(k,\pi)}$ be the random symmetric $kn\times kn$ matrix $(A_{ij}X_{ij})_{1\le i < j \le n}$, where $X_{ij}$ are independent samples from $\pi$. We prove that $$\mathbb{E} \|A^{(k,\pi)}\|\lesssim \max_{i}\sqrt{\sum_j A_{ij}^2}+\max_{ij}|A_{ij}|\sqrt{\log (kn)}.$$ This result can be viewed as an extension of existing spectral bounds on random matrices with independent entries, providing further instances where the multiplicative $\sqrt{\log n}$ factor in the Non-Commutative Khintchine inequality can be removed. As a direct application of our result, we prove an upper bound of $2(1+\epsilon)\sqrt{\Delta}+O(\sqrt{\log(kn)})$ on the new eigenvalues for random $k$-lifts of a fixed $G = (V,E)$ with $|V| = n$ and maximum degree $\Delta$, compared to the previous result of $O(\sqrt{\Delta\log(kn)})$ by Oliveira and the recent breakthrough by Bordenave and Collins which gives $2\sqrt{\Delta-1} + o(1)$ as $k\rightarrow\infty$ for $\Delta$-regular graph $G$.


title: Shape Analysis of HII Regions -- II. Synthetic Observations

abstract: The statistical shape analysis method developed for probing the link between physical parameters and morphologies of Galactic HII regions is applied here to a set of synthetic observations (SOs) of a numerically modelled HII region. The systematic extraction of HII region shape, presented in the first paper of this series, allows for a quantifiable confirmation of the accuracy of the numerical simulation, with respect to the real observational counterparts of the resulting SOs. A further aim of this investigation is to determine whether such SOs can be used for direct interpretation of the observational data, in a future supervised classification scheme based upon HII region shape. The numerical HII region data was the result of photoionisation and radiation pressure feedback of a 34 Msun star, in a 1000 Msun cloud. The SOs analysed herein comprised four evolutionary snapshots (0.1, 0.2, 0.4 and 0.6 Myr), and multiple viewing projection angles. The shape analysis results provided conclusive evidence of the efficacy of the numerical simulations. When comparing the shapes of the synthetic regions to their observational counterparts, the SOs were grouped in amongst the Galactic HII regions by the hierarchical clustering procedure. There was also an association between the evolutionary distribution of regions and the respective groups. This suggested that the shape analysis method could be further developed for morphological classification of HII regions by using a synthetic data training set, with differing initial conditions of well-defined parameters.


title: Embed Me If You Can: A Geometric Perceptron

abstract: Solving geometric tasks involving point clouds by using machine learning is a challenging problem. Standard feed-forward neural networks combine linear or, if the bias parameter is included, affine layers and activation functions. Their geometric modeling is limited, which motivated the prior work introducing the multilayer hypersphere perceptron (MLHP). Its constituent part, i.e., the hypersphere neuron, is obtained by applying a conformal embedding of Euclidean space. By virtue of Clifford algebra, it can be implemented as the Cartesian dot product of inputs and weights. If the embedding is applied in a manner consistent with the dimensionality of the input space geometry, the decision surfaces of the model units become combinations of hyperspheres and make the decision-making process geometrically interpretable for humans. Our extension of the MLHP model, the multilayer geometric perceptron (MLGP), and its respective layer units, i.e., geometric neurons, are consistent with the 3D geometry and provide a geometric handle of the learned coefficients. In particular, the geometric neuron activations are isometric in 3D, which is necessary for rotation and translation equivariance. When classifying the 3D Tetris shapes, we quantitatively show that our model requires no activation function in the hidden layers other than the embedding to outperform the vanilla multilayer perceptron. In the presence of noise in the data, our model is also superior to the MLHP.


title: Deep learning reconstruction of digital breast tomosynthesis images for accurate breast density and patient-specific radiation dose estimation

abstract: The two-dimensional nature of mammography makes estimation of the overall breast density challenging, and estimation of the true patient-specific radiation dose impossible. Digital breast tomosynthesis (DBT), a pseudo-3D technique, is now commonly used in breast cancer screening and diagnostics. Still, the severely limited 3rd dimension information in DBT has not been used, until now, to estimate the true breast density or the patient-specific dose. This study proposes a reconstruction algorithm for DBT based on deep learning specifically optimized for these tasks. The algorithm, which we name DBToR, is based on unrolling a proximal-dual optimization method. The proximal operators are replaced with convolutional neural networks and prior knowledge is included in the model. This extends previous work on a deep learning-based reconstruction model by providing both the primal and the dual blocks with breast thickness information, which is available in DBT. Training and testing of the model were performed using virtual patient phantoms from two different sources. Reconstruction performance, and accuracy in estimation of breast density and radiation dose, were estimated, showing high accuracy (density <+/-3%; dose <+/-20%) without bias, significantly improving on the current state-of-the-art. This work also lays the groundwork for developing a deep learning-based reconstruction algorithm for the task of image interpretation by radiologists.


title: Angular distribution of extensive air showers by TEL array under GELATICA experiment

abstract: Extensive Air Showers (EAS) arrival direction distribution is studied by means of a 4-detector installation in Telavi (TEL array), which is a node of GELATICA net in Georgia. The description of EAS arrival zenith angle distribution within the spheric layer model of the atmosphere and exponential absorption of showers with the air path is used. It is shown that the variation of zenith angles upper cutoff boundary allows a stable estimation of showers absorption path.


title: Understanding the Dynamics of Information Flow During Disaster Response Using Absorbing Markov Chains

abstract: This paper aims to derive a quantitative model to evaluate the impact of information flow on the effectiveness of disaster response. At the core of the model is a specialized absorbing Markov chain that models the process of delivering federal assistance to the community while considering stakeholder interactions and information flow uncertainty. Using the proposed model, the probability of community satisfaction is computed to reflect the effectiveness of disaster response. A hypothetical example is provided to demonstrate the applicability and interpretability of the derived quantitative model. Practically, the research provides governmental stakeholders interpretable insights for evaluating the impact of information flow on their disaster response effectiveness so that critical stakeholders can be targeted proactive actions for enhanced disaster response.


title: Relativistic quantum motion of an electron in spinning cosmic string spacetime in the presence of uniform magnetic field and Aharonov-Bohm potential

abstract: In this manuscript, we study the relativistic quantum mechanics of an electron in external fields in the spinning cosmic string spacetime. We obtain the Dirac equation, write the first and second-order equations from it, and then we solve these equations for bound states. We show that there are bound states solutions for the first-order equation Dirac. For the second-order equation, we show that its wave functions are given in terms of the Kummer functions and we determine the energies of the particle. We examine the behavior of the energies as a function of the physical parameters of the model, such as rotation, curvature, magnetic field, Aharonov-Bohm flux, and quantum numbers. Our study reveals that both curvature and rotation influence more intensely when these parameters have values smaller than 0.3. We also find that, depending on the values of these parameters, there are energy nonpermissible levels.


title: STRATIFY: a comprehensive and versatile MATLAB code for a multilayered sphere

abstract: We present a computer code for calculating near- and far-field electromagnetic properties of multilayered spheres. STRATIFY is one-of-a-kind open-source package which allows for the efficient calculation of electromagnetic near-field, energy density, total electromagnetic energy, radiative and non-radiative decay rates of a dipole emitter located in any (non-absorbing) shell (including a host medium), and fundamental cross-sections of a multilayered sphere, all within a single program. The developed software is typically more than 50 times faster than freely available packages based on boundary-element-method. Because of its speed and broad applicability, our package is a valuable tool for analysis of numerous light scattering problems, including, but not limited to fluorescence enhancement, upconversion, downconversion, second harmonic generation, surface enhanced Raman spectroscopy. The software is available for download from GitLab https://gitlab.com/iliarasskazov/stratify


title: On the Feasibility of Perfect Resilience with Local Fast Failover

abstract: In order to provide a high resilience and to react quickly to link failures, modern computer networks support fully decentralized flow rerouting, also known as local fast failover. In a nutshell, the task of a local fast failover algorithm is to pre-define fast failover rules for each node using locally available information only. These rules determine for each incoming link from which a packet may arrive and the set of local link failures (i.e., the failed links incident to a node), on which outgoing link a packet should be forwarded. Ideally, such a local fast failover algorithm provides a perfect resilience deterministically: a packet emitted from any source can reach any target, as long as the underlying network remains connected. Feigenbaum et al. (ACM PODC 2012) and also Chiesa et al. (IEEE/ACM Trans. Netw. 2017) showed that it is not always possible to provide perfect resilience. Interestingly, not much more is known currently about the feasibility of perfect resilience. This paper revisits perfect resilience with local fast failover, both in a model where the source can and cannot be used for forwarding decisions. We first derive several fairly general impossibility results: By establishing a connection between graph minors and resilience, we prove that it is impossible to achieve perfect resilience on any non-planar graph; furthermore, while planarity is necessary, it is also not sufficient for perfect resilience. On the positive side, we show that graph families closed under link subdivision allow for simple and efficient failover algorithms which simply skip failed links. We demonstrate this technique by deriving perfect resilience for outerplanar graphs and related scenarios, as well as for scenarios where the source and target are topologically close after failures.


title: Limited Lookahead Policies for the Control of Discrete-Event Systems: A Tutorial

abstract: Some problems in discrete-event systems (DES) model large, time-varying state spaces with complex legal languages. To address these problems, Chung et al. introduced limited lookahead policies (LLP) to provide online supervisory control for discrete-event systems. This seminal paper, along with an addendum of technical results, provided the field with a series of very important and powerful results, but in a notationally- and conceptually-dense manner. In this tutorial, we present Chung et al.'s problem formulation for online control and unravel the formal definitions and proofs from their original work with the aim of making the ideas behind limited lookahead accessible to all DES researchers. Finally, we introduce the Air Traffic Control problem as an example of an online control problem and demonstrate the synthesis of LLP supervisors.


title: A note on some reduction formulas for the incomplete beta function and the Lerch transcedent

abstract: We derive new reduction formulas for the incomplete beta function and the Lerch transcendent in terms of elementary functions. As an application, we calculate some new integrals. Also, we use these reduction formulas to test the performance of the algorithms devoted to the numerical evaluation of the incomplete beta function.


title: Between Broadway and the Hudson: A Bijection of Corridor Paths

abstract: We present a substantial generalization of the equinumeracy of grand Dyck paths and Dyck-path prefixes, constrained within a band. The number of constrained paths starting at level $i$ and ending in a window of size $2j+2$ is equal to the number starting at level $j$ and ending in a window of size $2i+2$ centered around the same point. A new encoding of lattice paths provides a bijective proof.


title: Evolution and global charge conservation for polarization singularities emerging from nonhermitian degeneracies

abstract: Core concepts in singular optics, especially the polarization singularity, have rapidly penetrated the surging fields of topological and nonhermitian photonics. For open photonic structures with degeneracies in particular, the polarization singularity would inevitably encounter another sweeping concept of Berry phase. Several investigations have discussed, in an inexplicit way, the connections between both concepts, hinting at that nonzero topological charges for far-field polarizations on a loop is inextricably linked to its nontrivial Berry phase when degeneracies are enclosed. In this work, we reexamine the seminal photonic crystal slab that supports the fundamental two-level nonhermitian degeneracies. Regardless of the invariance of nontrivial Berry phase for different loops enclosing both exceptional points, we demonstrate that the associated polarization fields exhibit topologically inequivalent patterns that are characterized by variant topological charges, including even the trivial scenario of zero charge. It is further revealed that for both bands, the seemingly complex evolutions of polarizations are bounded by the global charge conservation, with extra points of circular polarizations playing indispensable roles. This indicates that tough not directly associated with any local charges, the invariant Berry phase is directly linked to the globally conserved charge, the physical principles underlying which have all been further clarified by a modified Berry-Dennis model. Our work can potentially trigger an avalanche of studies to explore subtle interplays between Berry phase and all sorts of optical singularities, shedding new light on subjects beyond photonics that are related to both Berry phase and singularities.


title: Towards Expedited Impedance Tuning of a Robotic Prosthesis for Personalized Gait Assistance by Reinforcement Learning Control

abstract: Personalizing medical devices such as lower limb wearable robots is challenging. While the initial feasibility of automating the process of knee prosthesis control parameter tuning has been demonstrated in a principled way, the next critical issue is to improve tuning efficiency and speed it up for the human user, in clinic settings, while maintaining human safety. We, therefore, propose a policy iteration with constraint embedded (PICE) method as an innovative solution to the problem under the framework of reinforcement learning. Central to PICE is the use of a projected Bellman equation with a constraint of assuring positive semidefiniteness of performance values during policy evaluation. Additionally, we developed both online and offline PICE implementations that provide additional flexibility for the designer to fully utilize measurement data, either from on-policy or off-policy, to further improve PICE tuning efficiency. Our human subject testing showed that the PICE provided effective policies with significantly reduced tuning time. For the first time, we also experimentally evaluated and demonstrated the robustness of the deployed policies by applying them to different tasks and users. Putting it together, our new way of problem solving has been effective as PICE has demonstrated its potential toward truly automating the process of control parameter tuning for robotic knee prosthesis users.


title: Reserve Price Optimization for First Price Auctions

abstract: The display advertising industry has recently transitioned from second- to first-price auctions as its primary mechanism for ad allocation and pricing. In light of this, publishers need to re-evaluate and optimize their auction parameters, notably reserve prices. In this paper, we propose a gradient-based algorithm to adaptively update and optimize reserve prices based on estimates of bidders' responsiveness to experimental shocks in reserves. Our key innovation is to draw on the inherent structure of the revenue objective in order to reduce the variance of gradient estimates and improve convergence rates in both theory and practice. We show that revenue in a first-price auction can be usefully decomposed into a \emph{demand} component and a \emph{bidding} component, and introduce techniques to reduce the variance of each component. We characterize the bias-variance trade-offs of these techniques and validate the performance of our proposed algorithm through experiments on synthetic data and real display ad auctions data from Google ad exchange.


title: Achieving robustness in classification using optimal transport with hinge regularization

abstract: Adversarial examples have pointed out Deep Neural Networks vulnerability to small local noise. It has been shown that constraining their Lipschitz constant should enhance robustness, but make them harder to learn with classical loss functions. We propose a new framework for binary classification, based on optimal transport, which integrates this Lipschitz constraint as a theoretical requirement. We propose to learn 1-Lipschitz networks using a new loss that is an hinge regularized version of the Kantorovich-Rubinstein dual formulation for the Wasserstein distance estimation. This loss function has a direct interpretation in terms of adversarial robustness together with certifiable robustness bound. We also prove that this hinge regularized version is still the dual formulation of an optimal transportation problem, and has a solution. We also establish several geometrical properties of this optimal solution, and extend the approach to multi-class problems. Experiments show that the proposed approach provides the expected guarantees in terms of robustness without any significant accuracy drop. The adversarial examples, on the proposed models, visibly and meaningfully change the input providing an explanation for the classification.


title: Morita Equivalence of Graph and Ultragraph Leavitt Path Algebras

abstract: The primary purpose of this thesis is to show every ultragraph Leavitt path algebra is Morita equivalent, as a ring, to a graph Leavitt path algebra. Takeshi Katsura, Paul Muhly, Aidan Sims, and Mark Tomforde showed every ultragraph $C^{*}$-algebra is Morita equivalent, in the $C^{*}$-sense, to a graph $C^{*}$-algebra; our result is an algebraic analog of this fact. Further, we will use our result to give an alternate proof for established conditions which guarantee the simplicity of an ultragraph Leavitt path algebra over a field.


title: Squeezing-enhanced communication without a phase reference

abstract: We study the problem of transmitting classical information using quantum Gaussian states on a family of phase-noise channels with a finite decoherence time, such that the phase-reference is lost after $m$ consecutive uses of the transmission line. This problem is relevant for long-distance communication in free space and optical fiber, where phase noise is typically considered as a limiting factor. The Holevo capacity of these channels is always attained with photon-number encodings, challenging with current technology. Hence for coherent-state encodings the optimal rate depends only on the total-energy distribution and we provide upper and lower bounds for all $m$, the latter attainable at low energies with on/off modulation and photodetection. We generalize this lower bound to squeezed-coherent encodings, exhibiting for the first time to our knowledge an unconditional advantage with respect to any coherent encoding for $m=1$ and a considerable advantage with respect to its direct coherent counterpart for $m>1$. This advantage is robust with respect to moderate attenuation, and persists in a regime where Fock encodings with up to two-photon states are also suboptimal. Finally, we show that the use of part of the energy to establish a reference frame is sub-optimal even at large energies. Our results represent a key departure from the case of phase-covariant Gaussian channels and constitute a proof-of-principle of the advantages of using non-classical, squeezed light in a motivated communication setting.


title: Predicting inhibitors for SARS-CoV-2 RNA-dependent RNA polymerase using machine learning and virtual screening

abstract: Global coronavirus disease pandemic (COVID-19) caused by newly identified SARS- CoV-2 coronavirus continues to claim the lives of thousands of people worldwide. The unavailability of specific medications to treat COVID-19 has led to drug repositioning efforts using various approaches, including computational analyses. Such analyses mostly rely on molecular docking and require the 3D structure of the target protein to be available. In this study, we utilized a set of machine learning algorithms and trained them on a dataset of RNA-dependent RNA polymerase (RdRp) inhibitors to run inference analyses on antiviral and anti-inflammatory drugs solely based on the ligand information. We also performed virtual screening analysis of the drug candidates predicted by machine learning models and docked them against the active site of SARS- CoV-2 RdRp, a key component of the virus replication machinery. Based on the ligand information of RdRp inhibitors, the machine learning models were able to identify candidates such as remdesivir and baloxavir marboxil, molecules with documented activity against RdRp of the novel coronavirus. Among the other identified drug candidates were beclabuvir, a non-nucleoside inhibitor of the hepatitis C virus (HCV) RdRp enzyme, and HCV protease inhibitors paritaprevir and faldaprevir. Further analysis of these candidates using molecular docking against the SARS-CoV-2 RdRp revealed low binding energies against the enzyme active site. Our approach also identified anti-inflammatory drugs lupeol, lifitegrast, antrafenine, betulinic acid, and ursolic acid to have potential activity against SARS-CoV-2 RdRp. We propose that the results of this study are considered for further validation as potential therapeutic options against COVID-19.


title: Observation of a higher-order topological bound state in the continuum

abstract: Higher-order topological insulators are a recently discovered class of materials that can possess zero-dimensional localized states regardless of the dimension of the lattice. Here, we experimentally demonstrate that the topological corner-localized modes of higher-order topological insulators can be symmetry protected bound states in the continuum; these states do not hybridize with the surrounding bulk states of the lattice even in the absence of a bulk bandgap. As such, this class of structures has potential applications in confining and controlling light in systems that do not support a complete photonic bandgap.


title: Attentive WaveBlock: Complementarity-enhanced Mutual Networks for Unsupervised Domain Adaptation in Person Re-identification and Beyond

abstract: Unsupervised domain adaptation (UDA) for person re-identification is challenging because of the huge gap between the source and target domain. A typical self-training method is to use pseudo-labels generated by clustering algorithms to iteratively optimize the model on the target domain. However, a drawback to this is that noisy pseudo-labels generally cause trouble in learning. To address this problem, a mutual learning method by dual networks has been developed to produce reliable soft labels. However, as the two neural networks gradually converge, their complementarity is weakened and they likely become biased towards the same kind of noise. This paper proposes a novel light-weight module, the Attentive WaveBlock (AWB), which can be integrated into the dual networks of mutual learning to enhance the complementarity and further depress noise in the pseudo-labels. Specifically, we first introduce a parameter-free module, the WaveBlock, which creates a difference between features learned by two networks by waving blocks of feature maps differently. Then, an attention mechanism is leveraged to enlarge the difference created and discover more complementary features. Furthermore, two kinds of combination strategies, i.e. pre-attention and post-attention, are explored. Experiments demonstrate that the proposed method achieves state-of-the-art performance with significant improvements on multiple UDA person re-identification tasks. We also prove the generality of the proposed method by applying it to vehicle re-identification and image classification tasks. Our codes and models are available at https://github.com/WangWenhao0716/Attentive-WaveBlock.


title: Recurrent Neural Networks for Handover Management in Next-Generation Self-Organized Networks

abstract: In this paper, we discuss a handover management scheme for Next Generation Self-Organized Networks. We propose to extract experience from full protocol stack data, to make smart handover decisions in a multi-cell scenario, where users move and are challenged by deep zones of an outage. Traditional handover schemes have the drawback of taking into account only the signal strength from the serving, and the target cell, before the handover. However, we believe that the expected Quality of Experience (QoE) resulting from the decision of target cell to handover to, should be the driving principle of the handover decision. In particular, we propose two models based on multi-layer many-to-one LSTM architecture, and a multi-layer LSTM AutoEncoder (AE) in conjunction with a MultiLayer Perceptron (MLP) neural network. We show that using experience extracted from data, we can improve the number of users finalizing the download by 18%, and we can reduce the time to download, with respect to a standard event-based handover benchmark scheme. Moreover, for the sake of generalization, we test the LSTM Autoencoder in a different scenario, where it maintains its performance improvements with a slight degradation, compared to the original scenario.


title: MMA Regularization: Decorrelating Weights of Neural Networks by Maximizing the Minimal Angles

abstract: The strong correlation between neurons or filters can significantly weaken the generalization ability of neural networks. Inspired by the well-known Tammes problem, we propose a novel diversity regularization method to address this issue, which makes the normalized weight vectors of neurons or filters distributed on a hypersphere as uniformly as possible, through maximizing the minimal pairwise angles (MMA). This method can easily exert its effect by plugging the MMA regularization term into the loss function with negligible computational overhead. The MMA regularization is simple, efficient, and effective. Therefore, it can be used as a basic regularization method in neural network training. Extensive experiments demonstrate that MMA regularization is able to enhance the generalization ability of various modern models and achieves considerable performance improvements on CIFAR100 and TinyImageNet datasets. In addition, experiments on face verification show that MMA regularization is also effective for feature learning. Code is available at: https://github.com/wznpub/MMA_Regularization.


title: Unresolved Binary Exoplanet Host Stars Fit as Single Stars: Effects on the Stellar Parameters

abstract: In this work we quantify the effect of an unresolved companion star on the derived stellar parameters of the primary star if a blended spectrum is fit assuming the star is single. Fitting tools that determine stellar parameters from spectra typically fit for a single star, but we know that up to half of all exoplanet host stars may have one or more companion stars. We use high-resolution spectra of planet host stars in the Kepler field from the California-Kepler Survey to create simulated binaries; we select 8 stellar pairs and vary the contribution of the secondary star, then determine stellar parameters with SpecMatch-Emp and compare them to the parameters derived for the primary star alone. We find that in most cases the effective temperature, surface gravity, metallicity, and stellar radius derived from the composite spectrum are within 2-3 $\sigma$ of the values determined from the unblended spectrum, but the deviations depend on the properties of the two stars. Relatively bright companion stars that are similar to the primary star have the largest effect on the derived parameters; in these cases the stellar radii can be overestimated by up to 60%. We find that metallicities are generally underestimated, with values up to 8 times smaller than the typical uncertainty in [Fe/H]. Our study shows that follow-up observations are necessary to detect or set limits on stellar companions of planetary host stars so that stellar (and planet) parameters are as accurate as possible.


title: Dipole-dipole interaction driven antiblockade of two Rydberg atoms

abstract: Resonant laser excitation of multiple Rydberg atoms are prohibited, leading to Rydberg blockade, when the long-range van der Waals interactions are stronger than the laser-atom coupling. Rydberg blockade can be violated, i.e. simultaneous excitation of more than one Rydberg atoms, by off-resonant laser excitation, causing an excitation antiblockade. Rydberg antiblockade gives rise to strongly correlated many-body dynamics and spin-orbit coupling, and also finds quantum computation applications. Instead of commonly used van der Waals interactions, we investigate antiblockade dynamics of two Rydberg atoms interacting via dipole-dipole exchange interactions. We study typical situations in current Rydberg atoms experiments, where different types of dipole-dipole interactions can be achieved by varying Rydberg state couplings. Effective Hamiltonian governing underlying antiblockade dynamics is derived. We illustrate that geometric gates can be realized with the Rydberg antiblockade which is robust against decay of Rydberg states. Our study may stimulate new experimental and theoretical exploration of quantum optics and strongly interacting many-body dynamics with Rydberg antiblockade driven by dipole-dipole interactions.


title: Age-structured estimation of COVID-19 ICU demand from low quality data

abstract: We sample aggravated cases following age-structured probabilities from confirmed cases and use ICU occupation data to find a subnotification factor. A logistic fit is then employed to project the progression of the COVID-19 epidemic with plateau scenarios taken from locations that have reached this stage. Finally, the logistic curve found is corrected by the subnotification factor and sampled to project the future demand for ICU beds.


title: A Comparative Study of U-Net Topologies for Background Removal in Histopathology Images

abstract: During the last decade, the digitization of pathology has gained considerable momentum. Digital pathology offers many advantages including more efficient workflows, easier collaboration as well as a powerful venue for telepathology. At the same time, applying Computer-Aided Diagnosis (CAD) on Whole Slide Images (WSIs) has received substantial attention as a direct result of the digitization. The first step in any image analysis is to extract the tissue. Hence, background removal is an essential prerequisite for efficient and accurate results for many algorithms. In spite of the obvious discrimination for human operators, the identification of tissue regions in WSIs could be challenging for computers, mainly due to the existence of color variations and artifacts. Moreover, some cases such as alveolar tissue types, fatty tissues, and tissues with poor staining are difficult to detect. In this paper, we perform experiments on U-Net architecture with different network backbones (different topologies) to remove the background as well as artifacts from WSIs in order to extract the tissue regions. We compare a wide range of backbone networks including MobileNet, VGG16, EfficientNet-B3, ResNet50, ResNext101 and DenseNet121. We trained and evaluated the network on a manually labeled subset of The Cancer Genome Atlas (TCGA) Dataset. EfficientNet-B3 and MobileNet by almost 99% sensitivity and specificity reached the best results.


title: Kotani's Theorem for the Fourier Transform

abstract: In 1991, Shinichi Kotani proved a theorem giving a sufficient condition to conclude that a function $f(x)$ on ${\mathbb Z}^d$ decays like $|x|^{-(d-2)}$ for large $x$, assuming that its Fourier transform $\hat f(k)$ is such that $|k|^{2}\hat f(k)$ is well behaved for $k$ near zero. The proof was not published. We prove an extension of Kotani's Theorem, based on Kotani's unpublished proof.


title: Direct and inverse problems for the matrix Sturm-Liouville operator with the general self-adjoint boundary conditions

abstract: The matrix Sturm-Liouville operator on a finite interval with the boundary conditions in the general self-adjoint form and with the singular potential from the class $W_2^{-1}$ is studied. This operator generalizes Sturm-Liouville operators on geometrical graphs. We investigate structural and asymptotical properties of the spectral data (eigenvalues and weight matrices) of this operator. Furthermore, we prove the uniqueness of recovering the operator from its spectral data, by using the method of spectral mappings.


title: Gravitational baryogenesis of cosmological constant dominated universe

abstract: In this paper, we have studied the gravitational baryogenesis of isotropic and homogeneous universe in the frame-work of general relativity. We investigate an exact and new solution of Einstein's field equations for FRW metric. Our solution represents a transitioning model of the universe which was expanding in decelerated mode and it transit in accelerated mode after dominance of cosmological constant $\Lambda$. We observe that gravitational baryogenesis occurs in the derived universe and derived baryon entropy ratio is in good agreement with its observational value.


title: Privacy Adversarial Network: Representation Learning for Mobile Data Privacy

abstract: The remarkable success of machine learning has fostered a growing number of cloud-based intelligent services for mobile users. Such a service requires a user to send data, e.g. image, voice and video, to the provider, which presents a serious challenge to user privacy. To address this, prior works either obfuscate the data, e.g. add noise and remove identity information, or send representations extracted from the data, e.g. anonymized features. They struggle to balance between the service utility and data privacy because obfuscated data reduces utility and extracted representation may still reveal sensitive information. This work departs from prior works in methodology: we leverage adversarial learning to a better balance between privacy and utility. We design a \textit{representation encoder} that generates the feature representations to optimize against the privacy disclosure risk of sensitive information (a measure of privacy) by the \textit{privacy adversaries}, and concurrently optimize with the task inference accuracy (a measure of utility) by the \textit{utility discriminator}. The result is the privacy adversarial network (\systemname), a novel deep model with the new training algorithm, that can automatically learn representations from the raw data. Intuitively, PAN adversarially forces the extracted representations to only convey the information required by the target task. Surprisingly, this constitutes an implicit regularization that actually improves task accuracy. As a result, PAN achieves better utility and better privacy at the same time! We report extensive experiments on six popular datasets and demonstrate the superiority of \systemname compared with alternative methods reported in prior work.


title: Gapped boundaries and string-like excitations in (3+1)d gauge models of topological phases

abstract: We study lattice Hamiltonian realisations of (3+1)d Dijkgraaf-Witten theory with gapped boundaries. In addition to the bulk loop-like excitations, the Hamiltonian yields bulk dyonic string-like excitations that terminate at gapped boundaries. Using a tube algebra approach, we classify such excitations and derive the corresponding representation theory. Via a dimensional reduction argument, we relate this tube algebra to that describing (2+1)d boundary point-like excitations at interfaces between two gapped boundaries. Such point-like excitations are well known to be encoded into a bicategory of module categories over the input fusion category. Exploiting this correspondence, we define a bicategory that encodes the string-like excitations ending at gapped boundaries, showing that it is a sub-bicategory of the centre of the input bicategory of group-graded 2-vector spaces. In the process, we explain how gapped boundaries in (3+1)d can be labelled by so-called pseudo-algebra objects over this input bicategory.

categories: cond-mat.str-el hep-th math-ph math.MP quant-ph


title: Fast increased fidelity approximate Gibbs samplers for Bayesian Gaussian process regression

abstract: The use of Gaussian processes (GPs) is supported by efficient sampling algorithms, a rich methodological literature, and strong theoretical grounding. However, due to their prohibitive computation and storage demands, the use of exact GPs in Bayesian models is limited to problems containing at most several thousand observations. Sampling requires matrix operations that scale at $\mathcal{O}(n^3),$ where $n$ is the number of unique inputs. Storage of individual matrices scales at $\mathcal{O}(n^2),$ and can quickly overwhelm the resources of most modern computers. To overcome these bottlenecks, we develop a sampling algorithm using $\mathcal{H}$ matrix approximation of the matrices comprising the GP posterior covariance. These matrices can approximate the true conditional covariance matrix within machine precision and allow for sampling algorithms that scale at $\mathcal{O}(n \ \mbox{log}^2 n)$ time and storage demands scaling at $\mathcal{O}(n \ \mbox{log} \ n).$ We also describe how these algorithms can be used as building blocks to model higher dimensional surfaces at $\mathcal{O}(d \ n \ \mbox{log}^2 n)$, where $d$ is the dimension of the surface under consideration, using tensor products of one-dimensional GPs. Though various scalable processes have been proposed for approximating Bayesian GP inference when $n$ is large, to our knowledge, none of these methods show that the approximation's Kullback-Leibler divergence to the true posterior can be made arbitrarily small and may be no worse than the approximation provided by finite computer arithmetic. We describe $\mathcal{H}-$matrices, give an efficient Gibbs sampler using these matrices for one-dimensional GPs, offer a proposed extension to higher dimensional surfaces, and investigate the performance of this fast increased fidelity approximate GP, FIFA-GP, using both simulated and real data sets.


title: Faraday rotation measures of northern-hemisphere pulsars using CHIME/Pulsar

abstract: Using commissioning data from the first year of operation of the Canadian Hydrogen Intensity Mapping Experiment's (CHIME) Pulsar backend system, we conduct a systematic analysis of the Faraday Rotation Measure (RM) of the northern hemisphere pulsars detected by CHIME. We present 55 new RMs as well as obtain improved RM uncertainties for 25 further pulsars. CHIME's low observing frequency and wide bandwidth between 400-800 MHz contribute to the precision of our measurements, whereas the high cadence observation provide extremely high signal-to-noise co-added data. Our results represent a significant increase of the pulsar RM census, particularly regarding the northern hemisphere. These new RMs are for sources that are located in the Galactic plane out to 10 kpc, as well as off the plane to a scale height of ~16 kpc. This improved knowledge of the Faraday sky will contribute to future Galactic large-scale magnetic structure and ionosphere modelling.


title: Quantitative global-local mixing for accessible skew products

abstract: We study global-local mixing for a family of accessible skew products with an exponentially mixing base and non-compact fibers, preserving an infinite measure. For a dense set of almost periodic global observables, we prove rapid mixing; and for a dense set of global observables vanishing at infinity, we prove polynomial mixing. More generally, we relate the speed of mixing to the "low frequency behaviour" of the spectral measure associated to our global observables. Our strategy relies on a careful choice of the spaces of observables and on the study of a family of twisted transfer operators.


title: Impedances and Instabilities

abstract: The concepts of wake fields and impedance are introduced to describe the electromagnetic interaction of a bunch of charged particles with its environment in an particle accelerator. The wake fields can act back on the beam and lead to instabilities, which may limit the achievable current per bunch, the total current, or even both. Some typical examples are used to illustrate the wake function and its basic properties. Then the frequency-domain view of the wake field or impedance is explained, and basic properties of the impedance are derived. The impedance of a cavity mode is illustrated using an equivalent circuit model. The relation of wake field effects to important beam parameters is treated in the rigid beam approximation. Several examples are employed to illustrate the impact of the geometry and the material properties of the vacuum chamber on the impedance. Finally, a basic introduction to beam instabilities based on a head tail model of the beam is given.


title: Intelligent Surfaces for 6G Wireless Networks: A Survey of Optimization and Performance Analysis Techniques

abstract: This paper surveys the optimization frameworks and performance analysis methods for large intelligent surfaces (LIS), which have been emerging as strong candidates to support the sixth-generation wireless physical platforms (6G). Due to their ability to adjust the behavior of interacting electromagnetic (EM) waves through intelligent manipulations of the reflections phase shifts, LIS have shown promising merits at improving the spectral efficiency of wireless networks. In this context, researchers have been recently exploring LIS technology in depth as a means to achieve programmable, virtualized, and distributed wireless network infrastructures. From a system level perspective, LIS have also been proven to be a low-cost, green, sustainable, and energy-efficient solution for 6G systems. This paper provides a unique blend that surveys the principles of operation of LIS, together with their optimization and performance analysis frameworks. The paper first introduces the LIS technology and its physical working principle. Then, it presents various optimization frameworks that aim to optimize specific objectives, namely, maximizing energy efficiency, sum-rate, secrecy-rate, and coverage. The paper afterwards discusses various relevant performance analysis works including capacity analysis, the impact of hardware impairments on capacity, uplink/downlink data rate analysis, and outage probability. The paper further presents the impact of adopting the LIS technology for positioning applications. Finally, we identify numerous exciting open challenges for LIS-aided 6G wireless networks, including resource allocation problems, hybrid radio frequency/visible light communication (RF-VLC) systems, health considerations, and localization.


title: How simplifying and flexible is the simplifying assumption in pair-copula constructions -- analytic answers in dimension three and a glimpse beyond

abstract: Motivated by the increasing popularity and the seemingly broad applicability of pair-copula constructions underlined by numerous publications in the last decade, in this contribution we tackle the unavoidable question on how flexible and simplifying the commonly used `simplifying assumption' is from an analytic perspective and provide answers to two related open questions posed by Nagler and Czado in 2016. Aiming at a simplest possible setup for deriving the main results we first focus on the three-dimensional setting. We prove that the family of simplified copulas is flexible in the sense that it is dense in the set of all three-dimensional co\-pulas with respect to the uniform metric $d_\infty$ - considering stronger notions of convergence like the one induced by the metric $D_1$, by weak conditional convergence, by total variation, or by Kullback-Leibler divergence, however, the family even turn out to be nowhere dense and hence insufficient for any kind of flexible approximation. Furthermore, returning to $d_\infty$ we show that the partial vine copula is never the optimal simplified copula approximation of a given, non-simplified copula $C$, and derive examples illustrating that the corresponding approximation error can be strikingly large and extend to more than 28\% of the diameter of the metric space. Moreover, the mapping $\psi$ assigning each three-dimensional copula its unique partial vine copula turns out to be discontinuous with respect to $d_\infty$ (but continuous with respect to $D_1$ and to weak conditional convergence), implying a surprising sensitivity of partial vine copula approximations. The afore-mentioned main results concerning $d_\infty$ are then extended to the general multivariate setting.


title: Data and Incentives

abstract: "Big data" gives markets access to previously unmeasured characteristics of individual agents. Policymakers must decide whether and how to regulate the use of this data. We study how new data affects incentives for agents to exert effort in settings such as the labor market, where an agent's quality is initially unknown but is forecast from an observable outcome. We show that measurement of a new covariate has a systematic effect on the average effort exerted by agents, with the direction of the effect determined by whether the covariate is informative about long-run quality or about a shock to short-run outcomes. For a class of covariates satisfying a statistical property we call strong homoskedasticity, this effect is uniform across agents. More generally, new measurements can impact agents unequally, and we show that these distributional effects have a first-order impact on social welfare.


title: Parallel-in-Time Simulation of Power Converters Using Multirate PDEs

abstract: This paper presents an efficient numerical algorithm for the simulation of pulse-width modulated power converters via parallelization in time domain. The method applies the multirate partial differential equation approach on the coarse grid of the (two-grid) parallel-in-time algorithm Parareal. Performance of the proposed approach is illustrated via its application to a DC-DC converter.


title: Interpretable, similarity-driven multi-view embeddings from high-dimensional biomedical data

abstract: Similarity-driven multi-view linear reconstruction (SiMLR) is an algorithm that exploits inter-modality relationships to transform large scientific datasets into smaller, more well-powered and interpretable low-dimensional spaces. SiMLR contributes a novel objective function for identifying joint signal, regularization based on sparse matrices representing prior within-modality relationships and an implementation that permits application to joint reduction of large data matrices, each of which may have millions of entries. We demonstrate that SiMLR outperforms closely related methods on supervised learning problems in simulation data, a multi-omics cancer survival prediction dataset and multiple modality neuroimaging datasets. Taken together, this collection of results shows that SiMLR may be applied with default parameters to joint signal estimation from disparate modalities and may yield practically useful results in a variety of application domains.


title: Unique determination of the shape of a scattering screen from a passive measurement

abstract: We consider the problem of fixed frequency acoustic scattering from a sound-soft flat screen. More precisely the obstacle is restricted to a two-dimensional plane and interacting with a arbitrary incident wave, it scatters acoustic waves to three-dimensional space. The model is particularly relevant in the study and design of reflecting sonars and antennas, cases where one cannot assume that the incident wave is a plane wave. Our main result is that given the plane where the screen is located, the far-field pattern produced by any single arbitrary incident wave determines the exact shape of the screen, as long as it is not antisymmetric with respect to the plane. This holds even for screens whose shape is an arbitrary simply connected smooth domain. This is in contrast to earlier work where the incident wave had to be a plane wave, or more recent work where only polygonal scatterers are determined.


title: Avoiding Side Effects in Complex Environments

abstract: Reward function specification can be difficult. Rewarding the agent for making a widget may be easy, but penalizing the multitude of possible negative side effects is hard. In toy environments, Attainable Utility Preservation (AUP) avoided side effects by penalizing shifts in the ability to achieve randomly generated goals. We scale this approach to large, randomly generated environments based on Conway's Game of Life. By preserving optimal value for a single randomly generated reward function, AUP incurs modest overhead while leading the agent to complete the specified task and avoid many side effects. Videos and code are available at https://avoiding-side-effects.github.io/.


title: Efficient democratic decisions via nondeterministic proportional consensus

abstract: Are there voting methods which (i) give everyone, including minorities, an equal share of effective power even if voters act strategically, (ii) promote consensus rather than polarization and inequality, and (iii) do not favour the status quo or rely too much on chance? We show the answer is yes by describing two nondeterministic voting methods, one based on automatic bargaining over lotteries, the other on conditional commitments to approve compromise options. Our theoretical analysis and agent-based simulation experiments suggest that with these, majorities cannot consistently suppress minorities as with deterministic methods, proponents of the status quo cannot block decisions as in consensus-based approaches, the resulting aggregate welfare is comparable to existing methods, and average randomness is lower than for other nondeterministic methods.


title: Imaging Spatially Extended Objects with Interferometers: Mosaicking and the Short Spacing Correction

abstract: Interferometry is a powerful technique for making sensitive, high-fidelity images of the sky, but is limited in its ability to measure extended or diffuse emission. Better images of extended astronomical objects can be obtained by mosaicking together many pointings of the interferometer array. Even better images can be obtained by combining these data with data from a single-dish telescope. This lecture explains commonly practiced techniques for obtaining and analyzing these observations, and the theory behind them.


title: Aplica\c{c}\~oes do Setor John da Gravidade de Horndeski nos Cen\'arios de Brana Negra e Rela\c{c}\~ao de viscosidade/entropia, Mundo Brana e Cosmologia (In Portuguese)

abstract: In this thesis, we approached the applications of the John sector of the Horndeski gravity in the scenarios of black brane, braneworld and cosmology. Specially, we study AdS/CFT correspondence in the context of the Horndeski gravity to evaluate the viscosity/entropy density ratio, as a function of the fixed exponent and parameters of theory. In the braneworld, we address the issue of finding braneworld solutions in a five-dimensional Horndeski gravity and discuss gravity localization. Finally, in the cosmological scenario we finding solutions by using first-order formalism in the Horndeski gravity that constrains the superpotential and implies that no free choice of scalar potential is allowed. The scalar field equation of state reflects well the cosmological scenario at the present time and allows us to conclude that it can simulate the dark energy in the Horndeski gravity.


title: Upper Bounding Rainbow Connection Number by Forest Number

abstract: A path in an edge-colored graph is rainbow if no two edges of it are colored the same, and the graph is rainbow-connected if there is a rainbow path between each pair of its vertices. The minimum number of colors needed to rainbow-connect a graph $G$ is the rainbow connection number of $G$, denoted by $\text{rc}(G)$. A simple way to rainbow-connect a graph $G$ is to color the edges of a spanning tree with distinct colors and then re-use any of these colors to color the remaining edges of $G$. This proves that $\text{rc}(G) \le |V(G)|-1$. We ask whether there is a stronger connection between tree-like structures and rainbow coloring than that is implied by the above trivial argument. For instance, is it possible to find an upper bound of $t(G) -1$ for $\text{rc}(G)$, where $t(G)$ is the number of vertices in the largest induced tree of $G$? The answer turns out to be negative, as there are counter-examples that show that even $c\cdot t(G)$ is not an upper bound for $\text{rc}(G))$ for any given constant $c$. In this work we show that if we consider the forest number $f(G)$, the number of vertices in a maximum induced forest of $G$, instead of $t(G)$, then surprisingly we do get an upper bound. More specifically, we prove that $\text{rc}(G) \leq f(G) + 2$. Our result indicates a stronger connection between rainbow connection and tree-like structures than that was suggested by the simple spanning tree based upper bound.


title: A model for redistributing heat over the surface of irradiated spider companions

abstract: Spider pulsars are binary systems containing an energetic millisecond pulsar that intensely irradiates a closely orbiting low-mass companion. Modelling their companion's optical light curves is essential to the study of the orbital properties of the binary, including the determination of the pulsar mass, characterising the pulsar wind and the star itself. We aim to generalise the traditional direct heating model of irradiation, whereby energy deposited by the pulsar wind into the stellar envelope is locally re-emitted, by introducing heat redistribution via diffusion and convection within the outer stellar envelope. We approximate the irradiated stellar envelope as a two-dimensional shell. This allows us to propose an effective equation of energy conservation that can be solved at a reduced computational cost. We then implement this model in the \texttt{Icarus} software and use evidence sampling to determine the most likely convection and diffusion laws for the light curve of the redback companion of PSR J2215+5135. Redistribution effects concentrate near the terminator line of pulsar irradiation, and can create apparent hot and cold spots. Among the models tested for PSR J2215+5135, we find that all models with heat redistribution are more likely than symmetric direct heating. The best-fitting redistribution model involves diffusion together with a uniformly rotating envelope. However, we caution that all models still present serious systematic effects, and that prior knowledge from pulsar timing, spectroscopy and distance are key to determine with certainty the most accurate redistribution law. We propose an extension of the direct heating framework that allows for exploring a variety of heat redistribution effects. Future work is necessary to determine the relevant laws from first principles and empirically using complementary observations.


title: Stanza: A Nonlinear State Space Model for Probabilistic Inference in Non-Stationary Time Series

abstract: Time series with long-term structure arise in a variety of contexts and capturing this temporal structure is a critical challenge in time series analysis for both inference and forecasting settings. Traditionally, state space models have been successful in providing uncertainty estimates of trajectories in the latent space. More recently, deep learning, attention-based approaches have achieved state of the art performance for sequence modeling, though often require large amounts of data and parameters to do so. We propose Stanza, a nonlinear, non-stationary state space model as an intermediate approach to fill the gap between traditional models and modern deep learning approaches for complex time series. Stanza strikes a balance between competitive forecasting accuracy and probabilistic, interpretable inference for highly structured time series. In particular, Stanza achieves forecasting accuracy competitive with deep LSTMs on real-world datasets, especially for multi-step ahead forecasting.


title: Basis and Dimension of Exponential Vector Space

abstract: Exponential vector space [shortly \emph{evs}] is an algebraic order extension of vector space in the sense that every evs contains a vector space and conversely every vector space can be embedded into such a structure. This evs structure consists of a semigroup structure, a scalar multiplication and a partial order. In this paper we have developed the concepts of basis and dimension of an evs by introducing the ideas of orderly independent set and generating set with the help of partial order and algebraic operations. We have found that like vector space, an evs does not contain basis always. We have established a necessary and sufficient condition for an evs to have a basis. It was shown that equality of dimension is an evs property but the converse is not true. We have studied the dimension of subevs and found that every evs contains a subevs with all possible lower dimensions. Lastly we have computed basis and dimension of some evs which help us to explore the theory of basis by creating counter examples in different aspects.


title: Multi-Agent Reinforcement Learning in Stochastic Networked Systems

abstract: We study multi-agent reinforcement learning (MARL) in a stochastic network of agents. The objective is to find localized policies that maximize the (discounted) global reward. In general, scalability is a challenge in this setting because the size of the global state/action space can be exponential in the number of agents. Scalable algorithms are only known in cases where dependencies are static, fixed and local, e.g., between neighbors in a fixed, time-invariant underlying graph. In this work, we propose a Scalable Actor Critic framework that applies in settings where the dependencies can be non-local and stochastic, and provide a finite-time error bound that shows how the convergence rate depends on the speed of information spread in the network. Additionally, as a byproduct of our analysis, we obtain novel finite-time convergence results for a general stochastic approximation scheme and for temporal difference learning with state aggregation, which apply beyond the setting of MARL in networked systems.


title: Exact solution of the Floquet-PXP cellular automaton

abstract: We study the dynamics of a bulk deterministic Floquet model, the Rule 201 synchronous one-dimensional reversible cellular automaton (RCA201). The system corresponds to a deterministic, reversible, and discrete version of the PXP model, whereby a site flips only if both its nearest neighbours are unexcited. We show that the RCA201/Floquet-PXP model exhibits ballistic propagation of interacting quasiparticles - or solitons - corresponding to the domain walls between non-trivial three-fold vacuum states. Starting from the quasiparticle picture, we find the exact matrix product state form of the non-equilibrium stationary state for a range of boundary conditions, including both periodic and stochastic. We discuss further implications of the integrability of the model.


title: Robust Identification of Differential Equations by Numerical Techniques from a Single Set of Noisy Observation

abstract: We propose robust methods to identify underlying Partial Differential Equation (PDE) from a given set of noisy time dependent data. We assume that the governing equation is a linear combination of a few linear and nonlinear differential terms in a prescribed dictionary. Noisy data make such identification particularly challenging. Our objective is to develop methods which are robust against a high level of noise, and to approximate the underlying noise-free dynamics well. We first introduce a Successively Denoised Differentiation (SDD) scheme to stabilize the amplified noise in numerical differentiation. SDD effectively denoises the given data and the corresponding derivatives. Secondly, we present two algorithms for PDE identification: Subspace pursuit Time evolution error (ST) and Subspace pursuit Cross-validation (SC). Our general strategy is to first find a candidate set using the Subspace Pursuit (SP) greedy algorithm, then choose the best one via time evolution or cross validation. ST uses multi-shooting numerical time evolution and selects the PDE which yields the least evolution error. SC evaluates the cross-validation error in the least squares fitting and picks the PDE that gives the smallest validation error. We present a unified notion of PDE identification error to compare the objectives of related approaches. We present various numerical experiments to validate our methods. Both methods are efficient and robust to noise.


title: GR decompositions and their relations to Cholesky-like factorizations

abstract: For a given matrix, we are interested in computing GR decompositions $A=GR$, where $G$ is an isometry with respect to given scalar products. The orthogonal QR decomposition is the representative for the Euclidian scalar product. For a signature matrix, a respective factorization is given as the hyperbolic QR decomposition. Considering a skew-symmetric matrix leads to the symplectic QR decomposition. The standard approach for computing GR decompositions is based on the successive elimination of subdiagonal matrix entries. For the hyperbolic and symplectic case, this approach does in general not lead to a satisfying numerical accuracy. An alternative approach computes the QR decomposition via a Cholesky factorization, but also has bad stability. It is improved by repeating the procedure a second time. In the same way, the hyperbolic and the symplectic QR decomposition are related to the $LDL^T$ and a skew-symmetric Cholesky-like factorization. We show that methods exploiting this connection can provide better numerical stability than elimination-based approaches.


title: Self-consistent tomography and measurement-device independent cryptography

abstract: A recurring problem in quantum mechanics is to estimate either the state of a quantum system or the measurement operator applied to it. If we wish to estimate both, then the difficulty is that the state and the measurement always appear together: to estimate the state, we must use a measurement; to estimate the measurement operator, we must use a state. The data of such quantum estimation experiments come in the form of measurement frequencies. Ideally, the measured average frequencies can be attributed to an average state and an average measurement operator. If this is not the case, we have correlated state-preparation-and-measurement (SPAM) errors. We extend some tests developed to detect such correlated errors to apply to a cryptographic scenario in which two parties trust their individual states but not the measurement performed on the joint state.


title: Generalization error in high-dimensional perceptrons: Approaching Bayes error with convex optimization

abstract: We consider a commonly studied supervised classification of a synthetic dataset whose labels are generated by feeding a one-layer neural network with random iid inputs. We study the generalization performances of standard classifiers in the high-dimensional regime where $\alpha=n/d$ is kept finite in the limit of a high dimension $d$ and number of samples $n$. Our contribution is three-fold: First, we prove a formula for the generalization error achieved by $\ell_2$ regularized classifiers that minimize a convex loss. This formula was first obtained by the heuristic replica method of statistical physics. Secondly, focussing on commonly used loss functions and optimizing the $\ell_2$ regularization strength, we observe that while ridge regression performance is poor, logistic and hinge regression are surprisingly able to approach the Bayes-optimal generalization error extremely closely. As $\alpha \to \infty$ they lead to Bayes-optimal rates, a fact that does not follow from predictions of margin-based generalization error bounds. Third, we design an optimal loss and regularizer that provably leads to Bayes-optimal generalization error.

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


title: ScoreGAN: A Fraud Review Detector based on Multi Task Learning of Regulated GAN with Data Augmentation

abstract: The promising performance of Deep Neural Networks (DNNs) in text classification, has attracted researchers to use them for fraud review detection. However, the lack of trusted labeled data has limited the performance of the current solutions in detecting fraud reviews. The Generative Adversarial Network (GAN) as a semi-supervised method has demonstrated to be effective for data augmentation purposes. The state-of-the-art solutions utilize GANs to overcome the data scarcity problem. However, they fail to incorporate the behavioral clues in fraud generation. Additionally, state-of-the-art approaches overlook the possible bot-generated reviews in the dataset. Finally, they also suffer from a common limitation in scalability and stability of the GAN, slowing down the training procedure. In this work, we propose ScoreGAN for fraud review detection that makes use of both review text and review rating scores in the generation and detection process. Scores are incorporated through Information Gain Maximization (IGM) into the loss function for three reasons. One is to generate score-correlated reviews based on the scores given to the generator. Second, the generated reviews are employed to train the discriminator, so the discriminator can correctly label the possible bot-generated reviews through joint representations learned from the concatenation of GLobal Vector for Word representation (GLoVe) extracted from the text and the score. Finally, it can be used to improve the stability and scalability of the GAN. Results show that the proposed framework outperformed the existing state-of-the-art framework, namely FakeGAN, in terms of AP by 7\%, and 5\% on the Yelp and TripAdvisor datasets, respectively.


title: Superconducting radio-frequency cavity fault classification using machine learning at Jefferson Laboratory

abstract: We report on the development of machine learning models for classifying C100 superconducting radio-frequency (SRF) cavity faults in the Continuous Electron Beam Accelerator Facility (CEBAF) at Jefferson Lab. CEBAF is a continuous-wave recirculating linac utilizing 418 SRF cavities to accelerate electrons up to 12 GeV through 5-passes. Of these, 96 cavities (12 cryomodules) are designed with a digital low-level RF system configured such that a cavity fault triggers waveform recordings of 17 RF signals for each of the 8 cavities in the cryomodule. Subject matter experts (SME) are able to analyze the collected time-series data and identify which of the eight cavities faulted first and classify the type of fault. This information is used to find trends and strategically deploy mitigations to problematic cryomodules. However manually labeling the data is laborious and time-consuming. By leveraging machine learning, near real-time (rather than post-mortem) identification of the offending cavity and classification of the fault type has been implemented. We discuss performance of the ML models during a recent physics run. Results show the cavity identification and fault classification models have accuracies of 84.9% and 78.2%, respectively.


title: A Primer on Large Intelligent Surface (LIS) for Wireless Sensing in an Industrial Setting

abstract: One of the beyond-5G developments that is often highlighted is the integration of wireless communication and radio sensing. This paper addresses the potential of communication-sensing integration of Large Intelligent Surfaces (LIS) in an exemplary Industry 4.0 scenario. Besides the potential for high throughput and efficient multiplexing of wireless links, an LIS can offer a high-resolution rendering of the propagation environment. This is because, in an indoor setting, it can be placed in proximity to the sensed phenomena, while the high resolution is offered by densely spaced tiny antennas deployed over a large area. By treating an LIS as a radio image of the environment, we develop sensing techniques that leverage the usage of computer vision combined with machine learning. We test these methods for a scenario where we need to detect whether an industrial robot deviates from a predefined route. The results show that the LIS-based sensing offers high precision and has a high application potential in indoor industrial environments.


title: Convergence of an iterative scheme for the Monge-Amp\`ere eigenvalue problem with general initial data

abstract: In this note, we revisit an iterative scheme, due to Abedin and Kitagawa (Inverse Iteration for the Monge-Amp\`ere Eigenvalue Problem, Proc. Amer. Math. Soc. 148 (2020), no. 11, 4875--4886), to solve the Monge-Amp\`ere eigenvalue problem on a general bounded convex domain. Using a nonlinear integration by parts, we show that the scheme converges for all convex initial data having finite and nonzero Rayleigh quotient to a nonzero Monge-Amp\`ere eigenfunction.


title: Tidal disruption events in the first billion years of a galaxy

abstract: Accretion of stars on massive black holes (MBHs) can feed MBHs and generate tidal disruption events (TDEs). We introduce a new physically motivated model to self-consistently treat TDEs in cosmological simulations, and apply it to the assembly of a galaxy with final mass $3\times 10^{10}\mathrm{\, M}_\odot$ at $z=6$. This galaxy exhibits a TDE rate of $\sim 10^{-5}\mathrm{\, yr}^{-1}$, consistent with local observations but already in place when the Universe was one billion year old. A fraction of the disrupted stars participate in the growth of MBHs, dominating it until the MBH reaches mass $\sim 5 \times 10^5 \mathrm{\, M}_\odot$, but their contribution then becomes negligible compared to gas. TDEs could be a viable mechanism to grow light MBH seeds, but fewer TDEs are expected when the MBH becomes sufficiently massive to reach the luminosity of, and be detected as, an active galactic nucleus. Galaxy mergers bring multiple MBHs in the galaxy, resulting in an enhancement of the global TDE rate in the galaxy by $\sim 1$ order of magnitude during $100\mathrm{\, Myr}$ around mergers. This enhancement is not on the central MBH, but caused by the presence of MBHs in the infalling galaxies. This is the first self consistent study of TDEs in a cosmological environment and highlights that accretion of stars and TDEs are a natural process occurring in a Milky~Way-mass galaxy at early cosmic times.


title: Optimally Deceiving a Learning Leader in Stackelberg Games

abstract: Recent results in the ML community have revealed that learning algorithms used to compute the optimal strategy for the leader to commit to in a Stackelberg game, are susceptible to manipulation by the follower. Such a learning algorithm operates by querying the best responses or the payoffs of the follower, who consequently can deceive the algorithm by responding as if his payoffs were much different than what they actually are. For this strategic behavior to be successful, the main challenge faced by the follower is to pinpoint the payoffs that would make the learning algorithm compute a commitment so that best responding to it maximizes the follower's utility, according to his true payoffs. While this problem has been considered before, the related literature only focused on the simplified scenario in which the payoff space is finite, thus leaving the general version of the problem unanswered. In this paper, we fill in this gap, by showing that it is always possible for the follower to compute (near-)optimal payoffs for various scenarios about the learning interaction between leader and follower.


title: Exploring Category-Agnostic Clusters for Open-Set Domain Adaptation

abstract: Unsupervised domain adaptation has received significant attention in recent years. Most of existing works tackle the closed-set scenario, assuming that the source and target domains share the exactly same categories. In practice, nevertheless, a target domain often contains samples of classes unseen in source domain (i.e., unknown class). The extension of domain adaptation from closed-set to such open-set situation is not trivial since the target samples in unknown class are not expected to align with the source. In this paper, we address this problem by augmenting the state-of-the-art domain adaptation technique, Self-Ensembling, with category-agnostic clusters in target domain. Specifically, we present Self-Ensembling with Category-agnostic Clusters (SE-CC) -- a novel architecture that steers domain adaptation with the additional guidance of category-agnostic clusters that are specific to target domain. These clustering information provides domain-specific visual cues, facilitating the generalization of Self-Ensembling for both closed-set and open-set scenarios. Technically, clustering is firstly performed over all the unlabeled target samples to obtain the category-agnostic clusters, which reveal the underlying data space structure peculiar to target domain. A clustering branch is capitalized on to ensure that the learnt representation preserves such underlying structure by matching the estimated assignment distribution over clusters to the inherent cluster distribution for each target sample. Furthermore, SE-CC enhances the learnt representation with mutual information maximization. Extensive experiments are conducted on Office and VisDA datasets for both open-set and closed-set domain adaptation, and superior results are reported when comparing to the state-of-the-art approaches.


title: Learning a Unified Sample Weighting Network for Object Detection

abstract: Region sampling or weighting is significantly important to the success of modern region-based object detectors. Unlike some previous works, which only focus on "hard" samples when optimizing the objective function, we argue that sample weighting should be data-dependent and task-dependent. The importance of a sample for the objective function optimization is determined by its uncertainties to both object classification and bounding box regression tasks. To this end, we devise a general loss function to cover most region-based object detectors with various sampling strategies, and then based on it we propose a unified sample weighting network to predict a sample's task weights. Our framework is simple yet effective. It leverages the samples' uncertainty distributions on classification loss, regression loss, IoU, and probability score, to predict sample weights. Our approach has several advantages: (i). It jointly learns sample weights for both classification and regression tasks, which differentiates it from most previous work. (ii). It is a data-driven process, so it avoids some manual parameter tuning. (iii). It can be effortlessly plugged into most object detectors and achieves noticeable performance improvements without affecting their inference time. Our approach has been thoroughly evaluated with recent object detection frameworks and it can consistently boost the detection accuracy. Code has been made available at \url{https://github.com/caiqi/sample-weighting-network}.


title: Cosmological analogies in the search for new physics in high-energy collisions

abstract: In this paper, analogies between multiparticle production in high-energy collisions and the time evolution of the early universe are discussed. A common explanation is put forward under the assumption of an unconventional early state: a rapidly expanding universe before recombination (last scattering surface), followed by the CMB, later evolving up to present days, versus the formation of hidden/dark states in hadronic collisions followed by a conventional QCD parton shower yielding final-state particles. In particular, long-range angular correlations are considered pointing out deep connections between the two physical cases potentially useful for the discovery of new physics.


title: Transferring and Regularizing Prediction for Semantic Segmentation

abstract: Semantic segmentation often requires a large set of images with pixel-level annotations. In the view of extremely expensive expert labeling, recent research has shown that the models trained on photo-realistic synthetic data (e.g., computer games) with computer-generated annotations can be adapted to real images. Despite this progress, without constraining the prediction on real images, the models will easily overfit on synthetic data due to severe domain mismatch. In this paper, we novelly exploit the intrinsic properties of semantic segmentation to alleviate such problem for model transfer. Specifically, we present a Regularizer of Prediction Transfer (RPT) that imposes the intrinsic properties as constraints to regularize model transfer in an unsupervised fashion. These constraints include patch-level, cluster-level and context-level semantic prediction consistencies at different levels of image formation. As the transfer is label-free and data-driven, the robustness of prediction is addressed by selectively involving a subset of image regions for model regularization. Extensive experiments are conducted to verify the proposal of RPT on the transfer of models trained on GTA5 and SYNTHIA (synthetic data) to Cityscapes dataset (urban street scenes). RPT shows consistent improvements when injecting the constraints on several neural networks for semantic segmentation. More remarkably, when integrating RPT into the adversarial-based segmentation framework, we report to-date the best results: mIoU of 53.2%/51.7% when transferring from GTA5/SYNTHIA to Cityscapes, respectively.


title: Quasicrystalline electronic states in twisted bilayers and the effects of interlayer and sublattice symmetries

abstract: We study the electronic structure of quasicrystals composed of incommensurate stacks of atomic layers. We consider two systems: a pair of square lattices with a relative twist angle of $\theta=45^\circ$ and a pair of hexagonal lattices with a relative twist angle of $\theta=30^\circ$, with various interlayer interaction strengths. This constitutes every two-dimensional bilayer quasicrystal system. We investigate the resonant coupling governing the quasicrystalline order in each quasicrystal symmetry, and calculate the quasi-band dispersion. The resonant interaction emerges in bilayer quasicrystals if all the dominant interlayer interactions occur between the atomic orbitals that have the same magnetic quantum number. Thus, not only the quasicrystal composed of the widely studied graphene, but also those composed of transition metal dichalcogenides will exhibit the quasicrystalline states. We find that some quasicrystalline states, which are usually obscured by decoupled monolayer states, are more prominent, i.e., "exposed", in the systems with strong interlayer interaction. We also show that we can switch the states between quasicrystalline configuration and its layer components, by turning on and off the interlayer symmetry.


title: Primes in arithmetic progressions to large moduli I: Fixed residue classes

abstract: We prove new mean value theorems for primes in arithmetic progressions to moduli larger than $x^{1/2}$. Our main result shows that the primes are equidistributed for a fixed residue class over all moduli of size $x^{1/2+\delta}$ with a 'convenient sized' factor. As a consequence, the expected asymptotic holds for all but $O(\delta Q)$ moduli $q\sim Q=x^{1/2+\delta}$ and we get results for moduli as large as $x^{11/21}$. Our proof extends previous techniques of Bombieri, Fouvry, Friedlander and Iwaniec by incorporating new ideas inspired by amplification methods. We combine these with techniques of Zhang and Polymath tailored to our application. In particular, we ultimately rely on exponential sum bounds coming from the spectral theory of automorphic forms (the Kuznetsov trace formula) or from algebraic geometry (Weil and Deligne style estimates).


title: Spectral Image Segmentation with Global Appearance Modeling

abstract: We introduce a new spectral method for image segmentation that incorporates long range relationships for global appearance modeling. The approach combines two different graphs, one is a sparse graph that captures spatial relationships between nearby pixels and another is a dense graph that captures pairwise similarity between all pairs of pixels. We extend the spectral method for Normalized Cuts to this setting by combining the transition matrices of Markov chains associated with each graph. We also derive an efficient method for sparsifying the dense graph of appearance relationships. This leads to a practical algorithm for segmenting high-resolution images. The resulting method can segment challenging images without any filtering or pre-processing.


title: Dynamically Stable Infinite-Width Limits of Neural Classifiers

abstract: Recent research has been focused on two different approaches to studying neural networks training in the limit of infinite width (1) a mean-field (MF) and (2) a constant neural tangent kernel (NTK) approximations. These two approaches have different scaling of hyperparameters with the width of a network layer and as a result, different infinite-width limit models. We propose a general framework to study how the limit behavior of neural models depends on the scaling of hyperparameters with network width. Our framework allows us to derive scaling for existing MF and NTK limits, as well as an uncountable number of other scalings that lead to a dynamically stable limit behavior of corresponding models. However, only a finite number of distinct limit models are induced by these scalings. Each distinct limit model corresponds to a unique combination of such properties as boundedness of logits and tangent kernels at initialization or stationarity of tangent kernels. Existing MF and NTK limit models, as well as one novel limit model, satisfy most of the properties demonstrated by finite-width models. We also propose a novel initialization-corrected mean-field limit that satisfies all properties noted above, and its corresponding model is a simple modification for a finite-width model.


title: Normal Reflection Subgroups

abstract: We study normal reflection subgroups of complex reflection groups. Our point of view leads to a refinement of a theorem of Orlik and Solomon to the effect that the generating function for fixed-space dimension over a reflection group is a product of linear factors involving generalized exponents. Our refinement gives a uniform proof and generalization of a recent theorem of the second author.


title: Flexsealing BGP Against Route Leaks: Peerlock Active Measurement and Analysis

abstract: BGP route leaks frequently precipitate serious disruptions to interdomain routing. These incidents have plagued the Internet for decades while deployment and usability issues cripple efforts to mitigate the problem. Peerlock, introduced in 2016, addresses route leaks with a new approach. Peerlock enables filtering agreements between transit providers to protect their own networks without the need for broad cooperation or a trust infrastructure. We outline the Peerlock system and one variant, Peerlock-lite, and conduct live Internet experiments to measure their deployment on the control plane. Our measurements find evidence for significant Peerlock protection between Tier 1 networks in the peering clique, where 48% of potential Peerlock filters are deployed, and reveal that many other networks also deploy filters against Tier 1 leaks. To guide further deployment, we also quantify Peerlock's impact on route leaks both at currently observed levels and under hypothetical future deployment scenarios via BGP simulation. These experiments reveal present Peerlock deployment restricts Tier 1 leak export to 10% or fewer networks for 40% of simulated leaks. Strategic additional Peerlock-lite deployment at all large ISPs (fewer than 1% of all networks), in tandem with Peerlock within the peering clique as deployed, completely mitigates 80% of simulated Tier 1 route leaks.


title: Helium abundances and its radial gradient from the spectra of HII regions and ring nebulae of the Milky Way

abstract: We determine the radial abundance gradient of helium in the disc of the Galaxy from published spectra of 19 $\text{H}\thinspace \text{II}$ regions and ring nebulae surrounding massive O stars. We revise the Galactocentric distances of the objects considering {\it Gaia} DR2 parallaxes and determine the physical conditions and the ionic abundance of He$^{+}$ in a homogeneous way, using between 3 and 10 $\text{He}\thinspace \text{I}$ recombination lines in each object. We estimate the total He abundance of the nebulae and its radial abundance gradient using four different ICF(He) schemes. The slope of the gradient is always negative and weakly dependent on the ICF(He) scheme, especially when only the objects with log($\eta$) $<$ 0.9 are considered. The slope values go from $-$0.0078 to $-$0.0044 dex kpc$^{-1}$, consistent with the predictions of chemical evolution models of the Milky Way and chemodynamical simulations of disc galaxies. Finally, we estimate the abundance deviations of He, O and N in a sample of ring nebulae around Galactic WR stars, finding a quite similar He overabundance of about +0.24 $\pm$ 0.11 dex in three stellar ejecta ring nebulae.


title: Topologically twisted index of $T[SU(N)]$ at large $N$

abstract: We compute, in the large $N$ limit, the topologically twisted index of the 3d $T[SU(N)]$ theory, namely the partition function on $\Sigma_{\mathfrak{g}} \times S^1$, with a topological twist on the Riemann surface $\Sigma_{\mathfrak{g}}$. To provide an expression for this quantity, we take advantage of some recent results obtained for five dimensional quiver gauge theories. In case of a universal twist, we correctly reproduce the entropy of the universal black hole that can be embedded in the holographically dual solution.


title: Quantization of the charge in Coulomb plus harmonic potential

abstract: We consider two models where the wave equation can be reduced to the effective Schr\"odinger equation whose potential contains both harmonic and the Coulomb terms, $\omega^{2}r^{2}-a/r$. The equation reduces to the biconfluent Heun's equation, and we find that the charge as well as the energy must be quantized and state dependent. We also find that two quantum numbers are necessary to count radial degrees of freedom and suggest that this is a general feature of differential equation with higher singularity like the Heun's equation.


title: Online Learning in Iterated Prisoner's Dilemma to Mimic Human Behavior

abstract: As an important psychological and social experiment, the Iterated Prisoner's Dilemma (IPD) treats the choice to cooperate or defect as an atomic action. We propose to study the behaviors of online learning algorithms in the Iterated Prisoner's Dilemma (IPD) game, where we investigate the full spectrum of reinforcement learning agents: multi-armed bandits, contextual bandits and reinforcement learning. We evaluate them based on a tournament of iterated prisoner's dilemma where multiple agents can compete in a sequential fashion. This allows us to analyze the dynamics of policies learned by multiple self-interested independent reward-driven agents, and also allows us study the capacity of these algorithms to fit the human behaviors. Results suggest that considering the current situation to make decision is the worst in this kind of social dilemma game. Multiples discoveries on online learning behaviors and clinical validations are stated, as an effort to connect artificial intelligence algorithms with human behaviors and their abnormal states in neuropsychiatric conditions.


title: Asymptotic Errors for Teacher-Student Convex Generalized Linear Models (or : How to Prove Kabashima's Replica Formula)

abstract: There has been a recent surge of interest in the study of asymptotic reconstruction performance in various cases of generalized linear estimation problems in the teacher-student setting, especially for the case of i.i.d standard normal matrices. Here, we go beyond these matrices, and prove an analytical formula for the reconstruction performance of convex generalized linear models with rotationally-invariant data matrices with arbitrary bounded spectrum, rigorously confirming, under suitable assumptions, a conjecture originally derived using the replica method from statistical physics. The proof is achieved by leveraging on message passing algorithms and the statistical properties of their iterates, allowing to characterize the asymptotic empirical distribution of the estimator. For sufficiently strongly convex problems, we show that the two-layer vector approximate message passing algorithm (2-MLVAMP) converges, where the convergence analysis is done by checking the stability of an equivalent dynamical system, which gives the result for such problems. We then show that, under a concentration assumption, an analytical continuation may be carried out to extend the result to convex (non-strongly) problems. We illustrate our claim with numerical examples on mainstream learning methods such as sparse logistic regression and linear support vector classifiers, showing excellent agreement between moderate size simulation and the asymptotic prediction.

categories: stat.ML cond-mat.dis-nn cs.IT cs.LG math.IT math.PR


title: Finding Wombling Boundaries in LHC Data with Voronoi and Delaunay Tessellations

abstract: We address the problem of finding a wombling boundary in point data generated by a general Poisson point process, a specific example of which is an LHC event sample distributed in the phase space of a final state signature, with the wombling boundary created by some new physics. We discuss the use of Voronoi and Delaunay tessellations of the point data for estimating the local gradients and investigate methods for sharpening the boundaries by reducing the statistical noise. The outcome from traditional wombling algorithms is a set of boundary cell candidates with relatively large gradients, whose spatial properties must then be scrutinized in order to construct the boundary and evaluate its significance. Here we propose an alternative approach where we simultaneously form and evaluate the significance of all possible boundaries in terms of the total gradient flux. We illustrate our method with several toy examples of both straight and curved boundaries with varying amounts of signal present in the data.

categories: hep-ph hep-ex physics.comp-ph physics.data-an


title: Gauge Principle and Gauge Invariance in Quantum Two-Level Systems

abstract: The quantum Rabi model is a widespread description for the coupling between a two-level system and a quantized single mode of an electromagnetic resonator. Issues about this model's gauge invariance have been raised. These issues become evident when the light-matter interaction reaches the so-called ultrastrong coupling regime. Recently, a modified quantum Rabi model able to provide gauge-invariant physical results in any interaction regime was introduced [Nature Physics 15, 803 (2019)]. Here we provide an alternative derivation of this result, based on the implementation in two-state systems of the gauge principle, which is the principle from which all the fundamental interactions in quantum field theory are derived. The adopted procedure can be regarded as the two-site version of the general method used to implement the gauge principle in lattice gauge theories. Applying this method, we also obtain the gauge-invariant quantum Rabi model for asymmetric two-state systems, and the multi-mode gauge-invariant quantum Rabi model beyond the dipole approximation.


title: The agglomeration and dispersion dichotomy of human settlements on Earth

abstract: Human settlements on Earth are scattered in a multitude of shapes, sizes and spatial arrangements. These patterns are often not random but a result of complex geographical, cultural, economic and historical processes that have profound human and ecological impacts. However, little is known about the global distribution of these patterns and the spatial forces that creates them. This study analyses human settlements from high-resolution satellite imagery and provides a global classification of spatial patterns. We find two emerging classes, namely agglomeration and dispersion. In the former, settlements are fewer than expected based on the predictions of scaling theory, while an unexpectedly high number of settlements characterizes the latter. Our global classification of spatial patterns correlates with some urban outcomes, such as the amount of CO2 emitted for transportation, providing insights into the relationship between land use patterns and socio-economic and environmental indicators. To explain the observed spatial patterns, we also propose a model that combines two agglomeration forces and simulates human settlements' historical growth. Our results show that our model accurately matches the observed global classification (F1: 0.73), helps to understand and estimate the growth of human settlements and, in turn, the distribution and physical dynamics of all human settlements on Earth, from small villages to cities.


title: The Limits of Resolution and Dose for Aberration-Corrected Electron Tomography

abstract: Aberration-corrected electron microscopy can resolve the smallest atomic bond-lengths in nature. However, the high-convergence angles that enable spectacular resolution in 2D have unknown 3D resolution limits for all but the smallest objects ($< \sim$8nm). We show aberration-corrected electron tomography offers new limits for 3D imaging by measuring several focal planes at each specimen tilt. We present a theoretical foundation for aberration-corrected electron tomography by establishing analytic descriptions for resolution, sampling, object size, and dose---with direct analogy to the Crowther-Klug criterion. Remarkably, aberration-corrected scanning transmission electron tomography can measure complete 3D specimen structure of unbounded object sizes up to a specified cutoff resolution. This breaks the established Crowther limit when tilt increments are twice the convergence angle or smaller. Unprecedented 3D resolution is achievable across large objects. Atomic 3D imaging (1$\unicode{xC5}$) is allowed across extended objects larger than depth-of-focus (e.g. $>$ 20nm) using available microscopes and modest specimen tilting ($<$ 3$^\circ$). Furthermore, aberration-corrected tomography follows the rule of dose-fractionation where a specified total dose can be divided among tilts and defoci.