title: Computability models over categories

abstract: Generalising slightly the notions of a strict computability model and of a simulation between them, which were elaborated by Longley and Normann, we define canonical computability models over categories and appropriate Set-valued functors on them. We study the canonical total computability model over a category, and the partial one over a category with pullbacks. Our notions and results are generalised to categories with a base of computability, connecting Rosolini's theory of dominions with the theory of computability models.


title: Deep learned SVT: Unrolling singular value thresholding to obtain better MSE

abstract: Affine rank minimization problem is the generalized version of low rank matrix completion problem where linear combinations of the entries of a low rank matrix are observed and the matrix is estimated from these measurements. We propose a trainable deep neural network by unrolling a popular iterative algorithm called the singular value thresholding (SVT) algorithm to perform this generalized matrix completion which we call Learned SVT (LSVT). We show that our proposed LSVT with fixed layers (say T) reconstructs the matrix with lesser mean squared error (MSE) compared with that incurred by SVT with fixed (same T) number of iterations and our method is much more robust to the parameters which need to be carefully chosen in SVT algorithm.


title: Quantum algorithm for doubling the amplitude of the search problem's solution states

abstract: In this paper we present a quantum algorithm which increases the amplitude of the states corresponding to the solutions of the search problem by a factor of almost two.


title: The Antiferromagnetic Character of the Quantum Phase Transition in the Hubbard Model on the Honeycomb Lattice

abstract: We provide a unified, comprehensive treatment of all operators that contribute to the anti-ferromagnetic, ferromagnetic, and charge-density-wave structure factors and order parameters of the hexagonal Hubbard Model. We use the Hybrid Monte Carlo algorithm to perform a systematic, carefully controlled analysis in the temporal Trotter error and of the thermodynamic limit. We expect our findings to improve the consistency of Monte Carlo determinations of critical exponents. We perform a data collapse analysis and determine the critical exponent $\beta=0.898(37)$ for the semimetal-Mott insulator transition in the hexagonal Hubbard Model. Our methods are applicable to a wide range of lattice theories of strongly correlated electrons.


title: Large enhancement of Edelstein effect in Weyl semimetals from Fermi-arc surface states

abstract: One remarkable feature of Weyl semimetals is the manifestation of their topological nature in the form of the Fermi-arc surface states. In a recent calculation by \cite{Johansson2018}, the current-induced spin polarization or Edelstein effect has been predicted, within the semiclassical Boltzmann theory, to be strongly amplified in a Weyl semimetal TaAs due to the existence of the Fermi arcs. Motivated by this result, we calculate the Edelstein response of an effective model for an inversion-symmetry-breaking Weyl semimetal in the presence of an interface using linear response theory. The scatterings from scalar impurities are included and the vertex corrections are computed within the self-consistent ladder approximation. At chemical potentials close to the Weyl points, we find the surface states have a much stronger response near the interface than the bulk states by about one to two orders of magnitude. At higher chemical potentials, the surface states' response near the interface decreases to be about the same order of magnitude as the bulk states' response. We attribute this phenomenon to the decoupling between the Fermi arc states and bulk states at energies close to the Weyl points. The surface states which are effectively dispersing like a one-dimensional chiral fermion become nearly nondissipative. This leads to a large surface vertex correction and, hence, a strong enhancement of the surface states' Edelstein response.


title: Eliminating Thurston obstructions and controlling dynamics on curves

abstract: Every Thurston map $f\colon S^2\rightarrow S^2$ on a $2$-sphere $S^2$ induces a pull-back operation on Jordan curves $\alpha\subset S^2\setminus P_f$, where $P_f$ is the postcritical set of $f$. Here the isotopy class $[f^{-1}(\alpha)]$ (relative to $P_f$) only depends on the isotopy class $[\alpha]$. We study this operation for Thurston maps with four postcritical points. In this case a Thurston obstruction for the map $f$ can be seen as a fixed point of the pull-back operation. We show that if a Thurston map $f$ with a hyperbolic orbifold and four postcritical points has a Thurston obstruction, then one can "blow up" suitable arcs in the underlying $2$-sphere and construct a new Thurston map $\widehat f$ for which this obstruction is eliminated. We prove that no other obstruction arises and so $\widehat f$ is realized by a rational map. In particular, this allows for the combinatorial construction of a large class of rational Thurston maps with four postcritical points. We also study the dynamics of the pull-back operation under iteration. We exhibit a subclass of our rational Thurston maps with four postcritical points for which we can give positive answer to the global curve attractor problem.


title: Massless Preheating and Electroweak Vacuum Metastability

abstract: Current measurements of Standard Model parameters suggest that the electroweak vacuum is metastable. This metastability has important cosmological implications, because large fluctuations in the Higgs field could trigger vacuum decay in the early universe. For the false vacuum to survive, interactions which stabilize the Higgs during inflation -- e.g., inflaton-Higgs interactions or non-minimal couplings to gravity -- are typically necessary. However, the post-inflationary preheating dynamics of these same interactions could also trigger vacuum decay, thereby recreating the problem we sought to avoid. This dynamics is often assumed catastrophic for models exhibiting scale invariance since these generically allow for unimpeded growth of fluctuations. In this paper, we examine the dynamics of such "massless preheating" scenarios and show that the competing threats to metastability can nonetheless be balanced to ensure viability. We find that fully accounting for both the backreaction from particle production and the effects of perturbative decays reveals a large number of disjoint "islands of (meta)stability" over the parameter space of couplings. Ultimately, the interplay among Higgs-stabilizing interactions plays a significant role, leading to a sequence of dynamical phases that effectively extend the metastable regions to large Higgs-curvature couplings.


abstract: We calculate the NNLO QCD corrections to diphoton production with an additional jet at the LHC. Our calculation represents the first NNLO-accurate prediction for the transverse momentum distribution of the diphoton system. The improvement in the accuracy of the theoretical prediction is significant, by a factor of up to four relative to NLO QCD. Our calculation is exact except for the finite remainder of the two-loop amplitude which is included at leading color. The numerical impact of this approximated contribution is small. The results of this work are expected to further our understanding of the Higgs boson sector and of the behavior of higher-order corrections to LHC processes.


title: Development, validation and clinical usefulness of a prognostic model for relapse in relapsing-remitting multiple sclerosis

abstract: Prognosis on the occurrence of relapses in individuals with Relapsing-Remitting Multiple Sclerosis (RRMS), the most common subtype of Multiple Sclerosis (MS), could support individualized decisions and disease management and could be helpful for efficiently selecting patients in future randomized clinical trials. There are only three previously published prognostic models on this, all of them with important methodological shortcomings. We aim to present the development, internal validation, and evaluation of the potential clinical benefit of a prognostic model for relapses for individuals with RRMS using real world data. We followed seven steps to develop and validate the prognostic model. Finally, we evaluated the potential clinical benefit of the developed prognostic model using decision curve analysis. We selected eight baseline prognostic factors: age, sex, prior MS treatment, months since last relapse, disease duration, number of prior relapses, expanded disability status scale (EDSS), and gadolinium enhanced lesions. We also developed a web application where the personalized probabilities to relapse within two years are calculated automatically. The optimism-corrected c-statistic is 0.65 and the optimism-corrected calibration slope was 0.92. The model appears to be clinically useful between the range 15% and 30% of the threshold probability to relapse. The prognostic model we developed offers several advantages in comparison to previously published prognostic models on RRMS. Importantly, we assessed the potential clinical benefit to better quantify the clinical impact of the model. Our web application, once externally validated in the future, could be used by patients and doctors to calculate the individualized probability to relapse within two years and to inform the management of their disease.


title: VICEROY: GDPR-/CCPA-compliant Enforcement of Verifiable Accountless Consumer Requests

abstract: Recent data protection regulations (such as GDPR and CCPA) grant consumers various rights, including the right to access, modify or delete any personal information collected about them (and retained) by a service provider. To exercise these rights, one must submit a verifiable consumer request proving that the collected data indeed pertains to them. This action is straightforward for consumers with active accounts with a service provider at the time of data collection, since they can use standard (e.g., password-based) means of authentication to validate their requests. However, a major conundrum arises from the need to support consumers without accounts to exercise their rights. To this end, some service providers began requiring such accountless consumers to reveal and prove their identities (e.g., using government-issued documents, utility bills, or credit card numbers) as part of issuing a verifiable consumer request. While understandable as a short-term cure, this approach is cumbersome and expensive for service providers as well as privacy-invasive for consumers. Consequently, there is a strong need to provide better means of authenticating requests from accountless consumers. To achieve this, we propose VICEROY, a privacy-preserving and scalable framework for producing proofs of data ownership, which form a basis for verifiable consumer requests. Building upon existing web techniques and features, VICEROY allows accountless consumers to interact with service providers, and later prove that they are the same person in a privacy-preserving manner, while requiring minimal changes for both parties. We design and implement VICEROY with emphasis on security/privacy, deployability and usability. We also thoroughly assess its practicality via extensive experiments.


title: Efficient Spiking Neural Networks with Radix Encoding

abstract: Spiking neural networks (SNNs) have advantages in latency and energy efficiency over traditional artificial neural networks (ANNs) due to its event-driven computation mechanism and replacement of energy-consuming weight multiplications with additions. However, in order to reach accuracy of its ANN counterpart, it usually requires long spike trains to ensure the accuracy. Traditionally, a spike train needs around one thousand time steps to approach similar accuracy as its ANN counterpart. This offsets the computation efficiency brought by SNNs because longer spike trains mean a larger number of operations and longer latency. In this paper, we propose a radix encoded SNN with ultra-short spike trains. In the new model, the spike train takes less than ten time steps. Experiments show that our method demonstrates 25X speedup and 1.1% increment on accuracy, compared with the state-of-the-art work on VGG-16 network architecture and CIFAR-10 dataset.


title: The Greedy Algorithm is \emph{not} Optimal for On-Line Edge Coloring

