title: The backward Euler-Maruyama method for invariant measures of stochastic differential equations with super-linear coefficients

abstract: The backward Euler-Maruyama (BEM) method is employed to approximate the invariant measure of stochastic differential equations, where both the drift and the diffusion coefficient are allowed to grow super-linearly. The existence and uniqueness of the invariant measure of the numerical solution generated by the BEM method are proved and the convergence of the numerical invariant measure to the underlying one is shown. Simulations are provided to illustrate the theoretical results and demonstrate the application of our results in the area of system control.


title: Gauged $\mathcal{N}=3,D=4$ Supergravity: a new web of marginally connected vacua

abstract: We analyze the vacuum structure of $\mathcal{N}=3,D=4$ supergravity coupled to 9 vector multiplets with gauge group ${\rm SO}(3)\times {\rm SU}(3)$. Aside from the central $\mathcal{N}=3$ AdS$_4$ vacuum at the origin, on which the supermultiplet structure reproduces the massless sector of M-theory compactified on $\mathrm{N^{0,1,0}}$, we find a rich structure of AdS$_4$ vacua preserving $\mathcal{N}=0,1,2,3$ supersymmetry. These new vacua are arranged in a manifold spanned by scalar fields corresponding to exactly marginal deformations of the dual CFT. This manifold has the form $T^3/K$, where $K$ is a discrete subgroup of the gauge group: $\mathcal{N}=3,2$ and $1$ vacua correspond, respectively, to a point, a line and a surface in the three-dimensional vacuum manifold. We study RG flows from the central $\mathcal{N}=3$ vacuum and elaborate on the possible higher dimensional origin of the new vacua. For the reader's convenience we also provide a review of the embedding tensor formulation of $D=4$, $\mathcal{N}=3$ gauged supergravities. In particular we provide formulas involving the fermion shift tensors and mass matrices in $\mathcal{N}=3$ theories, which can be applied to a generic gauging.


title: Polarized x-rays constrain the disk-jet geometry in the black hole x-ray binary Cygnus X-1

abstract: A black hole x-ray binary (XRB) system forms when gas is stripped from a normal star and accretes onto a black hole, which heats the gas sufficiently to emit x-rays. We report a polarimetric observation of the XRB Cygnus X-1 using the Imaging x-ray Polarimetry Explorer. The electric field position angle aligns with the outflowing jet, indicating that the jet is launched from the inner x-ray emitting region. The polarization degree is (4.01+-0.20)% at 2 to 8 kiloelectronvolts, implying that the accretion disk is viewed closer to edge-on than the binary orbit. The observations reveal that hot x-ray emitting plasma is spatially extended in a plane perpendicular to the jet axis, not parallel to the jet.


title: Two-sided Robustly Testable Codes

abstract: We show that the tensor product of two random linear codes is robustly testable with high probability. This implies that one can obtain pairs of linear codes such that their product and the product of their dual codes are simultaneously robustly testable. Such two-sided robustly testable codes (with a much weaker form of robustness) were the key ingredient in the recent constructions of asymptotically good quantum LDPC codes, which ensured their linear minimum distance. We hope that the existence of such codes with a stronger form of robustness, shown here, can be used to simplify the proofs and provide better distance bounds in these constructions. We also give new very simple examples of non-robustly testable codes. We show that if the parity-checks of two codes are mutually orthogonal, then their product is not robustly testable. In particular, this implies that the product of a code with its dual can never be robustly testable. We also study a property of a collection of linear codes called product-expansion, which can be viewed as a coboundary expansion of the cochain complex naturally associated with the product of these codes. We show that this property is related with the robust testability and the agreement testability of the products of codes.


title: The Lieb-Oxford Lower Bounds on the Coulomb Energy, Their Importance to Electron Density Functional Theory, and a Conjectured Tight Bound on Exchange

abstract: Lieb and Oxford (1981) derived rigorous lower bounds, in the form of local functionals of the electron density, on the indirect part of the Coulomb repulsion energy. The greatest lower bound for a given electron number N depends monotonically upon N, and the N-> infinity limit is a bound for all N. These bounds have been shown to apply to the exact density functionals for the exchange- and exchange-correlation energies that must be approximated for an accurate and computationally efficient description of atoms, molecules, and solids. A tight bound on the exact exchange energy has been derived therefrom for two-electron ground states, and is conjectured to apply to all spin-unpolarized electronic ground states. Some of these and other exact constraints have been used to construct two generations of non-empirical density functionals beyond the local density approximation: the Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA), and the strongly constrained and appropriately normed (SCAN) meta-GGA.


title: Distributions of gas and small and large grains in the LkH$\alpha\,330$ disk trace a young planetary system

abstract: [abridged] We present new scattered light and millimeter observations of the protoplanetary disk around LkH$\alpha\,330$, using SPHERE/VLT and ALMA, respectively. The scattered-light SPHERE observations reveal an asymmetric ring at around 45au from the star in addition to two spiral arms with similar radial launching points at around 90au. The millimeter observations from ALMA (resolution of 0.06''$\times$0.04'') mainly show an asymmetric ring located at 110au from the star. In addition to this asymmetry, there are two faint symmetric rings at 60au and 200au. The $^{12}$CO, $^{13}$CO, and C$^{18}$O lines seem to be less abundant in the inner disk (these observations have a resolution of 0.16''$\times$0.11''). The $^{13}$CO peaks at a location similar to the inner ring observed with SPHERE, suggesting that this line is optically thick and traces variations of disk temperature instead of gas surface-density variations, while the C$^{18}$O peaks slightly further away at around 60au. We compare our observations with hydrodynamical simulations that include gas and dust evolution, and conclude that a 10$M_{\rm{Jup}}$ mass planet at 60au and in an eccentric orbit ($e=0.1$) can qualitatively explain most of the observed structures. A planet in a circular orbit leads to a much narrower concentration in the millimeter emission, while a planet in a more eccentric orbit leads to a very eccentric cavity as well. In addition, the outer spiral arm launched by the planet changes its pitch angle along the spiral due to the eccentricity and when it interacts with the vortex, potentially appearing in observations as two distinct spirals. Our observations and models show that LkH$\alpha\,330$ is an interesting target to search for (eccentric-) planets while they are still embedded in their parental disk, making it an excellent candidate for studies on planet-disk interaction.


title: Noise Estimation in Gaussian Process Regression

abstract: We develop a computational procedure to estimate the covariance hyperparameters for semiparametric Gaussian process regression models with additive noise. Namely, the presented method can be used to efficiently estimate the variance of the correlated error, and the variance of the noise based on maximizing a marginal likelihood function. Our method involves suitably reducing the dimensionality of the hyperparameter space to simplify the estimation procedure to a univariate root-finding problem. Moreover, we derive bounds and asymptotes of the marginal likelihood function and its derivatives, which are useful to narrowing the initial range of the hyperparameter search. Using numerical examples, we demonstrate the computational advantages and robustness of the presented approach compared to traditional parameter optimization.


title: Thompson Sampling Efficiently Learns to Control Diffusion Processes

abstract: Diffusion processes that evolve according to linear stochastic differential equations are an important family of continuous-time dynamic decision-making models. Optimal policies are well-studied for them, under full certainty about the drift matrices. However, little is known about data-driven control of diffusion processes with uncertain drift matrices as conventional discrete-time analysis techniques are not applicable. In addition, while the task can be viewed as a reinforcement learning problem involving exploration and exploitation trade-off, ensuring system stability is a fundamental component of designing optimal policies. We establish that the popular Thompson sampling algorithm learns optimal actions fast, incurring only a square-root of time regret, and also stabilizes the system in a short time period. To the best of our knowledge, this is the first such result for Thompson sampling in a diffusion process control problem. We validate our theoretical results through empirical simulations with real parameter matrices from two settings of airplane and blood glucose control. Moreover, we observe that Thompson sampling significantly improves (worst-case) regret, compared to the state-of-the-art algorithms, suggesting Thompson sampling explores in a more guarded fashion. Our theoretical analysis involves characterization of a certain optimality manifold that ties the local geometry of the drift parameters to the optimal control of the diffusion process. We expect this technique to be of broader interest.

categories: cs.LG cs.AI cs.SY eess.SY math.DS math.OC


title: German AI Start-Ups and AI Ethics: Using A Social Practice Lens for Assessing and Implementing Socio-Technical Innovation

abstract: Within the current AI ethics discourse, there is a gap in empirical research on understanding how AI practitioners understand ethics and socially organize to operationalize ethical concerns, particularly in the context of AI start-ups. This gap intensifies the risk of a disconnect between scholarly research, innovation, and application. This risk materializes acutely as mounting pressures to identify and mitigate the potential harms of AI systems have created an urgent need to assess and implement socio-technical innovation for fairness, accountability, and transparency. Building on social practice theory, we address this need via a framework that allows AI researchers, practitioners, and regulators to systematically analyze existing cultural understandings, histories, and social practices of ethical AI to define appropriate strategies for effectively implementing socio-technical innovations. Our contributions are threefold: 1) we introduce a practice-based approach for understanding ethical AI; 2) we present empirical findings from our study on the operationalization of ethics in German AI start-ups to underline that AI ethics and social practices must be understood in their specific cultural and historical contexts; and 3) based on our empirical findings, we suggest that ethical AI practices can be broken down into principles, needs, narratives, materializations, and cultural genealogies to form a useful backdrop for considering socio-technical innovations.


title: Mitigating Data Heterogeneity in Federated Learning with Data Augmentation

abstract: Federated Learning (FL) is a prominent framework that enables training a centralized model while securing user privacy by fusing local, decentralized models. In this setting, one major obstacle is data heterogeneity, i.e., each client having non-identically and independently distributed (non-IID) data. This is analogous to the context of Domain Generalization (DG), where each client can be treated as a different domain. However, while many approaches in DG tackle data heterogeneity from the algorithmic perspective, recent evidence suggests that data augmentation can induce equal or greater performance. Motivated by this connection, we present federated versions of popular DG algorithms, and show that by applying appropriate data augmentation, we can mitigate data heterogeneity in the federated setting, and obtain higher accuracy on unseen clients. Equipped with data augmentation, we can achieve state-of-the-art performance using even the most basic Federated Averaging algorithm, with much sparser communication.


