title: An Agent-Based Modelling Approach to Brain Drain

abstract: The phenomenon of brain drain, that is the emigration of highly skilled people, has many undesirable effects, particularly for developing countries. In this study, an agent-based model is developed to understand the dynamics of such emigration. We hypothesise that skilled people's emigration decisions are based on several factors including the overall economic and social difference between the home and host countries, people's ability and capacity to obtain good jobs and start a life abroad, and the barriers of moving abroad. Furthermore, the social network of individuals also plays a significant role. The model is validated using qualitative and quantitative pattern matching with real-world observations. Sensitivity and uncertainty analyses are performed in addition to several scenario analyses. Linear and random forest response surface models are created to provide quick predictions on the number of emigrants as well as to understand the effect sizes of individual parameters. Overall, the study provides an abstract model where brain drain dynamics can be explored. Findings from the simulation outputs show that future socioeconomic state of the country is more important than the current state, lack of barriers results in a large number of emigrants, and network effects ensue compounding effects on emigration. Upon further development and customisation, future versions can assist in the decision-making of social policymakers regarding brain drain.


title: Clustering multilayer graphs with missing nodes

abstract: Relationship between agents can be conveniently represented by graphs. When these relationships have different modalities, they are better modelled by multilayer graphs where each layer is associated with one modality. Such graphs arise naturally in many contexts including biological and social networks. Clustering is a fundamental problem in network analysis where the goal is to regroup nodes with similar connectivity profiles. In the past decade, various clustering methods have been extended from the unilayer setting to multilayer graphs in order to incorporate the information provided by each layer. While most existing works assume - rather restrictively - that all layers share the same set of nodes, we propose a new framework that allows for layers to be defined on different sets of nodes. In particular, the nodes not recorded in a layer are treated as missing. Within this paradigm, we investigate several generalizations of well-known clustering methods in the complete setting to the incomplete one and prove some consistency results under the Multi-Layer Stochastic Block Model assumption. Our theoretical results are complemented by thorough numerical comparisons between our proposed algorithms on synthetic data, and also on real datasets, thus highlighting the promising behaviour of our methods in various settings.


title: Of Moments and Matching: A Game-Theoretic Framework for Closing the Imitation Gap

abstract: We provide a unifying view of a large family of previous imitation learning algorithms through the lens of moment matching. At its core, our classification scheme is based on whether the learner attempts to match (1) reward or (2) action-value moments of the expert's behavior, with each option leading to differing algorithmic approaches. By considering adversarially chosen divergences between learner and expert behavior, we are able to derive bounds on policy performance that apply for all algorithms in each of these classes, the first to our knowledge. We also introduce the notion of moment recoverability, implicit in many previous analyses of imitation learning, which allows us to cleanly delineate how well each algorithmic family is able to mitigate compounding errors. We derive three novel algorithm templates (AdVIL, AdRIL, and DAeQuIL) with strong guarantees, simple implementation, and competitive empirical performance.


title: High-dimensional estimation of quadratic variation based on penalized realized variance

abstract: In this paper, we develop a penalized realized variance (PRV) estimator of the quadratic variation (QV) of a high-dimensional continuous It\^{o} semimartingale. We adapt the principle idea of regularization from linear regression to covariance estimation in a continuous-time high-frequency setting. We show that under a nuclear norm penalization, the PRV is computed by soft-thresholding the eigenvalues of realized variance (RV). It therefore encourages sparsity of singular values or, equivalently, low rank of the solution. We prove our estimator is minimax optimal up to a logarithmic factor. We derive a concentration inequality, which reveals that the rank of PRV is -- with a high probability -- the number of non-negligible eigenvalues of the QV. Moreover, we also provide the associated non-asymptotic analysis for the spot variance. We suggest an intuitive data-driven bootstrap procedure to select the shrinkage parameter. Our theory is supplemented by a simulation study and an empirical application. The PRV detects about three-five factors in the equity market, with a notable rank decrease during times of distress in financial markets. This is consistent with most standard asset pricing models, where a limited amount of systematic factors driving the cross-section of stock returns are perturbed by idiosyncratic errors, rendering the QV -- and also RV -- of full rank.


title: On the Complexity of Equilibrium Computation in First-Price Auctions

abstract: We consider the problem of computing a (pure) Bayes-Nash equilibrium in the first-price auction with continuous value distributions and discrete bidding space. We prove that when bidders have independent subjective prior beliefs about the value distributions of the other bidders, computing an $\varepsilon$-equilibrium of the auction is PPAD-complete, and computing an exact equilibrium is FIXP-complete. We also provide an efficient algorithm for solving a special case of the problem, for a fixed number of bidders and available bids.


title: Moshpit SGD: Communication-Efficient Decentralized Training on Heterogeneous Unreliable Devices

abstract: Training deep neural networks on large datasets can often be accelerated by using multiple compute nodes. This approach, known as distributed training, can utilize hundreds of computers via specialized message-passing protocols such as Ring All-Reduce. However, running these protocols at scale requires reliable high-speed networking that is only available in dedicated clusters. In contrast, many real-world applications, such as federated learning and cloud-based distributed training, operate on unreliable devices with unstable network bandwidth. As a result, these applications are restricted to using parameter servers or gossip-based averaging protocols. In this work, we lift that restriction by proposing Moshpit All-Reduce - an iterative averaging protocol that exponentially converges to the global average. We demonstrate the efficiency of our protocol for distributed optimization with strong theoretical guarantees. The experiments show 1.3x speedup for ResNet-50 training on ImageNet compared to competitive gossip-based strategies and 1.5x speedup when training ALBERT-large from scratch using preemptible compute nodes.


title: Learning ABCs: Approximate Bijective Correspondence for isolating factors of variation with weak supervision

abstract: Representational learning forms the backbone of most deep learning applications, and the value of a learned representation is intimately tied to its information content regarding different factors of variation. Finding good representations depends on the nature of supervision and the learning algorithm. We propose a novel algorithm that utilizes a weak form of supervision where the data is partitioned into sets according to certain inactive (common) factors of variation which are invariant across elements of each set. Our key insight is that by seeking correspondence between elements of different sets, we learn strong representations that exclude the inactive factors of variation and isolate the active factors that vary within all sets. As a consequence of focusing on the active factors, our method can leverage a mix of set-supervised and wholly unsupervised data, which can even belong to a different domain. We tackle the challenging problem of synthetic-to-real object pose transfer, without pose annotations on anything, by isolating pose information which generalizes to the category level and across the synthetic/real domain gap. The method can also boost performance in supervised settings, by strengthening intermediate representations, as well as operate in practically attainable scenarios with set-supervised natural images, where quantity is limited and nuisance factors of variation are more plentiful.


title: Detectability of gravitational higher order modes in the third-generation era

abstract: Detection of higher order modes of gravitational waves in third-generation (3G) ground-based detectors such as Cosmic Explorer and Einstein Telescope is explored. Using the astrophysical population of binary black holes based on events reported in the second gravitational wave catalog by Laser Interferometer Gravitational Wave Observatory (LIGO) and Virgo (GWTC-2), in conjunction with the Madau-Dickinson model for redshift evolution of the binary black hole mergers, we assess the detectability of these higher order modes using a network consisting of three third-generation detectors. We find that the two subleading modes [(3,3) and (4,4)] can be detected in approximately 30% of the population with a network signal-to-noise ratio of 3 or more, and for nearly 10% of the sources, the five leading modes will be detectable. Besides, a study concerning the effect of binary's mass ratio and its orbital inclination with the observer's line-of-sight in detecting various modes is presented. For a few selected events of the LIGO-Virgo catalog, we identify the modes that would have been detected if a third-generation detector was operational when these events were recorded. We also compute the detectability of higher modes by Voyager and find that only $\sim$ 6 and 2% of the detectable population will have an associated detection of (3,3) and (4,4) modes, respectively. Observing these higher order modes in the 3G era would have a huge impact on the science possible with these detectors ranging from astrophysics and cosmology to testing strong-field gravity.


title: Spatiotemporal Imaging of Thickness-Induced Band Bending Junctions

abstract: Van der Waals materials exhibit naturally passivated surfaces and can form versatile heterostructures, enabling observation of carrier transport mechanisms not seen in three-dimensional materials. Here we report observation of a "band bending junction", a new type of semiconductor homojunction whose surface potential landscape depends solely on a difference in thickness between the two semiconductor regions atop a buried heterojunction interface. Using MoS2 on Au to form a buried heterojunction interface, we find that lateral surface potential differences can arise in MoS2 from the local extent of vertical band bending in thin and thick MoS2 regions. Using scanning ultrafast electron microscopy, we examine the spatiotemporal dynamics of photogenerated charge carriers and find that lateral carrier separation is enabled by a band bending junction, which is confirmed with semiconductor transport simulations. Band bending junctions may therefore enable new electronic and optoelectronic devices in Van der Waals materials that rely on thickness variations rather than doping to separate charge carriers.


title: Anycost GANs for Interactive Image Synthesis and Editing

abstract: Generative adversarial networks (GANs) have enabled photorealistic image synthesis and editing. However, due to the high computational cost of large-scale generators (e.g., StyleGAN2), it usually takes seconds to see the results of a single edit on edge devices, prohibiting interactive user experience. In this paper, we take inspirations from modern rendering software and propose Anycost GAN for interactive natural image editing. We train the Anycost GAN to support elastic resolutions and channels for faster image generation at versatile speeds. Running subsets of the full generator produce outputs that are perceptually similar to the full generator, making them a good proxy for preview. By using sampling-based multi-resolution training, adaptive-channel training, and a generator-conditioned discriminator, the anycost generator can be evaluated at various configurations while achieving better image quality compared to separately trained models. Furthermore, we develop new encoder training and latent code optimization techniques to encourage consistency between the different sub-generators during image projection. Anycost GAN can be executed at various cost budgets (up to 10x computation reduction) and adapt to a wide range of hardware and latency requirements. When deployed on desktop CPUs and edge devices, our model can provide perceptually similar previews at 6-12x speedup, enabling interactive image editing. The code and demo are publicly available: https://github.com/mit-han-lab/anycost-gan.


title: Super-Resolution DOA Estimation for Wideband Signals using Arbitrary Linear Arrays without Focusing Matrices

