title: A Performance Guarantee for Spectral Clustering

abstract: The two-step spectral clustering method, which consists of the Laplacian eigenmap and a rounding step, is a widely used method for graph partitioning. It can be seen as a natural relaxation to the NP-hard minimum ratio cut problem. In this paper we study the central question: when is spectral clustering able to find the global solution to the minimum ratio cut problem? First we provide a condition that naturally depends on the intra- and inter-cluster connectivities of a given partition under which we may certify that this partition is the solution to the minimum ratio cut problem. Then we develop a deterministic two-to-infinity norm perturbation bound for the the invariant subspace of the graph Laplacian that corresponds to the $k$ smallest eigenvalues. Finally by combining these two results we give a condition under which spectral clustering is guaranteed to output the global solution to the minimum ratio cut problem, which serves as a performance guarantee for spectral clustering.


title: On invariant (co)homology of a group

abstract: There are different notions of homology and cohomology that can be defined for a group with an action of another group by group automorphisms. In this paper we address three natural questions that arise in this context. Namely, the relation of these notions with the usual (co)homology of a semidirect product, the interpretation of the first homology group as some kind of abelianization and the classification of (invariant) group extensions.


title: Diversity Improves Speed and Accuracy in Social Networks

abstract: How does temporally structured private and social information shape collective decisions? To address this question we consider a network of rational agents who independently accumulate private evidence that triggers a decision upon reaching a threshold. When seen by the whole network, the first agent's choice initiates a wave of new decisions; later decisions have less impact. In heterogeneous networks, first decisions are made quickly by impulsive individuals who need little evidence to make a choice, but, even when wrong, can reveal the correct options to nearly everyone else. We conclude that groups comprised of diverse individuals can make more efficient decisions than homogenous ones.


title: Computing Dense and Sparse Subgraphs of Weakly Closed Graphs

abstract: A graph $G$ is weakly $\gamma$-closed if every induced subgraph of $G$ contains one vertex $v$ such that for each non-neighbor $u$ of $v$ it holds that $|N(u)\cap N(v)|<\gamma$. The weak closure $\gamma(G)$ of a graph, recently introduced by Fox et al. [SIAM J. Comp. 2020], is the smallest number such that $G$ is weakly $\gamma$-closed. This graph parameter is never larger than the degeneracy (plus one) and can be significantly smaller. Extending the work of Fox et al. [SIAM J. Comp. 2020] on clique enumeration, we show that several problems related to finding dense subgraphs, such as the enumeration of bicliques and $s$-plexes, are fixed-parameter tractable with respect to $\gamma(G)$. Moreover, we show that the problem of determining whether a weakly $\gamma$-closed graph $G$ has a subgraph on at least $k$ vertices that belongs to a graph class $\mathcal{G}$ which is closed under taking subgraphs admits a kernel with at most $\gamma k^2$ vertices. Finally, we provide fixed-parameter algorithms for Independent Dominating Set and Dominating Clique when parameterized by $\gamma+k$ where $k$ is the solution size.


title: Joint Access Point Selection and Interference Cancellation for Cell-Free Massive MIMO

abstract: Cell-Free Massive MIMO is a highly promising approach to enhance network capacity by moving a large number of distributed access points (AP) closer to mobile users while utilizing simple matched filtering and conjugate beamforming. Recent work using minimum mean-squared-error (MMSE) receiver that suppress multi-user interference (MUI) shows significant capacity increase, but at the cost of high computational complexity and residual MUI enhancement. We propose a significantly lower complexity adaptive approach where central processing unit (CPU) removes MUI without amplifying the residual interference. It does so dynamically by using available knowledge of channel estimates to perform joint process of combining selected strongest AP signals for each user and subtracting the sum of interference estimates from other users at the same time. We provide signalto-interference plus noise-ratio (SINR) and complexity analyses backed by numerical results to show the superiority of this approach compared with the state-of-the-art techniques.


title: Measuring segregation characteristics of industrially relevant granular mixtures: Part I -- A continuum model approach

abstract: We present a method to estimate the segregation parameter, $S,$ a key input in a continuum transport model of particulate flows. $S$ is determined by minimizing the difference between measured and model-predicted concentration profiles. To validate the approach, we conduct discrete element method simulations of size-bidisperse mixtures in quasi-2D bounded heap flow; the resulting data show that $S$ calculated from concentration profiles is consistent with the directly measured value. The method's accuracy depends critically on the velocity profile during filling, but only weakly on the diffusion coefficient. When the velocity profile is nominally spanwise invariant, the error between estimated and measured $S$ is $10\%$. This method is intended for practical application (described in Part II), so we restrict characterization of the velocity profile to that which can be readily determined experimentally, and explore the sensitivity of concentration profiles to variation of the gap between the sidewalls of the heap.


title: High Resolution Spectral Metrology Leveraging Topologically Enhanced Optical Activity in Fibers

abstract: Optical rotation, a form of optical activity, is a phenomenon employed in various metrological applications and industries including chemical, food, and pharmaceutical. In naturally-occurring, as well as structured media, the integrated effect is, however, typically small. Here, we demonstrate that, by exploiting the inherent and stable spin-orbit interaction of orbital angular momentum fiber modes, giant, scalable optical activity can be obtained, and that we can use this effect to realize a new type of wavemeter by exploiting its optical rotary dispersion. The device we construct provides for an instantaneous wavelength-measurement technique with high resolving power $R =3.4 \times 10^6$ (i.e. resolution $< 0.3 pm$ at $1-{\mu}m$ wavelengths) and can also detect spectral bandwidths of known lineshapes with high sensitivity.


title: Vector Balancing in Lebesgue Spaces

abstract: A tantalizing conjecture in discrete mathematics is the one of Koml\'os, suggesting that for any vectors $\mathbf{a}_1,\ldots,\mathbf{a}_n \in B_2^m$ there exist signs $x_1, \dots, x_n \in \{ -1,1\}$ so that $\|\sum_{i=1}^n x_i\mathbf{a}_i\|_\infty \le O(1)$. It is a natural extension to ask what $\ell_q$-norm bound to expect for $\mathbf{a}_1,\ldots,\mathbf{a}_n \in B_p^m$. We prove that, for $2 \le p \le q \le \infty$, such vectors admit fractional colorings $x_1, \dots, x_n \in [-1,1]$ with a linear number of $\pm 1$ coordinates so that $\|\sum_{i=1}^n x_i\mathbf{a}_i\|_q \leq O(\sqrt{\min(p,\log(2m/n))}) \cdot n^{1/2-1/p+ 1/q}$, and that one can obtain a full coloring at the expense of another factor of $\frac{1}{1/2 - 1/p + 1/q}$. In particular, for $p \in (2,3]$ we can indeed find signs $\mathbf{x} \in \{ -1,1\}^n$ with $\|\sum_{i=1}^n x_i\mathbf{a}_i\|_\infty \le O(n^{1/2-1/p} \cdot \frac{1}{p-2})$. Our result generalizes Spencer's theorem, for which $p = q = \infty$, and is tight for $m = n$. Additionally, we prove that for any fixed constant $\delta>0$, in a centrally symmetric body $K \subseteq \mathbb{R}^n$ with measure at least $e^{-\delta n}$ one can find such a fractional coloring in polynomial time. Previously this was known only for a small enough constant -- indeed in this regime classical nonconstructive arguments do not apply and partial colorings of the form $\mathbf{x} \in \{ -1,0,1\}^n$ do not necessarily exist.


title: From unsteady to quasi-steady dynamics in the streamwise-oscillating cylinder wake

abstract: The flow around a cylinder oscillating in the streamwise direction with a frequency, $f_f$, much lower than the shedding frequency, $f_s$, has been relatively less studied than the case when these frequencies have the same order of magnitude, or the transverse oscillation configuration. In this study, Particle Image Velocimetry and Koopman Mode Decomposition are used to investigate the streamwise-oscillating cylinder wake for forcing frequencies $f_f/f_s \sim 0.04-0.2$ and mean Reynolds number, $Re_0 = 900$. The amplitude of oscillation is such that the instantaneous Reynolds number remains above the critical value for vortex shedding at all times. Characterization of the wake reveals a range of phenomena associated with the interaction of the two frequencies, including modulation of both the amplitude and frequency of the wake structure by the forcing. Koopman analysis reveals a frequency spreading of Koopman modes. A scaling parameter and associated transformation are developed to relate the unsteady, or forced, dynamics of a system to that of a quasi-steady, or unforced, system. For the streamwise-oscillating cylinder, it is shown that this transformation leads to a Koopman Mode Decomposition similar to that of the unforced system.


title: Adaptive Superresolution in Deconvolution of Sparse Peaks

abstract: The aim of this paper is to investigate superresolution in deconvolution driven by sparsity priors. The observed signal is a convolution of an original signal with a continuous kernel.With the prior knowledge that the original signal can be considered as a sparse combination of Dirac delta peaks, we seek to estimate the positions and amplitudes of these peaks by solving a finite dimensional convex problem on a computational grid. Because, the support of the original signal may or may not be on this grid, by studying the discrete deconvolution of sparse peaks using L1-norm sparsity prior, we confirm recent observations that canonically the discrete reconstructions will result in multiple peaks at grid points adjacent to the location of the true peak. Owning to the complexity of this problem, we analyse carefully the de-convolution of single peaks on a grid and gain a strong insight about the dependence of the reconstructed magnitudes on the exact peak location. This in turn allows us to infer further information on recovering the location of the exact peaks i.e. to perform super-resolution. We analyze in detail the possible cases that can appear and based on our theoretical findings, we propose an self-driven adaptive grid approach that allows to perform superresolution in one-dimensional and multi-dimensional spaces. With the view that the current study can provide a further step in the development of more robust algorithms for the detection of single molecules in fluorescence microscopy or identification of characteristic frequencies in spectral analysis, we demonstrate how the proposed approach can recover sparse signals using simulated clusters of point sources (peaks) of low-resolution in one and two-dimensional spaces.


title: Multilevel Digital Contact Tracing

