title: Realizing the "fictitious" beam splitter -- A stationary implementation of semi-counterfactual interaction-free imaging

abstract: Based on quantum counterfactual interaction-free measurement, we propose an implementation scheme for a beam splitter with anomalous reflection and transmission properties that looks impossible at first glance. Our scheme is stationary without requiring switchable mirrors and polarization rotators. Using the scheme for imaging will ensure that the optical radiation received by the object being imaged can be arbitrarily low. Thus it enables applications such as stealthy night vision devices that can work without detectable ambient light, or being used as a hackware against some counterfactual quantum cryptographic protocols.


title: Flows in near algebras with applications to harnesses

abstract: We introduce one-way flows in near algebras and two-way flows in double near algebras with two interrelated multiplications. We establish parametric representations of the one-way and two-way flows in terms of a single element of the algebra that we call a flow generator. We indicate probabilistic applications of the one-way flows to a study of polynomial stochastic processes. We apply our results on the two-way flows to harnesses and quadratic harnesses in probability theory, generalizing some previous results.


title: Scalable, Proposal-free Instance Segmentation Network for 3D Pixel Clustering and Particle Trajectory Reconstruction in Liquid Argon Time Projection Chambers

abstract: Liquid Argon Time Projection Chambers (LArTPCs) are high resolution particle imaging detectors, employed by accelerator-based neutrino oscillation experiments for high precision physics measurements. While images of particle trajectories are intuitive to analyze for physicists, the development of a high quality, automated data reconstruction chain remains challenging. One of the most critical reconstruction steps is particle clustering: the task of grouping 3D image pixels into different particle instances that share the same particle type. In this paper, we propose the first scalable deep learning algorithm for particle clustering in LArTPC data using sparse convolutional neural networks (SCNN). Building on previous works on SCNNs and proposal free instance segmentation, we build an end-to-end trainable instance segmentation network that learns an embedding of the image pixels to perform point cloud clustering in a transformed space. We benchmark the performance of our algorithm on PILArNet, a public 3D particle imaging dataset, with respect to common clustering evaluation metrics. 3D pixels were successfully clustered into individual particle trajectories with 90% of them having an adjusted Rand index score greater than 92% with a mean pixel clustering efficiency and purity above 96%. This work contributes to the development of an end-to-end optimizable full data reconstruction chain for LArTPCs, in particular pixel-based 3D imaging detectors including the near detector of the Deep Underground Neutrino Experiment. Our algorithm is made available in the open access repository, and we share our Singularity software container, which can be used to reproduce our work on the dataset.


title: Discrete self-similarity in the formation of satellites for viscous cavity break-up

abstract: The breakup of a jet of a viscous fluid with viscosity $\mu _{1}$ immersed into another viscous fluid with viscosity $\mu _{2}$ is considered in the limit when the viscosity ratio $\lambda =\mu _{1}/\mu _{2}$ is close to zero. We show that, in this limit, a transition from ordinary continuous selfsimilarity to discrete selfsimilarity takes place as $\lambda $decreases. The result being that instead of a single point breakup, the rupture of the inner jet occurs through the appearance of an infinite sequence of filaments of decreasing size that will eventually produce infinite sequences of bubbles of the inner fluid inside the outer fluid. The transition can be understood as the result of a Hopf bifurcation in the system of equations modelling the physical problem.


title: Wasserstein Distances for Stereo Disparity Estimation

abstract: Existing approaches to depth or disparity estimation output a distribution over a set of pre-defined discrete values. This leads to inaccurate results when the true depth or disparity does not match any of these values. The fact that this distribution is usually learned indirectly through a regression loss causes further problems in ambiguous regions around object boundaries. We address these issues using a new neural network architecture that is capable of outputting arbitrary depth values, and a new loss function that is derived from the Wasserstein distance between the true and the predicted distributions. We validate our approach on a variety of tasks, including stereo disparity and depth estimation, and the downstream 3D object detection. Our approach drastically reduces the error in ambiguous regions, especially around object boundaries that greatly affect the localization of objects in 3D, achieving the state-of-the-art in 3D object detection for autonomous driving. Our code will be available at https://github.com/Div99/W-Stereo-Disp.


title: Search for high energy $\gamma$-rays from the direction of the candidate electromagnetic counterpart to the binary black hole merger gravitational-wave event S190521g

abstract: The gravitational-wave event S190521g -- a likely binary black hole merger in the accretion disk of an active galactic nucleus -- was accompanied by an optical counterpart. Such dense environments around luminous energy release regions are favourable for high energy $\gamma$-ray production. We report on a search for high energy $\gamma$-rays from the direction of the candidate electromagnetic counterpart to the S190521g event using publicly-available data of the Fermi-LAT space $\gamma$-ray telescope. No significant signal was found. We present upper limits on the spectral energy distribution of the source in the 100 MeV - 300 GeV energy range. We discuss the importance of studying S190521g-like transients in the context of cosmic ray acceleration, $\gamma$-ray and neutrino production in such sources.


title: Smooth Formulation for Isothermal Compositional Simulation with Improved Nonlinear Convergence

abstract: Compositional simulation is challenging, because of highly nonlinear couplings between multi-component flow in porous media with thermodynamic phase behavior. The coupled nonlinear system is commonly solved by the fully-implicit scheme. Various compositional formulations have been proposed. However, severe convergence issues of Newton solvers can arise under the conventional formulations. Crossing phase boundaries produces kinks in discretized equations, and subsequently causing oscillations or even divergence of Newton iterations. The objective of this work is to develop a smooth formulation that removes all the property switches and discontinuities associated with phase changes. We show that it can be very difficult and costly to smooth the conservation equations directly. Therefore, we first reformulate the coupled system, so that the discontinuities are transferred to the phase equilibrium model. In this way a single and concise non-smooth equation is achieved and then a smoothing approximation can be made. The new formulation with a smoothing parameter provides smooth transitions of variables across all the phase regimes. In addition, we employ a continuation method where the solution progressively evolves toward the target system. We evaluate the efficiency of the new smooth formulation and the continuation method using several complex problems. Compared to the standard natural formulation, the developed formulation and method exhibit superior nonlinear convergence behaviors. The continuation method leads to smooth and stable iterative performance, with a negligible impact on solution accuracy. Moreover, it works robustly for a wide range of flow conditions without parameter tuning.


title: Arithmetic properties of the sum of divisors

abstract: The divisor function $\sigma(n)$ denotes the sum of the divisors of the positive integer $n$. For a prime $p$ and $m \in \mathbb{N}$, the $p$-adic valuation of $m$ is the highest power of $p$ which divides $m$. Formulas for $\nu_{p}(\sigma(n))$ are established. For $p=2$, these involve only the odd primes dividing $n$. These expressions are used to establish the bound $\nu_{2}(\sigma(n)) \leq \lceil\log_{2}(n) \rceil$, with equality if and only if $n$ is the product of distinct Mersenne primes, and for an odd prime $p$, the bound is $\nu_{p}(\sigma(n)) \leq \lceil \log_{p}(n) \rceil$, with equality related to solutions of the Ljunggren-Nagell diophantine equation.


title: Non-Gaussian tail in the force distribution: A hallmark of correlated disorder in the host media of elastic objects

abstract: Inferring the nature of disorder in the media where elastic objects are nucleated is of crucial importance for many applications but remains a challenging basic-science problem. Here we propose a method to discern whether weak-point or strong-correlated disorder dominates based on characterizing the distribution of the interaction forces between objects mapped in large fields-of-view. We illustrate our proposal with the case-study system of vortex structures nucleated in type-II superconductors with different pinning landscapes. Interaction force distributions are computed from individual vortex positions imaged in thousands-vortices fields-of-view in a two-orders-of-magnitude-wide vortex-density range. Vortex structures nucleated in point-disordered media present Gaussian distributions of the interaction force components. In contrast, if the media have dilute and randomly-distributed correlated disorder, these distributions present non-Gaussian algebraically-decaying tails for large force magnitudes. We propose that detecting this deviation from the Gaussian behavior is a fingerprint of strong disorder, in our case originated from a dilute distribution of correlated pinning centers.


title: Mapping images into ordinal networks

abstract: An increasing abstraction has marked some recent investigations in network science. Examples include the development of algorithms that map time series data into networks whose vertices and edges can have different interpretations, beyond the classical idea of parts and interactions of a complex system. These approaches have proven useful for dealing with the growing complexity and volume of diverse data sets. However, the use of such algorithms is mostly limited to one-dimension data, and there has been little effort towards extending these methods to higher-dimensional data such as images. Here we propose a generalization for the ordinal network algorithm for mapping images into networks. We investigate the emergence of connectivity constraints inherited from the symbolization process used for defining the network nodes and links, which in turn allows us to derive the exact structure of ordinal networks obtained from random images. We illustrate the use of this new algorithm in a series of applications involving randomization of periodic ornaments, images generated by two-dimensional fractional Brownian motion and the Ising model, and a data set of natural textures. These examples show that measures obtained from ordinal networks (such as average shortest path and global node entropy) extract important image properties related to roughness and symmetry, are robust against noise, and can achieve higher accuracy than traditional texture descriptors extracted from gray-level co-occurrence matrices in simple image classification tasks.


title: Electron-phonon coupling and superconductivity in the doped topological-crystalline insulator (Pb$_{0.5}$Sn$_{0.5}$)$_{1-x}$In$_x$Te

abstract: We present a neutron scattering study of phonons in single crystals of (Pb$_{0.5}$Sn$_{0.5}$)$_{1-x}$In$_x$Te with $x=0$ (metallic, but nonsuperconducting) and $x=0.2$ (nonmetallic normal state, but superconducting). We map the phonon dispersions (more completely for $x=0$) and find general consistency with theoretical calculations, except for the transverse and longitudinal optical (TO and LO) modes at the Brillouin zone center. At low temperature, both modes are strongly damped but sit at a finite energy ($\sim4$ meV in both samples), shifting to higher energy at room temperature. These modes are soft due to a proximate structural instability driven by the sensitivity of Pb-Te and Sn-Te $p$-orbital hybridization to off-center displacements of the metal atoms. The impact of the soft optical modes on the low-energy acoustic modes is inferred from the low thermal conductivity, especially at low temperature. Given that the strongest electron-phonon coupling is predicted for the LO mode, which should be similar for both studied compositions, it is intriguing that only the In-doped crystal is superconducting. In addition, we observe elastic diffuse (Huang) scattering that is qualitatively explained by the difference in Pb-Te and Sn-Te bond lengths within the lattice of randomly distributed Pb and Sn sites. We also confirm the presence of anomalous diffuse low-energy atomic vibrations that we speculatively attribute to local fluctuations of individual Pb atoms between off-center sites.