abstract: We focus on developing an effective Direction Of Arrival (DOA) estimation method for wideband sources based on the gridless sparse concept. Previous coherent methods have been designed by dividing wideband frequencies into a few subbands which are transferred to a reference subband using focusing matrices. In this work, as opposed to the previous techniques, we propose a convex optimization problem that leads to an accurate wideband DOA estimation method with no need for any focusing matrix. Moreover, in this method, no initial DOA estimates are required and it can be used for any arbitrary linear arrays. Numerical simulations show that in comparison to some well-known techniques, the proposed method generates outstanding accuracy and better robustness to noise. The effectiveness of the method is also verified in presence of close adjacent sources.


title: Cosmic voids and induced hyperbolicity. II. Sensitivity to void/wall scales

abstract: Cosmic voids as typical under-density regions in the large scale Universe are known for their hyperbolic properties as an ability to deviate the photon beams. The under-density then is acting as the negative curvature in the hyperbolic spaces. The hyperbolicity of voids has to lead to distortion in the statistical analysis at galactic surveys. We reveal the sensitivity of the hyperbolicity and hence of the distortion with respect to the ratio of void/wall scales which are observable parameters. This provides a principal possibility to use the distortion in the galactic surveys in revealing the line-of sight number of cosmic voids and their characteristic scales.


title: Geodesic orbit metrics in a class of homogeneous bundles over quaternionic Stiefel manifolds

abstract: Geodesic orbit spaces (or g.o. spaces) are defined as those homogeneous Riemannian spaces $(M=G/H,g)$ whose geodesics are orbits of one-parameter subgroups of $G$. The corresponding metric $g$ is called a geodesic orbit metric. We study the geodesic orbit spaces of the form $(\Sp(n)/\Sp(n_1)\times \cdots \times \Sp(n_s), g)$, with $0<n_1+\cdots +n_s\leq n$. Such spaces include spheres, quaternionic Stiefel manifolds, Grassmann manifolds and quaternionic flag manifolds. The present work is a contribution to the study of g.o. spaces $(G/H,g)$ with $H$ semisimple.


title: Time granularity impact on propagation of disruptions in a system-of-systems simulation of infrastructure and business networks

abstract: System-of-systems (SoS) approach is often used for simulating disruptions to business and infrastructure system networks allowing for integration of several models into one simulation. However, the integration is frequently challenging as each system is designed individually with different characteristics, such as time granularity. Understanding the impact of time granularity on propagation of disruptions between businesses and infrastructure systems and finding the appropriate granularity for the SoS simulation remain as major challenges. To tackle these, we explore how time granularity, recovery time, and disruption size affect the propagation of disruptions between constituent systems of an SoS simulation. To address this issue, we developed a High Level Architecture (HLA) simulation of 3 networks and performed a series of simulation experiments. Our results revealed that time granularity and especially recovery time have huge impact on propagation of disruptions. Consequently, we developed a model for selecting an appropriate time granularity for an SoS simulation based on expected recovery time. Our simulation experiments show that time granularity should be less than 1.13 of expected recovery time. We identified some areas for future research centered around extending the experimental factors space.


title: Accurate Gaunt factors for non-relativistic quadrupole bremsstrahlung

abstract: The exact result for non-relativistic quadrupole bremsstrahlung in a Coulomb field was established only recently in Pradler & Semmelrock (2020). It requires the evaluation and integration of hypergeometric functions across a wide range of parameters and arguments, which, in practice, is unfeasible. Here we provide a highly accurate tabulation of the Gaunt factor for quadrupole radiation, its thermal average in a Maxwellian plasma, and the associated cooling function over the entire kinematically relevant range. In addition, we provide a simple approximate formula for the emission cross section which works to within a few percent accuracy for all practical purposes. The results can be applied to the scattering of electrons with themselves, for which quadrupole radiation is the dominant process.


title: The Hubble Tension as a Hint of Leptogenesis and Neutrino Mass Generation

abstract: The majoron, a neutrinophilic pseudo-Goldstone boson conventionally arising in the context of neutrino mass models, can damp neutrino free-streaming and inject additional energy density into neutrinos prior to recombination. The combination of these effects for an eV-scale mass majoron has been shown to ameliorate the outstanding $H_0$ tension, however only if one introduces additional dark radiation at the level of $\Delta N_{\rm eff} \sim 0.5$. We show here that models of low-scale leptogenesis can naturally source this dark radiation by generating a primordial population of majorons from the decays of GeV-scale sterile neutrinos in the early Universe. Using a posterior predictive distribution conditioned on Planck2018+BAO data, we show that the value of $H_0$ observed by the SH$_0$ES collaboration is expected to occur at the level of $\sim 10\%$ in the primordial majoron cosmology (to be compared with $\sim 0.1\%$ in the case of $\Lambda$CDM). This insight provides an intriguing connection between the neutrino mass mechanism, the baryon asymmetry of the Universe, and the discrepant measurements of $H_0$.


title: Algorithmically solving the Tadpole Problem

abstract: The extensive computer-aided search applied in [arXiv:2010.10519] to find the minimal charge sourced by the fluxes that stabilize all the (flux-stabilizable) moduli of a smooth K3xK3 compactification uses differential evolutionary algorithms supplemented by local searches. We present these algorithms in detail and show that they can also solve our minimization problem for other lattices. Our results support the Tadpole Conjecture: The minimal charge grows linearly with the dimension of the lattice and, for K3xK3, this charge is larger than allowed by tadpole cancelation. Even if we are faced with an NP-hard lattice-reduction problem at every step in the minimization process, we find that differential evolution is a good technique for identifying the regions of the landscape where the fluxes with the lowest tadpole can be found. We then design a "Spider Algorithm," which is very efficient at exploring these regions and producing large numbers of minimal-tadpole configurations.


title: Reconstructing the Last Major Merger of the Milky Way with the H3 Survey

abstract: Several lines of evidence suggest the Milky Way underwent a major merger at z~2 with a galaxy known as Gaia-Sausage-Enceladus (GSE). Here we use H3 Survey data to argue that GSE entered the Galaxy on a retrograde orbit based on a population of highly retrograde stars with chemistry similar to the largely radial GSE debris. We present the first tailored, high-resolution N-body simulations of the merger. From a grid of ~500 simulations we find a GSE with $M_{*}=5\times10^{8}\ M_{\odot}, M_{\rm{DM}}=2\times10^{11} M_{\odot}$ (a 2.5:1 total mass merger) best matches the H3 data. This simulation shows the retrograde GSE stars are stripped from its outer disk early in the merger before the orbit loses significant angular momentum. Despite being selected purely on angular momenta and radial distributions, this simulation reproduces and explains the following empirical phenomena: (i) the elongated, triaxial shape of the inner halo (axis ratios $10:7.9:4.5$), whose major axis is at ~35{\deg} to the plane and connects GSE's apocenters, (ii) the Hercules-Aquila Cloud & the Virgo Overdensity, which arise due to apocenter pile-up, (iii) the 2 Gyr lag between the quenching of GSE and the truncation of the age distribution of the in-situ halo, which tracks the 2 Gyr gap between the first and final GSE pericenters. We make the following predictions: (i) the inner halo has a "double-break" density profile with breaks at both ~15-18 kpc and 30 kpc, coincident with the GSE apocenters, (ii) the outer halo has retrograde streams containing ~10% of GSE stars awaiting discovery at >30 kpc. The retrograde (radial) GSE debris originates from its outer (inner) disk -- exploiting this trend we reconstruct the stellar metallicity gradient of GSE ($-0.04\pm0.01$ dex $r_{\rm{50}}^{-1}$). These simulations imply GSE delivered ~20% of the Milky Way's present-day dark matter and ~50% of its stellar halo. (ABRIDGED)


title: Constraining the Epoch of Reionization With Highly Dispersed Fast Radio Bursts

abstract: The period in which hydrogen in the intergalactic medium (IGM) is ionized, known as the Epoch of Reionization (EoR) is still poorly understood. The timing and duration of the EoR is expected to be governed by the underlying astrophysics. Furthermore, most models of reionization predict a correlation between the density and ionization field. Here we consider using the mean dispersion measure (DM) of high redshift Fast Radio Bursts (FRBs) as a probe of the underlying astrophysics and morphology of the EoR. To do this, we forecast observational scenarios by building mock data sets of non-repeating FRBs between redshifts $8\leq z \leq 10$. It is assumed that all FRBs have accompanying spectroscopic redshift measurements. We find that samples of 100 high redshift FRBs, in the above mentioned narrow redshift range, can rule out uncorrelated reionization at $68\%$ credibility, while larger samples, $\geq 10^4$ FRBs, can rule out uncorrelated reionization at $95\%$ credibility. We also find 100 high redshift FRBs can rule out scenarios where the Universe is entirely neutral at $z = 10$ with $68\%$ credibility. Further with $\geq 10^5$ FRBs, we can constrain the duration $\Delta z$ of reionization (duration between mean ionized fraction 0.25 to 0.75) to $\Delta z = 2.0^{+0.5}_{-0.4}$, and the midpoint of reionization to $z = 7.8^{+0.4}_{-0.2}$ at $95\%$ credibility.


title: Freeze-In Dark Matter within the Seesaw mechanism

abstract: We show that the minimal Type-I Seesaw mechanism can successfully account for the observed dark matter abundance in the form of a keV sterile neutrino. This population can be produced by the decay of the heavier neutral leptons, with masses above the Higgs mass scale, while they are in thermal equilibrium in the early Universe (freeze-in). Moreover, the implementation of the relevant phenomenological constraints (relic abundance, indirect detection and structure formation) on this model automatically selects a region of the parameter space featuring an approximate lepton number symmetry.


title: Heavy decaying dark matter at future neutrino radio telescopes

abstract: In the next decades, ultra-high-energy neutrinos in the EeV energy range will be potentially detected by next-generation neutrino telescopes. Although their primary goals are to observe cosmogenic neutrinos and to gain insight into extreme astrophysical environments, they can also indirectly probe the nature of dark matter. In this paper, we study the projected sensitivity of up-coming neutrino radio telescopes, such as RNO-G, GRAND and IceCube-gen2 radio array, to decaying dark matter scenarios. We investigate different dark matter decaying channels and masses, from $10^7$ to $10^{15}$ GeV. By assuming the observation of cosmogenic or newborn pulsar neutrinos, we forecast conservative constraints on the lifetime of heavy dark matter particles. We find that these limits are competitive with and highly complementary to previous multi-messenger analyses.


