title: CoFiNet: Reliable Coarse-to-fine Correspondences for Robust Point Cloud Registration

abstract: We study the problem of extracting correspondences between a pair of point clouds for registration. For correspondence retrieval, existing works benefit from matching sparse keypoints detected from dense points but usually struggle to guarantee their repeatability. To address this issue, we present CoFiNet - Coarse-to-Fine Network which extracts hierarchical correspondences from coarse to fine without keypoint detection. On a coarse scale and guided by a weighting scheme, our model firstly learns to match down-sampled nodes whose vicinity points share more overlap, which significantly shrinks the search space of a consecutive stage. On a finer scale, node proposals are consecutively expanded to patches that consist of groups of points together with associated descriptors. Point correspondences are then refined from the overlap areas of corresponding patches, by a density-adaptive matching module capable to deal with varying point density. Extensive evaluation of CoFiNet on both indoor and outdoor standard benchmarks shows our superiority over existing methods. Especially on 3DLoMatch where point clouds share less overlap, CoFiNet significantly outperforms state-of-the-art approaches by at least 5% on Registration Recall, with at most two-third of their parameters.


title: Cloaking a nanolaser

abstract: Light emitters are bound to strongly interact with light through enhanced absorption and scattering, which imposes limitations on the design and performance of photonic devices such as solar cells, nanoantennas, and (nano) lasers. Overcoming these limitations forces the use of ineffective nonreciprocity approaches or separation of radiation and scattering in the frequency or time domain. A design that combines the properties of an efficient emitter in one state and the property of being invisible in another state is vital for various applications. In this work, we propose a nanolaser design based on a semiconductor nanoparticle with gain coated by a phase transition material (Sb2S3), switchable between lasing and cloaking (nonscattering) states at the same operating frequency without change in pumping. The operation characteristics of the nanolaser are rigorously investigated. The designed nanolaser can operate with optical or electric pumping and possesses attributes of a thresholdless laser due to the high beta-factor and strong Purcell enhancement in the strongly confined Mie resonance mode. We design a reconfigurable metasurface composed of lasing-cloaking metaatoms that can switch from lasing to a nonscattering state in a reversible manner.


title: Adversarial Online Learning with Variable Plays in the Pursuit-Evasion Game: Theoretical Foundations and Application in Connected and Automated Vehicle Cybersecurity

abstract: We extend the adversarial/non-stochastic multi-play multi-armed bandit (MPMAB) to the case where the number of arms to play is variable. The work is motivated by the fact that the resources allocated to scan different critical locations in an interconnected transportation system change dynamically over time and depending on the environment. By modeling the malicious hacker and the intrusion monitoring system as the attacker and the defender, respectively, we formulate the problem for the two players as a sequential pursuit-evasion game. We derive the condition under which a Nash equilibrium of the strategic game exists. For the defender side, we provide an exponential-weighted based algorithm with sublinear pseudo-regret. We further extend our model to heterogeneous rewards for both players, and obtain lower and upper bounds on the average reward for the attacker. We provide numerical experiments to demonstrate the effectiveness of a variable-arm play.


title: Nested Fermi surfaces and correlated electronic phases in hole-doped semiconductor quantum wells

abstract: We demonstrate the existence of novel interaction effects in hole-doped semiconductor quantum wells which are connected to dramatic changes in the Fermi surface geometry occurring upon variation of the doping. We present band structure calculations showing that quantum wells formed in $p$-type cubic semiconductors develop nested Fermi surfaces at a critical hole density set by the width $d$ of the quantum well $k_F\sim \pi/d$. Nesting gives rise to competing superconducting and charge or spin density wave order, which we analyze using the perturbative renormalization group method. The correlated phases may be created or destroyed by tuning the hole density towards or away from the critical density. Our results establish $p$-type semiconductor quantum wells as a platform for novel correlated phases, which may be precisely controlled using electrostatic gating and external magnetic fields.


title: Search for an anomalous excess of charged-current quasi-elastic $\nu_e$ interactions with the MicroBooNE experiment using Deep-Learning-based reconstruction

abstract: We present a measurement of the $\nu_e$-interaction rate in the MicroBooNE detector that addresses the observed MiniBooNE anomalous low-energy excess (LEE). The approach taken isolates neutrino interactions consistent with the kinematics of charged-current quasi-elastic (CCQE) events. The topology of such signal events has a final state with 1 electron, 1 proton, and 0 mesons ($1e1p$). Multiple novel techniques are employed to identify a $1e1p$ final state, including particle identification that use two methods of deep-learning-based image identification, and event isolation using a boosted decision-tree ensemble trained to recognize two-body scattering kinematics. This analysis selects 25 $\nu_e$-candidate events in the reconstructed neutrino energy range of 200--1200\,MeV, while $29.0 \pm 1.9_\text{(sys)} \pm 5.4_\text{(stat)}$ are predicted when using $\nu_\mu$ CCQE interactions as a constraint. We use a simplified model to translate the MiniBooNE LEE observation into a prediction for a $\nu_e$ signal in MicroBooNE. A $\Delta \chi^2$ test statistic, based on the combined Neyman--Pearson $\chi^2$ formalism, is used to define frequentist confidence intervals for the LEE signal strength. Using this technique, in the case of no LEE signal, we expect this analysis to exclude a normalization factor of 0.75 (0.98) times the median MiniBooNE LEE signal strength at 90\% ($2\sigma$) confidence level, while the MicroBooNE data yield an exclusion of 0.25 (0.38) times the median MiniBooNE LEE signal strength at 90\% ($2\sigma$) confidence


title: A Controlled Experiment of Different Code Representations for Learning-Based Bug Repair

abstract: Training a deep learning model on source code has gained significant traction recently. Since such models reason about vectors of numbers, source code needs to be converted to a code representation before vectorization. Numerous approaches have been proposed to represent source code, from sequences of tokens to abstract syntax trees. However, there is no systematic study to understand the effect of code representation on learning performance. Through a controlled experiment, we examine the impact of various code representations on model accuracy and usefulness in deep learning-based program repair. We train 21 different generative models that suggest fixes for name-based bugs, including 14 different homogeneous code representations, four mixed representations for the buggy and fixed code, and three different embeddings. We assess if fix suggestions produced by the model in various code representations are automatically patchable, meaning they can be transformed to a valid code that is ready to be applied to the buggy code to fix it. We also conduct a developer study to qualitatively evaluate the usefulness of inferred fixes in different code representations. Our results highlight the importance of code representation and its impact on learning and usefulness. Our findings indicate that (1) while code abstractions help the learning process, they can adversely impact the usefulness of inferred fixes from a developer's point of view; this emphasizes the need to look at the patches generated from the practitioner's perspective, which is often neglected in the literature, (2) mixed representations can outperform homogeneous code representations, (3) bug type can affect the effectiveness of different code representations; although current techniques use a single code representation for all bug types, there is no single best code representation applicable to all bug types.


title: Multifidelity multilevel Monte Carlo to accelerate approximate Bayesian parameter inference for partially observed stochastic processes

abstract: Models of stochastic processes are widely used in almost all fields of science. Theory validation, parameter estimation, and prediction all require model calibration and statistical inference using data. However, data are almost always incomplete observations of reality. This leads to a great challenge for statistical inference because the likelihood function will be intractable for almost all partially observed stochastic processes. This renders many statistical methods, especially within a Bayesian framework, impossible to implement. Therefore, computationally expensive likelihood-free approaches are applied that replace likelihood evaluations with realisations of the model and observation process. For accurate inference, however, likelihood-free techniques may require millions of expensive stochastic simulations. To address this challenge, we develop a new method based on recent advances in multilevel and multifidelity. Our approach combines the multilevel Monte Carlo telescoping summation, applied to a sequence of approximate Bayesian posterior targets, with a multifidelity rejection sampler to minimise the number of computationally expensive exact simulations required for accurate inference. We present the derivation of our new algorithm for likelihood-free Bayesian inference, discuss practical implementation details, and demonstrate substantial performance improvements. Using examples from systems biology, we demonstrate improvements of more than two orders of magnitude over standard rejection sampling techniques. Our approach is generally applicable to accelerate other sampling schemes, such as sequential Monte Carlo, to enable feasible Bayesian analysis for realistic practical applications in physics, chemistry, biology, epidemiology, ecology and economics.


title: Internal Structure of the 2019 April 2 CME

abstract: We present the first analysis of internal coronal mass ejection (CME) structure observed very close to the Sun by the Wide-field Imager for Solar PRobe (WISPR) instrument on board Parker Solar Probe (PSP). The transient studied here is a CME observed during PSP's second perihelion passage on 2019 April 2, when PSP was only 40 R_sun from the Sun. The CME was also well observed from 1 au by the STEREO-A spacecraft, which tracks the event all the way from the Sun to 1 au. However, PSP/WISPR observes internal structure not apparent in the images from 1 au. In particular, two linear features are observed, one bright and one dark. We model these features as two loops within the CME flux rope channel. The loops can be interpreted as bundles of field lines, with the brightness of the bright loop indicative of lots of mass being loaded into those field lines, and with the dark loop being devoid of such mass loading. It is possible that these loops are actually representative of two independent flux rope structures within the overall CME outline.


title: On the modeling of $t\bar{t}W^\pm$ signatures at the LHC

abstract: We discuss our recent progress on improving the theoretical description of the $pp \to t\bar{t}W^\pm$ process with particular focus on fiducial signatures that include top-quark decays. We employ a wide range of techniques from parton shower matching, narrow width approximation to full off-shell computations to assess the importance of the different modeling approaches.


title: Vertical Injection and Wideband Grating Coupler Based on Asymmetric Grating Trenches

abstract: A Silicon-on-insulator (SOI) perfectly vertical fibre-to-chip grating coupler is proposed and designed based on engineered subwavelength structures. The high directionality of the coupler is achieved by implementing step gratings to realize asymmetric diffraction and by applying effective index variation with auxiliary ultra-subwavelength gratings. The proposed structure is numerically analysed by using two-dimensional Finite Difference Time Domain (2D FDTD) method and achieves 76% (-1.19 dB) coupling efficiency and 39 nm 1-dB bandwidth.


title: Search for $B^0$ meson decays into $\Lambda$ and missing energy with a hadronic tagging method at Belle