title: Neural Subgraph Matching

abstract: Subgraph matching is the problem of determining the presence and location(s) of a given query graph in a large target graph. Despite being an NP-complete problem, the subgraph matching problem is crucial in domains ranging from network science and database systems to biochemistry and cognitive science. However, existing techniques based on combinatorial matching and integer programming cannot handle matching problems with both large target and query graphs. Here we propose NeuroMatch, an accurate, efficient, and robust neural approach to subgraph matching. NeuroMatch decomposes query and target graphs into small subgraphs and embeds them using graph neural networks. Trained to capture geometric constraints corresponding to subgraph relations, NeuroMatch then efficiently performs subgraph matching directly in the embedding space. Experiments demonstrate NeuroMatch is 100x faster than existing combinatorial approaches and 18% more accurate than existing approximate subgraph matching methods.


title: Feasibility Analysis for the Problem of Active Near Field/Far Field Acoustic Pattern Synthesis in Free Space and Shallow Water Environments

abstract: In this paper, a detailed sensitivity analysis of the active manipulation scheme for scalar Helmholtz fields proposed in our previous works, in both free space and constant-depth homogeneous ocean environments, is presented. We apply the method of moments (MoM) together with Tikhonov regularization with Morozov discrepancy principle to investigate the effects of problem parameters variations on the accuracy and feasibility of the proposed active field control strategy. We discuss the feasibility of the active scheme (power budget and control accuracy) as a function of the frequency, the distance between the control region and the active source, the mutual distance between the control regions, and the size of the control region. Process error is considered as well to investigate the possibility of an accurate active control in the presence of manufacturing noise. The numerical simulations show the accuracy of the active field control scheme and indicate some challenges and limitations for its physical implementations.


title: An analogue of Amitsur's property for the ring of pseudo-differential operators

abstract: Let R be a ring with a derivation \delta. In this paper, we prove that an analogue of Amitsur's property holds for left T-nilpotent radideals of pseudo-differential operator rings R((x^{-1}; \delta)), where R is a delta-compatible ring. As a direct consequence of this fact, we obtain an alternative characterization of the prime radical of R((x^{-1}; \delta)).


title: Kortweg de-Vries solitons on electrified liquid jets

abstract: The propagation of axisymmetric waves on the surface of a liquid jet under the action of a radial electric field is considered. The jet is assumed to be inviscid and perfectly conducting, and a field is set up by placing the jet concentrically inside a perfectly cylindrical tube whose wall is maintained at a constant potential. A nontrivial interaction arises between the hydrodynamics and the electric field in the annulus, resulting in the formation of electrocapillary waves. The main objective of the present study is to describe nonlinear aspects of such axisymmetric waves in the weakly nonlinear regime which is valid for long waves relative to the undisturbed jet radius. This is found to be possible if two conditions hold: the outer electrode radius is not too small, and the applied electric field is sufficiently strong. Under these conditions long waves are shown to be dispersive and a weakly nonlinear theory can be developed to describe the evolution of the disturbances. The canonical system that arises is the Kortweg de-Vries equation with coefficients that vary as the electric field and the electrode radius are varied. Interestingly, the coefficient of the highest order third derivative term does not change sign and remains strictly positive, whereas the coefficient $\alpha$ of the nonlinear term can change sign for certain values of the parameters. This finding implies that solitary electrocapillary waves are possible; there are waves of elevation for $\alpha>0$ and of depression for $\alpha<0$. Regions in parameter space are identified where such waves are found.


title: Domain Adaptation for Ultrasound Beamforming

abstract: Ultrasound B-Mode images are created from data obtained from each element in the transducer array in a process called beamforming. The beamforming goal is to enhance signals from specified spatial locations, while reducing signal from all other locations. On clinical systems, beamforming is accomplished with the delay-and-sum (DAS) algorithm. DAS is efficient but fails in patients with high noise levels, so various adaptive beamformers have been proposed. Recently, deep learning methods have been developed for this task. With deep learning methods, beamforming is typically framed as a regression problem, where clean, ground-truth data is known, and usually simulated. For in vivo data, however, it is extremely difficult to collect ground truth information, and deep networks trained on simulated data underperform when applied to in vivo data, due to domain shift between simulated and in vivo data. In this work, we show how to correct for domain shift by learning deep network beamformers that leverage both simulated data, and unlabeled in vivo data, via a novel domain adaption scheme. A challenge in our scenario is that domain shift exists both for noisy input, and clean output. We address this challenge by extending cycle-consistent generative adversarial networks, where we leverage maps between synthetic simulation and real in vivo domains to ensure that the learned beamformers capture the distribution of both noisy and clean in vivo data. We obtain consistent in vivo image quality improvements compared to existing beamforming techniques, when applying our approach to simulated anechoic cysts and in vivo liver data.


title: FleCSPH: The Next Generation FleCSIble Parallel Computational Infrastructure for Smoothed Particle Hydrodynamics

abstract: FleCSPH is a smoothed particle hydrodynamics simulation tool, based on the compile-time configurable framework FleCSI. The asynchronous distributed tree topology combined with a fast multipole method allows FleCSPH to efficiently compute hydrodynamics and long range particle-particle interactions. FleCSPH provides initial data generators, particle relaxation techniques, and standard evolution drivers, which can be easily modified and extended to user-specific setups. Data input/output uses the H5part format, compatible with modern visualization software.


title: Text Recognition -- Real World Data and Where to Find Them

abstract: We present a method for exploiting weakly annotated images to improve text extraction pipelines. The approach uses an arbitrary end-to-end text recognition system to obtain text region proposals and their, possibly erroneous, transcriptions. The proposed method includes matching of imprecise transcription to weak annotations and edit distance guided neighbourhood search. It produces nearly error-free, localised instances of scene text, which we treat as "pseudo ground truth" (PGT). We apply the method to two weakly-annotated datasets. Training with the extracted PGT consistently improves the accuracy of a state of the art recognition model, by 3.7~\% on average, across different benchmark datasets (image domains) and 24.5~\% on one of the weakly annotated datasets.


title: Eventual regularization for the 3D Muskat problem: Lipschitz for finite time implies global existence

abstract: We show that for any fixed Lipschitz constant $L$, there is a time $T^*<\infty$ depending only on $L$ such that if $f:[0,T^*]\times \mathbb{R}^{2}\to [0,1]$ is a classical solution of the stable Muskat problem with $||\nabla_x f||_{L^\infty}\leq L$, then $f$ extends uniquely to a classical solution defined for all time. Our proof is short, quantitative, and explicit.


title: Cost-sensitive Multi-class AdaBoost for Understanding Driving Behavior with Telematics

abstract: Powered with telematics technology, insurers can now capture a wide range of data, such as distance traveled, how drivers brake, accelerate or make turns, and travel frequency each day of the week, to better decode driver's behavior. Such additional information helps insurers improve risk assessments for usage-based insurance (UBI), an increasingly popular industry innovation. In this article, we explore how to integrate telematics information to better predict claims frequency. For motor insurance during a policy year, we typically observe a large proportion of drivers with zero claims, a less proportion with exactly one claim, and far lesser with two or more claims. We introduce the use of a cost-sensitive multi-class adaptive boosting (AdaBoost) algorithm, which we call SAMME.C2, to handle such imbalances. To calibrate SAMME.C2 algorithm, we use empirical data collected from a telematics program in Canada and we find improved assessment of driving behavior with telematics relative to traditional risk variables. We demonstrate our algorithm can outperform other models that can handle class imbalances: SAMME, SAMME with SMOTE, RUSBoost, and SMOTEBoost. The sampled data on telematics were observations during 2013-2016 for which 50,301 are used for training and another 21,574 for testing. Broadly speaking, the additional information derived from vehicle telematics helps refine risk classification of drivers of UBI.


title: Statistical Physics of Epidemic on Network Predictions for SARS-CoV-2 Parameters

abstract: The SARS-CoV-2 pandemic has necessitated mitigation efforts around the world. We use only reported deaths in the two weeks after the first death to determine infection parameters, in order to make predictions of hidden variables such as the time dependence of the number of infections. Early deaths are sporadic and discrete so the use of network models of epidemic spread is imperative, with the network itself a crucial random variable. Location-specific population age distributions and population densities must be taken into account when attempting to fit these events with parametrized models. These characteristics render naive Bayesian model comparison impractical as the networks have to be large enough to avoid finite-size effects. We reformulated this problem as the statistical physics of independent location-specific `balls' attached to every model in a six-dimensional lattice of 56448 parametrized models by elastic springs, with model-specific `spring constants' determined by the stochasticity of network epidemic simulations for that model. The distribution of balls then determines all Bayes posterior expectations. Important characteristics of the contagion are determinable: the fraction of infected patients that die ($0.017\pm 0.009$), the expected period an infected person is contagious ($22 \pm 6$ days) and the expected time between the first infection and the first death ($25 \pm 8$ days) in the US. The rate of exponential increase in the number of infected individuals is $0.18\pm 0.03$ per day, corresponding to 65 million infected individuals in one hundred days from a single initial infection, which fell to 166000 with even imperfect social distancing effectuated two weeks after the first recorded death. The fraction of compliant socially-distancing individuals matters less than their fraction of social contact reduction for altering the cumulative number of infections.


title: Natural Emergence of Heterogeneous Strategies in Artificially Intelligent Competitive Teams