title: Early Science with the Large Millimeter Telescope: Constraining the Gas Fraction of a Compact Quiescent Galaxy at z=1.883

abstract: We present constraints on the dust continuum flux and inferred gas content of a gravitationally lensed massive quiescent galaxy at $z$=1.883$\pm$0.001 using AzTEC 1.1mm imaging with the Large Millimeter Telescope. MRG-S0851 appears to be a prototypical massive compact quiescent galaxy, but has evidence that it experienced a centrally concentrated rejuvenation event in the last 100 Myr (see Akhshik et al. 2020). This galaxy is undetected in the AzTEC image but we calculate an upper limit on the millimeter flux and use this to estimate the H$_2$ mass limit via an empirically calibrated relation that assumes a constant molecular gas-to-dust ratio of 150. We constrain the 3$\sigma$ upper limit of the H$_2$ fraction from the dust continuum in MRG-S0851 to be ${M_{H_2}/M_{\star}}$ $\leq$ 6.8%. MRG-S0851 has a low gas fraction limit with a moderately low sSFR owing to the recent rejuvenation episode, which together results in a relatively short depletion time of $<$0.6 Gyr if no further H$_2$ gas is accreted. Empirical and analytical models both predict that we should have detected molecular gas in MRG-S0851, especially given the rejuvenation episode; this suggests that cold gas and/or dust is rapidly depleted in at least some early quiescent galaxies.


title: Doped Mott Insulators Break $\mathbb Z_2$ Symmetry of a Fermi Liquid: Stability of Strongly Coupled Fixed Points

abstract: Because Fermi liquids are inherently non-interacting states of matter, all electronic levels below the chemical potential are doubly occupied. Consequently, the simplest way of breaking Fermi liquid theory is to engineer a model in which some of those states are singly occupied keeping time-reversal invariance intact. We show that breaking an overlooked local-in-momentum space $\mathbb Z_2$ symmetry of a Fermi liquid does precisely this. As a result, while the Mott transition from a Fermi liquid is correctly believed to obtain without the breaking of any continuous symmetry, a discrete symmetry is broken. This symmetry breaking serves as an organizing principle for Mott physics whether it arises from the tractable Hatsugai-Kohmoto (HK) model or the intractable Hubbard model. That both are controlled by the same fixed point we establish through a renormalization group analysis. An experimental manifestation of this fixed point is the onset of particle-hole asymmetry, a widely observed phenomenon in strongly correlated systems. Theoretically, the singly-occupied region of the spectrum gives rise to a surface of zeros of the single-particle Green function, denoted as the Luttinger surface. Using K-homology, we show that the Bott topological invariant guarantees the stability of this surface to local perturbations. Our proof demonstrates that the strongly coupled fixed point only corresponds to those Luttinger surfaces with co-dimension $p+1$ with $p$ odd. We conclude that the Hubbard and HK models both lie in the same high temperature universality class and are controlled by the broken $\mathbb Z_2$ symmetry quartic fixed point.


title: Strong lens modelling: comparing and combining Bayesian neural networks and parametric profile fitting

abstract: The vast quantity of strong galaxy-galaxy gravitational lenses expected by future large-scale surveys necessitates the development of automated methods to efficiently model their mass profiles. For this purpose, we train an approximate Bayesian convolutional neural network (CNN) to predict mass profile parameters and associated uncertainties, and compare its accuracy to that of conventional parametric modelling for a range of increasingly complex lensing systems. These include standard smooth parametric density profiles, hydrodynamical EAGLE galaxies and the inclusion of foreground mass structures, combined with parametric sources and sources extracted from the Hubble Ultra Deep Field. In addition, we also present a method for combining the CNN with traditional parametric density profile fitting in an automated fashion, where the CNN provides initial priors on the latter's parameters. On average, the CNN achieved errors 19 $\pm$ 22 per cent lower than the traditional method's blind modelling. The combination method instead achieved 27 $\pm$ 11 per cent lower errors over the blind modelling, reduced further to 37 $\pm$ 11 per cent when the priors also incorporated the CNN-predicted uncertainties, with errors also 17 $\pm$ 21 per cent lower than the CNN by itself. While the CNN is undoubtedly the fastest modelling method, the combination of the two increases the speed of conventional fitting alone by factors of 1.73 and 1.19 with and without CNN-predicted uncertainties, respectively. This, combined with greatly improved accuracy, highlights the benefits one can obtain through combining neural networks with conventional techniques in order to achieve an efficient automated modelling approach.


title: On the Zero Point Constant of the Bolometric Correction Scale

abstract: Arbitrariness attributed to the zero point constant of the $V$ band bolometric corrections ($BC_V$) and its relation to "bolometric magnitude of a star ought to be brighter than its visual magnitude" and "bolometric corrections must always be negative" was investigated. The falsehood of the second assertion became noticeable to us after IAU 2015 General Assembly Resolution B2, where the zero point constant of bolometric magnitude scale was decided to have a definite value $C_{Bol}(W)= 71.197~425~...$~. Since the zero point constant of the $BC_V$ scale could be written as $C_2=C_{Bol}-C_V$, where $C_V$ is the zero point constant of the visual magnitudes in the basic definition $BC_V=M_{Bol}-M_V=m_{bol}-m_V$, and $C_{Bol}>C_V$, the zero point constant ($C_2$) of the $BC_V$ scale cannot be arbitrary anymore; rather, it must be a definite positive number obtained from the two definite positive numbers. The two conditions $C_2>0$ and $0<BC_V<C_2$ are also sufficient for $L_V<L$, a similar case to negative $BC_V$ numbers, which means that "bolometric corrections are not always negative". In sum it becomes apparent that the first assertion is misleading causing one to understand bolometric corrections must always be negative, which is not necessarily true.


title: Revisiting the Kepler field with TESS: Improved ephemerides using TESS 2min data

abstract: Up to date planet ephemerides are becoming increasingly important as exoplanet science moves from detecting exoplanets to characterising their architectures and atmospheres in depth. In this work ephemerides are updated for 22 Kepler planets and 4 Kepler planet candidates, constituting all Kepler planets and candidates with sufficient signal to noise in the TESS 2min dataset. A purely photometric method is utilised here to allow ephemeris updates for planets even when they do not posses significant radial velocity data. The obtained ephemerides are of very high precision and at least seven years 'fresher' than archival ephemerides. In particular, significantly reduced period uncertainties for Kepler-411d, Kepler-538b and the candidates K00075.01/K00076.01 are reported. O-C diagrams were generated for all objects, with the most interesting ones discussed here. Updated TTV fits of five known multiplanet systems with significant TTVs were also attempted (Kepler-18, Kepler-25, Kepler-51, Kepler-89, and Kepler-396), however these suffered from the comparative scarcity and dimness of these systems in TESS. Despite these difficulties, TESS has once again shown itself to be an incredibly powerful follow-up instrument as well as a planet-finder in its own right. Extension of the methods used in this paper to the 30min-cadence TESS data and TESS extended mission has the potential to yield updated ephemerides of hundreds more systems in the future.


title: Topological sectors, dimer correlations and monomers from the transfer-matrix solution of the dimer model

abstract: We solve the classical square-lattice dimer model with periodic boundaries and in the presence of a field $\boldsymbol{t}$ that couples to the (vector) flux, by diagonalizing a modified version of Lieb's transfer matrix. After deriving the torus partition function in the thermodynamic limit, we show how the configuration space divides into 'topological sectors' corresponding to distinct values of the flux. Additionally, we demonstrate in general that expectation values are $\boldsymbol{t}$-independent at leading order, and obtain explicit expressions for dimer occupation numbers, dimer-dimer correlation functions and the monomer distribution function. The last of these is expressed as a Toeplitz determinant, whose asymptotic behavior for large monomer separation is tractable using the Fisher-Hartwig conjecture. Our results reproduce those previously obtained using Pfaffian techniques.


title: Implications of the $z>5$ Lyman-$\alpha$ forest for the 21-cm power spectrum from the epoch of reionization

abstract: Our understanding of the intergalactic medium at redshifts $z=5$-$6$ has improved considerably in the last few years due to the discovery of quasars with $z>6$ that enable Lyman-$\alpha$ forest studies at these redshifts. A realisation from this has been that hydrogen reionization could end much later than previously thought, so that large "islands" of cold, neutral hydrogen could exist in the IGM at redshifts $z=5$-$6$. By using radiative transfer simulations of the IGM, we consider the implications of the presence of these neutral hydrogen islands for the 21-cm power spectrum signal and its potential detection by experiments such as HERA, SKA, LOFAR, and MWA. In contrast with previous models of the 21-cm signal, we find that thanks to the late end of reionization the 21-cm power in our simulation continues to be as high as $\Delta^2_{21}=10~\mathrm{mK}^2$ at $k\sim 0.1~h/$cMpc at $z=5$-$6$. This value of the power spectrum is several orders of magnitude higher than that in the conventional models considered in the literature for these redshifts. Such high values of the 21-cm power spectrum should be detectable by HERA and SKA1-LOW in $\sim 1000$ hours, assuming optimistic foreground subtraction. This redshift range is also attractive due to relatively low sky temperature and potentially greater abundance of multiwavelength data.


title: Bayesian probability for leptonic CP phase with strong CP prior

abstract: In a world devoid of axions, the smallness of the strong CP phase can have implications for leptonic CP violation being probed by neutrino experiments. For example, if nature adopted an axionless solution to the strong CP problem, the same symmetries that set the strong CP phase to zero at the tree-level, may also set the leptonic $CP$ phases to zero (mod \pi). This automatically happens in the left-right symmetric model when it is extended minimally to solve the strong CP problem by imposition of both P and CP. In the Nelson Barr solution the needed symmetries can be assigned so that leptonic CP violation is not generated. In the minimal left-right symmetric model with P (and not CP), where the strong CP problem remains, leptonic CP phases radiatively generate the strong CP phase in one loop, and therefore they may be absent (negligibly small) in large regions of parameter space. All these results motivate us to consider a Bayesian prior for leptonic Dirac CP phase \delta_{CP} of the PMNS matrix that has delta function like peaks at CP conserving values of 0 and \pi on top of a uniform distribution. We evaluate the posterior probability distribution for \delta_{CP} using the current global fit to neutrino experiments and find significant enhancements for the probability that \delta_{CP} is at or negligibly close to \pi. We also provide useful tables for the posterior probability considering present and future experimental sensitivities.