abstract: We present a search for the decays of $B^0$ mesons into a final state containing a $\Lambda$ baryon and missing energy. These results are obtained from a $711\text{fb}^{-1}$ data sample that contains $772 \times 10^6$ $B{\kern 0.18em}\overline{\kern -0.18em B}$ pairs and was collected near the $\Upsilon(4S)$ resonance with the Belle detector at the KEKB asymmetric-energy $e^+e^-$ collider. We use events in which one $B$ meson is fully reconstructed in a hadronic decay mode and require the remainder of the event to consist of only a single $\Lambda$. No evidence for these decays is found and we set $90\%$ confidence level upper limits on the branching fractions in the range $2.1$-$3.8\times 10^{-5}$. This measurement provides the world's most restrictive limits, with implications for baryogenesis and dark matter production.


title: Using the Climate App to learn about Planetary Habitability and Climate Change

abstract: Simple climate models have been around for more than a century but have recently come back into fashion: they are useful for explaining global warming and the habitability of extrasolar planets. The Climate App (https://www.climateapp.ca) is an interactive web-based application that describes the radiative transfer governing planetary climate. The App is currently available in French and English and is suitable for teaching high-school through college students, or public outreach. The beginner version can be used to explore the greenhouse effect and planetary albedo, sufficient for explaining anthropogenic climate change, the Faint Young Sun Paradox, the habitability of TRAPPIST planets and other simple scenarios. There is also an advanced option with more atmospheric layers and incorporating the absorption and scattering of shortwave radiation for students and educators wishing a deeper dive into atmospheric radiative transfer. A number of pedagogical activities are being beta tested and rolled out.


title: Shear-Driven Flow in an Elliptical Enclosure Generated by an Inner Rotating Circular Cylinder

abstract: Shear-driven flow between a rotating cylinder and a stationary elliptical enclosure is studied in this paper. Two-dimensional time-dependent Navier Stokes equations are solved using a meshless method where interpolations are done with Polyharmonic Spline Radial Basis Functions. The fluid flow is analyzed for various aspect ratios of the ellipse and eccentric placements of the inner cylinder. Contour plots of vorticity with streamlines, plots of non-dimensional torque, and the angle of eye of the primary vortex are presented in the paper for Reynolds numbers between 200 and 2000. Formation of Moffatt like vortices in the wide-gap region of the model is observed and some benchmark data are provided for various cases that are simulated.


title: The density distributions of cosmic structures: impact of the local environment on weak-lensing convergence

abstract: Whilst the underlying assumption of the Friedman-Lema\^itre-Robertson-Walker (FLRW) cosmological model is that matter is homogeneously distributed throughout the universe, gravitational influences over the life of the universe have resulted in mass clustered on a range of scales. Hence we expect that, in our inhomogeneous universe, the view of an observer will be influenced by the location and local environment. Here we analyse the one-point probability distribution functions and angular power spectra of weak-lensing (WL) convergence and magnification numerically to investigate the influence of our local environment on WL statistics in relativistic $N$-body simulations. To achieve this, we numerically solve the null geodesic equations which describe the propagation of light bundles backwards in time from today, and develop a ray-tracing algorithm, and from these calculate various WL properties. Our findings demonstrate how cosmological observations of large-scale structure through WL can be impacted by the locality of the observer. We also calculate the constraints on the cosmological parameters as a function of redshift from the theoretical and numerical study of the angular power spectrum of WL convergence. This study concludes the minimal redshift for the constraint on the parameter $\Omega_m$ ($H_0$) is $z \sim 0.2$ $(z \sim 0.6 )$ beyond which the local environment's effect is negligible and the data from WL surveys are more meaningful above that redshift. The outcomes of this study will have direct consequences for future surveys, where percent-level-precision is necessary.


title: Teardown and feasibility study of IronKey -- the most secure USB Flash drive

abstract: There are many solutions for protecting user data on USB Flash drives. However, the family of IronKey devices was designed with the highest security expectations. They are definitely standing above others by being certified to FIPS 140-2 Level 3 and also claimed as certified by NATO for Top-Secret use. Many encrypted USB drives had been evaluated and found insecure, however, no public research on IronKey devices was made. This feasibility study fills the gap by looking inside the IronKey family of devices. As a result the users of the IronKey devices could be assured about the real level of the security protection they get. Several generations of devices from IronKey family and competitors are teared down, their hardware solutions discussed and evaluated for possible attacks. Some potential flaws are exposed and those findings are likely to stimulate further research into specific solutions aimed to protect user data.


title: Connect-the-Dots: Bridging Semantics between Words and Definitions via Aligning Word Sense Inventories

abstract: Word Sense Disambiguation (WSD) aims to automatically identify the exact meaning of one word according to its context. Existing supervised models struggle to make correct predictions on rare word senses due to limited training data and can only select the best definition sentence from one predefined word sense inventory (e.g., WordNet). To address the data sparsity problem and generalize the model to be independent of one predefined inventory, we propose a gloss alignment algorithm that can align definition sentences (glosses) with the same meaning from different sense inventories to collect rich lexical knowledge. We then train a model to identify semantic equivalence between a target word in context and one of its glosses using these aligned inventories, which exhibits strong transfer capability to many WSD tasks. Experiments on benchmark datasets show that the proposed method improves predictions on both frequent and rare word senses, outperforming prior work by 1.2% on the All-Words WSD Task and 4.3% on the Low-Shot WSD Task. Evaluation on WiC Task also indicates that our method can better capture word meanings in context.


title: BioGrad: Biologically Plausible Gradient-Based Learning for Spiking Neural Networks

abstract: Spiking neural networks (SNN) are delivering energy-efficient, massively parallel, and low-latency solutions to AI problems, facilitated by the emerging neuromorphic chips. To harness these computational benefits, SNN need to be trained by learning algorithms that adhere to brain-inspired neuromorphic principles, namely event-based, local, and online computations. Yet, the state-of-the-art SNN training algorithms are based on backprop that does not follow the above principles. Due to its limited biological plausibility, the application of backprop to SNN requires non-local feedback pathways for transmitting continuous-valued errors, and relies on gradients from future timesteps. The introduction of biologically plausible modifications to backprop has helped overcome several of its limitations, but limits the degree to which backprop is approximated, which hinders its performance. We propose a biologically plausible gradient-based learning algorithm for SNN that is functionally equivalent to backprop, while adhering to all three neuromorphic principles. We introduced multi-compartment spiking neurons with local eligibility traces to compute the gradients required for learning, and a periodic "sleep" phase to further improve the approximation to backprop during which a local Hebbian rule aligns the feedback and feedforward weights. Our method achieved the same level of performance as backprop with multi-layer fully connected SNN on MNIST (98.13%) and the event-based N-MNIST (97.59%) datasets. We deployed our learning algorithm on Intel's Loihi to train a 1-hidden-layer network for MNIST, and obtained 93.32% test accuracy while consuming 400 times less energy per training sample than BioGrad on GPU. Our work shows that optimal learning is feasible in neuromorphic computing, and further pursuing its biological plausibility can better capture the benefits of this emerging computing paradigm.


title: The Physical Parameters of Four WC-type Wolf-Rayet Stars in the Large Magellanic Cloud: Evidence of Evolution

abstract: We present a spectral analysis of four LMC WC-type Wolf-Rayet (WR) stars (BAT99-8, BAT99-9, BAT99-11, and BAT99-52) to shed light on two evolutionary questions surrounding massive stars. The first is: are WO-type WR stars more oxygen enriched than the WC-type stars, indicating further chemical evolution, or are the strong high-excitation oxygen lines in the WO-type stars an indication of higher temperatures. This study will act as a baseline for answering the question of where WO-type stars fall in WR evolution. Each star's spectrum, extending from 1100~\AA\ to 25000~\AA, was modeled using \cmfgen\ to determine the star's physical properties such as luminosity, mass-loss rate, and chemical abundances. The oxygen abundance is a key evolutionary diagnostic, and with higher resolution data and an improved stellar atmosphere code, we found the oxygen abundance to be up to a factor of 5 lower than previous studies. The second evolutionary question revolves around the formation of WR stars: Do they evolve by themselves or is a close companion star necessary for their formation? Using our derived physical parameters, we compared our results to the Geneva single-star evolutionary models and the BPASS binary evolutionary models. We found that both the Geneva solar metallicity models and BPASS LMC metallicity models are in agreement with the four WC-type stars, while the Geneva LMC metallicity models are not. Therefore, these four WC4 stars could have been formed either via binary or single-star evolution.


title: Learning-Augmented $k$-means Clustering

abstract: $k$-means clustering is a well-studied problem due to its wide applicability. Unfortunately, there exist strong theoretical limits on the performance of any algorithm for the $k$-means problem on worst-case inputs. To overcome this barrier, we consider a scenario where "advice" is provided to help perform clustering. Specifically, we consider the $k$-means problem augmented with a predictor that, given any point, returns its cluster label in an approximately optimal clustering up to some, possibly adversarial, error. We present an algorithm whose performance improves along with the accuracy of the predictor, even though na\"{i}vely following the accurate predictor can still lead to a high clustering cost. Thus if the predictor is sufficiently accurate, we can retrieve a close to optimal clustering with nearly optimal runtime, breaking known computational barriers for algorithms that do not have access to such advice. We evaluate our algorithms on real datasets and show significant improvements in the quality of clustering.


title: Third Order Modal Exceptional Degeneracy in Waveguides with Glide-Time Symmetry

abstract: The dispersion of a three-way waveguide is engineered to exhibit exceptional modal characteristics. Two coupled waveguides with Parity-Time (PT) symmetry have been previously demonstrated to exhibit second order exceptional points of degeneracy (EPDs). In this work, we introduce and investigate a particular class of EPDs, applicable from radio frequency to optical wavelengths, whereby three coupled waveguides satisfy Glide-Time (GT) symmetry to exhibit a third order modal degeneracy with a real-valued wavenumber. GT symmetry involves glide symmetry of lossless/gainless components of the waveguide in addition to changing the sign of passive/active elements while applying a glide symmetry operation. This GT-symmetry condition allows three Floquet-Bloch eigenmodes of the structure to coalesce to a real-valued wavenumber at a single frequency, in addition of having one branch of the dispersion diagram with a purely real wavenumber. The proposed scheme may have applications including but not limited to distributed amplifiers, radiating arrays, and sensors, from radio frequency to optics.


title: Towards Robust Bisimulation Metric Learning

abstract: Learned representations in deep reinforcement learning (DRL) have to extract task-relevant information from complex observations, balancing between robustness to distraction and informativeness to the policy. Such stable and rich representations, often learned via modern function approximation techniques, can enable practical application of the policy improvement theorem, even in high-dimensional continuous state-action spaces. Bisimulation metrics offer one solution to this representation learning problem, by collapsing functionally similar states together in representation space, which promotes invariance to noise and distractors. In this work, we generalize value function approximation bounds for on-policy bisimulation metrics to non-optimal policies and approximate environment dynamics. Our theoretical results help us identify embedding pathologies that may occur in practical use. In particular, we find that these issues stem from an underconstrained dynamics model and an unstable dependence of the embedding norm on the reward signal in environments with sparse rewards. Further, we propose a set of practical remedies: (i) a norm constraint on the representation space, and (ii) an extension of prior approaches with intrinsic rewards and latent space regularization. Finally, we provide evidence that the resulting method is not only more robust to sparse reward functions, but also able to solve challenging continuous control tasks with observational distractions, where prior methods fail.


title: Automated Evaluation of Web Site Accessibility Using A Dynamic Accessibility Measurement Crawler

abstract: Achieving accessibility compliance is extremely important for many government agencies and businesses who wish to improve services for their consumers. With the growing reliance on dynamic web applications many organizations are finding it difficult to implement accessibility standards, often due to the inability of current automated testing tools to test the stateful environments created by dynamic web applications. In this paper, we present mathematical foundations and theory for the Demodocus framework and prototype, and outline its approach to using web science, web crawling,and accessibility testing to automatically navigate and test interactive content for accessibility. Our approach simulates the page interactions of users with and without disabilities, and compares graphs of reachable states from these simulations to determine both the accessibility and the difficulty of content access for these different users.


title: Provable Lifelong Learning of Representations

abstract: In lifelong learning, tasks (or classes) to be learned arrive sequentially over time in arbitrary order. During training, knowledge from previous tasks can be captured and transferred to subsequent ones to improve sample efficiency. We consider the setting where all target tasks can be represented in the span of a small number of unknown linear or nonlinear features of the input data. We propose a lifelong learning algorithm that maintains and refines the internal feature representation. We prove that for any desired accuracy on all tasks, the dimension of the representation remains close to that of the underlying representation. The resulting sample complexity improves significantly on existing bounds. In the setting of linear features, our algorithm is provably efficient and the sample complexity for input dimension $d$, $m$ tasks with $k$ features up to error $\epsilon$ is $\tilde{O}(dk^{1.5}/\epsilon+km/\epsilon)$. We also prove a matching lower bound for any lifelong learning algorithm that uses a single task learner as a black box. We complement our analysis with an empirical study, including a heuristic lifelong learning algorithm for deep neural networks. Our method performs favorably on challenging realistic image datasets compared to state-of-the-art continual learning methods.


title: Tight Concentrations and Confidence Sequences from the Regret of Universal Portfolio

abstract: A classic problem in statistics is the estimation of the expectation of random variables from samples. This gives rise to the tightly connected problems of deriving concentration inequalities and confidence sequences, that is confidence intervals that hold uniformly over time. Previous work has shown how to easily convert the regret guarantee of an online betting algorithm into a time-uniform concentration inequality. In this paper, we show that we can go even further: We show that the regret of universal portfolio algorithms give rise to new implicit time-uniform concentrations and state-of-the-art empirically calculated confidence sequences. In particular, our numerically obtained confidence sequences can never be vacuous, even with a single sample, and satisfy the law of iterated logarithm.

categories: stat.ML cs.IT cs.LG math.IT math.ST stat.ME stat.TH


title: An Algebraic Approach to Degenerate Bernoulli Numbers

abstract: In this work we study the properties of a new algebraic variant of the degenerate Bernoulli polynomial $\tilde{\beta}_{k}(m,x)$ and study the corresponding degenerate Bernoulli number $\tilde{\beta}_{k}(m,1)=m^{k}\beta_{k}(1/m)$, where $\beta_{k}(\lambda), \lambda\neq 0$ is the standard degenerate Bernoulli number. Our approach relies on a new algebraic framework for generating functions and the action of a symbolic evaluation function on powers of polynomials. We show that $\tilde{\beta}_{k}(m,x)$ displays surprising links with other mathematical objects (such as Circulant matrices and Galois fields) and enjoys many interesting algebraic, symmetric and dynamical properties that could be deployed to perform efficient computations.


title: Determination of hybrid charmonium meson masses

abstract: We report initial results from our study of the masses and decay constants of the lightest multiplet of charmonium-like hybrid mesons. We obtain precise measurements of the $1^{-+}$ state through the use of a variational basis and a large number of configurations at three lattice spacings. We use staggered fermion operators using configurations generated with the HISQ action with 2+1+1 dynamical flavours. The mixing of the vector hybrid with the $J/\psi$ is examined and a preliminary bound on the vector hybrid decay constant is presented.


title: Breather and interacting soliton and periodic waves for modified KdV equation

abstract: We present the discovery of a class of exact spatially localized as well as periodic wave solutions within the framework of the modified Korteweg-de Vries equation. This class comprises breather and interacting soliton solutions as well as interacting periodic wave solutions. The functional forms of these solutions include a joint parameter which can take both positive and negative values of unity. It is found that the existence of those closed form solutions depend strongly on whether the cubic nonlinearity parameter should be considered positive or negative. The derived wave structures show interesting properties that may find practical applications.


title: Gradient Ricci solitons carrying a closed conformal vector field

abstract: We show that a complete gradient Ricci soliton $(M^n,\,g)$ with constant scalar curvature and a non-parallel closed conformal vector field is isometric to either the Euclidean space, or an Euclidean sphere, or negatively Einstein warped product of the real line with a complete non-positively Einstein manifold. Moreover, we show that a K\"ahler gradient Ricci soliton of real dimension $\ge 4$, with a non-parallel closed conformal real vector field is Ricci-flat (Calabi-Yau) and is flat in dimension 4. Finally, we show that a Ricci soliton whose associated 1-form is harmonic, has constant scalar curvature.


title: Evidence for equilibrium excitons and exciton condensation in monolayer WTe2

abstract: A single monolayer of the layered semimetal WTe2 behaves as a two-dimensional topological insulator, with helical conducting edge modes surrounding a bulk state that becomes insulating at low temperatures. Here we present evidence that the bulk state has a very unusual nature, containing electrons and holes bound by Coulomb attraction (excitons) that spontaneously form in thermal equilibrium. On cooling from room temperature to 100 K the conductivity develops a V-shaped dependence on electrostatic doping, while the chemical potential develops a ~43 meV step at the neutral point. These features are much sharper than is possible in an independent-electron picture, but they can be largely accounted for by positing that some of the electrons and holes are paired in equilibrium. Our calculations from first principles show that the exciton binding energy is larger than 100 meV and the radius as small as 4 nm, explaining their formation at high temperature and doping levels. Below 100 K more strongly insulating behavior is seen, suggesting that a charge-ordered state forms. The observed absence of charge density waves in this state appears surprising within an excitonic insulator picture, but we show that it can be explained by the symmetries of the exciton wave function. Monolayer WTe2 therefore presents an exceptional combination of topological properties and strong correlations over a wide temperature range.


title: Negative string tension of higher-charge Schwinger model via digital quantum simulation

abstract: We study some properties of generalized global symmetry for the charge-$q$ Schwinger model in the Hamiltonian formalism, which is the $(1+1)$-dimensional quantum electrodynamics with a charge-$q$ Dirac fermion. This model has the $\mathbb{Z}_q$ $1$-form symmetry, which is a remnant of the electric $U(1)$ $1$-form symmetry in the pure Maxwell theory. It is known that if we put the theory on closed space, then the Hilbert space is decomposed into $q$ distinct sectors, called universes and some states with higher energy density do not decay to the ground state due to the selection rule of the $1$-form symmetry. Even with open boundaries, we can observe the stability of such states by seeing a negative string tension behavior, meaning that opposite charges repel with each other. We develop a method based on the adiabatic state preparation to see this feature with digital quantum simulation and confirm it using a classical simulator of quantum devices. Especially, we measure local energy density and see how it jumps between inside and outside of insertion of the probe charges. We explicitly see that the energy density inside is lower than the outside one. This is a clear signature of the negative string tension.


title: Gemini/GMOS Transmission Spectroscopy of the Grazing Planet Candidate WD 1856+534 b

abstract: WD 1856+534 b is a Jupiter-sized, cool giant planet candidate transiting the white dwarf WD 1856+534. Here, we report an optical transmission spectrum of WD 1856+534 b obtained from ten transits using the Gemini Multi-Object Spectrograph. This system is challenging to observe due to the faintness of the host star and the short transit duration. Nevertheless, our phase-folded white light curve reached a precision of 0.12 %. WD 1856+534 b provides a unique transit configuration compared to other known exoplanets: the planet is $8\times$ larger than its star and occults over half of the stellar disc during mid-transit. Consequently, many standard modeling assumptions do not hold. We introduce the concept of a `limb darkening corrected, time-averaged transmission spectrum' and propose that this is more suitable than $(R_{\mathrm{p}, \lambda} / R_{\mathrm{s}})^2$ for comparisons to atmospheric models for planets with grazing transits. We also present a modified radiative transfer prescription. Though the transmission spectrum shows no prominent absorption features, it is sufficiently precise to constrain the mass of WD 1856+534 b to be > 0.84 M$_\mathrm{J}$ (to $2 \, \sigma$ confidence), assuming a clear atmosphere and a Jovian composition. High-altitude cloud decks can allow lower masses. WD 1856+534 b could have formed either as a result of common envelope evolution or migration under the Kozai-Lidov mechanism. Further studies of WD 1856+534 b, alongside new dedicated searches for substellar objects around white dwarfs, will shed further light on the mysteries of post-main sequence planetary systems.


title: Electromagnetic form factors of baryons in nuclear medium

abstract: The electromagnetic structure of the baryons is modified in the nuclear medium. The modifications can be inferred from the comparison between the electromagnetic form factors in medium with the respective form factor in vacuum. Of particular interest is the ratio between the electric and magnetic form factors in medium ($G_E^*/G_M^*$) and vacuum ($G_E/G_M$) of the octet baryon. The deviation of the double ratios ($G_E^*/G_M^*)/(G_E/G_M)$ from unity measures the impact of the medium modification of the electromagnetic structure in a nuclear medium. Measurements of the double ratios $(G_E^*/G_M^*)/(G_E/G_M)$ for different nuclear densities may become available in a near future using the polarization-transfer method developed at Jefferson Lab. We present estimates of the double ratios of octet baryons based on a covariant constituent quark model, which takes into account pion cloud excitations of the baryon cores,for different nuclear densities. Our results manifest different features, namely, enhancement or quenching depending on the baryon flavor content.


title: Quality, Speed, and Scale: three key attributes to measure the performance of near-term quantum computers

abstract: Defining the right metrics to properly represent the performance of a quantum computer is critical to both users and developers of a computing system. In this white paper, we identify three key attributes for quantum computing performance: quality, speed, and scale. Quality and scale are measured by quantum volume and number of qubits, respectively. We propose a speed benchmark, using an update to the quantum volume experiments that allows the measurement of Circuit Layer Operations Per Second (CLOPS) and identify how both classical and quantum components play a role in improving performance. We prescribe a procedure for measuring CLOPS and use it to characterize the performance of some IBM Quantum systems.


title: Eigencurve: Optimal Learning Rate Schedule for SGD on Quadratic Objectives with Skewed Hessian Spectrums

abstract: Learning rate schedulers have been widely adopted in training deep neural networks. Despite their practical importance, there is a discrepancy between its practice and its theoretical analysis. For instance, it is not known what schedules of SGD achieve best convergence, even for simple problems such as optimizing quadratic objectives. In this paper, we propose Eigencurve, the first family of learning rate schedules that can achieve minimax optimal convergence rates (up to a constant) for SGD on quadratic objectives when the eigenvalue distribution of the underlying Hessian matrix is skewed. The condition is quite common in practice. Experimental results show that Eigencurve can significantly outperform step decay in image classification tasks on CIFAR-10, especially when the number of epochs is small. Moreover, the theory inspires two simple learning rate schedulers for practical applications that can approximate eigencurve. For some problems, the optimal shape of the proposed schedulers resembles that of cosine decay, which sheds light to the success of cosine decay for such situations. For other situations, the proposed schedulers are superior to cosine decay.


title: Mining frequency-based sequential trajectory co-clusters

abstract: Co-clustering is a specific type of clustering that addresses the problem of finding groups of objects without necessarily considering all attributes. This technique has shown to have more consistent results in high-dimensional sparse data than traditional clustering. In trajectory co-clustering, the methods found in the literature have two main limitations: first, the space and time dimensions have to be constrained by user-defined thresholds; second, elements (trajectory points) are clustered ignoring the trajectory sequence, assuming that the points are independent among them. To address the limitations above, we propose a new trajectory co-clustering method for mining semantic trajectory co-clusters. It simultaneously clusters the trajectories and their elements taking into account the order in which they appear. This new method uses the element frequency to identify candidate co-clusters. Besides, it uses an objective cost function that automatically drives the co-clustering process, avoiding the need for constraining dimensions. We evaluate the proposed approach using real-world a publicly available dataset. The experimental results show that our proposal finds frequent and meaningful contiguous sequences revealing mobility patterns, thereby the most relevant elements.


abstract: An invertible matrix is called a Perron similarity if one of its columns and the corresponding row of its inverse are both nonnegative or both nonpositive. Such matrices are of relevance and import in the study of the nonnegative inverse eigenvalue problem. In this work, Kronecker products of Perron similarities are examined and used to to construct ideal Perron similarities all of whose rows are extremal.


title: A Bayesian Receiver with Improved Complexity-Reliability Trade-off in Massive MIMO Systems

abstract: The stringent requirements on reliability and processing delay in the fifth-generation ($5$G) cellular networks introduce considerable challenges in the design of massive multiple-input-multiple-output (M-MIMO) receivers. The two main components of an M-MIMO receiver are a detector and a decoder. To improve the trade-off between reliability and complexity, a Bayesian concept has been considered as a promising approach that enhances classical detectors, e.g. minimum-mean-square-error detector. This work proposes an iterative M-MIMO detector based on a Bayesian framework, a parallel interference cancellation scheme, and a decision statistics combining concept. We then develop a high performance M-MIMO receiver, integrating the proposed detector with a low complexity sequential decoding for polar codes. Simulation results of the proposed detector show a significant performance gain compared to other low complexity detectors. Furthermore, the proposed M-MIMO receiver with sequential decoding ensures one order magnitude lower complexity compared to a receiver with stack successive cancellation decoding for polar codes from the 5G New Radio standard.


title: Regularisation in Nonperturbative Extensions of Effective Field Theory

abstract: The process of renormalisation in nonperturbative Hamiltonian Effective Field Theory (HEFT) is examined in the $\Delta$-resonance scattering channel. As an extension of effective field theory incorporating the L\"uscher formalism, HEFT provides a bridge between the infinite-volume scattering data of experiment and the finite-volume spectrum of energy eigenstates in lattice QCD. HEFT also provides phenomenological insight into the basis-state composition of the finite-volume eigenstates via the state eigenvectors. The Hamiltonian matrix is made finite through the introduction of finite-range regularisation. The extent to which the established features of this regularisation scheme survive in HEFT is examined. In a single-channel $\pi N$ analysis, fits to experimental phase shifts withstand large variations in the regularisation parameter, $\Lambda$, providing an opportunity to explore the sensitivity of the finite-volume spectrum and state composition on the regulator. While the L\"uscher formalism ensures the eigenvalues are insensitive to $\Lambda$ variation in the single-channel case, the eigenstate composition varies with $\Lambda$; the admission of short distance interactions diminishes single-particle contributions to the states. In the two-channel $\pi N$, $\pi \Delta$ analysis, $\Lambda$ is restricted to a small range by the experimental data. Here the inelasticity is particularly sensitive to variations in $\Lambda$ and its associated parameter set. This sensitivity is also manifest in the finite-volume spectrum for states near the opening of the $\pi \Delta$ scattering channel. Finally, HEFT has the unique ability to describe the quark-mass dependence of the finite-volume eigenstates. The robust nature of this capability is presented and used to confront current state-of-the-art lattice QCD calculations.


title: A C-Band Fully Polarimetric Automotive Synthetic Aperture Radar

abstract: Due to the rapid increase in 76 GHz automotive spectrum use in recent years, wireless interference is becoming a legitimate area of concern. However, the recent rise in interest of automated vehicles (AVs) has also spurred new growth and adoption of low frequency vehicle-to-everything (V2X) communications in and around the 5.8 GHz unlicensed bands, opening the possibility for communications spectrum reuse in the form of joint radar-communications (JRC). In this work, we present a low frequency 5.9 GHz side-looking polarimetric synthetic aperture radar (SAR) for automotive use, utilizing a ranging waveform in a common low frequency V2X communications band. A synthetic aperture technique is employed to address the angular resolution concerns commonly associated with radars at lower frequencies. Three side-looking fully polarimetric SAR images in various urban scenes are presented and discussed to highlight the unique opportunities for landmark inference afforded through measurement of co- and cross-polarized scattering.


title: Metrics of research impact in astronomy: Predicting later impact from metrics measured 10-15 years after the PhD

abstract: This paper calibrates how metrics derivable from the SAO/NASA Astrophysics Data System can be used to estimate the future impact of astronomy research careers and thereby to inform decisions on resource allocation such as job hires and tenure decisions. Three metrics are used, citations of refereed papers, citations of all publications normalized by the numbers of co-authors, and citations of all first-author papers. Each is individually calibrated as an impact predictor in the book Kormendy (2020), "Metrics of Research Impact in Astronomy" (Astron Soc Pac Conference Series Monograph 8, San Francisco). How this is done is reviewed in the first half of this paper. Then, I show that averaging results from three metrics produces more accurate predictions. Average prediction machines are constructed for different cohorts of 1990-2007 PhDs and used to postdict 2017 impact from metrics measured 10, 12, and 15 years after the PhD. The time span over which prediction is made ranges from 0 years for 2007 PhDs to 17 years for 1990 PhDs using metrics measured 10 years after the PhD. Calibration is based on perceived 2017 impact as voted by 22 experienced astronomers for 510 faculty members at 17 highly-ranked university astronomy departments world-wide. Prediction machinery reproduces voted impact estimates with an RMS uncertainty of 1/8 of the dynamic range for people in the study sample. The aim of this work is to lend some of the rigor that is normally used in scientific research to the difficult and subjective job of judging people's careers.

categories: astro-ph.IM astro-ph.EP astro-ph.SR cs.DL


title: Data-driven decomposition of brain dynamics with principal component analysis in different types of head impacts

abstract: Strain and strain rate are effective traumatic brain injury predictors. Kinematics-based models estimating these metrics suffer from significant different distributions of both kinematics and the injury metrics across head impact types. To address this, previous studies focus on the kinematics but not the injury metrics. We have previously shown the kinematic features vary largely across head impact types, resulting in different patterns of brain deformation. This study analyzes the spatial distribution of brain deformation and applies principal component analysis (PCA) to extract the representative patterns of injury metrics (maximum principal strain (MPS), MPS rate (MPSR) and MPSXMPSR) in four impact types (simulation, football, mixed martial arts and car crashes). We apply PCA to decompose the patterns of the injury metrics for all impacts in each impact type, and investigate the distributions among brain regions using the first principal component (PC1). Furthermore, we developed a deep learning head model (DLHM) to predict PC1 and then inverse-transform to predict for all brain elements. PC1 explained >80% variance on the datasets. Based on PC1 coefficients, the corpus callosum and midbrain exhibit high variance on all datasets. We found MPSXMPSR the most sensitive metric on which the top 5% of severe impacts further deviates from the mean and there is a higher variance among the severe impacts. Finally, the DLHM reached mean absolute errors of <0.018 for MPS, <3.7 (1/s) for MPSR and <1.1 (1/s) for MPSXMPSR, much smaller than the injury thresholds. The brain injury metric in a dataset can be decomposed into mean components and PC1 with high explained variance. The brain dynamics decomposition enables better interpretation of the patterns in brain injury metrics and the sensitivity of brain injury metrics across impact types. The decomposition also reduces the dimensionality of DLHM.


title: Forecasting with a Panel Tobit Model

abstract: We use a dynamic panel Tobit model with heteroskedasticity to generate forecasts for a large cross-section of short time series of censored observations. Our fully Bayesian approach allows us to flexibly estimate the cross-sectional distribution of heterogeneous coefficients and then implicitly use this distribution as prior to construct Bayes forecasts for the individual time series. In addition to density forecasts, we construct set forecasts that explicitly target the average coverage probability for the cross-section. We present a novel application in which we forecast bank-level loan charge-off rates for small banks.


title: Object-Aware Regularization for Addressing Causal Confusion in Imitation Learning

abstract: Behavioral cloning has proven to be effective for learning sequential decision-making policies from expert demonstrations. However, behavioral cloning often suffers from the causal confusion problem where a policy relies on the noticeable effect of expert actions due to the strong correlation but not the cause we desire. This paper presents Object-aware REgularizatiOn (OREO), a simple technique that regularizes an imitation policy in an object-aware manner. Our main idea is to encourage a policy to uniformly attend to all semantic objects, in order to prevent the policy from exploiting nuisance variables strongly correlated with expert actions. To this end, we introduce a two-stage approach: (a) we extract semantic objects from images by utilizing discrete codes from a vector-quantized variational autoencoder, and (b) we randomly drop the units that share the same discrete code together, i.e., masking out semantic objects. Our experiments demonstrate that OREO significantly improves the performance of behavioral cloning, outperforming various other regularization and causality-based methods on a variety of Atari environments and a self-driving CARLA environment. We also show that our method even outperforms inverse reinforcement learning methods trained with a considerable amount of environment interaction.


title: Vertex distortion detects the unknot

abstract: The first two authors introduced vertex distortion and showed that the vertex distortion of the unknot is trivial. It was conjectured that the vertex distortion of a knot is trivial if and only if the knot is trivial. We will use Denne-Sullivan's bound on Gromov distortion to bound the vertex distortion of nontrivial lattice knots. We will then conclude that trivial vertex distortion implies the unknot, proving the conjecture. Additionally, the first conjecture in vertex distortion's debut article is proven and a vertex distortion calculator is given.


title: ScaleCert: Scalable Certified Defense against Adversarial Patches with Sparse Superficial Layers

abstract: Adversarial patch attacks that craft the pixels in a confined region of the input images show their powerful attack effectiveness in physical environments even with noises or deformations. Existing certified defenses towards adversarial patch attacks work well on small images like MNIST and CIFAR-10 datasets, but achieve very poor certified accuracy on higher-resolution images like ImageNet. It is urgent to design both robust and effective defenses against such a practical and harmful attack in industry-level larger images. In this work, we propose the certified defense methodology that achieves high provable robustness for high-resolution images and largely improves the practicality for real adoption of the certified defense. The basic insight of our work is that the adversarial patch intends to leverage localized superficial important neurons (SIN) to manipulate the prediction results. Hence, we leverage the SIN-based DNN compression techniques to significantly improve the certified accuracy, by reducing the adversarial region searching overhead and filtering the prediction noises. Our experimental results show that the certified accuracy is increased from 36.3% (the state-of-the-art certified detection) to 60.4% on the ImageNet dataset, largely pushing the certified defenses for practical use.


title: On Computing the Hyperparameter of Extreme Learning Machines: Algorithm and Application to Computational PDEs, and Comparison with Classical and High-Order Finite Elements

abstract: We consider the use of extreme learning machines (ELM) for computational partial differential equations (PDE). In ELM the hidden-layer coefficients in the neural network are assigned to random values generated on $[-R_m,R_m]$ and fixed, where $R_m$ is a user-provided constant, and the output-layer coefficients are trained by a linear or nonlinear least squares computation. We present a method for computing the optimal value of $R_m$ based on the differential evolution algorithm. The presented method enables us to illuminate the characteristics of the optimal $R_m$ for two types of ELM configurations: (i) Single-Rm-ELM, in which a single $R_m$ is used for generating the random coefficients in all the hidden layers, and (ii) Multi-Rm-ELM, in which multiple $R_m$ constants are involved with each used for generating the random coefficients of a different hidden layer. We adopt the optimal $R_m$ from this method and also incorporate other improvements into the ELM implementation. In particular, here we compute all the differential operators involving the output fields of the last hidden layer by a forward-mode auto-differentiation, as opposed to the reverse-mode auto-differentiation in a previous work. These improvements significantly reduce the network training time and enhance the ELM performance. We systematically compare the computational performance of the current improved ELM with that of the finite element method (FEM), both the classical second-order FEM and the high-order FEM with Lagrange elements of higher degrees, for solving a number of linear and nonlinear PDEs. It is shown that the current improved ELM far outperforms the classical FEM. Its computational performance is comparable to that of the high-order FEM for smaller problem sizes, and for larger problem sizes the ELM markedly outperforms the high-order FEM.

categories: physics.comp-ph cs.LG cs.NA math.NA physics.flu-dyn


title: Data-Driven Representations for Testing Independence: Modeling, Analysis and Connection with Mutual Information Estimation

abstract: This work addresses testing the independence of two continuous and finite-dimensional random variables from the design of a data-driven partition. The empirical log-likelihood statistic is adopted to approximate the sufficient statistics of an oracle test against independence (that knows the two hypotheses). It is shown that approximating the sufficient statistics of the oracle test offers a learning criterion for designing a data-driven partition that connects with the problem of mutual information estimation. Applying these ideas in the context of a data-dependent tree-structured partition (TSP), we derive conditions on the TSP's parameters to achieve a strongly consistent distribution-free test of independence over the family of probabilities equipped with a density. Complementing this result, we present finite-length results that show our TSP scheme's capacity to detect the scenario of independence structurally with the data-driven partition as well as new sampling complexity bounds for this detection. Finally, some experimental analyses provide evidence regarding our scheme's advantage for testing independence compared with some strategies that do not use data-driven representations.


title: Newtonian Mechanics Based Transient Stability PART V: Inner-group Machine

abstract: This paper analyzes the mechanisms of the inner-group machine. It is first clarified that the inner-group machine is created from the difference between the equivalent system and the original system. The inner-group machine stability is analyzed based on the machine paradigms. In particular, strict correlation between the inner-group machine trajectory and the inner-group machine transient energy conversion is established through the I-CR system modeling. Then, the transient characteristics of the inner-group machine are analyzed. It is clarified that the inner-group motions might be inseparable or separable, and the inner-group machine DLP will occur later than the EDLP and IDLP. Simulation results show that the severity of the original system cannot be simply evaluated through its equivalent system once any inner-group motion becomes fierce.


title: A novel multiobjective evolutionary algorithm based on decomposition and multi-reference points strategy

abstract: Many real-world optimization problems such as engineering design can be eventually modeled as the corresponding multiobjective optimization problems (MOPs) which must be solved to obtain approximate Pareto optimal fronts. Multiobjective evolutionary algorithm based on decomposition (MOEA/D) has been regarded as a significantly promising approach for solving MOPs. Recent studies have shown that MOEA/D with uniform weight vectors is well-suited to MOPs with regular Pareto optimal fronts, but its performance in terms of diversity usually deteriorates when solving MOPs with irregular Pareto optimal fronts. In this way, the solution set obtained by the algorithm can not provide more reasonable choices for decision makers. In order to efficiently overcome this drawback, we propose an improved MOEA/D algorithm by virtue of the well-known Pascoletti-Serafini scalarization method and a new strategy of multi-reference points. Specifically, this strategy consists of the setting and adaptation of reference points generated by the techniques of equidistant partition and projection. For performance assessment, the proposed algorithm is compared with existing four state-of-the-art multiobjective evolutionary algorithms on benchmark test problems with various types of Pareto optimal fronts. According to the experimental results, the proposed algorithm exhibits better diversity performance than that of the other compared algorithms. Finally, our algorithm is applied to two real-world MOPs in engineering optimization successfully.


title: The top quark chromomagnetic dipole moment in the SM from the 4-body vertex function

abstract: A new proposal to compute the anomalous chromomagnetic dipole moment of the top quark, $\hat{\mu}_t$, in the Standard Model is presented. On the basis of the 5-dimensional effective Lagrangian operator that characterizes the quantum-loop induced chromodipolar vertices $gt\bar{t}$ and $ggt\bar{t}$, the $\hat{\mu}_t$ anomaly is derived via radiative correction at the 1-loop level from the non-Abelian 4-body vertex function $ggt\bar{t}$. We evaluate $\hat{\mu}_t(s)$ as a function of the energy scale $s=\pm E^2$, for $E=[10,1000]$ GeV, taking into account the running of the quark masses and alpha strong through the $\overline{\mathrm{MS}}$ scheme. In particular, we find that at the typical energy scale $E=m_Z$ for high-energy physics, similarly to $\alpha_s(m_Z^2)$, $\alpha(m_Z^2)$ and $s_W(m_Z^2)$, the spacelike evaluation yields $\hat{\mu}_t(-m_Z^2)$ $=$ $-0.025$$+$$0.00384i$ and the timelike $\hat{\mu}_t(m_Z^2)$ $=$ $-0.0318$$-$$0.0106i$. This Re$\thinspace\hat{\mu}_t(-m_Z^2)$ $=$ $-0.025$ from $ggt\bar{t}$ is even closer to the experimental central value $\hat{\mu}_t^\mathrm{Exp}=$ $-0.024$, than that coming from the known 3-body vertex function $gt\bar{t}$, $-0.0224$. Once again, the Im$\thinspace\hat{\mu}_t(-m_Z^2)$ part is due to the contribution of virtual charged currents, just like in the $gt\bar{t}$ case. We can infer that the spacelike prediction is the favored one.


title: Practical quantum access network over a 10 Gbit/s Ethernet passive optical network

abstract: Quantum key distribution (QKD) provides an information-theoretically secure method to share keys between legitimate users. To achieve large-scale deployment of QKD, it should be easily scalable and cost-effective. The infrastructure construction of quantum access network (QAN) expands network capacity and the integration between QKD and classical optical communications reduces the cost of channel. Here, we present a practical downstream QAN over a 10 Gbit/s Ethernet passive optical network (10G-EPON), which can support up to 64 users. In the full coexistence scheme using the single feeder fiber structure, the co-propagation of QAN and 10G-EPON signals with 9 dB attenuation is achieved over 21 km fiber, and the secure key rate for each of 16 users reaches 1.5 kbps. In the partial coexistence scheme using the dual feeder fiber structure, the combination of QAN and full-power 10G-EPON signals is achieved over 11 km with a network capacity of 64-user. The practical QAN over the 10G-EPON in our work implements an important step towards the achievement of large-scale QKD infrastructure.


title: Constrained Optimization Involving Nonconvex $\ell_p$ Norms: Optimality Conditions, Algorithm and Convergence

abstract: This paper investigates the optimality conditions for characterizing the local minimizers of the constrained optimization problems involving an $\ell_p$ norm ($0<p<1$) of the variables, which may appear in either the objective or the constraint. This kind of problems have strong applicability to a wide range of areas since usually the $\ell_p$ norm can promote sparse solutions. However, the nonsmooth and non-Lipschtiz nature of the $\ell_p$ norm often cause these problems difficult to analyze and solve. We provide the calculation of the subgradients of the $\ell_p$ norm and the normal cones of the $\ell_p$ ball. For both problems, we derive the first-order necessary conditions under various constraint qualifications. We also derive the sequential optimality conditions for both problems and study the conditions under which these conditions imply the first-order necessary conditions. We point out that the sequential optimality conditions can be easily satisfied for iteratively reweighted algorithms and show that the global convergence can be easily derived using sequential optimality conditions.


title: Improving Cell-Free Massive MIMO Detection Performance via Expectation Propagation

abstract: Cell-free (CF) massive multiple-input multiple-output (M-MIMO) technology plays a prominent role in the beyond fifth-generation (5G) networks. However, designing a high performance CF M-MIMO detector is a challenging task due to the presence of pilot contamination which appears when the number of pilot sequences is smaller than the number of users. This work proposes a CF M-MIMO detector referred to as CF expectation propagation (CF-EP) that incorporates the pilot contamination when calculating the posterior belief. The simulation results show that the proposed detector achieves significant improvements in terms of the bit-error rate and sum spectral efficiency performances as compared to the ones of the state-of-the-art CF detectors.


title: On the variational properties of the prescribed Ricci curvature functional

abstract: We study the prescribed Ricci curvature problem for homogeneous metrics. Given a (0,2)-tensor field $T$, this problem asks for solutions to the equation $\mathrm{Ric}(g)=cT$ for some constant $c$. Our approach is based on examining global properties of the scalar curvature functional whose critical points are solutions to this equation. We produce conditions for a general homogeneous space under which it has a global maximum. Finally, we study the behavior of the functional in specific examples to illustrate our result.


title: White-light Continuum Observation of the Off-limb Loops of the SOL2017-09-10 X8.2 Flare: Temporal and Spatial Variations

abstract: Observations of the Sun's off-limb white-light (WL) flares offer rare opportunities to study the energy release and transport mechanisms in flare loops. One of the best such events was SOL2017-09-10, an X8.2 flare that occurred near the Sun's west limb on 2017 September 10 and produced a WL loop system lasting more than 60 minutes and reaching an altitude higher than 30 Mm. The event was well observed by a suite of ground- and space-based instruments, including the Solar Dynamics Observatory/Helioseismic and Magnetic Imager (SDO/HMI) that captured its off-limb loops in WL continuum near Fe I 6173 A, and the Atmospheric Imager Assembly (SDO/AIA) that observed its ultraviolet (UV) and extreme-ultraviolet (EUV) counterparts. We found quasi-periodic pulsations in the WL and UV emissions at the flare loop-top with a period around 8.0 min. Each pulsation appears to have an EUV counterpart that occurs earlier in time and higher in altitude. Despite many similarities in the WL and UV images and light curves, the WL flux at the loop-top continues to grow for about 16 minutes while the UV fluxes gradually decay. We discuss the implication of these unprecedented observations on the understanding of the enigmatic off-limb WL flare emission mechanisms.


title: Temporal Knowledge Distillation for On-device Audio Classification

abstract: Improving the performance of on-device audio classification models remains a challenge given the computational limits of the mobile environment. Many studies leverage knowledge distillation to boost predictive performance by transferring the knowledge from large models to on-device models. However, most lack a mechanism to distill the essence of the temporal information, which is crucial to audio classification tasks, or similar architecture is often required. In this paper, we propose a new knowledge distillation method designed to incorporate the temporal knowledge embedded in attention weights of large transformer-based models into on-device models. Our distillation method is applicable to various types of architectures, including the non-attention-based architectures such as CNNs or RNNs, while retaining the original network architecture during inference. Through extensive experiments on both an audio event detection dataset and a noisy keyword spotting dataset, we show that our proposed method improves the predictive performance across diverse on-device architectures.


title: Angular-spectrum-dependent interference

abstract: Optical interference is not only a fundamental phenomenon that has enabled new theories of light to be derived but it has also been used in interferometry for the measurement of small displacements, refractive index changes and surface irregularities. In a two-beam interferometer, variations in the interference fringes are used as a diagnostic for anything that causes the optical path difference (OPD) to change; therefore, for a specified OPD, greater variation in the fringes indicates better measurement sensitivity. Here, we introduce and experimentally validate an interesting optical interference phenomenon that uses photons with a structured frequency-angular spectrum, which are generated from a spontaneous parametric down-conversion process in a nonlinear crystal. This interference phenomenon is manifested as interference fringes that vary much more rapidly with increasing OPD than the corresponding fringes for equal-inclination interference; the phenomenon is parameterised using an equivalent wavelength, which under our experimental conditions is 29.38 nm or about 1/27 of the real wavelength. This phenomenon not only enriches the knowledge with regard to optical interference but also offers promise for applications in interferometry.


title: Search for the charmed baryonium and dibaryon structures via the QCD sum rules

abstract: In the present work, we construct eight six-quark currents to study the $\Sigma_c\bar{\Sigma}_c$ baryonium and $\Sigma_c\Sigma_c$ dibaryon states via the QCD sum rules. For either $\Sigma_c\bar{\Sigma}_c$ baryonium or $\Sigma_c\Sigma_c$ dibaryon states, we construct four currents with the $J^P=1^-, 1^+, 0^+$ and $0^-$. Except the current $1^+$ of $\Sigma_c\Sigma_c$, we find the Borel windows for the other seven. There are two possible $\Sigma_c\bar{\Sigma}_c$ molecular states or compact six-quark states with the $J^P=1^-$ and $0^-$. The results of the current $0^+$ of $\Sigma_c\Sigma_c$ show a possible $\Sigma_c\Sigma_c$ molecular state with $20$ MeV binding energy. We find two $\Sigma_c\bar{\Sigma}_c$ baryonium resonances states with the $J^P=1^+, 0^+$ and two $\Sigma_c\Sigma_c$ dibaryon resonances states with the $J^P=0^-, 1^-$, respectively. Pole residues of the seven states are also calculated.


title: Uncertainty regions of observables and state-independent uncertainty relations

abstract: The optimal state-independent lower bounds for the sum of variances or deviations of observables are of significance for the growing number of experiments that reach the uncertainty limited regime. We present a framework for computing the tight uncertainty relations of variance or deviation via determining the uncertainty regions, which are formed by the tuples of two or more of quantum observables in random quantum states induced from the uniform Haar measure on the purified states. From the analytical formulae of these uncertainty regions, we present state-independent uncertainty inequalities satisfied by the sum of variances or deviations of two, three and arbitrary many observables, from which experimentally friend entanglement detection criteria are derived for bipartite and tripartite systems.


title: A simple numerical experiment on the dust temperature bias for Lyman break galaxies at $z\gtrsim 5$

abstract: Some studies suggest that the dust temperatures ($T_\mathrm{d}$) in high-redshift ($z\gtrsim 5$) Lyman break galaxies (LBGs) are high. However, possible observational bias in $T_\mathrm{d}$ is yet to be understood. Thus, we perform a simple test using random realizations of LBGs with various stellar masses, dust temperatures, and dust-to-stellar mass ratios, and examine how the sample detected by ALMA is biased in terms of $T_\mathrm{d}$. We show that ALMA tends to miss high-$T_\mathrm{d}$ objects even at total dust luminosity. LBGs are, however, basically selected by the stellar UV luminosity. The dust-temperature bias in a UV-selected sample is complicated because of the competing effects between high $T_\mathrm{d}$ and low dust abundance. For ALMA Band 6, there is no tendency of high-$T_\mathrm{d}$ LBGs being more easily detected in our experiment. Thus, we suggest that the observed trend of high $T_\mathrm{d}$ in $z\gtrsim 5$ LBGs is real. We also propose that the 450 $\mu$m band is useful in further clarifying the dust temperatures. To overcome the current shallowness of 450 $\mu$m observations, we examine a future Antarctic 30-m class telescope with a suitable atmospheric condition for wavelengths $\lesssim 450~\mu$m, where the detection is not confusion-limited. We find that, with this telescope, an $L_\mathrm{IR}$-selected sample with $\log(L_\mathrm{IR}/\mathrm{L}_{\odot})>11$ is constructed for $z\gtrsim 5$, and detection in the intermediate-$M_\star$ (stellar mass) range [$9<\log (M_\star /M_{\odot})<9.5$] is much improved, especially at high $T_\mathrm{d}$.


title: Numerical and Statistical Analysis of Aliquot Sequences

abstract: We present a variety of numerical data related to the growth of terms in aliquot sequences, iterations of the function $s(n) = \sigma(n) - n$. First, we compute the geometric mean of the ratio $s_k(n)/s_{k-1}(n)$ of $k$th iterates for $n \leq 2^{37}$ and $k=1,\dots,10.$ Second, we extend the computation of numbers not in the range of $s(n)$ (called untouchable) by Pollack and Pomerance to the bound of $2^{40}$ and use these data to compute the geometric mean of the ratio of consecutive terms limited to terms in the range of $s(n).$ Third, we give an algorithm to compute $k$-untouchable numbers ($k-1$st iterates of $s(n)$ but not $k$th iterates) along with some numerical data. Finally, inspired by earlier work of Devitt, we estimate the growth rate of terms in aliquot sequences using a Markov chain model based on data extracted from thousands of sequences.


title: Relationship Oriented Affordance Learning through Manipulation Graph Construction

abstract: In this paper, we propose Manipulation Relationship Graph (MRG), a novel affordance representation which captures the underlying manipulation relationships of an arbitrary scene. To construct such a graph from raw visual observations, a deep nerual network named AR-Net is introduced. It consists of an Attribute module and a Context module, which guide the relationship learning at object and subgraph level respectively. We quantitatively validate our method on a novel manipulation relationship dataset named SMRD. To evaluate the performance of the proposed model and representation, both visual perception and physical manipulation experiments are conducted. Overall, AR-Net along with MRG outperforms all baselines, achieving the success rate of 88.89% on task relationship recognition (TRR) and 73.33% on task completion (TC)


title: A Linear Bayesian Learning Receiver Scheme for Massive MIMO Systems

abstract: Much stringent reliability and processing latency requirements in ultra-reliable-low-latency-communication (URLLC) traffic make the design of linear massive multiple-input-multiple-output (M-MIMO) receivers becomes very challenging. Recently, Bayesian concept has been used to increase the detection reliability in minimum-mean-square-error (MMSE) linear receivers. However, the latency processing time is a major concern due to the exponential complexity of matrix inversion operations in MMSE schemes. This paper proposes an iterative M-MIMO receiver that is developed by using a Bayesian concept and a parallel interference cancellation (PIC) scheme, referred to as a linear Bayesian learning (LBL) receiver. PIC has a linear complexity as it uses a combination of maximum ratio combining (MRC) and decision statistic combining (DSC) schemes to avoid matrix inversion operations. Simulation results show that the bit-error-rate (BER) and latency processing performances of the proposed receiver outperform the ones of MMSE and best Bayesian-based receivers by minimum $2$ dB and $19$ times for various M-MIMO system configurations.


title: Closing the Gap Between Time-Domain Multi-Channel Speech Enhancement on Real and Simulation Conditions

abstract: The deep learning based time-domain models, e.g. Conv-TasNet, have shown great potential in both single-channel and multi-channel speech enhancement. However, many experiments on the time-domain speech enhancement model are done in simulated conditions, and it is not well studied whether the good performance can generalize to real-world scenarios. In this paper, we aim to provide an insightful investigation of applying multi-channel Conv-TasNet based speech enhancement to both simulation and real data. Our preliminary experiments show a large performance gap between the two conditions in terms of the ASR performance. Several approaches are applied to close this gap, including the integration of multi-channel Conv-TasNet into the beamforming model with various strategies, and the joint training of speech enhancement and speech recognition models. Our experiments on the CHiME-4 corpus show that our proposed approaches can greatly reduce the speech recognition performance discrepancy between simulation and real data, while preserving the strong speech enhancement capability in the frontend.


title: Polyhedral realizations for $B(\infty)$ and extended Young diagrams, Young walls of type ${\rm A}^{(1)}_{n-1}$, ${\rm C}^{(1)}_{n-1}$, ${\rm A}^{(2)}_{2n-2}$, ${\rm D}^{(2)}_{n}$

abstract: The crystal bases are quite useful combinatorial tools to study the representations of quantized universal enveloping algebras $U_q(\mathfrak{g})$. The polyhedral realization for $B(\infty)$ is a combinatorial description of the crystal base, which is defined as an image of embedding $\Psi_{\iota}:B(\infty)\hookrightarrow \mathbb{Z}^{\infty}_{\iota}$, where $\iota$ is an infinite sequence of indices and $\mathbb{Z}^{\infty}_{\iota}$ is an infinite $\mathbb{Z}$-lattice with a crystal structure associated with $\iota$. It is a natural problem to find an explicit form of the polyhedral realization ${\rm Im}(\Psi_{\iota})$. In this article, supposing that $\mathfrak{g}$ is of affine type ${\rm A}^{(1)}_{n-1}$, ${\rm C}^{(1)}_{n-1}$, ${\rm A}^{(2)}_{2n-2}$ or ${\rm D}^{(2)}_{n}$ and $\iota$ satisfies the condition of `adaptedness', we describe ${\rm Im}(\Psi_{\iota})$ by using several combinatorial objects such as extended Young diagrams and Young walls.


title: Modified Epidemic Diffusive Process on the Apollonian Network

abstract: We present an analysis of an epidemic spreading process on the Apollonian network that can describe an epidemic spreading in a non-sedentary population. The modified diffusive epidemic process was employed in this analysis in a computational context by means of the Monte Carlo method. Our model has been useful for modeling systems closer to reality consisting of two classes of individuals: susceptible (A) and infected (B). The individuals can diffuse in a network according to constant diffusion rates $D_{A}$ and $D_{B}$, for the classes A and B, respectively, and obeying three diffusive regimes, i.e., $D_{A}<D_{B}$, $D_{A}=D_{B}$ and $D_{A}>D_{B}$. Into the same site $i$, the reaction occurs according to the dynamical rule based on Gillespie's algorithm. Finite-size scaling analysis has shown that our model exhibit continuous phase transition to an absorbing state with a set of critical exponents given by $\beta/\nu=0.66(1)$, $1/\nu=0.46(2)$, and $\gamma/\nu=-0.24(2)$ common to every investigated regime. In summary, the continuous phase transition, characterized by this set of critical exponents, does not have the same exponents of the Mean-Field universality class in both regular lattices and complex networks.


title: Separating Long-Form Speech with Group-Wise Permutation Invariant Training

abstract: Multi-talker conversational speech processing has drawn many interests for various applications such as meeting transcription. Speech separation is often required to handle overlapped speech that is commonly observed in conversation. Although the original utterancelevel permutation invariant training-based continuous speech separation approach has proven to be effective in various conditions, it lacks the ability to leverage the long-span relationship of utterances and is computationally inefficient due to the highly overlapped sliding windows. To overcome these drawbacks, we propose a novel training scheme named Group-PIT, which allows direct training of the speech separation models on the long-form speech with a low computational cost for label assignment. Two different speech separation approaches with Group-PIT are explored, including direct long-span speech separation and short-span speech separation with long-span tracking. The experiments on the simulated meeting-style data demonstrate the effectiveness of our proposed approaches, especially in dealing with a very long speech input.


title: SOAT: A Scene- and Object-Aware Transformer for Vision-and-Language Navigation

abstract: Natural language instructions for visual navigation often use scene descriptions (e.g., "bedroom") and object references (e.g., "green chairs") to provide a breadcrumb trail to a goal location. This work presents a transformer-based vision-and-language navigation (VLN) agent that uses two different visual encoders -- a scene classification network and an object detector -- which produce features that match these two distinct types of visual cues. In our method, scene features contribute high-level contextual information that supports object-level processing. With this design, our model is able to use vision-and-language pretraining (i.e., learning the alignment between images and text from large-scale web data) to substantially improve performance on the Room-to-Room (R2R) and Room-Across-Room (RxR) benchmarks. Specifically, our approach leads to improvements of 1.8% absolute in SPL on R2R and 3.7% absolute in SR on RxR. Our analysis reveals even larger gains for navigation instructions that contain six or more object references, which further suggests that our approach is better able to use object features and align them to references in the instructions.


title: Physically Explainable CNN for SAR Image Classification

abstract: Integrating the special electromagnetic characteristics of Synthetic Aperture Radar (SAR) in deep neural networks is essential in order to enhance the explainability and physics awareness of deep learning. In this paper, we first propose a novel physically explainable convolutional neural network for SAR image classification, namely physics guided and injected learning (PGIL). It comprises three parts: (1) explainable models (XM) to provide prior physics knowledge, (2) physics guided network (PGN) to encode the knowledge into physics-aware features, and (3) physics injected network (PIN) to adaptively introduce the physics-aware features into classification pipeline for label prediction. A hybrid Image-Physics SAR dataset format is proposed for evaluation, with both Sentinel-1 and Gaofen-3 SAR data being experimented. The results show that the proposed PGIL substantially improve the classification performance in case of limited labeled data compared with the counterpart data-driven CNN and other pre-training methods. Additionally, the physics explanations are discussed to indicate the interpretability and the physical consistency preserved in the predictions. We deem the proposed method would promote the development of physically explainable deep learning in SAR image interpretation field.


title: On branching laws of Speh representations

abstract: In this paper, we consider the branching law of the Speh representation $\mathrm{Sp}(\pi,n+l)$ of $\mathrm{GL}_{2n+2l}$ with respect to the block diagonal subgroup $\mathrm{GL}_n\times\mathrm{GL}_{n+2l}$ for any irreducible generic representation $\pi$ of $\mathrm{GL}_2$ over any $p$-adic field. We use the Shalika model of $\mathrm{Sp}(\pi,n)$ to construct certain zeta integrals, which were defined by Ginzburg and Kaplan independently, and study them. Finally, using these zeta integrals, we obtain a nonzero $\mathrm{GL}_n\times\mathrm{GL}_{n+2l}$-map from $\mathrm{Sp}(\pi,n+l)$ to $\tau\boxtimes\tau^\vee\chi_\pi\times\mathrm{Sp}(\pi, l)$ for any irreducible representation $\tau$ of $\mathrm{GL}_n$. These results form part of the local theory of the Miyawaki lifting for unitary groups.


title: Tail contributions to gravitational conservative dynamics

abstract: We compute tail contributions to the conservative dynamics of a generic self-gravitating system, for every multipole order, of either electric and magnetic parity. Such contributions arise when gravitational radiation is backscattered by the static curvature generated by the source itself and reabsorbed. Tail processes are relevant for the study of compact binary dynamics starting from the fourth post-Newtonian order. Within this context, we resolve a discrepancy recently arised at fifth post-Newtonian order between Effective Field Theory and self-force results.


title: Image Comes Dancing with Collaborative Parsing-Flow Video Synthesis

abstract: Transferring human motion from a source to a target person poses great potential in computer vision and graphics applications. A crucial step is to manipulate sequential future motion while retaining the appearance characteristic.Previous work has either relied on crafted 3D human models or trained a separate model specifically for each target person, which is not scalable in practice.This work studies a more general setting, in which we aim to learn a single model to parsimoniously transfer motion from a source video to any target person given only one image of the person, named as Collaborative Parsing-Flow Network (CPF-Net). The paucity of information regarding the target person makes the task particularly challenging to faithfully preserve the appearance in varying designated poses. To address this issue, CPF-Net integrates the structured human parsing and appearance flow to guide the realistic foreground synthesis which is merged into the background by a spatio-temporal fusion module. In particular, CPF-Net decouples the problem into stages of human parsing sequence generation, foreground sequence generation and final video generation. The human parsing generation stage captures both the pose and the body structure of the target. The appearance flow is beneficial to keep details in synthesized frames. The integration of human parsing and appearance flow effectively guides the generation of video frames with realistic appearance. Finally, the dedicated designed fusion network ensure the temporal coherence. We further collect a large set of human dancing videos to push forward this research field. Both quantitative and qualitative results show our method substantially improves over previous approaches and is able to generate appealing and photo-realistic target videos given any input person image. All source code and dataset will be released at https://github.com/xiezhy6/CPF-Net.


title: Uniform Concentration Bounds toward a Unified Framework for Robust Clustering

abstract: Recent advances in center-based clustering continue to improve upon the drawbacks of Lloyd's celebrated $k$-means algorithm over $60$ years after its introduction. Various methods seek to address poor local minima, sensitivity to outliers, and data that are not well-suited to Euclidean measures of fit, but many are supported largely empirically. Moreover, combining such approaches in a piecemeal manner can result in ad hoc methods, and the limited theoretical results supporting each individual contribution may no longer hold. Toward addressing these issues in a principled way, this paper proposes a cohesive robust framework for center-based clustering under a general class of dissimilarity measures. In particular, we present a rigorous theoretical treatment within a Median-of-Means (MoM) estimation framework, showing that it subsumes several popular $k$-means variants. In addition to unifying existing methods, we derive uniform concentration bounds that complete their analyses, and bridge these results to the MoM framework via Dudley's chaining arguments. Importantly, we neither require any assumptions on the distribution of the outlying observations nor on the relative number of observations $n$ to features $p$. We establish strong consistency and an error rate of $O(n^{-1/2})$ under mild conditions, surpassing the best-known results in the literature. The methods are empirically validated thoroughly on real and synthetic datasets.


title: Diversity Matters When Learning From Ensembles

abstract: Deep ensembles excel in large-scale image classification tasks both in terms of prediction accuracy and calibration. Despite being simple to train, the computation and memory cost of deep ensembles limits their practicability. While some recent works propose to distill an ensemble model into a single model to reduce such costs, there is still a performance gap between the ensemble and distilled models. We propose a simple approach for reducing this gap, i.e., making the distilled performance close to the full ensemble. Our key assumption is that a distilled model should absorb as much function diversity inside the ensemble as possible. We first empirically show that the typical distillation procedure does not effectively transfer such diversity, especially for complex models that achieve near-zero training error. To fix this, we propose a perturbation strategy for distillation that reveals diversity by seeking inputs for which ensemble member outputs disagree. We empirically show that a model distilled with such perturbed samples indeed exhibits enhanced diversity, leading to improved performance.


title: Training Wasserstein GANs without gradient penalties

abstract: We propose a stable method to train Wasserstein generative adversarial networks. In order to enhance stability, we consider two objective functions using the $c$-transform based on Kantorovich duality which arises in the theory of optimal transport. We experimentally show that this algorithm can effectively enforce the Lipschitz constraint on the discriminator while other standard methods fail to do so. As a consequence, our method yields an accurate estimation for the optimal discriminator and also for the Wasserstein distance between the true distribution and the generated one. Our method requires no gradient penalties nor corresponding hyperparameter tuning and is computationally more efficient than other methods. At the same time, it yields competitive generators of synthetic images based on the MNIST, F-MNIST, and CIFAR-10 datasets.


title: Voxelization based packing analysis for discrete element simulations of non-spherical particles

abstract: A voxelization based post-processing algorithm is proposed to analyze the packing of non-spherical particle assemblies simulated using the Discrete Element Method. Voxelization of the particle data allows for isolating the geometric features of the granular assembly in various spatial sub-domains (2D surface or 3D region) and investigate the localized packing behaviour. Analyzing the local packing behaviour enables determining con-fined influences such as the wall-effect, stacking behaviour, local expansion/contraction and localized loading. The efficacy of the proposed technique to analyze practical granular assemblies is demonstrated through the packing analysis of different assemblies of superquadric cubes, superquadric ellipsoidals, and multi-spherical coffee beans.


title: A New Semi-Empirical Model for Cosmic Ray Muon Flux Estimation

abstract: Cosmic ray muons have emerged as a non-conventional high-energy radiation probe to monitor dense and large objects. Muons are the most abundant cosmic radiation on Earth, however, their flux at sea level is approximately 10,000 min^-1m^-2 much less than that of induced radiation. In addition, cosmic ray muon flux depends on not only various natural conditions (e.g., zenith angle, altitude, or solar activities) but also the geometric characteristic of detectors. Since the low muon flux typically results in long measurement times, an accurate estimation of measurable muon counts is important for muon applications. Here we propose a simple and versatile semi-empirical model to improve the accuracy in muon flux estimation at all zenith angles by incorporating the geometric parameters of detectors, and we name this the Effective Solid Angle model. To demonstrate the functionality, our model is compared with i) the cosine-squared, ii) PARMA model, and iii) Monte-Carlo simulations, and iv) measurements. Our results show that the muon count rate estimation capability is significantly improved resulting in increasing a mean C/E from 0.7 to 0.95. By selecting an appropriate intensity correlation, the model can be easily extended to estimate muon flux at various altitude and underground level.