abstract: Nearly three decades ago, Bar-Noy, Motwani and Naor showed that no online edge-coloring algorithm can edge color a graph optimally. Indeed, their work, titled "the greedy algorithm is optimal for on-line edge coloring", shows that the competitive ratio of $2$ of the na\"ive greedy algorithm is best possible online. However, their lower bound required bounded-degree graphs, of maximum degree $\Delta = O(\log n)$, which prompted them to conjecture that better bounds are possible for higher-degree graphs. While progress has been made towards resolving this conjecture for restricted inputs and arrivals or for random arrival orders, an answer for fully general \emph{adversarial} arrivals remained elusive. We resolve this thirty-year-old conjecture in the affirmative, presenting a $(1.9+o(1))$-competitive online edge coloring algorithm for general graphs of degree $\Delta = \omega(\log n)$ under vertex arrivals. At the core of our results, and of possible independent interest, is a new online algorithm which rounds a fractional bipartite matching $x$ online under vertex arrivals, guaranteeing that each edge $e$ is matched with probability $(1/2+c)\cdot x_e$, for a constant $c>0.027$.


title: The equivariant Hilbert series of the canonical ring of Fermat curves

abstract: We consider a Fermat curve $F_n:x^n+y^n+z^n=1$ over an algebraically closed field $k$ of characteristic $p\geq0$ and study the action of the automorphism group $G=\left(\mathbb{Z}/n\mathbb{Z}\times\mathbb{Z}/n\mathbb{Z}\right)\rtimes S_3$ on the canonical ring $R=\bigoplus H^0(F_n,\Omega_{F_n}^{\otimes m})$ when $p>3$, $p\nmid n$ and $n-1$ is not a power of $p$. In particular, we explicitly determine the classes $[H^0(F_n,\Omega_{F_n}^{\otimes m})]$ in the Grothendieck group $K_0(G,k)$ of finitely generated $k[G]$-modules, describe the respective equivariant Hilbert series $H_{R,G}(t)$ as a rational function, and use our results to write a program in Sage that computes $H_{R,G}(t)$ for an arbitrary Fermat curve.


title: Fine Selmer groups of congruent $p$-adic Galois representations

abstract: We compare the Pontryagin duals of fine Selmer groups of two congruent $p$-adic Galois representations over admissible pro-$p$, $p$-adic Lie extensions $K_\infty$ of number fields $K$. We prove that in several natural settings the $\pi$-primary submodules of the Pontryagin duals are pseudo-isomorphic over the Iwasawa algebra; if the coranks of the fine Selmer groups are not equal, then we can still prove inequalities between the $\mu$-invariants. In the special case of a $\mathbb{Z}_p$-extension $K_\infty/K$, we also compare the Iwasawa $\lambda$-invariants of the fine Selmer groups, even in situations where the $\mu$-invariants are non-zero. Finally, we prove similar results for certain abelian non-$p$-extensions.


title: Thank you BART! Rewarding Pre-Trained Models Improves Formality Style Transfer

abstract: Scarcity of parallel data causes formality style transfer models to have scarce success in preserving content. We show that fine-tuning pre-trained language (GPT-2) and sequence-to-sequence (BART) models boosts content preservation, and that this is possible even with limited amounts of parallel data. Augmenting these models with rewards that target style and content -- the two core aspects of the task -- we achieve a new state-of-the-art.


title: People construct simplified mental representations to plan

abstract: One of the most striking features of human cognition is the capacity to plan. Two aspects of human planning stand out: its efficiency and flexibility. Efficiency is especially impressive because plans must often be made in complex environments, and yet people successfully plan solutions to myriad everyday problems despite having limited cognitive resources. Standard accounts in psychology, economics, and artificial intelligence have suggested human planning succeeds because people have a complete representation of a task and then use heuristics to plan future actions in that representation. However, this approach generally assumes that task representations are fixed. Here, we propose that task representations can be controlled and that such control provides opportunities to quickly simplify problems and more easily reason about them. We propose a computational account of this simplification process and, in a series of pre-registered behavioral experiments, show that it is subject to online cognitive control and that people optimally balance the complexity of a task representation and its utility for planning and acting. These results demonstrate how strategically perceiving and conceiving problems facilitates the effective use of limited cognitive resources.


title: A Monte Carlo method for solving the electromagnetic scattering problem in dielectric bodies

abstract: In this work, we develop a novel Monte Carlo method for solving the electromagnetic scattering problem. The method is based on a formal solution of the scattering problem as a modified Born series whose coefficients are found by a conformal transformation. The terms of the Born series are approximated by sampling random elements of its matrix representation, computed by the Method of Moments. Unlike other techniques as the Fast Multiple Method, this Monte Carlo method does not require communications between processors, which makes it suitable for large parallel executions.


title: Plot and Rework: Modeling Storylines for Visual Storytelling

abstract: Writing a coherent and engaging story is not easy. Creative writers use their knowledge and worldview to put disjointed elements together to form a coherent storyline, and work and rework iteratively toward perfection. Automated visual storytelling (VIST) models, however, make poor use of external knowledge and iterative generation when attempting to create stories. This paper introduces PR-VIST, a framework that represents the input image sequence as a story graph in which it finds the best path to form a storyline. PR-VIST then takes this path and learns to generate the final story via an iterative training process. This framework produces stories that are superior in terms of diversity, coherence, and humanness, per both automatic and human evaluations. An ablation study shows that both plotting and reworking contribute to the model's superiority.


title: Higgs rare and exotic decays

abstract: Several recent searches for exotic and rare decays of the Standard Model Higgs boson with the ATLAS and CMS detectors are presented. The searches are performed on $\sqrt{s}$=13 TeV proton-proton collisions data collected at the LHC between 2015 and 2018. The topics covered include searches for the Higgs boson decays into two pseudoscalars, $H\rightarrow aa$, in three different final states, $2b2\mu$, $4b$ and $2\mu2\tau$, search for lepton-flavour violating Higgs decays, $H\rightarrow \mu\tau / e\tau$, and search for a rare $H\rightarrow \ell\ell \gamma$ decay. ATLAS presents evidence for the $H\rightarrow \ell\ell \gamma$ rare decay, amounting to an observed significance of $3.2\sigma$.


title: Radiative corrections of order $O(\alpha E_e/m_N)$ to Sirlin's radiative corrections of order $O(\alpha/\pi)$, induced by the hadronic structure of the neutron

abstract: We investigate the contributions of the hadronic structure of the neutron to radiative $O(\alpha E_e/m_N)$ corrections (or the inner $O(\alpha E_e/m_N)$ RC) to the neutron beta decay, where $\alpha$, $E_e$ and $m_N$ are the fine-structure constant, the electron energy and the nucleon mass, respectively. We perform the calculation within the effective quantum field theory of strong low-energy pion-nucleon interactions described by the linear $\sigma$-model with chiral $SU(2) \times SU(2)$ symmetry and electroweak hadron-hadron, hadron-lepton and lepton-lepton interactions for the electron-lepton family with $SU(2)_L \times U(1)_Y$ symmetry of the Standard Electroweak Theory (Ivanov et al., Phys. Rev. D99, 093006 (2019)). We show that after renormalization, carried out in accordance with Sirlin's prescription (Sirlin, Phys. Rev. 164, 1767 (1967)), the inner $O(\alpha E_e/m_N)$ RC are of the order of a few parts of $10^{-5} - 10^{-4}$. This agrees well with the results obtained in (Ivanov et al., Phys. Rev. D99, 093006 (2019)).


title: A complete ab initio view of Orbach and Raman spin-lattice relaxation in a Dysprosium coordination compound

abstract: The unique electronic and magnetic properties of Lanthanides molecular complexes place them at the forefront of the race towards high-temperature single-ion magnets and magnetic quantum bits. The design of compounds of this class has so far been almost exclusively driven by static crystal field considerations, with emphasis on increasing the magnetic anisotropy barrier. This guideline has now reached its maximum potential and new progress can only come from a deeper understanding of spin-phonon relaxation mechanisms. In this work we compute relaxation times fully ab initio and unveil the nature of all spin-phonon relaxation mechanisms, namely Orbach and Raman pathways, in a prototypical Dy single-ion magnet. Computational predictions are in agreement with the experimental determination of spin relaxation time and crystal field anisotropy, and show that Raman relaxation, dominating at low temperature, is triggered by low-energy phonons and little affected by further engineering of crystal field axiality. A comprehensive analysis of spin-phonon coupling mechanism reveals that molecular vibrations beyond the ion's first coordination shell can also assume a prominent role in spin relaxation through an electrostatic polarization effect. Therefore, this work shows the way forward in the field by delivering a novel and complete set of chemically-sound design rules tackling every aspect of spin relaxation at any temperature


title: Studying the $b\rightarrow s \ell^+\ell^-$ anomalies and $(g-2)_{\mu}$ in $R$-parity violating MSSM framework with the inverse seesaw mechanism

abstract: Inspired by the recent experimental results which show deviations from the standard model (SM) predictions of $b\rightarrow s \ell^+\ell^-$ transitions, we study the $R$-parity violating minimal supersymmetric standard model (RPV-MSSM) extended by the inverse seesaw mechanism. The trilinear $R$-parity violating terms, together with the chiral mixing of sneutrinos, induce the loop contributions to the $b\rightarrow s \ell^+\ell^-$ anomaly. We study the parameter space of the single-parameter scenario $C^{\rm NP}_{9,\mu}=-C^{\rm NP}_{10,\mu}=C_{\rm V}$ and the double-parameter scenario $(C_{\rm V},C_{\rm U})$, respectively, constrained by other experimental data such as $B_s-\bar{B}_s$ mixing, $B\rightarrow X_s \gamma$ decay, the lepton flavour violating decays, etc. Both the single-parameter and the double-parameter scenario can resolve the long existing muon anomalous magnetic moment problem as well, and allow the anomalous $t\rightarrow cg$ process to reach the sensitivity at the Future Circular hadron-hadron Collider (FCC-hh).


title: On the enumeration of plane bipolar posets and transversal structures