title: Semantic preservation for a type directed translation scheme of Featherweight Go

abstract: Featherweight Go (FG) is a minimal core calculus that includes essential Go features such as overloaded methods and interface types. The most straightforward semantic description of the dynamic behavior of FG programs is to resolve method calls based on run-time type information. A more efficient approach is to apply a type-directed translation scheme where interface-values are replaced by dictionaries that contain concrete method definitions. Thus, method calls can be resolved by a simple lookup of the method definition in the dictionary. Establishing that the target program obtained via the type-directed translation scheme preserves the semantics of the original FG program is an important task. To establish this property we employ logical relations that are indexed by types to relate source and target programs. We provide rigorous proofs and give a detailed discussion of the many subtle corners that we have encountered including the need for a step index due to recursive interfaces and method definitions.


title: Measuring Class-Imbalance Sensitivity of Deterministic Performance Evaluation Metrics

abstract: The class-imbalance issue is intrinsic to many real-world machine learning tasks, particularly to the rare-event classification problems. Although the impact and treatment of imbalanced data is widely known, the magnitude of a metric's sensitivity to class imbalance has attracted little attention. As a result, often the sensitive metrics are dismissed while their sensitivity may only be marginal. In this paper, we introduce an intuitive evaluation framework that quantifies metrics' sensitivity to the class imbalance. Moreover, we reveal an interesting fact that there is a logarithmic behavior in metrics' sensitivity meaning that the higher imbalance ratios are associated with the lower sensitivity of metrics. Our framework builds an intuitive understanding of the class-imbalance impact on metrics. We believe this can help avoid many common mistakes, specially the less-emphasized and incorrect assumption that all metrics' quantities are comparable under different class-imbalance ratios.


title: Fast calibration of weak FARIMA models

abstract: In this paper, we investigate the asymptotic properties of Le Cam's one-step estimator for weak Fractionally AutoRegressive Integrated Moving-Average (FARIMA) models. For these models, noises are uncorrelated but neither necessarily independent nor martingale differences errors. We show under some regularity assumptions that the one-step estimator is strongly consistent and asymptotically normal with the same asymptotic variance as the least squares estimator. We show through simulations that the proposed estimator reduces computational time compared with the least squares estimator. An application for providing remotely computed indicators for time series is proposed.


title: Mnemonic: A Parallel Subgraph Matching System for Streaming Graphs

abstract: Finding patterns in large highly connected datasets is critical for value discovery in business development and scientific research. This work focuses on the problem of subgraph matching on streaming graphs, which provides utility in a myriad of real-world applications ranging from social network analysis to cybersecurity. Each application poses a different set of control parameters, including the restrictions for a match, type of data stream, and search granularity. The problem-driven design of existing subgraph matching systems makes them challenging to apply for different problem domains. This paper presents Mnemonic, a programmable system that provides a high-level API and democratizes the development of a wide variety of subgraph matching solutions. Importantly, Mnemonic also delivers key data management capabilities and optimizations to support real-time processing on long-running, high-velocity multi-relational graph streams. The experiments demonstrate the versatility of Mnemonic, as it outperforms several state-of-the-art systems by up to two orders of magnitude.


title: Uncertainty and bias of cosmology and astrophysical population model from statistical dark sirens

abstract: Gravitational-wave (GW) radiation from a coalescing compact binary is a standard siren as the luminosity distance of each event can be directly measured from the amplitude of the signal. One possibility to constrain cosmology using the GW siren is to perform statistical inference on a population of binary black hole (BBH) events. In essence, this statistical method can be viewed as follows. We can modify the shape of the distribution of observed BBH events by changing cosmological parameters until it eventually matches the distribution constructed from an astrophysical population model, thereby allowing us to determine the cosmological parameters. In this work, we derive the Cram\'er-Rao bound for both cosmological parameters and those governing the astrophysical population model from this statistical dark siren method by examining the Fisher information contained in the event distribution. Our study provides analytical insights and enables fast yet accurate estimations of the statistical accuracy of dark siren cosmology. Furthermore, we consider the bias in cosmology due to unmodeled substructures in the merger rate and the mass distribution. We find a $1\%$ deviation in the astrophysical model can lead to a more than $1\%$ error in the Hubble constant. This could limit the accuracy of dark siren cosmology when there are more than $10^4$ BBH events detected.


title: Gauge and Scalar Fields on $\mathbb{CP}^2$: A Gauge-invariant Analysis II. The measure for gauge fields and a 4d WZW theory

abstract: We consider the volume of the gauge orbit space for gauge fields on four-dimensional complex projective space. The analysis uses a parametrization of gauge fields where gauge transformations act homogeneously on the fields, facilitating a manifestly gauge-invariant analysis. The volume element contains a four-dimensional Wess-Zumino-Witten (WZW) action for a hermitian matrix-valued field. There is also a mass-like term for certain components of the gauge field. We discuss how the mass term could be related to results from lattice simulations as well as Schwinger-Dyson equations. We argue for a kinematic regime where the Yang-Mills theory can be approximated by the 4d-WZW theory. The result is suggestive of the instanton liquid picture of QCD. Further it is also indicative of the mechanism for confinement being similar to what happens in two dimensions.


title: Elliptic curves, ACM bundles and Ulrich bundles on prime Fano threefolds

abstract: Let $X$ be any smooth prime Fano threefold of degree $2g-2$ in $\mathbb{P}^{g+1}$, with $g \in \{3,\ldots,10,12\}$. We prove that for any integer $d$ satisfying $\left\lfloor \frac{g+3}{2} \right\rfloor \leq d \leq g+3$ the Hilbert scheme parametrizing smooth irreducible elliptic curves of degree $d$ in $X$ is nonempty and has a component of dimension $d$, which is furthermore reduced except for the case when $(g,d)=(4,3)$ and $X$ is contained in a singular quadric. Consequently, we deduce that the moduli space of rank--two slope--stable $ACM$ bundles $\mathcal{F}_d$ on $X$ such that $\det(\mathcal{F}_d)=\mathcal{O}_X(1)$, $c_2(\mathcal{F}_d)\cdot \mathcal{O}_X(1)=d$ and $h^0(\mathcal{F}_d(-1))=0$ is nonempty and has a component of dimension $2d-g-2$, which is furthermore reduced except for the case when $(g,d)=(4,3)$ and $X$ is contained in a singular quadric. This completes the classification of rank-two $ACM$ bundles on prime Fano threefolds. Secondly, we prove that for every $h \in \mathbb{Z}^+$ the moduli space of stable Ulrich bundles $\mathcal{E}$ of rank $2h$ and determinant $\mathcal{O}_X(3h)$ on $X$ is nonempty and has a reduced component of dimension $h^2(g+3)+1$; this result is optimal in the sense that there are no other Ulrich bundles occurring on $X$. This in particular shows that any prime Fano threefold is Ulrich wild.


title: Measuring Gender Bias in Educational Videos: A Case Study on YouTube

abstract: Students are increasingly using online materials to learn new subjects or to supplement their learning process in educational institutions. Issues regarding gender bias have been raised in the context of formal education and some measures have been proposed to mitigate them. However, online educational materials in terms of possible gender bias and stereotypes which may appear in different forms are yet to be investigated in the context of search bias in a widely-used search platform. As a first step towards measuring possible gender bias in online platforms, we have investigated YouTube educational videos in terms of the perceived gender of their narrators. We adopted bias measures for ranked search results to evaluate educational videos returned by YouTube in response to queries related to STEM (Science, Technology, Engineering, and Mathematics) and NON-STEM fields of education. Gender is a research area by itself in social sciences which is beyond the scope of this work. In this respect, for annotating the perceived gender of the narrator of an instructional video we used only a crude classification of gender into Male, and Female. Then, for analysing perceived gender bias we utilised bias measures that have been inspired by search platforms and further incorporated rank information into our analysis. Our preliminary results demonstrate that there is a significant bias towards the male gender on the returned YouTube educational videos, and the degree of bias varies when we compare STEM and NON-STEM queries. Finally, there is a strong evidence that rank information might affect the results.


title: Improvements to Pan-STARRS1 Astrometry: II. Corrections for Differential Chromatic Refraction

abstract: In a previous paper, we applied the Gaia DR2 catalog to improve the astrometric accuracy of about 1.7 billion objects in Pan-STARRS1 Data Release 2 (PS1 DR2). We report here on further improvements made by utilizing Gaia EDR3 and correcting for the effects of differential chromatic refraction (DCR) in declination. We extend the correction algorithm in Paper 1 by iteratively subtracting color- and declination-dependent PS1/Gaia EDR3 declination residuals. We determine the astrometric improvement for ~440 million reference objects that are point-like and cross-match to Gaia EDR3. For this set of objects, Gaia EDR3 provides a ~3% improvement in PS1 astrometry over Gaia DR2, and DCR corrections provide an additional ~5% improvement. DCR corrections increase substantially for objects observed away from the zenith. DCR corrections lead to an astrometric improvement of ~30% for blue objects (0<g-i<1) that are 50{\deg} away from the zenith. The amplitude of systematic astrometric errors from these effects is substantially reduced to less than 1 mas for objects with PS1 colors in the range 0 < g-i < 4.5, which makes this a useful astrometric reference catalog in fields where there are few Gaia stars. The improved astrometric data will be available through the Mikulski Archive for Space Telescopes PS1 catalog interfaces.


title: Pushed-to-pulled front transitions: continuation, speed scalings, and hidden monotonicity

abstract: We analyze the transition between pulled and pushed fronts both analytically and numerically from a model-independent perspective. Based on minimal conceptual assumptions, we show that pushed fronts bifurcate from a branch of pulled fronts with an effective speed correction that scales quadratically in the bifurcation parameter. Strikingly, we find that in this general context without assumptions on comparison principles, the pulled front loses stability and gives way to a pushed front when monotonicity in the leading edge is lost. Our methods rely on far-field core decompositions that identify explicitly asymptotics in the leading edge of the front. We show how the theoretical construction can be directly implemented to yield effective algorithms that determine spreading speeds and bifurcation points with exponentially small error in the domain size. Example applications considered here include an extended Fisher-KPP equation, a Fisher-Burgers equation, negative taxis in combination with logistic population growth, an autocatalytic reaction, and a Lotka-Volterra model.