title: Mode structure and orbital angular momentum of spatiotemporal optical vortex (STOV) pulses

abstract: We identify a class of modal solutions for spatio-temporal optical vortex (STOV) electromagnetic pulses propagating in dispersive media with orbital angular momentum (OAM) orthogonal to propagation. We find that symmetric STOVs in vacuum can carry half-integer intrinsic orbital angular momentum (OAM); for general asymmetric STOVs in a dispersive medium, the OAM is quantized in integer multiples of a parameter that depends on the STOV symmetry and the group velocity dispersion. Our results suggest that STOVs propagating in dispersive media are accompanied by a polariton-like quasiparticle. The modal theory is in excellent agreement with measurements of free space propagation of STOVs.


title: Quantum routing with fast reversals

abstract: We present methods for implementing arbitrary permutations of qubits under interaction constraints. Our protocols make use of previous methods for rapidly reversing the order of qubits along a path. Given nearest-neighbor interactions on a path of length $n$, we show that there exists a constant $\epsilon \approx 0.034$ such that the quantum routing time is at most $(1-\epsilon)n$, whereas any swap-based protocol needs at least time $n-1$. This represents the first known quantum advantage over swap-based routing methods and also gives improved quantum routing times for realistic architectures such as grids. Furthermore, we show that our algorithm approaches a quantum routing time of $2n/3$ in expectation for uniformly random permutations, whereas swap-based protocols require time $n$ asymptotically. Additionally, we consider sparse permutations that route $k \le n$ qubits and give algorithms with quantum routing time at most $n/3 + O(k^2)$ on paths and at most $2r/3 + O(k^2)$ on general graphs with radius $r$.


title: Better SGD using Second-order Momentum

abstract: We develop a new algorithm for non-convex stochastic optimization that finds an $\epsilon$-critical point in the optimal $O(\epsilon^{-3})$ stochastic gradient and Hessian-vector product computations. Our algorithm uses Hessian-vector products to "correct" a bias term in the momentum of SGD with momentum. This leads to better gradient estimates in a manner analogous to variance reduction methods. In contrast to prior work, we do not require excessively large batch sizes, and are able to provide an adaptive algorithm whose convergence rate automatically improves with decreasing variance in the gradient estimates. We validate our results on a variety of large-scale deep learning architectures and benchmarks tasks.


title: Robustness of Noisy Quantum Networks

abstract: Quantum networks are a new paradigm of complex networks, allowing us to harness networked quantum technologies and to develop a quantum internet. But how robust is a quantum network when its links and nodes start failing? We show that quantum networks based on typical noisy quantum-repeater nodes are prone to discontinuous phase transitions with respect to the random loss of operating links and nodes, abruptly compromising the connectivity of the network, and thus significantly limiting the reach of its operation. Furthermore, we determine the critical quantum-repeater efficiency necessary to avoid this catastrophic loss of connectivity as a function of the network topology, the network size, and the distribution of entanglement in the network. In particular, our results indicate that a scale-free topology is a crucial design principle to establish a robust large-scale quantum internet.

categories: quant-ph cond-mat.dis-nn math-ph math.MP


title: Currents, charges and algebras in exceptional generalised geometry

abstract: A classical $E_{d(d)}$-invariant Hamiltonian formulation of world-volume theories of half-BPS p-branes in type IIb and eleven-dimensional supergravity is proposed, extending known results to $d \leq 6$. It consists of a Hamiltonian, characterised by a generalised metric, and a current algebra constructed s.t. it reproduces the $E_{d(d)}$ generalised Lie derivative. $E_{d(d)}$-covariance necessitates the introduction of so-called charges, specifying the type of p-brane and the choice of section. For p>2, currents of p-branes are generically non-geometric due to the imposition of U-duality, e.g. the M5-currents contain coordinates associated to the M2-momentum. A derivation of the $E_{d(d)}$-invariant current algebra from a canonical Poisson structure is in general not possible. At most, one can derive a current algebra associated to para-Hermitian exceptional geometry. The membrane in the SL(5)-theory is studied in detail. It is shown that in a generalised frame the current algebra is twisted by the generalised fluxes. As a consistency check, the double dimensional reduction from membranes in M-theory to strings in type IIa string theory is performed. Many features generalise to p-branes in SL(p+3) generalised geometries that form building blocks for the $E_{d(d)}$-invariant currents.


title: The HOSTS survey: evidence for an extended dust disk and constraints on the presence of giant planets in the Habitable Zone of $\beta$ Leo

abstract: The young (50-400 Myr) A3V star $\beta$ Leo is a primary target to study the formation history and evolution of extrasolar planetary systems as one of the few stars with known hot ($\sim$1600$^\circ$K), warm ($\sim$600$^\circ$K), and cold ($\sim$120$^\circ$K) dust belt components. In this paper, we present deep mid-infrared measurements of the warm dust brightness obtained with the Large Binocular Telescope Interferometer (LBTI) as part of its exozodiacal dust survey (HOSTS). The measured excess is 0.47\%$\pm$0.050\% within the central 1.5 au, rising to 0.81\%$\pm$0.026\% within 4.5 au, outside the habitable zone of $\beta$~Leo. This dust level is 50 $\pm$ 10 times greater than in the solar system's zodiacal cloud. Poynting-Robertson drag on the cold dust detected by Spitzer and Herschel under-predicts the dust present in the habitable zone of $\beta$~Leo, suggesting an additional delivery mechanism (e.g.,~comets) or an additional belt at $\sim$5.5 au. A model of these dust components is provided which implies the absence of planets more than a few Saturn masses between $\sim$5 au and the outer belt at $\sim$40 au. We also observationally constrain giant planets with the LBTI imaging channel at 3.8~$\mu$m wavelength. Assuming an age of 50 Myr, any planet in the system between approximately 5 au to 50 au must be less than a few Jupiter masses, consistent with our dust model. Taken together, these observations showcase the deep contrasts and detection capabilities attainable by the LBTI for both warm exozodiacal dust and giant exoplanets in or near the habitable zone of nearby stars.


title: Topology and geometry under the nonlinear electromagnetic spotlight

abstract: For many materials, a precise knowledge of their dispersion spectra is insufficient to predict their ordered phases and physical responses. Instead, these materials are classified by the geometrical and topological properties of their wavefunctions. A key challenge is to identify and implement experiments that probe or control these quantum properties. In this review, we describe recent progress in this direction, focusing on nonlinear electromagnetic responses that arise directly from quantum geometry and topology. We give an overview of the field by discussing new theoretical ideas, groundbreaking experiments, and the novel materials that drive them. We conclude by discussing how these techniques can be combined with new device architectures to uncover, probe, and ultimately control novel quantum phases with emergent topological and correlated properties.

categories: cond-mat.mtrl-sci cond-mat.mes-hall cond-mat.str-el cond-mat.supr-con


title: Extraction of the Sivers function from SIDIS, Drell-Yan, and $W^\pm/Z$ boson production data with TMD evolution

abstract: We perform a global fit of the available polarized Semi-Inclusive Deep Inelastic Scattering (SIDIS), polarized pion-induced Drell-Yan (DY) and $W^\pm/Z$ boson production data at N$^3$LO and NNLO accuracy of the Transverse Momentum Dependent (TMD) evolution, and extract the Sivers function for $u$, $d$, $s$ and for sea quarks. The Qiu-Sterman function is determined in a model independent way via the operator product expansion from the extracted Sivers function. The analysis is supplemented by additional studies, such as the estimation of applicability region, the impact of the unpolarized distributions' uncertainties, the universality of the Sivers functions, positivity constraints, the significance of the sign-change relation, and the comparison with the existing extractions


title: Super-resolution Method for Coherent DOA Estimation of Multiple Wideband Sources

abstract: We focus on coherent direction of arrival estimation of wideband sources based on spatial sparsity. This area of research is encountered in many applications such as passive radar, sonar, mining, and communication problems, in which an increasing attention has been devoted to improving the estimation accuracy and robustness to noise. By the development of super-resolution algorithms, narrowband direction of arrival estimation based on gridless sparse algorithms and atomic norm minimization has already been addressed. In this paper, a superresolution based method is proposed for coherent direction of arrival estimation of multiple wideband sources. First, unlike the conventional coherent methods, we develop a new focusing method to map the subband with the largest center frequency to the other ones, which leads to an accurate method with no requirement for initial estimates for DOAs. Then, we introduce an atomic norm problem by defining a new set of atoms and exploiting the signal joint sparsity of different frequency subbands in a continuous spatial domain. This problem is then cast as a semidefinite program, which leads to implementing a new coherent direction of arrival estimation method with higher resolution and more robustness to noise. Our method needs only one single snapshot for each frequency subband, leading to a small number of snapshots for the received wideband signal compared to the other coherent DOA techniques. Numerical simulations show the outperformance of the proposed method compared to the conventional ones.


title: Simulation study of Q statistic with constant weights for testing and estimation of heterogeneity of standardized mean differences in meta-analysis

abstract: Cochran's $Q$ statistic is routinely used for testing heterogeneity in meta-analysis. Its expected value is also used for estimation of between-study variance $\tau^2$. Cochran's $Q$, or $Q_{IV}$, uses estimated inverse-variance weights which makes approximating its distribution rather complicated. As an alternative, we are investigating a new $Q$ statistic, $Q_F$, whose constant weights use only the studies' effective sample sizes. For standardized mean difference as the measure of effect, we study, by simulation, approximations to distributions of $Q_{IV}$ and $Q_F$, as the basis for tests of heterogeneity and for new point and interval estimators of the between-study variance $\tau^2$. These include new DerSimonian-Kacker (2007)-type moment estimators based on the first moment of $Q_F$, and novel median-unbiased estimators of $\tau^2$.


title: Manipulating phonons of a trapped-ion system using optical tweezers