abstract: Multi agent strategies in mixed cooperative-competitive environments can be hard to craft by hand because each agent needs to coordinate with its teammates while competing with its opponents. Learning based algorithms are appealing but many scenarios require heterogeneous agent behavior for the team's success and this increases the complexity of the learning algorithm. In this work, we develop a competitive multi agent environment called FortAttack in which two teams compete against each other. We corroborate that modeling agents with Graph Neural Networks and training them with Reinforcement Learning leads to the evolution of increasingly complex strategies for each team. We observe a natural emergence of heterogeneous behavior amongst homogeneous agents when such behavior can lead to the team's success. Such heterogeneous behavior from homogeneous agents is appealing because any agent can replace the role of another agent at test time. Finally, we propose ensemble training, in which we utilize the evolved opponent strategies to train a single policy for friendly agents.


title: Resistance distance, Kirchhoff index, and Kemeny's constant in flower graphs

abstract: We obtain a general formula for the resistance distance (or effective resistance) between any pair of nodes in a general family of graphs which we call flower graphs. Flower graphs are obtained from identifying nodes of multiple copies of a given base graph in a cyclic way. We apply our general formula to two specific families of flower graphs, where the base graph is either a complete graph or a cycle. We also obtain bounds on the Kirchhoff index and Kemeny's constant of general flower graphs using our formula for resistance. For flower graphs whose base graph is a complete graph or a cycle, we obtain exact, closed form expressions for the Kirchhoff index and Kemeny's constant.


title: Mass spectrometric analyses of high performance polymers to assess radiopurity as ultra low background materials for rare event physics detectors

abstract: The mass concentrations of Th-232 and U-238 in several strong, unfilled, high performance polymers are reported as measures of their radiopurity. Highly radiopure polymers are required as dielectric materials in the construction of rare event physics detectors, in order to minimize background signals arising from the detector materials themselves. New data are reported for carefully sourced samples of polyetheretherketone (PEEK), polyetherketoneketone (PEKK), a polyamide imide (PAI, branded as Torlon) and polybenzimidazole (PBI). Data for solid polyetherimide (PEI) are also discussed and new data for PEI in the form of a commercial 3D printing filament stock are reported. Strong high performance polymers PEKK, PBI, PAI and PEI were found with levels for Th-232 and U-238 that are below one milliBq/kg, including the PEI 3D printing filament. Specifically, for Th-232 and U-238 respectively, in microBq/kg, we found values of 149 and 184 for Arkema Kepstan PEKK 6002 flake; 69 and 2250 for Solvay Ketaspire PEEK flake; 346 and 291 for PBI Performance Products low metals grade 100 mesh PBI powder; 66 and 105 for Drake Plastics PAI Torlon T-4200 pellets; 401 and 285 for Drake Plastics cured PAI rod; 32 and 41 for Ensinger PEI Ultem 1000 solid; and 15 and 85 microBq/kg, for ThermaX PEI Ultem 1010 filament material. These results were all obtained using a novel dry ashing method in crucibles constructed of ultra low background (ULB) electroformed copper. Samples were spiked with T-229h and U-233 as internal standards prior to ashing and determinations were made by inductively coupled plasma mass spectrometry (ICP-MS). Radiopurity is displayed graphically relative to numerical measures of mechanical strength for these and several other polymers.


title: Multiwavelength Analysis of the Variability of the Blazar 3C 273

abstract: We present multiwavelength light curves and polarimetric data of the Flat Spectrum Radio Quasar 3C 273 over 8 years. The wavelength range of our data set extends from radio to gamma-rays. We found that the optical emission in this source is dominated by the accretion disk during the entire time-frame of study. We additionally find that in contrast with the observed behaviour in other blazars, 3C 273 does not show a correlation between the gamma-ray spectral index and the gamma-ray luminosity. Finally, we identified an anti-correlation between the 15 GHz and V-band light curves for the time-range $JD_{245} = 4860 - 5760$, which we speculate is the consequence of the inner part of the accretion disk falling into the black hole, followed by the ejection of a component into the jet.


title: An Evaluation of Two Commercial Deep Learning-Based Information Retrieval Systems for COVID-19 Literature

abstract: The COVID-19 pandemic has resulted in a tremendous need for access to the latest scientific information, primarily through the use of text mining and search tools. This has led to both corpora for biomedical articles related to COVID-19 (such as the CORD-19 corpus (Wang et al., 2020)) as well as search engines to query such data. While most research in search engines is performed in the academic field of information retrieval (IR), most academic search engines$\unicode{x2013}$though rigorously evaluated$\unicode{x2013}$are sparsely utilized, while major commercial web search engines (e.g., Google, Bing) dominate. This relates to COVID-19 because it can be expected that commercial search engines deployed for the pandemic will gain much higher traction than those produced in academic labs, and thus leads to questions about the empirical performance of these search tools. This paper seeks to empirically evaluate two such commercial search engines for COVID-19, produced by Google and Amazon, in comparison to the more academic prototypes evaluated in the context of the TREC-COVID track (Roberts et al., 2020). We performed several steps to reduce bias in the available manual judgments in order to ensure a fair comparison of the two systems with those submitted to TREC-COVID. We find that the top-performing system from TREC-COVID on bpref metric performed the best among the different systems evaluated in this study on all the metrics. This has implications for developing biomedical retrieval systems for future health crises as well as trust in popular health search engines.


title: Multi-image Super Resolution of Remotely Sensed Images using Residual Feature Attention Deep Neural Networks

abstract: Convolutional Neural Networks (CNNs) have been consistently proved state-of-the-art results in image Super-Resolution (SR), representing an exceptional opportunity for the remote sensing field to extract further information and knowledge from captured data. However, most of the works published in the literature have been focusing on the Single-Image Super-Resolution problem so far. At present, satellite based remote sensing platforms offer huge data availability with high temporal resolution and low spatial resolution. In this context, the presented research proposes a novel residual attention model (RAMS) that efficiently tackles the multi-image super-resolution task, simultaneously exploiting spatial and temporal correlations to combine multiple images. We introduce the mechanism of visual feature attention with 3D convolutions in order to obtain an aware data fusion and information extraction of the multiple low-resolution images, transcending limitations of the local region of convolutional operations. Moreover, having multiple inputs with the same scene, our representation learning network makes extensive use of nestled residual connections to let flow redundant low-frequency signals and focus the computation on more important high-frequency components. Extensive experimentation and evaluations against other available solutions, either for single or multi-image super-resolution, have demonstrated that the proposed deep learning-based solution can be considered state-of-the-art for Multi-Image Super-Resolution for remote sensing applications.


title: Breakdown signatures of the phenomenological Lindblad master equation in the strong optomechanical coupling regime

abstract: The Lindblad form of the master equation has proven to be one of the most convenient ways to describe the impact of an environment interacting with a quantum system of interest. For single systems the jump operators characterizing these interactions usually take simple forms with a clear interpretation. However, for coupled systems these operators take significantly different forms and the full dynamics cannot be described by jump operators acting on the individual subsystems only. In this work, we investigate the differences between a common phenomenological model for the master equation and the more rigorous dressed-state master equation for optomechanical systems. We provide an analytical method to obtain the absorption spectrum of the system for both models and show the breakdown of the phenomenological model in both the bad cavity and the ultra-strong coupling limit. We present a careful discussion of the indirect dephasing of the optical cavity in both models and its role in the differences of their predicted absorption spectra. Our work provides a simple experimental test to determine whether the simpler phenomenological model can be used to describe the system and is a step forward toward a better understanding of the role of the coupling between subsystems for open-quantum-system dynamics.


title: Cycle-StarNet: Bridging the gap between theory and data by leveraging large datasets

abstract: The advancements in stellar spectroscopy data acquisition have made it necessary to accomplish similar improvements in efficient data analysis techniques. Current automated methods for analyzing spectra are either (a) data-driven, which requires prior knowledge of stellar parameters and elemental abundances, or (b) based on theoretical synthetic models that are susceptible to the gap between theory and practice. In this study, we present a hybrid generative domain adaptation method that turns simulated stellar spectra into realistic spectra by applying unsupervised learning to large spectroscopic surveys. We apply our technique to the APOGEE H-band spectra at R=22,500 and the Kurucz synthetic models. As a proof of concept, two case studies are presented. The first of which is the calibration of synthetic data to become consistent with observations. To accomplish this, synthetic models are morphed into spectra that resemble observations, thereby reducing the gap between theory and observations. Fitting the observed spectra shows an improved average reduced $\chi_R^2$ from 1.97 to 1.22, along with a reduced mean residual from 0.16 to -0.01 in normalized flux. The second case study is the identification of the elemental source of missing spectral lines in the synthetic modelling. A mock dataset is used to show that absorption lines can be recovered when they are absent in one of the domains. This method can be applied to other fields, which use large data sets and are currently limited by modelling accuracy. The code used in this study is made publicly available on github.

categories: astro-ph.SR astro-ph.GA astro-ph.IM physics.data-an stat.ML


title: The DG Products of Peeva and Srinivasan Coincide

abstract: Consider the ideal $(x_1 , \dotsc , x_n)^d \subseteq k[x_1 , \dotsc , x_n]$, where $k$ is any field. This ideal can be resolved by both the $L$-complexes of Buchsbaum and Eisenbud, and the Eliahou-Kervaire resolution. Both of these complexes admit the structure of an associative DG algebra, and it is a question of Peeva as to whether these DG structures coincide in general. In this paper, we construct an isomorphism of complexes between the aforementioned complexes that is also an isomorphism of algebras with their respective products, thus giving an affirmative answer to Peeva's question.


title: Two dimensional electron gas in a non-Euclidean space

abstract: A charged particle in the presence of a magnetic field is studied in the position dependent operator formalism. Instead of a quantum harmonic oscillator, the solution of the resulting Schr\"odinger-like equation is the one for the Morse oscillator. The anharmonicity that shows up naturally from the theory is analogous to the corrections introduced by relativistic ones.The degeneracy of the spin up and down levels is lifted due to the non-Euclidean space. It is shown that the Fermi energy and the total ground state energy of the 2DEG is also modified.


title: Interpreting Stellar Spectra with Unsupervised Domain Adaptation

