title: Galileo, Ignoramus: Mathematics versus Philosophy in the Scientific Revolution

abstract: I offer a revisionist interpretation of Galileo's role in the history of science. My overarching thesis is that Galileo lacked technical ability in mathematics, and that this can be seen as directly explaining numerous aspects of his life's work. I suggest that it is precisely because he was bad at mathematics that Galileo was keen on experiment and empiricism, and eagerly adopted the telescope. His reliance on these hands-on modes of research was not a pioneering contribution to scientific method, but a last resort of a mind ill equipped to make a contribution on mathematical grounds. Likewise, it is precisely because he was bad at mathematics that Galileo expounded at length about basic principles of scientific method. "Those who can't do, teach." The vision of science articulated by Galileo was less original than is commonly assumed. It had long been taken for granted by mathematicians, who, however, did not stop to pontificate about such things in philosophical prose because they were too busy doing advanced scientific work.


title: An Antisymmetric Berry Frictional Force At Equilibrium in the Presence of Spin-Orbit Coupling

abstract: We analytically calculate the electronic friction tensor for a molecule near a metal surface in the case that the electronic Hamiltonian is complex-valued, e.g. the case that there is spin-orbit coupling and/or an external magnetic field. In such a case, even at equilibrium, we show that the friction tensor is not symmetric. Instead, the tensor is the real-valued sum of one positive definite tensor (corresponding to dissipation) plus one antisymmetric tensor (corresponding to a Berry pseudomagnetic force). Moreover, we find that this Berry force can be much larger than the dissipational force, suggesting the possibility of strongly spin-polarized chemicurrents or strongly spin-dependent rate constants for systems with spin-orbit coupling.


title: Li-Yau type inequality for curves in any codimension

abstract: For immersed curves in Euclidean space of any codimension we establish a Li--Yau type inequality that gives a lower bound of the (normalized) bending energy in terms of multiplicity. The obtained inequality is optimal for any codimension and any multiplicity except for the case of planar closed curves with odd multiplicity; in this remaining case we discover a hidden algebraic obstruction and indeed prove an exhaustive non-optimality result. The proof is mainly variational and involves Langer--Singer's classification of elasticae and Andr\'{e}'s algebraic-independence theorem for certain hypergeometric functions. We also discuss applications to elastic flows, networks, and knots.


title: Eight fold quantum Hall phases in a time reversal symmetry broken tight binding model

abstract: We consider a time reversal symmetry (TRS) broken Kane-Mele model superimposed with Haldane model and chart out the phase diagram using spin Chern number to investigate the fate of quantum anomalous Hall insulator (QAHI) and quantum spin Hall insulator (QSHI) phases. Interestingly, in addition to QSHI and QAHI phase, the phase diagram unveils quantum anomalous spin Hall insulator (QASHI) phase where only one spin sector is topological. We also find multicritical points where three / four topological phase boundaries coalesce. These topological phases are protected by an effective TRS and a composite anti-unitary particle-hole symmetry leading to remarkable properties of edge modes. We find spin-selective, spin-polarized and spin-neutral edge transport in QASHI, QSHI and QAHI phases respectively. Our study indicates that the robustness of the topological phase mainly depends on the spin gap which does not necessarily vanish at the Dirac points across a topological phase transition. We believe that our proposals can be tested in near future using recent experimental advancements in solid state and cold atomic systems.


title: Neural Architecture Search as Program Transformation Exploration