abstract: We propose an experimental architecture where an array of optical tweezers affords site-dependent control over the confining potential of a conventional radio-frequency ion trap. The site-dependent control enables programmable manipulation of phonon modes of ions, with many potential applications in quantum information processing (QIP) and thermodynamics. We describe protocols for programming the array of optical tweezers to attain a set of target phonon modes with high accuracy. We propose applications of such controls in simulating quantum thermodynamics of a particle of programmable effective mass via Jarzynski's equality and improving the efficiency of sympathetic cooling and quantum logic gates in a multi-species ion system of disparate masses. We discuss the required optical parameters in a realistic ion trap system and potential adverse effects of optical tweezers in QIP. Our scheme extends the utility of trapped-ions as a platform for quantum computation and simulation.


title: Do biological constraints impair dendritic computation?

abstract: Computations on the dendritic trees of neurons have important constraints. Voltage dependent conductances in dendrites are not similar to arbitrary direct-current generation, they are the basis for dendritic nonlinearities and they do not allow converting positive currents into negative currents. While it has been speculated that the dendritic tree of a neuron can be seen as a multi-layer neural network and it has been shown that such an architecture could be computationally strong, we do not know if that computational strength is preserved under these biological constraints. Here we simulate models of dendritic computation with and without these constraints. We find that dendritic model performance on interesting machine learning tasks is not hurt by these constraints but may benefit from them. Our results suggest that single real dendritic trees may be able to learn a surprisingly broad range of tasks.


title: The Climate Extended Risk Model (CERM)

abstract: This paper addresses estimates of climate risk embedded within a bank credit portfolio. The proposed Climate Extended Risk Model (CERM) adapts well known credit risk models and makes it possible to calculate incremental credit losses on a loan portfolio that are rooted into physical and transition risks. The paper provides detailed description of the model hypotheses and steps.


title: Counting in Uncountably Categorical Pseudofinite Structures

abstract: We show that every definable subset of an uncountably categorical pseudofinite structure has pseudofinite cardinality which is polynomial (over the rationals) in the size of any strongly minimal subset, with the degree of the polynomial equal to the Morley rank of the subset. From this fact, we show that classes of finite structures whose ultraproducts all satisfy the same uncountably categorical theory are polynomial $R$-mecs as well as $N$-dimensional asymptotic classes, where $N$ is the Morley rank of the theory.


title: External gas accretion provides a fresh gas supply to the active S0 galaxy NGC 5077

abstract: In early type galaxies, externally accreted gas is thought to be the main source of gas replenishment at late times. We use MUSE integral field spectroscopy data to study the active S0 galaxy NGC 5077, known to have disturbed dynamics, indicative of a past external interaction. We confirm the presence of a stellar kinematically distinct core with a diameter of 2.8 kpc, counter-rotating with respect to the main stellar body of the galaxy. We find that the counter-rotating core consists of an old stellar population, not significantly different from the rest of the galaxy. The ionised gas is strongly warped and extends out to 6.5 kpc in the polar direction and in a filamentary structure. The gas dynamics is complex, with significant changes in the position angle as a function of radius. The ionised gas line ratios are consistent with LINER excitation by the AGN both in the nucleus and at kiloparsec scales. We discover a nuclear outflow with projected velocity V ~ 400 km/s, consistent with a hollow outflow cone intersecting the plan of the sky. The properties of the misaligned gas match predictions from numerical simulations of misaligned gas infall after a gas-rich merger. The warp and change in the gas orientation as a function of radius are consistent with gas relaxation due to stellar torques, that are stronger at small radii where the gas aligns faster than in the outer regions, driving gas to the nucleus. The stellar and gas dynamics indicate that NGC 5077 has had at least two external interactions, one that resulted in the formation of the counter-rotating core followed by late time external gas accretion. NGC 5077 illustrates the importance of external interactions in the replenishment of the galaxy gas reservoir and the nuclear gas content available for black hole fuelling.


title: A Deep Learning Approach to Mapping Irrigation: IrrMapper-U-Net

abstract: Accurate maps of irrigation are essential for understanding and managing water resources. We present a new method of mapping irrigation and demonstrate its accuracy for the state of Montana from years 2000-2019. The method is based off of an ensemble of convolutional neural networks that use reflectance information from Landsat imagery to classify irrigated pixels, that we call IrrMapper-U-Net. The methodology does not rely on extensive feature engineering and does not condition the classification with land use information from existing geospatial datasets. The ensemble does not need exhaustive hyperparameter tuning and the analysis pipeline is lightweight enough to be implemented on a personal computer. Furthermore, the proposed methodology provides an estimate of the uncertainty associated with classification. We evaluated our methodology and the resulting irrigation maps using a highly accurate novel spatially-explicit ground truth data set, using county-scale USDA surveys of irrigation extent, and using cadastral surveys. We found that that our method outperforms other methods of mapping irrigation in Montana in terms of overall accuracy and precision. We found that our method agrees better statewide with the USDA National Agricultural Statistics Survey estimates of irrigated area compared to other methods, and has far fewer errors of commission in rainfed agriculture areas. The method learns to mask clouds and ignore Landsat 7 scan-line failures without supervision, reducing the need for preprocessing data. This methodology has the potential to be applied across the entire United States and for the complete Landsat record.


title: Remember What You Want to Forget: Algorithms for Machine Unlearning

abstract: We study the problem of unlearning datapoints from a learnt model. The learner first receives a dataset $S$ drawn i.i.d. from an unknown distribution, and outputs a model $\widehat{w}$ that performs well on unseen samples from the same distribution. However, at some point in the future, any training datapoint $z \in S$ can request to be unlearned, thus prompting the learner to modify its output model while still ensuring the same accuracy guarantees. We initiate a rigorous study of generalization in machine unlearning, where the goal is to perform well on previously unseen datapoints. Our focus is on both computational and storage complexity. For the setting of convex losses, we provide an unlearning algorithm that can unlearn up to $O(n/d^{1/4})$ samples, where $d$ is the problem dimension. In comparison, in general, differentially private learning (which implies unlearning) only guarantees deletion of $O(n/d^{1/2})$ samples. This demonstrates a novel separation between differential privacy and machine unlearning.


title: Finding Efficient Trade-offs in Multi-Fidelity Response Surface Modeling

abstract: In the context of optimization approaches to engineering applications, time-consuming simulations are often utilized which can be configured to deliver solutions for various levels of accuracy, commonly referred to as different fidelity levels. It is common practice to train hierarchical surrogate models on the objective functions in order to speed-up the optimization process. These operate under the assumption that there is a correlation between the high- and low-fidelity versions of the problem that can be exploited to cheaply gain information. In the practical scenario where the computational budget has to be allocated between multiple fidelities, limited guidelines are available to help make that division. In this paper we evaluate a range of different choices for a two-fidelity setup that provide helpful intuitions about the trade-off between evaluating in high- or low-fidelity. We present a heuristic method based on subsampling from an initial Design of Experiments (DoE) to find a suitable division of the computational budget between the fidelity levels. This enables the setup of multi-fidelity optimizations which utilize the available computational budget efficiently, independent of the multi-fidelity model used.


title: CLAIMED, a visual and scalable component library for Trusted AI

abstract: Deep Learning models are getting more and more popular but constraints on explainability, adversarial robustness and fairness are often major concerns for production deployment. Although the open source ecosystem is abundant on addressing those concerns, fully integrated, end to end systems are lacking in open source. Therefore we provide an entirely open source, reusable component framework, visual editor and execution engine for production grade machine learning on top of Kubernetes, a joint effort between IBM and the University Hospital Basel. It uses Kubeflow Pipelines, the AI Explainability360 toolkit, the AI Fairness360 toolkit and the Adversarial Robustness Toolkit on top of ElyraAI, Kubeflow, Kubernetes and JupyterLab. Using the Elyra pipeline editor, AI pipelines can be developed visually with a set of jupyter notebooks.


title: Three-dimensional flux vacua from IIB on co-calibrated G2 orientifolds

abstract: We derive the 3D N=1 superpotential for the closed string sector of type IIB supergravity on toroidal O5 orientifolds with co-calibrated G2 structure and RR background flux. We find that such compactifications can provide full closed string moduli stabilization on supersymmetric AdS$_3$ vacua, and once we include brane-supersymmetry-breaking we also find indication for the existence of classical 3D de Sitter solutions. The latter however are rather difficult to reconcile with the shape moduli stabilization and flux quantization. We also discuss the possibility of achieving scale separation in AdS$_3$ and dS$_3$ vacua, but such effects seems to be hindered by the geometric flux quantization.


title: Cosmological Constraints from Galaxy Cluster Sparsity, Cluster Gas Mass Fraction and Baryon Acoustic Oscillations Data

abstract: In recent years, the availability of large, complete cluster samples has enabled numerous cosmological parameter inference analyses using cluster number counts. These have provided constraints on the cosmic matter density $\Omega_m$ and the amplitude of matter density fluctuations $\sigma_8$ alternative to those obtained from other standard probes. However, systematics uncertainties, such as the mass calibration bias and selection effects, may still significantly affect these data analyses. Hence, it is timely to explore other proxies of galaxy cluster cosmology that can provide cosmological constraints complementary to those obtained from cluster number counts. Here, we use measurements of the cluster sparsity from weak lensing mass estimates of the LC$^2$-{\it single} and HSC-XXL cluster catalogs to infer constraints on a flat $\Lambda$CDM model. The cluster sparsity has the advantage of being insensitive to selection and mass calibration bias. On the other hand, it primarily constrains a degenerate combination of $\Omega_m$ and $\sigma_8$ (along approximately constant curves of $S_8=\sigma_8\sqrt{\Omega_m/0.3}$), and to less extent the reduced Hubble parameter $h$. Hence, in order to break the internal parameter degeneracies we perform a combined likelihood analysis of cluster sparsities with cluster gas mass fraction measurements and BAO data. We find marginal constraints that are competitive with those from other standard cosmic probes: $\Omega_m=0.316\pm 0.013$, $\sigma_8=0.757\pm 0.067$ (corresponding to $S_8=0.776\pm 0.064$) and $h=0.696\pm 0.017$ at $1\sigma$. Moreover, assuming a conservative Gaussian prior on the mass bias of gas mass fraction data, we find a lower limit on the gas depletion factor $Y_{b,500c}\gtrsim 0.89$.