title: A detailed analysis of the Gl 486 planetary system

abstract: The Gl 486 system consists of a very nearby, relatively bright, weakly active M3.5 V star at just 8 pc with a warm transiting rocky planet of about 1.3 R_Terra and 3.0 M_Terra that is ideal for both transmission and emission spectroscopy and for testing interior models of telluric planets. To prepare for future studies, we collected light curves of seven new transits observed with the CHEOPS space mission and new radial velocities obtained with MAROON-X/Gemini North and CARMENES/Calar Alto telescopes, together with previously published spectroscopic and photometric data from the two spectrographs and TESS. We also performed interferometric observations with the CHARA Array and new photometric monitoring with a suite of smaller telescopes. From interferometry, we measure a limb-darkened disc angular size of the star Gl 486. Together with a corrected Gaia EDR3 parallax, we obtain a stellar radius. We also measure a stellar rotation period at P_rot ~ 49.9 d, an upper limit to its XUV (5-920 AA) flux with new Hubble/STIS data, and, for the first time, a variety of element abundances (Fe, Mg, Si, V, Sr, Zr, Rb) and C/O ratio. Besides, we impose restrictive constraints on the presence of additional components, either stellar or substellar, in the system. With the input stellar parameters and the radial-velocity and transit data, we determine the radius and mass of the planet Gl 486 b at R_p = 1.343+/0.063 R_Terra and M_p = 3.00+/-0.13 M_Terra. From the planet parameters and the stellar element abundances, we infer the most probable models of planet internal structure and composition, which are consistent with a relatively small metallic core with respect to the Earth, a deep silicate mantle, and a thin volatile upper layer. With all these ingredients, we outline prospects for Gl 486 b atmospheric studies, especially with forthcoming James Webb Space Telescope observations (abridged).


title: Model Optimization in Imbalanced Regression

abstract: Imbalanced domain learning aims to produce accurate models in predicting instances that, though underrepresented, are of utmost importance for the domain. Research in this field has been mainly focused on classification tasks. Comparatively, the number of studies carried out in the context of regression tasks is negligible. One of the main reasons for this is the lack of loss functions capable of focusing on minimizing the errors of extreme (rare) values. Recently, an evaluation metric was introduced: Squared Error Relevance Area (SERA). This metric posits a bigger emphasis on the errors committed at extreme values while also accounting for the performance in the overall target variable domain, thus preventing severe bias. However, its effectiveness as an optimization metric is unknown. In this paper, our goal is to study the impacts of using SERA as an optimization criterion in imbalanced regression tasks. Using gradient boosting algorithms as proof of concept, we perform an experimental study with 36 data sets of different domains and sizes. Results show that models that used SERA as an objective function are practically better than the models produced by their respective standard boosting algorithms at the prediction of extreme values. This confirms that SERA can be embedded as a loss function into optimization-based learning algorithms for imbalanced regression scenarios.


title: Hyperparameter Importance of Quantum Neural Networks Across Small Datasets

abstract: As restricted quantum computers are slowly becoming a reality, the search for meaningful first applications intensifies. In this domain, one of the more investigated approaches is the use of a special type of quantum circuit - a so-called quantum neural network -- to serve as a basis for a machine learning model. Roughly speaking, as the name suggests, a quantum neural network can play a similar role to a neural network. However, specifically for applications in machine learning contexts, very little is known about suitable circuit architectures, or model hyperparameters one should use to achieve good learning performance. In this work, we apply the functional ANOVA framework to quantum neural networks to analyze which of the hyperparameters were most influential for their predictive performance. We analyze one of the most typically used quantum neural network architectures. We then apply this to $7$ open-source datasets from the OpenML-CC18 classification benchmark whose number of features is small enough to fit on quantum hardware with less than $20$ qubits. Three main levels of importance were detected from the ranking of hyperparameters obtained with functional ANOVA. Our experiment both confirmed expected patterns and revealed new insights. For instance, setting well the learning rate is deemed the most critical hyperparameter in terms of marginal contribution on all datasets, whereas the particular choice of entangling gates used is considered the least important except on one dataset. This work introduces new methodologies to study quantum machine learning models and provides new insights toward quantum model selection.


title: Reconstructed pCT Images Using Monte Carlo Simulations of a Scintillating Glass Detector

abstract: The high cost and low image quality traditionally associated with proton computed tomography (pCT) have prevented it from seeing significant use in clinical settings. A cheap, compact, high-density scintillating glass detector capable of being attached to existing proton therapy gantries may help address these concerns. The design of the detector allows for use in conjunction with single-proton counting reconstruction algorithms, as well as beam-based algorithms that do not resolve individual protons within an accelerator bunch. This study presents quantitative reconstructed images of proton stopping power from Monte Carlo generated pCT scans using the radiation transport code MCNP6, demonstrating the feasibility of proton imaging using this detector design. Relative error and contrast have been examined and compared for images reconstructed using two reconstruction algorithms: a standard filtered backprojection algorithm to act as a benchmark, and a variant of a pCT algorithm which utilizes the concept of distance-driven binning.


title: Evolution of domain structure with electron doping in ferroelectric thin films

abstract: To minimize their electrostatic energy, insulating ferroelectric films tend to break up into nanoscale ``Kittel'' domains of opposite polarization that are separated by uncharged 180$^\circ$ domain walls. Here, I report on self-consistent solutions of coupled Landau-Ginzburg-Devonshire and Schr\"odinger equations for an electron-doped ferroelectric thin film. The model is based on LaAlO$_3$/SrTiO$_3$ interfaces in which the SrTiO$_3$ substrate is made ferroelectric by cation substitution or strain. I find that electron doping destabilizes the Kittel domains. As the two-dimensional electron density $n_\mathrm{2D}$ increases, there is a smooth crossover to a zigzag domain wall configuration. The domain wall is positively charged, but is compensated by the electron gas, which attaches itself to the domain wall and screens depolarizing fields. The domain wall approaches a flat head-to-head configuration in the limit of perfect screening. The polarization profile may be manipulated by an external bias voltage and the electron gas may be switched between surfaces of the ferroelectric film.


title: Modelling Populations of Interaction Networks via Distance Metrics

abstract: Network data arises through observation of relational information between a collection of entities. Recent work in the literature has independently considered when (i) one observes a sample of networks, connectome data in neuroscience being a ubiquitous example, and (ii) the units of observation within a network are edges or paths, such as emails between people or a series of page visits to a website by a user, often referred to as interaction network data. The intersection of these two cases, however, is yet to be considered. In this paper, we propose a new Bayesian modelling framework to analyse such data. Given a practitioner-specified distance metric between observations, we define families of models through location and scale parameters, akin to a Gaussian distribution, with subsequent inference of model parameters providing reasoned statistical summaries for this non-standard data structure. To facilitate inference, we propose specialised Markov chain Monte Carlo (MCMC) schemes capable of sampling from doubly-intractable posterior distributions over discrete and multi-dimensional parameter spaces. Through simulation studies we confirm the efficacy of our methodology and inference scheme, whilst its application we illustrate via an example analysis of a location-based social network (LSBN) data set.


title: Martingales on Principal Fiber Bundles

abstract: Let $P(M,G)$ be a principal fiber bundle, let $\omega$ be a connection form on $P(M,G)$, and consider a projectable connection $\nabla^{P}$ on $P(M,G)$. The aim of this work is to determine the $\nabla^{P}$-martingales in $P(M,G)$. Our results allow establishing new characterizations of harmonic maps from Riemannian manifolds to principal fiber bundles.


title: Search for resonant and nonresonant production of pairs of dijet resonances in proton-proton collisions at $\sqrt{s}$ = 13 TeV

abstract: A search for pairs of dijet resonances with the same mass is conducted in final states with at least four jets. Results are presented separately for the case where the four jet production proceeds via an intermediate resonant state and for nonresonant production. The search uses a data sample corresponding to an integrated luminosity of 138 fb$^{-1}$ collected by the CMS detector in proton-proton collisions at $\sqrt{s}$ = 13 TeV. Model-independent limits, at 95% confidence level, are reported on the production cross section of four-jet and dijet resonances. These first LHC limits on resonant pair production of dijet resonances via high mass intermediate states are applied to a signal model of diquarks that decay into pairs of vector-like quarks, excluding diquark masses below 7.6 TeV for a particular model scenario. There are two events in the tails of the distributions, each with a four-jet mass of 8 TeV and an average dijet mass of 2 TeV, resulting in local and global significances of 3.9 and 1.6 standard deviations, respectively, if interpreted as a signal. The nonresonant search excludes pair production of top squarks with masses between 0.50 TeV to 0.77 TeV, with the exception of a small interval between 0.52 and 0.58 TeV, for supersymmetric $R$-parity-violating decays to quark pairs, significantly extending previous limits. Here, the most significant excess above the predicted background occurs at an average dijet mass of 0.95 TeV, for which the local and global significances are 3.6 and 2.5 standard deviations, respectively.


title: Spin Torque Generated by Valley Hall Effect in WSe2

abstract: Monolayer transition metal dichalcogenides are promising materials for spintronics due to their robust spin-valley locked valence states, enabling efficient charge-to-spin conversion via valley Hall effect with non-equilibrium spins possessing long spin diffusion lengths of hundreds of nanometers. In this work, we show that the injection of a pure valley current, induced by valley Hall effect in a WSe2 monolayer, imparts a spin torque on the magnetization of an overlaid Fe or CoFe in a tunneling structure. The torque efficiency is found to be comparable to that in conventional perpendicular magnetic tunnel junctions and can be further optimized with valley Hall angle in WSe2. The valley nature of the spin torque gives rise to out-of-plane damping-like torques in a current-in-plane configuration, vanishing charge transport perpendicular-to-the-plane as well as torque efficiency tunable through gating.


title: Understanding RowHammer Under Reduced Wordline Voltage: An Experimental Study Using Real DRAM Devices