abstract: Digital contact tracing plays a crucial role in alleviating an outbreak, and designing multilevel digital contact tracing for a country is an open problem due to the analysis of large volumes of temporal contact data. We develop a multilevel digital contact tracing framework that constructs dynamic contact graphs from the proximity contact data. Prominently, we introduce the edge label of the contact graph as a binary circular contact queue, which holds the temporal social interactions during the incubation period. After that, our algorithm prepares the direct and indirect (multilevel) contact list for a given set of infected persons from the contact graph. Finally, the algorithm constructs the infection pathways for the trace list. We implement the framework and validate the contact tracing process with synthetic and real-world data sets. In addition, analysis reveals that for COVID-19 close contact parameters, the framework takes reasonable space and time to create the infection pathways. Our framework can apply to any epidemic spreading by changing the algorithm's parameters.


title: On the Performance of Direct Shaping Codes

abstract: In this work, we study a recently proposed direct shaping code for flash memory. This rate-1 code is designed to reduce the wear for SLC (one bit per cell) flash by minimizing the average fraction of programmed cells when storing structured data. Then we describe an adaptation of this algorithm that provides data shaping for MLC (two bits per cell) flash memory. It makes use of a page-dependent cost model and is designed to be compatible with the standard procedure of row-by-row, page-based, wordline programming. We also give experimental results demonstrating the performance of MLC data shaping codes when applied to English and Chinese language text. We then study the potential error propagation properties of direct shaping codes when used in a noisy flash device. In particular, we model the error propagation as a biased random walk in a multidimensional space. We prove an upper bound on the error propagation probability and propose an algorithm that can numerically approach a lower bound. Finally, we study the asymptotic performance of direct shaping codes. We prove that the SLC direct shaping code is suboptimal in the sense that it can only achieve the minimum average cost for a rate-1 code under certain conditions on the source distribution.


title: Electric conductivity in finite-density SU(2) lattice gauge theory with dynamical fermions

abstract: We study the dependence of the electric conductivity on chemical potential in finite-density $SU(2)$ gauge theory with $N_f = 2$ flavours of rooted staggered sea quarks, in combination with Wilson-Dirac and Domain Wall valence quarks. The pion mass is reasonably small with $m_{\pi}/m_{\rho} \approx 0.4$. We concentrate in particular on the vicinity of the chiral crossover, where we find the low-frequency electric conductivity to be most sensitive to small changes in fermion density. Working in the low-density QCD-like regime with spontaneously broken chiral symmetry, we obtain an estimate of the first nontrivial coefficient $c(T)$ of the expansion of conductivity $\sigma(T,\mu) = \sigma(T,0) \left(1 + c(T) (\mu/T)^2 + O(\mu^4)\right)$ in powers of $\mu$, which has rather weak temperature dependence and takes its maximal value $c(T) \approx 0.10 \pm 0.07$ around the critical temperature. At larger densities and lower temperatures, the conductivity quickly grows towards the diquark condensation phase, and also becomes closer to the free quark result. As a by-product of our study we confirm the conclusions of previous studies with heavier pion that for $SU(2)$ gauge theory the ratio of crossover temperature to pion mass $T_c/m_{\pi} \approx 0.4$ at $\mu=0$ is significantly smaller than in real QCD.


title: Constraining early-time dust formation in core-collapse supernovae

abstract: There is currently a severe discrepancy between theoretical models of dust formation in core-collapse supernovae (CCSNe), which predict $\gtrsim 0.01$ M$_\odot$ of ejecta dust forming within $\sim 1000$ days, and observations at these epochs, which infer much lower masses. We demonstrate that, in the optically thin case, these low dust masses are robust despite significant observational and model uncertainties. For a sample of 11 well-observed CCSNe, no plausible model reaches carbon dust masses above $10^{-4}$ M$_\odot$, or silicate masses above $\sim 10^{-3}$ M$_\odot$. Optically thick models can accommodate larger dust masses, but the dust must be clumped and have a low ($<0.1$) covering fraction to avoid conflict with data at optical wavelengths. These values are insufficient to reproduce the observed infrared fluxes, and the required covering fraction varies not only between SNe but between epochs for the same object. The difficulty in reconciling large dust masses with early-time observations of CCSNe, combined with well-established detections of comparably large dust masses in supernova remnants, suggests that a mechanism for late-time dust formation is necessary.


title: Thin-shell wormholes in Rastall gravity

abstract: We constructed thin-shell wormholes using the well-known "cut and paste" technique for static black holes sourced by an anisotropic fluid within the context of the Rastall gravity. Using the generalized Lanczos equations we explored the energy conditions at the wormhole's throat along with the traversability condition. We determined the stability regions when the wormhole configurations are surrounded by different bulk sources and obtained that the stability regions are considerably modified due to the non-trivial dependence of the Rastall parameter in the effective potential energy.


title: Magnetic effects of QCD parameters from finite energy sum rules

abstract: One of the advantages of the finite energy sum rules is the fact that every operator in the operator product expansion series can be selected individually by the use of an appropriate kernel function which removes other operator poles. This characteristic is maintained by QCD systems in the presence of external homogeneous magnetic field, providing interesting information about the magnetic evolution of QCD and hadronic parameters. In this work finite energy sum rules are applied on QCD in the light quark sector, combining axial and pseudoscalar channels in the presence of an external homogeneous magnetic field, obtaining the magnetic evolution of the light quark masses, pion mass, the pion decay constant, the gluon condensate and the continuum hadronic threshold.


title: Learning Local Complex Features using Randomized Neural Networks for Texture Analysis

abstract: Texture is a visual attribute largely used in many problems of image analysis. Currently, many methods that use learning techniques have been proposed for texture discrimination, achieving improved performance over previous handcrafted methods. In this paper, we present a new approach that combines a learning technique and the Complex Network (CN) theory for texture analysis. This method takes advantage of the representation capacity of CN to model a texture image as a directed network and uses the topological information of vertices to train a randomized neural network. This neural network has a single hidden layer and uses a fast learning algorithm, which is able to learn local CN patterns for texture characterization. Thus, we use the weighs of the trained neural network to compose a feature vector. These feature vectors are evaluated in a classification experiment in four widely used image databases. Experimental results show a high classification performance of the proposed method when compared to other methods, indicating that our approach can be used in many image analysis problems.


title: D-dimensional self-gravitating lattice gas in general relativity

abstract: Using a lattice equation of state combined with the D-dimensional Tolman-Oppenheimer-Volkoff equation and the Friedmann equations, we investigate the possibility of the formation of compact objects as well as the time evolution of the scale factor and the density profile of a self-gravitating material cluster. The numerical results show that in a $2+1$ dimensional spacetime, the mass is independent of the central pressure. Hence, the formation of only compact objects with a finite constant mass similar to the white dwarf is possible. However, in a $3+1$ dimensional spacetime, self-gravity leads to the formation of compact objects with a large gap of mass and the corresponding phase diagram has the same structure as the one for Neutron Star. The results also show that beyond certain critical central pressure, the star is unstable against gravitational collapse, and it may end in a black hole. Analysis of spacetimes of higher dimensions shows that gravity has the stronger effect in $3+1$ dimensions. Numerical solutions of the Friedmann equations show that the effect of the curvature of spacetime increases with increasing temperature, but decreases with increasing dimensionality beyond $D=3$.


title: Nuclear Surface Acoustic Resonance with Spin-Rotation Coupling

abstract: We show that, under an appropriate out-of-plane static magnetic field, nuclear spins in a thin specimen on a surface acoustic wave (SAW) cavity can be resonantly excited and detected through spin-rotation coupling. Since such a SAW cavity can have the quality factor as high as $10^{4}$ and the mode volume as small as $10^{-2}$ mm$^{3}$ the signal-to-noise ratio in detecting the resonance is estimated to be quite high. We argue that detecting nuclear spin resonance of a single flake of an atomically-thin layer of two-dimensional semiconductor, which has so far been beyond hope with the conventional inductive method, can be a realistic target with the proposed scheme.


title: Transformations between deep neural networks

abstract: We propose to test, and when possible establish, an equivalence between two different artificial neural networks by attempting to construct a data-driven transformation between them, using manifold-learning techniques. In particular, we employ diffusion maps with a Mahalanobis-like metric. If the construction succeeds, the two networks can be thought of as belonging to the same equivalence class. We first discuss transformation functions between only the outputs of the two networks; we then also consider transformations that take into account outputs (activations) of a number of internal neurons from each network. In general, Whitney's theorem dictates the number of measurements from one of the networks required to reconstruct each and every feature of the second network. The construction of the transformation function relies on a consistent, intrinsic representation of the network input space. We illustrate our algorithm by matching neural network pairs trained to learn (a) observations of scalar functions; (b) observations of two-dimensional vector fields; and (c) representations of images of a moving three-dimensional object (a rotating horse). The construction of such equivalence classes across different network instantiations clearly relates to transfer learning. We also expect that it will be valuable in establishing equivalence between different Machine Learning-based models of the same phenomenon observed through different instruments and by different research groups.


title: Finding Equilibrium in Multi-Agent Games with Payoff Uncertainty

abstract: We study the problem of finding equilibrium strategies in multi-agent games with incomplete payoff information, where the payoff matrices are only known to the players up to some bounded uncertainty sets. In such games, an ex-post equilibrium characterizes equilibrium strategies that are robust to the payoff uncertainty. When the game is one-shot, we show that in zero-sum polymatrix games, an ex-post equilibrium can be computed efficiently using linear programming. We further extend the notion of ex-post equilibrium to stochastic games, where the game is played repeatedly in a sequence of stages and the transition dynamics are governed by an Markov decision process (MDP). We provide sufficient condition for the existence of an ex-post Markov perfect equilibrium (MPE). We show that under bounded payoff uncertainty, the value of any two-player zero-sum stochastic game can be computed up to a tight value interval using dynamic programming.


title: A successful search for intervening 21 cm HI absorption in galaxies at 0.4 < z <1.0 with the Australian Square Kilometre Array Pathfinder (ASKAP)

