title: The Massive and Distant Clusters of WISE Survey IX: High Radio Activity in a Merging Cluster

abstract: We present a multi-wavelength investigation of the radio galaxy population in the galaxy cluster MOO J1506+5137 at $z$=1.09$\pm$0.03, which in previous work we identified as having multiple complex radio sources. The combined dataset used in this work includes data from the Low-Frequency Array Two-metre Sky Survey (LoTSS), NSF's Karl G. Jansky Very Large Array (VLA), the Robert C. Byrd Green Bank Telescope (GBT), the Spitzer Space Telescope, and the Dark Energy Camera Legacy Survey (DECaLS). We find that there are five radio sources which are all located within 500 kpc ($\sim$1$^{\prime}$) of the cluster center and have radio luminosities $P_{\mathrm{1.4GHz}}$ > 1.6$\times$10$^{24}$ W Hz$^{-1}$. The typical host galaxies are among the highest stellar mass galaxies in the cluster. The exceptional radio activity among the massive galaxy population appears to be linked to the dynamical state of the cluster. The galaxy distribution suggests an ongoing merger, with a subgroup found to the northwest of the main cluster. Further, two of the five sources are classified as bent-tail sources with one being a potential wide-angle tail (WAT)/hybrid morphology radio source (HyMoRS) indicating a dynamic environment. The cluster also lies in a region of the mass-richness plane occupied by other merging clusters in the Massive and Distant Clusters of WISE Survey (MaDCoWS). The data suggest that during the merger phase radio activity can be dramatically enhanced, which would contribute to the observed trend of increased radio activity in clusters with increasing redshift.


title: Reproducing sub-millimetre galaxy number counts with cosmological hydrodynamic simulations

abstract: Matching the number counts of high-$z$ sub-millimetre-selected galaxies (SMGs) has been a long standing problem for galaxy formation models. In this paper, we use 3D dust radiative transfer to model the sub-mm emission from galaxies in the SIMBA cosmological hydrodynamic simulations, and compare predictions to the latest single-dish observational constraints on the abundance of 850$\mathrm{\mu m}$-selected sources. We find good agreement with the shape of the integrated 850$\mathrm{\mu m}$ luminosity function, and the normalisation is within 0.25 dex at $> 3 \; \mathrm{mJy}$, unprecedented for a fully cosmological hydrodynamic simulation, along with good agreement in the redshift distribution of bright SMGs. The agreement is driven primarily by SIMBA's good match to infrared measures of the star formation rate (SFR) function between $z = 2-4$ at high SFRs. Also important is the self-consistent on-the-fly dust model in SIMBA, which predicts, on average, higher dust masses (by up to a factor of 2.5) compared to using a fixed dust-to-metals ratio of 0.3. We construct a lightcone to investigate the effect of far-field blending, and find that 52% of sources are blends of multiple components, which makes a small contribution to the normalisation of the bright-end of the number counts. We provide new fits to the 850$\mathrm{\mu m}$ luminosity as a function of SFR and dust mass. Our results demonstrate that exotic solutions to the discrepancy between sub-mm counts in simulations and observations, such as a top-heavy IMF, are unnecessary, and that sub-millimetre-bright phases are a natural consequence of massive galaxy evolution.


title: Improved $(g-2)_\mu$ Measurements and Supersymmetry

abstract: The electroweak (EW) sector of the Minimal Supersymmetric Standard Model (MSSM) can account for a variety of experimental data. The lighest supersymmetric particle (LSP), which we take as the lightest neutralino, $\tilde \chi_1^0$, can account for the observed Dark Matter (DM) content of the universe via coannihilation with the next-to-LSP (NLSP), while being in agreement with negative results from Direct Detection (DD) experiments. Owing to relatively small production cross-sections a comparably light EW sector of the MSSM is also in agreement with the unsuccessful searches at the LHC. Most importantly, the EW sector of the MSSM can account for the persistent $3-4\,\sigma$ discrepancy between the experimental result for the anomalous magnetic moment of the muon, $(g-2)_\mu$, and its Standard Model (SM) prediction. Under the assumption that the $\tilde \chi_1^0$ provides the full DM relic abundance we first analyze which mass ranges of neutralinos, charginos and scalar leptons are in agreement with all experimental data, including relevant LHC searches. We find an upper limit of $\sim 600$ GeV for the LSP and NLSP masses. In a second step we assume that the new result of the Run 1 of the ``MUON G-2'' collaboration at Fermilab yields a precision comparable to the existing experimental result with the same central value. We analyze the potential impact of the combination of the Run 1 data with the existing $(g-2)_\mu$ data on the allowed MSSM parameter space. We find that in this case the upper limits on the LSP and NLSP masses are substantially reduced by roughly $100$ GeV. This would yield improved upper limits on these masses of $\sim 500$ GeV. In this way, a clear target could be set for future LHC EW searches, as well as for future high-energy $e^+e^-$ colliders, such as the ILC or CLIC.


title: Relative Arbitrage Opportunities in $N$ Investors and Mean-Field Regimes

abstract: The relative arbitrage portfolio, formulated in Stochastic Portfolio Theory (SPT), outperforms a benchmark portfolio over a given time-horizon with probability one. This paper analyzes the market behavior and optimal investment strategies to attain relative arbitrage both in the $N$ investors and mean field regimes under some market conditions. An investor competes with a benchmark of market and peer investors, expecting to outperform the benchmark and minimizing the initial capital. With market price of risk processes depending on the market portfolio and investors, we develop a systematic way to solve multi-agent optimization problem within SPT's framework. The objective can be characterized by the smallest nonnegative continuous solution of a Cauchy problem. By a modification in the structure of the extended mean field game with common noise and its notion of the uniqueness of Nash equilibrium, we show a unique equilibrium in $N$-player games and mean field games with mild conditions on the equity market.


title: Informing dark matter direct detection limits with the ARTEMIS simulations

abstract: Dark matter (DM) direct detection experiments aim to place constraints on the DM--nucleon scattering cross-section and the DM particle mass. These constraints depend sensitively on the assumed local DM density and velocity distribution function. While astrophysical observations can inform the former (in a model-dependent way), the latter is not directly accessible with observations. Here we use the high-resolution ARTEMIS cosmological hydrodynamical simulation suite of 42 Milky Way-mass halos to explore the spatial and kinematical distributions of the DM in the solar neighbourhood, and we examine how these quantities are influenced by substructures, baryons, the presence of dark discs, as well as general halo-to-halo scatter (cosmic variance). We also explore the accuracy of the standard Maxwellian approach for modelling the velocity distribution function. We find significant halo-to-halo scatter in the density and velocity functions which, if propagated through the standard halo model for predicting the DM detection limits, implies a significant scatter about the typically quoted limit. We also show that, in general, the Maxwellian approximation works relatively well for simulations that include the important gravitational effects of baryons, but is less accurate for collisionless (DM-only) simulations. Given the significant halo-to-halo scatter in quantities relevant for DM direct detection, we advocate propagating this source of uncertainty through in order to derive conservative DM detection limits.


title: Dissecting power of intersection of two context-free languages

abstract: We say that a language $L$ is \emph{constantly growing} if there is a constant $c$ such that for every word $u\in L$ there is a word $v\in L$ with $\vert u\vert<\vert v\vert\leq c+\vert u\vert$. We say that a language $L$ is \emph{geometrically growing} if there is a constant $c$ such that for every word $u\in L$ there is a word $v\in L$ with $\vert u\vert<\vert v\vert\leq c\vert u\vert$. Given two infinite languages $L_1,L_2$, we say that $L_1$ \emph{dissects} $L_2$ if $\vert L_2\setminus L_1\vert=\infty$ and $\vert L_1\cap L_2\vert=\infty$. In 2013, it was shown that for every constantly growing language $L$ there is a regular language $R$ such that $R$ dissects $L$. In the current article we show how to dissect a geometrically growing language by a homomorphic image of intersection of two context-free languages. Consider three alphabets $\Gamma$, $\Sigma$, and $\Theta$ such that $\vert \Sigma\vert=1$ and $\vert \Theta\vert=4$. We prove that there are context-free languages $M_1,M_2\subseteq \Theta^*$, an erasing alphabetical homomorphism $\pi:\Theta^*\rightarrow \Sigma^*$, and a nonerasing alphabetical homomorphism $\varphi : \Gamma^*\rightarrow \Sigma^*$ such that: If $L\subseteq \Gamma^*$ is a geometrically growing language then there is a regular language $R\subseteq \Theta^*$ such that $\varphi^{-1}\left(\pi\left(R\cap M_1\cap M_2\right)\right)$ dissects the language $L$.


title: Atomic scale model and electronic structure of Cu$_2$O/CH$_3$NH$_3$PbI$_3$ interfaces in perovskite solar cells

abstract: Cuprous oxide has been conceived as a potential alternative to traditional organic hole transport layers in hybrid halide perovskite-based solar cells. Device simulations predict record efficiencies using this semiconductor, but experimental results do not yet show this trend. More detailed knowledge about the Cu$_2$O/perovskite interface is mandatory to improve the photoconversion efficiency. Using density functional theory calculations, here we study the interfaces of CH$_3$NH$_3$PbI$_3$ with Cu$_2$O to assess their influence on device performance. Several atomistic models of these interfaces are provided for the first time, considering different compositions of the interface atomic planes. The interface electronic properties are discussed on the basis of the optimal theoretical situation, but in connection with the experimental realizations and device simulations. It is shown that the formation of vacancies in the Cu$_2$O terminating planes is essential to eliminate dangling bonds and trap states. The four interface models that fulfill this condition present a band alignment favorable for photovoltaic conversion. Energy of adhesion, and charge transfer across the interfaces are also studied. The termination of CH$_3$NH$_3$PbI$_3$ in PbI$_2$ atomic planes seems optimal to maximize the photoconversion efficiency.


title: C19-TraNet: an empirical, global index-case transmission network of SARS-CoV-2