abstract: Improving the performance of deep neural networks (DNNs) is important to both the compiler and neural architecture search (NAS) communities. Compilers apply program transformations in order to exploit hardware parallelism and memory hierarchy. However, legality concerns mean they fail to exploit the natural robustness of neural networks. In contrast, NAS techniques mutate networks by operations such as the grouping or bottlenecking of convolutions, exploiting the resilience of DNNs. In this work, we express such neural architecture operations as program transformations whose legality depends on a notion of representational capacity. This allows them to be combined with existing transformations into a unified optimization framework. This unification allows us to express existing NAS operations as combinations of simpler transformations. Crucially, it allows us to generate and explore new tensor convolutions. We prototyped the combined framework in TVM and were able to find optimizations across different DNNs, that significantly reduce inference time - over 3$\times$ in the majority of cases. Furthermore, our scheme dramatically reduces NAS search time. Code is available at~\href{https://github.com/jack-willturner/nas-as-program-transformation-exploration}{this https url}.


title: Shape- and element-sensitive reconstruction of periodic nanostructures with grazing incidence X-ray fluorescence analysis and machine learning

abstract: The characterization of nanostructured surfaces with sensitivity in the sub-nm range is of high importance for the development of current and next generation integrated electronic circuits. Modern transistor architectures for e.g. FinFETs are realized by lithographic fabrication of complex, well ordered nanostructures. Recently, a novel characterization technique based on X-ray fluorescence measurements in grazing incidence geometry has been proposed for such applications. This technique uses the X-ray standing wave field, arising from an interference between incident and the reflected radiation, as a nanoscale sensor for the dimensional and compositional parameters of the nanostructure. The element sensitivity of the X-ray fluorescence technique allows for a reconstruction of the spatial element distribution using a finite-element method. Due to a high computational time, intelligent optimization methods employing machine learning algorithms are essential for a timely provision of results. Here, a sampling of the probability distributions by Bayesian optimization is not only fast, it also provides an initial estimate of the parameter uncertainties and sensitivities. The high sensitivity of the method requires a precise knowledge of the material parameters in the modeling of the dimensional shape provided that some physical properties of the material are known or determined beforehand. The unknown optical constants were extracted from an unstructured but otherwise identical layer system by means of soft X-ray reflectometry. The spatial distribution profiles of the different elements contained in the grating structure were compared to scanning electron and atomic force microscopy and the influence of carbon surface contamination on the modeling results were discussed.


title: 3D-1D coupling on non conforming meshes via three-field optimization based domain decomposition

abstract: A new numerical approach is proposed for the simulation of coupled three-dimensional and one-dimensional elliptic equations (3D-1D coupling) arising from dimensionality reduction of 3D-3D problems with thin inclusions. The method is based on a well posed mathematical formulation and results in a numerical scheme with high robustness and flexibility in handling geometrical complexities. This is achieved by means of a three-field approach to split the 1D problems from the bulk 3D problem, and then resorting to the minimization of a properly designed functional to impose matching conditions at the interfaces. Thanks to the structure of the functional, the method allows the use of independent meshes for the various subdomains.


title: Modeling Dynamic User Interests: A Neural Matrix Factorization Approach

abstract: In recent years, there has been significant interest in understanding users' online content consumption patterns. But, the unstructured, high-dimensional, and dynamic nature of such data makes extracting valuable insights challenging. Here we propose a model that combines the simplicity of matrix factorization with the flexibility of neural networks to efficiently extract nonlinear patterns from massive text data collections relevant to consumers' online consumption patterns. Our model decomposes a user's content consumption journey into nonlinear user and content factors that are used to model their dynamic interests. This natural decomposition allows us to summarize each user's content consumption journey with a dynamic probabilistic weighting over a set of underlying content attributes. The model is fast to estimate, easy to interpret and can harness external data sources as an empirical prior. These advantages make our method well suited to the challenges posed by modern datasets. We use our model to understand the dynamic news consumption interests of Boston Globe readers over five years. Thorough qualitative studies, including a crowdsourced evaluation, highlight our model's ability to accurately identify nuanced and coherent consumption patterns. These results are supported by our model's superior and robust predictive performance over several competitive baseline methods.


title: Semantically-Conditioned Negative Samples for Efficient Contrastive Learning

abstract: Negative sampling is a limiting factor w.r.t. the generalization of metric-learned neural networks. We show that uniform negative sampling provides little information about the class boundaries and thus propose three novel techniques for efficient negative sampling: drawing negative samples from (1) the top-$k$ most semantically similar classes, (2) the top-$k$ most semantically similar samples and (3) interpolating between contrastive latent representations to create pseudo negatives. Our experiments on CIFAR-10, CIFAR-100 and Tiny-ImageNet-200 show that our proposed \textit{Semantically Conditioned Negative Sampling} and Latent Mixup lead to consistent performance improvements. In the standard supervised learning setting, on average we increase test accuracy by 1.52\% percentage points on CIFAR-10 across various network architectures. In the knowledge distillation setting, (1) the performance of student networks increase by 4.56\% percentage points on Tiny-ImageNet-200 and 3.29\% on CIFAR-100 over student networks trained with no teacher and (2) 1.23\% and 1.72\% respectively over a \textit{hard-to-beat} baseline (Hinton et al., 2015).


title: Cockpit: A Practical Debugging Tool for the Training of Deep Neural Networks

abstract: When engineers train deep learning models, they are very much 'flying blind'. Commonly used methods for real-time training diagnostics, such as monitoring the train/test loss, are limited. Assessing a network's training process solely through these performance indicators is akin to debugging software without access to internal states through a debugger. To address this, we present Cockpit, a collection of instruments that enable a closer look into the inner workings of a learning machine, and a more informative and meaningful status report for practitioners. It facilitates the identification of learning phases and failure modes, like ill-chosen hyperparameters. These instruments leverage novel higher-order information about the gradient distribution and curvature, which has only recently become efficiently accessible. We believe that such a debugging tool, which we open-source for PyTorch, is a valuable help in troubleshooting the training process. By revealing new insights, it also more generally contributes to explainability and interpretability of deep nets.


title: Unleashing the Power of Contrastive Self-Supervised Visual Models via Contrast-Regularized Fine-Tuning

abstract: Contrastive self-supervised learning (CSL) has attracted increasing attention for model pre-training via unlabeled data. The resulted CSL models provide instance-discriminative visual features that are uniformly scattered in the feature space. During deployment, the common practice is to directly fine-tune CSL models with cross-entropy, which however may not be the best strategy in practice. Although cross-entropy tends to separate inter-class features, the resulting models still have limited capability for reducing intra-class feature scattering that exists in CSL models. In this paper, we investigate whether applying contrastive learning to fine-tuning would bring further benefits, and analytically find that optimizing the contrastive loss benefits both discriminative representation learning and model optimization during fine-tuning. Inspired by these findings, we propose Contrast-regularized tuning (Core-tuning), a new approach for fine-tuning CSL models. Instead of simply adding the contrastive loss to the objective of fine-tuning, Core-tuning further applies a novel hard pair mining strategy for more effective contrastive fine-tuning, as well as smoothing the decision boundary to better exploit the learned discriminative feature space. Extensive experiments on image classification and semantic segmentation verify the effectiveness of Core-tuning.


title: Optimisation of an active heat engine

abstract: Optimisation of heat engines at the micro-scale has applications in biological and artificial nano-technology, and stimulates theoretical research in non-equilibrium statistical physics. Here we consider non-interacting overdamped particles confined by an external harmonic potential, in contact either with a thermal reservoir or with a stochastic self-propulsion force (active Ornstein-Uhlenbeck model). A cyclical machine is produced by periodic variation of the parameters of the potential and of the noise. An exact mapping between the passive and the active model allows us to define the effective temperature $T_{eff}(t)$ which is meaningful for the thermodynamic performance of the engine. We show that $T_{eff}(t)$ is different from all other known active temperatures, typically used in static situations. The mapping allows us to optimise the active engine, whatever are the values of the persistence time or self-propulsion velocity. In particular - through linear irreversible thermodynamics (small amplitude of the cycle) - we give explicit formula for the optimal cycle period and phase delay (between the two modulated parameters, stiffness and temperature) achieving maximum power with Curzon-Ahlborn efficiency. In the quasi-static limit, the formula for $T_{eff}(t)$ simplifies and coincides with a recently proposed temperature for stochastic thermodynamics, bearing a compact expression for the maximum efficiency. A point - overlooked in recent literature - is made about the difficulty in defining efficiency without a consistent definition of effective temperature.


title: Intelligent Software Web Agents: A Gap Analysis

abstract: Semantic web technologies have shown their effectiveness, especially when it comes to knowledge representation, reasoning, and data integration. However, the original semantic web vision, whereby machine readable web data could be automatically actioned upon by intelligent software web agents, has yet to be realised. In order to better understand the existing technological opportunities and challenges, in this paper we examine the status quo in terms of intelligent software web agents, guided by research with respect to requirements and architectural components, coming from the agents community. We use the identified requirements to both further elaborate on the semantic web agent motivating use case scenario, and to summarise different perspectives on the requirements from the semantic web agent literature. We subsequently propose a hybrid semantic web agent architecture, and use the various components and subcomponents in order to provide a focused discussion in relation to existing semantic web standards and community activities. Finally, we highlight open research opportunities and challenges and take a broader perspective of the research by discussing the potential for intelligent software web agents as an enabling technology for emerging domains, such as digital assistants, cloud computing, and the internet of things.


title: Stable Bloch oscillations and Landau-Zener tunneling in a non-Hermitian $\cal{PT}$-symmetric flat band lattice

abstract: This article aims to study the existence of stable Bloch oscillations and Landau-Zener tunneling in a non-Hermitian system when exposed to external fields. We investigate a non-Hermitian $\cal{PT}$-symmetric diamond chain network and its transport dynamics in two different situations, namely in a flat band case and a non-flat band case. The considered system does not support unbroken-$\cal{PT}$ phase or completely real eigenspectra in any of the parametric regions in both the flat and non-flat band cases. In the flat band case, up to a critical value of the gain-loss parameter, the bands are found to be gapless or inseparable, and for other values the bands are isolated. Considering the non-flat band case, all the bands are found to be complex dispersive and are also isolated. In the case of completely broken $\cal{PT}$ phase, we look upon the possibility to have stable dynamics or Bloch oscillations upon the application of external fields like synthetic electric field. In particular, when the complex bands are isolated, we point out that the Landau-Zener tunneling induced by the synthetic electric field can enable Bloch oscillations. The amplitude of these Bloch oscillations is large and persists for a long propagation distance which reveals that super Bloch oscillations can be observed in the broken $\cal{PT}$ phase of the system. We also report the amplified Bloch oscillations which pave the way towards controlling transport phenomena in non-Hermitian systems.


title: Model-based Prediction and Optimal Control of Pandemics by Non-pharmaceutical Interventions

abstract: A model-based signal processing framework is proposed for pandemic trend forecasting and control, by using non-pharmaceutical interventions (NPI) at regional and country levels worldwide. The control objective is to prescribe quantifiable NPI strategies at different levels of stringency, which balance between human factors (such as new cases and death rates) and cost of intervention per region/country. Due to infrastructural disparities and differences in priorities of regions and countries, strategists are given the flexibility to weight between different NPIs and to select the desired balance between the human factor and overall NPI cost. The proposed framework is based on a \textit{finite-horizon optimal control} (FHOC) formulation of the bi-objective problem and the FHOC is numerically solved by using an ad hoc \textit{extended Kalman filtering/smoothing} framework for optimal NPI estimation and pandemic trend forecasting. The algorithm enables strategists to select the desired balance between the human factor and NPI cost with a set of weights and parameters. The parameters of the model are partially selected by epidemiological facts from COVID-19 studies, and partially trained by using machine learning techniques. The developed algorithm is applied on ground truth data from the Oxford COVID-19 Government Response Tracker project, which has categorized and quantified the regional responses to the pandemic for more than 300 countries and regions worldwide, since January 2020. The dataset was used for NPI-based prediction and prescription during the XPRIZE Pandemic Response Challenge.


title: Enhancing into the codec: Noise Robust Speech Coding with Vector-Quantized Autoencoders

abstract: Audio codecs based on discretized neural autoencoders have recently been developed and shown to provide significantly higher compression levels for comparable quality speech output. However, these models are tightly coupled with speech content, and produce unintended outputs in noisy conditions. Based on VQ-VAE autoencoders with WaveRNN decoders, we develop compressor-enhancer encoders and accompanying decoders, and show that they operate well in noisy conditions. We also observe that a compressor-enhancer model performs better on clean speech inputs than a compressor model trained only on clean speech.


title: Understanding population annealing Monte Carlo simulations

abstract: Population annealing is a recent addition to the arsenal of the practitioner in computer simulations in statistical physics and beyond that is found to deal well with systems with complex free-energy landscapes. Above all else, it promises to deliver unrivaled parallel scaling qualities, being suitable for parallel machines of the biggest calibre. Here we study population annealing using as the main example the two-dimensional Ising model which allows for particularly clean comparisons due to the available exact results and the wealth of published simulational studies employing other approaches. We analyze in depth the accuracy and precision of the method, highlighting its relation to older techniques such as simulated annealing and thermodynamic integration. We introduce intrinsic approaches for the analysis of statistical and systematic errors, and provide a detailed picture of the dependence of such errors on the simulation parameters. The results are benchmarked against canonical and parallel tempering simulations.


title: Pre-processing of Galaxies in the Early Stages of Cluster Formation in Abell 1882 at $z$=0.139

abstract: A rare opportunity to distinguish internal and environmental effects on galaxy evolution is afforded by "Supergroups", systems which are rich and massive but include several comparably rich substructures, surrounded by filaments. We present here a multiwavelength photometric and spectroscopic study of the galaxy population in the Supergroup Abell 1882 at $z$=0.139, combining new data from the MMT and Hectospec with archival results from the Galaxy And Mass Assembly (GAMA) survey, the Sloan Digital Sky Survey (SDSS), NED, the Gemini Multi-Object Spectrograph (GMOS) and Galaxy Evolution Explorer (GALEX). These provide spectroscopic classifications for 526 member galaxies, across wide ranges of local density and velocity dispersion. We identify three prominent filaments along which galaxies seem to be entering the Supergroup (mostly in E-W directions). Abell 1882 has a well-populated red sequence, containing most galaxies with stellar mass $>$ $10^{10.5}$ $M_{sun}$, and a pronounced color-density relation even within its substructures. Thus, galaxy evolution responds to the external environment as strongly in these unrelaxed systems as we find in rich and relaxed clusters. From these data, local density remains the primary factor with a secondary role for distance from the inferred center of the entire structure's potential well. Effects on star formation, as traced by optical and near-UV colors, depend on galaxy mass. We see changes in lower-mass galaxies (M $<$ $10^{10.5}$ $M_{sun}$) at four times the virial radius of major substructures, while the more massive $NUV$ Green Valley galaxies show low levels of star formation within two virial radii. Suppression of star formation ("quenching") occurs in the infall regions of these structures even before galaxies enter the denser group environment.


title: Improved LP-based Approximation Algorithms for Facility Location with Hard Capacities

abstract: The Capacitated Facility Location (CFL), a long-standing classic problem with intriguing approximability and literature dated back to the 90s, is considered. Following the open question posted in [Williamson and Shmoys, 2011] and the notable work due to [An et al., FOCS~2014], we present an LP-based approximation algorithm with a guarantee of $(10+\sqrt{67})/2 \approx 9.0927$, a significant improvement upon the previous LP-based ratio of $288$ due to An et al. in 2014. Our contribution for this part is a simple and elegant rounding algorithm that brings clear insights for the MFN relaxation and the CFL problem. For CFL with cardinality facility cost (CFL-CFC), we present an LP-based $4$-approximation algorithm, which improves upon the decades-old ratio of 5 due to Levi et al. that ages up since 2004. Prior to our work, it was not clear whether or not LP-based methods can be used to provide a guarantee better than 5 for the CFL problem, even for restricted versions of this problem, for which natural LPs are already known to have small integrality gaps. Our rounding algorithm provides the first affirmative answer on the case with cadinality facility cost.


title: TerraWatt: Sustaining Sustainable Computing of Containers in Containers

abstract: Each day the world inches closer to a climate catastrophe and a sustainability revolution. To avoid the former and achieve the latter we must transform our use of energy. Surprisingly, today's growing problem is that there is too much wind and solar power generation at the wrong times and in the wrong places. We argue for the construction of TerraWatt: a geographically-distributed, large-scale, zero-carbon compute infrastructure using renewable energy and older hardware. Delivering zero-carbon compute for general cloud workloads is challenging due to spatiotemporal power variability. We describe the systems challenges in using intermittent renewable power at scale to fuel such an older, decentralized compute infrastructure.


title: Exciton Dynamics in Conjugated Polymers

abstract: This is a review of exciton dynamics in conjugated polymers. Exciton dynamics encompass multiple time and length scales. Ultrafast femtosecond processes are intrachain and involve a quantum mechanical correlation of the exciton and nuclear degrees of freedom. In contrast, post-picosecond processes involve the incoherent Forster transfer of excitons between polymer chains. Exciton dynamics is also strongly determined by the spatial and temporal disorder that is ubiquitous in conjugated polymers. Since excitons are delocalized over hundreds of atoms, a theoretical understanding of these processes is only realistically possible by employing suitably parametrized coarse-grained exciton-phonon models. Moreover, to correctly account for ultrafast processes, the exciton and phonon modes must be treated on the same quantum mechanical basis and the Ehrenfest approximation must be abandoned. This further implies that sophisticated numerical techniques must be employed to solve these models. We begin by describing the energetic and spatial distribution of excitons in disordered polymer systems. Next, we discuss ultrafast intrachain exciton decoherence caused by exciton-phonon entanglement, which leads to fluorescence depolarization on the timescale of 10-fs. Interactions of the polymer with its environment causes the stochastic relaxation and localization of high-energy delocalized excitons onto chromophores. The coupling of excitons with torsional modes also leads to various dynamical processes. On sub-ps timescales it causes exciton density localization and local polymer planarization, while on post-ps timescales stochastic torsional fluctuations cause exciton diffusion along the polymer chain. Finally, we describe a first-principles, Forster-type model of intrachain exciton transfer and diffusion, whose starting point is a realistic description of the donor and acceptor chromophores.


title: Multi-access Coded Caching Scheme with Linear Sub-packetization using PDAs

abstract: We consider multi-access coded caching problem introduced by Hachem et.al., where each user has access to $L$ neighboring caches in a cyclic wrap-around fashion. We focus on the deterministic schemes for a specific class of multi-access coded caching problem based on the concept of PDA. We construct new PDAs which specify the delivery scheme for the specific class of multi-access coded caching problem discussed in this paper. For the proposed scheme, the coding gain is larger than that of the state-of-the-art while the sub-packetization level varies only linearly with the number of users. Hence, we achieve a lower transmission rate with the least sub-packetization level compared to the existing schemes.


title: Implementing the MULTI-VP coronal model in EUHFORIA: test case results and comparisons with the WSA coronal model

abstract: In this study, we focus on improving EUHFORIA (European Heliospheric Forecasting Information Asset), a recently developed 3D MHD space weather prediction tool. EUHFORIA consists of two parts, covering two spatial domains; the solar corona and the inner heliosphere. For the first part, the semi-empirical Wang-Sheeley-Arge (WSA) model is used by default, which employs the Potential Field Source Surface (PFSS) and Schatten Current Sheet (SCS) models to provide the necessary solar wind plasma and magnetic conditions above the solar surface, at 0.1 AU, that serve as boundary conditions for the inner heliospheric part. Herein, we present the first results of the implementation of an alternative coronal model in EUHFORIA, the so-called MULTI-VP model. We compared the output of the default coronal model with the output from MULTI-VP at the inner boundary of the heliospheric domain of EUHFORIA in order to understand differences between the two models, before they propagate to Earth. We also compared the performance of WSA+EUHFORIA-heliosphere and MULTI-VP+EUHFORIA-heliosphere against in situ observations at Earth. In the frame of this study, we considered two different high-speed stream cases, one during a period of low solar activity and one during a period of high solar activity. We also employed two different magnetograms, i.e., GONG and WSO. Our results show that the choice of both the coronal model and the magnetogram play an important role on the accuracy of the solar wind prediction. However, it is not clear which component plays the most important role for the modeled results obtained at Earth. A statistical analysis with an appropriate number of simulations is needed to confirm our findings.


title: Extending superposed harmonic initial data to higher spin

abstract: Numerical simulations of binary black holes are accompanied by an initial spurious burst of gravitational radiation (called `junk radiation') caused by a failure of the initial data to describe a snapshot of an inspiral that started at an infinite time in the past. A previous study showed that the superposed harmonic (SH) initial data gives rise to significantly smaller junk radiation. However, it is difficult to construct SH initial data for black holes with dimensionless spin $\chi\gtrsim0.7$. We here provide a class of spatial coordinate transformations that extend SH to higher spin. The new spatial coordinate system, which we refer to as superposed modified harmonic (SMH), is characterized by a continuous parameter -- Kerr-Schild and harmonic spatial coordinates are only two special cases of this new gauge. We compare SMH with the superposed Kerr-Schild (SKS) initial data by evolving several binary black hole systems with $\chi=0.8$ and $0.9$. We find that the new initial data still leads to less junk radiation and only small changes of black hole parameters (e.g. mass and spin). We also find that the volume-weighted constraint violations for the new initial data converge with resolution during the junk stage $(t\lesssim700M)$, which means there are fewer high-frequency components in waveforms at outer regions.


title: Differences in the spatial landscape of urban mobility: gender and socioeconomic perspectives

abstract: Many of our routines and activities are linked to our ability to move; be it commuting to work, shopping for groceries, or meeting friends. Yet, factors that limit the individuals' ability to fully realise their mobility needs will ultimately affect the opportunities they can have access to (e.g., cultural activities, professional interactions). One important aspect frequently overlooked in human mobility studies is how gender-centred issues can amplify other sources of mobility disadvantages (e.g., socioeconomic inequalities), unevenly affecting the pool of opportunities men and women have access to. In this work, we leverage on a combination of computational, statistical, and information-theoretical approaches to investigate the existence of systematic discrepancies in the mobility diversity (i.e., the diversity of travel destinations) of (1) men and women from different socioeconomic backgrounds, and (2) work and non-work travels. Our analysis is based on datasets containing multiple instances of large-scale, official, travel surveys carried out in three major metropolitan areas in South America: Medell\'in and Bogot\'a in Colombia, and S\~ao Paulo in Brazil. Our results indicate the presence of general discrepancies in the urban mobility diversities related to the gender and socioeconomic characteristics of the individuals. Lastly, this paper sheds new light on the possible origins of gender-level human mobility inequalities, contributing to the general understanding of disaggregated patterns in human mobility.


title: Hidden regular variation for point processes and the single/multiple large point heuristic

abstract: We consider regular variation for marked point processes with independent heavy-tailed marks and prove a single large point heuristic: the limit measure is concentrated on the cone of point measures with one single point. We then investigate successive hidden regular variation removing the cone of point measures with at most $k$ points, $k\geq 1$, and prove a multiple large point phenomenon: the limit measure is concentrated on the cone of point measures with $k+1$ points. We show how these results imply hidden regular variation in Skorokhod space of the associated risk process. Finally, we provide an application to risk theory in a reinsurance model where the $k$ largest claims are covered and we study the asymptotic behavior of the residual risk.


title: Optimizing Inference Performance of Transformers on CPUs

abstract: The Transformer architecture revolutionized the field of natural language processing (NLP). Transformers-based models (e.g., BERT) power many important Web services, such as search, translation, question-answering, etc. While enormous research attention is paid to the training of those models, relatively little efforts are made to improve their inference performance. This paper comes to address this gap by presenting an empirical analysis of scalability and performance of inferencing a Transformer-based model on CPUs. Focusing on the highly popular BERT model, we identify key components of the Transformer architecture where the bulk of the computation happens, and propose three optimizations to speed them up. The optimizations are evaluated using the inference benchmark from HuggingFace, and are shown to achieve the speedup of up to x2.37. The considered optimizations do not require any changes to the implementation of the models nor affect their accuracy.


title: MetaGrad: Adaptation using Multiple Learning Rates in Online Learning

abstract: We provide a new adaptive method for online convex optimization, MetaGrad, that is robust to general convex losses but achieves faster rates for a broad class of special functions, including exp-concave and strongly convex functions, but also various types of stochastic and non-stochastic functions without any curvature. We prove this by drawing a connection to the Bernstein condition, which is known to imply fast rates in offline statistical learning. MetaGrad further adapts automatically to the size of the gradients. Its main feature is that it simultaneously considers multiple learning rates, which are weighted directly proportional to their empirical performance on the data using a new meta-algorithm. We provide three versions of MetaGrad. The full matrix version maintains a full covariance matrix and is applicable to learning tasks for which we can afford update time quadratic in the dimension. The other two versions provide speed-ups for high-dimensional learning tasks with an update time that is linear in the dimension: one is based on sketching, the other on running a separate copy of the basic algorithm per coordinate. We evaluate all versions of MetaGrad on benchmark online classification and regression tasks, on which they consistently outperform both online gradient descent and AdaGrad.


title: Chiral magnetic effect reveals the topology of gauge fields in heavy-ion collisions

abstract: The topological structure of vacuum is the cornerstone of non-Abelian gauge theories describing strong and electroweak interactions within the standard model of particle physics. However, transitions between different topological sectors of the vacuum (believed to be at the origin of the baryon asymmetry of the Universe) have never been observed directly. An experimental observation of such transitions in Quantum Chromodynamics (QCD) has become possible in heavy-ion collisions, where the chiral magnetic effect converts the chiral asymmetry (generated by topological transitions in hot QCD matter) into an electric current, under the presence of the magnetic field produced by the colliding ions. The Relativistic Heavy Ion Collider program on heavy-ion collisions such as the Zr-Zr and Ru-Ru isobars, thus has the potential to uncover the topological structure of vacuum in a laboratory experiment. This discovery would have far-reaching implications for the understanding of QCD, the origin of the baryon asymmetry in the present-day Universe, and for other areas, including condensed matter physics.

categories: hep-ph cond-mat.str-el hep-th nucl-ex nucl-th


title: A Generative Model for Hallucinating Diverse Versions of Super Resolution Images

abstract: Traditionally, the main focus of image super-resolution techniques is on recovering the most likely high-quality images from low-quality images, using a one-to-one low- to high-resolution mapping. Proceeding that way, we ignore the fact that there are generally many valid versions of high-resolution images that map to a given low-resolution image. We are tackling in this work the problem of obtaining different high-resolution versions from the same low-resolution image using Generative Adversarial Models. Our learning approach makes use of high frequencies available in the training high-resolution images for preserving and exploring in an unsupervised manner the structural information available within these images. Experimental results on the CelebA dataset confirm the effectiveness of the proposed method, which allows the generation of both realistic and diverse high-resolution images from low-resolution images.


title: Corrections to the hadron resonance gas from lattice QCD and their effect on fluctuation-ratios at finite density

abstract: The hadron resonance gas (HRG) model is often believed to correctly describe the confined phase of QCD. This assumption is the basis of many phenomenological works on QCD thermodynamics and of the analysis of hadron yields in relativistic heavy ion collisions. We use first-principle lattice simulations to calculate corrections to the ideal HRG. Namely, we determine the sub-leading fugacity expansion coefficients of the grand canonical free energy, receiving contributions from processes like kaon-kaon or baryon-baryon scattering. We achieve this goal by performing a two dimensional scan on the imaginary baryon number chemical potential ($\mu_B$) - strangeness chemical potential ($\mu_S$) plane, where the fugacity expansion coefficients become Fourier coefficients. We carry out a continuum limit estimation of these coefficients by performing lattice simulations with temporal extents of $N_\tau=8,10,12$ using the 4stout-improved staggered action. We then use the truncated fugacity expansion to extrapolate ratios of baryon number and strangeness fluctuations and correlations to finite chemical potentials. Evaluating the fugacity expansion along the crossover line, we reproduce the trend seen in the experimental data on net-proton fluctuations by the STAR collaboration.


title: Do-calculus enables estimation of causal effects in partially observed biomolecular pathways

abstract: Estimating causal queries, such as changes in protein abundance in response to a perturbation, is a fundamental task in the analysis of biomolecular pathways. The estimation requires experimental measurements on the pathway components. However, in practice many pathway components are left unobserved (latent) because they are either unknown, or difficult to measure. Latent variable models (LVMs) are well-suited for such estimation. Unfortunately, LVM-based estimation of causal queries can be inaccurate when parameters of the latent variables are not uniquely identified, or when the number of latent variables is misspecified. This has limited the use of LVMs for causal inference in biomolecular pathways. In this manuscript, we propose a general and practical approach for LVM-based estimation of causal queries. We prove that, despite the challenges above, LVM-based estimators of causal queries are accurate if the queries are identifiable according to Pearl's do-calculus, and describe an algorithm for its estimation. We illustrate the breadth and the practical utility of this approach for estimating causal queries in four synthetic and two experimental case studies, where structures of biomolecular pathways challenge the existing methods for causal query estimation. The code and the data documenting all the case studies are available at \url{https://github.com/srtaheri/LVMwithDoCalculus}


title: An efficient method for goal-oriented linear Bayesian optimal experimental design: Application to optimal sensor placemen

abstract: Optimal experimental design (OED) plays an important role in the problem of identifying uncertainty with limited experimental data. In many applications, we seek to minimize the uncertainty of a predicted quantity of interest (QoI) based on the solution of the inverse problem, rather than the inversion model parameter itself. In these scenarios, we develop an efficient method for goal-oriented optimal experimental design (GOOED) for large-scale Bayesian linear inverse problem that finds sensor locations to maximize the expected information gain (EIG) for a predicted QoI. By deriving a new formula to compute the EIG, exploiting low-rank structures of two appropriate operators, we are able to employ an online-offline decomposition scheme and a swapping greedy algorithm to maximize the EIG at a cost measured in model solutions that is independent of the problem dimensions. We provide detailed error analysis of the approximated EIG, and demonstrate the efficiency, accuracy, and both data- and parameter-dimension independence of the proposed algorithm for a contaminant transport inverse problem with infinite-dimensional parameter field.


title: How to Make Traversable Wormholes: Eternal AdS$_4$ Wormholes from Coupled CFT's

abstract: We construct an eternal traversable wormhole connecting two asymptotically $\text{AdS}_4$ regions. The wormhole is dual to the ground state of a system of two identical holographic CFT's coupled via a single low-dimension operator. The coupling between the two CFT's leads to negative null energy in the bulk, which supports a static traversable wormhole. As the ground state of a simple Hamiltonian, it may be possible to make these wormholes in the lab or on a quantum computer.


title: The art of molecular computing: whence and whither

abstract: An astonishingly diverse biomolecular circuitry orchestrates the functioning machinery underlying every living cell. These biomolecules and their circuits have been engineered not only for various industrial applications but also to perform other atypical functions that they were not evolved for - including computation. Various kinds of computational challenges, such as solving NP-complete problems with many variables, logical computation, neural network operations, and cryptography, have all been attempted through this unconventional computing paradigm. In this review, we highlight key experiments across three different eras of molecular computation, beginning with molecular solutions, transitioning to logic circuits and ultimately, more complex molecular networks. We also discuss a variety of applications of molecular computation, from solving NP-hard problems to self-assembled nanostructures for delivering molecules, and provide a glimpse into the exciting potential that molecular computing holds for the future.


title: Replica Symmetry Breaking and Phase Transitions in a PT Symmetric Sachdev-Ye-Kitaev Model

abstract: We show that the low temperature phase of a conjugate pair of uncoupled, quantum chaotic, nonhermitian systems such as the Sachdev-Ye-Kitaev (SYK) model or the Ginibre ensemble of random matrices are dominated by replica symmetry breaking configurations with a nearly flat free energy that terminates in a first order phase transition. In the case of the SYK model, we show explicitly that the spectrum of the effective replica theory has a gap. These features are strikingly similar to those induced by wormholes in the gravity path integral which suggests a close relation between both configurations. For a non-chaotic SYK, the results are qualitatively different: the spectrum is gapless in the low temperature phase and there is an infinite number of second order phase transitions unrelated to the restoration of replica symmetry.


title: Downstream depolarization in the Sausage relic: a 1-4 GHz Very Large Array study

abstract: Radio relics are elongated sources related to shocks driven by galaxy cluster merger events. Although these objects are highly polarized at GHz frequencies ($\gtrsim 20\%$), high-resolution studies of their polarization properties are still lacking. We present the first high-resolution and high-sensitivity polarimetry study of the merging galaxy cluster CIZA J2242.8+5301 in the 1-4 GHz frequency band. We use the $QU$-fitting approach to model the Stokes $I$, $Q$ and $U$ emission, obtaining best-fit intrinsic polarization fraction ($p_0$), intrinsic polarization angle ($\chi_0$), Rotation Measure (RM) and wavelength-dependent depolarization ($\sigma_{\rm RM}$) maps of the cluster. Our analysis focuses on the northern relic (RN). For the first time in a radio relic, we observe a decreasing polarization fraction in the downstream region. Our findings are possibly explained by geometrical projections and/or by decreasing of the magnetic field anisotropy towards the cluster center. From the amount of depolarization of the only detected background radio galaxy, we estimate a turbulent magnetic field strength of $B_{\rm turb}\sim5.6~\mu$Gauss in the relic. Finally, we observe Rotation Measure fluctuations of about 30 rad m$^{-2}$ around at the median value of 140.8 rad m$^{-2}$ at the relic position.


title: Deep Reinforcement Learning for Backup Strategies against Adversaries

abstract: Many defensive measures in cyber security are still dominated by heuristics, catalogs of standard procedures, and best practices. Considering the case of data backup strategies, we aim towards mathematically modeling the underlying threat models and decision problems. By formulating backup strategies in the language of stochastic processes, we can translate the challenge of finding optimal defenses into a reinforcement learning problem. This enables us to train autonomous agents that learn to optimally support planning of defense processes. In particular, we tackle the problem of finding an optimal backup scheme in the following adversarial setting: Given $k$ backup devices, the goal is to defend against an attacker who can infect data at one time but chooses to destroy or encrypt it at a later time, potentially also corrupting multiple backups made in between. In this setting, the usual round-robin scheme, which always replaces the oldest backup, is no longer optimal with respect to avoidable exposure. Thus, to find a defense strategy, we model the problem as a hybrid discrete-continuous action space Markov decision process and subsequently solve it using deep deterministic policy gradients. We show that the proposed algorithm can find storage device update schemes which match or exceed existing schemes with respect to various exposure metrics.


title: Discrete-Time Consensus Networks: Scalability, Grounding and Countermeasures

abstract: We investigate the disruption of discrete-time consensus problems via grounding. Loosely speaking, grounding a network occurs if the state of one agent no longer responds to inputs from other agents and/or changes its dynamics. Then, the agent becomes a leader or a so-called stubborn agent. The disruption of the agent can be caused by internal faults, safety protocols or due to an external malicious attack. In this paper we investigate how grounding affects expander graph families that usually exhibit good scaling properties with increasing network size. It is shown that the algebraic connectivity and eigenratio of the network decrease due to the grounding causing the performance and scalability of the network to deteriorate, even to the point of losing consensusability. We then present possible countermeasures to such disruptions and discuss their practicality and limitations. In particular, for a specific countermeasure of deliberately grounding additional nodes, we investigate extensively how to select additional nodes to ground and how many nodes we need to ground to recover the consensus performance. Our findings are supported by a wide range of numerical simulations.


title: An Overview of Recommender Systems and Machine Learning in Feature Modeling and Configuration

abstract: Recommender systems support decisions in various domains ranging from simple items such as books and movies to more complex items such as financial services, telecommunication equipment, and software systems. In this context, recommendations are determined, for example, on the basis of analyzing the preferences of similar users. In contrast to simple items which can be enumerated in an item catalog, complex items have to be represented on the basis of variability models (e.g., feature models) since a complete enumeration of all possible configurations is infeasible and would trigger significant performance issues. In this paper, we give an overview of a potential new line of research which is related to the application of recommender systems and machine learning techniques in feature modeling and configuration. In this context, we give examples of the application of recommender systems and machine learning and discuss future research issues.


title: ReLU Neural Networks of Polynomial Size for Exact Maximum Flow Computation

abstract: This paper studies the expressive power of artificial neural networks with rectified linear units. In order to study them as a model of real-valued computation, we introduce the concept of Max-Affine Arithmetic Programs and show equivalence between them and neural networks concerning natural complexity measures. We then use this result to show that two fundamental combinatorial optimization problems can be solved with polynomial-size neural networks. First, we show that for any undirected graph with $n$ nodes, there is a neural network (with fixed weights and biases) of size $\mathcal{O}(n^3)$ that takes the edge weights as input and computes the value of a minimum spanning tree of the graph. Second, we show that for any directed graph with $n$ nodes and $m$ arcs, there is a neural network of size $\mathcal{O}(m^2n^2)$ that takes the arc capacities as input and computes a maximum flow. Our results imply that these two problems can be solved with strongly polynomial time algorithms that solely uses affine transformations and maxima computations, but no comparison-based branchings.


title: Neutron scattering studies on spin fluctuations in Sr$_2$RuO$_4$

abstract: The magnetic excitations in Sr$_2$RuO$_4$ are studied by polarized and unpolarized neutron scattering experiments as a function of temperature. At the scattering vector of the Fermi-surface nesting with a half-integer out-of-plane component, there is no evidence for the appearance of a resonance excitation in the superconducting phase. The body of existing data indicates weakening of the scattered intensity in the nesting spectrum to occur at very low energies. The nesting signal persists up to 290 K but is strongly reduced. In contrast, a quasi-ferromagnetic contribution maintains its strength and still exhibits a finite width in momentum space.


title: Griffiths positivity for Bismut curvature and its behaviour along Hermitian Curvature Flows

abstract: In this note we study a positivity notion for the curvature of the Bismut connection; more precisely, we study the notion of \emph{Bismut-Griffiths-positivity} for complex Hermitian non-K\"ahler manifolds. Since the K\"ahler-Ricci flow preserves and regularizes the usual Griffiths positivity we investigate the behaviour of the Bismut-Griffiths-positivity under the action of the Hermitian curvature flows. In particular we study two concrete classes of examples, namely, linear Hopf manifolds and six-dimensional Calabi-Yau solvmanifolds with holomorphically-trivial canonical bundle. From these examples we identify some HCFs which do not preserve Bismut-Griffiths-non-negativity.


title: UAVs Path Deviation Attacks: Survey and Research Challenges

abstract: Recently, Unmanned Aerial Vehicles (UAVs) are employed for a plethora of civilian applications. Such flying vehicles can accomplish tasks under the pilot's eyesight within the range of a remote controller, or autonomously according to a certain pre-loaded path configuration. Different path deviation attacks can be performed by malicious users against UAVs. We classify such attacks and the relative defenses based on the UAV's flight mode, i.e., (i) First Person View (FPV), (ii) civilian Global Navigation Satellite System based (GNSS), and (iii) GNSS "plus" auxiliary technologies (GNSS+), and on the multiplicity, i.e., (i) Single UAV, and (ii) Multiple UAVs. We found that very little has been done to secure the FPV flight mode against path deviation. In GNSS mode, spoofing is the most worrisome attack. The best defense against spoofing seems to be redundancy, such as adding vision chips to single UAV or using multiple arranged UAVs. No specific attacks and defenses have been found in literature for GNSS+ or for UAVs moving in group without a pre-ordered arrangement. These aspects require further investigation.


title: Topological surface currents accessed through reversible hydrogenation of the three-dimensional bulk

abstract: Hydrogen, the smallest and most abundant element in nature, can be efficiently incorporated within a solid and drastically modify its electronic state - it has been known to induce novel magnetoelectric effects in complex perovskites and modulate insulator-to-metal transition in a correlated Mott oxide. Here we demonstrate that hydrogenation resolves an outstanding challenge in chalcogenide classes of three-dimensional (3D) topological insulators and magnets - the control of intrinsic bulk conduction that denies access to quantum surface transport. With electrons donated by a reversible binding of H+ ions to Te(Se) chalcogens, carrier densities are easily changed by over 10^20 cm^-3, allowing tuning the Fermi level into the bulk bandgap to enter surface/edge current channels. The hydrogen-tuned topological materials are stable at room temperature and tunable disregarding bulk size, opening a breadth of platforms for harnessing emergent topological states.


title: Performance of the CMS Zero Degree Calorimeters in pPb collisions at the LHC

abstract: The two Zero Degree Calorimeters (ZDCs) of the CMS experiment are located at $\pm 140~$m from the collision point and detect neutral particles in the $|\eta| > 8.3$ pseudorapidity region. This paper presents a study on the performance of the ZDC in the 2016 pPb run. The response of the detectors to ultrarelativistic neutrons is studied using in-depth Monte Carlo simulations. A method of signal extraction based on template fits is presented, along with a dedicated calibration procedure. A deconvolution technique for the correction of overlapping collision events is discussed.


title: Crack occurrence in bodies with gradient polyconvex energies

abstract: Energy minimality selects among possible configurations of a continuous body with and without cracks those compatible with assigned boundary conditions of Dirichlet-type. Crack paths are described in terms of curvature varifolds so that we consider both \textquotedblleft phase" (cracked or non-cracked) and crack orientation. The energy considered is gradient polyconvex: it accounts for relative variations of second-neighbor surfaces and pressure-confinement effects. We prove the existence of minimizers for such an energy. They are pairs of deformations and varifolds. The former ones are taken to be $SBV$ maps satisfying an impenetrability condition. Their jump set is constrained to be in the varifold support.


title: Generalizations of Samuel's criteria for a ring to be a unique factorization domain

abstract: We give several criteria for a ring to be a UFD including generalizations of some criteria due to P. Samuel. These criteria are applied to construct, for any field k, (1) a Z-graded non-noetherian rational UFD of dimension three over k, and (2) k-affine rational UFDs defined by trinomial relations.


title: A High Pressure Time Projection Chamber with Optical Readout

abstract: Measurements of proton-nucleus scattering and high resolution neutrino-nucleus interaction imaging are key to reduce neutrino oscillation systematic uncertainties in future experiments. A High Pressure Time Projection Chamber (HPTPC) prototype has been constructed and operated at Royal Holloway University of London and CERN as a first step in the development of a HPTPC capable of performing these measurements as part of a future long-baseline neutrino oscillation experiment such as the Deep Underground Neutrino Experiment. In this paper we describe the design and operation of the prototype HPTPC with an argon based gas mixture. We report on the successful hybrid charge and optical readout, using four CCD cameras, of signals from Am-241 sources.


title: A Novel Nano Tomography Setup for Material Science and Engineering Applications

abstract: In a comprehensive study on several samples we demonstrate for our laboratory-based computed tomography system resolutions down to 150nm. The achieved resolution is validated by imaging com-mon test structures in 2D and Fourier Shell Correlation of 3D volumes. As representative application examples from nowadays material research, we show metallization processes in multilayer integrated circuits, ageing in lithium battery electrodes, and volumetric of metallic sub-micrometer fillers of com-posites. Thus, our laboratory system provides the unique possibility to image non-destructively struc-tures in the range of hundred nanometers, even for high density materials.


title: Bayesian Quadrature on Riemannian Data Manifolds

abstract: Riemannian manifolds provide a principled way to model nonlinear geometric structure inherent in data. A Riemannian metric on said manifolds determines geometry-aware shortest paths and provides the means to define statistical models accordingly. However, these operations are typically computationally demanding. To ease this computational burden, we advocate probabilistic numerical methods for Riemannian statistics. In particular, we focus on Bayesian quadrature (BQ) to numerically compute integrals over normal laws on Riemannian manifolds learned from data. In this task, each function evaluation relies on the solution of an expensive initial value problem. We show that by leveraging both prior knowledge and an active exploration scheme, BQ significantly reduces the number of required evaluations and thus outperforms Monte Carlo methods on a wide range of integration problems. As a concrete application, we highlight the merits of adopting Riemannian geometry with our proposed framework on a nonlinear dataset from molecular dynamics.


title: Explicit Basis Function Kernel Methods for Cloud Segmentation in Infrared Sky Images

abstract: Photovoltaic systems are sensitive to cloud shadow projection, which needs to be forecasted to reduce the noise impacting the intra-hour forecast of global solar irradiance. We present a comparison between different kernel discriminative models for cloud detection. The models are solved in the primal formulation to make them feasible in real-time applications. The performances are compared using the j-statistic. The infrared cloud images have been preprocessed to remove debris, which increases the performance of the analyzed methods. The use of neighboring features of the pixels also leads to a performance improvement. Discriminative models solved in the primal yield a dramatically lower computing time along with high performance in the segmentation.


title: Non-local reparametrization action in coupled Sachdev-Ye-Kitaev models

abstract: We continue the investigation of coupled Sachdev-Ye-Kitaev(SYK) models without Schwarzian action dominance. Like the original SYK, at large N and low energies these models have an approximate reparametrization symmetry. However, the dominant action for reparametrizations is non-local due to the presence of irrelevant local operator with small conformal dimension. We semi-analytically study different thermodynamic properties and the 4-point function and demonstrate that they significantly differ from the Schwarzian prediction. However, the residual entropy and maximal chaos exponent are the same as in Majorana SYK. We also discuss chain models and finite N corrections.


title: A Critical Look at the Consistency of Causal Estimation With Deep Latent Variable Models

abstract: Using deep latent variable models in causal inference has attracted considerable interest recently, but an essential open question is their ability to yield consistent causal estimates. While they have demonstrated promising results and theory exists on some simple model formulations, we also know that causal effects are not even identifiable in general with latent variables. We investigate this gap between theory and empirical results with analytical considerations and extensive experiments under multiple synthetic and real-world data sets, using the causal effect variational autoencoder (CEVAE) as a case study. While CEVAE seems to work reliably under some simple scenarios, it does not estimate the causal effect correctly with a misspecified latent variable or a complex data distribution, as opposed to its original motivation. Hence, our results show that more attention should be paid to ensuring the correctness of causal estimates with deep latent variable models.


title: Sharing pizza in n dimensions

abstract: We introduce and prove the $n$-dimensional Pizza Theorem: Let $\mathcal{H}$ be a hyperplane arrangement in $\mathbb{R}^{n}$. If $K$ is a measurable set of finite volume, the {pizza quantity} of $K$ is the alternating sum of the volumes of the regions obtained by intersecting $K$ with the arrangement $\mathcal{H}$. We prove that if $\mathcal{H}$ is a Coxeter arrangement different from $A_{1}^{n}$ such that the group of isometries $W$ generated by the reflections in the hyperplanes of $\mathcal{H}$ contains the map $-\mathrm{id}$, and if $K$ is a translate of a convex body that is stable under $W$ and contains the origin, then the pizza quantity of $K$ is equal to zero. Our main tool is an induction formula for the pizza quantity involving a subarrangement of the restricted arrangement on hyperplanes of $\mathcal{H}$ that we call the {even restricted arrangement}. More generally, we prove that for a class of arrangements that we call {even} (this includes the Coxeter arrangements above) and for a {sufficiently symmetric} set $K$, the pizza quantity of $K+a$ is polynomial in $a$ for $a$ small enough, for example if $K$ is convex and $0\in K+a$. We get stronger results in the case of balls, more generally, convex bodies bounded by quadratic hypersurfaces. For example, we prove that the pizza quantity of the ball centered at $a$ having radius $R\geq\|a\|$ vanishes for a Coxeter arrangement $\mathcal{H}$ with $|\mathcal{H}|-n$ an even positive integer. We also prove the Pizza Theorem for the surface volume: When $\mathcal{H}$ is a Coxeter arrangement and $|\mathcal{H}| - n$ is a nonnegative even integer, for an $n$-dimensional ball the alternating sum of the $(n-1)$-dimensional surface volumes of the regions is equal to zero.


title: Robust White Matter Hyperintensity Segmentation on Unseen Domain

abstract: Typical machine learning frameworks heavily rely on an underlying assumption that training and test data follow the same distribution. In medical imaging which increasingly begun acquiring datasets from multiple sites or scanners, this identical distribution assumption often fails to hold due to systematic variability induced by site or scanner dependent factors. Therefore, we cannot simply expect a model trained on a given dataset to consistently work well, or generalize, on a dataset from another distribution. In this work, we address this problem, investigating the application of machine learning models to unseen medical imaging data. Specifically, we consider the challenging case of Domain Generalization (DG) where we train a model without any knowledge about the testing distribution. That is, we train on samples from a set of distributions (sources) and test on samples from a new, unseen distribution (target). We focus on the task of white matter hyperintensity (WMH) prediction using the multi-site WMH Segmentation Challenge dataset and our local in-house dataset. We identify how two mechanically distinct DG approaches, namely domain adversarial learning and mix-up, have theoretical synergy. Then, we show drastic improvements of WMH prediction on an unseen target domain.


title: Coupled Sachdev-Ye-Kitaev models without Schwartzian dominance

abstract: We argue that in certain class of coupled Sachdev-Ye-Kitaev(SYK) models the low energy physics at large N is governed by a non-local action rather than the Schwartzian action. We present a partial analytic and extensive numerical evidence for this. We find that these models are maximally chaotic and have the same residual entropy as Majorana SYK. However, thermodynamic quantities, such as heat capacity and diffusion constant are different.


title: Barriers for recent methods in geodesic optimization

abstract: We study a class of optimization problems including matrix scaling, matrix balancing, multidimensional array scaling, operator scaling, and tensor scaling that arise frequently in theory and in practice. Some of these problems, such as matrix and array scaling, are convex in the Euclidean sense, but others such as operator scaling and tensor scaling are geodesically convex on a different Riemannian manifold. Trust region methods, which include box-constrained Newton's method, are known to produce high precision solutions very quickly for matrix scaling and matrix balancing (Cohen et. al., FOCS 2017, Allen-Zhu et. al. FOCS 2017), and result in polynomial time algorithms for some geodesically convex problems like operator scaling (Garg et. al. STOC 2018, B\"urgisser et. al. FOCS 2019). One is led to ask whether these guarantees also hold for multidimensional array scaling and tensor scaling. We show that this is not the case by exhibiting instances with exponential diameter bound: we construct polynomial-size instances of 3-dimensional array scaling and 3-tensor scaling whose approximate solutions all have doubly exponential condition number. Moreover, we study convex-geometric notions of complexity known as margin and gap, which are used to bound the running times of all existing optimization algorithms for such problems. We show that margin and gap are exponentially small for several problems including array scaling, tensor scaling and polynomial scaling. Our results suggest that it is impossible to prove polynomial running time bounds for tensor scaling based on diameter bounds alone. Therefore, our work motivates the search for analogues of more sophisticated algorithms, such as interior point methods, for geodesically convex optimization that do not rely on polynomial diameter bounds.


title: Development of sieve of Eratosthenes and sieve of Sundaram's proof

abstract: We make two algorithms that generate all prime numbers up to a given limit, they are a development of sieve of Eratosthenes algorithm, we use two formulas to achieve this development, where all the multiples of prime number 2 are eliminated in the first formula, and all the multiples of prime numbers 2 and 3 are eliminated in the second formula. Using the first algorithm we proof sieve of Sundaram's algorithm, then we improve it to be more efficient prime generating algorithm. We will show the difference in performance between all the algorithms we will make and sieve of Eratosthenes algorithm in terms of run time.


title: Higher algebra of $A_\infty$ and $\Omega B As$-algebras in Morse theory I

abstract: Elaborating on works by Abouzaid and Mescher, we prove that for a Morse function on a smooth compact manifold, its Morse cochain complex can be endowed with an $\Omega B As$-algebra structure by counting moduli spaces of perturbed Morse gradient trees. This rich structure descends to its already known $A_\infty$-algebra structure. We then introduce the notion of $\Omega B As$-morphism between two $\Omega B As$-algebras and prove that given two Morse functions, one can construct an $\Omega B As$-morphism between their associated $\Omega B As$-algebras by counting moduli spaces of two-colored perturbed Morse gradient trees. This morphism induces a standard $A_\infty$-morphism between the induced $A_\infty$-algebras. We work with integer coefficients, and provide to this extent a detailed account on the sign conventions for $A_\infty$ (resp. $\Omega B As$)-algebras and $A_\infty$ (resp. $\Omega B As$)-morphisms, using polytopes (resp. moduli spaces) which explicitly realize the dg-operadic objects encoding them. Our proofs also involve showing at the level of polytopes that an $\Omega B As$-morphism between $\Omega B As$-algebras naturally induces an $A_\infty$-morphism between $A_\infty$-algebras. This paper comes in particular with a short survey on operads, $A_\infty$-algebras and $A_\infty$-morphisms, the associahedra and the multiplihedra. All the details on transversality, gluing maps, signs and orientations for the moduli spaces defining the algebraic structures on the Morse cochains are thorougly carried out. It moreover lays the basis for a second article in which we solve the problem of finding a satisfactory homotopic notion of higher morphisms between $A_\infty$-algebras and between $\Omega B As$-algebras, and show how this higher algebra of $A_\infty$ and $\Omega B As$-algebras naturally arises in the context of Morse theory.


title: Responsibility and verification: Importance value in temporal logics

abstract: We aim at measuring the influence of the nondeterministic choices of a part of a system on its ability to satisfy a specification. For this purpose, we apply the concept of Shapley values to verification as a means to evaluate how important a part of a system is. The importance of a component is measured by giving its control to an adversary, alone or along with other components, and testing whether the system can still fulfill the specification. We study this idea in the framework of model-checking with various classical types of linear-time specification, and propose several ways to transpose it to branching ones. We also provide tight complexity bounds in almost every case.


title: Leveraging Artificial Intelligence to Analyze the COVID-19 Distribution Pattern based on Socio-economic Determinants

abstract: The spatialization of socioeconomic data can be used and integrated with other sources of information to reveal valuable insights. Such data can be utilized to infer different variations, such as the dynamics of city dwellers and their spatial and temporal variability. This work focuses on such applications to explore the underlying association between socioeconomic characteristics of different geographical regions in Dublin, Ireland, and the number of confirmed COVID cases in each area. Our aim is to implement a machine learning approach to identify demographic characteristics and spatial patterns. Spatial analysis was used to describe the pattern of interest in Electoral Divisions (ED), which are the legally defined administrative areas in the Republic of Ireland for which population statistics are published from the census data. We used the most informative variables of the census data to model the number of infected people in different regions at ED level. Seven clusters detected by implementing an unsupervised neural network method. The distribution of people who have contracted the virus was studied.


title: End-to-end Audio-visual Speech Recognition with Conformers

abstract: In this work, we present a hybrid CTC/Attention model based on a ResNet-18 and Convolution-augmented transformer (Conformer), that can be trained in an end-to-end manner. In particular, the audio and visual encoders learn to extract features directly from raw pixels and audio waveforms, respectively, which are then fed to conformers and then fusion takes place via a Multi-Layer Perceptron (MLP). The model learns to recognise characters using a combination of CTC and an attention mechanism. We show that end-to-end training, instead of using pre-computed visual features which is common in the literature, the use of a conformer, instead of a recurrent network, and the use of a transformer-based language model, significantly improve the performance of our model. We present results on the largest publicly available datasets for sentence-level speech recognition, Lip Reading Sentences 2 (LRS2) and Lip Reading Sentences 3 (LRS3), respectively. The results show that our proposed models raise the state-of-the-art performance by a large margin in audio-only, visual-only, and audio-visual experiments.


title: MG1-688432: A Peculiar Variable System

abstract: The short period variable star MG1-688432 has been discovered to exhibit occasional extremely high energy optical outbursts as high as 10^31 joules. Outbursts are typically of several hours duration. These events are often highly structured, resembling sequential associated releases of energy. Twenty years of time sequence photometry is presented, indicating a basic sinusoidal light curve of mean period 6.65d, with some phase shifting and long-term temporal trends in amplitude and mean brightness. Spectroscopy reveals a peculiar star, best resembling a K3 subgiant that has evolved off the main sequence moderately red-ward of the giant branch. Spectroscopic and radial velocity analyses indicate a binary system orbiting its barycenter with an unseen companion to the K3IV primary. This is not an eclipsing system with the inclination of the orbit precluding eclipse by the secondary. The system is at a distance of 1.5kpc and analysis of GAIA observations leads to the conclusion that the HR diagram position of MG1-688432 is established by an intrinsic feature of the system, most likely either the stellar evolutionary state of the observed star or the presence of small (non-gray) dust within the system. Two mechanisms that might give rise to the system are 1) impacts with tidally disrupted planetary debris, and 2) magnetically induced chromospheric activity. An intriguing idea that requires further investigation suggests that the unseen companion is perhaps a white dwarf star which has encountered a planet and tidally shredded it to produce a debris and dust veil that modulates the brightness of the primary.


title: Leveraging Artificial Intelligence to Analyze Citizens' Opinions on Urban Green Space

abstract: Continued population growth and urbanization is shifting research to consider the quality of urban green space over the quantity of these parks, woods, and wetlands. The quality of urban green space has been hitherto measured by expert assessments, including in-situ observations, surveys, and remote sensing analyses. Location data platforms, such as TripAdvisor, can provide people's opinion on many destinations and experiences, including UGS. This paper leverages Artificial Intelligence techniques for opinion mining and text classification using such platform's reviews as a novel approach to urban green space quality assessments. Natural Language Processing is used to analyze contextual information given supervised scores of words by implementing computational analysis. Such an application can support local authorities and stakeholders in their understanding of and justification for future investments in urban green space.


title: Lattice QCD equation of state at finite chemical potential from an alternative expansion scheme

abstract: Taylor expansion of the equation of state of QCD suffers from shortcomings at chemical potentials $\mu_B \geq (2-2.5)T$. First, one faces difficulties inherent in performing such an expansion with a limited number of coefficients; second, higher order coefficients determined from lattice calculations suffer from a poor signal-to-noise ratio. In this work, we present a novel scheme for extrapolating the equation of state of QCD to finite, real chemical potential that can extend its reach further than previous methods. We present continuum extrapolated lattice results for the new expansion coefficients and show the thermodynamic observables up to $\mu_B/T\le3.5$.


title: Numerical investigation of Mach number consistent Roe solvers for the Euler equations of gas dynamics

abstract: While traditional approaches to prevent the carbuncle phenomenon in gas dynamics simulations increase the viscosity on entropy and shear waves near shocks, it was quite recently suggested to instead decrease the viscosity on the acoustic waves for low Mach numbers. The goal is to achieve what, in this paper, we call Mach number consistency: for all waves, the numerical viscosity decreases with the same order of the Mach number when the Mach number tends to zero. We take the simple approach that was used for the proof of concept together with the simple model for the increased numerical viscosity on linear waves and investigate the possibilities of combining both in an adaptive manner while locally maintaining Mach number consistency.


title: Qualifying Software Engineers Undergraduates in DevOps -- Challenges of Introducing Technical and Non-technical Concepts in a Project-oriented Course

abstract: The constant changes in the software industry, practices, and methodologies impose challenges to teaching and learning current software engineering concepts and skills. DevOps is particularly challenging because it covers technical concepts, such as pipeline automation, and non-technical ones, such as team roles and project management. The present study investigates a course setup to introduce these concepts to software engineering undergraduates. We designed the course by employing coding to associate DevOps concepts to Agile, Lean, and Open source practices and tools. We present the main aspects of this project-oriented DevOps course, with 240 students enrolled in it since its first offering in 2016. We conducted an empirical study, with both a quantitative and qualitative analysis, to evaluate this project-oriented course setup. We collected the data from the projects repository and students perceptions from a questionnaire. We mined 148 repositories (corresponding to 72 projects) and obtained 86 valid responses to the questionnaire. We also mapped the concepts which are more challenging to students learn from experience. The results evidence that first-hand experience facilitates the comprehension of DevOps concepts and enriches classes discussions. We present a set of lessons learned, which may help professors better design and conduct project-oriented courses to cover DevOps concepts.


title: Potential singularity formation of incompressible axisymmetric Euler equations with degenerate viscosity coefficients

abstract: In this paper, we present strong numerical evidences that the incompressible axisymmetric Euler equations with degenerate viscosity coefficients and smooth initial data of finite energy develop a potential finite-time locally self-similar singularity at the origin. An important feature of this potential singularity is that the solution develops a two-scale traveling wave that travels towards the origin. The two-scale feature is characterized by the scaling property that the center of the traveling wave is located at a ring of radius $O((T-t)^{1/2})$ surrounding the symmetry axis while the thickness of the ring collapses at a rate $O(T-t)$. The driving mechanism for this potential singularity is due to an antisymmetric vortex dipole that generates a strong shearing layer in both the radial and axial velocity fields. Without the viscous regularization, the $3$D Euler equations develop a sharp front and some shearing instability in the far field. On the other hand, the Navier-Stokes equations with a constant viscosity coefficient regularize the two-scale solution structure and do not develop a finite-time singularity for the same initial data.


title: Vortex structures and electron beam dynamics in magnetized plasma

abstract: The formation of vortex structures at reflection of electron beam from the double layer of the Jupiter ionosphere is investigated in this paper. And also the influence of these vortex structures on the formation of dense upward electron fluxes, accelerated by the double layer potential along the Io flux tube is studied. Then a phase transition to the cyclotron superradiance mode becomes possible for these electrons. The conditions of the vortex perturbations formation are considered. The nonlinear equation is found that describes the vortex dynamics of electrons and its consequences are studied.


title: Bayesian Uncertainty Estimation of Learned Variational MRI Reconstruction

abstract: Recent deep learning approaches focus on improving quantitative scores of dedicated benchmarks, and therefore only reduce the observation-related (aleatoric) uncertainty. However, the model-immanent (epistemic) uncertainty is less frequently systematically analyzed. In this work, we introduce a Bayesian variational framework to quantify the epistemic uncertainty. To this end, we solve the linear inverse problem of undersampled MRI reconstruction in a variational setting. The associated energy functional is composed of a data fidelity term and the total deep variation (TDV) as a learned parametric regularizer. To estimate the epistemic uncertainty we draw the parameters of the TDV regularizer from a multivariate Gaussian distribution, whose mean and covariance matrix are learned in a stochastic optimal control problem. In several numerical experiments, we demonstrate that our approach yields competitive results for undersampled MRI reconstruction. Moreover, we can accurately quantify the pixelwise epistemic uncertainty, which can serve radiologists as an additional resource to visualize reconstruction reliability.


title: What helped, and what did not? An Evaluation of the Strategies to Improve Continuous Integration

abstract: Continuous integration (CI) is a widely used practice in modern software engineering. Unfortunately, it is also an expensive practice - Google and Mozilla estimate their CI systems in millions of dollars. There are a number of techniques and tools designed to or having the potential to save the cost of CI or expand its benefit - reducing time to feedback. However, their benefits in some dimensions may also result in drawbacks in others. They may also be beneficial in other scenarios where they are not designed to help. In this paper, we perform the first exhaustive comparison of techniques to improve CI, evaluating 14 variants of 10 techniques using selection and prioritization strategies on build and test granularity. We evaluate their strengths and weaknesses with 10 different cost and time-tofeedback saving metrics on 100 real-world projects. We analyze the results of all techniques to understand the design decisions that helped different dimensions of benefit. We also synthesized those results to lay out a series of recommendations for the development of future research techniques to advance this area.


title: Triangulating metric surfaces

abstract: We prove that any length metric space homeomorphic to a surface may be decomposed into non-overlapping convex triangles of arbitrarily small diameter. This generalizes a previous result of Alexandrov--Zalgaller for surfaces of bounded curvature.


title: Numerical analysis of a model of two phase compressible fluid flow

abstract: We consider a model of a binary mixture of two immiscible compressible fluids. We propose a numerical scheme and discuss its basic properties: Stability, consistency, convergence. The convergence is established via the method of generalized weak solutions combined with the weak-strong uniqueness principle.


title: Oxygen magnetic polarization, nodes in spin density, and zigzag spin order in oxides

abstract: Recent DFT calculations for Ba2CoO4 (BCO) and neutron scattering experiments for SrRuO3 (SRO) have shown that oxygen develops a magnetic polarization. Moreover, DFT calculations for these compounds also unveiled unexpected nodes in the spin density, both along Co-O and Ru-O. For BCO, the overall antiferromagnetic state in its triangular lattice contains unusual zigzag spin patterns. Here, using simple model calculations supplemented by DFT we explain and extend these results. We predict that ligands that in principle should be spinless, such as O$^{2-}$, will develop a net polarization when they act as electronic bridges between transition metal (TM) spins ferromagnetically ordered, regardless of the number of intermediate ligand atoms. The reason is the hybridization between atoms and mobility of the electrons with spins opposite to those of the closest TM atoms. Moreover, for bonds with TMs antiferromagnetically ordered, counterintuitively our calculations show that oxygens should also have a net magnetization for the super-super-exchange cases TM-O-O-TM while for only one oxygen, as in Cu-O-Cu, the O-polarization should cancel. Our simple model also allows us to explain the presence of nodes based on the antibonding character of the dominant singly occupied molecular orbitals along the TM-O bonds. Finally, the zigzag pattern order becomes the ground state mainly due to the influence of the Hubbard $U$, that creates the moments, in combination with a robust easy-axis anisotropy that suppresses the competing 120$^{\circ}$ degree antiferromagnetic order of a triangular lattice. Our predictions are generic and should be applicable to any other compound with characteristics similar to those of BCO and SRO.


title: Some simplifications in the proof of the Sims conjecture

abstract: We prove an elementary lemma concerning primitive amalgams and use it to greatly simplify the proof of the Sims conjecture in the case of almost simple groups.


title: Online voluntary mentoring: Optimising the assignment of students and mentors

abstract: After the closure of the schools in Hungary from March 2020 due to the pandemic, many students were left at home with no or not enough parental help for studying, and in the meantime some people had more free time and willingness to help others in need during the lockdown. In this paper we describe the optimisation aspects of a joint NGO project for allocating voluntary mentors to students using a web-based coordination mechanism. The goal of the project has been to form optimal pairs and study groups by taking into the preferences and the constraints of the participants. In this paper we present the optimisation concept, and the integer programming techniques used for solving the allocation problems. Furthermore, we conducted computation simulations on real and generated data for evaluate the performance of this dynamic matching scheme under different parameter settings.


title: Computations of the Comodule Structures of the Chow rings of Flag Varieties

abstract: Let $G$ be a connected reductive group, and $G/B$ be its flag variety. Let $\pi:G\to G/B$ be the natural projection. In this paper, we developed an algorithm to describe the map $\pi^* :\operatorname{CH}^*(G/B;\mathbb{F}_p)\longrightarrow \operatorname{CH}^*(G;\mathbb{F}_p)$ in terms of Schubert cells. Taking advantage of the Pieri rule, we give an explicit formula for $A$-type, $C$-type, $G_2$, $F_4$ of the cohomology map $\pi^* :\operatorname{CH}^*(G/B;\mathbb{F}_p)\longrightarrow \operatorname{CH}^*(G;\mathbb{F}_p)$, and some partial result of $\pi^*$ is given for $E_6$ and $E_7$. Denote the group action map $\mu:G\times G/B\to G/B$, we also give an explicit formula for $A$-type, $C$-type, $G_2$, $F_4$ of the cohomology map $\mu^*: \operatorname{CH}^*(G/B;\mathbb{F}_p)\longrightarrow \operatorname{CH}^*(G\times G/B;\mathbb{F}_p)$.


title: Proof complexity of positive branching programs

abstract: We investigate the proof complexity of systems based on positive branching programs, i.e. non-deterministic branching programs (NBPs) where, for any 0-transition between two nodes, there is also a 1-transition. Positive NBPs compute monotone Boolean functions, just like negation-free circuits or formulas, but constitute a positive version of (non-uniform) NL, rather than P or NC1, respectively. The proof complexity of NBPs was investigated in previous work by Buss, Das and Knop, using extension variables to represent the dag-structure, over a language of (non-deterministic) decision trees, yielding the system eLNDT. Our system eLNDT+ is obtained by restricting their systems to a positive syntax, similarly to how the 'monotone sequent calculus' MLK is obtained from the usual sequent calculus LK by restricting to negation-free formulas. Our main result is that eLNDT+ polynomially simulates eLNDT over positive sequents. Our proof method is inspired by a similar result for MLK by Atserias, Galesi and Pudl\'ak, that was recently improved to a bona fide polynomial simulation via works of Je\v{r}\'abek and Buss, Kabanets, Kolokolova and Kouck\'y. Along the way we formalise several properties of counting functions within eLNDT+ by polynomial-size proofs and, as a case study, give explicit polynomial-size poofs of the propositional pigeonhole principle.


title: A model for traffic incident prediction using emergency braking data

abstract: This article presents a model for traffic incident prediction. Specifically, we address the fundamental problem of data scarcity in road traffic accident prediction by training our model on emergency braking events instead of accidents. Based on relevant risk factors for traffic accidents and corresponding data categories, we evaluate different options for preprocessing sparse data and different Machine Learning models. Furthermore, we present a prototype implementing a traffic incident prediction model for Germany based on emergency braking data from Mercedes-Benz vehicles as well as weather, traffic and road data, respectively. After model evaluation and optimisation, we found that a Random Forest model trained on artificially balanced (under-sampled) data provided the highest classification accuracy of 85% on the original imbalanced data. Finally, we present our conclusions and discuss further work; from gathering more data over a longer period of time to build stronger classification systems, to addition of internal factors such as the driver's visual and cognitive attention.


title: The role of the electric Bond number in the stability of pasta phases

abstract: The stability of pasta phases in cylindrical and spherical Wigner-Seitz (W-S) cells is examined. The electric Bond number is introduced as the ratio of electric and surface energies. In the case of a charged rod in vacuum, other kinds of instabilities appear in addition to the well known Plateau- Rayleigh mode. For the case of a rod confined in a W-S cell the variety of unstable modes is reduced. It comes from the virial theorem, which bounds the value of the Bond number from above and reduces the role played by electric forces. A similar analysis is done for the spherical W-S cell, where it appears that the inclusion of the virial theorem stabilizes all of the modes.


title: A multiscale model of vascular function in chronic thromboembolic pulmonary hypertension

abstract: Chronic thromboembolic pulmonary hypertension (CTEPH) is caused by recurrent or unresolved pulmonary thromboemboli, leading to perfusion defects and increased arterial wave reflections. CTEPH treatment aims to reduce pulmonary arterial pressure and reestablish adequate lung perfusion, yet patients with distal lesions are inoperable by standard surgical intervention. Instead, these patients undergo balloon pulmonary angioplasty (BPA), a multi-session, minimally invasive surgery that disrupts the thromboembolic material within the vessel lumen using a catheter balloon. However, there still lacks an integrative, holistic tool for identifying optimal target lesions for treatment. To address this insufficiency, we simulate CTEPH hemodynamics and BPA therapy using a multiscale fluid dynamics model. The large pulmonary arterial geometry is derived from a computed tomography (CT) image, whereas a fractal tree represents the small vessels. We model ring- and web-like lesions, common in CTEPH, and simulate normotensive conditions and four CTEPH disease scenarios; the latter includes both large artery lesions and vascular remodeling. BPA therapy is simulated by simultaneously reducing lesion severity in three locations. Our predictions mimic severe CTEPH, manifested by an increase in mean proximal pulmonary arterial pressure above 20 mmHg and prominent wave reflections. Both flow and pressure decrease in vessels distal to the lesions and increase in unobstructed vascular regions. We use the main pulmonary artery (MPA) pressure, a wave reflection index, and a measure of flow heterogeneity to select optimal target lesions for BPA. In summary, this study provides a multiscale, image-to-hemodynamics pipeline for BPA therapy planning for inoperable CTEPH patients.


title: Codifferentials and Quasidifferentials of the Expectation of Nonsmooth Random Integrands and Two-Stage Stochastic Programming

abstract: This work is devoted to an analysis of exact penalty functions and optimality conditions for nonsmooth two-stage stochastic programming problems. To this end, we first study the co-/quasi-differentiability of the expectation of nonsmooth random integrands and obtain explicit formulae for its co- and quasidifferential under some natural assumptions on the integrand. Then we analyse exact penalty functions for a variational reformulation of two-stage stochastic programming problems and obtain sufficient conditions for the global exactness of these functions with two different penalty terms. In the end of the paper, we combine our results on the co-/quasi-differentiability of the expectation of nonsmooth random integrands and exact penalty functions to derive optimality conditions for nonsmooth two-stage stochastic programming problems in terms of codifferentials.


title: Thermodynamic Uncertainty Relation Bounds the Extent of Anomalous Diffusion

abstract: In a finite system driven out of equilibrium by a constant external force the thermodynamic uncertainty relation (TUR) bounds the variance of the conjugate current variable by the thermodynamic cost of maintaining the non-equilibrium stationary state. Here we highlight a new facet of the TUR by showing that it also bounds the time-scale on which a finite system can exhibit anomalous kinetics. In particular, we demonstrate that the TUR bounds subdiffusion in a single file confined to a ring as well as a dragged Gaussian polymer chain even when detailed balance is satisfied. Conversely, the TUR bounds the onset of superdiffusion in the active comb model. Remarkably, the fluctuations in a comb model evolving from a steady state behave anomalously as soon as detailed balance is broken. Our work establishes a link between stochastic thermodynamics and the field of anomalous dynamics that will fertilize further investigations of thermodynamic consistency of anomalous diffusion models.


title: Adversarial Branch Architecture Search for Unsupervised Domain Adaptation

abstract: Unsupervised Domain Adaptation (UDA) is a key issue in visual recognition, as it allows to bridge different visual domains enabling robust performances in the real world. To date, all proposed approaches rely on human expertise to manually adapt a given UDA method (e.g. DANN) to a specific backbone architecture (e.g. ResNet). This dependency on handcrafted designs limits the applicability of a given approach in time, as old methods need to be constantly adapted to novel backbones. Existing Neural Architecture Search (NAS) approaches cannot be directly applied to mitigate this issue, as they rely on labels that are not available in the UDA setting. Furthermore, most NAS methods search for full architectures, which precludes the use of pre-trained models, essential in a vast range of UDA settings for reaching SOTA results. To the best of our knowledge, no prior work has addressed these aspects in the context of NAS for UDA. Here we tackle both aspects with an Adversarial Branch Architecture Search for UDA (ABAS): i. we address the lack of target labels by a novel data-driven ensemble approach for model selection; and ii. we search for an auxiliary adversarial branch, attached to a pre-trained backbone, which drives the domain alignment. We extensively validate ABAS to improve two modern UDA techniques, DANN and ALDA, on three standard visual recognition datasets (Office31, Office-Home and PACS). In all cases, ABAS robustly finds the adversarial branch architectures and parameters which yield best performances.