title: No-Regret Learning in Unknown Games with Correlated Payoffs

abstract: We consider the problem of learning to play a repeated multi-agent game with an unknown reward function. Single player online learning algorithms attain strong regret bounds when provided with full information feedback, which unfortunately is unavailable in many real-world scenarios. Bandit feedback alone, i.e., observing outcomes only for the selected action, yields substantially worse performance. In this paper, we consider a natural model where, besides a noisy measurement of the obtained reward, the player can also observe the opponents' actions. This feedback model, together with a regularity assumption on the reward function, allows us to exploit the correlations among different game outcomes by means of Gaussian processes (GPs). We propose a novel confidence-bound based bandit algorithm GP-MW, which utilizes the GP model for the reward function and runs a multiplicative weight (MW) method. We obtain novel kernel-dependent regret bounds that are comparable to the known bounds in the full information setting, while substantially improving upon the existing bandit results. We experimentally demonstrate the effectiveness of GP-MW in random matrix games, as well as real-world problems of traffic routing and movie recommendation. In our experiments, GP-MW consistently outperforms several baselines, while its performance is often comparable to methods that have access to full information feedback.


title: Specification and Optimal Reactive Synthesis of Run-time Enforcement Shields

abstract: A system with sporadic errors (SSE) is a controller which produces high quality output but it may occasionally violate a critical requirement REQ(I,O). A run-time enforcement shield is a controller which takes (I,O) (coming from SSE) as its input, and it produces a corrected output O' which guarantees the invariance of requirement REQ(I,O'). Moreover, the output sequence O' must deviate from O "as little as possible" to maintain the quality. In this paper, we give a method for logical specification of shields using formulas of logic Quantified Discrete Duration Calculus (QDDC). The specification consists of a correctness requirement REQ as well as a hard deviation constraint HDC which must both be mandatorily and invariantly satisfied by the shield. Moreover, we also use quantitative optimization to give a shield which minimizes the expected value of cumulative deviation in an H-optimal fashion. We show how tool DCSynth implementing soft requirement guided synthesis can be used for automatic synthesis of shields from a given specification. Next, we give logical formulas specifying several notions of shields including the k-Stabilizing shield of Bloem "et al." as well as the Burst-error shield of Wu "et al.", and a new e,d-shield. Shields can be automatically synthesized for all these specifications using the tool DCSynth. We give experimental results showing the performance of our shield synthesis tool in relation to previous work. We also compare the performance of the shields synthesized under diverse hard deviation constraints in terms of their expected deviation and the worst case burst-deviation latency.


title: Extremely Weak Supervised Image-to-Image Translation for Semantic Segmentation

abstract: Recent advances in generative models and adversarial training have led to a flourishing image-to-image (I2I) translation literature. The current I2I translation approaches require training images from the two domains that are either all paired (supervised) or all unpaired (unsupervised). In practice, obtaining paired training data in sufficient quantities is often very costly and cumbersome. Therefore solutions that employ unpaired data, while less accurate, are largely preferred. In this paper, we aim to bridge the gap between supervised and unsupervised I2I translation, with application to semantic image segmentation. We build upon pix2pix and CycleGAN, state-of-the-art seminal I2I translation techniques. We propose a method to select (very few) paired training samples and achieve significant improvements in both supervised and unsupervised I2I translation settings over random selection. Further, we boost the performance by incorporating both (selected) paired and unpaired samples in the training process. Our experiments show that an extremely weak supervised I2I translation solution using only one paired training sample can achieve a quantitative performance much better than the unsupervised CycleGAN model, and comparable to that of the supervised pix2pix model trained on thousands of pairs.


title: On the Order Type of Scattered Context-Free Orderings

abstract: We show that if a context-free grammar generates a language whose lexicographic ordering is well-ordered of type less than $\omega^2$, then its order type is effectively computable.


title: Distance Geometry and Data Science

abstract: Data are often represented as graphs. Many common tasks in data science are based on distances between entities. While some data science methodologies natively take graphs as their input, there are many more that take their input in vectorial form. In this survey we discuss the fundamental problem of mapping graphs to vectors, and its relation with mathematical programming. We discuss applications, solution methods, dimensional reduction techniques and some of their limits. We then present an application of some of these ideas to neural networks, showing that distance geometry techniques can give competitive performance with respect to more traditional graph-to-vector mappings.


title: Parallel Concatenation of Non-Binary Linear Random Fountain Codes with Maximum Distance Separable Codes

abstract: The performance and the decoding complexity of a novel coding scheme based on the concatenation of maximum distance separable (MDS) codes and linear random fountain codes are investigated. Differently from Raptor codes (which are based on a serial concatenation of a high-rate outer block code and an inner Luby-transform code), the proposed coding scheme can be seen as a parallel concatenation of a MDS code and a linear random fountain code, both operating on the same finite field. Upper and lower bounds on the decoding failure probability under maximum-likelihood (ML) decoding are developed. It is shown how, for example, the concatenation of a $(15,10)$ Reed-Solomon (RS) code and a linear random fountain code over a finite field of order $16$, $\mathbb {F}_{16}$, brings to a decoding failure probability $4$ orders of magnitude lower than the one of a linear random fountain code for the same receiver overhead in a channel with a erasure probability of $\epsilon=5\cdot10^{-2}$. It is illustrated how the performance of the novel scheme approaches that of an idealized fountain code for higher-order fields and moderate erasure probabilities. An efficient decoding algorithm is developed for the case of a (generalized) RS code.


title: Flag-transitive block designs and unitary groups

abstract: In this article, we study $2$-designs with $\gcd(r, \lambda)=1$ admitting a flag-transitive automorphism group. The automorphism groups of these designs are point-primitive of almost simple or affine type. We determine all pairs $(\mathcal{D}, G)$, where $\mathcal{D}$ is a $2$-design with $\gcd(r, \lambda)=1$ and $G$ is a flag-transitive almost simple automorphism group of $\mathcal{D}$ whose socle is $X=\mathrm{PSU}(n, q)$ with $(n, q)\neq (3, 2)$ and prove that such a design belongs to one of the two infinite families of Hermitian unitals and Witt-Bose-Shrikhande spaces, or it is isomorphic to a design with parameters $(6, 3, 2)$, $(7, 3, 1)$, $(8, 4, 3)$, $(10, 6, 5)$, $(11, 5, 2)$ or $(28, 7, 2)$.


title: Enhanced energy transfer via graphene-coated wire surface plasmons

abstract: This work analyzes the electromagnetic energy transfer rate between donor and acceptor quantum emitters close to a graphene-coated wire. We discuss the modification of the energy transfer rate when the emitters are interfaced via surface plasmon (SP) environments. All of the notable effects on the spatial dependence of the energy transfer are highlighted and discussed in terms of SP propagation characteristics. Our results show that a dramatically enhancement of the energy transfer occur when the graphene wire SPs are excited. Moreover, different dipole moment orientations influence differently this enhancement. As a consequence of the quasi-one-dimensional graphene wire SPs, we found that the normalized energy transfer rate reaches a maximum value at a donor-acceptor distance which is twice the value corresponding to its two-dimensional counterpart consisting of a single graphene sheet or a flat graphene waveguide. In particular, we provide a simplified model that reproduces the main features of the numerical results.


title: Inclusive productions of $\Upsilon(1S,2S,3S)$ and $\chi_b(1P,2P,3P)$ via the Higgs boson decay

abstract: In this paper, we carry out the complete $\mathcal O(\alpha\alpha_s^{2})$-order study on the inclusive productions of $\Upsilon(nS)$ and $\chi_b(nP)$ ($n=1,2,3$) via the Standard Model Higgs boson decay, within the framework of nonrelativistic QCD. The feeddown effects via the higher excited states are found to be substantial. The color-octet $^3S_1^{[8]}$ state related processes consisting of $H^0 \to b\bar{b}[^3S_1^{[8]}]+g$ and $H^0 \to b\bar{b}[^3S_1^{[8]}]+Q+\bar{Q}$ ($Q=c,b$) play a vital role in the predictions on the decay widths. Moreover, our newly calculated next-to-leading order QCD corrections to $H^0 \to b\bar{b}[^3S_1^{[8]}]+g$ can enhance its leading-order result by 3-4 times, subsequently magnifying the total $^3S_1^{[8]}$ contributions by about $40\%$. Such a remarkable enhancement will to a large extent influence the phenomenological conclusions. For the color-singlet $^3P_J^{[1]}$ state, in addition to $H^{0} \to b\bar{b}[^3P_J^{[1]}]+b+\bar{b}$, the newly introduced light hadrons associated process, $H^{0} \to b\bar{b}[^3P_J^{[1]}]+g+g$, can also provide non-negligible contributions, especially for $^3P_2^{[1]}$. Summing up all the contributions, we have $\mathcal B_{H^0 \to \chi_b(nP)+X} \sim 10^{-6}-10^{-5}$ and $\mathcal B_{H^0 \to \Upsilon(nS)+X} \sim 10^{-5}-10^{-4}$, which meets marginally nowadays LHC experimental data and can help in understanding the heavy quarkonium production mechanism as well as the Yukawa couplings.


title: Knowledge representation and diagnostic inference using Bayesian networks in the medical discourse

abstract: For the diagnostic inference under uncertainty Bayesian networks are investigated. The method is based on an adequate uniform representation of the necessary knowledge. This includes both generic and experience-based specific knowledge, which is stored in a knowledge base. For knowledge processing, a combination of the problem-solving methods of concept-based and case-based reasoning is used. Concept-based reasoning is used for the diagnosis, therapy and medication recommendation and evaluation of generic knowledge. Exceptions in the form of specific patient cases are processed by case-based reasoning. In addition, the use of Bayesian networks allows to deal with uncertainty, fuzziness and incompleteness. Thus, the valid general concepts can be issued according to their probability. To this end, various inference mechanisms are introduced and subsequently evaluated within the context of a developed prototype. Tests are employed to assess the classification of diagnoses by the network.


title: Evolution of 3-dimensional Shape of Passively Evolving and Star-forming Galaxies at $z<1$