abstract: RowHammer is a circuit-level DRAM vulnerability, where repeatedly activating and precharging a DRAM row, and thus alternating the voltage of a row's wordline between low and high voltage levels, can cause bit flips in physically nearby rows. Recent DRAM chips are more vulnerable to RowHammer: with technology node scaling, the minimum number of activate-precharge cycles to induce a RowHammer bit flip reduces and the RowHammer bit error rate increases. Therefore, it is critical to develop effective and scalable approaches to protect modern DRAM systems against RowHammer. To enable such solutions, it is essential to develop a deeper understanding of the RowHammer vulnerability of modern DRAM chips. However, even though the voltage toggling on a wordline is a key determinant of RowHammer vulnerability, no prior work experimentally demonstrates the effect of wordline voltage (VPP) on the RowHammer vulnerability. Our work closes this gap in understanding. This is the first work to experimentally demonstrate on 272 real DRAM chips that lowering VPP reduces a DRAM chip's RowHammer vulnerability. We show that lowering VPP 1) increases the number of activate-precharge cycles needed to induce a RowHammer bit flip by up to 85.8% with an average of 7.4% across all tested chips and 2) decreases the RowHammer bit error rate by up to 66.9% with an average of 15.2% across all tested chips. At the same time, reducing VPP marginally worsens a DRAM cell's access latency, charge restoration, and data retention time within the guardbands of system-level nominal timing parameters for 208 out of 272 tested chips. We conclude that reducing VPP is a promising strategy for reducing a DRAM chip's RowHammer vulnerability without requiring modifications to DRAM chips.


title: Self Tuned Criticality: Controlling a neuron near its bifurcation point via temporal correlations

abstract: Previous work showed that the collective activity of large neuronal networks can be tamed to remain near its critical point by a feedback control that maximizes the temporal correlations of the mean-field fluctuations. Since such correlations behave similarly near instabilities across nonlinear dynamical systems, it is expected that the principle should control also low dimensional dynamical systems exhibiting continuous or discontinuous bifurcations from fixed points to limit cycles. Here we present numerical evidence that the dynamics of a single neuron can be controlled in the vicinity of its bifurcation point. The approach is tested in two models: a 2D generic excitable map and the paradigmatic FitzHugh-Nagumo neuron model. The results show that in both cases, the system can be self-tuned to its bifurcation point by modifying the control parameter according to the first coefficient of the autocorrelation function.


title: Nodal solutions for the weighted biharmonic equation with critical exponential growth

abstract: In this paper, we deal with the logarithmic weighted fourth order elliptic equation in the unit disk of $B\subset\R^{4}$: $$\displaystyle(P_{\lambda})~~\Delta(w(x) \Delta u) = \lambda\ f(x,u) \quad\mbox{ in }\quad B, \quad u=\frac{\partial u}{\partial n} \quad\mbox{ on } \quad\partial B,$$ where the nonlinearity $f$ is assumed to have exponential critical growth in view of Adam's type inequalities. By using the constrained minimization in Nehari set combined with the quantitative deformation lemma and degree theory, we prove the existence of nodal solutions to the problem $(P_{\lambda})$ .


title: The $\mu$-neutral fractional multi-delayed differential equations

abstract: The $\mu$-neutral linear fractional multi-delayed differential nonhomogeneous system with noncommutative coefficient matrices is introduced. The novel $\mu$-neutral multi-delayed perturbation of Mittag-Leffler type matrix function is proposed. Based on this, an explicit solution to the system is investigated step by step. The existence uniqueness of solutions to $\mu$-neutral nonlinear fractional multi-delayed differential system is obtained with regard to the supremum norm. The notion of stability analysis in the sense of solutions to the described system is discussed on the grounds of the fixed point approach.


title: Folding rotationally symmetrical tableaux via webs

abstract: Rectangular standard Young tableaux with 2 or 3 rows are in bijection with $U_q(\mathfrak{sl}_2)$-webs and $U_q(\mathfrak{sl}_3)$-webs respectively. When $W$ is a web with a reflection symmetry, the corresponding tableau $T_W$ has a rotational symmetry. Folding $T_W$ transforms it into a domino tableau $D_W$. We study the relationships between these correspondences. For 2-row tableaux, folding a rotationally symmetric tableau corresponds to "literally folding" the web along its axis of symmetry. For $3$-row tableaux, we give simple algorithms, which provide direct bijective maps between symmetrical webs and domino tableaux (in both directions). These details of these algorithms reflect the intuitive idea that $D_W$ corresponds to "$W$ modulo symmetry".


title: Quantitative bounds for product of simplices in subsets of the unit cube

abstract: For each $1\leq i \le n$, let $k_i\geq 1$ and let $\Delta_i$ be a set of vertices of a non-degenerate simplex of $k_i+1$ points in $\mathbb{R}^{k_i+1}$. If $A\subseteq [0,1]^{k_1+1}\times \cdots \times [0,1]^{k_n+1}$ is a Lebesgue measurable set of measure at least $\delta$, we show that there exists an interval $I=I(\Delta_1,\ldots, \Delta_n,A)$ of length at least $\exp(-\delta^{-C(\Delta_1,\ldots, \Delta_n)})$ such that for each $\lambda\in I$, the set $A$ contains $\Delta'_1\times \cdots \times \Delta'_n$, where each $\Delta_i'$ is an isometric copy of $\lambda\Delta_i$. This is a quantitative improvement of a result by Lyall and Magyar. Our proof relies on harmonic analysis. The main ingredient in the proof are cancellation estimates for forms similar to multilinear singular integrals associated with $n$-partite $n$-regular hypergraphs.


title: Feasibility study of clinical target volume definition for soft-tissue sarcoma using muscle fiber orientations derived from diffusion tensor imaging

abstract: Objective: Soft-tissue sarcoma spreads preferentially along muscle fibers. We explore the utility of deriving muscle fiber orientations from diffusion tensor MRI (DT-MRI) for defining the boundary of the clinical target volume in muscle tissue. Approach: We recruited eight healthy volunteers to acquire MR images of the left and right thigh. The imaging session consisted of (a) two MRI spin-echo-based scans, T1- and T2-weighted; (b) a diffusion weighted (DW) spin-echo-based scan using an echo planar acquisition with fat suppression. The thigh muscles were auto-segmented using CNN. DT-MRI data was used as a geometry encoding input to solve the anisotropic Eikonal equation with Hamiltonian Fast-Marching method. The isosurfaces of the solution modeled the CTV boundary. Main results: The auto-segmented muscles of the thigh agreed with manually delineated with the Dice score ranging from 0.8 to 0.94 for different muscles. Anisotropy of the isosurfaces was compared across muscles with different anatomical orientations within a thigh, between muscles in left and right thighs of each subject, and between different subjects. Analysis showed a high degree of consistency across all comparisons. The distance from the GTV to the isosurface and the eigenvalues ratio are two controlling parameters for the extent and shape of the CTV. Significance: Our feasibility study with healthy volunteers shows the promise of using muscle fiber orientations derived from diffusion weighted MRI data for automated generation of anisotropic CTV boundary in soft tissue sarcoma. Our contribution is significant as it is expected to lead to the improvements in the treatment outcomes of soft-tissue sarcoma patients undergoing radiotherapy and decrease amputation rate for a subset of patients. We expect such improvements to have a strong positive impact for the cancer centers with small volume of sarcoma patients.


title: Capabilities and Limits of the Unitary Coupled-cluster Approach with Generalized Two-body Cluster Operators

abstract: Unitary cluster expansions of the electronic wavefunction have recently gained much interest because of their use in conjunction with quantum algorithms. In this contribution, we investigate some aspects of an ansatz using generalized two-body excitations operators, which has been considered in some recent works on quantum algorithms for quantum chemistry. Our numerical results show that in particular two-body operators with effective particle-hole excitation level of one in connection with the usual particle-hole double excitation operators lead to a very accurate yet compact representation of the wavefunction. Generalized two-body operators with effective excitation rank zero have a considerably less pronounced effect. We compare to standard and unitary coupled-cluster expansions and show that the above mentioned approach matches or even surpasses the accuracy of expansions with three-body particle-hole excitations, in particular at the onset of strong correlation. A downside of the approach is that it is rather difficult to rigorously converge it to its variational minimum.


title: Building Blocks for Network-Accelerated Distributed File Systems

abstract: High-performance clusters and datacenters pose increasingly demanding requirements on storage systems. If these systems do not operate at scale, applications are doomed to become I/O bound and waste compute cycles. To accelerate the data path to remote storage nodes, remote direct memory access (RDMA) has been embraced by storage systems to let data flow from the network to storage targets, reducing overall latency and CPU utilization. Yet, this approach still involves CPUs on the data path to enforce storage policies such as authentication, replication, and erasure coding. We show how storage policies can be offloaded to fully programmable SmartNICs, without involving host CPUs. By using PsPIN, an open-hardware SmartNIC, we show latency improvements for writes (up to 2x), data replication (up to 2x), and erasure coding (up to 2x), when compared to respective CPU- and RDMA-based alternatives.


title: Watkins's conjecture for quadratic twists of Elliptic Curves with Prime Power Conductor

abstract: Watkins' conjecture asserts that the rank of an elliptic curve is upper bounded by the $2$-adic valuation of its modular degree. We show that this conjecture is satisfied when $E$ is any quadratic twist of an elliptic curve with rational $2$-torsion and prime power conductor. Furthermore, we give a lower bound of the congruence number for elliptic curves of the form $y^2=x^3-dx$, with $d$ a biquadratefree integer.


title: Event-Case Correlation for Process Mining using Probabilistic Optimization