abstract: We discuss how to achieve mapping from large sets of imperfect simulations and observational data with unsupervised domain adaptation. Under the hypothesis that simulated and observed data distributions share a common underlying representation, we show how it is possible to transfer between simulated and observed domains. Driven by an application to interpret stellar spectroscopic sky surveys, we construct the domain transfer pipeline from two adversarial autoencoders on each domains with a disentangling latent space, and a cycle-consistency constraint. We then construct a differentiable pipeline from physical stellar parameters to realistic observed spectra, aided by a supplementary generative surrogate physics emulator network. We further exemplify the potential of the method on the reconstructed spectra quality and to discover new spectral features associated to elemental abundances.


title: Examining COVID-19 Forecasting using Spatio-Temporal Graph Neural Networks

abstract: In this work, we examine a novel forecasting approach for COVID-19 case prediction that uses Graph Neural Networks and mobility data. In contrast to existing time series forecasting models, the proposed approach learns from a single large-scale spatio-temporal graph, where nodes represent the region-level human mobility, spatial edges represent the human mobility based inter-region connectivity, and temporal edges represent node features through time. We evaluate this approach on the US county level COVID-19 dataset, and demonstrate that the rich spatial and temporal information leveraged by the graph neural network allows the model to learn complex dynamics. We show a 6% reduction of RMSLE and an absolute Pearson Correlation improvement from 0.9978 to 0.998 compared to the best performing baseline models. This novel source of information combined with graph based deep learning approaches can be a powerful tool to understand the spread and evolution of COVID-19. We encourage others to further develop a novel modeling paradigm for infectious disease based on GNNs and high resolution mobility data.


title: Efficient Conformal Prediction via Cascaded Inference with Expanded Admission

abstract: In this paper, we present a novel approach for conformal prediction (CP), in which we aim to identify a set of promising prediction candidates -- in place of a single prediction. This set is guaranteed to contain a correct answer with high probability, and is well-suited for many open-ended classification tasks. In the standard CP paradigm, the predicted set can often be unusably large and also costly to obtain. This is particularly pervasive in settings where the correct answer is not unique, and the number of total possible answers is high. We first expand the CP correctness criterion to allow for additional, inferred "admissible" answers, which can substantially reduce the size of the predicted set while still providing valid performance guarantees. Second, we amortize costs by conformalizing prediction cascades, in which we aggressively prune implausible labels early on by using progressively stronger classifiers -- again, while still providing valid performance guarantees. We demonstrate the empirical effectiveness of our approach for multiple applications in natural language processing and computational chemistry for drug discovery.


title: Arid or Cloudy: Characterizing the Atmosphere of the super-Earth 55 Cancri e using High-Resolution Spectroscopy

abstract: The nearby super-Earth 55 Cnc e orbits a bright (V = 5.95 mag) star with a period of ~ 18 hours and a mass of ~ 8 Earth masses. Its atmosphere may be water-rich and have a large scale-height, though attempts to characterize it have yielded ambiguous results. Here we present a sensitive search for water and TiO in its atmosphere at high spectral resolution using the Gemini North telescope and the GRACES spectrograph. We combine observations with previous observations from Subaru and CFHT, improving the constraints on the presence of water vapor. We adopt parametric models with an updated planet radius based on recent measurements, and use a cross-correlation technique to maximize sensitivity. Our results are consistent with atmospheres that are cloudy or contain minimal amounts of water and TiO. Using these parametric models, we rule out a water-rich atmosphere (VMR >= 0.1%) with a mean molecular weight of <= 15 g/mol at a 3 sigma confidence level, improving on the previous limit by a significant margin. For TiO, we rule out a mean molecular weight of <= 5 g/mol with a 3 sigma confidence level for a VMR greater than 10^-8; for a VMR of greater than 10^-7, the limit rises to a mean molecular weight of <= 10 g/mol. We can rule out low mean-molecular-weight chemical equilibrium models both including and excluding TiO/VO at very high confidence levels (> 10 sigma). Overall, our results are consistent with an atmosphere with a high mean molecular weight and/or clouds, or no atmosphere.


title: Ruelle resonances from cohomological equations

abstract: These notes are based on lectures given by the author at the Summer School on Teichm\"uller dynamics, mapping class groups and applications in Grenoble, France, in June 2018 and at the Oberwolfach Seminar on Anisotropic Spaces and their Applications to Hyperbolic and Parabolic Systems in June 2019. We derive results about the so-called Ruelle resonances and the asymptotics of correlations for several classes of systems from known results on cohomological equations and invariant distributions for the respective unstable vector fields. In particular, we consider pseudo-Anosov diffeomorphisms on surfaces of higher genus, for horocycle flows on surfaces of constant negative curvature and for partially hyperbolic automorphisms of Heisenberg 3-dimensional nilmanfolds. Ruelles resonances for pseudo-Anosov maps with applications to the cohomological equation for their unstable translation flows was recently studied in depth by F. Faure, S. Gou\"ezel and E. Lanneau [FGL] by methods based on the analysis of the transfer operator of the pseudo-Anosov map. Ruelle resonances for geodesic flows on hyperbolic compact manifolds of any dimension and of partially hyperbolic automorphisms of Heisenberg 3-dimensional nilmanfolds are studied by general results of Dyatlov, Faure and Guillarmou [DFG] and Faure and Tsujii [FT15] based on methods of semi-classical analysis. These works do not derive results on cohomological equations for unstable flows or horospherical foliations of these systems.


title: Multi-Fidelity Bayesian Optimization via Deep Neural Networks

abstract: Bayesian optimization (BO) is a popular framework to optimize black-box functions. In many applications, the objective function can be evaluated at multiple fidelities to enable a trade-off between the cost and accuracy. To reduce the optimization cost, many multi-fidelity BO methods have been proposed. Despite their success, these methods either ignore or over-simplify the strong, complex correlations across the fidelities, and hence can be inefficient in estimating the objective function. To address this issue, we propose Deep Neural Network Multi-Fidelity Bayesian Optimization (DNN-MFBO) that can flexibly capture all kinds of complicated relationships between the fidelities to improve the objective function estimation and hence the optimization performance. We use sequential, fidelity-wise Gauss-Hermite quadrature and moment-matching to fulfill a mutual information-based acquisition function, which is computationally tractable and efficient. We show the advantages of our method in both synthetic benchmark datasets and real-world applications in engineering design.


title: Toward $e^+e^-\to \pi^+\pi^-$ annihilation inspired by higher $\rho$ mesonic states around 2.2 GeV

abstract: Very recently, the {\it BABAR} Collaboration indicated that there exist an explicit enhancement structure near 2.2 GeV when focusing on the $e^+e^-\to\pi^+\pi^-$ process again, which inspires our interest in studying the production of higher $\rho$ mesonic states. Since the branching ratio of $\pi^+\pi^-$ channel of $D-$wave $\rho$ states are much smaller than $S-$wave states, we choose $\rho(1900)$ and $\rho(2150)$ as the intermediate states in $e^+e^-\to\pi^+\pi^-$, where $\rho(1900)$ and $\rho(2150)$ are treated as $\rho(3S)$ and $\rho(4S)$ states, respectively. Our result indicates that the $BABAR$'s data of $e^+e^-\to\pi^+\pi^-$ around 2 GeV can be depicted well, which shows that this enhancement structure near 2.2 GeV existing in $e^+e^-\to\pi^+\pi^-$ can be due to the contribution from two $\rho$ mesons, $\rho(1900)$ and $\rho(2150)$. Additionally, this conclusion can be enforced by the consistence of the extracted values of $\Gamma_{e^+e^-}\mathcal{B}(\pi^+\pi^-)$ of $\rho(1900)$ and $\rho(2150)$ in the whole fitting processes and the corresponding theoretical calculations. The present study of $e^+e^-\to\pi^+\pi^-$ data may provide valuable information to establish the $\rho$ meson family.


title: X-ray Photon-Counting Data Correction through Deep Learning

abstract: X-ray photon-counting detectors (PCDs) are drawing an increasing attention in recent years due to their low noise and energy discrimination capabilities. The energy/spectral dimension associated with PCDs potentially brings great benefits such as for material decomposition, beam hardening and metal artifact reduction, as well as low-dose CT imaging. However, X-ray PCDs are currently limited by several technical issues, particularly charge splitting (including charge sharing and K-shell fluorescence re-absorption or escaping) and pulse pile-up effects which distort the energy spectrum and compromise the data quality. Correction of raw PCD measurements with hardware improvement and analytic modeling is rather expensive and complicated. Hence, here we proposed a deep neural network based PCD data correction approach which directly maps imperfect data to the ideal data in the supervised learning mode. In this work, we first establish a complete simulation model incorporating the charge splitting and pulse pile-up effects. The simulated PCD data and the ground truth counterparts are then fed to a specially designed deep adversarial network for PCD data correction. Next, the trained network is used to correct separately generated PCD data. The test results demonstrate that the trained network successfully recovers the ideal spectrum from the distorted measurement within $\pm6\%$ relative error. Significant data and image fidelity improvements are clearly observed in both projection and reconstruction domains.


title: Using electric fields to induce patterning in leaky dielectric fluids in a rod-annular geometry

abstract: The stability and axisymmetric deformation of two immiscible, viscous, perfect or leaky dielectric fluids confined in the annulus between two concentric cylinders are studied in the presence of radial electric fields. The fields are set up by imposing a constant voltage potential difference between the inner and outer cylinders. We derive a set of equations for the interface in the long-wavelength approximation which retains the essential physics of the system and allows for interfacial deformations to be as large as the annular gap hence accounting for possible touchdown at the inner or outer electrode. The effects of the electric parameters are evaluated initially by performing a linear stability analysis which shows excellent agreement with the linear theory of the full axisymmetric problem in the appropriate long wavelength regime. The nonlinear interfacial dynamics are investigated by carrying out direct numerical simulations of the derived long wave models, both in the absence and presence of electric fields. For non-electrified thin layer flows (i.e. one of the layers thin relative to the other) the long-time dynamics agree with the lubrication approximation results found in literature. When the liquid layers have comparable thickness our results demonstrate the existence of both finite time and infinite time singularities (asymptotic touching solutions) in the system. It is shown that a two-side touching solution is possible for both the non-electrified and perfect dielectric cases, while only one-side touching is found in the case of leaky dielectric liquids, where the flattened interface shape resembles the pattern solutions found in literature. Meanwhile the finite-time singular solution agrees qualitatively with the experiments of~\cite{Reynolds1965}.