abstract: We show that plane bipolar posets (i.e., plane bipolar orientations with no transitive edge) and transversal structures can be set in correspondence to certain (weighted) models of quadrant walks, via suitable specializations of a bijection due to Kenyon, Miller, Sheffield and Wilson. We then derive exact and asymptotic counting results. In particular we prove (computationally and then bijectively) that the number of plane bipolar posets on $n+2$ vertices equals the number of plane permutations of size $n$. Regarding transversal structures, for each $v\geq 0$ we consider $t_n(v)$ the number of such structures with $n+4$ vertices and weight $v$ per quadrangular inner face (the case $v=0$ corresponds to having only triangular inner faces). We obtain a recurrence to compute $t_n(v)$, and an asymptotic formula that for $v=0$ gives $t_n(0)\sim c\ \!(27/2)^nn^{-1-\pi/\mathrm{arccos}(7/8)}$ for some $c>0$, which also ensures that the associated generating function is not D-finite.


title: Information-theoretic Evolution of Model Agnostic Global Explanations

abstract: Explaining the behavior of black box machine learning models through human interpretable rules is an important research area. Recent work has focused on explaining model behavior locally i.e. for specific predictions as well as globally across the fields of vision, natural language, reinforcement learning and data science. We present a novel model-agnostic approach that derives rules to globally explain the behavior of classification models trained on numerical and/or categorical data. Our approach builds on top of existing local model explanation methods to extract conditions important for explaining model behavior for specific instances followed by an evolutionary algorithm that optimizes an information theory based fitness function to construct rules that explain global model behavior. We show how our approach outperforms existing approaches on a variety of datasets. Further, we introduce a parameter to evaluate the quality of interpretation under the scenario of distributional shift. This parameter evaluates how well the interpretation can predict model behavior for previously unseen data distributions. We show how existing approaches for interpreting models globally lack distributional robustness. Finally, we show how the quality of the interpretation can be improved under the scenario of distributional shift by adding out of distribution samples to the dataset used to learn the interpretation and thereby, increase robustness. All of the datasets used in our paper are open and publicly available. Our approach has been deployed in a leading digital marketing suite of products.


title: On the absolute convergence of automorphic Dirichlet series