abstract: Process mining supports the analysis of the actual behavior and performance of business processes using event logs. % such as, e.g., sales transactions recorded by an ERP system. An essential requirement is that every event in the log must be associated with a unique case identifier (e.g., the order ID of an order-to-cash process). In reality, however, this case identifier may not always be present, especially when logs are acquired from different systems or extracted from non-process-aware information systems. In such settings, the event log needs to be pre-processed by grouping events into cases -- an operation known as event correlation. Existing techniques for correlating events have worked with assumptions to make the problem tractable: some assume the generative processes to be acyclic, while others require heuristic information or user input. Moreover, %these techniques' primary assumption is that they abstract the log to activities and timestamps, and miss the opportunity to use data attributes. % In this paper, we lift these assumptions and propose a new technique called EC-SA-Data based on probabilistic optimization. The technique takes as inputs a sequence of timestamped events (the log without case IDs), a process model describing the underlying business process, and constraints over the event attributes. Our approach returns an event log in which every event is associated with a case identifier. The technique allows users to incorporate rules on process knowledge and data constraints flexibly. The approach minimizes the misalignment between the generated log and the input process model, maximizes the support of the given data constraints over the correlated log, and the variance between activity durations across cases. Our experiments with various real-life datasets show the advantages of our approach over the state of the art.


title: Extremal graph realizations and graph Laplacian eigenvalues

abstract: For a regular polyhedron (or polygon) centered at the origin, the coordinates of the vertices are eigenvectors of the graph Laplacian for the skeleton of that polyhedron (or polygon) associated with the first (non-trivial) eigenvalue. In this paper, we generalize this relationship. For a given graph, we study the eigenvalue optimization problem of maximizing the first (non-trivial) eigenvalue of the graph Laplacian over non-negative edge weights. We show that the spectral realization of the graph using the eigenvectors corresponding to the solution of this problem, under certain assumptions, is a centered, unit-distance graph realization that has maximal total variance. This result gives a new method for generating unit-distance graph realizations and is based on convex duality. A drawback of this method is that the dimension of the realization is given by the multiplicity of the extremal eigenvalue, which is typically unknown prior to solving the eigenvalue optimization problem. Our results are illustrated with a number of examples.


title: When Does Re-initialization Work?

abstract: Re-initializing a neural network during training has been observed to improve generalization in recent works. Yet it is neither widely adopted in deep learning practice nor is it often used in state-of-the-art training protocols. This raises the question of when re-initialization works, and whether it should be used together with regularization techniques such as data augmentation, weight decay and learning rate schedules. In this work, we conduct an extensive empirical comparison of standard training with a selection of re-initialization methods to answer this question, training over 15,000 models on a variety of image classification benchmarks. We first establish that such methods are consistently beneficial for generalization in the absence of any other regularization. However, when deployed alongside other carefully tuned regularization techniques, re-initialization methods offer little to no added benefit for generalization, although optimal generalization performance becomes less sensitive to the choice of learning rate and weight decay hyperparameters. To investigate the impact of re-initialization methods on noisy data, we also consider learning under label noise. Surprisingly, in this case, re-initialization significantly improves upon standard training, even in the presence of other carefully tuned regularization techniques.


title: Limitations of the NTK for Understanding Generalization in Deep Learning

abstract: The ``Neural Tangent Kernel'' (NTK) (Jacot et al 2018), and its empirical variants have been proposed as a proxy to capture certain behaviors of real neural networks. In this work, we study NTKs through the lens of scaling laws, and demonstrate that they fall short of explaining important aspects of neural network generalization. In particular, we demonstrate realistic settings where finite-width neural networks have significantly better data scaling exponents as compared to their corresponding empirical and infinite NTKs at initialization. This reveals a more fundamental difference between the real networks and NTKs, beyond just a few percentage points of test accuracy. Further, we show that even if the empirical NTK is allowed to be pre-trained on a constant number of samples, the kernel scaling does not catch up to the neural network scaling. Finally, we show that the empirical NTK continues to evolve throughout most of the training, in contrast with prior work which suggests that it stabilizes after a few epochs of training. Altogether, our work establishes concrete limitations of the NTK approach in understanding generalization of real networks on natural datasets.


title: Measuring the Effect of Training Data on Deep Learning Predictions via Randomized Experiments

abstract: We develop a new, principled algorithm for estimating the contribution of training data points to the behavior of a deep learning model, such as a specific prediction it makes. Our algorithm estimates the AME, a quantity that measures the expected (average) marginal effect of adding a data point to a subset of the training data, sampled from a given distribution. When subsets are sampled from the uniform distribution, the AME reduces to the well-known Shapley value. Our approach is inspired by causal inference and randomized experiments: we sample different subsets of the training data to train multiple submodels, and evaluate each submodel's behavior. We then use a LASSO regression to jointly estimate the AME of each data point, based on the subset compositions. Under sparsity assumptions ($k \ll N$ datapoints have large AME), our estimator requires only $O(k\log N)$ randomized submodel trainings, improving upon the best prior Shapley value estimators.


title: Deep Partial Least Squares for Empirical Asset Pricing

abstract: We use deep partial least squares (DPLS) to estimate an asset pricing model for individual stock returns that exploits conditioning information in a flexible and dynamic way while attributing excess returns to a small set of statistical risk factors. The novel contribution is to resolve the non-linear factor structure, thus advancing the current paradigm of deep learning in empirical asset pricing which uses linear stochastic discount factors under an assumption of Gaussian asset returns and factors. This non-linear factor structure is extracted by using projected least squares to jointly project firm characteristics and asset returns on to a subspace of latent factors and using deep learning to learn the non-linear map from the factor loadings to the asset returns. The result of capturing this non-linear risk factor structure is to characterize anomalies in asset returns by both linear risk factor exposure and interaction effects. Thus the well known ability of deep learning to capture outliers, shed lights on the role of convexity and higher order terms in the latent factor structure on the factor risk premia. On the empirical side, we implement our DPLS factor models and exhibit superior performance to LASSO and plain vanilla deep learning models. Furthermore, our network training times are significantly reduced due to the more parsimonious architecture of DPLS. Specifically, using 3290 assets in the Russell 1000 index over a period of December 1989 to January 2018, we assess our DPLS factor model and generate information ratios that are approximately 1.2x greater than deep learning. DPLS explains variation and pricing errors and identifies the most prominent latent factors and firm characteristics.

categories: q-fin.PR cs.LG q-fin.PM q-fin.RM stat.ML


title: An interval-valued recursive estimation framework for linearly parameterized systems

abstract: This paper proposes a recursive interval-valued estimation framework for identifying the parameters of linearly parameterized systems which may be slowly time-varying. It is assumed that the model error (which may consist in measurement noise or model mismatch or both) is unknown but lies at each time instant in a known interval. In this context, the proposed method relies on bounding the error generated by a given reference point-valued recursive estimator, for example, the well-known recursive least squares algorithm. We discuss the trade-off between computational complexity and tightness of the estimated parametric interval.


title: High-accuracy Rb$_{2}^+$ interaction potentials based on coupled cluster calculations