abstract: We have used the Australian Square Kilometre Array Pathfinder (ASKAP) radio telescope to search for intervening 21 cm neutral hydrogen (HI) absorption along the line of sight to 53 bright radio continuum sources. Our observations are sensitive to HI column densities typical of Damped Lyman Alpha absorbers (DLAs) in cool gas with an HI spin temperature below about 300-500 K. The six-dish Boolardy Engineering Test Array (BETA) and twelve-antenna Early Science array (ASKAP-12) covered a frequency range corresponding to redshift $0.4<z<1.0$ and $0.37<z<0.77$ respectively for the HI line. Fifty of the 53 radio sources observed have reliable optical redshifts, giving a total redshift path $\Delta z$ = 21.37. This was a spectroscopically-untargeted survey, with no prior assumptions about the location of the lines in redshift space. Four intervening HI lines were detected, two of them new. In each case, the estimated HI column density lies above the DLA limit for HI spin temperatures above 50-80 K, and we estimate a DLA number density at redshift $z\sim0.6$ of $n(z)=0.19\substack{+0.15 \\ -0.09}$. This value lies somewhat above the general trend of $n(z)$ with redshift seen in optical DLA studies. Although the current sample is small, it represents an important proof of concept for the much larger 21cm First Large Absorption Survey in HI (FLASH) project to be carried out with the full 36-antenna ASKAP telescope, probing a total redshift path $\Delta z\sim\,50,000$.


title: Broadband and giant nonreciprocity at the subwavelength scale in magnetoplasmonic materials

abstract: We unveil a previously overlooked wave propagation regime in magnetized plasmonic (gyrotropic) materials with comparable plasma and cyclotron frequencies, which enables a giant and broadband (nondispersive) nonreciprocal response. We show that this effect is due to a natural form of dispersion compensation that ultimately originates from the subtle implications of the principle of causality for gyrotropic plasmonic media, which allows the existence of a low-loss frequency window with anomalous nonmonotonic dispersion for the extraordinary mode. This is in stark contrast with conventional nongyrotropic passive materials, for which the frequency derivative of the permittivity dispersion function is always positive in low-loss regions. These findings pave the way for superior nonreciprocal components in terms of bandwidth of operation and compactness, with orders-of-magnitude reductions in size. As a relevant example, we consider indium antimonide (InSb) to theoretically demonstrate a deeply subwavelength, broadband, THz isolator operating under moderate magnetic bias and at room temperature.


title: Detecting entanglement of unknown continuous variable states with random measurements

abstract: We explore the possibility of entanglement detection in continuous variable systems by entanglement witnesses based on covariance matrices, constructible from random homodyne measurements. We propose new linear constraints characterizing the entanglement witnesses based on second moments, and use them in a semidefinite program providing the optimal entanglement test for given random measurements. We test the method on the class of squeezed vacuum states and study the efficiency of entanglement detection in general unknown covariance matrices.


title: Deep or Simple Models for Semantic Tagging? It Depends on your Data [Experiments]

abstract: Semantic tagging, which has extensive applications in text mining, predicts whether a given piece of text conveys the meaning of a given semantic tag. The problem of semantic tagging is largely solved with supervised learning and today, deep learning models are widely perceived to be better for semantic tagging. However, there is no comprehensive study supporting the popular belief. Practitioners often have to train different types of models for each semantic tagging task to identify the best model. This process is both expensive and inefficient. We embark on a systematic study to investigate the following question: Are deep models the best performing model for all semantic tagging tasks? To answer this question, we compare deep models against "simple models" over datasets with varying characteristics. Specifically, we select three prevalent deep models (i.e. CNN, LSTM, and BERT) and two simple models (i.e. LR and SVM), and compare their performance on the semantic tagging task over 21 datasets. Results show that the size, the label ratio, and the label cleanliness of a dataset significantly impact the quality of semantic tagging. Simple models achieve similar tagging quality to deep models on large datasets, but the runtime of simple models is much shorter. Moreover, simple models can achieve better tagging quality than deep models when targeting datasets show worse label cleanliness and/or more severe imbalance. Based on these findings, our study can systematically guide practitioners in selecting the right learning model for their semantic tagging task.


title: Chiral symmetry restoration in RQED at finite temperature in the supercritical coupling regime

abstract: We explore the conditions for chiral symmetry breaking in reduced (or pseudo) quantum electrodynamics at finite temperature in connection with graphene and other 2D-materials with an underlying Dirac behavior of the charge carriers. By solving the corresponding Schwinger-Dyson equation in rainbow approximation, in a Landau-like gauge and neglecting wavefunction renormalization effects, we find the need of the coupling to exceed a critical value $\alpha_c$ in order for chiral symmetry to be broken, in agreement with known results from other groups. In this supercritical regime, we add the effects of a thermal bath at temperature $T$ and find the critical values of this parameter that leads to chiral symmetry restoration.


title: Flavored axions and the flavor problem

abstract: A Peccei-Quinn~(PQ) symmetry is proposed in order to generate in the Standard Model~(SM) quark sector a realistic mass matrix ansatz with five texture-zeros. Limiting our analysis to Hermitian mass matrices, we show that this requires a minimum of 4 Higgs doublets. This model allows assigning values close to 1 for several Yukawa couplings, giving insight into the origin of the mass scales in the SM. Since the PQ charges are non-universal, the model features Flavor-Changing Neutral Currents~(FCNC) at the tree level. From the analytical expressions for the FCNC we report the allowed region in the parameter space obtained from the measurements of branching ratios of semileptonic meson decays.


title: Fully Nonlinear Equations with Applications to Grad Equations in Plasma Physics

abstract: In this paper we generalize an equation studied by Mossino and Temam, to the fully nonlinear case. This equation arises in plasma physics as an approximation to Grad equations, which were introduced by Harold Grad, to model the behavior of plasma confined in a toroidal vessel. We prove existence of a $W^{2,p}$-viscosity solution and regularity up to $C^{1,\alpha}(\overline{\Omega})$ for any $\alpha<1$(we improve this regularity near the boundary). The difficulty of this problem lays on a right hand side which involves the measure of the superlevel sets, making the problem nonlocal.


title: Evolving Graphical Planner: Contextual Global Planning for Vision-and-Language Navigation

abstract: The ability to perform effective planning is crucial for building an instruction-following agent. When navigating through a new environment, an agent is challenged with (1) connecting the natural language instructions with its progressively growing knowledge of the world; and (2) performing long-range planning and decision making in the form of effective exploration and error correction. Current methods are still limited on both fronts despite extensive efforts. In this paper, we introduce the Evolving Graphical Planner (EGP), a model that performs global planning for navigation based on raw sensory input. The model dynamically constructs a graphical representation, generalizes the action space to allow for more flexible decision making, and performs efficient planning on a proxy graph representation. We evaluate our model on a challenging Vision-and-Language Navigation (VLN) task with photorealistic images and achieve superior performance compared to previous navigation architectures. For instance, we achieve a 53% success rate on the test split of the Room-to-Room navigation task through pure imitation learning, outperforming previous navigation architectures by up to 5%.


title: Convex Hulls for Graphs of Quadratic Functions With Unit Coefficients: Even Wheels and Complete Split Graphs

abstract: We study the convex hull of the graph of a quadratic function $f(\mathbf{x})=\sum_{ij\in E}x_ix_j$, where the sum is over the edge set of a graph $G$ with vertex set $\{1,\dots,n\}$. Using an approach proposed by Gupte et al. (Discrete Optimization $\textbf{36}$, 2020, 100569), we investigate minimal extended formulations using additional variables $y_{ij}$, $1\leq i<j\leq n$, representing the products $x_ix_j$. The basic idea is to identify a set of facets of the Boolean Quadric Polytope which is sufficient for characterizing the convex hull for the given graph. Our main results are extended formulations for the cases that the underlying graph $G$ is either an even wheel or a complete split graph.


title: Hardware Implementation of Deep Network Accelerators Towards Healthcare and Biomedical Applications

abstract: The advent of dedicated Deep Learning (DL) accelerators and neuromorphic processors has brought on new opportunities for applying both Deep and Spiking Neural Network (SNN) algorithms to healthcare and biomedical applications at the edge. This can facilitate the advancement of medical Internet of Things (IoT) systems and Point of Care (PoC) devices. In this paper, we provide a tutorial describing how various technologies including emerging memristive devices, Field Programmable Gate Arrays (FPGAs), and Complementary Metal Oxide Semiconductor (CMOS) can be used to develop efficient DL accelerators to solve a wide variety of diagnostic, pattern recognition, and signal processing problems in healthcare. Furthermore, we explore how spiking neuromorphic processors can complement their DL counterparts for processing biomedical signals. The tutorial is augmented with case studies of the vast literature on neural network and neuromorphic hardware as applied to the healthcare domain. We benchmark various hardware platforms by performing a sensor fusion signal processing task combining electromyography (EMG) signals with computer vision. Comparisons are made between dedicated neuromorphic processors and embedded AI accelerators in terms of inference latency and energy. Finally, we provide our analysis of the field and share a perspective on the advantages, disadvantages, challenges, and opportunities that various accelerators and neuromorphic processors introduce to healthcare and biomedical domains.


title: Search for the lepton flavor violating decay $\tau$ $\to$ 3$\mu$ in proton-proton collisions at $\sqrt{s} =$ 13 TeV

abstract: Results are reported from a search for the lepton flavor violating decay $\tau$ $\to$ 3$\mu$ in proton-proton collisions at $\sqrt{s} =$ 13 TeV. The data sample corresponds to an integrated luminosity of 33.2 fb$^{-1}$ recorded by the CMS experiment at the LHC in 2016. The search exploits $\tau$ leptons produced in both W boson and heavy-flavor hadron decays. No significant excess above the expected background is observed. An upper limit on the branching fraction $\mathcal{B}(\tau$ $\to$ 3$\mu)$ of 8.0 $\times$ 10$^{-8}$ at 90% confidence level is obtained, with an expected upper limit of 6.9 $\times$ 10$^{-8}$.


title: New cosmological bounds on axions in the XENON1T window

abstract: Motivated by a possible $\sim$ eV-mass solar axion explanation to excess events recently detected by the XENON1T experiment, I revisit and update cosmological constraints on axions in this mass range. I find that of the allowed XENON1T mass window (0.1 - 4.1 eV for DFSZ axions and 46 - 56 eV for KSVZ axions), only 0.1 - 0.35 eV remains viable at 95% confidence given current cosmological probes. If a 0.35 eV DFSZ axion existed, it would be detectable at ${\sim}7\sigma$ via two independent physical effects with the next-generation CMB-S4 experiment. Conversely, even a combination of CMB-S4 with future DESI measurements falls just short of guaranteeing a 0.1 eV-mass axion can be detected or ruled out. A future limit of $\Delta N_{\rm eff}<0.027$ could rule out any generic axion-like particle across a wide range of masses as long as the reheating temperature is not too low, or alternatively, a future cosmological detection of such an axion-like particle could become the tightest existing observational lower bound on the reheating temperature.