abstract: Using the HST/ACS $I_{\rm F814W}$-band data, we investigated distribution of apparent axial ratios of $\sim21000$ galaxies with $M_{V}<-20$ at $0.2<z<1.0$ in the COSMOS field as a function of stellar mass, specific star formation rate (sSFR), and redshift. We statistically estimated intrinsic 3-dimensional shapes of these galaxies by fitting the axial-ratio distribution with triaxial ellipsoid models characterized by face-on (middle-to-long) and edge-on (short-to-long) axial ratios $B/A$ and $C/A$. We found that the transition from thin disk to thick spheroid occurs at $\Delta$MS $\sim-1$ dex, i.e., 10 times lower sSFR than that of the main sequence for galaxies with $M_{\rm star} = 10^{10}$--$10^{11} M_{\odot}$ at $0.2<z<1.0$. Furthermore, the intrinsic thickness ($C/A$) of passively evolving galaxies with $M_{\rm star}=10^{10}$--$10^{11}M_{\odot}$ significantly decreases with time from $C/A \sim 0.40$ -- $0.50$ at $z\sim 0.8$ to $C/A\sim0.33$ -- $0.37$ at $z\sim0.4$, while those galaxies with $M_{\rm star}>10^{11}M_{\odot}$ have $C/A\sim0.5$ irrespective of redshift. On the other hand, star-forming galaxies on the main sequence with $10^{9.5}$--$10^{11}M_{\odot}$ show no significant evolution in their shape at $0.2<z<1.0$, but their thickness depends on stellar mass;more massive star-forming galaxies tend to have lower $C/A$ (thinner shape) than low-mass ones. These results suggest that some fraction of star-forming galaxies with a thin disk, which started to appear around $z\sim1$, quench their star formation without violent morphological change, and these newly added quiescent galaxies with a relatively thin shape cause the significant evolution in the axial-ratio distribution of passively evolving galaxies with $M_{\rm star}<10^{11}M_{\odot}$ at $z<1$.


title: Topological Recursion in The Ramond Sector

abstract: We investigate supereigenvalue models in the Ramond sector and their recursive structure. We prove that the free energy truncates at quadratic order in Grassmann coupling constants, and consider super loop equations of the models with the assumption that the 1/N expansion makes sense. Subject to this assumption, we obtain the associated genus-zero algebraic curve with two ramification points (one regular and the other irregular) and also the supersymmetric partner polynomial equation. Starting with these polynomial equations, we present a recursive formalism that computes all the correlation functions of these models. Somewhat surprisingly, correlation functions obtained from the new recursion formalism have no poles at the irregular ramification point due to a supersymmetric correction -- the new recursion may lead us to a further development of supersymmetric generalizations of the Eynard-Orantin topological recursion.


title: An Automated Engineering Assistant: Learning Parsers for Technical Drawings

abstract: From a set of technical drawings and expert knowledge, we automatically learn a parser to interpret such a drawing. This enables automatic reasoning and learning on top of a large database of technical drawings. In this work, we develop a similarity based search algorithm to help engineers and designers find or complete designs more easily and flexibly. This is part of an ongoing effort to build an automated engineering assistant. The proposed methods make use of both neural methods to learn to interpret images, and symbolic methods to learn to interpret the structure in the technical drawing and incorporate expert knowledge.


title: Monosynaptic inference via finely-timed spikes

abstract: Observations of finely-timed spike relationships in population recordings have been used to support partial reconstruction of neural microcircuit diagrams. In this approach, fine-timescale components of paired spike train interactions are isolated and subsequently attributed to synaptic parameters. Recent perturbation studies strengthen the case for such an inference, yet the complete set of measurements needed to calibrate statistical models are unavailable. To address this gap, we study features of pairwise spiking in a large-scale in vivo dataset where presynaptic neurons were explicitly decoupled from network activity by juxtacellular stimulation. We then construct biophysical models of paired spike trains to reproduce the observed phenomenology of in vivo monosynaptic interactions, including both fine-timescale spike-spike correlations and firing irregularity. A key characteristic of these models is that the paired neurons are coupled by rapidly-fluctuating background inputs. We quantify a monosynapse's causal effect by comparing the postsynaptic train with its counterfactual, when the monosynapse is removed. Subsequently, we develop statistical techniques for estimating this causal effect from the pre- and post-synaptic spike trains. A particular focus is the justification and application of a nonparametric separation of timescale principle to implement synaptic inference. Using simulated data generated from the biophysical models, we characterize the regimes in which the estimators accurately identify the monosynaptic effect. A secondary goal is to initiate a critical exploration of neurostatistical assumptions in terms of biophysical mechanisms, particularly with regards to the challenging but arguably fundamental issue of fast, unobservable nonstationarities in background dynamics.


title: Investigating the high-temperature thermoelectric properties of n-type rutile TiO$_2$

abstract: Transition metal oxides are considered promising thermoelectric materials for harvesting high-temperature waste heat due to their stability, abundance and low toxicity. Despite their typically strong ionic character, they can exhibit surprisingly high power factors $\sigma S^2$, as in n-type SrTiO$_3$ for instance. Thus, it is worth examining other transition metal oxides that might surpass the performances of SrTiO$_3$. This theoretical paper investigates the thermoelectric properties of n-type rutile TiO$_2$, which is the most stable phase of titanium oxide up to 2000 K. The electronic structure is obtained through ab initio calculations, while the prominent features of strong electron-phonon interaction and defects states are modelled using a small number of parameters. The theoretical results are compared with a wealth of experimental data from the literature, yielding very good agreements over a wide range of carrier concentrations. This validates the hypothesis of band conduction in rutile TiO$_2$ and allows the prediction of the high-temperature thermoelectric properties.


title: Field-based mechanical phenotyping of cereal crops to assess lodging resistance

abstract: Plant mechanical failure, also known as lodging, is the cause of significant and unpredictable yield losses in cereal crops. Lodging occurs in two distinct failure modes - stalk lodging and root lodging. Despite the prevalence and detrimental impact of lodging on crop yields, there is little consensus on how to phenotype plants in the field for lodging resistance and thus breed for mechanically resilient plants. This review provides an overview of field-based mechanical testing approaches to assess stalk and root lodging resistance. These approaches are placed in the context of future perspectives. Best practices and recommendations for acquiring field-based mechanical phenotypes of plants are also presented.


title: $f$-distance of knotoids and protein structure

abstract: Recent studies classify the topology of proteins by analysing the distribution of their projections using knotoids. The approximation of this distribution depends on the number of projection directions that are sampled. Here we investigate the relation between knotoids differing only by small perturbations of the direction of projection. Since such knotoids are connected by at most a single forbidden move, we characterise forbidden moves in terms of equivariant band attachment between strongly invertible knots and of strand passages between $\theta$-curves. This allows for the determination of the optimal sample size needed to produce a well approximated knotoid distribution. Based on that and on topological properties of the distribution, we propose a numerical measure for the determination of deeply knotted proteins that does not require the computationally expensive method of subchain analysis.


title: Default Disambiguation for Online Parsers

abstract: Since composed grammars are often ambiguous, grammar composition requires a mechanism for dealing with ambiguity: either ruling it out by using delimiters (which are awkward to work with), or by using disambiguation operators to filter a parse forest down to a single parse tree (where, in general, we cannot be sure that we have covered all possible parse forests). In this paper, we show that default disambiguation, which is inappropriate for batch parsing, works well for online parsing, where it can be overridden by the user if necessary. We extend language boxes -- a delimiter-based algorithm atop incremental parsing -- in such a way that default disambiguation can automatically insert, remove, or resize, language boxes, leading to the automatic language boxes algorithm. The nature of the problem means that default disambiguation cannot always match a user's intention. However, our experimental evaluation shows that automatic language boxes behave acceptably in 98.8% of tests involving compositions of real-world programming languages.


title: A note on alternating direction method of multipliers with generalized augmented terms for constrained sparse least absolute deviation

abstract: This technical note is an ancillary material for our research paper (Itoh and Parente, 2019). We discuss an alternating direction method of multipliers with generalized augmented terms (ADMM-GAT) and introduce a generalized residual balance technique for efficiently employing ADMM-GAT. These techniques are applied to least absolute deviation and its constrained version and their algorithmic details are presented. These algorithms are used for the implementation of the method described in (Itoh and Parente, 2019).


title: On the Right Path: A Modal Logic for Supervised Learning

abstract: Formal learning theory formalizes the process of inferring a general result from examples, as in the case of inferring grammars from sentences when learning a language. Although empirical evidence suggests that children can learn a language without responding to the correction of linguistic mistakes, the importance of Teacher in many other paradigms is significant. Instead of focusing only on learner(s), this work develops a general framework---the supervised learning game (SLG)---to investigate the interaction between Teacher and Learner. In particular, our proposal highlights several interesting features of the agents: on the one hand,Learner may make mistakes in the learning process, and she may also ignore the potential relation between different hypotheses; on the other hand, Teacher is able to correct Learner's mistakes, eliminate potential mistakes and point out the facts ignored by Learner. To reason about strategies in this game, we develop a modal logic of supervised learning (SLL). Broadly, this work takes a small step towards studying the interaction between graph games, logics and formal learning theory.


title: From Bipedal Walking to Quadrupedal Locomotion: Full-Body Dynamics Decomposition for Rapid Gait Generation

abstract: This paper systematically decomposes a quadrupedal robot into bipeds to rapidly generate walking gaits and then recomposes these gaits to obtain quadrupedal locomotion. We begin by decomposing the full-order, nonlinear and hybrid dynamics of a three-dimensional quadrupedal robot, including its continuous and discrete dynamics, into two bipedal systems that are subject to external forces. Using the hybrid zero dynamics (HZD) framework, gaits for these bipedal robots can be rapidly generated (on the order of seconds) along with corresponding controllers. The decomposition is achieved in such a way that the bipedal walking gaits and controllers can be composed to yield dynamic walking gaits for the original quadrupedal robot --- the result is the rapid generation of dynamic quadruped gaits utilizing the full-order dynamics. This methodology is demonstrated through the rapid generation (3.96 seconds on average) of four stepping-in-place gaits and one diagonally symmetric ambling gait at 0.35 m/s on a quadrupedal robot --- the Vision 60, with 36 state variables and 12 control inputs --- both in simulation and through outdoor experiments. This suggested a new approach for fast quadrupedal trajectory planning using full-body dynamics, without the need for empirical model simplification, wherein methods from dynamic bipedal walking can be directly applied to quadrupeds.