abstract: This work discusses a protocol for constructing highly accurate potential energy curves (PECs) for the lowest two states of Rb$_{2}^+$, i.e. $X\,{}^2{\Sigma}{_g^+}$ and $(1) {}^2\Sigma{_u^+}$, using an additivity scheme based on coupled-cluster theory. The approach exploits the findings of our previous work [J. Schnabel, L. Cheng and A. K\"ohn, J. Chem. Phys. 155, 124101 (2021)] to avoid the unphysical repulsive long-range barrier occurring for symmetric molecular ions when perturbative estimates of higher-order cluster operators are employed. Furthermore, care was taken to reproduce the physically correct exchange splitting of the $X {}^2{\Sigma}{_g^+}$ and $(1) {}^2{\Sigma}{_u^+}$ PECs. The accuracy of our computational approach is benchmarked for ionization energies of Rb and for spectroscopic constants as well as vibrational levels of the $a {}^3{\Sigma}{_u^+}$ triplet state of Rb\textsubscript{2}. We study high-level correlation contributions, high-level relativistic effects and inner-shell correlation contributions and find very good agreement with experimental reference values for the atomic ionization potential and the binding energy of Rb$_{2}$ in the $a\,{}^3{\Sigma}{_u^+}$ triplet state. Our final best estimate for the binding energy of the Rb$_{2}^+$ $X {}^2{\Sigma}{_g^+}$ state including zero-point vibrational contributions is $D_0 = 6179\,\mathrm{cm}^{-1}$ with an estimated error bound of $\mathcal{O}(\pm 30\,\mathrm{cm}^{-1})$. This value is smaller than the experimentally inferred lower bond of $D_0\ge 6307.5\,\mathrm{cm}^{-1}$ [Bellos et al., Phys. Rev. A 87, 012508 (2013)] and will require further investigation. For the $(1) {}^2{\Sigma}{_u^+}$ state a shallow potential with $D_0 = 78.4\,\mathrm{cm}^{-1}$ and an error bound of $\pm 9\,\mathrm{cm}^{-1}$ is computed.


title: Pattern bounds for principal specializations of $\beta$-Grothendieck Polynomials

abstract: There has been recent interest in lower bounds for the principal specializations of Schubert polynomials $\nu_w := \mathfrak S_w(1,\dots,1)$. We prove a conjecture of Yibo Gao in the setting of $1243$-avoiding permutations that gives a lower bound for $\nu_w$ in terms of the permutation patterns contained in $w$. We extended this result to principal specializations of $\beta$-Grothendieck polynomials $\nu^{(\beta)}_w := \mathfrak G^{(\beta)}_w(1,\dots,1)$ by restricting to the class of vexillary $1243$-avoiding permutations. Our methods are bijective, offering a combinatorial interpretation of the coefficients $c_w$ and $c^{(\beta)}_w$ appearing in these conjectures.


title: Propagation of chaos for maxima of particle systems with mean-field drift interaction

abstract: We study the asymptotic behavior of the normalized maxima of real-valued diffusive particles with mean-field drift interaction. Our main result establishes propagation of chaos: in the large population limit, the normalized maxima behave as those arising in an i.i.d. system where each particle follows the associated McKean--Vlasov limiting dynamics. Because the maximum depends on all particles, our result does not follow from classical propagation of chaos, where convergence to an i.i.d. limit holds for any fixed number of particles but not all particles simultaneously. The proof uses a change of measure argument that depends on a delicate combinatorial analysis of the iterated stochastic integrals appearing in the chaos expansion of the Radon--Nikodym density.


title: flow-based clustering and spectral clustering: a comparison

abstract: We propose and study a novel graph clustering method for data with an intrinsic network structure. Similar to spectral clustering, we exploit an intrinsic network structure of data to construct Euclidean feature vectors. These feature vectors can then be fed into basic clustering methods such as k-means or Gaussian mixture model (GMM) based soft clustering. What sets our approach apart from spectral clustering is that we do not use the eigenvectors of a graph Laplacian to construct the feature vectors. Instead, we use the solutions of total variation minimization problems to construct feature vectors that reflect connectivity between data points. Our motivation is that the solutions of total variation minimization are piece-wise constant around a given set of seed nodes. These seed nodes can be obtained from domain knowledge or by simple heuristics that are based on the network structure of data. Our results indicate that our clustering methods can cope with certain graph structures that are challenging for spectral clustering methods.


title: On the nodal structures of random fields -- a decade of results

abstract: We survey a decade worth of work pertaining to the nodal structures of random fields, with emphasis on the transformative techniques that shaped the field.


title: Hidden-nucleons neural-network quantum states for the nuclear many-body problem

abstract: We generalize the hidden-fermion family of neural network quantum states to encompass both continuous and discrete degrees of freedom and solve the nuclear many-body Schr\"odinger equation in a systematically improvable fashion. We demonstrate that adding hidden nucleons to the original Hilbert space considerably augments the expressivity of the neural-network architecture compared to the Slater-Jastrow ansatz. The benefits of explicitly encoding in the wave function point symmetries such as parity and time-reversal are also discussed. Leveraging on improved optimization methods and sampling techniques, the hidden-nucleon ansatz achieves an accuracy comparable to the numerically-exact hyperspherical harmonic method in light nuclei and to the auxiliary field diffusion Monte Carlo in $^{16}$O. Thanks to its polynomial scaling with the number of nucleons, this method opens the way to highly-accurate quantum Monte Carlo studies of medium-mass nuclei.


title: Stochastic Online Learning with Feedback Graphs: Finite-Time and Asymptotic Optimality

abstract: We revisit the problem of stochastic online learning with feedback graphs, with the goal of devising algorithms that are optimal, up to constants, both asymptotically and in finite time. We show that, surprisingly, the notion of optimal finite-time regret is not a uniquely defined property in this context and that, in general, it is decoupled from the asymptotic rate. We discuss alternative choices and propose a notion of finite-time optimality that we argue is \emph{meaningful}. For that notion, we give an algorithm that admits quasi-optimal regret both in finite-time and asymptotically.


title: Direct link between disorder, mobility and magnetoresistance in topological semimetals

abstract: The extent to which disorder influences the properties of topological semimetals remains an open question and is relevant to both the understanding of topological states and the use of topological materials in practical applications. Here, we achieve unmatched and systematic control of point defect concentrations in the prototypical Dirac semimetal Cd$_3$As$_2$ to gain important insight into the role of disorder on electron transport behavior. We find that arsenic vacancies introduce localized states near the Fermi level and strongly influence the electron mobility. Reducing arsenic vacancies by changing the As/Cd flux ratio used during deposition results in an increase in the magnetoresistance from 200%-1000% and an increase in mobility from 5000-18,000 cm$^2$/Vs. However, the degree of linear magnetoresistance, which has previously been linked to disorder, is found here to correlate inversely with measures of disorder, including disorder potential and disorder correlation lengths. This finding yields important new information in the quest to identify the origin of linear magnetoresistance in a wider range of materials.


title: Magnetic ground states of honeycomb lattice Wigner crystals

abstract: Lattice Wigner crystal states stabilized by long-range Coulomb interactions have recently been realized in two-dimensional moir\'e materials. We employ large-scale unrestricted Hartree-Fock techniques to unveil the magnetic phase diagrams of honeycomb lattice Wigner crystals. For the three lattice filling factors with the largest charge gaps, n = 2/3, 1/2, 1/3, the magnetic phase diagrams contain multiple phases, including ones with non-collinear and non-coplanar spin arrangements. We discuss magnetization evolution with external magnetic field, which has potential as an experimental signature of exotic spin states. Our theoretical results could potentially be validated in moir\'e materials formed from group VI transition metal dichalcogenide twisted homobilayers.


title: Learning from Positive and Negative Examples: New Proof for Binary Alphabets

abstract: One of the most fundamental problems in computational learning theory is the the problem of learning a finite automaton $A$ consistent with a finite set $P$ of positive examples and with a finite set $N$ of negative examples. By consistency, we mean that $A$ accepts all strings in $P$ and rejects all strings in $N$. It is well known that this problem is NP-complete. In the literature, it is stated that this NP-hardness holds even in the case of a binary alphabet. As a standard reference for this theorem, the work of Gold from 1978 is either cited or adapted. But as a crucial detail, the work of Gold actually considered Mealy machines and not deterministic finite state automata (DFAs) as they are considered nowadays. As Mealy automata are equipped with an output function, they can be more compact than DFAs which accept the same language. We show that the adaptions of Gold's construction for Mealy machines stated in the literature have some issues and give a new construction for DFAs with a binary alphabet ourselves.


title: Large-scale Hydrodynamical Shocks as the Smoking Gun Evidence for a Bar in M31

abstract: The formation and evolutionary history of M31 are closely related to its dynamical structures, which remain unclear due to its high inclination. Gas kinematics could provide crucial evidence for the existence of a rotating bar in M31. Using the position-velocity diagram of [OIII] and HI, we are able to identify clear sharp velocity jump (shock) features with a typical amplitude over 100 km/s in the central region of M31 (4.6 kpc X 2.3 kpc, or 20 arcmin X 10 arcmin). We also simulate gas morphology and kinematics in barred M31 potentials and find that the bar-induced shocks can produce velocity jumps similar to those in [OIII]. The identified shock features in both [OIII] and HI are broadly consistent, and they are found mainly on the leading sides of the bar/bulge, following a hallmark pattern expected from the bar-driven gas inflow. Shock features on the far side of the disk are clearer than those on the near side, possibly due to limited data coverage on the near side, as well as obscuration by the warped gas and dust layers. Further hydrodynamical simulations with more sophisticated physics are desired to fully understand the observed gas features and to better constrain the parameters of the bar in M31.


title: DNA: Proximal Policy Optimization with a Dual Network Architecture

abstract: This paper explores the problem of simultaneously learning a value function and policy in deep actor-critic reinforcement learning models. We find that the common practice of learning these functions jointly is sub-optimal, due to an order-of-magnitude difference in noise levels between these two tasks. Instead, we show that learning these tasks independently, but with a constrained distillation phase, significantly improves performance. Furthermore, we find that the policy gradient noise levels can be decreased by using a lower \textit{variance} return estimate. Whereas, the value learning noise level decreases with a lower \textit{bias} estimate. Together these insights inform an extension to Proximal Policy Optimization we call \textit{Dual Network Architecture} (DNA), which significantly outperforms its predecessor. DNA also exceeds the performance of the popular Rainbow DQN algorithm on four of the five environments tested, even under more difficult stochastic control settings.


title: Intention-Aware Navigation in Crowds with Extended-Space POMDP Planning

abstract: This paper presents a hybrid online Partially Observable Markov Decision Process (POMDP) planning system that addresses the problem of autonomous navigation in the presence of multi-modal uncertainty introduced by other agents in the environment. As a particular example, we consider the problem of autonomous navigation in dense crowds of pedestrians and among obstacles. Popular approaches to this problem first generate a path using a complete planner (e.g., Hybrid A*) with ad-hoc assumptions about uncertainty, then use online tree-based POMDP solvers to reason about uncertainty with control over a limited aspect of the problem (i.e. speed along the path). We present a more capable and responsive real-time approach enabling the POMDP planner to control more degrees of freedom (e.g., both speed AND heading) to achieve more flexible and efficient solutions. This modification greatly extends the region of the state space that the POMDP planner must reason over, significantly increasing the importance of finding effective roll-out policies within the limited computational budget that real time control affords. Our key insight is to use multi-query motion planning techniques (e.g., Probabilistic Roadmaps or Fast Marching Method) as priors for rapidly generating efficient roll-out policies for every state that the POMDP planning tree might reach during its limited horizon search. Our proposed approach generates trajectories that are safe and significantly more efficient than the previous approach, even in densely crowded dynamic environments with long planning horizons.


title: SynWMD: Syntax-aware Word Mover's Distance for Sentence Similarity Evaluation

abstract: Word Mover's Distance (WMD) computes the distance between words and models text similarity with the moving cost between words in two text sequences. Yet, it does not offer good performance in sentence similarity evaluation since it does not incorporate word importance and fails to take inherent contextual and structural information in a sentence into account. An improved WMD method using the syntactic parse tree, called Syntax-aware Word Mover's Distance (SynWMD), is proposed to address these two shortcomings in this work. First, a weighted graph is built upon the word co-occurrence statistics extracted from the syntactic parse trees of sentences. The importance of each word is inferred from graph connectivities. Second, the local syntactic parsing structure of words is considered in computing the distance between words. To demonstrate the effectiveness of the proposed SynWMD, we conduct experiments on 6 textual semantic similarity (STS) datasets and 4 sentence classification datasets. Experimental results show that SynWMD achieves state-of-the-art performance on STS tasks. It also outperforms other WMD-based methods on sentence classification tasks.


title: High-order harmonic generation in three-dimensional Weyl semimetals

abstract: In this paper, the nonlinear interaction of Weyl semimetal (WSM) with a strong driving electromagnetic wave-field is investigated. In the scope of the structure-gauge invariant low-energy nonlinear electrodynamic theory, the polarization-resolved high-order harmonic generation spectra in WSM are analyzed. The obtained results show that the spectra in WSM are completely different compared to 2D graphene case. In particular, at the non-collinear arrangement of the electric and Weyl nodes' momentum separation vectors, an anomalous harmonics are generated which are polarised perpendicular to the pump wave electric field. The intensities of anomalous harmonics are quadratically dependent on the momentum space separation of the Weyl nodes. If the right and the left Weyl fermions are merged, we have a 4-component trivial massless Dirac fermion and, as a consequence, the anomalous harmonics vanish. In contrast to the anomalous harmonics, the intensities of normal harmonics do not depend on the Weyl nodes' momentum separation vector, and the harmonics spectra resemble the picture for a massless 3D Dirac fermion.


title: Semisimple Field Theories Detect Stable Diffeomorphism

abstract: Extending the work of the first author, we introduce a notion of semisimple topological field theory in arbitrary even dimension and show that such field theories necessarily lead to stable diffeomorphism invariants. The main result of this paper is a proof that this 'upper bound' is optimal: To this end, we introduce and study a class of semisimple topological field theories, generalizing the well known finite gauge theories constructed by Dijkgraaf-Witten, Freed and Quinn. We show that manifolds satisfying a certain finiteness condition -- including 4-manifolds with finite fundamental group -- are indistinguishable to these field theories if and only if they are stably diffeomorphic. Hence, such generalized Dijkgraaf-Witten theories provide the strongest semisimple TFT invariants possible. These results hold for a large class of ambient tangential structures. We discuss a number of applications, including the constructions of unoriented 4-dimensional semisimple field theories which can distinguish unoriented smooth structure and oriented higher-dimensional semisimple field theories which can distinguish certain exotic spheres. Along the way, we show that dimensional reductions of generalized Dijkgraaf-Witten theories are again generalized Dijkgraaf-Witten theories, we utilize ambidexterity in the rational setting, and we develop techniques related to the $\infty$-categorical Moebius inversion principle of Galvez-Carrillo--Kock--Tonks.


title: Communication-Efficient Federated Learning With Data and Client Heterogeneity

abstract: Federated Learning (FL) enables large-scale distributed training of machine learning models, while still allowing individual nodes to maintain data locally. However, executing FL at scale comes with inherent practical challenges: 1) heterogeneity of the local node data distributions, 2) heterogeneity of node computational speeds (asynchrony), but also 3) constraints in the amount of communication between the clients and the server. In this work, we present the first variant of the classic federated averaging (FedAvg) algorithm which, at the same time, supports data heterogeneity, partial client asynchrony, and communication compression. Our algorithm comes with a rigorous analysis showing that, in spite of these system relaxations, it can provide similar convergence to FedAvg in interesting parameter regimes. Experimental results in the rigorous LEAF benchmark on setups of up to $300$ nodes show that our algorithm ensures fast convergence for standard federated tasks, improving upon prior quantized and asynchronous approaches.