abstract: Let $F(s)=\sum_{n=1}^{\infty}\frac{a_n}{n^s}$ be a Dirichlet series in the axiomatically defined class ${\mathfrak A}^{\#}$ . The class ${\mathfrak A}^{\#}$ is known to contain the extended Selberg class ${\mathcal S}^{\#}$, as well as all the $L$-functions of automorphic forms on $GL_n/K$, where $K$ is a number field. Let $d$ be the degree of $F(s)$. We show that $\sum_{n<X}|a_n|=\Omega(X^{\frac{1}{2}+\frac{1}{2d}})$, and hence, that the abscissa of absolute convergence of $\sigma_a$ of $F(s)$ must satisfy $\sigma_a\ge 1/2+1/2d$.


title: A Hypothesis Testing Approach to Nonstationary Source Separation

abstract: The extraction of nonstationary signals from blind and semi-blind multivariate observations is a recurrent problem. Numerous algorithms have been developed for this problem, which are based on the exact or approximate joint diagonalization of second or higher order cumulant matrices/tensors of multichannel data. While a great body of research has been dedicated to joint diagonalization algorithms, the selection of the diagonalized matrix/tensor set remains highly problem-specific. Herein, various methods for nonstationarity identification are reviewed and a new general framework based on hypothesis testing is proposed, which results in a classification/clustering perspective to semi-blind source separation of nonstationary components. The proposed method is applied to noninvasive fetal ECG extraction, as case study.


title: Many-body theory of positron binding in polyatomic molecules

abstract: Positrons bind to molecules leading to vibrational excitation and spectacularly enhanced annihilation. Whilst positron binding energies have been measured via resonant annihilation spectra for $\sim$90 molecules in the past two decades, an accurate \emph{ab initio} theoretical description has remained elusive. Of the molecules studied experimentally, calculations exist for only 6, and for these, standard quantum chemistry approaches have proved severely deficient, agreeing with experiment to at best 25% accuracy for polar molecules, and failing to predict binding in nonpolar molecules. The mechanisms of binding are not understood. Here, we develop a many-body theory of positron-molecule interactions and uncover the role of strong many-body correlations including polarization of the electron cloud, screening of the positron-electron Coulomb interaction by molecular electrons, and crucially, the unique non-perturbative process of virtual-positronium formation (where a molecular electron temporarily tunnels to the positron): they dramatically enhance binding in polar molecules and enable binding in nonpolars. We also elucidate the role of individual molecular orbitals, highlighting the importance of electronic $\pi$ bonds. Overall, we calculate binding energies in agreement with experiment (to within 1% in cases), and we predict binding in formamide and nucleobases. As well as supporting resonant annihilation experiments and positron-based molecular spectroscopy, the approach can be extended to positron scattering and annihilation $\gamma$ spectra in molecules and condensed matter, to provide fundamental insight and predictive capability required to properly interpret materials science diagnostics, develop antimatter-based technologies (including positron traps, beams and positron emission tomography), and understand positrons in the galaxy.

categories: physics.atom-ph cond-mat.str-el physics.atm-clus physics.comp-ph quant-ph


title: Thompson Sampling for Gaussian Entropic Risk Bandits

abstract: The multi-armed bandit (MAB) problem is a ubiquitous decision-making problem that exemplifies exploration-exploitation tradeoff. Standard formulations exclude risk in decision making. Risknotably complicates the basic reward-maximising objectives, in part because there is no universally agreed definition of it. In this paper, we consider an entropic risk (ER) measure and explore the performance of a Thompson sampling-based algorithm ERTS under this risk measure by providing regret bounds for ERTS and corresponding instance dependent lower bounds.


title: Ready When You Are: Efficient Condition Variables via Delegated Condition Evaluation

abstract: Multi-thread applications commonly utilize condition variables for communication between threads. Condition variables allow threads to block and wait until a certain condition holds, and also enable threads to wake up their blocked peers notifying them about a change to the state of shared data. Quite often such notifications are delivered to all threads, while only a small number of specific threads is interested in it. This results in so-called futile wakeups, where threads receiving the notification wake up and resume their execution only to realize that the condition they are waiting for does not hold and they need to wait again. Those wakeups cause numerous context switches, increase lock contention and cache pressure, translating into lots of wasted computing cycles and energy. In this work, we propose to delegate conditions on which threads are waiting to the thread sending notifications. This enables the latter to evaluate the conditions and send the notification(s) only to the relevant thread(s), practically eliminating futile wakeups altogether. Our initial evaluation of this idea shows promising results, achieving 3-4x throughput improvement over legacy condition variables.


title: A remark on Gersten complex for Milnor $K$-theory

abstract: In this note, we consider the Gersten complex for Milnor $K$-theory over a regular local Henselian domain $S$ and prove that in degrees $\geq \dim S\geq 1$, the Gersten complex of an essentially smooth Henselian local $S$-scheme is exact.


title: The elastic vacuum

abstract: The quantum gravity vacuum must contain virtual fluctuations of black hole microstates. These extended-size fluctuations get `crushed' when a closed trapped surface forms, and turn into on-shell `fuzzball' states that resolve the information puzzle. We argue that these same fluctuations can get `stretched' by the anti-trapped surfaces in an expanding cosmology, and that this stretching generates vacuum energy. The stretching happen when the Hubble deceleration reduces quickly, which happens whenever the pressure drops quickly. We thus get an inflation-scale vacuum energy when the heavy GUTS particles become nonrelativistic, and again a small vacuum energy when the radiation phase turns to dust. The expansion law in the radiation phase does not allow stretching, in agreement with the observed irrelevance of vacuum energy in that phase. The extra energy induced when the radiation phase changes to dust may explain the tension in the Hubble constant between low and high redshift data.


title: BNNpriors: A library for Bayesian neural network inference with different prior distributions

abstract: Bayesian neural networks have shown great promise in many applications where calibrated uncertainty estimates are crucial and can often also lead to a higher predictive performance. However, it remains challenging to choose a good prior distribution over their weights. While isotropic Gaussian priors are often chosen in practice due to their simplicity, they do not reflect our true prior beliefs well and can lead to suboptimal performance. Our new library, BNNpriors, enables state-of-the-art Markov Chain Monte Carlo inference on Bayesian neural networks with a wide range of predefined priors, including heavy-tailed ones, hierarchical ones, and mixture priors. Moreover, it follows a modular approach that eases the design and implementation of new custom priors. It has facilitated foundational discoveries on the nature of the cold posterior effect in Bayesian neural networks and will hopefully catalyze future research as well as practical applications in this area.


title: Counterfactual Interventions Reveal the Causal Effect of Relative Clause Representations on Agreement Prediction

abstract: When language models process syntactically complex sentences, do they use their representations of syntax in a manner that is consistent with the grammar of the language? We propose AlterRep, an intervention-based method to address this question. For any linguistic feature of a given sentence, AlterRep generates counterfactual representations by altering how the feature is encoded, while leaving intact all other aspects of the original representation. By measuring the change in a model's word prediction behavior when these counterfactual representations are substituted for the original ones, we can draw conclusions about the causal effect of the linguistic feature in question on the model's behavior. We apply this method to study how BERT models of different sizes process relative clauses (RCs). We find that BERT variants use RC boundary information during word prediction in a manner that is consistent with the rules of English grammar; this RC boundary information generalizes to a considerable extent across different RC types, suggesting that BERT represents RCs as an abstract linguistic category.


title: A Spatio-Temporal Model for Predicting Wind Speeds in Southern California

abstract: The share of wind power in fuel mixes worldwide has increased considerably. The main ingredient when deriving wind power predictions are wind speed data; the closer to the wind farms, the better they forecast the power supply. The current paper proposes a hybrid model for predicting wind speeds at convenient locations. It is then applied to Southern California power price area. We build random fields with time series of gridded historical forecasts and actual wind speed observations. We estimate with ordinary kriging the spatial variability of the temporal parameters and derive predictions. The advantages of this work are twofold: (1) an accurate daily wind speed forecast at any location in the area and (2) a general method applicable to other markets.


title: Open-set Face Recognition for Small Galleries Using Siamese Networks

abstract: Face recognition has been one of the most relevant and explored fields of Biometrics. In real-world applications, face recognition methods usually must deal with scenarios where not all probe individuals were seen during the training phase (open-set scenarios). Therefore, open-set face recognition is a subject of increasing interest as it deals with identifying individuals in a space where not all faces are known in advance. This is useful in several applications, such as access authentication, on which only a few individuals that have been previously enrolled in a gallery are allowed. The present work introduces a novel approach towards open-set face recognition focusing on small galleries and in enrollment detection, not identity retrieval. A Siamese Network architecture is proposed to learn a model to detect if a face probe is enrolled in the gallery based on a verification-like approach. Promising results were achieved for small galleries on experiments carried out on Pubfig83, FRGCv1 and LFW datasets. State-of-the-art methods like HFCN and HPLS were outperformed on FRGCv1. Besides, a new evaluation protocol is introduced for experiments in small galleries on LFW.


title: On Gravitational Stefan-Boltzmann Law and Casimir Effect in FRW Universe

abstract: Both Stefan-Boltzmann law and the Casimir effect, in a universe described by the FRW metric with zero curvature, are calculated. These effects are described by Thermo Field Dynamics (TFD). The gravitational energy-momentum tensor is defined in the context of Teleparallel Equivalent to General Relativity (TEGR). Each of the two effects gives a consistent prediction with what is observed on a cosmological scale. One of the effect establishes a minimum range for the deceleration parameter. While another leads to the conclusion that a possible cosmological constant has a very small order of magnitude.


title: On the continuous dual Hahn process

abstract: In this note we extend the continuous dual Hahn process constructed by Corwin and Knizel on a finite time interval to the entire real line by taking a limit of a closely related Markov process. We also characterize this Markov processes by conditional means and variances under bidirectional conditioning.


title: VEGAS: A VST Early-type GAlaxy Survey.VI. The diffuse light in HCG 86 from the ultra-deep VEGAS images

abstract: Context. In this paper we present ultra deep images of the compact group of galaxies HCG 86 as part of the VEGAS survey. Aims. Our main goals are to estimate the amount of intragroup light (IGL), to study the light and color distributions in order to address the main formation process of the IGL component in groups of galaxies. Methods. We derived the azimuthally averaged surface brightness profiles in the g,r and i bands with g - r and r - i average colors and color profiles for all group members. By fitting the light distribution, we have extrapolated the contribution of the stellar halos plus the diffuse light from the brightest component of each galaxy. The results are compared with theoretical predictions. Results. The long integration time and wide area covered make our data deeper than previous literature studies of the IGL in compact groups of galaxies and allow us to produce an extended (~160 kpc) map of the IGL, down to a surface brightness level of about 30 mag/arcsec^2 in the g band. The IGL in HCG 86 is mainly in diffuse form and has average colors of g - r ~ 0.8 mag and r - i ~ 0.4 mag. The fraction of IGL in HCG 86 is ~ 16% of the total luminosity of the group, and this is consistent with estimates available for other compact groups and loose groups of galaxies of similar virial masses. A weak trend is present between the amount of IGL and the early-type to late-type galaxy ratio. Conclusions. By comparing the IGL fraction and colors with those predicted by simulations, the amount of IGL in HCG 86 would be the result of the disruption of satellites at an epoch of z ~ 0.4. At this redshift, observed colors are consistent with the scenario where the main contribution to the mass of the IGL comes from the intermediate-massive galaxies.


title: Dynamical mean-field theory of the Anderson-Hubbard model with local and non-local disorder in tensor formulation

abstract: To explore correlated electrons in the presence of local and non-local disorder, the Blackman-Esterling-Berk method for averaging over off-diagonal disorder is implemented into dynamical mean-field theory using tensor notation. The impurity model combining disorder and correlations is solved using the recently developed fork tensor-product state solver, which allows one to calculate the single particle spectral functions on the real-frequency axis. In the absence of off-diagonal hopping, we establish exact bounds of the spectral function of the non-interacting Bethe lattice with coordination number $Z$. In the presence of interaction, the Mott insulating paramagnetic phase of the one-band Hubbard model is computed at zero temperature in alloys with site- and off-diagonal disorder. When the Hubbard $U$ parameter is increased, transitions from an alloy band-insulator through a correlated metal into a Mott insulating phase are found to take place.


title: Gravitational waves from spectator Gauge-flation

abstract: We investigate the viability of inflation with a spectator sector comprised of non-Abelian gauge fields coupled through a higher order operator. We dub this model "spectator Gauge-flation". We study the predictions for the amplitude and tensor tilt of chiral gravitational waves and conclude that a slightly red-tilted tensor power spectrum is preferred $n_T=-{\cal O}(0.01)$. As with related models, the enhancement of chiral gravitational waves with respect to the single-field vacuum gravitational wave background is controlled by the parameter $\gamma=g^2 Q^2/H^2$, where $g$ is the gauge coupling, $H$ is the Hubble scale and $Q$ is the VEV of the $SU(2)$ sector. The requirement that the $SU(2)$ is a spectator sector leads to a maximum allowed value for $\gamma$, thereby constraining the possible amplification. In order to provide concrete predictions, we use an $\alpha$-attractor T-model potential for the inflaton sector. Potential observation of chiral gravitational waves with significantly tilted tensor spectra would then indicate the presence of additional couplings of the gauge fields to axions, like in the spectator axion-SU(2) model, or additional gauge field operators.


title: Multiparty Session Types for Safe Runtime Adaptation in an Actor Language (Extended version)

abstract: Human fallibility, unpredictable operating environments, and the heterogeneity of hardware devices are driving the need for software to be able to adapt as seen in the Internet of Things or telecommunication networks. Unfortunately, mainstream programming languages do not readily allow a software component to sense and respond to its operating environment, by discovering, replacing, and communicating with components that are not part of the original system design, while maintaining static correctness guarantees. In particular, if a new component is discovered at runtime, there is no guarantee that its communication behaviour is compatible with existing components. We address this problem by using multiparty session types with explicit connection actions, a type formalism used to model distributed communication protocols. By associating session types with software components, the discovery process can check protocol compatibility and, when required, correctly replace components without jeapordising safety. We present the design and implementation of EnsembleS, the first actor-based language with adaptive features and a static session type system, and apply it to a case study based on an adaptive DNS server. We formalise the type system of EnsembleS and prove the safety of well-typed programs, making essential use of recent advances in non-classical multiparty session types.


title: A Survey of Security Vulnerabilities in Ethereum Smart Contracts

abstract: Ethereum Smart Contracts based on Blockchain Technology (BT)enables monetary transactions among peers on a blockchain network independent of a central authorizing agency. Ethereum smart contracts are programs that are deployed as decentralized applications, having the building blocks of the blockchain consensus protocol. This enables consumers to make agreements in a transparent and conflict-free environment. However, there exist some security vulnerabilities within these smart contracts that are a potential threat to the applications and their consumers and have shown in the past to cause huge financial losses. In this study, we review the existing literature and broadly classify the BT applications. As Ethereum smart contracts find their application mostly in e-commerce applications, we believe these are more commonly vulnerable to attacks. In these smart contracts, we mainly focus on identifying vulnerabilities that programmers and users of smart contracts must avoid. This paper aims at explaining eight vulnerabilities that are specific to the application level of BT by analyzing the past exploitation case scenarios of these security vulnerabilities. We also review some of the available tools and applications that detect these vulnerabilities in terms of their approach and effectiveness. We also investigated the availability of detection tools for identifying these security vulnerabilities and lack thereof to identify some of them


title: Saddle point preconditioners for weak-constraint 4D-Var

abstract: Data assimilation algorithms combine information from observations and prior model information to obtain the most likely state of a dynamical system. The linearised weak-constraint four-dimensional variational assimilation problem can be reformulated as a saddle point problem, which admits more scope for preconditioners than the primal form. In this paper we design new terms which can be used within existing preconditioners, such as block diagonal and constraint-type preconditioners. Our novel preconditioning approaches: (i) incorporate model information, and (ii) are designed to target correlated observation error covariance matrices. To our knowledge (i) has not previously been considered for data assimilation problems. We develop new theory demonstrating the effectiveness of the new preconditioners within Krylov subspace methods. Linear and non-linear numerical experiments reveal that our new approach leads to faster convergence than existing state-of-the-art preconditioners for a broader range of problems than indicated by the theory alone. We present a range of numerical experiments performed in serial.


title: Sizable triple Higgs couplings in the 2HDM: Prospects for future $e^+e^-$ colliders

abstract: In the framework of the $\mathcal{CP}$ conserving Two Higgs Doublet Model (2HDM), type I and II, we study the triple Higgs couplings with at least one light $h$ Higgs boson that is identified by the 125 GeV Higgs boson. We define benchmark planes that exhibit large values of triple Higgs couplings, while being in agreement with all experimental and theoretical constraints. Finally, we analyze the impact of the triple Higgs couplings on the production cross section of two neutral Higgs bosons in two channels, $\sigma(e^+e^-\to h_i h_j Z)$ and $\sigma(e^+e^- \to h_i h_j \nu\bar{\nu})$ with $h_i h_j = hh, hH, HH, AA$. We show that the triple Higgs couplings have an important impact on these $e^+e^-$ production cross sections.


title: Do Context-Aware Translation Models Pay the Right Attention?

abstract: Context-aware machine translation models are designed to leverage contextual information, but often fail to do so. As a result, they inaccurately disambiguate pronouns and polysemous words that require context for resolution. In this paper, we ask several questions: What contexts do human translators use to resolve ambiguous words? Are models paying large amounts of attention to the same context? What if we explicitly train them to do so? To answer these questions, we introduce SCAT (Supporting Context for Ambiguous Translations), a new English-French dataset comprising supporting context words for 14K translations that professional translators found useful for pronoun disambiguation. Using SCAT, we perform an in-depth analysis of the context used to disambiguate, examining positional and lexical characteristics of the supporting words. Furthermore, we measure the degree of alignment between the model's attention scores and the supporting context from SCAT, and apply a guided attention strategy to encourage agreement between the two.


title: NICER Detection of Thermal X-ray Pulsations from the Massive Millisecond Pulsars PSR J0740+6620 and PSR J1614-2230

abstract: We report the detection of X-ray pulsations from the rotation-powered millisecond-period pulsars PSR J0740+6620 and PSR J1614-2230, two of the most massive neutron stars known, using observations with the Neutron Star Interior Composition Explorer (NICER). We also analyze X-ray Multi-Mirror Mission (XMM-Newton) data for both pulsars to obtain their time-averaged fluxes and study their respective X-ray fields. PSR J0740+6620 exhibits a broad double-peaked profile with a separation of ~0.4 in phase. PSR J1614-2230, on the other hand, has a broad single-peak profile. The broad modulations with soft X-ray spectra of both pulsars are indicative of thermal radiation from one or more small regions of the stellar surface. We show the NICER detections of X-ray pulsations for both pulsars and also discuss the phase relationship to their radio pulsations. In the case of PSR J0740+6620, this paper documents the data reduction performed to obtain the pulsation detection and prepare for pulse profile modeling analysis.


title: The Radius of PSR J0740+6620 from NICER and XMM-Newton Data

abstract: PSR J0740$+$6620 has a gravitational mass of $2.08\pm 0.07~M_\odot$, which is the highest reliably determined mass of any neutron star. As a result, a measurement of its radius will provide unique insight into the properties of neutron star core matter at high densities. Here we report a radius measurement based on fits of rotating hot spot patterns to Neutron Star Interior Composition Explorer (NICER) and X-ray Multi-Mirror (XMM-Newton) X-ray observations. We find that the equatorial circumferential radius of PSR J0740$+$6620 is $13.7^{+2.6}_{-1.5}$ km (68%). We apply our measurement, combined with the previous NICER mass and radius measurement of PSR J0030$+$0451, the masses of two other $\sim 2~M_\odot$ pulsars, and the tidal deformability constraints from two gravitational wave events, to three different frameworks for equation of state modeling, and find consistent results at $\sim 1.5-3$ times nuclear saturation density. For a given framework, when all measurements are included the radius of a $1.4~M_\odot$ neutron star is known to $\pm 4$% (68% credibility) and the radius of a $2.08~M_\odot$ neutron star is known to $\pm 5$%. The full radius range that spans the $\pm 1\sigma$ credible intervals of all the radius estimates in the three frameworks is $12.45\pm 0.65$ km for a $1.4~M_\odot$ neutron star and $12.35\pm 0.75$ km for a $2.08~M_\odot$ neutron star.


title: A NICER View of the Massive Pulsar PSR J0740+6620 Informed by Radio Timing and XMM-Newton Spectroscopy

abstract: We report on Bayesian estimation of the radius, mass, and hot surface regions of the massive millisecond pulsar PSR J0740$+$6620, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer X-ray Timing Instrument (NICER XTI) event data. We condition on informative pulsar mass, distance, and orbital inclination priors derived from the joint NANOGrav and CHIME/Pulsar wideband radio timing measurements of arXiv:2104.00880. We use XMM European Photon Imaging Camera spectroscopic event data to inform our X-ray likelihood function. The prior support of the pulsar radius is truncated at 16 km to ensure coverage of current dense matter models. We assume conservative priors on instrument calibration uncertainty. We constrain the equatorial radius and mass of PSR J0740$+$6620 to be $12.39_{-0.98}^{+1.30}$ km and $2.072_{-0.066}^{+0.067}$ M$_{\odot}$ respectively, each reported as the posterior credible interval bounded by the 16% and 84% quantiles, conditional on surface hot regions that are non-overlapping spherical caps of fully-ionized hydrogen atmosphere with uniform effective temperature; a posteriori, the temperature is $\log_{10}(T$ [K]$)=5.99_{-0.06}^{+0.05}$ for each hot region. All software for the X-ray modeling framework is open-source and all data, model, and sample information is publicly available, including analysis notebooks and model modules in the Python language. Our marginal likelihood function of mass and equatorial radius is proportional to the marginal joint posterior density of those parameters (within the prior support) and can thus be computed from the posterior samples.


title: Constraints on the dense matter equation of state and neutron star properties from NICER's mass-radius estimate of PSR J0740+6620 and multimessenger observations

abstract: In recent years our understanding of the dense matter equation of state (EOS) of neutron stars has significantly improved by analyzing multimessenger data from radio/X-ray pulsars, gravitational wave events, and from nuclear physics constraints. Here we study the additional impact on the EOS from the jointly estimated mass and radius of PSR J0740+6620, presented in Riley et al. (2021) by analyzing a combined dataset from X-ray telescopes NICER and XMM-Newton. We employ two different high-density EOS parameterizations: a piecewise-polytropic (PP) model and a model based on the speed of sound in a neutron star (CS). At nuclear densities these are connected to microscopic calculations of neutron matter based on chiral effective field theory interactions. In addition to the new NICER data for this heavy neutron star, we separately study constraints from the radio timing mass measurement of PSR J0740+6620, the gravitational wave events of binary neutron stars GW190425 and GW170817, and for the latter the associated kilonova AT2017gfo. By combining all these, and the NICER mass-radius estimate of PSR J0030+0451 we find the radius of a 1.4 solar mass neutron star to be constrained to the 95% credible ranges 12.33^{+0.76}_{-0.81} km (PP model) and 12.18^{+0.56}_{-0.79} km (CS model). In addition, we explore different chiral effective field theory calculations and show that the new NICER results provide tight constraints for the pressure of neutron star matter at around twice saturation density, which shows the power of these observations to constrain dense matter interactions at intermediate densities.


title: EASE: Extractive-Abstractive Summarization with Explanations

abstract: Current abstractive summarization systems outperform their extractive counterparts, but their widespread adoption is inhibited by the inherent lack of interpretability. To achieve the best of both worlds, we propose EASE, an extractive-abstractive framework for evidence-based text generation and apply it to document summarization. We present an explainable summarization system based on the Information Bottleneck principle that is jointly trained for extraction and abstraction in an end-to-end fashion. Inspired by previous research that humans use a two-stage framework to summarize long documents (Jing and McKeown, 2000), our framework first extracts a pre-defined amount of evidence spans as explanations and then generates a summary using only the evidence. Using automatic and human evaluations, we show that explanations from our framework are more relevant than simple baselines, without substantially sacrificing the quality of the generated summary.


title: Fast post-adiabatic waveforms in the time domain: Applications to compact binary coalescences in LIGO and Virgo

abstract: We present a computationally efficient (time-domain) multipolar waveform model for quasi-circular spin-aligned compact binary coalescences. The model combines the advantages of the numerical-relativity informed, effective-one-body (EOB) family of models with a post-adiabatic solution of the equations of motion for the inspiral part of the two-body dynamics. We benchmark this model against other state-of-the-art waveforms in terms of efficiency and accuracy. We find a speed-up of one to two orders of magnitude compared to the underlying time-domain EOB model for the total mass range $2 - 100 M_{\odot}$. More specifically, for a low total-mass system, such as a binary neutron star with equal masses of $1.4 M_{\odot}$, like GW170817, the computational speedup is around 100 times; for an event with total mass $\sim 40 M_\odot$ and mass ratio $\sim 3$, like GW190412, the speedup is by a factor of $\sim 20$, while for a binary system of comparable masses and total mass of $\sim 70 M_{\odot}$, like GW150914, it is by a factor of $\sim 10$. We demonstrate that the new model is extremely faithful to the underlying EOB model with unfaithfulness less than $0.01\%$ across the entire applicable region of parameter space. Finally, we present successful applications of this new waveform model to parameter estimation studies and tests of general relativity.


title: Laser-induced electron emission from Au nanowires: A probe for orthogonal polarizations

abstract: Photoelectron field emission, induced by femtosecond laser pulses focused on metallic nanotips, provides spatially coherent and temporally short electron pulses. The properties of the photoelectron yield give insight into both the material properties of the nanostructure and the exciting laser focus. Ultralong nanoribbons, grown as a single crystal attached to a metallic taper, are sources of electron field emission that have not yet been characterized. In this report, photoemission from gold nanoribbon samples is studied and compared to emission from tungsten and gold tips. We observe that the emission from sharp tips generally depends on one transverse component of the exciting laser field, while the emission of a blunted nanoribbon is found to be sensitive to both components. We propose that this property makes photoemission from nanoribbons a candidate for position-sensitive detection of the longitudinal field component in a tightly focused beam.


title: Spatial dynamics of a population in a heterogeneous environment

abstract: We consider a certain lattice branching random walk with on-site competition and in an environment which is heterogeneous at a macroscopic scale $1/\varepsilon$ in space and time. This can be seen as a model for the spatial dynamics of a biological population in a habitat which is heterogeneous at a large scale (mountains, temperature or precipitation gradient...). The model incorporates another parameter, $K$, which is a measure of the local population density. We study the model in the limit when first $\varepsilon\to 0$ and then $K\to\infty$. In this asymptotic regime, we show that the rescaled position of the front as a function of time converges to the solution of an explicit ODE. We further discuss the relation with another popular model of population dynamics, the Fisher-KPP equation, which arises in the limit $K\to\infty$. Combined with known results on the Fisher-KPP equation, our results show in particular that the limits $\varepsilon\to0$ and $K\to\infty$ do not commute in general. We conjecture that an interpolating regime appears when $\log K$ and $1/\varepsilon$ are of the same order.


title: Evaluating the Robustness of Self-Supervised Learning in Medical Imaging

abstract: Self-supervision has demonstrated to be an effective learning strategy when training target tasks on small annotated data-sets. While current research focuses on creating novel pretext tasks to learn meaningful and reusable representations for the target task, these efforts obtain marginal performance gains compared to fully-supervised learning. Meanwhile, little attention has been given to study the robustness of networks trained in a self-supervised manner. In this work, we demonstrate that networks trained via self-supervised learning have superior robustness and generalizability compared to fully-supervised learning in the context of medical imaging. Our experiments on pneumonia detection in X-rays and multi-organ segmentation in CT yield consistent results exposing the hidden benefits of self-supervision for learning robust feature representations.


title: Scaling Ensemble Distribution Distillation to Many Classes with Proxy Targets

abstract: Ensembles of machine learning models yield improved system performance as well as robust and interpretable uncertainty estimates; however, their inference costs may often be prohibitively high. \emph{Ensemble Distribution Distillation} is an approach that allows a single model to efficiently capture both the predictive performance and uncertainty estimates of an ensemble. For classification, this is achieved by training a Dirichlet distribution over the ensemble members' output distributions via the maximum likelihood criterion. Although theoretically principled, this criterion exhibits poor convergence when applied to large-scale tasks where the number of classes is very high. In our work, we analyze this effect and show that the Dirichlet log-likelihood criterion classes with low probability induce larger gradients than high-probability classes. This forces the model to focus on the distribution of the ensemble tail-class probabilities. We propose a new training objective that minimizes the reverse KL-divergence to a \emph{Proxy-Dirichlet} target derived from the ensemble. This loss resolves the gradient issues of Ensemble Distribution Distillation, as we demonstrate both theoretically and empirically on the ImageNet and WMT17 En-De datasets containing 1000 and 40,000 classes, respectively.


title: Automatic Non-Linear Video Editing Transfer

abstract: We propose an automatic approach that extracts editing styles in a source video and applies the edits to matched footage for video creation. Our Computer Vision based techniques considers framing, content type, playback speed, and lighting of each input video segment. By applying a combination of these features, we demonstrate an effective method that automatically transfers the visual and temporal styles from professionally edited videos to unseen raw footage. We evaluated our approach with real-world videos that contained a total of 3872 video shots of a variety of editing styles, including different subjects, camera motions, and lighting. We reported feedback from survey participants who reviewed a set of our results.


title: Phase Space Crystal Vibrations: Chiral Edge States with Preserved Time-reversal Symmetry

abstract: It was recently discovered that atoms subject to a time-periodic drive can give rise to a crystal structure in phase space. In this work, we point out the atom-atom interactions give rise to collective phonon excitations of phase-space crystal via a pairing interaction with intrinsically complex phases that can lead to a phonon Chern insulator, accompanied by topologically robust chiral transport along the edge of the phase-space crystal. This topological phase is realized even in scenarios where the time-reversal transformation is a symmetry, which is surprising because the breaking of time-reversal symmetry is a strict precondition for topological chiral transport in the standard setting of real-space crystals. Our work has also important implications for the dynamics of 2D charged particles in a strong magnetic field.


title: BERT Busters: Outlier Dimensions that Disrupt Transformers

abstract: Multiple studies have shown that Transformers are remarkably robust to pruning. Contrary to this received wisdom, we demonstrate that pre-trained Transformer encoders are surprisingly fragile to the removal of a very small number of features in the layer outputs (<0.0001% of model weights). In case of BERT and other pre-trained encoder Transformers, the affected component is the scaling factors and biases in the LayerNorm. The outliers are high-magnitude normalization parameters that emerge early in pre-training and show up consistently in the same dimensional position throughout the model. We show that disabling them significantly degrades both the MLM loss and the downstream task performance. This effect is observed across several BERT-family models and other popular pre-trained Transformer architectures, including BART, XLNet and ELECTRA; we also show a similar effect in GPT-2.


title: Structural formulas for matrix-valued orthogonal polynomials related to $2\times 2$ hypergeometric operators

abstract: We give some structural formulas for the family of matrix-valued orthogonal polynomials of size $2\times 2$ introduced by C. Calder\'on et al. in an earlier work, which are common eigenfunctions of a differential operator of hypergeometric type. Specifically, we give a Rodrigues formula that allows us to write this family of polynomials explicitly in terms of the classical Jacobi polynomials, and write, for the sequence of orthonormal polynomials, the three-term recurrence relation and the Christoffel-Darboux identity. We obtain a Pearson equation, which enables us to prove that the sequence of derivatives of the orthogonal polynomials is also orthogonal, and to compute a Rodrigues formula for these polynomials as well as a matrix-valued differential operator having these polynomials as eigenfunctions. We also describe the second-order differential operators of the algebra associated with the weight matrix.


title: Generalized LR-drawings of trees

abstract: The LR-drawing-method is a method of drawing an ordered rooted binary tree based on drawing one root-to-leaf path on a vertical line and attaching recursively obtained drawings of the subtrees on the left and right. In this paper, we study how to generalize this drawing-method to trees of higher arity. We first prove that (with some careful modifications) the proof of existence of a special root-to-leaf path transfers to trees of higher arity. Then we use such paths to obtain generalized LR-drawings of trees of arbitrary arity.


title: Omnimatte: Associating Objects and Their Effects in Video

abstract: Computer vision is increasingly effective at segmenting objects in images and videos; however, scene effects related to the objects -- shadows, reflections, generated smoke, etc -- are typically overlooked. Identifying such scene effects and associating them with the objects producing them is important for improving our fundamental understanding of visual scenes, and can also assist a variety of applications such as removing, duplicating, or enhancing objects in video. In this work, we take a step towards solving this novel problem of automatically associating objects with their effects in video. Given an ordinary video and a rough segmentation mask over time of one or more subjects of interest, we estimate an omnimatte for each subject -- an alpha matte and color image that includes the subject along with all its related time-varying scene elements. Our model is trained only on the input video in a self-supervised manner, without any manual labels, and is generic -- it produces omnimattes automatically for arbitrary objects and a variety of effects. We show results on real-world videos containing interactions between different types of subjects (cars, animals, people) and complex effects, ranging from semi-transparent elements such as smoke and reflections, to fully opaque effects such as objects attached to the subject.


title: Finite-dimensional representations of map superalgebras

abstract: We obtain a complete classification of all finite-dimensional irreducible modules over classical map superalgebras, provide formulas for their (super)characters and a description of their extension groups. Furthermore, we describe the block decomposition of the category of finite-dimensional modules for such map superalgebras. As an application, we specialize our results to the case of loop superalgebras in order to obtain a classification of finite-dimensional irreducible modules and block decomposition of the category of finite-dimensional modules over affine Lie superalgebras.


title: Improving exoplanet detection capabilities with the false inclusion probability. Comparison with other detection criteria in the context of radial velocities

abstract: Context. In exoplanet searches with radial velocity data, the most common statistical significance metrics are the Bayes factor and the false alarm probability (FAP). Both have proved useful, but do not directly address whether an exoplanet detection should be claimed. Furthermore, it is unclear which detection threshold should be taken and how robust the detections are to model misspecification. Aims. The present work aims at defining a detection criterion which conveys as precisely as possible the information needed to claim an exoplanet detection. We compare this new criterion to existing ones in terms of sensitivity and robustness. Methods. We define a significance metric called the false inclusion probability (FIP) based on the posterior probability of presence of a planet. Posterior distributions are computed with the nested sampling package Polychord. We show that for FIP and Bayes factor calculations, defining priors on linear parameters as Gaussian mixture models allows to significantly speed up computations. The performances of the FAP, Bayes factor and FIP are studied with simulations as well as analytical arguments. We compare the methods assuming the model is correct, then evaluate their sensitivity to the prior and likelihood choices. Results. Among other properties, the FIP offers ways to test the reliability of the significance levels, it is particularly efficient to account for aliasing and allows to exclude the presence of planets with a certain confidence. We find that, in our simulations, the FIP outperforms existing detection metrics. We show that planet detections are sensitive to priors on period and semi-amplitude and that letting free the noise parameters offers better performances than fixing a noise model based on a fit to ancillary indicators.

categories: astro-ph.EP astro-ph.IM stat.AP stat.CO stat.ME


title: Analytical Framework for Quantum Alternating Operator Ans\"atze

abstract: We develop a framework for analyzing layered quantum algorithms such as quantum alternating operator ans\"atze. Our framework relates quantum cost gradient operators, derived from the cost and mixing Hamiltonians, to classical cost difference functions that reflect cost function neighborhood structure. By considering QAOA circuits from the Heisenberg picture, we derive exact general expressions for expectation values as series expansions in the algorithm parameters, cost gradient operators, and cost difference functions. This enables novel interpretability and insight into QAOA behavior in various parameter regimes. For single-level QAOA1 we show the leading-order changes in the output probabilities and cost expectation value explicitly in terms of classical cost differences, for arbitrary cost functions. This demonstrates that, for sufficiently small positive parameters, probability flows from lower to higher cost states on average. By selecting signs of the parameters, we can control the direction of flow. We use these results to derive a classical random algorithm emulating QAOA1 in the small-parameter regime, i.e., that produces bitstring samples with the same probabilities as QAOA1 up to small error. For deeper QAOAp circuits we apply our framework to derive analogous and additional results in several settings. In particular we show QAOA always beats random guessing. We describe how our framework incorporates cost Hamiltonian locality for specific problem classes, including causal cone approaches, and applies to QAOA performance analysis with arbitrary parameters. We illuminate our results with a number of examples including applications to QUBO problems, MaxCut, and variants of MaxSat. We illustrate the application to QAOA circuits using mixing unitaries beyond the transverse-field mixer through two examples of constrained optimization, Max Independent Set and Graph Coloring.


title: Comment on "An axially symmetric spacetime with causality violation [Phys. Scr. 96 (07), 075208 (2021)], arXiv:2004.13769[gr-qc]"

abstract: Recently, an axially symmetric spacetime with causality violation is appeared in Ref. [1]. There, author presented a type III metric with vanishing expansion, shear and twist. The matter-energy represents a pure radial field with a negative cosmological constant. The spacetime is asymptotically anti-de Sitter space in the radial direction. Actually, this work with different for the metric was already published in Ref. [4]. Therefore, the current paper Ref. [1] is a duplicate one of the previous work.


title: A causal learning framework for the analysis and interpretation of COVID-19 clinical data

abstract: We present a workflow for clinical data analysis that relies on Bayesian Structure Learning (BSL), an unsupervised learning approach, robust to noise and biases, that allows to incorporate prior medical knowledge into the learning process and that provides explainable results in the form of a graph showing the causal connections among the analyzed features. The workflow consists in a multi-step approach that goes from identifying the main causes of patient's outcome through BSL, to the realization of a tool suitable for clinical practice, based on a Binary Decision Tree (BDT), to recognize patients at high-risk with information available already at hospital admission time. We evaluate our approach on a feature-rich COVID-19 dataset, showing that the proposed framework provides a schematic overview of the multi-factorial processes that jointly contribute to the outcome. We discuss how these computational findings are confirmed by current understanding of the COVID-19 pathogenesis. Further, our approach yields to a highly interpretable tool correctly predicting the outcome of 85% of subjects based exclusively on 3 features: age, a previous history of chronic obstructive pulmonary disease and the PaO2/FiO2 ratio at the time of arrival to the hospital. The inclusion of additional information from 4 routine blood tests (Creatinine, Glucose, pO2 and Sodium) increases predictive accuracy to 94.5%.


title: Actuarial strategy for pricing Asian options under a mixed fractional Brownian motion with jumps

abstract: The mixed fractional Brownian motion ($mfBm$) has become quite popular in finance, since it allows one to model long-range dependence and self-similarity while remaining, for certain values of the Hurst parameter, arbitrage-free. In the present paper, we propose approximate closed-form solutions for pricing arithmetic Asian options on an underlying described by the $mfBm$. Specifically, we consider both arithmetic Asian options and arithmetic Asian power options, and we obtain analytical formulas for pricing them based on a convenient approximation of the strike price. Both the standard $mfBm$ and the $mfBm$ with Poisson log-normally distributed jumps are taken into account.


title: Impacts of Water Latent Heat on the Thermal Structure of Ultra-Cool Objects: Brown Dwarfs and Free-Floating Planets

abstract: Brown dwarfs are essential targets for understanding planetary and sub-stellar atmospheres across a wide range of thermal and chemical conditions. As surveys continue to probe ever deeper, and as observing capabilities continue to improve, the number of known Y dwarfs -- the coldest class of sub-stellar objects, with effective temperatures below about 600 K -- is rapidly growing. Critically, this class of ultra-cool objects has atmospheric conditions that overlap with Solar System worlds and, as a result, tools and ideas developed from studying Earth, Jupiter, Saturn and other nearby worlds are well-suited for application to sub-stellar atmospheres. To that end, we developed a one-dimensional (vertical) atmospheric structure model for ultra-cool objects that includes moist adiabatic convection, as this is an important process for many Solar System planets. Application of this model across a range of effective temperatures (350, 300, 250, 200 K), metallicities ([M/H] of 0.0, 0.5, 0.7, 1.5), and gravities (log $g$ of 4.0, 4.5, 4.7, 5.0) demonstrates strong impacts of water latent heat release on simulated temperature-pressure profiles. At the highest metallicities, water vapor mixing ratios reach an Earth-like 3%, with associated major alterations to the thermal structure in the atmospheric regions where water condenses. Spectroscopic and photometric signatures of metallicity and moist convection should be readily detectable at near- and mid-infrared wavelengths, especially with James Webb Space Telescope observations, and can help indicate the formation history of an object.


title: Multiwavelength investigation of the candidate Galactic PeVatron MGRO J1908+06

abstract: The candidate PeVatron MGRO J1908+06, which shows a hard spectrum beyond 100 TeV, is one of the most peculiar $\gamma$-ray sources in the Galactic plane. Its complex morphology and some possible counterparts spatially related with the VHE emission region, preclude to distinguish between a hadronic or leptonic nature of the $\gamma$-ray emission. In this paper we illustrate a new multiwavelength analysis of MGRO J1908+06, with the aim to shed light on its nature and the origin of its ultra high-energy emission. We performed an analysis of the $^{12}$CO and $^{13}$CO molecular line emission demonstrating the presence of dense molecular clouds spatially correlated with the source region. We also analyzed 12-years of Fermi-LAT data between 10 GeV and 1 TeV finding a counterpart with a hard spectrum ($\Gamma \sim 1.6$). Our reanalysis of XMM-Newton data allowed us to put a more stringent constraint on the X-ray flux from this source. We demonstrate that a single accelerator cannot explain the whole set of multiwavelength data, regardless of whether it accelerates protons or electrons, but a 2-zone model is needed to explain the emission from MGRO J1908+06. The VHE emission seems most likely the superposition of a TeV PWN powered by PSR J1907+0602, in the southern part, and of the interaction between the supernova remnant G40.5-0.5 and the molecular clouds towards the northern region.


title: Real-time gravitational replicas: Low dimensional examples

abstract: We continue the study of real-time replica wormholes initiated in arXiv:2012.00828. Previously, we had discussed the general principles and had outlined a variational principle for obtaining stationary points of the real-time gravitational path integral. In the current work we present several explicit examples in low-dimensional gravitational theories where the dynamics is amenable to analytic computation. We demonstrate the computation of R\'enyi entropies in the cases of JT gravity and for holographic two-dimensional CFTs (using the dual gravitational dynamics). In particular, we explain how to obtain the large central charge result for subregions comprising of disjoint intervals directly from the real-time path integral.


title: Breaching the limit: formation of GW190521-like and IMBH mergers in young massive clusters

abstract: The LIGO-Virgo-Kagra collaboration (LVC) discovered recently GW190521, a gravitational wave (GW) source associated with the merger between two black holes (BHs) with mass $66$ M$_\odot$ and $>85$ M$_\odot$. GW190521 represents the first BH binary (BBH) merger with a primary mass falling in the "upper mass-gap" and the first leaving behind a $\sim 150$ M$_\odot$ remnant. So far, the LVC reported the discovery of four further mergers having a total mass $>100$ M$_\odot$, i.e. in the intermediate-mass black holes (IMBH) mass range. Here, we discuss results from a series of 80 $N$-body simulations of young massive clusters (YMCs) that implement relativistic corrections to follow compact object mergers. We discover the development of a GW190521-like system as the result of a 3rd-generation merger, and four IMBH-BH mergers with total mass $~(300-350)$ M$_\odot$. We show that these IMBH-BH mergers are low-frequency GW sources detectable with LISA and DECIGO out to redshift $z=0.01-0.1$ and $z>100$, and we discuss how their detection could help unravelling IMBH natal spins. For the GW190521 test case, we show that the 3rd-generation merger remnant has a spin and effective spin parameter that matches the $90\%$ credible interval measured for GW190521 better than a simpler double merger and comparably to a single merger. Due to GW recoil kicks, we show that retaining the products of these mergers require birth-sites with escape velocities $\gtrsim 50-100$ km s$^{-1}$, values typically attained in galactic nuclei and massive clusters with steep density profiles.


title: New Views of Black Holes from Computational Imaging

abstract: The unique challenges associated with imaging a black hole motivated the development of new computational imaging algorithms. As the Event Horizon Telescope continues to expand, these algorithms will need to evolve to keep pace with the increasingly demanding volume and dimensionality of the data.


title: Stellar Activity Manifesting at a One Year Alias Explains Barnard b as a False Positive

abstract: Barnard's star is among the most studied stars given its proximity to the Sun. It is often considered $the$ Radial Velocity (RV) standard for fully convective stars due to its RV stability and equatorial declination. Recently, an $M \sin i = 3.3 M_{\oplus}$ super-Earth planet candidate with a 233 day orbital period was announced by Ribas et al. (2018). New observations from the near-infrared Habitable-zone Planet Finder (HPF) Doppler spectrometer do not show this planetary signal. We ran a suite of experiments on both the original data and a combined original + HPF data set. These experiments include model comparisons, periodogram analyses, and sampling sensitivity, all of which show the signal at the proposed period of 233 days is transitory in nature. The power in the signal is largely contained within 211 RVs that were taken within a 1000 day span of observing. Our preferred model of the system is one which features stellar activity without a planet. We propose that the candidate planetary signal is an alias of the 145 day rotation period. This result highlights the challenge of analyzing long-term, quasi-periodic activity signals over multi-year and multi-instrument observing campaigns.


title: Optimal Virtual Network Embeddings for Tree Topologies

abstract: The performance of distributed and data-centric applications often critically depends on the interconnecting network. Applications are hence modeled as virtual networks, also accounting for resource demands on links. At the heart of provisioning such virtual networks lies the NP-hard Virtual Network Embedding Problem (VNEP): how to jointly map the virtual nodes and links onto a physical substrate network at minimum cost while obeying capacities. This paper studies the VNEP in the light of parameterized complexity. We focus on tree topology substrates, a case often encountered in practice and for which the VNEP remains NP-hard. We provide the first fixed-parameter algorithm for the VNEP with running time $O(3^r (s+r^2))$ for requests and substrates of $r$ and $s$ nodes, respectively. In a computational study our algorithm yields running time improvements in excess of 200x compared to state-of-the-art integer programming approaches. This makes it comparable in speed to the well-established ViNE heuristic while providing optimal solutions. We complement our algorithmic study with hardness results for the VNEP and related problems.


title: Pathfinder for a High Statistics Search for Missing Energy In Gamma Cascades

abstract: We investigate the feasibility of a high statistics experiment to search for invisible decay modes in nuclear gamma cascades using 200 kg of %36 Cs(Tl) scintillators that are presently available at Texas A\&M. The experiment aims to search for missing energy by robustly establishing the absence of a photon in a well identified gamma cascade. We report on the experimental demonstration of the energy resolution necessary for this search. Prior explorations of this detector concept focused on baryonically coupled physics that could be emitted in $E_2$ transitions. We point out that this protocol can also search for particles that are coupled to photons by searching for the conversion of a photon produced in a gamma cascade into a hidden particle. Examples of these processes include the oscillation of a photon into a hidden photon and the conversion of a photon into an axion-like-particle either in the presence of a magnetic field or via the Primakoff process. This proof-of-concept apparatus appears to have the ability to search for hitherto unconstrained baryonically coupled scalars and pseudoscalars produced in $E_0$ and $M_0$ transitions. If successfully implemented, this experiment serves as a pathfinder for a larger detector with greater containment that can thoroughly probe the existence of new particles with mass below 4 MeV that lie in the poorly constrained supernova ``trapping window'' that exists between 100 keV and 30 MeV.


title: Quantum Phases of Transition Metal Dichalcogenide Moir\'e Systems

abstract: Moir\'e systems provide a rich platform for studies of strong correlation physics. Recent experiments on hetero-bilayer transition metal dichalcogenide (TMD) Moir\'e systems are exciting in that they manifest a relatively simple model system of an extended Hubbard model on a triangular lattice. Inspired by the prospect of the hetero-TMD Moir\'e system's potential as a solid-state-based quantum simulator, we explore the extended Hubbard model on the triangular lattice using the density matrix renormalization group (DMRG). Specifically, we explore the two-dimensional phase space of the kinetic energy relative to the interaction strength $t/U$ and the further-range interaction strength $V_1/U$. We find competition between Fermi fluid, chiral spin liquid, spin density wave, and charge density wave. In particular, our finding of the optimal further-range interaction for the chiral correlation presents a tantalizing possibility.


title: Discrepancy between cosmological and electroweak observables in Higgs Inflation

abstract: In this work, we revisit the non-minimally coupled Higgs Inflation scenario and investigate its observational viability in light of the current Cosmic Microwave Background, Baryon Acoustic Oscillation and type Ia Supernovae data. We explore the effects of the Coleman-Weinberg approximation to the Higgs potential in the primordial universe, connecting the predictions for the Lagrangian parameters at inflationary scales to the electroweak observables through Renormalization Group methods at two-loop order. As the main result, we find that observations on the electroweak scale are in disagreement with the constraints obtained from the cosmological data sets used in the analysis. Specifically, an $\approx 8\sigma$-discrepancy between the inflationary parameters and the electroweak value of the top quark mass is found, which suggests that a significant deviation from the scenario analysed is required by the cosmological data.


title: Momentum space quantum Monte Carlo on twisted bilayer Graphene

abstract: We report an implementation of the momentum space quantum Monte Carlo (QMC) method on the interaction model for the twisted bilayer graphene (TBG) at integer fillings. The long-range Coulomb repulsion is treated exactly with the flat bands, spin and valley degrees of freedom of electrons taking into account. We prove the absence of the minus sign problem for QMC simulation at integer fillings when either the two valley or the two spin degrees of freedom are considered. By taking the realistic parameters of the twist angle and interlayer tunnelings into the simulation, we benchmark the QMC data with the exact band gap obtained at the chiral limit, to reveal the insulating ground states at the charge neutrality point (CNP). Then, with the exact Green's functions from QMC, we perform stochastic analytic continuation to obtain the first set of single-particle spectral function for the TBG model at CNP. Our momentum space QMC scheme therefore offers the controlled computation pathway for systematic investigation of the electronic states in realistic TBG model at various electron fillings.


title: Event rate predictions of strongly lensed gravitational waves with detector networks and more realistic templates

abstract: Strong lensing of gravitational waves (GWs) is attracting growing attention of the community. The event rates of lensed GWs by galaxies were predicted in numerous papers, which used some approximations to evaluate the GW strains detectable by a single detector. The joint-detection of GW signals by a network of instruments will increase the detecting ability of fainter and farther GW signals, which could increase the detection rate of the lensed GWs, especially for the 3rd generation detectors, e.g., Einstein Telescope (ET) and Cosmic Explorer (CE). Moreover, realistic GW templates will improve the accuracy of the prediction. In this work, we consider the detection of galaxy-scale lensed GW events under the 2nd, 2.5th, and 3rd generation detectors with the network scenarios and adopt the realistic templates to simulate GW signals. Our forecast is based on the Monte Carlo technique which enables us to take Earth's rotation into consideration. We find that the overall detection rate is improved, especially for the 3rd generation detector scenarios. More precisely, it increases by ~37% adopting realistic templates, and under network detection strategy, further increases by ~58% comparing with adoption of the realistic templates, and we estimate that the 3rd generation GW detectors will detect hundreds lensed events per year. The effect from the Earth's rotation is weakened in the detector network strategy.


title: Matrix-product state approach to the generalized nuclear pairing Hamiltonian

abstract: We show that from the point of view of the generalized pairing Hamiltonian, the atomic nucleus is a system with small entanglement and can thus be described efficiently using a 1D tensor network (matrix-product state) despite the presence of long-range interactions. The ground state can be obtained using the density-matrix renormalization group (DMRG) algorithm, which is accurate up to machine precision even for large nuclei, is numerically as cheap as the widely used BCS (Bardeen-Cooper-Schrieffer) approach, and does not suffer from any mean-field artifacts. We apply this framework to compute the even-odd mass differences of all known lead isotopes from $^{178}$Pb to $^{220}$Pb in a very large configuration space of 13 shells between the neutron magic numbers 82 and 184 (i.e., two major shells) and find good agreement with the experiment. We also treat pairing with non-zero angular momentum and determine the lowest excited states in the full configuration space of one major shell, which we demonstrate for the $N=126$, $Z\geq 82$ isotones. To demonstrate the capabilities of the method beyond low-lying excitations, we calculate the first 100 excited states of $^{208}$Pb with singlet pairing and the two-neutron removal spectral function of $^{210}$Pb, which relates to a two-neutron pickup experiment.


title: Interpolating Detailed Simulations of Kilonovae: Adaptive Learning and Parameter Inference Applications

abstract: Detailed radiative transfer simulations of kilonovae are difficult to apply directly to observations; they only sparsely cover simulation parameters, such as the mass, velocity, morphology, and composition of the ejecta. On the other hand, semianalytic models for kilonovae can be evaluated continuously over model parameters, but neglect important physical details which are not incorporated in the simulations, thus introducing systematic bias. Starting with a grid of 2D anisotropic simulations of kilonova light curves covering a wide range of ejecta properties, we apply adaptive-learning techniques to iteratively choose new simulations and produce high-fidelity surrogate models for those simulations. These surrogate models allow for continuous evaluation across model parameters while retaining the microphysical details about the ejecta. Using a new code for multimessenger inference, we demonstrate how to use our interpolated models to infer kilonova parameters. Comparing to inferences using simplified analytic models, we recover different ejecta properties. We discuss the implications of this analysis which is qualitatively consistent with similar previous work using detailed ejecta opacity calculations and which illustrates systematic challenges for kilonova modeling. An associated data and code release provides our interpolated light-curve models, interpolation implementation which can be applied to reproduce our work or extend to new models, and our multimessenger parameter inference engine.


title: SMURF: Self-Teaching Multi-Frame Unsupervised RAFT with Full-Image Warping

abstract: We present SMURF, a method for unsupervised learning of optical flow that improves state of the art on all benchmarks by $36\%$ to $40\%$ (over the prior best method UFlow) and even outperforms several supervised approaches such as PWC-Net and FlowNet2. Our method integrates architecture improvements from supervised optical flow, i.e. the RAFT model, with new ideas for unsupervised learning that include a sequence-aware self-supervision loss, a technique for handling out-of-frame motion, and an approach for learning effectively from multi-frame video data while still only requiring two frames for inference.


title: Kostka semigroups and generalized Dyck paths

abstract: We prove a conjecture of S. Gao-J. Kiers-G. Orelowitz-A. Yong which asserts the reducibility of certain generalized Dyck paths. This gives a strengthening, and new proof, for their Width Bound Theorem on the Hilbert basis of the Kostka semigroup.


title: Invariant manifolds in Hamiltonian systems with applications to the Earth-Moon system

abstract: Invariant manifolds are the skeleton of the chaotic dynamics in Hamiltonian systems. In Celestial Mechanics, for instance, these geometrical structures are applied to a multitude of physical and practical problems, such as to the description of the natural transport of asteroids and to the construction of trajectories for artificial satellites. In this work, we use efficient numerical methods to visually illustrate the influence of invariant manifolds, which are associated with specific equilibrium points and unstable periodic orbits, in the dynamical properties of Hamiltonian systems. First, we investigate an area-preserving version of the two-dimensional H\'enon map. Later, we focus our investigation on the motion of a body with negligible mass that moves due to the gravitational attraction of both the Earth and the Moon. As a model, we adopt the planar circular restricted three-body problem, a near-integrable Hamiltonian system with two degrees of freedom. For both systems, we show how the selected invariant manifolds and the structure of the phase space evolve alongside each other. Our results contribute to the understanding of the connection between dynamics and geometry in Hamiltonian systems.


title: Efficient Quasi-Geodesics on the Stiefel Manifold

abstract: Solving the so-called geodesic endpoint problem, i.e., finding a geodesic that connects two given points on a manifold, is at the basis of virtually all data processing operations, including averaging, clustering, interpolation and optimization. On the Stiefel manifold of orthonormal frames, this problem is computationally involved. A remedy is to use quasi-geodesics as a replacement for the Riemannian geodesics. Quasi-geodesics feature constant speed and covariant acceleration with constant (but possibly non-zero) norm. For a well-known type of quasi-geodesics, we derive a new representation that is suited for large-scale computations. Moreover, we introduce a new kind of quasi-geodesics that turns out to be much closer to the Riemannian geodesics.


title: Simplified Hartree-Fock Computations on Second-Row Atoms

abstract: Simplified Hartree-Fock computations are carried out on the atoms He through Ne, using orthonormalized basis functions for the 1s, 2s and 2p orbitals dependent on three parameters. Using Mathematica with the new Apple M1 chip, computations require about 0.005 seconds of CPU time. Approximate energies within 1% of the best H-F values are thereby obtained, with an order of magnitude less computational effort.


title: Chord Recognition- Music and Audio Information Retrieval

abstract: Music Information Retrieval (MIR) is a collaborative scientific study that help to build innovative information research themes, novel frameworks, and developing connected delivery mechanisms in addition to making the world's massive collection of music open for everyone. Modern rock music proved to be difficult to estimate tempo and chord recognition did not work. All of the findings indicate that modern rock and metal music can be analysed, despite its complexity, but that further research is needed in this area to make it useful. Using a neural network has been one of the simplest ways of dealing with it. The pitch class profile vector is used in the neural network method. Because the vector only contains 12 elements of semi-tone values, it is enough for chord recognition. Of course, there are other ways of achieving this work, most of them depend on pitch class profiling to transform the chord into a type that can be recognised, but the recognition process is time-consuming centred on extremely complicated and memory-intensive methods.