title: Integral Quadratic Forms Avoiding Arithmetic Progressions

abstract: For every positive integer k, it is shown that there exists a positive definite diagonal quaternary integral quadratic form that represents all positive integers except for precisely those which lie in k arithmetic progressions. For k=1, all forms with this property are determined.


title: The fastest relaxation rate of higher-dimensional Reissner--Nordstr\"{o}m black hole

abstract: In the eikonal regime, we analytically calculate quasinormal resonance frequencies for massless scalar perturbations of the higher-dimensional Reissner--Nordstr\"{o}m (RN) black holes. Remarkably, we find that the higher-dimensional RN black holes coupled with the massless scalar fields have the fastest relaxation rates in the Schwarzschild limit, this is qualitatively different from the four-dimensional case where the black hole with non-vanishing charge has the fastest relaxation rate.


title: On the restricted conformal group of the (1+n)-Einstein static universe

abstract: Explicit models for the restricted conformal group of the Einstein static universe of dimension greater than two and for its universal covering group are constructed. Based on these models, as an application we determine all oriented and time-oriented conformal Lorentz manifolds whose restricted conformal group has maximal dimension. They amount to the Einstein static universe itself and two countably infinite series of its compact quotients.


title: Informe-pais Brasil

abstract: Structural socioeconomic analysis of Brazil. All basic information about this South American country is gathered in a comprehensive outlook that includes the challenges Brazil faces, as well as their causes and posible economic solutions.


title: Construction of Machine Learned Force Fields with Quantum Chemical Accuracy: Applications and Chemical Insights

abstract: Highly accurate force fields are a mandatory requirement to generate predictive simulations. Here we present the path for the construction of machine learned molecular force fields by discussing the hierarchical pathway from generating the dataset of reference calculations to the construction of the machine learning model, and the validation of the physics generated by the model. We will use the symmetrized gradient-domain machine learning (sGDML) framework due to its ability to reconstruct complex high-dimensional potential-energy surfaces (PES) with high precision even when using just a few hundreds of molecular conformations for training. The data efficiency of the sGDML model allows using reference atomic forces computed with high-level wavefunction-based approaches, such as the $gold$ $standard$ coupled cluster method with single, double, and perturbative triple excitations (CCSD(T)). We demonstrate that the flexible nature of the sGDML framework captures local and non-local electronic interactions (e.g. H-bonding, lone pairs, steric repulsion, changes in hybridization states (e.g. $sp^2 \rightleftharpoons sp^3$), $n\to\pi^*$ interactions, and proton transfer) without imposing any restriction on the nature of interatomic potentials. The analysis of sGDML models trained for different molecular structures at different levels of theory (e.g. density functional theory and CCSD(T)) provides empirical evidence that a higher level of theory generates a smoother PES. Additionally, a careful analysis of molecular dynamics simulations yields new qualitative insights into dynamics and vibrational spectroscopy of small molecules close to spectroscopic accuracy.

categories: physics.chem-ph cond-mat.mtrl-sci physics.atm-clus physics.comp-ph physics.data-an


title: LISA Telescope: Phase noise due to pointing jitter

abstract: In a space based gravitational wave antenna like LISA, involving long light paths linking distant emitter/receiver spacecrafts, signal detection amounts to measuring the light-distance variationsthrough a phase change at the receiver. This is why spurious phase fluctuations due to various mechanical/thermal effects must be carefully studied. We consider here a possible pointing jitter in the light beam sent from the emitter. We show how the resulting phase noise depends on the quality of the wavefront due to the incident beam impinging on the telescope and due to the imperfections of the telescope itself. Namely, we numerically assess the crossed influence of various defects (aberrations and astigmatisms), inherent to a real telescope with pointing fluctuations.


title: New physics in $b\to s\ell\ell$ transitions at one loop

abstract: We investigate new-physics contributions to $b\to s \ell\ell$ transitions in the context of an effective field theory extension of the Standard Model, including operator mixing at one loop. We identify the few scenarios where a single Wilson coefficient, $C/\Lambda^2 \sim 1/{\rm TeV}^2$, induces a substantial shift in the lepton flavour universality ratios $R_K$ and $R_{K^*}$ at one loop, while evading $Z$-pole precision tests, collider bounds, and other flavour constraints. Good fits to the present data are achieved by a left-handed current operator with quark-flavour indices $(2,2)$ or $(3,3)$, hitherto overlooked. Interestingly, the running of the Standard Model Yukawa matrices gives the dominant effect for these scenarios. We match the favoured effective-theory scenarios to minimal, single-mediator models, which are subject to additional stringent constraints. Notably, we recognise three viable instances of a leptoquark with one coupling to fermions only. If the anomalies were confirmed, it appears that one-loop explanations have good prospects of being directly tested at the LHC.


title: From Farey fractions to the Klein quartic and beyond

abstract: In his 1878/79 paper "Ueber die transformation siebenter ordnung der elliptischen functionen", Klein produced his famous 14-sided polygon representing the Klein quartic, his Riemann surface of genus 3 which has PSL(2,7) as its automorphism group. The construction and method of side pairings are fairly complicated. By considering the Farey map modulo 7 we show how to obtain a fundamental polygon for Klein's surface using arithmetic. Now the side pairings are immediate and essentially the same as in Klein's paper. We also extend this idea from 7 to 11 as Klein attempted to do in his follow up paper "Ueber die transformation elfter ordnung der elliptischen functionen", in 1879.


title: Local Probability Conservation in Discrete Time Quantum Walks

abstract: We show that probability is locally conserved in discrete time quantum walks, corresponding to a particle evolving in discrete space and time. In particular, for a spatial structure represented by an arbitrary directed graph, and any unitary evolution of a particle which respects that locality structure, we can define probability currents which also respect the locality structure and which yield the correct final probability distribution.


title: Topological properties of bound-states-in-the-continuum in geometries with broken anisotropy-symmetry

abstract: Waveguiding structures made of anisotropic media support bound states in the continuum (BICs) that arise when the radiation channel of otherwise semi-leaky modes is suppressed. Hitherto, only structures with optical axes aligned in symmetric orientations inside the waveguide plane, where BICs appear as lines in the momentum-frequency dispersion diagram, have been considered. Here we address settings where such symmetry is broken and unveil a number of fundamental new features. Weak and strong symmetry-breaking regimes are identified, corresponding to azimuthal and polar optical axes orientation asymmetries, respectively. The azimuthal symmetry-breaking is found to still preserve the existence loci of BICs in the momentum-frequency dispersion diagram as lines. However, all possible BICs become interferometric, while the polarization separable states that occur in symmetric settings cease to exist. The polar symmetry-breaking has stronger effects and transforms the BICs' existence loci from lines to points, which correspond to full-vector states that exist at discrete values of the optical axis orientation for a given wavelength. Such transformation results in fundamental changes in the topological properties of the radiated field around the BICs.


title: A Standard Model explanation for the excess of electron-like events in MiniBooNE

abstract: We study the dependence of neutral current (NC) neutrino-induced $\pi^0$/photon production ($\nu_\mu + A \to \nu_\mu +1\pi^0 / \gamma + X$) on the atomic number of the target nucleus, A, at 4-momentum transfers relevant to the MiniBooNE experiment: $\Delta$ resonance mass region. Our conclusion is based on experimental data for photon-nucleus interactions from the A2 collaboration at the Mainz MAMI accelerator. We work in the approximation that decays of $\Delta$ resonance unaffected by its production channel, via photon or Z boson. $1\pi^0+X$ production scales as A$^{2/3}$, the surface area of the nucleus. Meanwhile the photons created in $\Delta$ decays will leave the nucleus, and that cross section will be proportional to the atomic number of the nucleus. Thus the ratio of photon production to $\pi^0$ production is proportional to A$^{1/3}$. For carbon $^{12}$C this factor is $\approx$2.3. MiniBooNE normalises the rate of photon production to the measured $\pi^0$ production rate. The reduced neutral pion production rate would yield at least twice as many photons as previously expected, thus significantly lowering the number of unexplained electron-like events.


title: Probing parton saturation with forward $Z^0$-boson production at small transverse momentum in p+p and p+A collisions

abstract: We calculate and compare the differential cross sections for forward $Z^0$-boson production at small transverse momentum, in proton-proton and proton-nucleus collisions, using both the collinear and dilute-dense factorization frameworks. In both cases, we implement a Sudakov resummation of the large logarithms generated by soft-gluon emissions, which is essential in order to describe the transverse momentum distribution of forward $Z^0$ bosons measured at the Tevatron and the LHC. We further compute the nuclear modification factor in the dilute-dense framework, hoping to single out signals of saturation effects at small values of $x$. Our predictions are compared with those obtained in the collinear factorization framework, using two different nuclear parton distribution functions.


title: Statistical and machine learning ensemble modelling to forecast sea surface temperature

abstract: In situ and remotely sensed observations have potential to facilitate data-driven predictive models for oceanography. A suite of machine learning models, including regression, decision tree and deep learning approaches were developed to estimate sea surface temperatures (SST). Training data consisted of satellite-derived SST and atmospheric data from The Weather Company. Models were evaluated in terms of accuracy and computational complexity. Predictive skill were assessed against observations and a state-of-the-art, physics-based model from the European Centre for Medium Weather Forecasting. Results demonstrated that by combining automated feature engineering with machine-learning approaches, accuracy comparable to existing state-of-the-art can be achieved. Models captured seasonal patterns in the data and qualitatively reproduce short-term variations driven by atmospheric forcing. Further, it demonstrated that machine-learning-based approaches can be used as transportable prediction tools for ocean variables -- the data-driven nature of the approach naturally integrates with automatic deployment frameworks, where model deployments are guided by data rather than user-parametrisation and expertise. The low computational cost of inference makes the approach particularly attractive for edge-based computing where predictive models could be deployed on low-power devices in the marine environment.


title: Learning Discrepancy Models From Experimental Data