title: Test Time Transform Prediction for Open Set Histopathological Image Recognition

abstract: Tissue typology annotation in Whole Slide histological images is a complex and tedious, yet necessary task for the development of computational pathology models. We propose to address this problem by applying Open Set Recognition techniques to the task of jointly classifying tissue that belongs to a set of annotated classes, e.g. clinically relevant tissue categories, while rejecting in test time Open Set samples, i.e. images that belong to categories not present in the training set. To this end, we introduce a new approach for Open Set histopathological image recognition based on training a model to accurately identify image categories and simultaneously predict which data augmentation transform has been applied. In test time, we measure model confidence in predicting this transform, which we expect to be lower for images in the Open Set. We carry out comprehensive experiments in the context of colorectal cancer assessment from histological images, which provide evidence on the strengths of our approach to automatically identify samples from unknown categories. Code is released at https://github.com/agaldran/t3po .


title: Deep Learning Models on CPUs: A Methodology for Efficient Training

abstract: GPUs have been favored for training deep learning models due to their highly parallelized architecture. As a result, most studies on training optimization focus on GPUs. There is often a trade-off, however, between cost and efficiency when deciding on how to choose the proper hardware for training. In particular, CPU servers can be beneficial if training on CPUs was more efficient, as they incur fewer hardware update costs and better utilizing existing infrastructure. This paper makes several contributions to research on training deep learning models using CPUs. First, it presents a method for optimizing the training of deep learning models on Intel CPUs and a toolkit called ProfileDNN, which we developed to improve performance profiling. Second, we describe a generic training optimization method that guides our workflow and explores several case studies where we identified performance issues and then optimized the Intel Extension for PyTorch, resulting in an overall 2x training performance increase for the RetinaNet-ResNext50 model. Third, we show how to leverage the visualization capabilities of ProfileDNN, which enabled us to pinpoint bottlenecks and create a custom focal loss kernel that was two times faster than the official reference PyTorch implementation.


title: On deformations of isolated singularity functions

abstract: We study multi-parameters deformations of isolated singularity function-germs on either a subanalytic set or a complex analytic spaces. We prove that if such a deformation has no coalescing of singular points, then it has constant topological type. This extends some classical results due to L\^e \& Ramanujam (1976) and Parusi\'nski (1999), as well as a recent result due to Jesus-Almeida and the first author. It also provides a sufficient condition for a one-parameter family of complex isolated singularity surfaces in $\C^3$ to have constant topological type. On the other hand, for complex isolated singularity families defined on an isolated determinantal singularity, we prove that $\mu$-constancy implies constant topological type.


title: Assessment of the Center of Inertia and Regional Inertia with Load Contribution via a Fully Data-Driven Method

abstract: This paper proposes a new comprehensive and fully data-driven methodology to estimate the center of inertia (COI) and the regional inertia, considering the displacement of the COI due to disturbances and load inertial contributions. The strategy uses the typicality-based data analysis (TDA) technique to detect the right pilot-bus that represents the COI. In the TDA, a compound of correlation and cosine similarities is implemented to approximate the actual distribution of the data and find the point (bus) closest to the mean which is elected as the pilot-bus. Then, the frequency response at the pilot-bus and the active power deviations are embedded into an autoregressive moving average exogenous input (ARMAX)-based approach to determine the regional inertia represented by an equivalent machine, whose inertia constant corresponds to the inertial contribution in the Region. % the This is done in order to perform the inertial estimation of the region using the frequency response of the pilot-bus and the active power of the region interconnections to model the region as and equivalent machine whose inertial constant is estimated using ARMAX identification technique. The methodology is tested using the IEEE 68-bus benchmark test system and an adapted version with aggregated dynamical loads, corroborating the method effectiveness.


title: Recent progress on bounds for sets with no three terms in arithmetic progression

abstract: This is the text accompanying my Bourbaki seminar on the work of Bloom and Sisask, Croot, Lev, and Pach, and Ellenberg and Gijswijt.


title: Understanding a Robot's Guiding Ethical Principles via Automatically Generated Explanations

abstract: The continued development of robots has enabled their wider usage in human surroundings. Robots are more trusted to make increasingly important decisions with potentially critical outcomes. Therefore, it is essential to consider the ethical principles under which robots operate. In this paper we examine how contrastive and non-contrastive explanations can be used in understanding the ethics of robot action plans. We build upon an existing ethical framework to allow users to make suggestions about plans and receive automatically generated contrastive explanations. Results of a user study indicate that the generated explanations help humans to understand the ethical principles that underlie a robot's plan.


title: Slowly cooling white dwarfs in NGC 6752

abstract: Recently, a new class of white dwarfs (``slowly cooling WDs'') has been identified in the globular cluster M13. The cooling time of these stars is increased by stable thermonuclear hydrogen burning in their residual envelope. These WDs are thought to be originated by horizontal branch (HB) stars populating the HB blue tail, which skipped the asymptotic giant branch phase. To further explore this phenomenon, we took advantage of deep photometric data acquired with the Hubble Space Telescope in the near-ultraviolet and investigate the bright portion of the WD cooling sequence in NGC 6752, another Galactic globular cluster with metallicity, age and HB morphology similar to M13. The normalized WD luminosity function derived in NGC 6752 turns out to be impressively similar to that observed in M13, in agreement with the fact that the stellar mass distribution along the HB of these two systems is almost identical. As in the case of M13, the comparison with theoretical predictions is consistent with $\sim 70\%$ of the investigated WDs evolving at slower rates than standard, purely cooling WDs. Thanks to its relatively short distance from Earth, NGC 6752 photometry reaches a luminosity one order of a magnitude fainter than the case of M13, allowing us to sample a regime where the cooling time delay, with respect to standard WD models, reaches $\sim 300$ Myr. The results presented in this paper provide new evidence for the existence of slowly cooling WDs and further support to the scenario proposing a direct causal connection between this phenomenon and the horizontal branch morphology of the host stellar cluster.


title: Arnold Tongues in Area-Preserving Maps

abstract: In the early 60's J. B. Keller and D. Levy discovered a fundamental property: the instability tongues in Mathieu-type equations lose sharpness with the addition of higher-frequency harmonics in the Mathieu potentials. 20 years later V. Arnold rediscovered a similar phenomenon on sharpness of Arnold tongues in circle maps (and rediscovered the result of Keller and Levy). In this paper we find a third class of objects where a similarly flavored behavior takes place: area-preserving maps of the cylinder. Speaking loosely, we show that periodic orbits of standard maps are extra fragile with respect to added drift (i.e. non-exactness) if the potential of the map is a trigonometric polynomial. That is, higher-frequency harmonics make periodic orbits more robust with respect to ``drift". The observation was motivated by the study of traveling waves in the discretized sine-Gordon equation which in turn models a wide variety of physical systems.


title: MPA: MultiPath++ Based Architecture for Motion Prediction

abstract: Autonomous driving technology is developing rapidly and nowadays first autonomous rides are being provided in city areas. This requires the highest standards for the safety and reliability of the technology. Motion prediction part of the general self-driving pipeline plays a crucial role in providing these qualities. In this work we present one of the solutions for Waymo Motion Prediction Challenge 2022 based on MultiPath++ ranked the 3rd as of May, 26 2022. Our source code is publicly available on GitHub.


title: Consistent circuits for indefinite causal order

abstract: Over the past decade, a number of quantum processes have been proposed which are logically consistent, yet feature a cyclic causal structure. However, there is no general formal method to construct a process with an exotic causal structure in a way that ensures, and makes clear why, it is consistent. Here we provide such a method, given by an extended circuit formalism. This only requires directed graphs endowed with Boolean matrices, which encode basic constraints on operations. Our framework (a) defines a set of elementary rules for checking the validity of any such graph, (b) provides a way of constructing consistent processes as a circuit from valid graphs, and (c) yields an intuitive interpretation of the causal relations within a process and an explanation of why they do not lead to inconsistencies. We display how several standard examples of exotic processes, including ones that violate causal inequalities, are among the class of processes that can be generated in this way; we conjecture that this class in fact includes all unitarily extendible processes.