title: Generating W states with braiding operators

abstract: Braiding operators can be used to create entangled states out of product states, thus establishing a correspondence between topological and quantum entanglement. This is well-known for maximally entangled Bell and GHZ states and their equivalent states under Stochastic Local Operations and Classical Communication, but so far a similar result for W states was missing. Here we use generators of extraspecial 2-groups to obtain the W state in a four-qubit space and partition algebras to generate the W state in a three-qubit space. We also present a unitary generalized Yang-Baxter operator that embeds the W$_n$ state in a $(2n-1)$-qubit space.


title: Deep Patch-based Human Segmentation

abstract: 3D human segmentation has seen noticeable progress in re-cent years. It, however, still remains a challenge to date. In this paper, weintroduce a deep patch-based method for 3D human segmentation. Wefirst extract a local surface patch for each vertex and then parameterizeit into a 2D grid (or image). We then embed identified shape descriptorsinto the 2D grids which are further fed into the powerful 2D Convolu-tional Neural Network for regressing corresponding semantic labels (e.g.,head, torso). Experiments demonstrate that our method is effective inhuman segmentation, and achieves state-of-the-art accuracy.


title: Local Lipschitz bounds for solutions to certain singular elliptic equations involving one-Laplacian

abstract: In this paper local Lipschitz regularity of weak solutions to certain singular elliptic equations involving one-Laplacian is studied. Equations treated here also contains another well-behaving elliptic operator such as $p$-Laplacian with $1<p<\infty$. The problem is that one-Laplacian is too singular on degenerate points, what is often called facet, which makes it difficult to obtain even Lipschitz regularity of weak solutions. This difficulty is overcome by making suitable approximation schemes, and by avoiding analysis on facet for approximated solutions. The key estimate is a local a priori uniform Lipschitz estimate for classical solutions to regularized equations, which is proved by Moser's iteration. Another local a priori uniform Lipschitz bounds can also be obtained by De Giorgi's truncation. Proofs of local Lipschitz estimates in this paper are rather classical and elementary in the sense that nonlinear potential estimates are not used at all.


title: Quasi-Periodic WaveNet: An Autoregressive Raw Waveform Generative Model with Pitch-dependent Dilated Convolution Neural Network

abstract: In this paper, a pitch-adaptive waveform generative model named Quasi-Periodic WaveNet (QPNet) is proposed to improve the limited pitch controllability of vanilla WaveNet (WN) using pitch-dependent dilated convolution neural networks (PDCNNs). Specifically, as a probabilistic autoregressive generation model with stacked dilated convolution layers, WN achieves high-fidelity audio waveform generation. However, the pure-data-driven nature and the lack of prior knowledge of audio signals degrade the pitch controllability of WN. For instance, it is difficult for WN to precisely generate the periodic components of audio signals when the given auxiliary fundamental frequency ($F_{0}$) features are outside the $F_{0}$ range observed in the training data. To address this problem, QPNet with two novel designs is proposed. First, the PDCNN component is applied to dynamically change the network architecture of WN according to the given auxiliary $F_{0}$ features. Second, a cascaded network structure is utilized to simultaneously model the long- and short-term dependencies of quasi-periodic signals such as speech. The performances of single-tone sinusoid and speech generations are evaluated. The experimental results show the effectiveness of the PDCNNs for unseen auxiliary $F_{0}$ features and the effectiveness of the cascaded structure for speech generation.


title: Constraints on circum-galactic media from Sunyaev-Zel'dovich effects and X-ray data

abstract: We use observational measurements of thermal and kinetic Sunyaev-Zel'dovich effects, as well as soft X-ray emission associated with galaxy groups to constrain the gas density and temperature in the circumgalactic medium (CGM) for dark matter halos with masses above $10^{12.5}M_{\odot}$. A number of generic models are used together with a Bayesian scheme to make model inferences. We find that gas with a single temperature component cannot provide a consistent model to match the observational data. A simple two-phase model assuming a hot component and an ionized warm component can accommodate all the three observations. The total amount of the gas in individual halos is inferred to be comparable to the universal baryon fraction corresponding to the halo mass. The inferred temperature of the hot component is comparable to the halo virial temperature. The fraction of the hot component increases from $(15-40)\%$ for $10^{12.5}{M}_\odot$ halos to $(40-60)\%$ for $10^{14.5}{M}_\odot$ halos, where the ranges reflect uncertainties in the assumed gas density profile. Our results suggest that a significant fraction of the halo gas is in a non-thermalized component with temperature much lower than the virial temperature.


title: A Computational Separation between Private Learning and Online Learning

abstract: A recent line of work has shown a qualitative equivalence between differentially private PAC learning and online learning: A concept class is privately learnable if and only if it is online learnable with a finite mistake bound. However, both directions of this equivalence incur significant losses in both sample and computational efficiency. Studying a special case of this connection, Gonen, Hazan, and Moran (NeurIPS 2019) showed that uniform or highly sample-efficient pure-private learners can be time-efficiently compiled into online learners. We show that, assuming the existence of one-way functions, such an efficient conversion is impossible even for general pure-private learners with polynomial sample complexity. This resolves a question of Neel, Roth, and Wu (FOCS 2019).


title: Observation of inverted band structure in topological Dirac-semimetal candidate CaAuAs

abstract: We have performed high-resolution angle-resolved photoemission spectroscopy of ternary pnictide CaAuAs which is predicted to be a three-dimensional topological Dirac semimetal (TDS). By accurately determining the bulk-band structure, we have revealed the coexistence of three-dimensional and quasi-two-dimensional Fermi surfaces with dominant hole carriers. The band structure around the Brillouin-zone center is characterized by an energy overlap between hole and electron pockets, in excellent agreement with first-principles band-structure calculations. This indicates the occurrence of bulk-band inversion, supporting the TDS state in CaAuAs. Because of the high tunability in the chemical composition besides the TDS nature, CaAuAs provides a precious opportunity for investigating the quantum phase transition from TDS to other exotic topological phases.


title: To Filter Prune, or to Layer Prune, That Is The Question

abstract: Recent advances in pruning of neural networks have made it possible to remove a large number of filters or weights without any perceptible drop in accuracy. The number of parameters and that of FLOPs are usually the reported metrics to measure the quality of the pruned models. However, the gain in speed for these pruned models is often overlooked in the literature due to the complex nature of latency measurements. In this paper, we show the limitation of filter pruning methods in terms of latency reduction and propose LayerPrune framework. LayerPrune presents a set of layer pruning methods based on different criteria that achieve higher latency reduction than filter pruning methods on similar accuracy. The advantage of layer pruning over filter pruning in terms of latency reduction is a result of the fact that the former is not constrained by the original model's depth and thus allows for a larger range of latency reduction. For each filter pruning method we examined, we use the same filter importance criterion to calculate a per-layer importance score in one-shot. We then prune the least important layers and fine-tune the shallower model which obtains comparable or better accuracy than its filter-based pruning counterpart. This one-shot process allows to remove layers from single path networks like VGG before fine-tuning, unlike in iterative filter pruning, a minimum number of filters per layer is required to allow for data flow which constraint the search space. To the best of our knowledge, we are the first to examine the effect of pruning methods on latency metric instead of FLOPs for multiple networks, datasets and hardware targets. LayerPrune also outperforms handcrafted architectures such as Shufflenet, MobileNet, MNASNet and ResNet18 by 7.3%, 4.6%, 2.8% and 0.5% respectively on similar latency budget on ImageNet dataset.


title: The compressible Euler equations in a physical vacuum: a comprehensive Eulerian approach

abstract: This article is concerned with the local well-posedness problem for the compressible Euler equations in gas dynamics. For this system we consider the free boundary problem which corresponds to a physical vacuum. Despite the clear physical interest in this system, the prior work on this problemis limited to Lagrangian coordinates, in high regularity spaces. Instead, the objective of the present work is to provide a new, fully Eulerian approach to this problem, which provides a complete, Hadamard style well-posedness theory for this problem in low regularity Sobolev spaces. In particular we give new proofs for both existence, uniqueness, and continuous dependence on the data with sharp, scale invariant energy estimates, and continuation criterion.


title: Hilbert series in the category of trees with contractions

abstract: We consider Hilbert series associated to modules over various categories of trees. Using the technology of Sam and Snowden, we show that these Hilbert series must be algebraic. We then apply these technical theorems to prove facts about certain natural generating functions associated to trees.


title: An Asymptotically Optimal Multi-Armed Bandit Algorithm and Hyperparameter Optimization

abstract: The evaluation of hyperparameters, neural architectures, or data augmentation policies becomes a critical model selection problem in advanced deep learning with a large hyperparameter search space. In this paper, we propose an efficient and robust bandit-based algorithm called Sub-Sampling (SS) in the scenario of hyperparameter search evaluation. It evaluates the potential of hyperparameters by the sub-samples of observations and is theoretically proved to be optimal under the criterion of cumulative regret. We further combine SS with Bayesian Optimization and develop a novel hyperparameter optimization algorithm called BOSS. Empirical studies validate our theoretical arguments of SS and demonstrate the superior performance of BOSS on a number of applications, including Neural Architecture Search (NAS), Data Augmentation (DA), Object Detection (OD), and Reinforcement Learning (RL).


title: A new data assimilation method of recovering turbulent flow field at high-Reynolds numbers for turbulence machine learning

abstract: This paper proposes a new data assimilation method for recovering high fidelity turbulent flow field around airfoil at high Reynolds numbers based on experimental data, which is called Proper Orthogonal Decomposition Inversion (POD-Inversion) data assimilation method. Aiming at the flows including shock wave discontinuities or separated flows at high angle of attack, the proposed method can reconstruct high-fidelity turbulent flow field combining with experimental distributed force coefficients. We firstly perform the POD analysis to the turbulent eddy viscosity fields computed by SA model and obtain the base POD modes. Then optimized the POD coefficients by global optimization algorithm coupling with the Navier-Stokes equations solver. The high-fidelity turbulent flied are recovered by several main modes, which can dramatically reduce the dimensions of the system. The effectiveness of the method is verified by the cases of transonic flow around the RAE2822 airfoil at high Reynolds numbers and the separated flow at high angles of attack. The results demonstrate that the proposed assimilation method can recover the turbulent flow field which optimally match the experimental data, and significantly reduce the error of pressure coefficients. The proposed data assimilation method can offer high-fidelity field data for turbulent model based on machine learning.