abstract: First principles modeling of physical systems has led to significant technological advances across all branches of science. For nonlinear systems, however, small modeling errors can lead to significant deviations from the true, measured behavior. Even in mechanical systems, where the equations are assumed to be well-known, there are often model discrepancies corresponding to nonlinear friction, wind resistance, etc. Discovering models for these discrepancies remains an open challenge for many complex systems. In this work, we use the sparse identification of nonlinear dynamics (SINDy) algorithm to discover a model for the discrepancy between a simplified model and measurement data. In particular, we assume that the model mismatch can be sparsely represented in a library of candidate model terms. We demonstrate the efficacy of our approach on several examples including experimental data from a double pendulum on a cart. We further design and implement a feed-forward controller in simulations, showing improvement with a discrepancy model.


title: Brown dwarfs and the minimum mass of stars

abstract: Stars form from large clouds of gas and dust that contract under their own gravity. The birth of a star occurred when a fusion reaction of hydrogen into helium has ignited in its core. The key variable that determines the formation of a star is mass. If the mass of the contracting cloud is below certain minimum value, instead of a star, a substelar object -- known as a brown dwarf -- will form. How much mass is required for a star to form? This article aims to answer this question by means of a simple heuristic argument. The found value is 0.016 solar masses, which is of the same order of magnitude as the accepted value 0.08 solar masses. This article may be useful as pedagogical material in an introductory undergraduate astronomy course.


title: A monolithic model for phase-field fracture and waves in solid-fluid media towards earthquakes