title: Multi-Armed Bandits with Local Differential Privacy

abstract: This paper investigates the problem of regret minimization for multi-armed bandit (MAB) problems with local differential privacy (LDP) guarantee. In stochastic bandit systems, the rewards may refer to the users' activities, which may involve private information and the users may not want the agent to know. However, in many cases, the agent needs to know these activities to provide better services such as recommendations and news feeds. To handle this dilemma, we adopt differential privacy and study the regret upper and lower bounds for MAB algorithms with a given LDP guarantee. In this paper, we prove a lower bound and propose algorithms whose regret upper bounds match the lower bound up to constant factors. Numerical experiments also confirm our conclusions.


title: Using classical Cepheids to study the far side of the Milky Way disk: I. Spectroscopic classification and the metallicity gradient

abstract: The structure, kinematics, and chemical composition of the far side of the Milky Way disk, beyond the bulge, are still to be revealed. Classical Cepheids (CCs) are young and luminous standard candles. We aim to use a well-characterized sample of these variable stars to study the present time properties of the far side of the Galactic disk. A sample of 45 Cepheid variable star candidates were selected from near infrared time series photometry obtained by the VVV survey. We characterized this sample using high quality near infrared spectra obtained with VLT/X-Shooter, deriving radial velocities and iron abundances for all the sample Cepheids. This allowed us to separate the CCs, which are metal rich and with kinematics consistent with the disk rotation, from type II Cepheids (T2Cs), which are more metal poor and with different kinematics. We estimated individual distances and extinctions using VVV photometry and period-luminosity relations, reporting the characterization of 30 CCs located on the far side of the Galactic disk, plus 8 T2Cs mainly located in the bulge region, of which 10 CCs and 4 T2Cs are new discoveries. This is the first sizeable sample of CCs in this distant region of our Galaxy that has been spectroscopically confirmed. We use their positions, kinematics, and metallicities to confirm that the general properties of the far disk are similar to those of the well-studied disk on the solar side of the Galaxy. In addition, we derive for the first time the radial metallicity gradient on the disk's far side. Considering all the CCs with $R_{\mathrm{GC}} < 17\,\rm{kpc}$, we measure a gradient with a slope of $-0.062 \, \mathrm{dex\, kpc^{-1}}$ and an intercept of $+0.59 \, \rm{dex}$, which is in agreement with previous determinations based on CCs on the near side of the disk.


title: Learning Embeddings for Image Clustering: An Empirical Study of Triplet Loss Approaches

abstract: In this work, we evaluate two different image clustering objectives, k-means clustering and correlation clustering, in the context of Triplet Loss induced feature space embeddings. Specifically, we train a convolutional neural network to learn discriminative features by optimizing two popular versions of the Triplet Loss in order to study their clustering properties under the assumption of noisy labels. Additionally, we propose a new, simple Triplet Loss formulation, which shows desirable properties with respect to formal clustering objectives and outperforms the existing methods. We evaluate all three Triplet loss formulations for K-means and correlation clustering on the CIFAR-10 image classification dataset.


title: Efficiency of learning vs. processing: Towards a normative theory of multitasking

abstract: A striking limitation of human cognition is our inability to execute some tasks simultaneously. Recent work suggests that such limitations can arise from a fundamental tradeoff in network architectures that is driven by the sharing of representations between tasks: sharing promotes quicker learning, at the expense of interference while multitasking. From this perspective, multitasking failures might reflect a preference for learning efficiency over multitasking capability. We explore this hypothesis by formulating an ideal Bayesian agent that maximizes expected reward by learning either shared or separate representations for a task set. We investigate the agent's behavior and show that over a large space of parameters the agent sacrifices long-run optimality (higher multitasking capacity) for short-term reward (faster learning). Furthermore, we construct a general mathematical framework in which rational choices between learning speed and processing efficiency can be examined for a variety of different task environments.


title: MAICRM : A general model for rapid simulation of hot dense plasmas

abstract: We propose a general model, Multi-Average Ion Collisional-Radiative Model (MAICRM), to rapid simulate the ionization and population distributions of hot dense plasmas. In MAICRM, the orbital occupation numbers of ions at the same charge stage are averaged and determined by the excitation and de-excitation processes; the populations of the average ions are determined by the ionization and recombination processes with the fixed orbital average occupation numbers in each ion. The calculated mean ionizations and charge state distributions of MAICRM are in general agreement with the other theoretical and experimental results especially for the mid- and high-density plasmas. Since MAICRM considers more detailed transitions and ionization balances than the average atom model and is faster than DCA/SCA models, this model has the advantage to be combined into hydrodynamic simulations.


title: On a thermostated Kac model with rescaling

abstract: We introduce a global thermostat on Kac's 1D model for the velocities of particles in a space-homogeneous gas subjected to binary collisions, also interacting with a (local) Maxwellian thermostat. The global thermostat rescales the velocities of all the particles, thus restoring the total energy of the system, which leads to an additional drift term in the corresponding nonlinear kinetic equation. We prove ergodicity for this equation, and show that its equilibrium distribution has a density that, depending on the parameters of the model, can exhibit heavy tails, and whose behaviour at the origin can range from being analytic, to being $C^k$, and even to blowing-up. Finally, we prove propagation of chaos for the associated $N$-particle system, with a uniform-in-time rate of order $N^{-\eta}$ in the squared $2$-Wasserstein metric, for an explicit $\eta \in (0, 1/3]$.


title: Free-Boundary Minimal Surfaces of Constant Kahler Angle in Complex Space Forms

abstract: In real space forms, Fraser and Schoen proved that a free-boundary minimal disk in a geodesic ball is totally geodesic. In this note, we consider free-boundary minimal surfaces $\Sigma$ (of any genus) in geodesic balls of complex space forms. In $\mathbb{CP}^2$, $\mathbb{C}^2$ and $\mathbb{CH}^2$, we show that if $\Sigma$ is Lagrangian, then $\Sigma$ is totally geodesic. In $\mathbb{CP}^n$, $\mathbb{C}^n$ and $\mathbb{CH}^n$ for $n \geq 2$, we show that if $\Sigma$ has K\"{a}hler angle $\pi/2$, then $\Sigma$ is superminimal.


title: Neutron Star Extreme Matter Observatory: A kilohertz-band gravitational-wave detector in the global network

abstract: Gravitational waves from coalescing neutron stars encode information about nuclear matter at extreme densities, inaccessible by laboratory experiments. The late inspiral is influenced by the presence of tides, which depend on the neutron star equation of state. Neutron star mergers are expected to often produce rapidly-rotating remnant neutron stars that emit gravitational waves. These will provide clues to the extremely hot post-merger environment. This signature of nuclear matter in gravitational waves contains most information in the 2-4 kHz frequency band, which is outside of the most sensitive band of current detectors. We present the design concept and science case for a neutron star extreme matter observatory (NEMO): a gravitational-wave interferometer optimized to study nuclear physics with merging neutron stars. The concept uses high circulating laser power, quantum squeezing and a detector topology specifically designed to achieve the high-frequency sensitivity necessary to probe nuclear matter using gravitational waves. Above one kHz, the proposed strain sensitivity is comparable to full third-generation detectors at a fraction of the cost. Such sensitivity changes expected event rates for detection of post-merger remnants from approximately one per few decades with two A+ detectors to a few per year, and potentially allows for the first gravitational-wave observations of supernovae, isolated neutron stars, and other exotica.


title: Confounding Ghost Channels and Causality: A New Approach to Causal Information Flows

abstract: Information theory provides a fundamental framework for the quantification of information flows through channels, formally Markov kernels. However, quantities such as mutual information and conditional mutual information do not necessarily reflect the causal nature of such flows. We argue that this is often the result of conditioning based on sigma algebras that are not associated with the given channels. We propose a version of the (conditional) mutual information based on families of sigma algebras that are coupled with the underlying channel. This leads to filtrations which allow us to prove a corresponding causal chain rule as a basic requirement within the presented approach.


title: Electromyogram (EMG) Removal by Adding Sources of EMG (ERASE) -- A novel ICA-based algorithm for removing myoelectric artifacts from EEG -- Part 1

abstract: Electroencephalographic (EEG) recordings are often contaminated by electromyographic (EMG) artifacts, especially when recording during movement. Existing methods to remove EMG artifacts include independent component analysis (ICA), and other high-order statistical methods. However, these methods can not effectively remove most of EMG artifacts. Here, we proposed a modified ICA model for EMG artifacts removal in the EEG, which is called EMG Removal by Adding Sources of EMG (ERASE). In this new approach, additional channels of real EMG from neck and head muscles (reference artifacts) were added as inputs to ICA in order to "force" the most power from EMG artifacts into a few independent components (ICs). The ICs containing EMG artifacts (the "artifact ICs") were identified and rejected using an automated procedure. Simulation results showed ERASE removed EMG artifacts from EEG significantly more effectively than conventional ICA. Subsequently, EEG was collected from 8 healthy participants while they moved their hands to test the realistic efficacy of this approach. Results showed that ERASE successfully removed EMG artifacts (on average, about 75% of EMG artifacts were removed when using real EMGs as reference artifacts) while preserving the expected EEG features related to movement. We also tested the ERASE procedure using simulated EMGs as reference artifacts (about 63% of EMG artifacts removed). Compared to conventional ICA, ERASE removed on average 26% more EMG artifacts from EEG. These results indicate that using additional real or simulated EMG sources can increase the effectiveness of ICA in removing EMG artifacts from EEG. Combined with automated artifact IC rejection, ERASE also minimizes potential user bias.


title: Prioritized Restreaming Algorithms for Balanced Graph Partitioning