abstract: Originating in Wuhan, the novel coronavirus, severe acute respiratory syndrome 2 (SARS-CoV-2), has astonished health-care systems across globe due to its rapid and simultaneous spread to the neighboring and distantly located countries. To gain the systems level understanding of the role of global transmission routes in the COVID-19 spread, in this study, we have developed the first, empirical, global, index-case transmission network of SARS-CoV-2 termed as C19-TraNet. We manually curated the travel history of country wise index-cases using government press releases, their official social media handles and online news reports to construct this C19-TraNet that is a spatio-temporal, sparse, growing network comprising of 187 nodes and 199 edges and follows a power-law degree distribution. To model the growing C19-TraNet, a novel stochastic scale free (SSF) algorithm is proposed that accounts for stochastic addition of both nodes as well as edges at each time step. A peculiar connectivity pattern in C19-TraNet is observed, characterized by a fourth degree polynomial growth curve, that significantly diverges from the average random connectivity pattern obtained from an ensemble of its 1,000 SSF realizations. Partitioning the C19-TraNet, using edge betweenness, it is found that most of the large communities are comprised of a heterogeneous mixture of countries belonging to different world regions suggesting that there are no spatial constraints on the spread of disease. This work characterizes the superspreaders that have very quickly transported the virus, through multiple transmission routes, to long range geographical locations alongwith their local neighborhoods.


title: Monotone Normality and Nabla-Products

abstract: Roitman's combinatorial principle $\Delta$ is equivalent to monotone normality of the nabla product, $\nabla (\omega +1)^\omega$. If $\{ X_n : n\in \omega\}$ is a family of metrizable spaces and $\nabla_n X_n$ is monotonically normal, then $\nabla_n X_n$ is hereditarily paracompact. Hence, if $\Delta$ holds then the box product $\square (\omega +1)^\omega$ is paracompact. Large fragments of $\Delta$ hold in $\mathsf{ZFC}$, yielding large subspaces of $\nabla (\omega+1)^\omega$ that are `really' monotonically normal. Countable nabla products of metrizable spaces which are respectively: arbitrary, of size $\le \mathfrak{c}$, or separable, are monotonically normal under respectively: $\mathfrak{b}=\mathfrak{d}$, $\mathfrak{d}=\mathfrak{c}$ or the Model Hypothesis. It is consistent and independent that $\nabla A(\omega_1)^\omega$ and $\nabla (\omega_1+1)^\omega$ are hereditarily normal (or hereditarily paracompact, or monotonically normal). In $\mathsf{ZFC}$ neither $\nabla A(\omega_2)^\omega$ nor $\nabla (\omega_2+1)^\omega$ is hereditarily normal.


title: The Carnegie Supernova Project-I: Correlation Between Type Ia Supernovae and Their Host Galaxies from Optical to Near-Infrared Bands

abstract: We present optical and near-infrared ($ugriYJH$) photometry of host galaxies of Type Ia supernovae (SN~Ia) observed by the \textit{Carnegie Supernova Project-I}. We determine host galaxy stellar masses and, for the first time, study their correlation with SN~Ia standardized luminosity across optical and near-infrared ($uBgVriYJH$) bands. In the individual bands, we find that SNe~Ia are more luminous in more massive hosts with luminosity offsets ranging between $-0.07 \pm0.03$ mag to $-0.15\pm0.04$ mag after light-curve standardization. The slope of the SN~Ia Hubble residual-host mass relation is negative across all $uBgVriYJH$ bands with values ranging between $-0.036\pm 0.025$ mag/dex to $-0.097\pm 0.027$ mag/dex -- implying that SNe~Ia in more massive galaxies are brighter than expected. The near-constant observed correlations across optical and near-infrared bands indicate that dust may not play a significant role in the observed luminosity offset--host mass correlation. We measure projected separations between SNe~Ia and their host centers, and find that SNe~Ia that explode beyond a projected 10 kpc have a $\rm 30\% \ to \ 50\%$ reduction of the dispersion in Hubble residuals across all bands -- making them a more uniform subset of SNe~Ia. Dust in host galaxies, peculiar velocities of nearby SN~Ia, or a combination of both may drive this result as the color excesses of SNe~Ia beyond 10 kpc are found to be generally lower than those interior, but there is also a diminishing trend of the dispersion as we exclude nearby events. We do not find that SN~Ia average luminosity varies significantly when they are grouped in various host morphological types. Host galaxy data from this work will be useful, in conjunction with future high-redshift samples, in constraining cosmological parameters.


title: Forecasting Precipitable Water Vapor Using LSTMs

abstract: Long-Short-Term-Memory (LSTM) networks have been used extensively for time series forecasting in recent years due to their ability of learning patterns over different periods of time. In this paper, this ability is applied to learning the pattern of Global Positioning System (GPS)-based Precipitable Water Vapor (PWV) measurements over a period of 4 hours. The trained model was evaluated on more than 1500 hours of recorded data. It achieves a root mean square error (RMSE) of 0.098 mm for a forecasting interval of 5 minutes in the future, and outperforms the naive approach for a lead-time of up to 40 minutes.


title: Dominate or Delete: Decentralized Competing Bandits in Serial Dictatorship

abstract: Online learning in a two-sided matching market, with demand side agents continuously competing to be matched with supply side (arms), abstracts the complex interactions under partial information on matching platforms (e.g. UpWork, TaskRabbit). We study the decentralized serial dictatorship setting, a two-sided matching market where the demand side agents have unknown and heterogeneous valuation over the supply side (arms), while the arms have known uniform preference over the demand side (agents). We design the first decentralized algorithm -- UCB with Decentralized Dominant-arm Deletion (UCB-D3), for the agents, that does not require any knowledge of reward gaps or time horizon. UCB-D3 works in phases, where in each phase, agents delete \emph{dominated arms} -- the arms preferred by higher ranked agents, and play only from the non-dominated arms according to the UCB. At the end of the phase, agents broadcast in a decentralized fashion, their estimated preferred arms through {\em pure exploitation}. We prove both, a new regret lower bound for the decentralized serial dictatorship model, and that UCB-D3 is order optimal.


title: Deep Involutive Generative Models for Neural MCMC