title: A critical survey on the kinetic assays of DNA polymerase fidelity from a new theoretical perspective

abstract: The high fidelity of DNA polymerase is critical for the faithful replication of genomic DNA. Several approaches were proposed to quantify the fidelity of DNA polymerase. Direct measurements of the error frequency of the replication products definitely give the true fidelity but turn out very hard to implement. Two biochemical kinetic approaches, the steady-state assay and the transient-state assay, were then suggested and widely adopted. In these assays, the error frequency is indirectly estimated by using the steady-state or the transient-state kinetic theory combined with the measured kinetic rates. However, whether these indirectly estimated fidelities are equivalent to the true fidelity has never been clarified theoretically, and in particular there are different strategies to quantify the proofreading efficiency of DNAP but often lead to inconsistent results. The reason for all these confusions is that it's mathematically challenging to formulate a rigorous and general theory of the true fidelity. Recently we have succeeded to establish such a theoretical framework. In this paper, we develop this theory to make a comprehensive examination on the theoretical foundation of the kinetic assays and the relation between fidelities obtained by different methods. We conclude that while the steady-state assay and the transient-state assay can always measure the true fidelity of exonuclease-deficient DNA polymerases, they only do so for exonuclease-efficient DNA polymerases conditionally (the proper way to use these assays to quantify the proofreading efficiency is also suggested). We thus propose a new kinetic approach, the single-molecule assay, which indirectly but precisely characterizes the true fidelity of either exonuclease-deficient or exonuclease-efficient DNA polymerases.


title: iRDRC: An Intelligent Real-time Dual-functional Radar-Communication System for Automotive Vehicles

abstract: This letter introduces an intelligent Real-time Dual-functional Radar-Communication (iRDRC) system for autonomous vehicles (AVs). This system enables an AV to perform both radar and data communications functions to maximize bandwidth utilization as well as significantly enhance safety. In particular, the data communications function allows the AV to transmit data, e.g., of current traffic, to edge computing systems and the radar function is used to enhance the reliability and reduce the collision risks of the AV, e.g., under bad weather conditions. The problem of the iRDRC is to decide when to use the communication mode or the radar mode to maximize the data throughput while minimizing the miss detection probability of unexpected events given the uncertainty of surrounding environment. To solve the problem, we develop a deep reinforcement learning algorithm that allows the AV to quickly obtain the optimal policy without requiring any prior information about the environment. Simulation results show that the proposed scheme outperforms baseline schemes in terms of data throughput, miss detection probability, and convergence rate.


title: Illuminating Mario Scenes in the Latent Space of a Generative Adversarial Network

abstract: Generative adversarial networks (GANs) are quickly becoming a ubiquitous approach to procedurally generating video game levels. While GAN generated levels are stylistically similar to human-authored examples, human designers often want to explore the generative design space of GANs to extract interesting levels. However, human designers find latent vectors opaque and would rather explore along dimensions the designer specifies, such as number of enemies or obstacles. We propose using state-of-the-art quality diversity algorithms designed to optimize continuous spaces, i.e. MAP-Elites with a directional variation operator and Covariance Matrix Adaptation MAP-Elites, to efficiently explore the latent space of a GAN to extract levels that vary across a set of specified gameplay measures. In the benchmark domain of Super Mario Bros, we demonstrate how designers may specify gameplay measures to our system and extract high-quality (playable) levels with a diverse range of level mechanics, while still maintaining stylistic similarity to human authored examples. An online user study shows how the different mechanics of the automatically generated levels affect subjective ratings of their perceived difficulty and appearance.


title: Towards Cross-Granularity Few-Shot Learning: Coarse-to-Fine Pseudo-Labeling with Visual-Semantic Meta-Embedding

abstract: Few-shot learning aims at rapidly adapting to novel categories with only a handful of samples at test time, which has been predominantly tackled with the idea of meta-learning. However, meta-learning approaches essentially learn across a variety of few-shot tasks and thus still require large-scale training data with fine-grained supervision to derive a generalized model, thereby involving prohibitive annotation cost. In this paper, we advance the few-shot classification paradigm towards a more challenging scenario, i.e., cross-granularity few-shot classification, where the model observes only coarse labels during training while is expected to perform fine-grained classification during testing. This task largely relieves the annotation cost since fine-grained labeling usually requires strong domain-specific expertise. To bridge the cross-granularity gap, we approximate the fine-grained data distribution by greedy clustering of each coarse-class into pseudo-fine-classes according to the similarity of image embeddings. We then propose a meta-embedder that jointly optimizes the visual- and semantic-discrimination, in both instance-wise and coarse class-wise, to obtain a good feature space for this coarse-to-fine pseudo-labeling process. Extensive experiments and ablation studies are conducted to demonstrate the effectiveness and robustness of our approach on three representative datasets.


title: Learning Object Depth from Camera Motion and Video Object Segmentation

abstract: Video object segmentation, i.e., the separation of a target object from background in video, has made significant progress on real and challenging videos in recent years. To leverage this progress in 3D applications, this paper addresses the problem of learning to estimate the depth of segmented objects given some measurement of camera motion (e.g., from robot kinematics or vehicle odometry). We achieve this by, first, introducing a diverse, extensible dataset and, second, designing a novel deep network that estimates the depth of objects using only segmentation masks and uncalibrated camera movement. Our data-generation framework creates artificial object segmentations that are scaled for changes in distance between the camera and object, and our network learns to estimate object depth even with segmentation errors. We demonstrate our approach across domains using a robot camera to locate objects from the YCB dataset and a vehicle camera to locate obstacles while driving.


title: Constraining Primordial Black Hole Fraction at the Galactic Centre using radio observational data

abstract: Recent gamma-ray and cosmic-ray observations have put strong constraints on the amount of primordial black holes (PBHs) in our universe. In this article, we use the archival radio data of the inner Galactic Centre to constrain the PBH to dark matter ratio for three different PBH mass distributions including monochromatic, log-normal and power-law. We show that the amount of PBHs only constitutes a very minor component of dark matter at the Galactic Centre for a large parameter space.


title: Charm fragmentation functions in the Nambu-Jona-Lasinio model

abstract: We derive the fragmentation function (FF), which describes the probability for a charm quark to emit a $D$ meson with a certain momentum fraction, in the Nambu-Jona-Lasinio (NJL) model. The corresponding elementary FF is calculated with the quark-meson coupling determined in the NJL model involving charm quarks. The FF in the infinite momentum frame is constructed through the jet process governed by the elementary FF, and then evolved to the charm scale, at which it is defined. To prepare the FF suitable for an analysis of $D$ meson production at CLEO, we further match the above FF to that in the finite momentum frame at one loop in QCD. It is shown that the charm quark FF including the finite momentum effects leads to theoretical results in agreement with the CLEO data.


title: Transport and deposition of dilute microparticles in turbulent thermal convection

abstract: We analyze the transport and deposition behavior of dilute microparticles in turbulent Rayleigh-B\'enard convection. Two-dimensional direct numerical simulations were carried out for the Rayleigh number ($Ra$) of $10^{8}$ and the Prandtl number ($Pr$) of 0.71 (corresponding to the working fluids of air). The Lagrangian point particle model was used to describe the motion of microparticles in the turbulence. Our results show that the suspended particles are homogeneously distributed in the turbulence for the Stokes number ($St$) less than $10^{-3}$, and they tend to cluster into bands for $10^{-3} \lesssim St \lesssim 10^{-2}$. At even larger $St$, the microparticles will quickly sediment in the convection. We also calculate the mean-square displacement (MSD) of the particle's trajectories. At short time intervals, the MSD exhibits a ballistic regime, and it is isotropic in vertical and lateral directions; at longer time intervals, the MSD reflects a confined motion for the particles, and it is anisotropic in different directions. We further obtained a phase diagram of the particle deposition positions on the wall, and we identified three deposition states depending on the particle's density and diameter. An interesting finding is that the dispersed particles preferred to deposit on the vertical wall where the hot plumes arise, which is verified by tilting the cell and altering the rotation direction of the large-scale circulation.


title: Capacity Improvement in Wideband Reconfigurable Intelligent Surface-Aided Cell-Free Network

abstract: Thanks to the strong ability against the inter-cell interference, cell-free network has been considered as a promising technique to improve the network capacity of future wireless systems. However, for further capacity enhancement, it requires to deploy more base stations (BSs) with high cost and power consumption. To address the issue, inspired by the recently proposed technique called reconfigurable intelligent surface (RIS), we propose the concept of RIS-aided cell-free network to improve the network capacity with low cost and power consumption. Then, for the proposed RIS-aided cell-free network in the typical wideband scenario, we formulate the joint precoding design problem at the BSs and RISs to maximize the network capacity. Due to the non-convexity and high complexity of the formulated problem, we develop an alternating optimization algorithm to solve this challenging problem. Note that most of the considered scenarios in existing works are special cases of the general scenario in this paper, and the proposed joint precoding framework can also serve as a general solution to maximize the capacity in most of existing RIS-aided scenarios. Finally, simulation results verify that, compared with the conventional cell-free network, the network capacity of the proposed scheme can be improved significantly.


title: Root estimation in Galton-Watson trees

abstract: Given only the free-tree structure of a tree, the root estimation problem asks if one can guess which of the free tree's nodes is the root of the original tree. We determine the maximum-likelihood estimator for the root of a free tree when the underlying tree is a size-conditioned Galton-Watson tree and calculate its probability of being correct.


title: Magnon-bipolar carrier drag thermopower in antiferromagnetic/ferromagnetic semiconductors: theoretical formulation and experimental evidence