title: Robust hybridization gap in a Kondo Insulator YbB${}_{12}$ probed by femtosecond optical spectroscopy

abstract: In heavy fermions the relaxation dynamics of photoexcited carriers has been found to be governed by the low energy indirect gap, E$_{g}$, resulting from hybridization between localized moments and conduction band electrons. Here, carrier relaxation dynamics in a prototype Kondo insulator YbB${}_{12}$ is studied over large range of temperatures and over three orders of magnitude. We utilize the intrinsic non-linearity of dynamics to quantitatively determine microscopic parameters, such as electron-hole recombination rate. The extracted value reveals that hybridization is accompanied by a strong charge transfer from localized 4f-levels. The results imply the presence of a hybridization gap up to temperatures of the order of E$_{g}$/k$_{B}\approx200$ K, which is extremely robust against electronic excitation. Finally, below 20 K the data reveal changes in the low energy electronic structure, attributed to short-range antiferromagnetic correlations between the localized levels.


title: A vertex scheme for two-phase flow in heterogeneous media

abstract: This paper presents the numerical solution of immiscible two-phase flows in porous media, obtained by a first-order finite element method equipped with mass-lumping and flux up-winding. The unknowns are the physical phase pressure and phase saturation. Our numerical experiments confirm that the method converges optimally for manufactured solutions. For both structured and unstructured meshes, we observe the high-accuracy wetting saturation profile that ensures minimal numerical diffusion at the front. Performing several examples of quarter-five spot problems in two and three dimensions, we show that the method can easily handle heterogeneities in the permeability field. Two distinct features that make the method appealing to reservoir simulators are: (i) maximum principle is satisfied, and (ii) mass balance is locally conserved.


title: Extremal points of Lorenz curves and applications to inequality analysis

abstract: We find the set of extremal points of Lorenz curves with fixed Gini index and compute the maximal $L^1$-distance between Lorenz curves with given values of their Gini coefficients. As an application we introduce a bidimensional index that simultaneously measures relative inequality and dissimilarity between two populations. This proposal employs the Gini indices of the variables and an $L^1$-distance between their Lorenz curves. The index takes values in a right-angled triangle, two of whose sides characterize perfect relative inequality-expressed by the Lorenz ordering between the underlying distributions. Further, the hypotenuse represents maximal distance between the two distributions. As a consequence, we construct a chart to, graphically, either see the evolution of (relative) inequality and distance between two income distributions over time or to compare the distribution of income of a specific population between a fixed time point and a range of years. We prove the mathematical results behind the above claims and provide a full description of the asymptotic properties of the plug-in estimator of this index. Finally, we apply the proposed bidimensional index to several real EU-SILC income datasets to illustrate its performance in practice.


title: Epistemic Signaling Games for Cyber Deception with Asymmetric Recognition

abstract: This study provides a model of cyber deception with asymmetric recognition represented by private beliefs. Signaling games, which are often used in existing works, are built on the implicit premise that the receiver's belief is public information. However, this assumption, which leads to symmetric recognition, is unrealistic in adversarial decision making. For a precise evaluation of risks arising from cognitive gaps, this paper proposes epistemic signaling games based on the Mertens-Zamir model, which explicitly quantifies players' asymmetric recognition. Equilibria of the games are analytically characterized with an interpretation.


title: The Effect of Network Topology on Credit Network Throughput

abstract: Credit networks rely on decentralized, pairwise trust relationships (channels) to exchange money or goods. Credit networks arise naturally in many financial systems, including the recent construct of payment channel networks in blockchain systems. An important performance metric for these networks is their transaction throughput. However, predicting the throughput of a credit network is nontrivial. Unlike traditional communication channels, credit channels can become imbalanced; they are unable to support more transactions in a given direction once the credit limit has been reached. This potential for imbalance creates a complex dependency between a network's throughput and its topology, path choices, and the credit balances (state) on every channel. Even worse, certain combinations of these factors can lead the credit network to deadlocked states where no transactions can make progress. In this paper, we study the relationship between the throughput of a credit network and its topology and credit state. We show that the presence of deadlocks completely characterizes a network's throughput sensitivity to different credit states. Although we show that identifying deadlocks in an arbitrary topology is NP-hard, we propose a peeling algorithm inspired by decoding algorithms for erasure codes that upper bounds the severity of the deadlock. We use the peeling algorithm as a tool to compare the performance of different topologies as well as to aid in the synthesis of topologies robust to deadlocks.


title: Interstellar Objects in the Solar System: 1. Isotropic Kinematics from the Gaia Early Data Release 3

abstract: 1I/'Oumuamua (or 1I) and 2I/Borisov (or 2I), the first InterStellar Objects (ISOs) discovered passing through the solar system, have opened up entirely new areas of exobody research. Finding additional ISOs and planning missions to intercept or rendezvous with these bodies will greatly benefit from knowledge of their likely orbits and arrival rates. Here, we use the local velocity distribution of stars from the Gaia Early Data Release 3 Catalogue of Nearby Stars and a standard gravitational focusing model to predict the velocity dependent flux of ISOs entering the solar system. With an 1I-type ISO number density of $\sim$0.1 AU$^{-3}$, we predict that a total of $\sim$6.9 such objects per year should pass within 1 AU of the Sun. There will be a fairly large high-velocity tail to this flux, with half of the incoming ISOs predicted to have a velocity at infinity, v$_{\infty}$, $>$ 40 km s$^{-1}$. Our model predicts that $\sim$92\% of incoming ISOs will be residents of the galactic thin disk, $\sim$6\% ($\sim$4 per decade) will be from the thick disk, $\sim$1 per decade will be from the halo and at most $\sim$3 per century will be unbound objects, ejected from our galaxy or entering the Milky Way from another galaxy. The rate of ISOs with very low v$_{\infty}$ $\lesssim$ 1.5 km s$^{-1}$ is so low in our model that any incoming very low velocity ISOs are likely to be previously lost solar system objects. Finally, we estimate a cometary ISO number density of $\sim$7 $\times$ 10$^{-5}$ AU$^{-3}$ for 2I type ISOs, leading to discovery rates for these objects possibly approaching once per decade with future telescopic surveys.


title: Decreasing Impatience

abstract: We characterize decreasing impatience, a common behavioral phenomenon in intertemporal choice. Discount factors that display decreasing impatience are characterized through a convexit y axiom for investments at fixed interest rates. Then we show that they are equivalent to a geometric average of generalized quasi-hype rbolic discount rates. Finally, they emerge through parimutuel preference aggregation of exponential discount factors.


title: Probing inflation with large-scale structure data: the contribution of information at small scales

abstract: Upcoming full-sky large-scale structure surveys such as Euclid can probe the primordial Universe. Using the specifications for the Euclid survey, we estimate the constraints on the inflation potential beyond slow-roll. We use mock Euclid and Planck data from fiducial cosmological models using the Wiggly Whipped Inflation (WWI) framework, which generates features in the primordial power spectrum. We include Euclid cosmic shear and galaxy clustering, with two setups (Conservative and Realistic) for the non-linear cut-off. We find that the addition of Euclid data gives an improvement in constraints in the WWI potential, with the Realistic setup providing marginal improvement over the Conservative for most models. This shows that Euclid may allow us to identify oscillations in the primordial spectrum present at intermediate to small scales.


title: Efficient generative modeling of protein sequences using simple autoregressive models

abstract: Generative models emerge as promising candidates for novel sequence-data driven approaches to protein design, and for the extraction of structural and functional information about proteins deeply hidden in rapidly growing sequence databases. Here we propose simple autoregressive models as highly accurate but computationally efficient generative sequence models. We show that they perform similarly to existing approaches based on Boltzmann machines or deep generative models, but at a substantially lower computational cost (by a factor between $10^2$ and $10^3$). Furthermore, the simple structure of our models has distinctive mathematical advantages, which translate into an improved applicability in sequence generation and evaluation. Within these models, we can easily estimate both the probability of a given sequence, and, using the model's entropy, the size of the functional sequence space related to a specific protein family. In the example of response regulators, we find a huge number of ca. $10^{68}$ possible sequences, which nevertheless constitute only the astronomically small fraction $10^{-80}$ of all amino-acid sequences of the same length. These findings illustrate the potential and the difficulty in exploring sequence space via generative sequence models.

categories: q-bio.BM cond-mat.dis-nn cond-mat.stat-mech q-bio.QM


title: Accelerating Second-Order Differential Dynamic Programming for Rigid-Body Systems

abstract: This letter presents a method to reduce the computational demands of including second-order dynamics sensitivity information into the Differential Dynamic Programming (DDP) trajectory optimization algorithm. An approach to DDP is developed where all the necessary derivatives are computed with the same complexity as in the iterative Linear Quadratic Regulator (iLQR). Compared to linearized models used in iLQR, DDP more accurately represents the dynamics locally, but it is not often used since the second-order derivatives of the dynamics are tensorial and expensive to compute. This work shows how to avoid the need for computing the derivative tensor by instead leveraging reverse-mode accumulation of derivative information to compute a key vector-tensor product directly. We also show how the structure of the dynamics can be used to further accelerate these computations in rigid-body systems. Benchmarks of this approach for trajectory optimization with multi-link manipulators show that the benefits of DDP can often be included without sacrificing evaluation time, and can be done in fewer iterations than iLQR.


title: An Almost Optimal Edit Distance Oracle

abstract: We consider the problem of preprocessing two strings $S$ and $T$, of lengths $m$ and $n$, respectively, in order to be able to efficiently answer the following queries: Given positions $i,j$ in $S$ and positions $a,b$ in $T$, return the optimal alignment of $S[i \mathinner{.\,.} j]$ and $T[a \mathinner{.\,.} b]$. Let $N=mn$. We present an oracle with preprocessing time $N^{1+o(1)}$ and space $N^{1+o(1)}$ that answers queries in $\log^{2+o(1)}N$ time. In other words, we show that we can query the alignment of every two substrings in almost the same time it takes to compute just the alignment of $S$ and $T$. Our oracle uses ideas from our distance oracle for planar graphs [STOC 2019] and exploits the special structure of the alignment graph. Conditioned on popular hardness conjectures, this result is optimal up to subpolynomial factors. Our results apply to both edit distance and longest common subsequence (LCS). The best previously known oracle with construction time and size $\mathcal{O}(N)$ has slow $\Omega(\sqrt{N})$ query time [Sakai, TCS 2019], and the one with size $N^{1+o(1)}$ and query time $\log^{2+o(1)}N$ (using a planar graph distance oracle) has slow $\Omega(N^{3/2})$ construction time [Long & Pettie, SODA 2021]. We improve both approaches by roughly a $\sqrt N$ factor.