abstract: We introduce deep involutive generative models, a new architecture for deep generative modeling, and use them to define Involutive Neural MCMC, a new approach to fast neural MCMC. An involutive generative model represents a probability kernel $G(\phi \mapsto \phi')$ as an involutive (i.e., self-inverting) deterministic function $f(\phi, \pi)$ on an enlarged state space containing auxiliary variables $\pi$. We show how to make these models volume preserving, and how to use deep volume-preserving involutive generative models to make valid Metropolis-Hastings updates based on an auxiliary variable scheme with an easy-to-calculate acceptance ratio. We prove that deep involutive generative models and their volume-preserving special case are universal approximators for probability kernels. This result implies that with enough network capacity and training time, they can be used to learn arbitrarily complex MCMC updates. We define a loss function and optimization algorithm for training parameters given simulated data. We also provide initial experiments showing that Involutive Neural MCMC can efficiently explore multi-modal distributions that are intractable for Hybrid Monte Carlo, and can converge faster than A-NICE-MC, a recently introduced neural MCMC technique.


title: Train and You'll Miss It: Interactive Model Iteration with Weak Supervision and Pre-Trained Embeddings

abstract: Our goal is to enable machine learning systems to be trained interactively. This requires models that perform well and train quickly, without large amounts of hand-labeled data. We take a step forward in this direction by borrowing from weak supervision (WS), wherein models can be trained with noisy sources of signal instead of hand-labeled data. But WS relies on training downstream deep networks to extrapolate to unseen data points, which can take hours or days. Pre-trained embeddings can remove this requirement. We do not use the embeddings as features as in transfer learning (TL), which requires fine-tuning for high performance, but instead use them to define a distance function on the data and extend WS source votes to nearby points. Theoretically, we provide a series of results studying how performance scales with changes in source coverage, source accuracy, and the Lipschitzness of label distributions in the embedding space, and compare this rate to standard WS without extension and TL without fine-tuning. On six benchmark NLP and video tasks, our method outperforms WS without extension by 4.1 points, TL without fine-tuning by 12.8 points, and traditionally-supervised deep networks by 13.1 points, and comes within 0.7 points of state-of-the-art weakly-supervised deep networks-all while training in less than half a second.


title: Frequency Shifts of photons emitted from geodesics of nonlinear electromagnetic black holes

abstract: We analyzed the frequency shifts of photons emitted from massive particles that move in stable geodesics around black holes with non-linear electrodynamics. The motion of the photon is modified by an effective metric when the Einstein gravity is coupled with non-linear electrodynamics. The kinematic shifts for photons are modified by the magnetic and electric factors that make the difference between the linear and non-linear electrodynamics. As an illustration, we present the frequency shifts of the photons for black holes with non-linear electrodynamics, two magnetic black holes, and two electric black holes, the Bardeen and Bronnikov black holes for the magnetic charge and the Born-Infeld and Dymnikova black holes for electric charge. These frequency shifts are compared with their linear electromagnetic counterpart.


title: Dynamical evolution of the inner asteroid belt

abstract: A determination of the dynamical evolution of the asteroid belt is difficult because the asteroid belt has evolved since the time of asteroid formation through mechanisms that include: (1) catastrophic collisions, (2) rotational disruption, (3) chaotic orbital evolution and (4) orbital evolution driven by Yarkovsky radiation forces. The timescales of these loss mechanisms are uncertain and there is a need for more observational constraints. In the inner main belt, the mean size of the non-family asteroids increases with increasing inclination. Here, we use that observation to show that all inner main belt asteroids originate from either the known families or from ghost families, that is, old families with dispersed orbital elements. We estimate that the average age of the asteroids in the ghost families is a factor of 1/3 less than the Yarkovsky orbital evolution timescale. However, this orbital evolution timescale is a long-term average that must allow for the collisional evolution of the asteroids and for stochastic changes in their spin directions. By applying these constraints on the orbital evolution timescales to the evolution of the size-frequency distribution of the Vesta asteroid family, we estimate that the age of this family is greater than 1.3 $Gyr$ and could be comparable with the age of the solar system. By estimating the number of ghost families, we calculate that the number of asteroids that are the root sources of the meteorites and the near-Earth asteroids that originate from the inner main belt is about 20.


title: Internal surface plasmon excitation as the root cause of laser-induced periodic plasma structure and self-organized nanograting formation in the volume of transparent dielectric

abstract: A computer simulation of dynamics of an optical discharge produced in the volume of a transparent dielectric (fused silica) by a focused femtosecond laser pulse was carried out taking into account the possibility of developing small-scale ionization-field instability. The presence of small foreign inclusions in the fused silica was taken into account with the model of a nanodispersed heterogeneous medium by using Maxwell Garnett formulas. The results of the calculations made it possible to reveal the previously unknown physical mechanism that determines the periodicity of the ordered plasma-field structure that is formed in each single breakdown pulse and is the root cause of the ordered volume nanograting formation in dielectric material exposed to a series of repeated pulse. Two main points are decisive in this mechanism: (i) the formation of a thin overcritical plasma layer at the breakdown wave front counter-propagated to the incident laser pulse and (ii) the excitation of the "internal surface plasmon" at this front, resulting in a rapid amplification of the corresponding spatial harmonic of random seed perturbations in the plasma and formation of a contrast structure with a period equal to the wavelength of the surface plasmon (0.7 of the wavelength in dielectric).


title: A Comparison of Uncertainty Estimation Approaches in Deep Learning Components for Autonomous Vehicle Applications

abstract: A key factor for ensuring safety in Autonomous Vehicles (AVs) is to avoid any abnormal behaviors under undesirable and unpredicted circumstances. As AVs increasingly rely on Deep Neural Networks (DNNs) to perform safety-critical tasks, different methods for uncertainty quantification have recently been proposed to measure the inevitable source of errors in data and models. However, uncertainty quantification in DNNs is still a challenging task. These methods require a higher computational load, a higher memory footprint, and introduce extra latency, which can be prohibitive in safety-critical applications. In this paper, we provide a brief and comparative survey of methods for uncertainty quantification in DNNs along with existing metrics to evaluate uncertainty predictions. We are particularly interested in understanding the advantages and downsides of each method for specific AV tasks and types of uncertainty sources.


title: Orlov's Theorem in the Smooth Proper Case

abstract: We extend Orlov's result that certain functors between derived categories of smooth projective varieties are Fourier--Mukai transforms to the case when those varieties are smooth and proper.


title: Beyond Quality and Quantity: Contact Distribution Encodes Frictional Strength

abstract: Classically, the quantity of contact area $A_R$ between two bodies is considered a proxy for the force of friction. However, bond density across the interface - quality of contact - is also relevant, and contemporary debate often centers around the relative importance of these two factors. In this work, we demonstrate that a third factor, often overlooked, plays a significant role in static frictional strength: the distribution of contact. We perform static friction measurements, $\mu$, on three pairs of solid blocks while imaging the contact plane. By using linear regression on hundreds of image-$\mu$ pairs, we are able to predict future friction measurements with 3 to 7 times better accuracy than existing benchmarks, including total quantity of contact area. Our model has no access to quality of contact, and we therefore conclude that a large portion of the interfacial state is encoded in the spatial distribution of contact, rather than its quality or quantity.


title: Application of Neuroevolution in Autonomous Cars

abstract: With the onset of Electric vehicles, and them becoming more and more popular, autonomous cars are the future in the travel/driving experience. The barrier to reaching level 5 autonomy is the difficulty in the collection of data that incorporates good driving habits and the lack thereof. The problem with current implementations of self-driving cars is the need for massively large datasets and the need to evaluate the driving in the dataset. We propose a system that requires no data for its training. An evolutionary model would have the capability to optimize itself towards the fitness function. We have implemented Neuroevolution, a form of genetic algorithm, to train/evolve self-driving cars in a simulated virtual environment with the help of Unreal Engine 4, which utilizes Nvidia's PhysX Physics Engine to portray real-world vehicle dynamics accurately. We were able to observe the serendipitous nature of evolution and have exploited it to reach our optimal solution. We also demonstrate the ease in generalizing attributes brought about by genetic algorithms and how they may be used as a boilerplate upon which other machine learning techniques may be used to improve the overall driving experience.


title: Bookworm continual learning: beyond zero-shot learning and continual learning

abstract: We propose bookworm continual learning(BCL), a flexible setting where unseen classes can be inferred via a semantic model, and the visual model can be updated continually. Thus BCL generalizes both continual learning (CL) and zero-shot learning (ZSL). We also propose the bidirectional imagination (BImag) framework to address BCL where features of both past and future classes are generated. We observe that conditioning the feature generator on attributes can actually harm the continual learning ability, and propose two variants (joint class-attribute conditioning and asymmetric generation) to alleviate this problem.


title: Initial fluctuations and power spectrum of flow anisotropies in relativistic heavy-ion collisions

abstract: Flow has emerged as a crucial probe for the properties of the thermalized medium produced in relativistic heavy-ion collisions. The evolution of initial state fluctuations leaves imprints on the power spectrum of flow coefficients. Therefore flow coefficients are a crucial probe of initial state fluctuations arising from the parton distributions of the colliding nuclei. This has a very strong correspondence with the physics of power spectrum of cosmic microwave background radiation (CMBR) anisotropies which directly probes initial inflationary fluctuations. Much work has been done to probe these interesting interconnections, in particular, in developing techniques for the measurements of higher flow coefficients. We present a short review of these developments. The effect of initial magnetic field on these features will also be reviewed. All this acquires special importance in view of upcoming electron-ion collider which will directly probe initial parton distribution of the colliding nucleus.


title: State Complexity of Permutation and Related Decision Problems on Alphabetical Pattern Constraints

abstract: We investigate the state complexity of the permutation operation, or the commutative closure, on Alphabetical Pattern Constraints (APC). This class corresponds to level $3/2$ of the Straubing-Th{\'e}rien Hierarchy and includes the finite, the piecewise-testable, or $\mathcal J$-trivial, and the $\mathcal R$-trivial and $\mathcal L$-trivial languages. We give a sharp state complexity bound expressed in terms of the longest strings in the unary projection languages of an associated finite language and which is already sharp for the subclass of finite languages. Additionally, for two subclasses, we give sharp bounds expressed in terms of the size of a recognizing input automaton and the size of the alphabet. Lastly, we investigate the inclusion and universality problem on APCs up to permutational equivalence, two problems known to be PSPACE-complete on APCs even for fixed alphabets in general, and show them to be decidable in polynomial time for fixed alphabets in this case.


title: An Analysis of Academic Performance of University Students in Namibia

abstract: Based on observations it seems that a considerable proportion of university students in Namibia need to enhance their academic performance in regard to be fully competitive in a globalised working environment. This specifically includes thorough understanding of mathematics, as it forms an integral part of modern knowledge disciplines such as information technology and the management sciences. Yet students struggle or fail to engage in relevant coursework, either because of an absence of adequate material, capable educators, their own will power, or a combination thereof. This study aims to investigate the critical factors that are related to academic performance of university students in Namibia.


title: Characteristic polynomials of products of Wigner matrices: finite-N results and Lyapunov universality

abstract: We compute the average characteristic polynomial of the hermitised product of $M$ real or complex Wigner matrices of size $N\times N$ and the average of the characteristic polynomial of a product of $M$ such Wigner matrices times the characteristic polynomial of the conjugate matrix. Surprisingly, the results agree with that of the product of $M$ real or complex Ginibre matrices at finite-$N$, which have i.i.d. Gaussian entries. For the latter the average characteristic polynomial yields the orthogonal polynomial for the singular values of the product matrix, whereas the product of the two characteristic polynomials involves the kernel of complex eigenvalues. This extends the result of Forrester and Gamburd for one characteristic polynomial of a single random matrix and only depends on the first two moments. In the limit $M\to\infty$ at fixed $N$ we determine the locations of the zeros of a single characteristic polynomial, rescaled as Lyapunov exponents by taking the logarithm of the $M$th root. The position of the $j$th zero agrees asymptotically for large-$j$ with the position of the $j$th Lyapunov exponent for products of Gaussian random matrices, hinting at the universality of the latter.


title: Designing Laboratories for Online Instruction using the iOLab Device

abstract: Scientific laboratories are among the most challenging course components to integrate into online instruction. Available technology restricts the design and nature of experiments and it can be hard to replicate the collaborative lab environment where frequent and immediate instructor feedback is the norm. Here we report on technological and pedagogical aspects of newly developed labs for online courses using the Interactive Online Lab (iOLab) device. We argue that this technology, coupled with an online course design emphasizing teamwork, targeted feedbacks, and self-regulation skills, provides a robust framework for students to do reliable, engaging, inquiry-based and hands-on labs outside the classroom. After describing the implementation and technology, we explain our lab objectives and how the labs were integrated into two introductory physics courses. We conclude with an example lab on kinematics.


title: Dynamic Constraint-based Influence Framework and its Application in Stochastic Modeling of Load Balancing

abstract: Components connected over a network influence each other and interact in various ways. Examples of such systems are networks of computing nodes, which the nodes interact by exchanging workload, for instance, for load balancing purposes. In this paper, we first study the Influence Model, a networked Markov chain framework, for modeling network interactions and discuss two key limitations of this model, which cause it to fall short in modeling constraint-based and dynamic interactions in networks. Next, we propose the Dynamic and Constraint-based Influence Model (DCIM) to alleviate the limitations. The DCIM extends the application of the Influence Model to more general network interaction scenarios. In this paper, the proposed DCIM is successfully applied to stochastic modeling of load balancing in networks of computing nodes allowing for prediction of the load distribution in the system, which is a novel application for the Influence Model. The DCIM is further used to identify the optimum workload distribution policy for load balancing in networked computing systems.


title: Real-Time Real Economic Activity: Entering and Exiting the Pandemic Recession of 2020

abstract: Entering and exiting the Pandemic Recession, I study the high-frequency real-activity signals provided by a leading nowcast, the ADS Index of Business Conditions produced and released in real time by the Federal Reserve Bank of Philadelphia. I track the evolution of real-time vintage beliefs and compare them to a later-vintage chronology. Real-time ADS plunges and then swings as its underlying economic indicators swing, but the ADS paths quickly converge to indicate a return to brisk positive growth by mid-May. I show, moreover, that the daily real activity path was highly correlated with the daily COVID-19 cases. Finally, I provide a comparative assessment of the real-time ADS signals provided when exiting the Great Recession.


title: Multiscale Microtubule Dynamics in Active Nematics

abstract: In microtubule-based active nematics, motor-driven extensile motion of microtubule bundles powers chaotic large-scale dynamics. We quantify the interfilament sliding motion both in isolated bundles and in a dense active nematic. The extension speed of an isolated microtubule pair is comparable to the molecular motor stepping speed. In contrast, the net extension in dense 2D active nematics is significantly slower; the interfilament sliding speeds are widely distributed about the average and the filaments exhibit both contractile and extensile relative motion. These measurements highlight the challenge of connecting the extension rate of isolated bundles to the multi-motor and multi-filament interactions present in a dense 2D active nematic. They also provide quantitative data that is essential for building multiscale models.


title: Hyperuniformity in cyclically driven glasses

abstract: We present a numerical investigation of the density fluctuations in a model glass under cyclic shear deformation. At low amplitude of shear, below yielding, the system reaches a steady absorbing state in which density fluctuations are suppressed revealing a clear fingerprint of hyperuniformity up to a finite length scale. The opposite scenario is observed above yielding, where the density fluctuations are strongly enhanced. We demonstrate that the transition to this state is accompanied by a spatial phase separation into two distinct hyperuniform regions, as a consequence of shear band formation above the yield amplitude.


title: Region-of-interest guided Supervoxel Inpainting for Self-supervision

abstract: Self-supervised learning has proven to be invaluable in making best use of all of the available data in biomedical image segmentation. One particularly simple and effective mechanism to achieve self-supervision is inpainting, the task of predicting arbitrary missing areas based on the rest of an image. In this work, we focus on image inpainting as the self-supervised proxy task, and propose two novel structural changes to further enhance the performance of a deep neural network. We guide the process of generating images to inpaint by using supervoxel-based masking instead of random masking, and also by focusing on the area to be segmented in the primary task, which we term as the region-of-interest. We postulate that these additions force the network to learn semantics that are more attuned to the primary task, and test our hypotheses on two applications: brain tumour and white matter hyperintensities segmentation. We empirically show that our proposed approach consistently outperforms both supervised CNNs, without any self-supervision, and conventional inpainting-based self-supervision methods on both large and small training set sizes.


title: A high-order well-balanced positivity-preserving moving mesh DG method for the shallow water equations with non-flat bottom topography

abstract: A rezoning-type adaptive moving mesh discontinuous Galerkin method is proposed for the numerical solution of the shallow water equations with non-flat bottom topography. The well-balance property is crucial to the simulation of perturbation waves over the lake-at-rest steady state such as waves on a lake or tsunami waves in the deep ocean. To ensure the well-balance and positivity-preserving properties, strategies are discussed in the use of slope limiting, positivity-preservation limiting, and data transferring between meshes. Particularly, it is suggested that a DG-interpolation scheme be used for the interpolation of both the flow variables and bottom topography from the old mesh to the new one and after each application of the positivity-preservation limiting on the water depth, a high-order correction be made to the approximation of the bottom topography according to the modifications in the water depth. Moreover, mesh adaptation based on the equilibrium variable and water depth is shown to give more desirable results than that based on the commonly used entropy function. Numerical examples in one and two spatial dimensions are presented to demonstrate the well-balance and positivity-preserving properties of the method and its ability to capture small perturbations of the lake-at-rest steady state.


title: Weakly $\omega$-Categorified Models of Algebraic Theories

abstract: We provide the expected constructions of weakly $\omega$-categorified models (in the sense of Bressie) of the theory of groups and quandles which arise by replacing the homotopies used to give equivalence relations in the theory of fundamental groups, fundamental quandles, and knot quandles with homotopies of all orders used as arrows of categorical dimensions one and greater, and discuss other related constructions of weakly $\omega$-categorifed algebras.


title: Interpretable Factorization for Neural Network ECG Models

abstract: The ability of deep learning (DL) to improve the practice of medicine and its clinical outcomes faces a looming obstacle: model interpretation. Without description of how outputs are generated, a collaborating physician can neither resolve when the model's conclusions are in conflict with his or her own, nor learn to anticipate model behavior. Current research aims to interpret networks that diagnose ECG recordings, which has great potential impact as recordings become more personalized and widely deployed. A generalizable impact beyond ECGs lies in the ability to provide a rich test-bed for the development of interpretive techniques in medicine. Interpretive techniques for Deep Neural Networks (DNNs), however, tend to be heuristic and observational in nature, lacking the mathematical rigor one might expect in the analysis of math equations. The motivation of this paper is to offer a third option, a scientific approach. We treat the model output itself as a phenomenon to be explained through component parts and equations governing their behavior. We argue that these component parts should also be "black boxes" --additional targets to interpret heuristically with clear functional connection to the original. We show how to rigorously factor a DNN into a hierarchical equation consisting of black box variables. This is not a subdivision into physical parts, like an organism into its cells; it is but one choice of an equation into a collection of abstract functions. Yet, for DNNs trained to identify normal ECG waveforms on PhysioNet 2017 Challenge data, we demonstrate this choice yields interpretable component models identified with visual composite sketches of ECG samples in corresponding input regions. Moreover, the recursion distills this interpretation: additional factorization of component black boxes corresponds to ECG partitions that are more morphologically pure.


title: Making DensePose fast and light

abstract: DensePose estimation task is a significant step forward for enhancing user experience computer vision applications ranging from augmented reality to cloth fitting. Existing neural network models capable of solving this task are heavily parameterized and a long way from being transferred to an embedded or mobile device. To enable Dense Pose inference on the end device with current models, one needs to support an expensive server-side infrastructure and have a stable internet connection. To make things worse, mobile and embedded devices do not always have a powerful GPU inside. In this work, we target the problem of redesigning the DensePose R-CNN model's architecture so that the final network retains most of its accuracy but becomes more light-weight and fast. To achieve that, we tested and incorporated many deep learning innovations from recent years, specifically performing an ablation study on 23 efficient backbone architectures, multiple two-stage detection pipeline modifications, and custom model quantization methods. As a result, we achieved $17\times$ model size reduction and $2\times$ latency improvement compared to the baseline model.


title: Is SGD a Bayesian sampler? Well, almost

abstract: Overparameterised deep neural networks (DNNs) are highly expressive and so can, in principle, generate almost any function that fits a training dataset with zero error. The vast majority of these functions will perform poorly on unseen data, and yet in practice DNNs often generalise remarkably well. This success suggests that a trained DNN must have a strong inductive bias towards functions with low generalisation error. Here we empirically investigate this inductive bias by calculating, for a range of architectures and datasets, the probability $P_{SGD}(f\mid S)$ that an overparameterised DNN, trained with stochastic gradient descent (SGD) or one of its variants, converges on a function $f$ consistent with a training set $S$. We also use Gaussian processes to estimate the Bayesian posterior probability $P_B(f\mid S)$ that the DNN expresses $f$ upon random sampling of its parameters, conditioned on $S$. Our main findings are that $P_{SGD}(f\mid S)$ correlates remarkably well with $P_B(f\mid S)$ and that $P_B(f\mid S)$ is strongly biased towards low-error and low complexity functions. These results imply that strong inductive bias in the parameter-function map (which determines $P_B(f\mid S)$), rather than a special property of SGD, is the primary explanation for why DNNs generalise so well in the overparameterised regime. While our results suggest that the Bayesian posterior $P_B(f\mid S)$ is the first order determinant of $P_{SGD}(f\mid S)$, there remain second order differences that are sensitive to hyperparameter tuning. A function probability picture, based on $P_{SGD}(f\mid S)$ and/or $P_B(f\mid S)$, can shed new light on the way that variations in architecture or hyperparameter settings such as batch size, learning rate, and optimiser choice, affect DNN performance.


title: Directed force propagation in semiflexible networks

abstract: We consider the propagation of tension along specific filament of a semiflexible filament network in response to the application of a point force using a combination of numerical simulations and analytic theory. We find the distribution of force within the network is highly heterogenous, with a small number of fibers supporting a significant fraction of the applied load over distances of multiple mesh sizes surrounding the point of force application. We suggest that these structures may be thought of as tensile force chains, whose structure we explore via simulation. We develop self-consistent calculations of the point-force response function and introduce a transfer matrix approach to explore the decay of tension (into bending) energy and the branching of tensile force chains in the network.


title: $L^2$-minimal extensions over Hermitian symmetric domains

abstract: In this paper, we study the $L^2$-minimal extension problem for polarized variations of Hodge structures over Hermitian symmetric domains. We are able to explicitly find the $L^2$-minimal extensions using a group-theoretic construction. In particular, this gives a construction without using $L^2$-estimates as in the Ohsawa-Takegoshi type extension theorems. The key ingredient is the Harish-Chandra embedding of Hermitian symmetric domains. The construction of holomorphic sections might be of independent interest since it gives a concrete description in the setting of Hermitian VHS.


title: Online learning with Corrupted context: Corrupted Contextual Bandits

abstract: We consider a novel variant of the contextual bandit problem (i.e., the multi-armed bandit with side-information, or context, available to a decision-maker) where the context used at each decision may be corrupted ("useless context"). This new problem is motivated by certain on-line settings including clinical trial and ad recommendation applications. In order to address the corrupted-context setting,we propose to combine the standard contextual bandit approach with a classical multi-armed bandit mechanism. Unlike standard contextual bandit methods, we are able to learn from all iteration, even those with corrupted context, by improving the computing of the expectation for each arm. Promising empirical results are obtained on several real-life datasets.


title: Uncovering socioeconomic gaps in mobility reduction during the COVID-19 pandemic using location data

abstract: Using smartphone location data from Colombia, Mexico, and Indonesia, we investigate how non-pharmaceutical policy interventions intended to mitigate the spread of the COVID-19 pandemic impact human mobility. In all three countries, we find that following the implementation of mobility restriction measures, human movement decreased substantially. Importantly, we also uncover large and persistent differences in mobility reduction between wealth groups: on average, users in the top decile of wealth reduced their mobility up to twice as much as users in the bottom decile. For decision-makers seeking to efficiently allocate resources to response efforts, these findings highlight that smartphone location data can be leveraged to tailor policies to the needs of specific socioeconomic groups, especially the most vulnerable.


title: A shell-model study of calcium isotopes towards their drip line

abstract: We report in this paper a study in terms of the nuclear shell model about the location of the calcium isotopes drip line. The starting point is considering the realistic two-body potential derived by Entem and Machleidt within chiral perturbation theory at next-to-next-to-next-to-leading order (N3LO), as well as a chiral three-body force at next-to-next-to-leading order (N2LO) whose structure and low-energy constants are consistent with the two-body potential. Then we construct the effective single-particle energies and residual interaction needed to diagonalize the shell-model Hamiltonian. The calculated two-neutron separation energies agree nicely with experiment until 56Ca, which is the heaviest isotope whose mass has been measured, and do not show any sign of two-neutron emission until 70Ca. We discuss the role of the choice of the model space in determining the neutron drip line, and also the dependence of the results on the parameters of the shell-model Hamiltonian.


title: Low-magnitude Seismicity with a Downhole Distributed Acoustic Sensing Array -- examples from the FORGE Geothermal Experiment

abstract: We show the capabilities of a downhole Distributed Acoustic Sensing (DAS) array in detecting, locating and characterizing low-magnitude earthquakes occurring in the vicinity of the Frontier Observatory for Research in Geothermal Energy (FORGE) site in Utah. 10.5 days of continuous data were acquired in a monitoring well at the FORGE geothermal site during the initial stimulation of an Enhanced Geothermal System in April-May 2019. Earthquake activity beneath Mineral Mountains, Utah also occurred within 10 km of the FORGE monitoring well. During the experiment, four events from those areas were cataloged by the University of Utah Seismograph Stations. Our processing of DAS data, including template matching, finds 82 earthquakes during that period, of which 16 are visible on the regional network. The magnitude of completeness obtained by DAS processing is better by at least M=0.5 than the dense surface array around the FORGE site. While a single vertical DAS array is limited in terms of event location due to its azimuthal ambiguity, multiple DAS wells or a combination of a downhole array with surface stations or near-surface horizontal DAS could jointly resolve locations. All detected events probably originated from the two active source areas and can be clustered into several distinct families.


title: Ionized and hot molecular outflows in the inner 500 pc of NGC1275

abstract: The role of feedback from Active Galactic Nuclei (AGN) in the evolution of galaxies is still not not fully understood, mostly due to the lack of observational constraints in the multi-phase gas kinematics on the ten to hundred parsec scales. We have used the Gemini Near-infrared Integral Field Spectrograph (NIFS) to map the molecular and ionized gas kinematics in the inner 900$\times$900 pc$^2$ of the Seyfert galaxy NGC1275 at a spatial resolution of $\sim$70 pc. From the fitting of the CO absorption bandheads in the K-band, we derive a stellar velocity dispersion of $265\pm26$ km s$^{-1}$, which implies a black hole mass of $M_{\rm SMBH}=1.1^{+0.9}_{-0.5}\times10^9$ M$_\odot$. We find hot ($T\gtrsim1000$ K) molecular and ionized outflows with velocities of up to 2 000 km s$^{-1}$ and mass outflow rates of $2.7\times10^{-2} {\rm M_\odot}$ yr$^{-1}$ and $1.6 {\rm M_\odot}$ yr$^{-1}$, respectively, in each of these gas phases. The kinetic power of the ionized outflows corresponds to only 0.05 per cent of the luminosity of the AGN of NGC 1275, indicating that they are not powerful enough to provide significant AGN feedback, but may be effective in redistributing the gas in the central region of the galaxy. The AGN driven outflows seem to be responsible for the shocks necessary to produce the observed H$_2$ and [Fe II] line emission.


title: DDPG++: Striving for Simplicity in Continuous-control Off-Policy Reinforcement Learning

abstract: This paper prescribes a suite of techniques for off-policy Reinforcement Learning (RL) that simplify the training process and reduce the sample complexity. First, we show that simple Deterministic Policy Gradient works remarkably well as long as the overestimation bias is controlled. This is contrast to existing literature which creates sophisticated off-policy techniques. Second, we pinpoint training instabilities, typical of off-policy algorithms, to the greedy policy update step; existing solutions such as delayed policy updates do not mitigate this issue. Third, we show that ideas in the propensity estimation literature can be used to importance-sample transitions from the replay buffer and selectively update the policy to prevent deterioration of performance. We make these claims using extensive experimentation on a set of challenging MuJoCo tasks. A short video of our results can be seen at https://tinyurl.com/scs6p5m .


title: Finite sets containing near-primitive roots

abstract: Fix $a \in \mathbb{Z}$, $a\notin \{0,\pm 1\}$. A simple argument shows that for each $\epsilon > 0$, and almost all (asymptotically 100% of) primes $p$, the multiplicative order of $a$ modulo $p$ exceeds $p^{\frac12-\epsilon}$. It is an open problem to show the same result with $\frac12$ replaced by any larger constant. We show that if $a,b$ are multiplicatively independent, then for almost all primes $p$, one of $a,b,ab, a^2b, ab^2$ has order exceeding $p^{\frac{1}{2}+\frac{1}{30}}$. The same method allows one to produce, for each $\epsilon > 0$, explicit finite sets $\mathcal{A}$ with the property that for almost all primes $p$, some element of $\mathcal{A}$ has order exceeding $p^{1-\epsilon}$. Similar results hold for orders modulo general integers $n$ rather than primes $p$.


title: What is Known About the Polarization of Starlight in the Southern Hole

abstract: Among the greatest mysteries in cosmology are the flatness problem, concerning the lack of curvature of the universe, and the homogeneity problem, questioning why the universe is almost isotropic despite having regions that are causally disconnected. These problems served as motivation for the theory of inflation, which suggests a period of exponential expansion in the early universe, and the inflationary origin of the universe can be traced by B-mode polarization. In an effort to better understand the potential foreground systematics, especially the levels of polarized dust emission, we queried the Heiles catalog to produce a list of starlight polarization data in the so-called "Southern Hole", which is an approximately $20\times20$ degree region centered at RA: $00^h12^m00^s$ and DEC: $-59{\deg}18'00''$ that is being examined by multiple CMB polarization experiments. Because magnetic field tends to dictate the orientation of dust grains, which in turn determines how starlight is polarized, starlight polarization can be used to trace magnetic fields. Therefore, to improve our understanding of the properties of this region, we used this catalog, along with Gaia data as tracers of the three-dimensional distribution of dust, as a potential indicator of magnetic field orientation throughout the galaxy in the Southern Hole region. We then analyzed these data with the hope that magnetic field data can be used to create a template to aid in subtracting away the contamination of CMB B-mode searches by polarized dust emission. While the results of the analysis are promising, we found that the currently available data are severely inadequate for the purpose of creating a template, thus demonstrating the need for improved and more uniform coverage of the Southern Hole when it comes to polarization measurements.


title: Likelihood Maximization and Moment Matching in Low SNR Gaussian Mixture Models

abstract: We derive an asymptotic expansion for the log likelihood of Gaussian mixture models (GMMs) with equal covariance matrices in the low signal-to-noise regime. The expansion reveals an intimate connection between two types of algorithms for parameter estimation: the method of moments and likelihood optimizing algorithms such as Expectation-Maximization (EM). We show that likelihood optimization in the low SNR regime reduces to a sequence of least squares optimization problems that match the moments of the estimate to the ground truth moments one by one. This connection is a stepping stone toward the analysis of EM and maximum likelihood estimation in a wide range of models. A motivating application for the study of low SNR mixture models is cryo-electron microscopy data, which can be modeled as a GMM with algebraic constraints imposed on the mixture centers. We discuss the application of our expansion to algebraically constrained GMMs, among other example models of interest.


title: The Galactic Plane Infrared Polarization Survey (GPIPS): Data Release 4

abstract: The Galactic Plane Infrared Polarization Survey (GPIPS) seeks to characterize the magnetic field in the dusty Galactic disk using near-infrared stellar polarimetry. All GPIPS observations were completed using the 1.83 m Perkins telescope and Mimir instrument. GPIPS observations surveyed 76 sq-deg of the northern Galactic plane, from Galactic longitudes 18 to 56 deg and latitudes -1 to +1 deg, in the H-band (1.6 um). Surveyed stars span 7th to 16th mag, resulting in nearly 10 million stars with measured linear polarizations. Of these stars, ones with m_H < 12.5 mag and polarization percentage uncertainties under 2% were judged to be high quality and number over one million. GPIPS data reveal plane-of-sky magnetic field orientations for numerous interstellar clouds for AV values to ~30 mag. The average sky separation of stars with m_H < 12.5 mag is about 30arcsec, or about 60 per Planck polarization resolution element. Matching to Gaia DR2 showed the brightest GPIPS stars are red giants with distances in the 0.6-7.5 kpc range. Polarization orientations are mostly parallel to the Galactic disk, with some zones showing significant orientation departures. Changes in orientations are stronger as a function of Galactic longitude than of latitude. Considered at 10 arcmin angular scales, directions that show the greatest polarization fractions and narrowest polarization position angle distributions are confined to about ten large, coherent structures that are not correlated with star forming clouds. The GPIPS polarimetric and photometric data products (Data Release 4 catalogs and images) are publicly available for over 13million stars.


title: Marcus Theory of Thermoelectricity in Molecular Junctions

abstract: Thermoelectric energy conversion is perhaps the most promising of the potential applications of molecular electronics. Ultimately, it is desirable for this technology to operate at around room temperature, and it is therefore important to consider the role of dissipative effects in these conditions. Here, we develop a theory of thermoelectricity which accounts for the vibrational coupling within the framework of Marcus theory. We demonstrate that the inclusion of lifetime broadening is necessary in the theoretical description of this phenomenon. We further show that the Seebeck coefficient and the power factor decrease with increasing reorganisation energy, and identify the optimal operating conditions in the case of non-zero reorganisation energy. Finally, with the aid of DFT calculations, we consider a prototypical fullerene-based molecular junction. We estimate the maximum power factor that can be obtained in this system, and confirm that C$_{60}$ is an excellent candidate for thermoelectric heat-to-energy conversion. This work provides general guidance that should be followed in order to achieve high-efficiency molecular thermoelectric materials.


title: Exclusive Diffraction at the LHC

abstract: In this report, we describe the most recent results on exclusive diffraction from the ATLAS, CMS, LHCb, TOTEM experiments at the LHC concerning exclusive pions, $J/\Psi$, $\Psi(2S)$, dilepton, diphoton, $WW$ productions and prospects concerning the search for anomalous couplings and axion-like particle production.


title: Compact sets with large projections and nowhere dense sumset

abstract: We answer a question of Banakh, Jab\l{}o\'nska and Jab\l{}o\'nski by showing that for $d\ge 2$ there exists a compact set $K \subseteq \mathbb{R}^d$ such that the projection of $K$ onto each hyperplane is of non-empty interior, but $K+K$ is nowhere dense. The proof relies on a random construction. A natural approach in the proofs is to construct such a $K$ in the unit cube with full projections, that is, such that the projections of $K$ agree with that of the unit cube. We investigate the generalization of these problems for projections onto various dimensional subspaces as well as for $\ell$-fold sumsets. We obtain numerous positive and negative results, but also leave open many interesting cases. We also show that in most cases if we have a specific example of such a compact set then actually the generic (in the sense of Baire category) compact set in a suitably chosen space is also an example. Finally, utilizing a computer-aided construction, we show that the compact set in the plane with full projections and nowhere dense sumset can be self-similar.


title: ATOM: Robustifying Out-of-distribution Detection Using Outlier Mining

abstract: Detecting out-of-distribution (OOD) inputs is critical for safely deploying deep learning models in an open-world setting. However, existing OOD detection solutions can be brittle in the open world, facing various types of adversarial OOD inputs. While methods leveraging auxiliary OOD data have emerged, our analysis on illuminative examples reveals a key insight that the majority of auxiliary OOD examples may not meaningfully improve or even hurt the decision boundary of the OOD detector, which is also observed in empirical results on real data. In this paper, we provide a theoretically motivated method, Adversarial Training with informative Outlier Mining (ATOM), which improves the robustness of OOD detection. We show that, by mining informative auxiliary OOD data, one can significantly improve OOD detection performance, and somewhat surprisingly, generalize to unseen adversarial attacks. ATOM achieves state-of-the-art performance under a broad family of classic and adversarial OOD evaluation tasks. For example, on the CIFAR-10 in-distribution dataset, ATOM reduces the FPR (at TPR 95%) by up to 57.99% under adversarial OOD inputs, surpassing the previous best baseline by a large margin.


title: Precision constraints for three-flavor neutrino oscillations from the full MINOS+ and MINOS data set

abstract: We report the final measurement of the neutrino oscillation parameters $\Delta m^2_{32}$ and $\sin^2\theta_{23}$ using all data from the MINOS and MINOS+ experiments. These data were collected using a total exposure of $23.76 \times 10^{20}$ protons on target producing $\nu_{mu}$ and $\overline{\nu_\mu}$ beams and 60.75 kt$\cdot$yr exposure to atmospheric neutrinos. The measurement of the disappearance of $\nu_{\mu}$ and the appearance of $\nu_e$ events between the Near and Far detectors yields $|\Delta m^2_{32}|=2.40^{+0.08}_{-0.09}~(2.45^{+0.07}_{-0.08}) \times 10^{-3}$ eV$^2$ and $\sin^2\theta_{23} = 0.43^{+0.20}_{-0.04} ~(0.42^{+0.07}_{-0.03})$ at 68% C.L. for Normal (Inverted) Hierarchy.


title: Probing BFKL dynamics at hadronic colliders in jet gap jet events

abstract: In this report, we give the Balitsky Fadin Kuraev Lipton formalism for jet gap jet events at hadronic colliders. We also discuss the case where in addition at least one proton is intact in the final state in diffractive events.


title: Automatic Integration

abstract: The purpose of this paper is to introduce the concept of the automatic integration and present a new way of approximating definite integrals using the automatic integration based on an associative algebra with zero divisors.


title: Reply to Comment on "Comparison of the Lateral Retention Forces on Sessile,Pendant, and Inverted Sessile Drops"

abstract: We address the issues raised in [R. Tadmor et al., Langmuir 2020, 36, 475-476]. In particular, we explain why we did not use Tadmor's theory to explain our results.


title: Hybrid Models for Learning to Branch

abstract: A recent Graph Neural Network (GNN) approach for learning to branch has been shown to successfully reduce the running time of branch-and-bound algorithms for Mixed Integer Linear Programming (MILP). While the GNN relies on a GPU for inference, MILP solvers are purely CPU-based. This severely limits its application as many practitioners may not have access to high-end GPUs. In this work, we ask two key questions. First, in a more realistic setting where only a CPU is available, is the GNN model still competitive? Second, can we devise an alternate computationally inexpensive model that retains the predictive power of the GNN architecture? We answer the first question in the negative, and address the second question by proposing a new hybrid architecture for efficient branching on CPU machines. The proposed architecture combines the expressive power of GNNs with computationally inexpensive multi-layer perceptrons (MLP) for branching. We evaluate our methods on four classes of MILP problems, and show that they lead to up to 26% reduction in solver running time compared to state-of-the-art methods without a GPU, while extrapolating to harder problems than it was trained on. The code for this project is publicly available at https://github.com/pg2455/Hybrid-learn2branch.


title: Simulating human interactions in supermarkets to measure the risk of COVID-19 contagion at scale

abstract: Taking the context of simulating a retail environment using agent based modelling, a theoretical model is presented that describes the probability distribution of customer "collisions" using a novel space transformation to the Torus $Tor^2$. A method for generating the distribution of customer paths based on historical basket data is developed. Finally a calculation of the number of simulations required for statistical significance is developed. An implementation of this modelling approach to run simulations on multiple store geometries at industrial scale is being developed with current progress detailed in the technical appendix.


title: Improving MF-DFA model with applications in precious metals market

abstract: With the aggravation of the global economic crisis and inflation, the precious metals with safe-haven function have become more popular. An improved MF-DFA method is proposed to analyze price fluctuations of the precious metals market. Based on the widely used multifractal detrended fluctuation analysis method (MF-DFA), we compare these two methods and find that the Bi-OSW-MF-DFA method possesses better efficiency. This article analyzes the degree of multifractality between spot gold market and spot silver market as well as their risks. From the numerical results and figures, it is found that two elements constitute the contributions in the formation of multifractality in time series and the risk of the spot silver market is higher than that of the spot gold market. This attempt could lead to a better understanding of complicated precious metals market.


title: A refined Gallai-Edmonds structure theorem for weighted matching polynomials

abstract: In this work, we prove a refinement of the Gallai-Edmonds structure theorem for weighted matching polynomials by Ku and Wong. Our proof uses a connection between matching polynomials and branched continued fractions. We also show how this is related to a modification by Sylvester of the classical Sturm's theorem on the number of zeros of a real polynomial in an interval. In addition, we obtain some other results about zeros of matching polynomials.


title: Mapping the South African health landscape in response to COVID-19

abstract: When the COVID-19 disease pandemic infiltrated the world, there was an immediate need for accurate information. As with any outbreak, the outbreak follows a clear trajectory, and subsequently, the supporting information for that outbreak needs to address the needs associated with that stage of the outbreak. At first, there was a need to inform the public of the information related to the initial situation related to the "who" of the COVID-19 disease. However, as time continued, the "where", "when" and "how to" related questions started to emerge in relation to the public healthcare system themselves. Questions surrounding the health facilities including COVID-19 hospital bed capacity, locations of designated COVID-19 facilities, and general information related to these facilities were not easily accessible to the general public. Furthermore, the available information was found to be outdated, fragmented across several platforms, and still had gaps in the data related to these facilities. To rectify this problem, a group of volunteers working on the covid19za project stepped in to assist. Each member leading a part of the project chose to focus on one of four problems related to the challenges associated with the Hospital information including: data quality, data completeness, data source validation and data visualisation capacity. As the project developed, so did the sophistication of the data, visualisation and core function of the project. The future prospects of this project relate to a Progressive Web Application that will avail this information for the public as well as healthcare workers through comprehensive mapping and data quality.


title: Exploring color-octet scalar parameter space in minimal $R$-symmetric models

abstract: In this work we study the collider phenomenology of color-octet scalars (sgluons) in minimal supersymmetric models endowed with a global continuous $R$ symmetry. We systematically catalog the significant decay channels of scalar and pseudoscalar sgluons and identify novel features that are natural in these models. These include decays in nonstandard diboson channels, such as to a gluon and a photon; three-body decays with considerable branching fractions; and long-lived particles with displaced vertex signatures. We also discuss the single and pair production of these particles and show that they can evade existing constraints from the Large Hadron Collider, to varying extents, in large regions of reasonable parameter space. We find, for instance, that a 725 GeV scalar and a 350 GeV or lighter pseudoscalar can still be accommodated in realistic scenarios.


title: Traditional and accelerated gradient descent for neural architecture search

abstract: In this paper we introduce two algorithms for neural architecture search (NASGD and NASAGD) following the theoretical work by two of the authors [5] which used the geometric structure of optimal transport to introduce the conceptual basis for new notions of traditional and accelerated gradient descent algorithms for the optimization of a function on a semi-discrete space. Our algorithms, which use the network morphism framework introduced in [2] as a baseline, can analyze forty times as many architectures as the hill climbing methods [2, 14] while using the same computational resources and time and achieving comparable levels of accuracy. For example, using NASGD on CIFAR-10, our method designs and trains networks with an error rate of 4.06 in only 12 hours on a single GPU.


title: Two-photon spontaneous emission in atomically thin plasmonic nanostructures

abstract: The ability to harness light-matter interactions at the few-photon level plays a pivotal role in quantum technologies. Single photons - the most elementary states of light - can be generated on-demand in atomic and solid state emitters. Two-photon states are also key quantum assets, but achieving them in individual emitters is challenging because their generation rate is much slower than competing one-photon processes. We demonstrate that atomically thin plasmonic nanostructures can harness two-photon spontaneous emission, resulting in giant far-field two-photon production, a wealth of resonant modes enabling tailored photonic and plasmonic entangled states, and plasmon-assisted single-photon creation orders of magnitude more efficient than standard one-photon emission. We unravel the two-photon spontaneous emission channels and show that their spectral line-shapes emerge from an intricate interplay between Fano and Lorentzian resonances. Enhanced two-photon spontaneous emission in two-dimensional nanostructures paves the way to an alternative efficient source of light-matter entanglement for on-chip quantum information processing and free-space quantum communications.


title: Recent results from the TOTEM collaboration at the LHC

abstract: We describe the most recent results from the TOTEM collaboration at the LHC, namely the elastic cross section measurements at a center-of-mass on 2.76, 7, 8 and 13 TeV. No structure or resonance is observed at high $t$ at high center-of-mass energies. A pure exponential form of $d \sigma/dt$ is excluded both at 8 and 13 TeV. Accessing the very low $t$ region allows measuring the $\rho$ parameter at 13 TeV.


title: Semi-discrete optimization through semi-discrete optimal transport: a framework for neural architecture search

abstract: In this paper we introduce a theoretical framework for semi-discrete optimization using ideas from optimal transport. Our primary motivation is in the field of deep learning, and specifically in the task of neural architecture search. With this aim in mind, we discuss the geometric and theoretical motivation for new techniques for neural architecture search (in a companion paper we show that algorithms inspired by our framework are competitive with contemporaneous methods). We introduce a Riemannian-like metric on the space of probability measures over a semi-discrete space $\mathbb{R}^d \times \mathcal{G}$ where $\mathcal{G}$ is a finite weighted graph. With such Riemmanian structure in hand, we derive formal expressions for the gradient flow of a relative entropy functional, as well as second order dynamics for the optimization of said energy. Then, with the aim of providing a rigorous motivation for the gradient flow equations derived formally, we also consider an iterative procedure known as minimizing movement scheme (i.e., Implicit Euler scheme, or JKO scheme) and apply it to the relative entropy with respect to a suitable cost function. For some specific choices of metric and cost, we rigorously show that the minimizing movement scheme of the relative entropy functional converges to the gradient flow process provided by the formal Riemannian structure. This flow coincides with a system of reaction-diffusion equations on $\mathbb{R}^d$.


title: BERTology Meets Biology: Interpreting Attention in Protein Language Models

abstract: Transformer architectures have proven to learn useful representations for protein classification and generation tasks. However, these representations present challenges in interpretability. In this work, we demonstrate a set of methods for analyzing protein Transformer models through the lens of attention. We show that attention: (1) captures the folding structure of proteins, connecting amino acids that are far apart in the underlying sequence, but spatially close in the three-dimensional structure, (2) targets binding sites, a key functional component of proteins, and (3) focuses on progressively more complex biophysical properties with increasing layer depth. We find this behavior to be consistent across three Transformer architectures (BERT, ALBERT, XLNet) and two distinct protein datasets. We also present a three-dimensional visualization of the interaction between attention and protein structure. Code for visualization and analysis is available at https://github.com/salesforce/provis.


title: Perception-Prediction-Reaction Agents for Deep Reinforcement Learning

abstract: We introduce a new recurrent agent architecture and associated auxiliary losses which improve reinforcement learning in partially observable tasks requiring long-term memory. We employ a temporal hierarchy, using a slow-ticking recurrent core to allow information to flow more easily over long time spans, and three fast-ticking recurrent cores with connections designed to create an information asymmetry. The \emph{reaction} core incorporates new observations with input from the slow core to produce the agent's policy; the \emph{perception} core accesses only short-term observations and informs the slow core; lastly, the \emph{prediction} core accesses only long-term memory. An auxiliary loss regularizes policies drawn from all three cores against each other, enacting the prior that the policy should be expressible from either recent or long-term memory. We present the resulting \emph{Perception-Prediction-Reaction} (PPR) agent and demonstrate its improved performance over a strong LSTM-agent baseline in DMLab-30, particularly in tasks requiring long-term memory. We further show significant improvements in Capture the Flag, an environment requiring agents to acquire a complicated mixture of skills over long time scales. In a series of ablation experiments, we probe the importance of each component of the PPR agent, establishing that the entire, novel combination is necessary for this intriguing result.


title: Evaluation of Embedded Platforms for Lower Limb Prosthesis with Visual Sensing Capabilities

abstract: Lower limb prosthesis can benefit from embedded systems capable of applying computer vision techniques to enhance autonomous control and context awareness for intelligent decision making. In order to fill in the gap of current literature of autonomous systems for prosthetic legs employing computer vision methods, we evaluate the performance capabilities of two off-the-shelf platforms, the Jetson TX2 and Raspberry Pi 3, by assessing their CPU load, memory usage, run time and classification accuracy for different image sizes and widespread computer vision algorithms. We make use of a dataset that we collected for terrain recognition using images from a camera mounted on a leg, which would enable context-awareness for lower limb prosthesis. We show that, given reasonably large images and an appropriate frame selection method, it is possible to identify the terrain that a subject is going to step on, with the possibility of reconstructing the surface and obtaining its inclination. This is part of a proposed system equipped with an embedded camera and inertial measurement unit to recognize different types of terrain and estimate the slope of the surface in front of the subject.


title: The ratio-cut polytope and K-means clustering

abstract: We introduce the ratio-cut polytope defined as the convex hull of ratio-cut vectors corresponding to all partitions of $n$ points in $\mathbb R^m$ into at most $K$ clusters. This polytope is closely related to the convex hull of the feasible region of a number of clustering problems such as K-means clustering and spectral clustering. We study the facial structure of the ratio-cut polytope and derive several types of facet-defining inequalities. We then consider the problem of K-means clustering and introduce a novel linear programming (LP) relaxation for it. Subsequently, we focus on the case of two clusters and derive sufficient conditions under which the proposed LP relaxation recovers the underlying clusters exactly. Namely, we consider the stochastic ball model, a popular generative model for K-means clustering, and we show that if the separation distance between cluster centers satisfies $\Delta > 1+\sqrt 3$, then the LP relaxation recovers the planted clusters with high probability. This is a major improvement over the only existing recovery guarantee for an LP relaxation of K-means clustering stating that recovery is possible with high probability if and only if $\Delta > 4$. Our numerical experiments indicate that the proposed LP relaxation significantly outperforms a popular semidefinite programming relaxation in recovering the planted clusters.


title: Riemannian Optimization on the Symplectic Stiefel Manifold

abstract: The symplectic Stiefel manifold, denoted by $\mathrm{Sp}(2p,2n)$, is the set of linear symplectic maps between the standard symplectic spaces $\mathbb{R}^{2p}$ and $\mathbb{R}^{2n}$. When $p=n$, it reduces to the well-known set of $2n\times 2n$ symplectic matrices. Optimization problems on $\mathrm{Sp}(2p,2n)$ find applications in various areas, such as optics, quantum physics, numerical linear algebra and model order reduction of dynamical systems. The purpose of this paper is to propose and analyze gradient-descent methods on $\mathrm{Sp}(2p,2n)$, where the notion of gradient stems from a Riemannian metric. We consider a novel Riemannian metric on $\mathrm{Sp}(2p,2n)$ akin to the canonical metric of the (standard) Stiefel manifold. In order to perform a feasible step along the antigradient, we develop two types of search strategies: one is based on quasi-geodesic curves, and the other one on the symplectic Cayley transform. The resulting optimization algorithms are proved to converge globally to critical points of the objective function. Numerical experiments illustrate the efficiency of the proposed methods.


title: Bypassing the computational bottleneck of quantum-embedding theories for strong electron correlations with machine learning

abstract: A cardinal obstacle to performing quantum-mechanical simulations of strongly-correlated matter is that, with the theoretical tools presently available, sufficiently-accurate computations are often too expensive to be ever feasible. Here we design a computational framework combining quantum-embedding (QE) methods with machine learning. This allows us to bypass altogether the most computationally-expensive components of QE algorithms, making their overall cost comparable to bare Density Functional Theory (DFT). We perform benchmark calculations of a series of actinide systems, where our method describes accurately the correlation effects, reducing by orders of magnitude the computational cost. We argue that, by producing a larger-scale set of training data, it will be possible to apply our method to systems with arbitrary stoichiometries and crystal structures, paving the way to virtually infinite applications in condensed matter physics, chemistry and materials science.


title: HypervolGAN: An efficient approach for GAN with multi-objective training function

abstract: Since the advent of generative adversarial networks (GANs), various loss functions have been developed and combined to constitute the overall training objective function, in order to improve model performance or for specific learning tasks. For instance, in image enhancement or restoration, there are often several criteria to consider such as signal-noise ratio, smoothness, structures and details. However, when the optimization goal has more than one adversarial loss, balancing multiple losses in the overall function becomes a challenging, critical and time-consuming problem. In this paper, we propose to tackle the problem by means of efficient multi-objective optimization. The proposed HypervolGAN adopts an adapted version of hypervolume maximization method to effectively define the multi-objective training function for GAN. We tested our proposed method on solving single image super-resolution problem. Experiments show that the proposed HypervolGAN is efficient in saving computational time and efforts for fine-tuning weights of various losses, and can generate enhanced samples that have better quality than results given by baseline GANs. The work explores the integration of adversarial learning and optimization techniques, which can benefit not only image processing but also a wide range of applications.


title: CheXpert++: Approximating the CheXpert labeler for Speed,Differentiability, and Probabilistic Output

abstract: It is often infeasible or impossible to obtain ground truth labels for medical data. To circumvent this, one may build rule-based or other expert-knowledge driven labelers to ingest data and yield silver labels absent any ground-truth training data. One popular such labeler is CheXpert, a labeler that produces diagnostic labels for chest X-ray radiology reports. CheXpert is very useful, but is relatively computationally slow, especially when integrated with end-to-end neural pipelines, is non-differentiable so can't be used in any applications that require gradients to flow through the labeler, and does not yield probabilistic outputs, which limits our ability to improve the quality of the silver labeler through techniques such as active learning. In this work, we solve all three of these problems with $\texttt{CheXpert++}$, a BERT-based, high-fidelity approximation to CheXpert. $\texttt{CheXpert++}$ achieves 99.81\% parity with CheXpert, which means it can be reliably used as a drop-in replacement for CheXpert, all while being significantly faster, fully differentiable, and probabilistic in output. Error analysis of $\texttt{CheXpert++}$ also demonstrates that $\texttt{CheXpert++}$ has a tendency to actually correct errors in the CheXpert labels, with $\texttt{CheXpert++}$ labels being more often preferred by a clinician over CheXpert labels (when they disagree) on all but one disease task. To further demonstrate the utility of these advantages in this model, we conduct a proof-of-concept active learning study, demonstrating we can improve accuracy on an expert labeled random subset of report sentences by approximately 8\% over raw, unaltered CheXpert by using one-iteration of active-learning inspired re-training. These findings suggest that simple techniques in co-learning and active learning can yield high-quality labelers under minimal, and controllable human labeling demands.


title: Dependence of the density of states outer measure on the potential for deterministic Schr\"odinger operators on graphs with applications to ergodic and random models

abstract: We continue our study of the dependence of the density of states measure and related spectral functions of Schr\"odinger operators on the potential. Whereas our earlier work focused on random Schr\"odinger operators, we extend these results to Schr\"odinger operators on infinite graphs with deterministic potentials and ergodic potentials, and improve our results for random potentials. In particular, we prove the Lipschitz continuity of the DOSm for random Schr\"odinger operators on the lattice, recovering results of \cite{kachkovskiy, shamis}. For our treatment of deterministic potentials, we first study the density of states outer measure (DOSoM), defined for all Schr\"odinger operators, and prove a deterministic result of the modulus of continuity of the DOSoM with respect to the potential. We apply these results to Schr\"odinger operators on the lattice $Z^d$ and the Bethe lattice. In the former case, we prove the Lipschitz continuity of the DOSoM, and in the latter case, we prove that the DOSoM is $\frac{1}{2}$-log-H\"older continuous. Our technique combines the abstract Lipschitz property of one-parameter families of self-adjoint operators with a new finite-range reduction that allows us to study the dependency of the DOSoM and related functions on only finitely-many variables and captures the geometry of the graph at infinity.


title: A Random Walk Model for Dark Matter Halo Concentrations

abstract: For idealized (spherical, smooth) dark matter halos described by single-parameter density profiles (such as the NFW profile) there exists a one-to-one mapping between the energy of the halo and the scale radius of its density profile. The energy therefore uniquely determines the concentration parameter of such halos. We exploit this fact to predict the concentrations of dark matter halos via a random walk in halo energy space. Given a full merger tree for a halo, the total internal energy of each halo in that tree is determined by summing the internal and orbital energies of progenitor halos. We show that, when calibrated, this model can accurately reproduce the mean of the concentration--mass relation measured in N-body simulations, and reproduces more of the scatter in that relation than previous models. We further test this model by examining both the autocorrelation of scale radii across time, and the correlations between halo concentration and spin, and comparing to results measured from cosmological N-body simulations. In both cases we find that our model closely matches the N-body results. Our model is implemented within the open source Galacticus toolkit.


title: The classification of symmetry protected topological phases of one-dimensional fermion systems

abstract: We introduce an index for symmetry protected topological (SPT) phases of infinite fermionic chains with an on-site symmetry given by a finite group $G$. This index takes values in $\mathbb{Z}_2 \times H^1(G,\mathbb{Z}_2) \times H^2(G, U(1)_{\mathfrak{p}})$ with a generalized Wall group law under stacking. We show that this index is an invariant of the classification of SPT phases. When the ground state is translation invariant and has reduced density matrices with uniformly bounded rank on finite intervals, we derive a fermionic matrix product representative of this state with on-site symmetry.

categories: math-ph cond-mat.stat-mech math.MP math.OA quant-ph


title: On the Relationship Between Probabilistic Circuits and Determinantal Point Processes

abstract: Scaling probabilistic models to large realistic problems and datasets is a key challenge in machine learning. Central to this effort is the development of tractable probabilistic models (TPMs): models whose structure guarantees efficient probabilistic inference algorithms. The current landscape of TPMs is fragmented: there exist various kinds of TPMs with different strengths and weaknesses. Two of the most prominent classes of TPMs are determinantal point processes (DPPs) and probabilistic circuits (PCs). This paper provides the first systematic study of their relationship. We propose a unified analysis and shared language for discussing DPPs and PCs. Then we establish theoretical barriers for the unification of these two families, and prove that there are cases where DPPs have no compact representation as a class of PCs. We close with a perspective on the central problem of unifying these tractable models.


title: Efficient 2D Tensor Network Simulation of Quantum Systems

abstract: Simulation of quantum systems is challenging due to the exponential size of the state space. Tensor networks provide a systematically improvable approximation for quantum states. 2D tensor networks such as Projected Entangled Pair States (PEPS) are well-suited for key classes of physical systems and quantum circuits. However, direct contraction of PEPS networks has exponential cost, while approximate algorithms require computations with large tensors. We propose new scalable algorithms and software abstractions for PEPS-based methods, accelerating the bottleneck operation of contraction and refactorization of a tensor subnetwork. We employ randomized SVD with an implicit matrix to reduce cost and memory footprint asymptotically. Further, we develop a distributed-memory PEPS library and study accuracy and efficiency of alternative algorithms for PEPS contraction and evolution on the Stampede2 supercomputer. We also simulate a popular near-term quantum algorithm, the Variational Quantum Eigensolver (VQE), and benchmark Imaginary Time Evolution (ITE), which compute ground states of Hamiltonians.


title: The stabilizers for the action of orthogonal similarity on symmetric matrices and orthogonal $*$-conjugacy on Hermitian matrices

abstract: We describe the recursive algorithmic procedure to compute the stabilizers of the group of complex orthogonal matrices with respect to the action of similarity on the set of all symmetric matrices. Futhermore, lower bounds for dimensions of the stabilizers for the action of orthogonal $*$-conjugation on Hermitian matrices are obtained. We also prove a result that completes the classification of normal forms of Hermitian matrices under orthogonal $*$-conjugation. A key step in our proof is to solve a certain block matrix equation with Toeplitz blocks. These results are then applied to provide a theorem on normal forms of the quadratic parts of flat complex points in a real codimension $2$ submanifold in a complex manifold.


title: Orthogonal polynomials and Hankel Determinants for certain Bernoulli and Euler Polynomials

abstract: Using continued fraction expansions of certain polygamma functions as a main tool, we find orthogonal polynomials with respect to the odd-index Bernoulli polynomials $B_{2k+1}(x)$ and the Euler polynomials $E_{2k+\nu}(x)$, for $\nu=0, 1, 2$. In the process we also determine the corresponding Jacobi continued fractions (or J-fractions) and Hankel determinants. In all these cases the Hankel determinants are polynomials in $x$ which factor completely over the rationals.


title: Why fractional derivatives with nonsingular kernels should not be used

abstract: In recent years, many papers discuss the theory and applications of new fractional-order derivatives that are constructed by replacing the singular kernel of the Caputo or Riemann-Liouville derivative by a non-singular (i.e., bounded) kernel. It will be shown here, through rigorous mathematical reasoning, that these non-singular kernel derivatives suffer from several drawbacks which should forbid their use. They fail to satisfy the fundamental theorem of fractional calculus since they do not admit the existence of a corresponding convolution integral of which the derivative is the left-inverse; and the value of the derivative at the initial time $t=0$ is always zero, which imposes an unnatural restriction on the differential equations and models where these derivatives can be used. For the particular cases of the so-called Caputo-Fabrizio and Atangana-Baleanu derivatives, it is shown that when this restriction holds the derivative can be simply expressed in terms of integer derivatives and standard Caputo fractional derivatives, thus demonstrating that these derivatives contain nothing new.


title: Non-adiabatic quantum interference effects and chaoticity in the ultracold Li + LiNa $\to$ Li$_2$ + Na reaction

abstract: Electronically non-adiabatic effects play an important role in many chemical reactions. How these effects manifest in cold and ultracold chemistry remain largely unexplored. Here, through first principles non-adiabatic quantum dynamics calculations of the Li + LiNa $\to$ Li$_2$ + Na chemical reaction, it is shown that non-adiabatic dynamics induces quantum interference effects that dramatically alter the ultracold rotationally resolved reaction rate coefficients. The interference effect arises from a conical intersection between the ground and an excited electronic state that is energetically accessible even for ultracold collisions. These unique interference effects might be exploited for quantum control applications as a quantum molecular switch. A statistical analysis of rotational populations of the Li$_2$ product reveals a Poisson distribution implying an underlying classically chaotic dynamics. The Poisson distribution is robust and amenable to experimental verification and appears to be a universal property of ultracold reactions involving alkali metal dimers.


title: Decay of two-dimensional quantum turbulence in binary Bose-Einstein condensates

abstract: We study two-dimensional quantum turbulence in miscible binary Bose-Einstein condensates in either a harmonic trap or a steep-wall trap through the numerical simulations of the Gross-Pitaevskii equations. The turbulence is generated through a Gaussian stirring potential. When the condensates have unequal intra-component coupling strengths or asymmetric trap frequencies, the turbulent condensates undergo a dramatic decay dynamics to an interlaced array of vortex-antidark structures, a quasi-equilibrium state, of like-signed vortices with an extended size of the vortex core. The time of formation of this state is shortened when the parameter asymmetry of the intra-component couplings or the trap frequencies are enhanced. The corresponding spectrum of the incompressible kinetic energy exhibits two noteworthy features: (i) a $k^{-3}$ power-law around the range of the wave number determined by the spin healing length (the size of the extended vortex-core) and (ii) a flat region around the range of the wave number determined by the density healing length. The latter is associated with the small scale phase fluctuation relegated outside the Thomas-Fermi radius and is more prominent as the strength of intercomponent interaction approaches the strength of intra-component interaction. We also study the impact of the inter-component interaction to the cluster formation of like-signed vortices in an elliptical steep-wall trap, finding that the inter-component coupling gives rise to the decay of the clustered configuration.

categories: cond-mat.quant-gas nlin.CD nlin.PS physics.flu-dyn