abstract: Coupling of rupture processes in solids with waves also propagating in fluids is a prominent phenomenon arising during tectonic earthquakes. It is executed here in a single `monolithic' model which can asymptotically capture both damageable solids (rocks) and (visco-)elastic fluids (outer core or oceans). Both ruptures on pre-existing lithospheric faults and a birth of new faults in compact rocks are covered by this model, together with emission and propagation of seismic waves, including, e.g., reflection of S-waves and refraction of P-waves on the solid-fluid interfaces. A robust, energy conserving, and convergent staggered FEM discretisation is devised. Using a rather simplified variant of such models for rupture, three computational experiments documenting the applicability of this approach are presented. Some extensions of the model towards more realistic geophysical modelling are outlined, too.


title: Dynamical systems on chain complexes and canonical minimal resolutions

abstract: We introduce notions of vector field and its (discrete time) flow on a chain complex. The resulting dynamical systems theory provides a set of tools with a broad range of applicability that allow, among others, to replace in a canonical way a chain complex with a "smaller" one of the same homotopy type. As applications we construct in an explicit, canonical, and symmetry-preserving fashion a minimal free resolution for every toric ring and every monomial ideal. Our constructions work in all characteristics and over any base field. A key subtle new point is that in certain finitely many positive characteristics (which depend on the object that is being resolved) a transcendental extension of the base field is produced before a resolution is obtained, while in all other characteristics the base field is kept unchanged. In the monomial case we show that such a transcendental base field extension cannot in general be avoided, and we conjecture that the same holds in the toric case.


title: Estimating maternal mortality using data from national civil registration vital statistics systems: A Bayesian hierarchical bivariate random walk model to estimate sensitivity and specificity of reporting

abstract: Civil registration vital statistics (CRVS) data are used to produce national estimates of maternal mortality, but are often subject to substantial reporting errors due to misclassification of maternal deaths. The accuracy of CRVS systems can be assessed by comparing CRVS-based counts of maternal and non-maternal deaths to those obtained from specialized studies, which are rigorous assessments of maternal mortality for a given country-period. We developed a Bayesian bivariate random walk model to estimate sensitivity and specificity of the reporting on maternal mortality in CRVS data, and associated CRVS adjustment factors. The model was fitted to a global data set of CRVS and specialized study data. Validation exercises suggest that the model performs well in terms of predicting CRVS-based proportions of maternal deaths for country-periods without specialized studies. The new model is used by the UN Maternal Mortality Inter-Agency Group to account for misclassification errors when estimating maternal mortality using CRVS data.


title: Semantic and Cognitive Tools to Aid Statistical Science: Replace Confidence and Significance by Compatibility and Surprise

abstract: Researchers often misinterpret and misrepresent statistical outputs. This abuse has led to a large literature on modification or replacement of testing thresholds and $P$-values with confidence intervals, Bayes factors, and other devices. Because the core problems appear cognitive rather than statistical, we review simple aids to statistical interpretations. These aids emphasize logical and information concepts over probability, and thus may be more robust to common misinterpretations than are traditional descriptions. We use the Shannon transform of the $P$-value $p$, also known as the binary surprisal or $S$-value $s=-\log_{2}(p)$, to measure the information supplied by the testing procedure, and to help calibrate intuitions against simple physical experiments like coin tossing. We also use tables or graphs of test statistics for alternative hypotheses, and interval estimates for different percentile levels, to thwart fallacies arising from arbitrary dichotomies. Finally, we reinterpret $P$-values and interval estimates in unconditional terms, which describe compatibility of data with the entire set of analysis assumptions. We illustrate these methods with a reanalysis of data from an existing record-based cohort study. In line with other recent recommendations, we advise that teaching materials and research reports discuss $P$-values as measures of compatibility rather than significance, compute $P$-values for alternative hypotheses whenever they are computed for null hypotheses, and interpret interval estimates as showing values of high compatibility with data, rather than regions of confidence. Our recommendations emphasize cognitive devices for displaying the compatibility of the observed data with various hypotheses of interest, rather than focusing on single hypothesis tests or interval estimates. We believe these simple reforms are well worth the minor effort they require.


title: Re-ID Driven Localization Refinement for Person Search

abstract: Person search aims at localizing and identifying a query person from a gallery of uncropped scene images. Different from person re-identification (re-ID), its performance also depends on the localization accuracy of a pedestrian detector. The state-of-the-art methods train the detector individually, and the detected bounding boxes may be sub-optimal for the following re-ID task. To alleviate this issue, we propose a re-ID driven localization refinement framework for providing the refined detection boxes for person search. Specifically, we develop a differentiable ROI transform layer to effectively transform the bounding boxes from the original images. Thus, the box coordinates can be supervised by the re-ID training other than the original detection task. With this supervision, the detector can generate more reliable bounding boxes, and the downstream re-ID model can produce more discriminative embeddings based on the refined person localizations. Extensive experimental results on the widely used benchmarks demonstrate that our proposed method performs favorably against the state-of-the-art person search methods.


title: A proof of Carleson's $\varepsilon^2$-conjecture

abstract: In this paper we provide a proof of the Carleson $\varepsilon^2$-conjecture. This result yields a characterization (up to exceptional sets of zero length) of the tangent points of a Jordan curve in terms of the finiteness of the associated Carleson $\varepsilon^2$-square function.


title: Code-Switched Language Models Using Neural Based Synthetic Data from Parallel Sentences

abstract: Training code-switched language models is difficult due to lack of data and complexity in the grammatical structure. Linguistic constraint theories have been used for decades to generate artificial code-switching sentences to cope with this issue. However, this require external word alignments or constituency parsers that create erroneous results on distant languages. We propose a sequence-to-sequence model using a copy mechanism to generate code-switching data by leveraging parallel monolingual translations from a limited source of code-switching data. The model learns how to combine words from parallel sentences and identifies when to switch one language to the other. Moreover, it captures code-switching constraints by attending and aligning the words in inputs, without requiring any external knowledge. Based on experimental results, the language model trained with the generated sentences achieves state-of-the-art performance and improves end-to-end automatic speech recognition.


title: To Aid Scientific Inference, Emphasize Unconditional Compatibility Descriptions of Statistics

abstract: All scientific interpretations of statistical outputs depend on background (auxiliary) assumptions that are rarely delineated or explicitly interrogated. These include not only the usual modeling assumptions, but also deeper assumptions about the data-generating mechanism that are implicit in conventional statistical interpretations yet are unrealistic in most health, medical and social research. We provide arguments and methods for reinterpreting statistics such as P-values and interval estimates in unconditional terms, which describe compatibility of observations with an entire set of underlying assumptions, rather than with a narrow target hypothesis conditional on the assumptions. Emphasizing unconditional interpretations helps avoid overconfident and misleading inferences in light of uncertainties about the assumptions used to arrive at the statistical results. These include not only mathematical assumptions, but also those about absence of systematic errors, protocol violations, and data corruption. Unconditional descriptions introduce assumption uncertainty directly into the primary statistical interpretations of results, rather than leaving it for the discussion of limitations after presentation of conditional interpretations. The unconditional approach does not entail different methods or calculations, only different interpretation of the usual results. We view use of unconditional description as a vital component of effective statistical training and presentation. By interpreting statistical outputs in unconditional terms, researchers can avoid making overconfident statements based on statistical outputs. Instead, reports should emphasize the compatibility of results with a range of plausible explanations, including assumption violations.


title: Influence of the Yarkovsky force on Jupiter Trojan asteroids

abstract: Aims. We investigate the influence of the Yarkovsky force on the long-term orbital evolution of Jupiter Trojan asteroids. Methods. Clones of the observed population with different sizes and different thermal properties were numerically integrated for 1 Gyr with and without the Yarkovsky effect. The escape rate of these objects from the Trojan region as well as changes in the libration amplitude, eccentricity, and inclination were used as a metric of the strength of the Yarkovsky effect on the Trojan orbits. Results. Objects with radii $R\leq$1 km are significantly influenced by the Yarkovsky force. The effect causes a depletion of these objects over timescales of a few hundred million years. As a consequence, we expect the size-frequency distribution of small Trojans to show a shallower slope than that of the currently observable population ($R$ $\gtrsim$ 1 km), with a turning point between $R$ = 100 m and $R$ = 1 km. The effect of the Yarkovsky acceleration on the orbits of Trojans depends on the sense of rotation in a complex way. The libration amplitude of prograde rotators decreases with time while the eccentricity increases. Retrograde rotators experience the opposite effect, which results in retrograde rotators being ejected faster from the 1:1 resonance region. Furthermore, for objects affected by the Yarkovsky force, we find indications that the effect tends to smooth out the differences in the orbital distribution between the two clouds.


title: Decoupling stochastic optimal control problems for efficient solution: insights from experiments across a wide range of noise regimes

abstract: We consider the problem of robotic planning under uncertainty in this paper. This problem may be posed as a stochastic optimal control problem, a solution to which is fundamentally intractable owing to the infamous "curse of dimensionality". Hence, we consider the extension of a "decoupling principle" that was recently proposed by some of the authors, wherein a nominal open-loop problem is solved followed by a linear feedback design around the open-loop, and which was shown to be near-optimal to second order in terms of a "small noise" parameter, to a much wider range of noise levels. Our empirical evidence suggests that this allows for tractable planning over a wide range of uncertainty conditions without unduly sacrificing performance.


title: Random concave functions on an equilateral lattice with periodic Hessians I: entropy and Laplacians

abstract: We show that a random concave function having a periodic hessian on an equilateral lattice has a quadratic scaling limit, if the average hessian of the function satisfies certain conditions. We consider the set of all concave functions $g$ on an equilateral lattice $\mathbb L$ that when shifted by an element of $n \mathbb L$, incur addition by a linear function (this condition is equivalent to the periodicity of the hessian of $g$). We identify this set, up to addition by a constant, with a convex polytope $P_n(s)$, where $s$ corresponds to the average hessian. We show that the $\ell_\infty$ diameter of $P_n(s)$ is bounded below by $c(s) n^2$, where $c(s)$ is a positive constant depending only on $s$. Our main result is that, for any $\epsilon_0 > 0$, the normalized Lebesgue measure of all points in $P_n(s)$ that are not contained in a $n^2$ dimensional cube $Q$ of sidelength $2 \epsilon_0 n^2$, centered at the unique (up to addition of a linear term) quadratic polynomial with hessian $s$, tends to $0$ as $n$ tends to $\infty$.


title: Pattern generation and symbolic dynamics in a nanocontact vortex oscillator

abstract: Harnessing chaos or intrinsic nonlinear behaviours from dynamical systems is a promising avenue for the development of unconventional information processing technologies. However, the exploitation of such features in spintronic devices has not been attempted despite the many theoretical and experimental evidence of nonlinear behaviour of the magnetization dynamics in nanomagnetic systems. Here, we propose a first step in that direction by unveiling and characterizing the patterns and symbolic dynamics originating from the nonlinear chaotic time-resolved electrical signals generated experimentally by a nanocontact vortex oscillator (NCVO). We use advanced filtering methods to dissociate nonlinear deterministic patterns from thermal fluctuations and show that the emergence of chaos results in the unpredictable alternation of simple oscillatory patterns controlled by the NCVO's core-polarity switching. With phase-space reconstruction techniques, we perform a symbolic analysis of the time series to assess the level of complexity and entropy generated in the chaotic regime. We find that at the centre of its incommensurate region, it can exhibit maximal entropy and complexity. This suggests that NCVOs are promising nonlinear nanoscale source of entropy that could be harnessed for information processing.


title: The dimension of the boundary of a Liouville quantum gravity metric ball

abstract: Let $\gamma \in (0,2)$, let $h$ be the planar Gaussian free field, and consider the $\gamma$-Liouville quantum gravity (LQG) metric associated with $h$. We show that the essential supremum of the Hausdorff dimension of the boundary of a $\gamma$-LQG metric ball with respect to the Euclidean (resp. $\gamma$-LQG) metric is $2 - \frac{\gamma}{d_\gamma}\left(\frac{2}{\gamma} + \frac{\gamma}{2} \right) + \frac{\gamma^2}{2d_\gamma^2}$ (resp. $d_\gamma-1$), where $d_\gamma$ is the Hausdorff dimension of the whole plane with respect to the $\gamma$-LQG metric. For $\gamma = \sqrt{8/3}$, in which case $d_{\sqrt{8/3}}=4$, we get that the essential supremum of Euclidean (resp. $\sqrt{8/3}$-LQG) dimension of a $\sqrt{8/3}$-LQG ball boundary is $5/4$ (resp. $3$). We also compute the essential suprema of the Euclidean and $\gamma$-LQG Hausdorff dimensions of the intersection of a $\gamma$-LQG ball boundary with the set of metric $\alpha$-thick points of the field $h$ for each $\alpha\in \mathbb R$. Our results show that the set of $\gamma/d_\gamma$-thick points on the ball boundary has full Euclidean dimension and the set of $\gamma$-thick points on the ball boundary has full $\gamma$-LQG dimension.


title: On the Maximal Parameter Range of Global Stability for a Nonlocal Thermostat Model

abstract: The global asymptotic stability of the unique steady state of a nonlinear scalar parabolic equation with a nonlocal boundary condition is studied. The equation describes the evolution of the temperature profile that is subject to a feedback control loop. It can be viewed as a model of a rudimentary thermostat, where a parameter controls the intensity of the heat flow in response to the magnitude of the deviation from the reference temperature at a boundary point. The system is known to undergo a Hopf bifurcation when the parameter exceeds a critical value. Results on the characterization of the maximal parameter range where the reference steady state is globally asymptotically stable are obtained by analyzing a closely related nonlinear Volterra integral equation. Its kernel is derived from the trace of a fundamental solution of a linear heat equation. A version of the Popov criterion is adapted and applied to the Volterra integral equation to obtain a sufficient condition for the asymptotic decay of its solutions.


title: A macroelement stabilization for multiphase poromechanics

abstract: Strong coupling between geomechanical deformation and multiphase fluid flow appears in a variety of geoscience applications. A common discretization strategy for these problems is a continuous Galerkin finite element scheme for the momentum balance equations and a finite volume scheme for the mass balance equations. When applied within a fully-implicit solution strategy, however, this discretization is not intrinsically stable. In the limit of small time steps or low permeabilities, spurious oscillations in the pressure field, i.e. checkerboarding, may be observed. Further, eigenvalues associated with the spurious modes will control the conditioning of the matrices and can dramatically degrade the convergence rate of iterative linear solvers. Here, we propose a stabilization technique in which the balance of mass equations are supplemented with stabilizing flux terms on a macroelement basis. The additional stabilization terms are dependent on a stabilization parameter. We identify an optimal value for this parameter using an analysis of the eigenvalue distribution of the macroelement Schur complement matrix. The resulting method is simple to implement and preserves the underlying sparsity pattern of the original discretization. Another appealing feature of the method is that mass is exactly conserved on macroelements, despite the addition of artificial fluxes. The efficacy of the proposed technique is demonstrated with several numerical examples.


title: Magnetic domain dynamics in an insulating quantum ferromagnet

abstract: The statistics and form of avalanches in a driven system reveal the nature of the underlying energy landscape and dynamics. In conventional metallic ferromagnets, eddy-current back action can dominate the dynamics. Here, we study Barkhausen noise in Li(Ho,Y)F4, an insulating Ising ferromagnet that cannot sustain eddy currents. For large avalanches at temperatures approaching the Curie point, we find a symmetric response free of drag effects. In the low temperature limit, drag effects contribute to the dynamics, which we link to enhanced pinning from local random fields that are enabled by the microscopic dipole-coupled Hamiltonian (the Ising model in transverse field).


title: Efficient Formulation of Ab Initio Quantum Embedding in Periodic Systems: Dynamical Mean-Field Theory

abstract: We present an efficient ab initio dynamical mean-field theory (DMFT) implementation for quantitative simulations in solids. Our DMFT scheme employs ab initio Hamiltonians defined for impurities comprising the full unit cell or a supercell of atoms and for realistic quantum chemical basis sets. We avoid double counting errors by using Hartree-Fock as the low-level theory. Intrinsic and projected atomic orbitals (IAO+PAO) are chosen as the local embedding basis, facilitating numerical bath truncation. Using an efficient integral transformation and coupled-cluster Green's function (CCGF) impurity solvers, we are able to handle embedded impurity problems with several hundred orbitals. We apply our ab initio DMFT approach to study a hexagonal boron nitride monolayer, crystalline silicon, and nickel oxide in the antiferromagnetic phase, with up to 104 and 78 impurity orbitals in spin-restricted and unrestricted cluster DMFT calculations and over 100 bath orbitals. We show that our scheme produces accurate spectral functions compared to both benchmark periodic coupled-cluster computations and experimental spectra.


title: Fine-Tuning Language Models from Human Preferences

abstract: Reward learning enables the application of reinforcement learning (RL) to tasks where reward is defined by human judgment, building a model of reward by asking humans questions. Most work on reward learning has used simulated environments, but complex information about values is often expressed in natural language, and we believe reward learning for language is a key to making RL practical and safe for real-world tasks. In this paper, we build on advances in generative pretraining of language models to apply reward learning to four natural language tasks: continuing text with positive sentiment or physically descriptive language, and summarization tasks on the TL;DR and CNN/Daily Mail datasets. For stylistic continuation we achieve good results with only 5,000 comparisons evaluated by humans. For summarization, models trained with 60,000 comparisons copy whole sentences from the input but skip irrelevant preamble; this leads to reasonable ROUGE scores and very good performance according to our human labelers, but may be exploiting the fact that labelers rely on simple heuristics.


title: NMR, magnetization, and heat capacity studies of the uniform spin-1/2 chain compound Bi6V3O16

abstract: We report the local (NMR) and bulk (magnetization and heat capacity) properties of the vanadium-based S=1/2 uniform spin chain compound Bi6V3O16(Bi4V2O10.66). In the low-temperature {\alpha} phase, the magnetic ions (V4+) are arranged in one-dimensional chains. The magnetic susceptibility shows a broad maximum around 50 K signifying a short-range magnetic order. Heat capacity measurements also reveal low-dimensional magnetism. The 51V magic angle spinning nuclear magnetic resonance measurements clearly show that the magnetic V4+ and nonmagnetic V5+ species are located on different crystallographic sites with no mixed occupation. The spin susceptibility calculated from the shift of the 51V NMR spectra reproduces the behavior observed in magnetic susceptibility and agrees well with the S=1/2 uniform spin chain model with J=113(5) K.


title: Coded Aperture Imaging in High-Energy Astrophysics

abstract: Hard X-ray and low-energy gamma-ray coded-aperture imaging instruments have been highly successful as high-energy surveyors and transient-source discoverers and trackers over the past decades. Albeit having relatively low sensitivity as compared to focussing instruments, coded-aperture telescopes still represent a very good choice for simultaneous, high cadence spectral measurements of individual point sources in large source fields. Here I present a review of the fundamentals of coded-aperture imaging instruments in high-energy astrophysics. Emphasis is on fundamental aspects of the technique, coded-mask instrument characteristics, and properties of the reconstructed images.


title: Efficient Implementation of Ab Initio Quantum Embedding in Periodic Systems: Density Matrix Embedding Theory

abstract: We describe an efficient quantum embedding framework for realistic ab initio density matrix embedding (DMET) calculations in solids. We discuss in detail the choice of orbitals and mapping to a lattice, treatment of the virtual space and bath truncation, and the lattice-to embedded integral transformation. We apply DMET in this ab initio framework to a hexagonal boron nitride monolayer, crystalline silicon, and nickel monoxide in the antiferromagnetic phase, using large embedded clusters with up to 300 embedding orbitals. We demonstrate our formulation of ab initio DMET in the computation of ground-state properties such as the total energy, equation of state, magnetic moment and correlation functions.


title: Dynamical behavior of black-hole spacetimes

abstract: The theory of General Relativity predicts that black holes are a natural constituent of our Universe. Before a black hole binary relaxes to a final state, it undergoes damped oscillations which completely characterize the final object. These oscillations are described by quasinormal modes which are present in many dissipative systems. Quasinormal modes are linked to the way a black hole spacetime responds to small fluctuations. Although the field equations admit a well defined future evolution of spacetime, some black holes, like the ones formed after a binary merger, possess a boundary, beyond which the field equations lose their predictive power. This boundary is called a Cauchy horizon and designates the region beyond which the spacetime can potentially be described in a highly non-unique manner. Strong cosmic censorship states that, generically, suitable initial data should be future inextendible beyond the Cauchy horizon, as a solution to the field equations. Recent studies indicate that the fate of Cauchy horizons, such as those found inside charged and rotating cosmological black holes, is intrinsically connected to the decay of perturbations exterior to the event horizon. As such, the validity of strong cosmic censorship is tied to how effectively the exterior damps fluctuations. In this thesis we will perform a quantitative stability analysis of Cauchy horizons in electrically charged black holes immersed in a de Sitter Universe. We will provide strong numerical evidence which designate that the (linear analog of) strong cosmic censorship is violated by near-extremal Reissner-Nordstr\"om-de Sitter black holes. Moreover, we shall investigate a superradiant instability in $d-$dimensional Reissner-Nordstr\"om-de Sitter black holes and show that the increment of dimensions decreases the timescale of the instability.


title: First-order system least squares finite-elements for singularly perturbed reaction-diffusion equations

abstract: We propose a new first-order-system least squares (FOSLS) finite-element discretization for singularly perturbed reaction-diffusion equations. Solutions to such problems feature layer phenomena, and are ubiquitous in many areas of applied mathematics and modelling. There is a long history of the development of specialized numerical schemes for their accurate numerical approximation. We follow a well-established practice of employing a priori layer-adapted meshes, but with a novel finite-element method that yields a symmetric formulation while also inducing a so-called "balanced" norm. We prove continuity and coercivity of the FOSLS weak form, present a suitable piecewise uniform mesh, and report on the results of numerical experiments that demonstrate the accuracy and robustness of the method.


title: Feature Pyramid Encoding Network for Real-time Semantic Segmentation

abstract: Although current deep learning methods have achieved impressive results for semantic segmentation, they incur high computational costs and have a huge number of parameters. For real-time applications, inference speed and memory usage are two important factors. To address the challenge, we propose a lightweight feature pyramid encoding network (FPENet) to make a good trade-off between accuracy and speed. Specifically, we use a feature pyramid encoding block to encode multi-scale contextual features with depthwise dilated convolutions in all stages of the encoder. A mutual embedding upsample module is introduced in the decoder to aggregate the high-level semantic features and low-level spatial details efficiently. The proposed network outperforms existing real-time methods with fewer parameters and improved inference speed on the Cityscapes and CamVid benchmark datasets. Specifically, FPENet achieves 68.0\% mean IoU on the Cityscapes test set with only 0.4M parameters and 102 FPS speed on an NVIDIA TITAN V GPU.


title: Fixed Point Analysis of Douglas-Rachford Splitting for Ptychography and Phase Retrieval

abstract: Douglas-Rachford Splitting (DRS) methods based on the proximal point algorithms for the Poisson and Gaussian log-likelihood functions are proposed for ptychography and phase retrieval. Fixed point analysis shows that the DRS iterated sequences are always bounded explicitly in terms of the step size and that the fixed points are attracting if and only if the fixed points are regular solutions. This alleviates two major drawbacks of the classical Douglas-Rachford algorithm: slow convergence when the feasibility problem is consistent and divergent behavior when the feasibility problem is inconsistent. Fixed point analysis also leads to a simple, explicit expression for the optimal step size in terms of the spectral gap of an underlying matrix. When applied to the challenging problem of blind ptychography, which seeks to recover both the object and the probe simultaneously, Alternating Minimization with the DRS inner loops, even with a far from optimal step size, converges geometrically under the nearly minimum conditions established in the uniqueness theory.


title: Diversity-enabled sweet spots in layered architectures and speed-accuracy trade-offs in sensorimotor control

abstract: Nervous systems sense, communicate, compute and actuate movement using distributed components with severe trade-offs in speed, accuracy, sparsity, noise and saturation. Nevertheless, brains achieve remarkably fast, accurate, and robust control performance due to a highly effective layered control architecture. Here we introduce a driving task to study how a mountain biker mitigates the immediate disturbance of trail bumps and responds to changes in trail direction. We manipulated the time delays and accuracy of the control input from the wheel as a surrogate for manipulating the characteristics of neurons in the control loop. The observed speed-accuracy trade-offs (SATs) motivated a theoretical framework consisting of layers of control loops with components having diverse speeds and accuracies within each physical level, such as nerve bundles containing axons with a wide range of sizes. Our model explains why the errors from two control loops -- one fast but inaccurate reflexive layer that corrects for bumps, and a planning layer that is slow but accurate -- are additive, and show how the errors in each control loop can be decomposed into the errors caused by the limited speeds and accuracies of the components. These results demonstrate that an appropriate diversity in the properties of neurons across layers helps to create "diversity-enabled sweet spots" (DESSs) so that both fast and accurate control is achieved using slow or inaccurate components.

categories: math.OC cs.IT cs.SY eess.SY math.IT q-bio.NC


title: Deep Model Reference Adaptive Control

abstract: We present a new neuroadaptive architecture: Deep Neural Network based Model Reference Adaptive Control (DMRAC). Our architecture utilizes the power of deep neural network representations for modeling significant nonlinearities while marrying it with the boundedness guarantees that characterize MRAC based controllers. We demonstrate through simulations and analysis that DMRAC can subsume previously studied learning based MRAC methods, such as concurrent learning and GP-MRAC. This makes DMRAC a highly powerful architecture for high-performance control of nonlinear systems with long-term learning properties.


title: Band spectra of periodic hybrid $\delta$-$\delta^\prime$ structures

abstract: We present a detailed study of a generalised one-dimensional Kronig-Penney model using $\delta\text{-}\delta'$ potentials. We analyse the band structure and the density of states in two situations. In the first case we consider an infinite array formed by identical $\delta\text{-}\delta'$ potentials standing at the linear lattice nodes. This case will be known throughout the paper as the one-species hybrid Dirac comb. We investigate the consequences of adding the $\delta'$ interaction to the Dirac comb by comparing the band spectra and the density of states of pure Dirac-$\delta$ combs and one-species hybrid Dirac combs. Secondly we study the quantum system that arises when the periodic potential is the one obtained from the superposition of two one-species hybrid Dirac combs displaced one with respect to the other and with different couplings. The latter will be known as the two-species hybrid Dirac comb. One of the most remarkable results is the appearance of a curvature change in the band spectrum when the $\delta'$ couplings are above a critical value.


title: ModelicaGym: Applying Reinforcement Learning to Modelica Models

abstract: This paper presents ModelicaGym toolbox that was developed to employ Reinforcement Learning (RL) for solving optimization and control tasks in Modelica models. The developed tool allows connecting models using Functional Mock-up Interface (FMI) toOpenAI Gym toolkit in order to exploit Modelica equation-based modelling and co-simulation together with RL algorithms as a functionality of the tools correspondingly. Thus, ModelicaGym facilitates fast and convenient development of RL algorithms and their comparison when solving optimal control problem for Modelicadynamic models. Inheritance structure ofModelicaGymtoolbox's classes and the implemented methods are discussed in details. The toolbox functionality validation is performed on Cart-Pole balancing problem. This includes physical system model description and its integration using the toolbox, experiments on selection and influence of the model parameters (i.e. force magnitude, Cart-pole mass ratio, reward ratio, and simulation time step) on the learning process of Q-learning algorithm supported with the discussion of the simulation results.


title: Graduated Non-Convexity for Robust Spatial Perception: From Non-Minimal Solvers to Global Outlier Rejection

abstract: Semidefinite Programming (SDP) and Sums-of-Squares (SOS) relaxations have led to certifiably optimal non-minimal solvers for several robotics and computer vision problems. However, most non-minimal solvers rely on least-squares formulations, and, as a result, are brittle against outliers. While a standard approach to regain robustness against outliers is to use robust cost functions, the latter typically introduce other non-convexities, preventing the use of existing non-minimal solvers. In this paper, we enable the simultaneous use of non-minimal solvers and robust estimation by providing a general-purpose approach for robust global estimation, which can be applied to any problem where a non-minimal solver is available for the outlier-free case. To this end, we leverage the Black-Rangarajan duality between robust estimation and outlier processes (which has been traditionally applied to early vision problems), and show that graduated non-convexity (GNC) can be used in conjunction with non-minimal solvers to compute robust solutions, without requiring an initial guess. Although GNC's global optimality cannot be guaranteed, we demonstrate the empirical robustness of the resulting robust non-minimal solvers in applications, including point cloud and mesh registration, pose graph optimization, and image-based object pose estimation (also called shape alignment). Our solvers are robust to 70-80% of outliers, outperform RANSAC, are more accurate than specialized local solvers, and faster than specialized global solvers. We also propose the first certifiably optimal non-minimal solver for shape alignment using SOS relaxation.


title: Simultaneous Segmentation and Recognition: Towards more accurate Ego Gesture Recognition

abstract: Ego hand gestures can be used as an interface in AR and VR environments. While the context of an image is important for tasks like scene understanding, object recognition, image caption generation and activity recognition, it plays a minimal role in ego hand gesture recognition. An ego hand gesture used for AR and VR environments conveys the same information regardless of the background. With this idea in mind, we present our work on ego hand gesture recognition that produces embeddings from RBG images with ego hands, which are simultaneously used for ego hand segmentation and ego gesture recognition. To this extent, we achieved better recognition accuracy (96.9%) compared to the state of the art (92.2%) on the biggest ego hand gesture dataset available publicly. We present a gesture recognition deep neural network which recognises ego hand gestures from videos (videos containing a single gesture) by generating and recognising embeddings of ego hands from image sequences of varying length. We introduce the concept of simultaneous segmentation and recognition applied to ego hand gestures, present the network architecture, the training procedure and the results compared to the state of the art on the EgoGesture dataset


title: Towards a Public Key Management Framework for Virtual Assets and Virtual Asset Service Providers

abstract: The recent FATF Recommendations defines virtual assets and virtual assets service providers (VASP), and requires under the Travel Rule that originating VASPs obtain and hold required and accurate originator information and required beneficiary information on virtual asset transfers. In this paper we discuss the notion of key ownership evidence as a core part of originator and beneficiary information required by the FATF Recommendation. We discuss approaches to securely communicate the originator and beneficiary information between VASPs, and review existing standards for public key certificates as applied to VASPs and virtual asset transfers. We propose the notion of a trust network of VASPs in which originator and beneficiary information, including key ownership information, can be exchanged securely while observing individual privacy requirements.


title: Assignment and Pricing of Shared Rides in Ride-Sourcing using Combinatorial Double Auctions

abstract: Transportation Network Companies employ dynamic pricing methods at periods of peak travel to incentivise driver participation and balance supply and demand for rides. Surge pricing multipliers are commonly used and are applied following demand and estimates of customer and driver trip valuations. Combinatorial double auctions have been identified as a suitable alternative, as they can achieve maximum social welfare in the allocation by relying on customers and drivers stating their valuations. A shortcoming of current models, however, is that they fail to account for the effects of trip detours that take place in shared trips and their impact on the accuracy of pricing estimates. To resolve this, we formulate a new shared-ride assignment and pricing algorithm using combinatorial double auctions. We demonstrate that this model is reduced to a maximum weighted independent set model, which is known to be APX-hard. A fast local search heuristic is also presented, which is capable of producing results that lie within 10% of the exact approach for practical implementations. Our proposed algorithm could be used as a fast and reliable assignment and pricing mechanism of ride-sharing requests to vehicles during peak travel times.


title: Multiple Galactic Sources with Emission Above 56 TeV Detected by HAWC

abstract: We present the first catalog of gamma-ray sources emitting above 56 and 100 TeV with data from the High Altitude Water Cherenkov (HAWC) Observatory, a wide field-of-view observatory capable of detecting gamma rays up to a few hundred TeV. Nine sources are observed above 56 TeV, all of which are likely Galactic in origin. Three sources continue emitting past 100 TeV, making this the highest-energy gamma-ray source catalog to date. We report the integral flux of each of these objects. We also report spectra for three highest-energy sources and discuss the possibility that they are PeVatrons.


title: Sample Efficient Policy Gradient Methods with Recursive Variance Reduction

abstract: Improving the sample efficiency in reinforcement learning has been a long-standing research problem. In this work, we aim to reduce the sample complexity of existing policy gradient methods. We propose a novel policy gradient algorithm called SRVR-PG, which only requires $O(1/\epsilon^{3/2})$ episodes to find an $\epsilon$-approximate stationary point of the nonconcave performance function $J(\boldsymbol{\theta})$ (i.e., $\boldsymbol{\theta}$ such that $\|\nabla J(\boldsymbol{\theta})\|_2^2\leq\epsilon$). This sample complexity improves the existing result $O(1/\epsilon^{5/3})$ for stochastic variance reduced policy gradient algorithms by a factor of $O(1/\epsilon^{1/6})$. In addition, we also propose a variant of SRVR-PG with parameter exploration, which explores the initial policy parameter from a prior probability distribution. We conduct numerical experiments on classic control problems in reinforcement learning to validate the performance of our proposed algorithms.


title: Class Feature Pyramids for Video Explanation

abstract: Deep convolutional networks are widely used in video action recognition. 3D convolutions are one prominent approach to deal with the additional time dimension. While 3D convolutions typically lead to higher accuracies, the inner workings of the trained models are more difficult to interpret. We focus on creating human-understandable visual explanations that represent the hierarchical parts of spatio-temporal networks. We introduce Class Feature Pyramids, a method that traverses the entire network structure and incrementally discovers kernels at different network depths that are informative for a specific class. Our method does not depend on the network's architecture or the type of 3D convolutions, supporting grouped and depth-wise convolutions, convolutions in fibers, and convolutions in branches. We demonstrate the method on six state-of-the-art 3D convolution neural networks (CNNs) on three action recognition (Kinetics-400, UCF-101, and HMDB-51) and two egocentric action recognition datasets (EPIC-Kitchens and EGTEA Gaze+).


title: Large e-retailer image dataset for visual search and product classification

abstract: Recent results of deep convolutional networks in visual recognition challenges open the path to a whole new set of disruptive user experiences such as visual search or recommendation. The list of companies offering this type of service is growing everyday but the adoption rate and the relevancy of results may vary a lot. We believe that the availability of large and diverse datasets is a necessary condition to improve the relevancy of such recommendation systems and facilitate their adoption. For that purpose, we wish to share with the community this dataset of more than 12M images of the 7M products of our online store classified into 5K categories. This original dataset is introduced in this article and several features are described. We also present some aspects of the winning solutions of our image classification challenge that was organized on the Kaggle platform around this set of images.


title: Investigating Phase Transition and Morphology of Bi-Te Thermoelectric System

abstract: The optimization of secondary phase in thermoelectric(TE) materials can {helps in improvisation of material's efficiency}. Being a potential {contender} for lower temperature TE application, bismuth telluride(Bi$_2$Te$_3$) nanoparticles were synthesized via different routes and profiles to optimize their pure single phase. Systematic characterizations were performed with the help of X-ray diffraction (XRD), Rietveld refinement and field effect-scanning electron microscopy(FE-SEM) for structural and morphological behavior, while TE properties such as Seebeck coefficient, electrical conductivity and power-factor were measured for the purest sample chosen. Rietveld refinement in the XRD pattern of the samples revealed that only a small amount ($\sim$ 1.6\%) of Bi$_2$Te$_3$ was formed in co-precipitation method, while the hydrothermal technique increases this phase with increment in synthesis duration. This work focused on the phase evolution of Bi$_2$Te$_3$ with increasing synthesis duration time at constant temperature and vice-versa. XRD and Rietveld refinement revealed that the hydrothermal technique (150 $^\circ$C for 48 hours) can synthesize purest samples (84\% Bi$_2$Te$_3$ phase in this case). {FE-SEM and Energy Dispersive X-ray analysis unveiled that the impure phases in the system {were} quantitatively reduced, and it supported by decline in atomic percentage of oxygen from 37\% to 11\%, in addition to this, it was also found that particle size was also decreased with increase in temperature.} The observed electrical conductivity of the chosen sample is $\sim$20 times greater, while Seebeck coefficient is $\sim$3 times lower than that of pure Bi$_2$Te$_3$ phase. The detailed analysis has generalized the growth mechanism in Bi$_2$Te$_3$ phase evolution by the diffusion of Bi into Te nanorods to fabricate hexagonal Bi$_2$Te$_3$.


title: Searching for Magnetostatic Modes in Spin-Polarized Atomic Hydrogen

abstract: We consider a possibility of the magnetostatic type spin waves driven by a long-range magnetic dipole interactions, to account for the peaks in the ESR spectra observed in our previous work. The Walker equation for magnetostatic modes is solved for a cylinder of atomic hydrogen, first in a uniform magnetic field and second in a linearly decreasing magnetic field. The frequency behavior of the solved modes with length of the cylinder and density of the gas is compared to experimental data. We conclude that magnetostatic modes are unlikely to account for the observed modulations of ESR spectra.


title: Simulations of the inner magnetospheric energetic electrons using theIMPTAM-VERB coupled model

abstract: In this study, we present initial results of the coupling between the Inner Magnetospheric Particle Transport and Acceleration Model (IMPTAM) and the Versatile Electron Radiation Belt (VERB-3D) code. IMPTAM traces electrons of 10 - 100 keV energies from the plasma sheet (L = 9 Re) to inner L-shell regions. The flux evolution modeled by IMPTAM is used at the low energy and outer L* computational boundaries of the VERB code (assuming a dipole approximation) to perform radiation belt simulations of energetic electrons. The model was tested on the March 17th, 2013 storm, for a six-day period. Four different simulations were performed and their results compared to satellites observations from Van Allen probes and GOES. The coupled IMPTAM-VERB model reproduces evolution and storm-time features of electron fluxes throughout the studied storm in agreement with the satellite data (within ~0.5 orders of magnitude). Including dynamics of the low energy population at L* = 6.6 increases fluxes closer to the heart of the belt and has a strong impact in the VERB simulations at all energies. However, inclusion of magnetopause losses leads to drastic flux decreases even below L* = 3. The dynamics of low energy electrons (max. 10s of keV) do not affect electron fluxes at energies >= 900 keV. Since the IMPTAM-VERB coupled model is only driven by solar wind parameters and the Dst and Kp indexes, it is suitable as a forecasting tool. In this study, we demonstrate that the estimation of electron dynamics with satellite-data-independent models is possible and very accurate.


title: Gradient Extension of Classical Material Models From Nuclear & Condensed Matter Scales to Earth & Cosmological Scales

abstract: The various mathematical models developed in the past to interpret the behavior of natural and manmade materials were based on observations and experiments made at that time. Classical laws (such as Newton's for gravity, Hooke's for elasticity, Navier-Stokes for fluidity, Fick's/Fourier's for diffusion/heat transfer, Coulomb's for electricity, as well as Maxwell's for electromagnetism and Einstein's for relativity) formed the basis of current technology and shaping of our civilization. The discovery of new phenomena with the aid of recently developed experimental probes have led to various modifications of these laws across disciplines and the scale spectrum: from subatomic and elementary particle physics to cosmology and from atomistic and nano/micro to macro/giga scales. The emergence of nanotechnology and the further advancement of space technology are ultimately connected with the design of novel tools for observation and measurements, as well as the development of new methods and approaches for quantification and understanding. The paper first reviews the author's previously developed weakly nonlocal or gradient models for elasticity, diffusion and plasticity within a unifying internal length gradient (ILG) framework. It then proposes a similar extension for fluids and Maxwell's equations of electromagnetism. Finally, it ventures a gradient modification of Newton's law of gravity and examines its implications to some problems of elementary particle physics, also relevant to cosmology. Along similar lines, it suggests an analogous extension of London's quantum mechanical potential to include both an "attractive" and a "repulsive" branch. It concludes with some comments on a fractional generalization of the ILG framework.


title: Ability of the e-TellTale sensor to detect flow features over wind turbine blades: flow separation/reattachment dynamics

abstract: Monitoring the flow features over wind turbine blades is a challenging task that has become more and more crucial. This paper is devoted to demonstrate the ability of the e-TellTale sensor to detect the flow separation/reattachment dynamics over wind turbine blades. This sensor is made of a strip with a strain gauge sensor at its base. The velocity field was acquired using TR- PIV measurements over an oscillating thick blade section equipped with an e-TellTale sensor. PIV images were post-processed to detect movements of the strip, which was compared to movements of flow. Results show good agreement between the measured velocity field and movements of the strip regarding the separation/reattachment dynamics.


title: Emergent Black Hole Dynamics in Critical Floquet Systems

abstract: While driven interacting quantum matter is generically subject to heating and scrambling, certain classes of systems evade this paradigm. We study such an exceptional class in periodically driven critical (1 + 1)-dimensional systems with a spatially modulated, but disorder-free time evolution operator. Instead of complete scrambling, the excitations of the system remain well-defined. Their propagation is analogous to the evolution along light cones in a curved space-time obtained by two Schwarzschild black holes. The Hawking temperature serves as an order parameter which distinguishes between heating and non-heating phases. Beyond a time scale determined by the inverse Hawking temperature, excitations are absorbed by the black holes resulting in a singular concentration of energy at their center. We obtain these results analytically within conformal field theory, capitalizing on a mapping to sine-square deformed field theories. Furthermore, by means of numerical calculations for an interacting XXZ spin-1/2 chain, we demonstrate that our findings survive lattice regularization.


title: X-ray Observations of a $z\sim6.2$ Quasar/Galaxy Merger

abstract: Quasars at early redshifts ($z > 6$) with companion galaxies offer unique insights into the growth and evolution of the first supermassive black holes. Here, we report on a 150 ks Chandra observation of PSO J308.0416$-$21.2339, a $z=6.23$ quasar with a merging companion galaxy identified in [C II] and rest-frame UV emission. With $72.3^{+9.6}_{-8.6}$ net counts, we find that PSO J308.0416$-$21.2339 is powerful ($L_X = 2.31^{+1.14}_{-0.76} \times 10^{45}\ \textrm{erg}\,\textrm{s}^{-1}\,\textrm{cm}^{-2}$ in rest-frame $2.0-10.0$ keV) yet soft (spectral power-law index $\Gamma=2.39^{+0.37}_{-0.36}$ and optical-to-X-ray slope $\alpha_{\rm OX} = -1.41 \pm 0.11$). In addition, we detect three hard-energy photons $2.''0$ to the west of the main quasar, cospatial with the brightest UV emission of the merging companion. As no soft-energy photons are detected in the same area, this is potentially indicative of a highly-obscured source. With conservative assumptions, and accounting for both background fluctuations and the extended wings of the quasar's emission, these photons only have a probability $P=0.021$ of happening by chance. If confirmed by deeper observations, this system is the first high redshift quasar and companion individually detected in X-rays and is likely a dual AGN.


title: Horizon Molecules in Causal Set Theory

abstract: We propose a new definition of "horizon molecules" in Causal Set Theory following pioneering work by Dou and Sorkin. The new concept applies for any causal horizon and its intersection with any spacelike hypersurface. In the continuum limit, as the discreteness scale tends to zero, the leading behaviour of the expected number of horizon molecules is shown to be the area of the horizon in discreteness units, up to a dimension dependent factor of order one. We also determine the first order corrections to the continuum value, and show how such corrections can be exploited to obtain further geometrical information about the horizon and the spacelike hypersurface from the causal set.


title: Induced Spins from Scattering Experiments of Initially Nonspinning Black Holes

abstract: When two relativistically boosted, nonspinning black holes pass by one another on a scattering trajectory, we might expect the tidal interaction to spin up each black hole. We present the first exploration of this effect, appearing at fourth post-Newtonian order, with full numerical relativity calculations. The basic setup for the calculations involves two free parameters: the initial boost of each black hole and the initial angle between the velocity vectors and a line connecting the centers of the black holes, with zero angle corresponding to a head-on trajectory. To minimize gauge effects, we measure final spins only if the black holes reach a final separation of at least $20M$. Fixing the initial boost, we find that as the initial angle decreases toward the scattering/nonscattering limit, the spin-up grows nonlinearly. In addition, as initial boosts are increased from $0.42c$ to $0.78c$, the largest observed final dimensionless spin on each black hole increases nonlinearly from $0.02$ to $0.20$. Based on these results, we conclude that much higher spin-ups may be possible with larger boosts, although achieving this will require improved numerical techniques.


title: Detailed study of the Milky Way globular cluster Laevens 3

abstract: We present a photometric and spectroscopic study of the Milky Way satellite Laevens 3. Using MegaCam/CFHT g and i photometry and Keck II/DEIMOS multi-object spectroscopy, we refine the structural and stellar properties of the system. The Laevens 3 colour-magnitude diagram shows that it is quite metal-poor, old (13.0 +-1.0 Gyr), and at a distance of 61.4+-1.0 kpc, partly based on two RR Lyrae stars. The system is faint (Mv = -2.8+0.2-0.3 mag) and compact (r = 11.4+-1.0 pc). From the spectroscopy, we constrain the systemic metallicity ([Fe/H]spectro = -1.8+-0.1 dex) but the metallicity and velocity dispersions are both unresolved. Using Gaia DR2, we infer a mean proper motion of (mu_alpha, mu_delta) = (0.51+-0.28, -0.83+-0.27) mas yr-1, which, combined with the system's radial velocity (<vr> = -70.2+-0.5 km s-1), translates into a halo orbit with a pericenter and apocenter of 40.7+5.6-14.7 and 85.6+17.2-5.9 kpc, respectively. Overall, Laevens 3 shares the typical properties of the Milky Way's outer halo globular clusters. Furthermore, we find that this system shows signs of mass-segregation which strengthens our conclusion that Laevens 3 is a globular cluster.


title: A Search for Neutrino Point-Source Populations in 7 Years of IceCube Data with Neutrino-count Statistics

abstract: The presence of a population of point sources in a dataset modifies the underlying neutrino-count statistics from the Poisson distribution. This deviation can be exactly quantified using the non-Poissonian template fitting technique, and in this work we present the first application this approach to the IceCube high-energy neutrino dataset. Using this method, we search in 7 years of IceCube data for point-source populations correlated with the disk of the Milky Way, the Fermi bubbles, the Schlegel, Finkbeiner, and Davis dust map, or with the isotropic extragalactic sky. No evidence for such a population is found in the data using this technique, and in the absence of a signal we establish constraints on population models with source count distribution functions that can be described by a power-law with a single break. The derived limits can be interpreted in the context of many possible source classes. In order to enhance the flexibility of the results, we publish the full posterior from our analysis, which can be used to establish limits on specific population models that would contribute to the observed IceCube neutrino flux.


title: ALMACAL VI: Molecular gas mass density across cosmic time via a blind search for intervening molecular absorbers

abstract: We are just starting to understand the physical processes driving the dramatic change in cosmic star-formation rate between $z\sim 2$ and the present day. A quantity directly linked to star formation is the molecular gas density, which should be measured through independent methods to explore variations due to cosmic variance and systematic uncertainties. We use intervening CO absorption lines in the spectra of mm-bright background sources to provide a census of the molecular gas mass density of the Universe. The data used in this work are taken from ALMACAL, a wide and deep survey utilizing the ALMA calibrator archive. While we report multiple Galactic absorption lines and one intrinsic absorber, no extragalactic intervening molecular absorbers are detected. However, thanks to the large redshift path surveyed ($\Delta z=182$), we provide constraints on the molecular column density distribution function beyond $z\sim 0$. In addition, we probe column densities of N(H$_2$) > 10$^{16}$ atoms~cm$^{-2}$, five orders of magnitude lower than in previous studies. We use the cosmological hydrodynamical simulation IllustrisTNG to show that our upper limits of $\rho ({\rm H}_2)\lesssim 10^{8.3} \text{M}_{\odot} \text{Mpc}^{-3}$ at $0 < z \leq 1.7$ already provide new constraints on current theoretical predictions of the cold molecular phase of the gas. These results are in agreement with recent CO emission-line surveys and are complementary to those studies. The combined constraints indicate that the present decrease of the cosmic star-formation rate history is consistent with an increasing depletion of molecular gas in galaxies compared to $z\sim 2$.


title: The $F$-term Problem and other Challenges of Stringy Quintessence

abstract: We attempt a systematic analysis of string-theoretic quintessence models as an alternative to metastable de Sitter vacua. It appears that, within the boundaries of what is known, large-volume type-IIB flux compactifications are preferred. Here the quintessence scalar is the ratio of certain 4-cycle volumes. It has already been noticed that the volume modulus, which must be stabilized, tends to remain too light. One may call this the "light volume problem". In addition, we identify an "$F$-term problem": The positive energy density of standard-model SUSY breaking is higher than the depth of all known negative contributions. We discuss what it would take to resolve these issues and comment on partially related challenges for axionic quintessence. In particular, large cancellations between positive and negative potential terms appear unavoidable in general. As a further challenge, one should then explain why a small de-tuning cannot be used to uplift into a deep slow-roll regime, violating de Sitter swampland conjectures.