title: The discovery of a highly accreting, radio-loud quasar at z=6.82

abstract: Radio sources at the highest redshifts can provide unique information on the first massive galaxies and black holes, the densest primordial environments, and the epoch of reionization. The number of astronomical objects identified at z>6 has increased dramatically over the last few years, but previously only three radio-loud (R2500>10) sources had been reported at z>6, with the most distant being a quasar at z=6.18. Here we present the discovery and characterization of P172+18, a radio-loud quasar at z=6.823. This source has an MgII-based black hole mass of ~3x10^8 Msun and is one of the fastest accreting quasars, consistent with super-Eddington accretion. The ionized region around the quasar is among the largest measured at these redshifts, implying an active phase longer than the average lifetime of the z>6 quasar population. From archival data, there is evidence that its 1.4 GHz emission has decreased by a factor of two over the last two decades. The quasar's radio spectrum between 1.4 and 3.0 GHz is steep (alpha=-1.31) and has a radio-loudness parameter R2500~90. A second steep radio source (alpha=-0.83) of comparable brightness to the quasar is only 23.1" away (~120 kpc at z=6.82; projection probability <2%), but shows no optical or near-infrared counterpart. Further follow-up is required to establish whether these two sources are physically associated.


title: PVG at WASSA 2021: A Multi-Input, Multi-Task, Transformer-Based Architecture for Empathy and Distress Prediction

abstract: Active research pertaining to the affective phenomenon of empathy and distress is invaluable for improving human-machine interaction. Predicting intensities of such complex emotions from textual data is difficult, as these constructs are deeply rooted in the psychological theory. Consequently, for better prediction, it becomes imperative to take into account ancillary factors such as the psychological test scores, demographic features, underlying latent primitive emotions, along with the text's undertone and its psychological complexity. This paper proffers team PVG's solution to the WASSA 2021 Shared Task on Predicting Empathy and Emotion in Reaction to News Stories. Leveraging the textual data, demographic features, psychological test score, and the intrinsic interdependencies of primitive emotions and empathy, we propose a multi-input, multi-task framework for the task of empathy score prediction. Here, the empathy score prediction is considered the primary task, while emotion and empathy classification are considered secondary auxiliary tasks. For the distress score prediction task, the system is further boosted by the addition of lexical features. Our submission ranked 1$^{st}$ based on the average correlation (0.545) as well as the distress correlation (0.574), and 2$^{nd}$ for the empathy Pearson correlation (0.517).


title: Theoretical realization of a two qubit quantum controlled-not logic gate and a single qubit Hadamard logic gate in the anti-Jaynes-Cummings model

abstract: We provide a theoretical scheme for realizing a Hadamard and a quantum controlled-NOT logic gates operations in the anti-Jaynes-Cummings interaction process. Standard Hadamard operation for a specified initial atomic state is achieved by setting a specific sum frequency and photon number in the anti-Jaynes-Cummings qubit state transition operation with the interaction component of the anti-Jaynes-Cummings Hamiltonian generating the state transitions. The quantum controlled-NOT logic gate is realized when a single atomic qubit defined in a two-dimensional Hilbert space is the control qubit and two non-degenerate and orthogonal polarized cavities defined in a two-dimensional Hilbert space make the target qubit. With precise choice of interaction time in the anti-Jaynes-Cummings qubit state transition operations defined in the anti-Jaynes-Cummings sub-space spanned by normalized but non-orthogonal basic qubit state vectors, we obtain ideal unit probabilities of success in the quantum controlled-NOT operations.


title: On a class of generalized Fermat equations of signature $(2,2n,3)$

abstract: We consider the Diophantine equation $7x^{2} + y^{2n} = 4z^{3}$. We determine all solutions to this equation for $n = 2, 3, 4$ and $5$. We formulate a Kraus type criterion for showing that the Diophantine equation $7x^{2} + y^{2p} = 4z^{3}$ has no non-trivial proper integer solutions for specific primes $p > 7$. We computationally verify the criterion for all primes $7 < p < 10^9$, $p \neq 13$. We use the symplectic method and quadratic reciprocity to show that the Diophantine equation $7x^{2} + y^{2p} = 4z^{3}$ has no non-trivial proper solutions for a positive proportion of primes $p$. In the paper \cite{ChDS} we consider the Diophantine equation $x^{2} +7y^{2n} = 4z^{3}$, determining all families of solutions for $n=2$ and $3$, as well as giving a (mostly) conjectural description of the solutions for $n=4$ and primes $n \geq 5$.


title: Total positive curvature and the equality case in the relative isoperimetric inequality outside convex domains

abstract: We settle the case of equality for the relative isoperimetric inequality outside any arbitrary convex set with not empty interior.


title: A nonparametric algorithm for optimal stopping based on robust optimization

abstract: Optimal stopping is a fundamental class of stochastic dynamic optimization problems with numerous applications in finance and operations management. We introduce a new approach for solving computationally-demanding stochastic optimal stopping problems with known probability distributions. The approach uses simulation to construct a robust optimization problem that approximates the stochastic optimal stopping problem to any arbitrary accuracy; we then solve the robust optimization problem to obtain near-optimal Markovian stopping rules for the stochastic optimal stopping problem. In this paper, we focus on designing algorithms for solving the robust optimization problems that approximate the stochastic optimal stopping problems. These robust optimization problems are challenging to solve because they require optimizing over the infinite-dimensional space of all Markovian stopping rules. We overcome this challenge by characterizing the structure of optimal Markovian stopping rules for the robust optimization problems. In particular, we show that optimal Markovian stopping rules for the robust optimization problems have a structure that is surprisingly simple and finite-dimensional. We leverage this structure to develop an exact reformulation of the robust optimization problem as a zero-one bilinear program over totally unimodular constraints. We show that the bilinear program can be solved in polynomial time in special cases, establish computational complexity results for general cases, and develop polynomial-time heuristics by relating the bilinear program to the maximal closure problem from graph theory. Numerical experiments demonstrate that our algorithms for solving the robust optimization problems are practical and can outperform state-of-the-art simulation-based algorithms in the context of widely-studied stochastic optimal stopping problems from high-dimensional option pricing.


title: On topological McAlister semigroups

abstract: In this paper we consider McAlister semigroups over arbitrary cardinals and investigate their algebraic and topological properties. We show that the group of automorphisms of a McAlister semigroup $\mathcal{M}_{\lambda}$ is isomorphic to the direct product $Sym(\lambda){\times}\mathbb{Z}_2$, where $Sym(\lambda)$ is the group of permutations of the cardinal $\lambda$. This fact correlates with the result of Mashevitzky, Schein and Zhitomirski which states that the group of automorphisms of the free inverse semigroup over a cardinal $\lambda$ is isomorphic to the wreath product of $Sym(\lambda)$ and $\mathbb{Z}_2$. Each McAlister semigroup admits a compact semigroup topology. Consequently, the Green's relations $\mathscr D$ and $\mathscr J$ coincide in McAlister semigroups. The latter fact complements results of Lawson. We showed that each non-zero element of a Hausdorff semitopological McAlister semigroup is isolated. This fact is an analogue of the result of Mesyan, Mitchell, Morayne and P\'{e}resse, who proved that each non-zero element of Hausdorff topological polycyclic monoid is isolated. Also, it follows that the free inverse semigroup over a singleton admits only the discrete Hausdorff shift-continuous topology. We proved that a Hausdorff locally compact semitopological semigroup $\mathcal{M}_1$ is either compact or discrete. This fact is similar to the result of Gutik, who showed that a Hausdorff locally compact semitopological polycyclic monoid $\mathcal{P}_1$ is either compact or discrete. However, this dichotomy does not hold for the semigroup $\mathcal{M}_2$. Moreover, $\mathcal{M}_2$ admits continuum many different Hausdorff locally compact inverse semigroup topologies.


title: Ensembles of Random SHAPs

abstract: Ensemble-based modifications of the well-known SHapley Additive exPlanations (SHAP) method for the local explanation of a black-box model are proposed. The modifications aim to simplify SHAP which is computationally expensive when there is a large number of features. The main idea behind the proposed modifications is to approximate SHAP by an ensemble of SHAPs with a smaller number of features. According to the first modification, called ER-SHAP, several features are randomly selected many times from the feature set, and Shapley values for the features are computed by means of "small" SHAPs. The explanation results are averaged to get the final Shapley values. According to the second modification, called ERW-SHAP, several points are generated around the explained instance for diversity purposes, and results of their explanation are combined with weights depending on distances between points and the explained instance. The third modification, called ER-SHAP-RF, uses the random forest for preliminary explanation of instances and determining a feature probability distribution which is applied to selection of features in the ensemble-based procedure of ER-SHAP. Many numerical experiments illustrating the proposed modifications demonstrate their efficiency and properties for local explanation.


title: Earth-skimming Ultra-high Energy Tau Neutrinos simulated with MonteCarlo method and CONEX code

abstract: This paper aims to study the feasibility of building an Earth-skimming cosmic tau neutrinos detector, with the aim of eventually identifying the ideal dimensions of a natural site (mountainvalley) for the detection, with the energy range to be determined (evidently, the highest possible numbers range from 1015 eV to 1020 eV), and possibly locate one such site in Algeria. First, a Monte Carlo simulation of the neutrino-[mountain]matter interaction as well as the resulting decay of the tau lepton is conducted to determine the optimal dimensions of the mountain as well as the location of the tau decay in the valley. Second, a CORSIKA (COsmic Ray Simulation for KAscade) [1] simulation with the CONEX option is conducted to track the evolution of the almost horizontal air shower born from the tau lepton. Among the particles produced in the shower are: electrons, muons, gammas, pions, etc). The study of the spatial distribution of these particles enables the discovery of the optimal width of the valley, and consequently, the distance at which to lay the detection network.


title: A Hybrid Controller for DOS-Resilient String-Stable Vehicle Platoons