abstract: Quantized spin-wave known as magnon, a bosonic quasiparticle, can drag electrons or holes via sd exchange interaction and boost the thermopower over the conventional diffusive thermopower. P-type magnon-drag thermopower has been observed in both ferromagnetic and antiferromagnetic metallic and degenerate semiconductors. However, it has been less reported for intrinsic or n-type magnetic semiconductors; therefore, the impact of magnon-bipolar carrier drag on thermopower has remained unexplored. Here, a theoretical model for magnon-bipolar carrier drag thermopower is derived based on the magnon-carrier interaction lifetimes. The model predicts that the bipolar carrier drag thermopower becomes independent of both the carrier and magnon relaxation times. A proof of concept experiment is presented that confirms this prediction. We also report the observation of magnon-carrier drag thermopower in n-type and intrinsic ferromagnetic semiconductors experimentally. The p-type antiferromagnetic MnTe is doped with different amounts of Cr to produce non-degenerate and n-type semiconductors of various carrier concentrations. Cr dopants have a donor nature and create ferromagnetic-antiferromagnetic clusters due to the Cr3+ oxidation state. Heat capacity measurements confirm the presence of magnons in Cr-doped MnTe. It is shown that the magnon-drag thermopower is significantly reduced for 3%-5% Cr-doped samples due to bipolar drag effects and becomes negative for 14% and 20% Cr-doped MnTe due to dominant magnon-electron drag thermopower.


title: Batch-level Experience Replay with Review for Continual Learning

abstract: Continual learning is a branch of deep learning that seeks to strike a balance between learning stability and plasticity. The CVPR 2020 CLVision Continual Learning for Computer Vision challenge is dedicated to evaluating and advancing the current state-of-the-art continual learning methods using the CORe50 dataset with three different continual learning scenarios. This paper presents our approach, called Batch-level Experience Replay with Review, to this challenge. Our team achieved the 1'st place in all three scenarios out of 79 participated teams. The codebase of our implementation is publicly available at https://github.com/RaptorMai/CVPR20_CLVision_challenge


title: Fast Real-time Counterfactual Explanations

abstract: Counterfactual explanations are considered, which is to answer {\it why the prediction is class A but not B.} Different from previous optimization based methods, an optimization-free Fast ReAl-time Counterfactual Explanation (FRACE) algorithm is proposed benefiting from the development of multi-domain image to image translation algorithms. Built from starGAN, a transformer is trained as a residual generator conditional on a classifier constrained under a proposal perturbation loss which maintains the content information of the query image, but just the class-specific semantic information is changed. The transformer can transfer the query image to any counterfactual class, and during inference, our explanation can be generated by it only within a forward time. It is fast and can satisfy the real-time practical application. Because of the adversarial training of GAN, our explanation is also more realistic compared to other counterparts. The experimental results demonstrate that our proposal is better than the existing state of the art in terms of quality and speed.


title: NeuralExplorer: State Space Exploration of Closed Loop Control Systems Using Neural Networks

abstract: In this paper, we propose a framework for performing state space exploration of closed loop control systems. Our approach involves approximating sensitivity and a newly introduced notion of inverse sensitivity by a neural network. We show how the approximation of sensitivity and inverse sensitivity can be used for computing estimates of the reachable set. We then outline algorithms for performing state space exploration by generating trajectories that reach a neighborhood. We demonstrate the effectiveness of our approach by applying it not only to standard linear and nonlinear dynamical systems, but also to nonlinear hybrid systems and also neural network based feedback control systems.


title: Event-based Signal Processing for Radioisotope Identification

abstract: This paper identifies the problem of unnecessary high power overhead of the conventional frame-based radioisotope identification process and proposes an event-based signal processing process to address the problem established. It also presents the design flow of the neuromorphic processor.


title: Fast Video Object Segmentation With Temporal Aggregation Network and Dynamic Template Matching

abstract: Significant progress has been made in Video Object Segmentation (VOS), the video object tracking task in its finest level. While the VOS task can be naturally decoupled into image semantic segmentation and video object tracking, significantly much more research effort has been made in segmentation than tracking. In this paper, we introduce "tracking-by-detection" into VOS which can coherently integrate segmentation into tracking, by proposing a new temporal aggregation network and a novel dynamic time-evolving template matching mechanism to achieve significantly improved performance. Notably, our method is entirely online and thus suitable for one-shot learning, and our end-to-end trainable model allows multiple object segmentation in one forward pass. We achieve new state-of-the-art performance on the DAVIS benchmark without complicated bells and whistles in both speed and accuracy, with a speed of 0.14 second per frame and J&F measure of 75.9% respectively.


title: The Kuramoto model in presence of additional interactions that break rotational symmetry

abstract: The Kuramoto model serves as a paradigm to study the phenomenon of spontaneous collective synchronization. We study here a nontrivial generalization of the Kuramoto model by including an interaction that breaks explicitly the rotational symmetry of the model. In an inertial frame (e.g., the laboratory frame), the Kuramoto model does not allow for a stationary state, that is, a state with time-independent value of the so-called Kuramoto (complex) synchronization order parameter $z\equiv re^{i\psi}$; Note that a time-independent $z$ implies $r$ and $\psi$ both time independent, with the latter fact corresponding to a state in which $\psi$ rotates at zero frequency (no rotation). In this backdrop, we ask: Does the introduction of the symmetry-breaking term suffice to allow for the existence of a stationary state in the laboratory frame? Compared to the original model, we reveal a rather rich phase diagram of the resulting model, with the existence of both stationary and standing wave phases. While in the former the synchronization order parameter $r$ has a long-time value that is time independent, one has in the latter an oscillatory behavior of the order parameter as a function of time that nevertheless yields a non-zero and time-independent time average. Our results are based on numerical integration of the dynamical equations as well as an exact analysis of the dynamics by invoking the so-called Ott-Antonsen ansatz that allows to derive a reduced set of time-evolution equations for the order parameter.


title: You Need to Calm Down: Calmness Regularity for a Class of Seminorm Optimization Problems

abstract: Compressed sensing involves solving a minimization problem with objective function $\Omega(\boldsymbol{x}) = \|\boldsymbol{x}\|_1$ and linear constraints $\boldsymbol{A} \boldsymbol{x} = \boldsymbol{b}$. Previous work has explored robustness to errors in $\boldsymbol{A}$ and $\boldsymbol{b}$ under special assumptions. Motivated by these results, we explore robustness to errors in $\boldsymbol{A}$ for a wider class of objective functions $\Omega$ and for a more general setting, where the solution may not be unique. Similar results for errors in $\boldsymbol{b}$ are known and easier to prove. More precisely, for a seminorm $\Omega(\boldsymbol{x})$ with a polyhedral unit ball, we prove that the set-valued map $S(\boldsymbol{A}) = \arg \min_{\boldsymbol{A} \boldsymbol{x} = \boldsymbol{b}} \Omega(\boldsymbol{x})$ is calm in $\boldsymbol{A}$, where calmness is a kind of local Lipschitz regularity.


title: A Unified Linear Speedup Analysis of Federated Averaging and Nesterov FedAvg

abstract: Federated learning (FL) learns a model jointly from a set of participating devices without sharing each other's privately held data. The characteristics of non-i.i.d. data across the network, low device participation, high communication costs, and the mandate that data remain private bring challenges in understanding the convergence of FL algorithms, particularly regarding how convergence scales with the number of participating devices. In this paper, we focus on Federated Averaging (FedAvg), one of the most popular and effective FL algorithms in use today, as well as its Nesterov accelerated variant, and conduct a systematic study of how their convergence scale with the number of participating devices under non-i.i.d. data and partial participation in convex settings. We provide a unified analysis that establishes convergence guarantees for FedAvg under strongly convex, convex, and overparameterized strongly convex problems. We show that FedAvg enjoys linear speedup in each case, although with different convergence rates and communication efficiencies. For strongly convex and convex problems, we also characterize the corresponding convergence rates for the Nesterov accelerated FedAvg algorithm, which are the first linear speedup guarantees for momentum variants of FedAvg in convex settings. Empirical studies of the algorithms in various settings have supported our theoretical results.


title: Discrete diffusion semigroups associated with Dunkl-Jacobi and exceptional Jacobi polynomials

abstract: Some weighted inequalities for the maximal operator with respect to the discrete diffusion semigroups associated with exceptional Jacobi and Dunkl-Jacobi polynomials are given. This setup allows to extend the corresponding results obtained for discrete heat semigroup recently to richer class of differential-difference operators.


title: How Does GAN-based Semi-supervised Learning Work?

abstract: Generative adversarial networks (GANs) have been widely used and have achieved competitive results in semi-supervised learning. This paper theoretically analyzes how GAN-based semi-supervised learning (GAN-SSL) works. We first prove that, given a fixed generator, optimizing the discriminator of GAN-SSL is equivalent to optimizing that of supervised learning. Thus, the optimal discriminator in GAN-SSL is expected to be perfect on labeled data. Then, if the perfect discriminator can further cause the optimization objective to reach its theoretical maximum, the optimal generator will match the true data distribution. Since it is impossible to reach the theoretical maximum in practice, one cannot expect to obtain a perfect generator for generating data, which is apparently different from the objective of GANs. Furthermore, if the labeled data can traverse all connected subdomains of the data manifold, which is reasonable in semi-supervised classification, we additionally expect the optimal discriminator in GAN-SSL to also be perfect on unlabeled data. In conclusion, the minimax optimization in GAN-SSL will theoretically output a perfect discriminator on both labeled and unlabeled data by unexpectedly learning an imperfect generator, i.e., GAN-SSL can effectively improve the generalization ability of the discriminator by leveraging unlabeled information.


title: Electric Currents through J-shaped and Non-J-shaped Flare Ribbons

abstract: Recently solar flares exhibiting a double J-shaped ribbons in the lower solar atmosphere have been paid increasing attention in the context of extending the two-dimensional standard flare model to three dimensions, as motivated by the spatial correlation between photospheric current channels and flare ribbons. Here we study the electric currents through the photospheric area swept by flare ribbons (termed synthesized ribbon area or SRA), with a sample of 71 two-ribbon flares, of which 36 are J-shaped. Electric currents flowing through one ribbon are highly correlated with those through the other, therefore belonging to the same current system. The non-neutrality factor of this current system is independent of the flare magnitude, implying that both direct and return currents participate in flares. J-shaped flares are distinct from non-J-shaped flares in the following aspects: 1) Electric current densities within J-shaped SRA are significantly smaller than those within non-J-shaped SRA, but J-shaped SRA and its associated magnetic flux are also significantly larger. 2) Electric currents through SRA are positively correlated with the flare magnitude, but J-shaped flares show stronger correlation than non-J-shaped flares. 3) The majority (75\%) of J-shaped flares are eruptive, while the majority (86\%) of non-J-shaped flares are confined; accordingly, hosting active regions of J-shaped flares are more likely to be sigmoidal than non-J-shaped flares. Thus, J-shaped flares constitute a distinct subset of two-ribbon flares, probably the representative of eruptive ones. Further, we found that combining SRA and its associated magnetic flux has the potential to differentiate eruptive from confined flares.