categories: astro-ph.IM astro-ph.HE hep-ex nucl-ex physics.ins-det


title: Differentially Private Federated Bayesian Optimization with Distributed Exploration

abstract: Bayesian optimization (BO) has recently been extended to the federated learning (FL) setting by the federated Thompson sampling (FTS) algorithm, which has promising applications such as federated hyperparameter tuning. However, FTS is not equipped with a rigorous privacy guarantee which is an important consideration in FL. Recent works have incorporated differential privacy (DP) into the training of deep neural networks through a general framework for adding DP to iterative algorithms. Following this general DP framework, our work here integrates DP into FTS to preserve user-level privacy. We also leverage the ability of this general DP framework to handle different parameter vectors, as well as the technique of local modeling for BO, to further improve the utility of our algorithm through distributed exploration (DE). The resulting differentially private FTS with DE (DP-FTS-DE) algorithm is endowed with theoretical guarantees for both the privacy and utility and is amenable to interesting theoretical insights about the privacy-utility trade-off. We also use real-world experiments to show that DP-FTS-DE achieves high utility (competitive performance) with a strong privacy guarantee (small privacy loss) and induces a trade-off between privacy and utility.


title: The ASTRID simulation: the evolution of Supermassive Black Holes

abstract: We present the evolution of black holes (BHs) and their relationship with their host galaxies in Astrid, a large-volume cosmological hydrodynamical simulation with box size 250 $h^{-1} \rm Mpc$ containing $2\times5500^3$ particles evolved to z=3. Astrid statistically models BH gas accretion and AGN feedback to their environments, applies a power-law distribution for BH seed mass $M_{\rm sd}$, uses a dynamical friction model for BH dynamics and executes a physical treatment of BH mergers. The BH population is broadly consistent with empirical constraints on the BH mass function, the bright end of the luminosity functions, and the time evolution of BH mass and accretion rate density. The BH mass and accretion exhibit a tight correlation with host stellar mass and star formation rate. We trace BHs seeded before z>10 down to z=3, finding that BHs carry virtually no imprint of the initial $M_{\rm sd}$ except those with the smallest $M_{\rm sd}$, where less than 50\% of them have doubled in mass. Gas accretion is the dominant channel for BH growth compared to BH mergers. With dynamical friction, Astrid predicts a significant delay for BH mergers after the first encounter of a BH pair, with a typical elapse time of about 200 Myrs. There are in total $4.5 \times 10^5$ BH mergers in Astrid at z>3, $\sim 10^3$ of which have X-ray detectable EM counterparts: a bright kpc scale dual AGN with $L_X>10^{43}$ erg/s. BHs with $M_{\rm BH} \sim 10^{7-8} M_{\odot}$ experience the most frequent mergers. Galaxies that host BH mergers are unbiased tracers of the overall $M_{\rm BH} - M_{*}$ relation. Massive ($>10^{11} M_{\odot}$) galaxies have a high occupation number (>10) of BHs, and hence host the majority of BH mergers.