abstract: Balanced graph partitioning is a critical step for many large-scale distributed computations with relational data. As graph datasets have grown in size and density, a range of highly-scalable balanced partitioning algorithms have appeared to meet varied demands across different domains. As the starting point for the present work, we observe that two recently introduced families of iterative partitioners---those based on restreaming and those based on balanced label propagation (including Facebook's Social Hash Partitioner)---can be viewed through a common modular framework of design decisions. With the help of this modular perspective, we find that a key combination of design decisions leads to a novel family of algorithms with notably better empirical performance than any existing highly-scalable algorithm on a broad range of real-world graphs. The resulting prioritized restreaming algorithms employ a constraint management strategy based on multiplicative weights, borrowed from the restreaming literature, while adopting notions of priority from balanced label propagation to optimize the ordering of the streaming process. Our experimental results consider a range of stream orders, where a dynamic ordering based on what we call ambivalence is broadly the most performative in terms of the cut quality of the resulting balanced partitions, with a static ordering based on degree being nearly as good.


title: Learning to learn generative programs with Memoised Wake-Sleep

abstract: We study a class of neuro-symbolic generative models in which neural networks are used both for inference and as priors over symbolic, data-generating programs. As generative models, these programs capture compositional structures in a naturally explainable form. To tackle the challenge of performing program induction as an 'inner-loop' to learning, we propose the Memoised Wake-Sleep (MWS) algorithm, which extends Wake Sleep by explicitly storing and reusing the best programs discovered by the inference network throughout training. We use MWS to learn accurate, explainable models in three challenging domains: stroke-based character modelling, cellular automata, and few-shot learning in a novel dataset of real-world string concepts.


title: The Sample Complexity of Best-$k$ Items Selection from Pairwise Comparisons

abstract: This paper studies the sample complexity (aka number of comparisons) bounds for the active best-$k$ items selection from pairwise comparisons. From a given set of items, the learner can make pairwise comparisons on every pair of items, and each comparison returns an independent noisy result about the preferred item. At any time, the learner can adaptively choose a pair of items to compare according to past observations (i.e., active learning). The learner's goal is to find the (approximately) best-$k$ items with a given confidence, while trying to use as few comparisons as possible. In this paper, we study two problems: (i) finding the probably approximately correct (PAC) best-$k$ items and (ii) finding the exact best-$k$ items, both under strong stochastic transitivity and stochastic triangle inequality. For PAC best-$k$ items selection, we first show a lower bound and then propose an algorithm whose sample complexity upper bound matches the lower bound up to a constant factor. For the exact best-$k$ items selection, we first prove a worst-instance lower bound. We then propose two algorithms based on our PAC best items selection algorithms: one works for $k=1$ and is sample complexity optimal up to a loglog factor, and the other works for all values of $k$ and is sample complexity optimal up to a log factor.


title: A group-theoretical classification of three-tone and four-tone harmonic chords

abstract: We classify three-tone and four-tone chords based on subgroups of the symmetric group acting on chords contained within a twelve-tone scale. The actions are inversion, major-minor duality, and augmented-diminished duality. These actions correspond to elements of symmetric groups, and also correspond directly to intuitive concepts in the harmony theory of music. We produce a graph of how these actions relate different seventh chords that suggests a concept of distance in the theory of harmony.


title: Effective equidistribution of horospherical flows in infinite volume rank one homogeneous spaces

abstract: We prove effective equidistribution of horospherical flows in $\operatorname{SO}(n,1)^\circ / \Gamma$ when $\Gamma$ is geometrically finite and the frame flow is exponentially mixing for the Bowen-Margulis-Sullivan measure. We also discuss settings in which such an exponential mixing result is known to hold. As part of the proof, we show that the Patterson-Sullivan measure satisfies some friendly-like properties when $\Gamma$ is geometrically finite.


title: Electromyogram (EMG) Removal by Adding Sources of EMG (ERASE) -- A novel ICA-based algorithm for removing myoelectric artifacts from EEG -- Part 2

abstract: Extraction of the movement-related high-gamma (80 - 160 Hz) in electroencephalogram (EEG) from traumatic brain injury (TBI) patients who have had hemicraniectomies, remains challenging due to a confounding bandwidth overlap with surface electromyogram (EMG) artifacts related to facial and head movements. In part 1, we described an augmented independent component analysis (ICA) approach for removal of EMG artifacts from EEG, and referred to as EMG Reduction by Adding Sources of EMG (ERASE). Here, we tested ERASE on EEG recorded from six TBI patients with hemicraniectomies while they performed a thumb flexion task. ERASE removed a mean of 52 +/- 12% (mean +/- S.E.M) (maximum 73%) of EMG artifacts. In contrast, conventional ICA removed a mean of 27 +/- 19\% (mean +/- S.E.M) of EMG artifacts from EEG. In particular, high-gamma synchronization was significantly improved in the contralateral hand motor cortex area within the hemicraniectomy site after ERASE was applied. We computed fractal dimension (FD) of EEG high-gamma on each channel. We found relative FD of high-gamma over hemicraniectomy after applying ERASE were strongly correlated to the amplitude of finger flexion force. Results showed that significant correlation coefficients across the electrodes related to thumb flexion averaged 0.76, while the coefficients across the homologous electrodes in non-hemicraniectomy areas were nearly 0. Across all subjects, an average of 83% of electrodes significantly correlated with force was located in the hemicraniectomy areas after applying ERASE. After conventional ICA, only 19% of electrodes with significant correlations were located in the hemicraniectomy. These results indicated that the new approach isolated electrophysiological features during finger motor activation while selectively removing confounding EMG artifacts.


title: Predicting Afrobeats Hit Songs Using Spotify Data

abstract: This study approached the Hit Song Science problem with the aim of predicting which songs in the Afrobeats genre will become popular among Spotify listeners. A dataset of 2063 songs was generated through the Spotify Web API, with the provided audio features. Random Forest and Gradient Boosting algorithms proved to be successful with approximately F1 scores of 86%.


title: Long range vortex configurations in generalized models with the Maxwell or Chern-Simons dynamics

abstract: In this work we deal with vortices in Maxwell-Higgs or Chern-Simons-Higgs models that engender long range tails. We find first order differential equations that support minimum energy solutions which solve the equations of motion. In the Maxwell scenario, we work with generalised magnetic permeabilities that lead to vortices described by solutions, magnetic field and energy density with power-law tails that extend farther than the standard exponential ones. We also find a manner to obtain a Chern-Simons model with the same scalar and magnetic field profiles of the Maxwell case. By doing so, we also find vortices with the aforementioned long range feature, which is also present in the electric field in the Chern-Simons model. The present results may motivate investigations on nonrelativistic models, in particular in the case involving Rydberg atoms, which are known to present long range interactions and relatively long lifetimes.


title: Loosely trapped surface and dynamically transversely trapping surface in Einstein-Maxwell system

abstract: We study the properties of the loosely trapped surface (LTS) and the dynamically transversely trapping surface (DTTS) in Einstein-Maxwell systems. These concepts of surfaces were proposed by the four of the present authors in order to characterize strong gravity regions. We prove the Penrose-like inequalities for the area of LTSs/DTTSs. Interestingly, although the naively expected upper bound for the area is that of the photon sphere of a Reissner-Nordstroem black hole with the same mass and charge, the obtained inequalities include corrections represented by the energy density or pressure/tension of electromagnetic fields. Due to this correction, the Penrose-like inequality for the area of LTSs is tighter than the naively expected one. We also evaluate the correction term numerically in the Majumdar-Papapetrou two-black-hole spacetimes.


title: Research on Annotation Rules and Recognition Algorithm Based on Phrase Window

abstract: At present, most Natural Language Processing technology is based on the results of Word Segmentation for Dependency Parsing, which mainly uses an end-to-end method based on supervised learning. There are two main problems with this method: firstly, the la-beling rules are complex and the data is too difficult to label, the workload of which is large; secondly, the algorithm cannot recognize the multi-granularity and diversity of language components. In order to solve these two problems, we propose labeling rules based on phrase windows, and designed corresponding phrase recognition algorithms. The labeling rule uses phrases as the minimum unit, di-vides sentences into 7 types of nestable phrase types, and marks the grammatical dependencies between phrases. The corresponding algorithm, drawing on the idea of identifying the target area in the image field, can find the start and end positions of various phrases in the sentence, and realize the synchronous recognition of nested phrases and grammatical dependencies. The results of the experiment shows that the labeling rule is convenient and easy to use, and there is no ambiguity; the algorithm is more grammatically multi-granular and diverse than the end-to-end algorithm. Experiments on the CPWD dataset improve the accuracy of the end-to-end method by about 1 point. The corresponding method was applied to the CCL2018 competition, and the first place in the Chinese Metaphor Sentiment Analysis Task.


title: Dual Mixup Regularized Learning for Adversarial Domain Adaptation

abstract: Recent advances on unsupervised domain adaptation (UDA) rely on adversarial learning to disentangle the explanatory and transferable features for domain adaptation. However, there are two issues with the existing methods. First, the discriminability of the latent space cannot be fully guaranteed without considering the class-aware information in the target domain. Second, samples from the source and target domains alone are not sufficient for domain-invariant feature extracting in the latent space. In order to alleviate the above issues, we propose a dual mixup regularized learning (DMRL) method for UDA, which not only guides the classifier in enhancing consistent predictions in-between samples, but also enriches the intrinsic structures of the latent space. The DMRL jointly conducts category and domain mixup regularizations on pixel level to improve the effectiveness of models. A series of empirical studies on four domain adaptation benchmarks demonstrate that our approach can achieve the state-of-the-art.


title: An Optimal Energy-Saving Home Energy Management Supporting User Comfort and Electricity Selling with Different Prices

abstract: In this study, we investigate the operation of an optimal home energy management system (HEMS) with integrated renewable energy system (RES) and energy storage system (ESS) supporting electricity selling functions. A multi-objective mixed integer nonlinear programming model, including RES, ESS, home appliances and the main grid, is proposed to optimize different and conflicting objectives which are energy cost, user comfort and PAR. The effect of different selling prices on the objectives is also considered in detail. We further develop a formula for the lower bound of energy cost to help residents or engineers quickly choose best parameters of RES and ESS for their homes during the installation process. The performance of our system is verified through extensive simulations under three different scenarios of normal, economic, and smart with different selling prices using real data, and simulation results are compared in terms of daily energy cost, PAR, user's convenience and consecutive waiting time to use appliances. Numerical results clearly show that the economic scenario achieves 51.6% reduction of daily energy cost compared to the normal scenario while sacrificing the user's convenience, PAR, and consecutive waiting time by 49%, 132%, and 1 hour, respectively. On the other hand, the smart scenario shows only slight degradation of user's convenience and PAR by 2% and 18%, respectively while achieving 46.4% reduction of daily energy cost and the same level of consecutive waiting time. Furthermore, our simulation results show that a decrease of selling prices has tiny impacts on PAR and user comfort even though the daily energy cost increases.


title: On the Efficient Evaluation of the Exchange Correlation Potential on Graphics Processing Unit Clusters

abstract: The predominance of Kohn-Sham density functional theory (KS-DFT) for the theoretical treatment of large experimentally relevant systems in molecular chemistry and materials science relies primarily on the existence of efficient software implementations which are capable of leveraging the latest advances in modern high performance computing (HPC). With recent trends in HPC leading towards in increasing reliance on heterogeneous accelerator based architectures such as graphics processing units (GPU), existing code bases must embrace these architectural advances to maintain the high-levels of performance which have come to be expected for these methods. In this work, we purpose a three-level parallelism scheme for the distributed numerical integration of the exchange-correlation (XC) potential in the Gaussian basis set discretization of the Kohn-Sham equations on large computing clusters consisting of multiple GPUs per compute node. In addition, we purpose and demonstrate the efficacy of the use of batched kernels, including batched level-3 BLAS operations, in achieving high-levels of performance on the GPU. We demonstrate the performance and scalability of the implementation of the purposed method in the NWChemEx software package by comparing to the existing scalable CPU XC integration in NWChem.


title: On the equidistribution of unstable curves for pseudo-Anosov diffeomorphisms of compact surfaces

abstract: Weprovethattheasymptoticsofergodicintegralsalonganinvariant foliation of a toral Anosov diffeomorphism, or of a pseudo-Anosov diffeomorphism on a compact orientable surface of higher genus, are determined (up to a logarithmic error) by the action of the diffeomorphism on the cohomology of the surface. As a consequence of our argument and of the results of Giulietti and Liverani [GL] on horospherical averages, toral Anosov diffeomorphisms have no Ruelle resonances in the open interval $(1,e^{h_{top}} )$.


title: Stability and breakup of liquid threads and annular layers in a corrugated tube

abstract: We study the stability and breakup behavior of an axisymmetric liquid thread which is surrounded by another viscous fluid layer through a long wave approximation. The two fluids are immiscible and confined in a concentrically placed cylindrical tube. The effect of the tube wall corrugation is taken into account in the model, which allows the access of the interaction between the wall shape and the thread interfacial dynamics. The linearized system is studied by the Floquet theory due to the presence of non-constant coefficients in the equation and the spectrum is computed numerically via the Fourier-Floquet-Hill method. The resulting features agree qualitatively with those obtained based on a lubrication model in the thin annulus limit, where the short-wave disturbances that would be stabilizing owing to the capillarity in the absence of wall corrugation, can excite some unstable long waves. Those results from the linear theory are also confirmed numerically by the direct numerical simulation on the evolution equations. Meanwhile, a transition on the dominant modes from the one for a tube without corrugation to the one with wall shape included is found. Finally the drop formation is shown in the nonlinear regime when the pinch off of the core thread is obtained. The annular film drainage regime is also approached slowly along with pinching when the tube wall is close to the thread interface. In addition, our results demonstrate the possibility to suppress pinching, depending on the averaged annular layer thickness and the variation in tube radius.


title: Determination of the phase coherence length of PdCoO$_2$ nanostructures

abstract: The two-dimensional layered compound PdCoO$_2$ is one of the best oxide conductors, providing an intriguing research arena opened by the long mean free path and the very high mobility of ~51000 cm2/Vs. These properties turn PdCoO$_2$ into a candidate material for nanoscale quantum devices. By exploring universal conductance fluctuations originating at nanoscale PdCoO$_2$ Hall-bar devices, we determined the phase coherence length of electron transport in c-axis oriented PdCoO$_2$ thin films to equal ~100 nm. The weak temperature dependence of the measured phase coherence length suggests that defect scattering at twin boundaries in the PdCoO$_2$ thin film governs phase breaking. These results suggest that phase coherent devices can be achieved by realizing the devices smaller than the size of twin domains, via refined microfabrication and suppression of twin boundaries.


title: Unidirectional spin transport of a spin-orbit-coupled atomic matter wave using a moving Dirac $\delta$-potential well

abstract: We study the transport of a spin-orbit-coupled atomic matter wave using a moving Dirac $\delta$-potential well. In a spin-orbit-coupled system, bound states can be formed in both ground and excited energy levels with a Dirac $\delta$ potential. Because Galilean invariance is broken in a spin-orbit-coupled system, moving of the potential will induce a velocity-dependent effective detuning. This induced detuning breaks the spin symmetry and makes the ground-state transporting channel be spin-$\uparrow$ ($\downarrow$) favored while makes the excited-state transporting channel be spin-$\downarrow$ ($\uparrow$) favored for a positive-direction (negative-direction) transporting. When the $\delta$-potential well moves at a small velocity, both the ground-state and the excited-state channels contribute to the transportation, and thus both the spin components can be efficiently transported. However, when the moving velocity of the $\delta$-potential well exceeds a critical value, the induced detuning is large enough to eliminate the excited bound state, and makes the ground bound state the only transporting channel, in which only the spin-$\uparrow$ ($\downarrow$) component can be efficiently transported in a positive (negative) direction. This work demonstrates a prototype of unidirectional spin transport.


title: All-electron Gaussian-based $G_0W_0$ for Valence and Core Excitation Energies of Periodic Systems

abstract: We describe an all-electron $G_0W_0$ implementation for periodic systems with $k$-point sampling implemented in a crystalline Gaussian basis. Our full-frequency $G_0W_0$ method relies on efficient Gaussian density fitting integrals and includes both analytic continuation and contour deformation schemes. Due to the compactness of Gaussian bases, no virtual state truncation is required as is seen in many plane-wave formulations. Finite size corrections are included by taking the $q \to 0$ limit of the Coulomb divergence. Using our implementation, we study quasiparticle energies and band structures across a range of systems including molecules, semiconductors, rare gas solids, and metals. We find that the $G_0W_0$ band gaps of traditional semiconductors converge rapidly with respect to the basis size, even for the conventionally challenging case of ZnO. Using correlation-consistent bases of polarized triple-zeta quality, we find the mean absolute relative error of the extrapolated $G_0W_0$@PBE band gaps to be only 5.2% when compared to experimental values. For core excitation binding energies (CEBEs), we find that $G_0W_0$ predictions improve significantly over those from DFT if the $G_0W_0$ calculations are started from hybrid functionals with a high percentage of exact exchange.


title: A Stochastic-Robust Approach for Resilient Microgrid Investment Planning Under Static and Transient Islanding Security Constraints

abstract: When planning the investment in Microgrids (MGs), usually static security constraints are included to ensure their resilience and ability to operate in islanded mode. However, unscheduled islanding events may trigger cascading disconnections of Distributed Energy Resources (DERs) inside the MG due to the transient response, leading to a partial or full loss of load. In this paper, a min-max-min, hybrid, stochastic-robust investment planning model is proposed to obtain a resilient MG considering both High-Impact-Low-Frequency (HILF) and Low-Impact-High-Frequency (LIHF) uncertainties. The HILF uncertainty pertains to the unscheduled islanding of the MG after a disastrous event, and the LIHF uncertainty relates to correlated loads and DER generation, characterized by a set of scenarios. The MG resilience under both types of uncertainty is ensured by incorporating static and transient islanding constraints into the proposed investment model. The inclusion of transient response constraints leads to a min-max-min problem with a non-linear dynamic frequency response model that cannot be solved directly by available optimization tools. Thus, in this paper, a three-stage solution approach is proposed to find the optimal investment plan. The performance of the proposed algorithm is tested on the CIGRE 18-node distribution network.


title: Using the Mark Weighted Correlation Functions to Improve the Constraints on Cosmological Parameters

abstract: We used the mark weighted correlation functions (MCFs), $W(s)$, to study the large scale structure of the Universe. We studied five types of MCFs with the weighting scheme $\rho^\alpha$, where $\rho$ is the local density, and $\alpha$ is taken as $-1,\ -0.5,\ 0,\ 0.5$, and 1. We found that different MCFs have very different amplitudes and scale-dependence. Some of the MCFs exhibit distinctive peaks and valleys that do not exist in the standard correlation functions. Their locations are robust against the redshifts and the background geometry, however it is unlikely that they can be used as ``standard rulers'' to probe the cosmic expansion history. Nonetheless we find that these features may be used to probe parameters related with the structure formation history, such as the values of $\sigma_8$ and the galaxy bias. Finally, after conducting a comprehensive analysis using the full shapes of the $W(s)$s and $W_{\Delta s}(\mu)$s, we found that, combining different types of MCFs can significantly improve the cosmological parameter constraints. Compared with using only the standard correlation function, the combinations of MCFs with $\alpha=0,\ 0.5,\ 1$ and $\alpha=0,\ -1,\ -0.5,\ 0.5,\ 1$ can improve the constraints on $\Omega_m$ and $w$ by $\approx30\%$ and $50\%$, respectively. We find highly significant evidence that MCFs can improve cosmological parameter constraints.


title: Curriculum learning for multilevel budgeted combinatorial problems

abstract: Learning heuristics for combinatorial optimization problems through graph neural networks have recently shown promising results on some classic NP-hard problems. These are single-level optimization problems with only one player. Multilevel combinatorial optimization problems are their generalization, encompassing situations with multiple players taking decisions sequentially. By framing them in a multi-agent reinforcement learning setting, we devise a value-based method to learn to solve multilevel budgeted combinatorial problems involving two players in a zero-sum game over a graph. Our framework is based on a simple curriculum: if an agent knows how to estimate the value of instances with budgets up to $B$, then solving instances with budget $B+1$ can be done in polynomial time regardless of the direction of the optimization by checking the value of every possible afterstate. Thus, in a bottom-up approach, we generate datasets of heuristically solved instances with increasingly larger budgets to train our agent. We report results close to optimality on graphs up to $100$ nodes and a $185 \times$ speedup on average compared to the quickest exact solver known for the Multilevel Critical Node problem, a max-min-max trilevel problem that has been shown to be at least $\Sigma_2^p$-hard.


title: The gem5 Simulator: Version 20.0+

abstract: The open-source and community-supported gem5 simulator is one of the most popular tools for computer architecture research. This simulation infrastructure allows researchers to model modern computer hardware at the cycle level, and it has enough fidelity to boot unmodified Linux-based operating systems and run full applications for multiple architectures including x86, Arm, and RISC-V. The gem5 simulator has been under active development over the last nine years since the original gem5 release. In this time, there have been over 7500 commits to the codebase from over 250 unique contributors which have improved the simulator by adding new features, fixing bugs, and increasing the code quality. In this paper, we give and overview of gem5's usage and features, describe the current state of the gem5 simulator, and enumerate the major changes since the initial release of gem5. We also discuss how the gem5 simulator has transitioned to a formal governance model to enable continued improvement and community support for the next 20 years of computer architecture research.


title: High-speed Millimeter-wave 5G/6G Image Transmission via Artificial Intelligence

abstract: Artificial Intelligence (AI) has been used to jointly optimize a mmWave Compressed Sensing (CS) for high-speed 5G/6G image transmission. Specifically, we have developed a Dictionary Learning Compressed Sensing neural Network (DL-CSNet) to realize three key functionalities: 1) to learn the dictionary basis of the images for transmission; 2) to optimize the Hadamard measurement matrix; and 3) to reconstruct the lossless images with the learned dictionary basis. A 94-GHz prototype has been built and up to one order of image transmission speed increase has been realized for letters ``A" to ``Z".