abstract: This paper deals with the design of resilient Cooperative Adaptive Cruise Control (CACC) for homogeneous vehicle platoons in which communication is vulnerable to Denial-of-Service (DOS) attacks. We consider DOS attacks as consecutive packet dropouts. We present a controller tuning procedure based on linear matrix inequalities (LMI) that maximizes the resiliency to DOS attacks, while guaranteeing performance and string stability. The design procedure returns controller gains and gives a lower bound on the maximum allowable number of successive packet dropouts. A numerical example is employed to illustrate the effectiveness of the proposed approach.


title: Predicting Kidney Transplant Survival using Multiple Feature Representations for HLAs

abstract: Kidney transplantation can significantly enhance living standards for people suffering from end-stage renal disease. A significant factor that affects graft survival time (the time until the transplant fails and the patient requires another transplant) for kidney transplantation is the compatibility of the Human Leukocyte Antigens (HLAs) between the donor and recipient. In this paper, we propose 4 new biologically-relevant feature representations for incorporating HLA information into machine learning-based survival analysis algorithms. We evaluate our proposed HLA feature representations on a database of over 100,000 transplants and find that they improve prediction accuracy by about 1%, modest at the patient level but potentially significant at a societal level. Accurate prediction of survival times can improve transplant survival outcomes, enabling better allocation of donors to recipients and reducing the number of re-transplants due to graft failure with poorly matched donors.


title: Temperature as perturbation in quantum mechanics

abstract: The perturbative approach was adopted to develop a temperature-dependent version of non-relativistic quantum mechanics in the limit of low-enough temperatures. A generalized, self-consistent Hamiltonian was therefore constructed for an arbitrary quantum-mechanical system in a way that the ground-state Hamiltonian turned out to be just a limiting case at absolute zero. The weak-coupling term connecting the system of interest and its immediate environment was accordingly treated as the perturbation. Applying the obtained generalized Hamiltonian to some typical quantum systems with exact zero-temperature solutions, including the free particle in a box, the free particle in vacuum, and the harmonic oscillator, up to the first order of self-consistency, therefore corrected their associated Hamiltonians, energy spectrums, and wavefunctions to be consistent with the low-temperature limit. Further investigation revealed some kind of quantum tunneling effect by a residual probability for the free particle in a box, as a chief consequence of thermally coupling to the reservoir. The possible effects of thermal environment on the main properties of the wavefunctions were also thoroughly examined and discussed.


title: Conservative Optimistic Policy Optimization via Multiple Importance Sampling

abstract: Reinforcement Learning (RL) has been able to solve hard problems such as playing Atari games or solving the game of Go, with a unified approach. Yet modern deep RL approaches are still not widely used in real-world applications. One reason could be the lack of guarantees on the performance of the intermediate executed policies, compared to an existing (already working) baseline policy. In this paper, we propose an online model-free algorithm that solves conservative exploration in the policy optimization problem. We show that the regret of the proposed approach is bounded by $\tilde{\mathcal{O}}(\sqrt{T})$ for both discrete and continuous parameter spaces.


title: The occurrence of riddled basins and blowout bifurcations in a parametric nonlinear system

abstract: In this paper, a two parameters family $F_{\beta_1,\beta_2}$ of maps of the plane living two different subspaces invariant is studied. We observe that, our model exhibits two chaotic attractors $A_i$, $i=0,1$, lying in these invariant subspaces and identify the parameters at which $A_i$ has a locally riddled basin of attraction or becomes a chaotic saddle. Then, the occurrence of riddled basin in the global sense is investigated in an open region of $\beta_1\beta_2$-plane. We semi-conjugate our system to a random walk model and define a fractal boundary which separates the basins of attraction of the two chaotic attractors, then we describe riddled basin in detail. We show that the model undergos a sequence of bifurcations: "a blowout bifurcation", "a bifurcation to normal repulsion" and "a bifurcation by creating a new chaotic attractor with an intermingled basin". Numerical simulations are presented graphically to confirm the validity of our results.


title: Translating declarative control elements to imperative using 'l-value redefinition graphs'

abstract: We focus on control constructs that allow programmers define actions to be performed when respective conditions are met without requiring the explicit evaluation and testing of conditions as part of an imperative algorithm. Such elements are commonly referred as declarative, not theoretically related to declarative languages. We introduce declarative constructs in the C++ language, presenting the translation method to standard C++. The innovative feature of our method is the accommodation of l-values involving arbitrary pointer / array expressions and objects, supporting immediate runtime evaluation upon content update even if such l-values bind to variant storage locations at runtime. To accomplish this we define 'l-value redefinition graphs', capturing storage binding dependencies among variables, being the floor-plan of our code generation and runtime management approach.


title: On Angular Speed Estimation of Rigid Bodies

abstract: The problem of estimating the angular speed of a solid body from attitude measurements is addressed. To solve this problem, we propose an observer whose dynamics are not constrained to evolve on any specific manifold. This drastically simplifies the analysis of the proposed observer. Using Lyapunov analysis, sufficient conditions for global asymptotic stability of a set wherein the estimation error is equal to zero are established. In addition, the proposed methodology is adapted to deal with angular speed estimation for systems evolving on the unit circle. The approach is illustrated through several numerical simulations.


title: Checkpointing SPAdes for Metagenome Assembly: Transparency versus Performance in Production

abstract: The SPAdes assembler for metagenome assembly is a long-running application commonly used at the NERSC supercomputing site. However, NERSC, like many other sites, has a 48-hour limit on resource allocations. The solution is to chain together multiple resource allocations in a single run, using checkpoint-restart. This case study provides insights into the "pain points" in applying a well-known checkpointing package (DMTCP: Distributed MultiThreaded CheckPointing) to long-running production workloads of SPAdes. This work has exposed several bugs and limitations of DMTCP, which were fixed to support the large memory and fragmented intermediate files of SPAdes. But perhaps more interesting for other applications, this work reveals a tension between the transparency goals of DMTCP and performance concerns due to an I/O bottleneck during the checkpointing process when supporting large memory and many files. Suggestions are made for overcoming this I/O bottleneck, which provides important "lessons learned" for similar applications.


title: The Cauchy problem for the fast $p-$Laplacian evolution equation. Characterization of the global Harnack principle and fine asymptotic behaviour

abstract: We study fine global properties of nonnegative solutions to the Cauchy Problem for the fast $p$-Laplacian evolution equation $u_t=\Delta_p u$ on the whole Euclidean space, in the so-called "good fast diffusion range" $\tfrac{2N}{N+1}<p<2$. It is well-known that non-negative solutions behave for large times as $\mathcal{B}$, the Barenblatt (or fundamental) solution, which has an explicit expression. We prove the so-called Global Harnack Principle (GHP), that is, precise global pointwise upper and lower estimates of nonnegative solutions in terms of $\mathcal{B}$. This can be considered the nonlinear counterpart of the celebrated Gaussian estimates for the linear heat equation. We characterize the maximal (hence optimal) class of initial data such that the GHP holds, by means of an integral tail condition, easy to check. The GHP is then used as a tool to analyze the fine asymptotic behavior for large times. For initial data that satisfy the same integral condition, we prove that the corresponding solutions behave like the Barenblatt with the same mass, uniformly in relative error. When the integral tail condition is not satisfied we show that both the GHP and the uniform convergence in relative error, do not hold anymore, and we provide also explicit counterexamples. We then prove a "generalized GHP", that is, pointwise upper and lower bounds in terms of explicit profiles with a tail different from $\mathcal{B}$. Finally, we derive sharp global quantitative upper bounds of the modulus of the gradient of the solution, and, when data are radially decreasing, we show uniform convergence in relative error for the gradients. To the best of our knowledge, analogous issues for the linear heat equation $p=2$, do not possess such clear answers, only partial results are known.


title: Robust Learning-Based Trajectory Planning for Emerging Mobility Systems

abstract: In this paper, we extend a framework that we developed earlier for coordination of connected and automated vehicles (CAVs) at a signal-free intersection to incorporate uncertainty. Using the possibly noisy observations of actual time trajectories and leveraging Gaussian process regression, we learn the bounded confidence intervals for deviations from the nominal trajectories of CAVs online. Incorporating these confidence intervals, we reformulate the trajectory planning as a robust coordination problem, the solution of which guarantees that constraints in the system are satisfied in the presence of bounded deviations from the nominal trajectories. We demonstrate the effectiveness of our extended framework through a numerical simulation.


title: Consistent Answers of Aggregation Queries using SAT Solvers

abstract: The framework of database repairs and consistent answers to queries is a principled approach to managing inconsistent databases. We describe the first system able to compute the consistent answers of general aggregation queries with the COUNT(A), COUNT(*), SUM(A), MIN(A), and MAX(A) operators, and with or without grouping constructs. Our system uses reductions to optimization versions of Boolean satisfiability (SAT) and then leverages powerful SAT solvers. We carry out an extensive set of experiments on both synthetic and real-world data that demonstrate the usefulness and scalability of this approach.


title: A fault-tolerant domain decomposition method based on space-filling curves

abstract: We propose a simple domain decomposition method for $d$-dimensional elliptic PDEs which involves an overlapping decomposition into local subdomain problems and a global coarse problem. It relies on a space-filling curve to create equally sized subproblems and to determine a certain overlap based on the one-dimensional ordering of the space-filling curve. Furthermore we employ agglomeration and a purely algebraic Galerkin discretization in the construction of the coarse problem. This way, the use of $d$-dimensional geometric information is avoided. The subproblems are dealt with in an additive, parallel way, which gives rise to a subspace correction type linear iteration and a preconditioner for the conjugate gradient method. To make the algorithm fault-tolerant we store on each processor, besides the data of the associated subproblem, a copy of the coarse problem and also the data of a fixed amount of neighboring subproblems with respect to the one-dimensional ordering of the subproblems induced by the space-filling curve. This redundancy then allows to restore the necessary data if processors fail during the computation. Theory supports that the convergence rate of such a linear iteration method stays the same in expectation, and only its order constant deteriorates slightly due to the faults. We observe this in numerical experiments for the preconditioned conjugate gradient method in slightly weaker form as well. Altogether, we obtain a fault-tolerant, parallel and efficient domain decomposition method based on space-filling curves which is especially suited for higher-dimensional elliptic problems.