title: Unusual magnetic and transport properties in HoMn$_6$Sn$_6$ kagome magnet

abstract: With intricate lattice structures, kagome materials are an excellent platform to study various fascinating topological quantum states. In particular, kagome materials, revealing large responses to external stimuli such as pressure or magnetic field, are subject to special investigation. Here, we study the kagome-net HoMn$_6$Sn$_6$ magnet that undergoes paramagnetic to ferrimagnetic transition (below 376 K) and reveals spin-reorientation transition below 200 K. In this compound, we observe the topological Hall effect and substantial contribution of anomalous Hall effect above 100 K. We unveil the pressure effects on magnetic ordering at a low magnetic field from the pressure tunable magnetization measurement. By utilizing high-resolution angle-resolved photoemission spectroscopy, Dirac-like dispersion at the high-symmetry point K is revealed in the vicinity of the Fermi level, which is well supported by the first-principles calculations, suggesting a possible Chern-gapped Dirac cone in this compound. Our investigation will pave the way to understand the magneto-transport and electronic properties of various rare-earth-based kagome magnets.


title: Divisibility of certain $\ell$-regular partitions by $2$

abstract: For a positive integer $\ell$, let $b_{\ell}(n)$ denote the number of $\ell$-regular partitions of a nonnegative integer $n$. Motivated by some recent conjectures of Keith and Zanello, we establish infinite families of congruences modulo $2$ for $b_3(n)$ and $b_{21}(n)$. We prove a specific case of a conjecture of Keith and Zanello on self-similarities of $b_3(n)$ modulo $2$. We next prove that the series $\sum_{n=0}^{\infty}b_9(2n+1)q^n$ is lacunary modulo arbitrary powers of $2$. We also prove that the series $\sum_{n=0}^{\infty}b_9(4n)q^n$ is lacunary modulo $2$.