title: Discretization-Aware Architecture Search

abstract: The search cost of neural architecture search (NAS) has been largely reduced by weight-sharing methods. These methods optimize a super-network with all possible edges and operations, and determine the optimal sub-network by discretization, \textit{i.e.}, pruning off weak candidates. The discretization process, performed on either operations or edges, incurs significant inaccuracy and thus the quality of the final architecture is not guaranteed. This paper presents discretization-aware architecture search (DA\textsuperscript{2}S), with the core idea being adding a loss term to push the super-network towards the configuration of desired topology, so that the accuracy loss brought by discretization is largely alleviated. Experiments on standard image classification benchmarks demonstrate the superiority of our approach, in particular, under imbalanced target network configurations that were not studied before.


title: Decentralized policy learning with partial observation and mechanical constraints for multiperson modeling

abstract: Extracting the rules of real-world multi-agent behaviors is a current challenge in various scientific and engineering fields. Biological agents independently have limited observation and mechanical constraints; however, most of the conventional data-driven models ignore such assumptions, resulting in lack of biological plausibility and model interpretability for behavioral analyses. Here we propose sequential generative models with partial observation and mechanical constraints in a decentralized manner, which can model agents' cognition and body dynamics, and predict biologically plausible behaviors. We formulate this as a decentralized multi-agent imitation-learning problem, leveraging binary partial observation and decentralized policy models based on hierarchical variational recurrent neural networks with physical and biomechanical penalties. Using real-world basketball and soccer datasets, we show the effectiveness of our method in terms of the constraint violations, long-term trajectory prediction, and partial observation. Our approach can be used as a multi-agent simulator to generate realistic trajectories using real-world data.