title: Long-Term Planning with Deep Reinforcement Learning on Autonomous Drones

abstract: In this paper, we study a long-term planning scenario that is based on drone racing competitions held in real life. We conducted this experiment on a framework created for "Game of Drones: Drone Racing Competition" at NeurIPS 2019. The racing environment was created using Microsoft's AirSim Drone Racing Lab. A reinforcement learning agent, a simulated quadrotor in our case, has trained with the Policy Proximal Optimization(PPO) algorithm was able to successfully compete against another simulated quadrotor that was running a classical path planning algorithm. Agent observations consist of data from IMU sensors, GPS coordinates of drone obtained through simulation and opponent drone GPS information. Using opponent drone GPS information during training helps dealing with complex state spaces, serving as expert guidance allows for efficient and stable training process. All experiments performed in this paper can be found and reproduced with code at our GitHub repository


title: Ferromagnetic Weyl Fermions in Two-Dimensional Layered Electride Gd$_2$C

abstract: Recently, two-dimensional layered electrides have emerged as a new class of materials which possess anionic electron layers in the interstitial spaces between cationic layers. Here, based on first-principles calculations, we discover a time-reversal-symmetry-breaking Weyl semimetal phase in a unique two-dimensional layered ferromagnetic (FM) electride Gd$_2$C. It is revealed that the crystal field mixes the interstitial electron states and Gd 5$d$ orbitals near the Fermi energy to form band inversions. Meanwhile, the FM order induces two spinful Weyl nodal lines (WNLs), which are converted into multiple pairs of Weyl nodes through spin-orbit coupling. Further, we not only identify Fermi-arc surface states connecting the Weyl nodes but also predict a large intrinsic anomalous Hall conductivity due to the Berry curvature produced by the gapped WNLs. Our findings demonstrate the existence of Weyl fermions in the room-temperature FM electride Gd$_2$C, therefore offering a new platform to investigate the intriguing interplay between electride materials and magnetic Weyl physics.


title: Evidence for $X(3872)\rightarrow J/\psi \pi^+\pi^-$ produced in single-tag two-photon interactions

abstract: We report the first evidence for $X(3872)$ production in two-photon interactions by tagging either the electron or the position in the final state, exploring the highly virtual photon region. The search is performed in $e^+e^- \rightarrow e^+e^-J/\psi\pi^+\pi^-$, using 825 fb$^{-1}$ of data collected by the Belle detector operated at the KEKB $e^+e^-$ collider. We observe three $X(3872)$ candidates with an expected background of $0.11\pm 0.10$ events, with a significance of 3.2$\sigma$. We obtain an estimated value for $\tilde{\Gamma}_{\gamma\gamma}{\cal B}(X(3872)\rightarrow J/\psi\pi^+\pi^-$) assuming the $Q^2$ dependence predicted by a $c\bar{c}$ meson model, where $-Q^2$ is the invariant mass-squared of the virtual photon. No $X(3915)\rightarrow J/\psi\pi^+\pi^-$ candidates are found.


title: Projective and Whittaker functors on category $\mathcal{O}$

abstract: We show that the Whittaker functor on a regular block of the BGG-category $\mathcal{O}$ of a semisimple complex Lie algebra can be obtained by composing a translation to the wall functor with Soergel and Mili\v{c}i\'{c}'s equivalence between the category of Whittaker modules and a singular block of $\mathcal{O}$. We show that the Whittaker functor is a quotient functor that commutes with all projective functors and endomorphisms between them.


title: From Heun Class Equations to Painlev\'e Equations

abstract: In the first part of our paper we discuss linear 2nd order differential equations in the complex domain, especially Heun class equations, that is, the Heun equation and its confluent cases. The second part of our paper is devoted to Painlev\'e I-VI equations. Our philosophy is to treat these families of equations in a unified way. This philosophy works especially well for Heun class equations. We discuss its classification into 5 supertypes, subdivided into 10 types (not counting trivial cases). We also introduce in a unified way deformed Heun class equations, which contain an additional nonlogarithmic singularity. We show that there is a direct relationship between deformed Heun class equations and all Painlev\'e equations. In particular, Painlev\'e equations can be also divided into 5 supertypes, and subdivided into 10 types. This relationship is not so easy to describe in a completely unified way, because the choice of the ''time variable'' may depend on the type. We describe unified treatments for several possible ''time variables''.


title: Differential operators on Hitchin variety

abstract: We introduce the notion of Hitchin variety over $\C$. Let $L$ be a holomorphic line bundle over a Hitchin variety $X$. We investigate the space of all global sections of sheaf of differential operators $\cat{D}^k (L)$ and symmetric powers of sheaf of first order differential operators $\cat{S}^k(\cat{D}^1 (L))$ over $X$ and show that for a projective Hithcin variety both the spaces are one dimensional. As an application, we show that the space $\cat{C}(L)$ of holomorphic connections on $L$ does not admit any non-constant regular function.


title: M-Evolve: Structural-Mapping-Based Data Augmentation for Graph Classification

abstract: Graph classification, which aims to identify the category labels of graphs, plays a significant role in drug classification, toxicity detection, protein analysis etc. However, the limitation of scale in the benchmark datasets makes it easy for graph classification models to fall into over-fitting and undergeneralization. To improve this, we introduce data augmentation on graphs (i.e. graph augmentation) and present four methods:random mapping, vertex-similarity mapping, motif-random mapping and motif-similarity mapping, to generate more weakly labeled data for small-scale benchmark datasets via heuristic transformation of graph structures. Furthermore, we propose a generic model evolution framework, named M-Evolve, which combines graph augmentation, data filtration and model retraining to optimize pre-trained graph classifiers. Experiments on six benchmark datasets demonstrate that the proposed framework helps existing graph classification models alleviate over-fitting and undergeneralization in the training on small-scale benchmark datasets, which successfully yields an average improvement of 3 - 13% accuracy on graph classification tasks.


title: Isomorphisms of $BV(\sigma)$ spaces

abstract: In this paper we investigate the relationship between the properties of a compact set $\sigma \subseteq \mathbb{C}$ and the structure of the space $BV(\sigma)$ of functions of bounded variation (in the sense of Ashton and Doust) defined on $\sigma$. For the subalgebras of absolutely continuous functions on $\sigma$, it is known that for certain classes of compact sets one obtains a Gelfand--Kolmogorov type result: the function spaces $AC(\sigma_1)$ and $AC(\sigma_2)$ are isomorphic if and only if the domain sets $\sigma_1$ and $\sigma_2$ are homeomorphic. Our main theorem is that in this case the isomorphism must extend to an isomorphism of the $BV(\sigma)$ spaces. An application is given to the spectral theory of $AC(\sigma)$ operators.


title: $AC(\sigma)$ spaces for polygonally inscribed curves