title: Classification Utility, Fairness, and Compactness via Tunable Information Bottleneck and R\'enyi Measures

abstract: Designing machine learning algorithms that are accurate yet fair, not discriminating based on any sensitive attribute, is of paramount importance for society to accept AI for critical applications. In this article, we propose a novel fair representation learning method termed the R\'enyi Fair Information Bottleneck Method (RFIB) which incorporates constraints for utility, fairness, and compactness (compression) of representation, and apply it to image and tabular data classification. A key attribute of our approach is that we consider - in contrast to most prior work - both demographic parity and equalized odds as fairness constraints, allowing for a more nuanced satisfaction of both criteria. Leveraging a variational approach, we show that our objectives yield a loss function involving classical Information Bottleneck (IB) measures and establish an upper bound in terms of two R\'enyi measures of order $\alpha$ on the mutual information IB term measuring compactness between the input and its encoded embedding. We study the influence of the $\alpha$ parameter as well as two other tunable IB parameters on achieving utility/fairness trade-off goals, and show that the $\alpha$ parameter gives an additional degree of freedom that can be used to control the compactness of the representation. Experimenting on three different image datasets (EyePACS, CelebA, and FairFace) and two tabular datasets (Adult and COMPAS), using both binary and categorical sensitive attributes, we show that on various utility, fairness, and compound utility/fairness metrics RFIB outperforms current state-of-the-art approaches.


title: Identifiability of deep generative models without auxiliary information

abstract: We prove identifiability of a broad class of deep latent variable models that (a) have universal approximation capabilities and (b) are the decoders of variational autoencoders that are commonly used in practice. Unlike existing work, our analysis does not require weak supervision, auxiliary information, or conditioning in the latent space. Specifically, we show that for a broad class of generative (i.e. unsupervised) models with universal approximation capabilities, the side information $u$ is not necessary: We prove identifiability of the entire generative model where we do not observe $u$ and only observe the data $x$. The models we consider match autoencoder architectures used in practice that leverage mixture priors in the latent space and ReLU/leaky-ReLU activations in the encoder, such as VaDE and MFC-VAE. Our main result is an identifiability hierarchy that significantly generalizes previous work and exposes how different assumptions lead to different "strengths" of identifiability, and includes certain "vanilla" VAEs with isotropic Gaussian priors as a special case. For example, our weakest result establishes (unsupervised) identifiability up to an affine transformation, and thus partially resolves an open problem regarding model identifiability raised in prior work. These theoretical results are augmented with experiments on both simulated and real data.


title: Detection of quasi-periodic micro-structure in three millisecond pulsars with the Large European Array for Pulsars

abstract: We report on the detection of quasi-periodic micro-structure in three millisecond pulsars (MSPs), PSRs J1022+1001, J2145-0750 and J1744-1134, using high time resolution data acquired with the Large European Array for Pulsars at a radio frequency of 1.4 GHz. The occurrence rate of quasi-periodic micro-structure is consistent among pulses with different peak flux densities. Using an auto-correlation analysis, we measure the periodicity and width of the micro-structure in these three pulsars. The detected micro-structure from PSRs J1022+1001 and J1744-1134 is often highly linearly polarised. In PSR J1022+1001, the linear polarisation position angles of micro-structure pulses are in general flat with a small degree of variation. Using these results, we further examine the frequency and rotational period dependency of micro-structure properties established in previous work, along with the angular beaming and temporal modulation models that explains the appearance of micro-structure. We also discuss a possible link of micro-structure to the properties of some of the recently discovered fast radio bursts which exhibit a very similar emission morphology.


title: Impediments to diffusion in quantum graphs: geometry-based upper bounds on the spectral gap

abstract: We derive several upper bounds on the spectral gap of the Laplacian with standard or Dirichlet vertex conditions on compact metric graphs. In particular, we obtain estimates based on the length of a shortest cycle (girth), diameter, total length of the graph, as well as further metric quantities introduced here for the first time, such as the avoidance diameter. Using known results about Ramanujan graphs, a class of expander graphs, we also prove that some of these metric quantities, or combinations thereof, do not to deliver any spectral bounds with the correct scaling.


title: Measure-Theoretically Mixing Subshifts of Minimal Word Complexity

abstract: We resolve a long-standing open question on the relationship between measure-theoretic dynamical complexity and symbolic complexity by establishing the exact word complexity at which measure-theoretic strong mixing manifests: For every superlinear $f : \mathbb{N} \to \mathbb{N}$, i.e. $f(q)/q \to \infty$, there exists a subshift admitting a (strongly) mixing of all orders probability measure with word complexity $p$ such that $p(q)/f(q) \to 0$. For a subshift with word complexity $p$ which is non-superlinear, i.e. $\liminf p(q)/q < \infty$, every ergodic probability measure is partially rigid.


title: FedER: Federated Learning through Experience Replay and Privacy-Preserving Data Synthesis

abstract: In the medical field, multi-center collaborations are often sought to yield more generalizable findings by leveraging the heterogeneity of patient and clinical data. However, recent privacy regulations hinder the possibility to share data, and consequently, to come up with machine learning-based solutions that support diagnosis and prognosis. Federated learning (FL) aims at sidestepping this limitation by bringing AI-based solutions to data owners and only sharing local AI models, or parts thereof, that need then to be aggregated. However, most of the existing federated learning solutions are still at their infancy and show several shortcomings, from the lack of a reliable and effective aggregation scheme able to retain the knowledge learned locally to weak privacy preservation as real data may be reconstructed from model updates. Furthermore, the majority of these approaches, especially those dealing with medical data, relies on a centralized distributed learning strategy that poses robustness, scalability and trust issues. In this paper we present a federated and decentralized learning strategy, FedER, that, exploiting experience replay and generative adversarial concepts, effectively integrates features from local nodes, providing models able to generalize across multiple datasets while maintaining privacy. FedER is tested on two tasks -- tuberculosis and melanoma classification -- using multiple datasets in order to simulate realistic non-i.i.d. medical data scenarios. Results show that our approach achieves performance comparable to standard (non-federated) learning and significantly outperforms state-of-the-art federated methods in their centralized (thus, more favourable) formulation. Code is available at https://github.com/perceivelab/FedER


title: The Capacity of 3 User Linear Computation Broadcast

abstract: The $K$ User Linear Computation Broadcast (LCBC) problem is comprised of $d$ dimensional data (from $\mathbb{F}_q$), that is fully available to a central server, and $K$ users, who require various linear computations of the data, and have prior knowledge of various linear functions of the data as side-information. The optimal broadcast cost is the minimum number of $q$-ary symbols to be broadcast by the server per computation instance, for every user to retrieve its desired computation. The reciprocal of the optimal broadcast cost is called the capacity. The main contribution of this paper is the exact capacity characterization for the $K=3$ user LCBC for all cases, i.e., for arbitrary finite fields $\mathbb{F}_q$, arbitrary data dimension $d$, and arbitrary linear side-informations and demands at each user. A remarkable aspect of the converse is that unlike the $2$ user LCBC whose capacity was determined previously, the entropic formulation (where the entropies of demands and side-informations are specified, but not their functional forms) is insufficient to obtain a tight converse for the $3$ user LCBC. Instead, the converse exploits functional submodularity. Notable aspects of achievability include a decomposition of the users' collective signal space into subspaces that allow different degrees of efficiency in broadcast cost, revealing a tradeoff that leads to a constrained water-filling solution. Random coding arguments are invoked to resolve compatibility issues that arise as each user has a different view of these subspaces, conditioned on its own side-information.


title: Better Incentives for Proof-of-Work

abstract: This work proposes a novel proof-of-work blockchain incentive scheme such that, barring exogenous motivations, following the protocol is guaranteed to be the optimal strategy for miners. Our blockchain takes the form of a directed acyclic graph, resulting in improvements with respect to throughput and speed. More importantly, for our blockchain to function, it is not expected that the miners conform to some presupposed protocol in the interest of the system's operability. Instead, our system works if miners act selfishly, trying to get the maximum possible rewards, with no consideration for the overall health of the blockchain.


title: WiNDS: An H$\alpha$ kinematics survey of nearby spiral galaxies -- Vertical perturbations in nearby disk-type galaxies

abstract: We present the Waves in Nearby Disk galaxies Survey (WiNDS) consisting of 40 nearby low inclination disk galaxies observed through H$\alpha$ high-resolution Fabry Perot interferometry. WiNDS consists of 12 new galaxy observations and 28 data archived observations obtained from different galaxy surveys. We derive two-dimensional line-of-sight velocity fields that are analyzed to identify the possible presence of vertical velocity flows in the galactic disks of these low-inclination late-type galaxies using velocity residual maps, derived from the subtraction of an axisymmetric rotation model to rotational velocity map. Large and globally coherent flows in the line-of-sight velocity of nearly face-on galaxies can be associated with large vertical displacement of the disk with respect to its mid-plane. Our goal is to characterize how frequent vertical perturbations, such as those observed in the Milky Way, arise in the Local Universe. Our currently available data have allowed us to identify 20$\%$ of WiNDS galaxies with strong velocity perturbations that are consistent with vertically perturbed galactic disks.


title: Theory of Radiation Pressure on a Diffractive Solar Sail

abstract: Solar sails propelled by radiation pressure enable space missions that cannot be achieved using chemical rockets alone. Significant in-space propulsion for missions such as a solar polar orbiter may be achieved with a sail that deviates sunlight at a large average angular direction. The momentum transfer efficiency of sunlight diffracted from a sun-facing diffractive sail comprised of periodic sawtooth prism elements is examined here. The spectrally averaged efficiency is found to approach that of a monolithic prism when the grating period is much longer than the peak of the solar spectrum. This idealized diffraction analysis predicts a greater transverse radiation pressure force compared to an idealized reflective sail. With modern optical design and fabrication techniques, diffractive solar sails may one day replace reflective sails.