title: Speed-of-sound imaging by differential phase contrast with angular compounding

abstract: We describe a technique to reveal speed-of-sound (SoS) variations within an echogenic sample. The technique uses the same receive data as standard pulse-echo imaging based on plane-wave compounding, and can be operated in parallel. Point-like scatterers randomly distributed throughout the sample serve as local probes of the downstream transmit-beam phase shifts caused by aberrating structures within the sample. Phase shifts are monitored in a differential manner, providing signatures of transverse gradients of the local sample SoS. The contrast of the signatures is augmented by a method of angular compounding, which provides ``focus" control of the image sharpness, which, in turn, enables a visual localization of aberrating inclusions within the sample on the fly. The localization can be performed in 2D when operated with standard B-mode imaging, or in 3D when operated with C-mode imaging. Finally, we present a wave-acoustic forward model that provides insight into the principle of differential phase contrast (DPC) imaging, and roughly recapitulates experimental results obtained with an elastography phantom. In particular, we demonstrate that our technique easily reveals relative SoS variations as small as 0.5\% in real time. Such imaging may ultimately be useful for clinical diagnosis of pathologies in soft tissue.


title: Size-regulated symmetry breaking in reaction-diffusion models of developmental transitions

abstract: The development of multicellular organisms proceeds through a series of morphogenetic and cell-state transitions, transforming homogeneous zygotes into complex adults by a process of self-organization. Many of these transitions are achieved by spontaneous symmetry breaking mechanisms, allowing cells and tissues to acquire pattern and polarity by virtue of local interactions without an upstream supply of information. The combined work of theory and experiment has elucidated how these systems break symmetry during developmental transitions. Given such transitions are multiple and their temporal ordering is crucial, an equally important question is how these developmental transitions are coordinated in time. Using a minimal mass-conserved substrate-depletion model for symmetry breaking as our case study, we elucidate mechanisms by which cells and tissues can couple reaction-diffusion driven symmetry breaking to the timing of developmental transitions, arguing that the dependence of patterning mode on system size may be a generic principle by which developing organisms measure time. By analyzing different regimes of our model, simulated on growing domains, we elaborate three distinct behaviours, allowing for clock-, timer-, or switch-like dynamics. By relating these behaviours to experimentally documented case studies of developmental timing, we provide a minimal conceptual framework to interrogate how developing organisms coordinate developmental transitions.


title: The LoCA Regret: A Consistent Metric to Evaluate Model-Based Behavior in Reinforcement Learning

abstract: Deep model-based Reinforcement Learning (RL) has the potential to substantially improve the sample-efficiency of deep RL. While various challenges have long held it back, a number of papers have recently come out reporting success with deep model-based methods. This is a great development, but the lack of a consistent metric to evaluate such methods makes it difficult to compare various approaches. For example, the common single-task sample-efficiency metric conflates improvements due to model-based learning with various other aspects, such as representation learning, making it difficult to assess true progress on model-based RL. To address this, we introduce an experimental setup to evaluate model-based behavior of RL methods, inspired by work from neuroscience on detecting model-based behavior in humans and animals. Our metric based on this setup, the Local Change Adaptation (LoCA) regret, measures how quickly an RL method adapts to a local change in the environment. Our metric can identify model-based behavior, even if the method uses a poor representation and provides insight in how close a method's behavior is from optimal model-based behavior. We use our setup to evaluate the model-based behavior of MuZero on a variation of the classic Mountain Car task.


title: Q-data and representation theory of untwisted quantum affine algebras

abstract: For a complex finite-dimensional simple Lie algebra $\mathfrak{g}$, we introduce the notion of Q-datum, which generalizes the notion of a Dynkin quiver with a height function from the viewpoint of Weyl group combinatorics. Using this notion, we develop a unified theory describing the twisted Auslander-Reiten quivers and the twisted adapted classes introduced in [O.-Suh, J. Algebra, 2019] with an appropriate notion of the generalized Coxeter elements. As a consequence, we obtain a combinatorial formula expressing the inverse of the quantum Cartan matrix of $\mathfrak{g}$, which generalizes the result of [Hernandez-Leclerc, J. Reine Angew. Math., 2015] in the simply-laced case. We also find several applications of our combinatorial theory of Q-data to the finite-dimensional representation theory of the untwisted quantum affine algebra of $\mathfrak{g}$. In particular, in terms of Q-data and the inverse of the quantum Cartan matrix, (i) we give an alternative description of the block decomposition results due to [Chari-Moura, Int. Math. Res. Not., 2005] and [Kashiwara-Kim-O.-Park, arXiv:2003.03265], (ii) we present a unified (partially conjectural) formula of the denominators of the normalized R-matrices between all the Kirillov-Reshetikhin modules, and (iii) we compute the invariants $\Lambda(V,W)$ and $\Lambda^\infty(V, W)$ introduced in [Kashiwara-Kim-O.-Park, Compos. Math., 2020] for each pair of simple modules $V$ and $W$.


title: Dynamic coherent perfect absorption in nonlinear metasurfaces

abstract: We propose a tunable coherent perfect absorber based on ultrathin nonlinear metasurfaces. The nonlinear metasurface is made of plasmonic nanoantennas coupled to an epsilon-near-zero material with a large optical nonlinearity. The coherent perfect absorption is achieved by controlling the relative phases of the input beams. We show that the optical response of the nonlinear metasurface can be tuned from a complete to a partial absorption by changing the intensity of the pump beam. The proposed nonlinear metasurface can be used to design optically tunable thermal emitters, modulators, and sensors.


title: Decay properties of beauty and charm mesons within Isgur-Wise function formalism

abstract: We investigate the decay properties of some beauty and charm mesons with a phenomenological potential model. First, we consider the nonrelativistic Hamiltonian of the mesonic system with Coulomb plus exponential terms and study the wave function and the energy of the system using the variational approach. Thereby, we compute the masses, the decay constants, the leptonic branching fractions of heavy-light mesons and the mixing mass parameter $\Delta {m_{{B_q}}}$. We study the radiative leptonic decay widths of ${D_s} \to \gamma \ell \bar \nu $, ${D^ - } \to \gamma \ell \bar \nu $ and the semileptonic decay widths of ${\bar B_{(s)}} \to {D_{(s)}}\ell \bar \nu $, ${\bar B_{(s)}} \to D_{(s)}^*\ell \bar \nu $. Using Isgur-Wise functions, we calculate the branching ratios of $B \to {D^{(*)}}\pi $ and two-body nonleptonic decay of $D \to K\pi $. Our results are consistent with other theoretical models and the experimental results.