abstract: For certain families of compact subsets of the plane, the isomorphism class of the algebra of absolutely continuous functions on a set is completely determined by the homeomorphism class of the set. This is analogous to the Gelfand--Kolmogorov theorem for $C(K)$ spaces. In this paper we define a family of compact sets comprising finite unions of convex curves and show that this family has the `Gelfand--Kolmogorov' property.


title: Graph Neural Networks for Massive MIMO Detection

abstract: In this paper, we innovately use graph neural networks (GNNs) to learn a message-passing solution for the inference task of massive multiple multiple-input multiple-output (MIMO) detection in wireless communication. We adopt a graphical model based on the Markov random field (MRF) where belief propagation (BP) yields poor results when it assumes a uniform prior over the transmitted symbols. Numerical simulations show that, under the uniform prior assumption, our GNN-based MIMO detection solution outperforms the minimum mean-squared error (MMSE) baseline detector, in contrast to BP. Furthermore, experiments demonstrate that the performance of the algorithm slightly improves by incorporating MMSE information into the prior.


title: On mixing enhancement by secondary baroclinic vorticity in shock-bubble interaction

abstract: To investigate the intrinsic mechanism for mixing enhancement by variable density behaviour, a canonical variable density (VD) mixing extracted from a supersonic streamwise vortex protocol, shock bubble interaction (SBI), is numerically studied and compared with a counterpart of passive scalar (PS) mixing. It is meaningful to observe that the maximum concentration decays much faster in VD SBI than in PS SBI regardless of the shock Mach number (Ma = 1.22 ~ 4). The quasi-Lamb-Oseen type velocity distribution in the PS SBI is found by analyzing the azimuthal velocity that stretches the bubble. Meanwhile, for the VD SBI, an additional stretching enhanced by the secondary baroclinic vorticity (SBV) production contributes to the faster-mixing decay. The underlying mechanism of the SBV-enhanced stretching is further revealed through the density and velocity difference between the shocked light bubble and the heavy ambient air. By combining the SBV-accelerated stretching model and the initial shock compression, a novel mixing time estimation for VD SBI is theoretically proposed by solving the advection-diffusion equation under a deformation field of an axisymmetric vortex with the additional SBV induced azimuthal velocity. Based on the mixing time model, a mixing enhancement number defined by the ratio of VD and PS mixing time further reveals the contribution from the variable-density effect, which implies a better control of density distribution for mixing enhancement in a supersonic streamwise vortex.


title: Nonlinear small-gain theorems for input-to-state stability of infinite interconnections

abstract: We consider infinite heterogeneous networks, consisting of input-to-state stable subsystems of possibly infinite dimension. We show that the network is input-to-state stable, provided that the gain operator satisfies a certain small-gain condition. We show that for finite networks of nonlinear systems this condition is equivalent to the so-called strong small-gain condition of the gain operator (and thus our results extend available results for finite networks), and for infinite networks with a linear gain operator they correspond to the condition that the spectral radius of the gain operator is less than one. We provide efficient criteria for input-to-state stability of infinite networks with linear gains, governed by linear and homogeneous gain operators, respectively.


title: Cascade Network with Guided Loss and Hybrid Attention for Two-view Geometry

abstract: In this paper, we are committed to designing a high-performance network for two-view geometry. We first propose a Guided Loss and theoretically establish the direct negative correlation between the loss and Fn-measure by dynamically adjusting the weights of positive and negative classes during training, so that the network is always trained towards the direction of increasing Fn-measure. By this way, the network can maintain the advantage of the cross-entropy loss while maximizing the Fn-measure. We then propose a hybrid attention block to extract feature, which integrates the bayesian attentive context normalization (BACN) and channel-wise attention (CA). BACN can mine the prior information to better exploit global context and CA can capture complex channel context to enhance the channel awareness of the network. Finally, based on our Guided Loss and hybrid attention block, a cascade network is designed to gradually optimize the result for more superior performance. Experiments have shown that our network achieves the state-of-the-art performance on benchmark datasets.


title: A fully noncommutative Painlev\'e II hierarchy: Lax pair and solutions related to Fredholm determinants

abstract: We consider Fredholm determinants of matrix convolution operators associated to matrix versions of the $n - $th Airy functions. Using the theory of integrable operators, we relate them to a fully noncommutative Painlev\'e II hierarchy, defined through a matrix valued version of the Lenard operators. In particular, the Riemann-Hilbert technique used to study these integrable operators allows to find a Lax pair for each member of the hierarchy. Finally, the coefficients of the Lax matrices are explicitely written in terms of these matrix valued Lenard operators and some solution of the hierarchy are written in terms of Fredholm determinants of the square of the matrix Airy convolution operators.


title: Dark Information in Black Hole with $\lambda\varphi$ Fluid

abstract: It has been shown that the nonthermal spectrum of Hawking radiation will lead to information-carrying correlations between emitted particles in the radiation. The mutual information carried by such correlations can not be locally observed and hence is dark. With dark information, the black hole information is conserved. In this paper, we look for the spherically symmetric black hole solution in a $\lambda\varphi$ fluid model and investigate the radiation spectrum and dark information of the black hole. The spacetime structure of this black hole is similar to that of the Schwarzschild one, while its horizon radius is decreased by the $\lambda\varphi$ fluid. By using the statistical mechanical method, the nonthermal radiation spectrum is calculated. This radiation spectrum is very different from the Schwarzschild case at its last stage because of the effect of the $\lambda\varphi$ fluid. The $\lambda\varphi$ fluid reduces the lifetime of the black hole, but increases the dark information of the Hawking radiation.


title: Scoring Interval Forecasts: Equal-Tailed, Shortest, and Modal Interval

abstract: We consider different types of predictive intervals and ask whether they are elicitable, i.e. are unique minimizers of a loss or scoring function in expectation. The equal-tailed interval is elicitable, with a rich class of suitable loss functions, though subject to translation invariance, or positive homogeneity and differentiability, the Winkler interval score becomes a unique choice. The modal interval also is elicitable, with a sole consistent scoring function, up to equivalence. However, the shortest interval fails to be elicitable relative to practically relevant classes of distributions. These results provide guidance in interval forecast evaluation and support recent choices of performance measures in forecast competitions.


title: Constitutive model for shear-thickening suspensions: Predictions for steady shear with superposed transverse oscillations

abstract: We recently developed a tensorial constitutive model for dense, shear-thickening particle suspensions that combines rate-independent microstructural evolution with a stress-dependent jamming threshold. This gives a good qualitative account for reversing flows, although it quantitatively over-estimates structural anisotropy [J.~J.~J.~Gillissen {\em et al.}, Phys. Rev. Lett. {\bf 123} (21), 214504 (2019)]. Here we use the model to predict the unjamming effect of superposed transverse oscillations on a steady shear flow in the thickened regime [N.~Y.~C.~Lin {\em et al.}, Proc.~Nat.~Acad.~Sci.~USA {\bf 113}, 10774 (2016)]. The model successfully reproduces the oscillation-mediated viscosity drop observed experimentally. We compare the time-dependent components of the stress and microstructure tensors to discrete-element simulations. Although the model correctly captures the main qualitative behaviour, it generally over-predicts the microstructural anisotropy in steady shear, and it under-predicts the number of particle contacts in oscillating shear. It also does not fully capture the correct variation in phase angle between the transverse component of the microstructure and the shear rate oscillations, as the amplitude of the latter is increased. These discrepancies suggest avenues for future improvements to the model.


title: Prime knots with arc index 12 up to 16 crossings

abstract: As a continuation of the previous works to tabulate the prime knots up to arc index 11, we provide the list of prime knots with arc index 12 up to 16 crossings and their minimal grid diagrams. There are 19,513 prime knots of arc index 12 up to 16 crossings.


title: Scaling behavior of density gradient accelerated mixing rate in shock bubble interaction

abstract: Variable-density mixing in shock bubble interaction, a canonical flow of Richtermyer-Meshkov instability, is studied by the high-resolution simulation. While the dissipation mainly controls the passive scalar mixing rate, an objective definition of variable-density mixing rate characterizing the macroscopic mixing formation is still lacking, and the fundamental behavior of mixing rate evolution is not yet well understood. Here, we first show that the variable-density mixing of shock bubble interaction is distinctly different from the previous observations in the passive scalar mixing. The widely-accepted hyperbolic conservation of the first moment of concentration in the scalar mixing, i.e., the conservation of the mean concentration, is violated in variable-density flows. We further combine the compositional transport equation and the divergence relation for the miscible flows to provide the evidence that the existence of density gradient accelerated mixing rate, decomposed by the accelerated dissipation term and redistributed diffusion term, contributes to the anomalous decrease or increase of the mean concentration depending on Atwood number. Further analyzing a number of simulations for the cylindrical or spherical bubbles under a broad range of shock Mach numbers, Reynolds numbers, and Peclet numbers, the density gradient accelerated mixing rate exhibits weak dependent on Peclet numbers, and identifies an Atwood number range with high mixing rate, which can be theoretically predicted based on the mode of hyperbolic conservation violation behavior.


title: Two-stage growth mode for lift-off mechanism in oblique shock-wave/jet interaction

abstract: The lift-off characteristics of supersonic streamwise vortex in oblique shock-wave/jet interaction (OS/JI for short), extracted from a wall-mount ramp injector in scramjet, is studied through Large-eddy simulation method. Shocked helium jet is deformed into a pair of streamwise vortex with a co-rotating companion vortex showing the lift-off characteristic immediately after shock. Based on the objective coordinate system in frame of oblique shock structure, it is found that the nature of three-dimensional lift-off structure of a shockinduced streamwise vortex is inherently and precisely controlled by a two-stage growth mode of structure kinetics of a shock bubble interaction (SBI for short). The striking similar structures between OS/JI and SBI support the proposition that the lift-off of streamwise vortex is the result of a underlying two-dimensional vortical motion. By considering the first stage impulsive linear growth rate, an improved vortex propagation model suitable for SBI is proposed and validated. The lift-off phenomena of both numerical OS/JI case in this paper and wall-mounted ramp injector cases in literature are well explained under the two-stage structure kinetics model of SBI. This model further predicts that for higher free stream Mach number (M > 5), increasing ramp compression shows little effect on elevating streamwise vortex while evident lift-off may occur for lower Mach number (M < 3.5), which offers the new way for preliminary design of streamwise vortex-based ramp injector in scramjet.


title: Measurement on the cosmic curvature using the Gaussian process method

abstract: Inflation predicts that the Universe is spatially flat. The Planck 2018 measurements of the cosmic microwave background anisotropy favour a spatially closed universe at more than 2$\sigma$ confidence level. We use model independent methods to study the issue of cosmic curvature. The method reconstructs the Hubble parameter $H(z)$ from cosmic chronometers data with the Gaussian process method. The distance modulus is then calculated with the reconstructed function $H(z)$ and fitted by type Ia supernovae data. Combining the cosmic chronometers and type Ia supernovae data, we obtain $\Omega_{k0}h^2=0.102\pm 0.066$ which is consistent with a spatially flat universe at the 2$\sigma$ confidence level. By adding the redshift space distortions data to the type Ia supernovae data with a proposed novel model independent method, we obtain $\Omega_{k0}h^2=0.117^{+0.058}_{-0.045}$ and no deviation from $\Lambda$CDM model is found.


title: Self-scattering in protoplanetary disks with dust settling

abstract: Scattering of re-emitted flux is considered to be at least partially responsible for the observed polarisation in the (sub-)millimetre wavelength range of several protoplanetary disks. Although the degree of polarisation produced by scattering is highly dependent on the dust model, early studies investigating this mechanism relied on the assumption of single grain sizes and simple density distribution of the dust. However, in the dense inner regions where this mechanism is usually most efficient, the existence of dust grains with sizes ranging from nanometres to millimetres has been confirmed. Additionally, the presence of gas forces larger grains to migrate vertically towards the disk midplane, introducing a dust segregation in the vertical direction. Using polarisation radiative transfer simulations, we analyse the dependence of the resulting scattered light polarisation at 350 $\mu$m, 850 $\mu$m, 1.3 mm, and 2 mm on various parameters describing protoplanetary disks, including the effect of dust grain settling. We find that the different disk parameters change the degree of polarisation mostly by affecting the anisotropy of the radiation field, the optical depth, or both. It is therefore very challenging to deduce certain disk parameter values directly from polarisation measurements alone. However, assuming a high dust albedo, it is possible to trace the transition from optically thick to optically thin disk regions. The degree of polarisation in most of the considered disk configurations is lower than what is found observationally, implying the necessity to revisit models that describe the dust properties and disk structure.


title: Shanks and Anderson-type acceleration techniques for systems of nonlinear equations

abstract: This paper examines a number of extrapolation and acceleration methods, and introduces a few modifications of the standard Shanks transformation that deal with general sequences. One of the goals of the paper is to lay out a general framework that encompasses most of the known acceleration strategies. The paper also considers the Anderson Acceleration method under a new light and exploits a connection with quasi-Newton methods, in order to establish local linear convergence results of a stabilized version of Anderson Acceleration method. The methods are tested on a number of problems, including a few that arise from nonlinear Partial Differential Equations.