title: Speech2Phone: A Novel and Efficient Method for Training Speaker Recognition Models

abstract: In this paper we present an efficient method for training models for speaker recognition using small or under-resourced datasets. This method requires less data than other SOTA (State-Of-The-Art) methods, e.g. the Angular Prototypical and GE2E loss functions, while achieving similar results to those methods. This is done using the knowledge of the reconstruction of a phoneme in the speaker's voice. For this purpose, a new dataset was built, composed of 40 male speakers, who read sentences in Portuguese, totaling approximately 3h. We compare the three best architectures trained using our method to select the best one, which is the one with a shallow architecture. Then, we compared this model with the SOTA method for the speaker recognition task: the Fast ResNet-34 trained with approximately 2,000 hours, using the loss functions Angular Prototypical and GE2E. Three experiments were carried out with datasets in different languages. Among these three experiments, our model achieved the second best result in two experiments and the best result in one of them. This highlights the importance of our method, which proved to be a great competitor to SOTA speaker recognition models, with 500x less data and a simpler approach.


title: An astrophysically motivated ranking criterion for low-latency electromagnetic follow-up of gravitational wave events

abstract: We investigate the properties of the host galaxies of compact binary mergers across cosmic time. To this end, we combine population synthesis simulations together with galaxy catalogues from the hydrodynamical cosmological simulation EAGLE to derive the properties of the host galaxies of binary neutron star (BNS), black hole-neutron star (BHNS) and binary black hole (BBH) mergers. Within this framework, we derive the host galaxy probability, i.e., the probability that a galaxy hosts a compact binary coalescence as a function of its stellar mass, star formation rate, $K_s$ magnitude and $B$ magnitude. This quantity is particularly important for low-latency searches of gravitational wave (GW) sources as it provides a way to rank galaxies lying inside the credible region in the sky of a given GW detection, hence reducing the number of viable host candidates. Furthermore, even if no electromagnetic counterpart is detected, the proposed ranking criterion can still be used to classify the galaxies contained in the error box. Our results show that massive galaxies (or equivalently galaxies with a high luminosity in $K_s$ band) have a higher probability of hosting BNS, BHNS, and BBH mergers. We provide the probabilities in a suitable format to be implemented in future low-latency searches.


title: EmbPred30: Assessing 30-days Readmission for Diabetic Patients using Categorical Embeddings

abstract: Hospital readmission is a crucial healthcare quality measure that helps in determining the level of quality of care that a hospital offers to a patient and has proven to be immensely expensive. It is estimated that more than $25 billion are spent yearly due to readmission of diabetic patients in the USA. This paper benchmarks existing models and proposes a new embedding based state-of-the-art deep neural network(DNN). The model can identify whether a hospitalized diabetic patient will be readmitted within 30 days or not with an accuracy of 95.2% and Area Under the Receiver Operating Characteristics(AUROC) of 97.4% on data collected from 130 US hospitals between 1999-2008. The results are encouraging with patients having changes in medication while admitted having a high chance of getting readmitted. Identifying prospective patients for readmission could help the hospital systems in improving their inpatient care, thereby saving them from unnecessary expenditures.


title: Plasma Temperature Measurement with a Silicon Photomultiplier (SiPM)

abstract: The temperature of a nonneutral plasma confined in a Penning-Malmberg trap can be determined by slowly lowering one side of the trap's electrostatic axial confinement barrier; the temperature is inferred from the rate at which particles escape the trap as a function of the barrier height. Often, the escaping particles are directed toward a microchannel plate (MCP), and the resulting amplified charge is collected on a phosphor screen. The screen is used for imaging the plasma, but can also be used as a Faraday cup FC for a temperature measurement. The sensitivity limit is then set by microphonic noise enhanced by the screen's high voltage bias. Alternately, a silicon photomultiplier (SiPM) can be employed to measure the charge via the light emitted from the phosphor screen. This decouples the signal from the microphonic noise and allows the temperature of colder and smaller plasmas to be measured than could be measured previously; this paper focusses on the advantages of a SiPM over a FC.


title: Anneal-path correction in flux qubits

abstract: Quantum annealers require accurate control and optimized operation schemes to reduce noise levels, in order to eventually demonstrate a computational advantage over classical algorithms. We study a high coherence four-junction capacitively shunted flux qubit (CSFQ), using dispersive measurements to extract system parameters and model the device. Josephson junction asymmetry inherent to the device causes a deleterious nonlinear cross-talk when annealing the qubit. We implement a nonlinear annealing path to correct the asymmetry in-situ, resulting in a substantial increase in the probability of the qubit being in the correct state given an applied flux bias. We also confirm the multi-level structure of our CSFQ circuit model by annealing it through small spectral gaps and observing quantum signatures of energy level crossings. Our results demonstrate an anneal-path correction scheme designed and implemented to improve control accuracy for high-coherence and high-control quantum annealers, which leads to an enhancement of success probability in annealing protocols.


title: Heavily Hydride-ion-doped 1111-type Iron-based Superconductors: Synthesis, Physical Properties and Electronic Structure

abstract: Iron-based superconductors have grown to be a new continent of high Tc superconductors comparable to cuprates. The optimal critical temperature (Tc) of 56 K in electron-doped 1111-type iron oxypnictides attracts considerable attention of physicists and chemists. Carrier doping is not only essential to induce superconductivity but also is a critical parameter that governs the electronic, magnetic, and crystallographic properties of ground states in high-Tc superconductors. Hydride ion (H-) which is an anionic state of hydrogen acts as an efficient electron donor in the 1111-type, leading to several important discoveries such as two-dome structure of superconducting phase and bipartite parent phase. This article summarizes the synthesis, physical properties, and electronic structure of the H- bearing 1111-type iron-based superconductors along with the relevant phenomena of other superconductors. We show several common characteristics of iron-based superconductors with high-Tc over 50 K and suggest a way to achieve higher Tc.


title: Convex Geometry and Duality of Over-parameterized Neural Networks

abstract: We develop a convex analytic approach to analyze finite width two-layer ReLU networks. We first prove that an optimal solution to the regularized training problem can be characterized as extreme points of a convex set, where simple solutions are encouraged via its convex geometrical properties. We then leverage this characterization to show that an optimal set of parameters yield linear spline interpolation for regression problems involving one dimensional or rank-one data. We also characterize the classification decision regions in terms of a kernel matrix and minimum $\ell_1$-norm solutions. This is in contrast to Neural Tangent Kernel which is unable to explain predictions of finite width networks. Our convex geometric characterization also provides intuitive explanations of hidden neurons as auto-encoders. In higher dimensions, we show that the training problem can be cast as a finite dimensional convex problem with infinitely many constraints. Then, we apply certain convex relaxations and introduce a cutting-plane algorithm to globally optimize the network. We further analyze the exactness of the relaxations to provide conditions for the convergence to a global optimum. Our analysis also shows that optimal network parameters can be also characterized as interpretable closed-form formulas in some practically relevant special cases.


title: Style Transfer for Light Field Photography

abstract: As light field images continue to increase in use and application, it becomes necessary to adapt existing image processing methods to this unique form of photography. In this paper we explore methods for applying neural style transfer to light field images. Feed-forward style transfer networks provide fast, high-quality results for monocular images, but no such networks exist for full light field images. Because of the size of these images, current light field data sets are small and are insufficient for training purely feed-forward style-transfer networks from scratch. Thus, it is necessary to adapt existing monocular style transfer networks in a way that allows for the stylization of each view of the light field while maintaining visual consistencies between views. Instead, the proposed method backpropagates the loss through the network, and the process is iterated to optimize (essentially overfit) the resulting stylization for a single light field image alone. The network architecture allows for the incorporation of pre-trained fast monocular stylization networks while avoiding the need for a large light field training set.


title: Distributed Weighted Least-squares Estimation for Networked Systems with Edge Measurements

abstract: This paper studies the problem of distributed weighted least-squares (WLS) estimation for an interconnected linear measurement network with additive noise. Two types of measurements are considered: self measurements for individual nodes, and edge measurements for the connecting nodes. Each node in the network carries out distributed estimation by using its own measurement and information transmitted from its neighbours. We study two distributed estimation algorithms: a recently proposed distributed WLS algorithm and the so-called Gaussian Belief Propagation (BP) algorithm. We first establish the equivalence of the two algorithms. We then prove a key result which shows that the information matrix is always generalised diagonally dominant, under some very mild condition. Using these two results and some known convergence properties of the Gaussian BP algorithm, we show that the aforementioned distributed WLS algorithm gives the globally optimal WLS estimate asymptotically. A bound on its convergence rate is also presented.


title: Focal-shape effects on the efficiency of the tunnel-ionization probe for extreme laser intensities

abstract: We examine the effect of laser focusing on the potential of a recently discussed scheme [M.F. Ciappina {\em et al}, Phys. Rev. A 99, 043405 (2019); Las. Phys. Lett. 17, 025301 (2020)] for in situ determination of ultra-high intensities of electromagnetic radiation delivered by multi-petawatt laser facilities. Using two model intensity distributions in the focus of a laser beam, we show how the resulting yields of highly charged ions generated in the process of multiple sequential tunneling of electrons from atoms, depend on the shape of these distributions. Our findings lead to the conclusion that an accurate extraction of the peak laser intensity can be made either in the near-threshold regime when the production of the highest charge state happens only in a small part of the laser focus close to the point where the intensity is maximal, or through the determination of the points where the ion yields of close charges become equal. We show that, for realistic parameters of the gas target, the number of ions generated in the central part of the focus in the threshold regime should be sufficient for a reliable measurement with highly sensitive time-of-flight detectors. Although positions of the intersection points generally depend on the focal shape, they can be used to localize the peak intensity value in a certain interval. Additionally to this analysis, we discuss the method in comparison to other recently proposed approaches for direct measurement of extreme laser intensities.


title: Device Heterogeneity in Federated Learning: A Superquantile Approach

abstract: We propose a federated learning framework to handle heterogeneous client devices which do not conform to the population data distribution. The approach hinges upon a parameterized superquantile-based objective, where the parameter ranges over levels of conformity. We present an optimization algorithm and establish its convergence to a stationary point. We show how to practically implement it using secure aggregation by interleaving iterations of the usual federated averaging method with device filtering. We conclude with numerical experiments on neural networks as well as linear models on tasks from computer vision and natural language processing.


title: Large data existence theory for three-dimensional unsteady flows of rate-type viscoelastic fluids with stress diffusion

abstract: We prove that there exists a weak solution to a system governing an unsteady flow of a viscoelastic fluid in three dimensions, for arbitrarily large time interval and data. The fluid is described by the incompressible Navier-Stokes equations for the velocity $v$, coupled with a diffusive variant of a combination of the Oldroyd-B and the Giesekus models for a tensor $\mathbb{B}$. By a proper choice of the constitutive relations for the Helmholtz free energy (which, however, is non-standard in the current literature, despite the fact that this choice is well motivated from the point of view of physics) and for the energy dissipation, we are able to prove that $\mathbb{B}$ enjoys the same regularity as $v$ in the classical three-dimensional Navier-Stokes equations. This enables us to handle any kind of objective derivative of $\mathbb{B}$, thus obtaining existence results for the class of diffusive Johnson-Segalman models as well. Moreover, using a suitable approximation scheme, we are able to show that $\mathbb{B}$ remains positive definite if the initial datum was a positive definite matrix (in a pointwise sense). We also show how the model we are considering can be derived from basic balance equations and thermodynamical principles in a natural way.


title: A classification algorithm for time-domain novelties in preparation for LSST alerts: Application to variable stars and transients detected with DECam in the Galactic Bulge

abstract: With the advent of the Large Synoptic Survey Telescope (LSST), time-domain astronomy will be faced with an unprecedented volume and rate of data. Real-time processing of variables and transients detected by such large-scale surveys is critical to identifying the more unusual events and allocating scarce follow-up resources efficiently. We develop an algorithm to identify these novel events within a given population of variable sources. We determine the distributions of magnitude changes (dm) over time intervals (dt) for a given passband f, pf(dm|dt), and use these distributions to compute the likelihood of a test source being consistent with the population, or an outlier. We demonstrate our algorithm by applying it to the DECam multi-band time-series data of more than 2000 variable stars identified by Saha et al. (2019) in the Galactic Bulge that are largely dominated by long-period variables and pulsating stars. Our algorithm discovers 18 outlier sources in the sample, including a microlensing event, a dwarf nova, and two chromospherically active RS CVn stars, as well as sources in the Blue Horizontal Branch region of the color-magnitude diagram without any known counterparts. We compare the performance of our algorithm for novelty detection with multivariate KDE and Isolation Forest on the simulated PLAsTiCC dataset. We find that our algorithm yields comparable results despite its simplicity. Our method provides an efficient way for flagging the most unusual events in a real-time alert-broker system.


title: Deep Learning and Statistical Models for Time-Critical Pedestrian Behaviour Prediction

abstract: The time it takes for a classifier to make an accurate prediction can be crucial in many behaviour recognition problems. For example, an autonomous vehicle should detect hazardous pedestrian behaviour early enough for it to take appropriate measures. In this context, we compare the switching linear dynamical system (SLDS) and a three-layered bi-directional long short-term memory (LSTM) neural network, which are applied to infer pedestrian behaviour from motion tracks. We show that, though the neural network model achieves an accuracy of 80%, it requires long sequences to achieve this (100 samples or more). The SLDS, has a lower accuracy of 74%, but it achieves this result with short sequences (10 samples). To our knowledge, such a comparison on sequence length has not been considered in the literature before. The results provide a key intuition of the suitability of the models in time-critical problems.


title: Search on Vertex-Transitive Graphs by Lackadaisical Quantum Walk

abstract: The lackadaisical quantum walk is a discrete-time, coined quantum walk on a graph with a weighted self-loop at each vertex. It uses a generalized Grover coin and the flip-flop shift, which makes it equivalent to Szegedy's quantum Markov chain. It has been shown that a lackadaisical quantum walk can improve spatial search on the complete graph, discrete torus, cycle, and regular complete bipartite graph. In this paper, we observe that these are all vertex-transitive graphs, and when there is a unique marked vertex, the optimal weight of the self-loop equals the degree of the loopless graph divided by the total number of vertices. We propose that this holds for all vertex-transitive graphs with a unique marked vertex. We present a number of numerical simulations supporting this hypothesis, including search on periodic cubic lattices of arbitrary dimension, strongly regular graphs, Johnson graphs, and the hypercube.


title: Enforcing Mean Reversion in State Space Models for Prawn Pond Water Quality Forecasting

abstract: The contribution of this study is a novel approach to introduce mean reversion in multi-step-ahead forecasts of state-space models. This approach is demonstrated in a prawn pond water quality forecasting application. The mean reversion constrains forecasts by gradually drawing them to an average of previously observed dynamics. This corrects deviations in forecasts caused by irregularities such as chaotic, non-linear, and stochastic trends. The key features of the approach include (1) it enforces mean reversion, (2) it provides a means to model both short and long-term dynamics, (3) it is able to apply mean reversion to select structural state-space components, and (4) it is simple to implement. Our mean reversion approach is demonstrated on various state-space models and compared with several time-series models on a prawn pond water quality dataset. Results show that mean reversion reduces long-term forecast errors by over 60% to produce the most accurate models in the comparison.


title: A recursion for a symmetric function generalization of the $q$-Dyson constant term identity

abstract: In 2000, Kadell gave an orthogonality conjecture for a symmetric function generalization of the $q$-Dyson constant term identity or the Zeilberger--Bressoud $q$-Dyson theorem. The non-zero part of Kadell's orthogonality conjecture is a constant term identity indexed by a weak composition $v=(v_1,\dots,v_n)$ in the case when only one $v_i\neq 0$. This conjecture was first proved by K\'{a}rolyi, Lascoux and Warnaar in 2015. They further formulated a closed-form expression for the above mentioned constant term in the case when all the parts of $v$ are distinct. Recently we obtain a recursion for this constant term provided that the largest part of $v$ occurs with multiplicity one in $v$. In this paper, we generalize our previous result to all compositions $v$.


title: Crossed product Leavitt path algebras

abstract: If $E$ is a directed graph and $K$ is a field, the Leavitt path algebra $L_K(E)$ of $E$ over $K$ is naturally graded by the group of integers $\mathbb Z.$ We formulate properties of the graph $E$ which are equivalent with $L_K(E)$ being a crossed product, a skew group ring, or a group ring with respect to this natural grading. We state this main result so that the algebra properties of $L_K(E)$ are also characterized in terms of the pre-ordered group properties of the Grothendieck $\mathbb Z$-group of $L_K(E)$. If $E$ has finitely many vertices, we characterize when $L_K(E)$ is strongly graded in terms of the properties of $K_0^\Gamma(L_K(E)).$ Our proof also provides an alternative to the known proof of the equivalence $L_K(E)$ is strongly graded if and only if $E$ has no sinks for a finite graph $E.$ We also show that, if unital, the algebra $L_K(E)$ is strongly graded and graded unit-regular if and only if $L_K(E)$ is a crossed product. In the process of showing the main result, we obtain conditions on a group $\Gamma$ and a $\Gamma$-graded division ring $K$ equivalent with the requirements that a $\Gamma$-graded matrix ring $R$ over $K$ is strongly graded, a crossed product, a skew group ring, or a group ring. We characterize these properties also in terms of the action of the group $\Gamma$ on the Grothendieck $\Gamma$-group $K_0^\Gamma(R).$


title: Benchmarking calculations with spectroscopic accuracy of excitation energies and wavelengths in sulfur-like tungsten

abstract: Atomic properties of S-like W are evaluated through a state-of-the-art method, namely, the multi-configuration Dirac-Hartree-Fock (MCDHF) method combined with the relativistic configuration interaction (RCI) approach. The level energies, wavelengths, and transition parameters involving the 88 lowest levels of W$^{+58}$ (W LIX) are calculated. We discuss in detail the relative importance of the valence- and core-valence electron correlation effects, the Breit interaction, the higher order retardation correction beyond the Breit interaction through the transverse photon interaction, and the quantum electrodynamical (QED) corrections. The present level energies are highly accurate, with uncertainties close to what can be achieved from spectroscopy. As such, they provide benchmark tests for other theoretical calculations of S-like W and should assist the spectroscopists in their assignment/identification of observed lines in complex spectra.


title: Residual finiteness for central pushouts

abstract: We prove that pushouts $A*_CB$ of residually finite-dimensional (RFD) $C^*$-algebras over central subalgebras are always residually finite-dimensional provided the fibers $A_p$ and $B_p$, $p\in \mathrm{spec}~C$ are RFD, recovering and generalizing results by Korchagin and Courtney-Shulman. This then allows us to prove that certain central pushouts of amenable groups have RFD group $C^*$-algebras. Along the way, we discuss the problem of when, given a central group embedding $H\le G$, the resulting $C^*$-algebra morphism is a continuous field: this is always the case for amenable $G$ but not in general.


title: Heavy Quark Energy Loss in the Quark-Gluon Plasma in the Moller theory

abstract: We study the energy loss of a heavy quark propagating in the Quark-Gluon Plasma (QGP) in the framework of the Moller theory, including possible large Coulomb logarithms as a perturbation to BDMPSZ bremsstrahlung, described in the Harmonic Oscillator (HO) approximation. We derive the analytical expression that describes the energy loss in the entire emitted gluon frequency region. In the small frequencies region, for angles larger than the dead cone angle, the energy loss is controlled by the BDMPSZ mechanism, while for larger frequencies it is described by N=1 term in the GLV opacity expansion. We estimate corresponding quenching rates for different values of the heavy quark path and different $m/E$ ratios.


title: Analysis of Lackadaisical Quantum Walks

abstract: The lackadaisical quantum walk is a quantum analogue of the lazy random walk obtained by adding a self-loop to each vertex in the graph. We analytically prove that lackadaisical quantum walks can find a unique marked vertex on any regular locally arc-transitive graph with constant success probability quadratically faster than the hitting time. This result proves several speculations and numerical findings in previous work, including the conjectures that the lackadaisical quantum walk finds a unique marked vertex with constant success probability on the torus, cycle, Johnson graphs, and other classes of vertex-transitive graphs. Our proof establishes and uses a relationship between lackadaisical quantum walks and quantum interpolated walks for any locally arc-transitive graph.


title: Lossless monochromation for electron microscopy with pulsed photoemission sources and rf cavities

abstract: Resonant radiofrequency cavities enable exquisite time-energy control of electron beams when synchronized with laser driven photoemission. We present a lossless monochromator design that exploits this fine control in the one-electron-per-pulse regime. The theoretically achievable maximum beam current on target is orders of magnitude greater than state-of-the-art monochromators for the same spatial and energy resolution. This improvement is the result of monochromating in the time domain, unconstrained by the transverse brightness of the electron source. We show analytically and confirm numerically that cavity parameters chosen to minimize energy spread perform the additional function of undoing the appreciable effect of chromatic aberration in the upstream optics. We argue that our design has significant applications in both ultra-fast and non-time-resolved microscopy, provided photoelectron sources of sufficiently small size and laser sources of sufficiently high repetition rate. Our design achieves in simulations more than two orders of magnitude reduction in beam energy spread, down to single digit meV. Overcoming the minimum probe-size limit that chromatic aberration imposes, our design clears a path for high-current, high-resolution electron beam applications at primary energies from single to hundreds of keV.


title: Paired Comparisons Modeling using t-Distribution with Bayesian Analysis

abstract: A paired comparison analysis is the simplest way to make comparative judgments between objects where objects may be goods, services or skills. For a set of problems, this technique helps to choose the most important problem to solve first and/or provides the solution that will be the most effective. This paper presents the theory of paired comparisons method and contributes to the paired comparisons models by developing a new model based on t-distribution. The developed model is illustrated using a data set of citations among four famous journals of Statistics. Using Bayesian analysis, the journals are ranked as JRSS-B --> Biometrika --> JASA --> Comm. in Stats.


title: Spectral Sparsification via Bounded-Independence Sampling

abstract: We give a deterministic, nearly logarithmic-space algorithm for mild spectral sparsification of undirected graphs. Given a weighted, undirected graph $G$ on $n$ vertices described by a binary string of length $N$, an integer $k\leq \log n$, and an error parameter $\epsilon > 0$, our algorithm runs in space $\tilde{O}(k\log (N\cdot w_{\mathrm{max}}/w_{\mathrm{min}}))$ where $w_{\mathrm{max}}$ and $w_{\mathrm{min}}$ are the maximum and minimum edge weights in $G$, and produces a weighted graph $H$ with $\tilde{O}(n^{1+2/k}/\epsilon^2)$ edges that spectrally approximates $G$, in the sense of Spielmen and Teng [ST04], up to an error of $\epsilon$. Our algorithm is based on a new bounded-independence analysis of Spielman and Srivastava's effective resistance based edge sampling algorithm [SS08] and uses results from recent work on space-bounded Laplacian solvers [MRSV17]. In particular, we demonstrate an inherent tradeoff (via upper and lower bounds) between the amount of (bounded) independence used in the edge sampling algorithm, denoted by $k$ above, and the resulting sparsity that can be achieved.


title: Graph Vertex Sampling with Arbitrary Graph Signal Hilbert Spaces

abstract: Graph vertex sampling set selection aims at selecting a set of ver-tices of a graph such that the space of graph signals that can be reconstructed exactly from those samples alone is maximal. In this context, we propose to extend sampling set selection based on spectral proxies to arbitrary Hilbert spaces of graph signals. Enabling arbitrary inner product of graph signals allows then to better account for vertex importance on the graph for a sampling adapted to the application. We first state how the change of inner product impacts sampling set selection and reconstruction, and then apply it in the context of geometric graphs to highlight how choosing an alternative inner product matrix can help sampling set selection and reconstruction.


title: Extremes of Censored and Uncensored Lifetimes in Survival Data

abstract: The i.i.d. censoring model for survival analysis assumes two independent sequences of i.i.d. positive random variables, $(T_i^*)_{1\le i\le n}$ and $(U_i)_{1\le i\le n}$. The data consists of observations on the random sequence $\big(T_i=\min(T_i^*,U_i)$ together with accompanying censor indicators. Values of $T_i$ with $T_i^*\le U_i$ are said to be uncensored, those with $T_i^*> U_i$ are censored. We assume that the distributions of the $T_i^*$ and $U_i$ are in the domain of attraction of the Gumbel distribution and obtain the asymptotic distributions, as sample size $n\to\infty$, of the maximum values of the censored and uncensored lifetimes in the data, and of statistics related to them. These enable us to examine questions concerning the possible existence of cured individuals in the population.


title: Gaussian unitary ensemble with jump discontinuities and the coupled Painlev\'e II and IV systems

abstract: We study the orthogonal polynomials and the Hankel determinants associated with Gaussian weight with two jump discontinuities. When the degree $n$ is finite, the orthogonal polynomials and the Hankel determinants are shown to be connected to the coupled Painlev\'e IV system. In the double scaling limit as the jump discontinuities tend to the edge of the spectrum and the degree $n$ grows to infinity, we establish the asymptotic expansions for the Hankel determinants and the orthogonal polynomials, which are expressed in terms of solutions of the coupled Painlev\'{e} II system. As applications, we re-derive the recently found Tracy-Widom type expressions for the gap probability of there being no eigenvalues in a finite interval near the the extreme eigenvalue of large Gaussian unitary ensemble and the limiting conditional distribution of the largest eigenvalue in Gaussian unitary ensemble by considering a thinned process.


title: Lightweight Online Separation of the Sound Source of Interest through BLSTM-Based Binary Masking

abstract: Online audio source separation has been an important part of auditory scene analysis and robot audition. The main type of technique to carry this out, because of its online capabilities, has been spatial filtering (or beamforming), where it is assumed that the location (mainly, the direction of arrival; DOA) of the source of interest (SOI) is known. However, these techniques suffer from considerable interference leakage in the final result. In this paper, we propose a two step technique: 1) a phase-based beamformer that provides, in addition to the estimation of the SOI, an estimation of the cumulative environmental interference; and 2) a BLSTM-based TF binary masking stage that calculates a binary mask that aims to separate the SOI from the cumulative environmental interference. In our tests, this technique provides a signal-to-interference ratio (SIR) above 20 dB with simulated data. Because of the nature of the beamformer outputs, the label permutation problem is handled from the beginning. This makes the proposed solution a lightweight alternative that requires considerably less computational resources (almost an order of magnitude) compared to current deep-learning based techniques, while providing a comparable SIR performance.


title: Attacks Which Do Not Kill Training Make Adversarial Learning Stronger

abstract: Adversarial training based on the minimax formulation is necessary for obtaining adversarial robustness of trained models. However, it is conservative or even pessimistic so that it sometimes hurts the natural generalization. In this paper, we raise a fundamental question---do we have to trade off natural generalization for adversarial robustness? We argue that adversarial training is to employ confident adversarial data for updating the current model. We propose a novel approach of friendly adversarial training (FAT): rather than employing most adversarial data maximizing the loss, we search for least adversarial (i.e., friendly adversarial) data minimizing the loss, among the adversarial data that are confidently misclassified. Our novel formulation is easy to implement by just stopping the most adversarial data searching algorithms such as PGD (projected gradient descent) early, which we call early-stopped PGD. Theoretically, FAT is justified by an upper bound of the adversarial risk. Empirically, early-stopped PGD allows us to answer the earlier question negatively---adversarial robustness can indeed be achieved without compromising the natural generalization.


title: Correspondence Analysis between the Location and the Leading Causes of Death in the United States

abstract: Correspondence Analysis analyzes two-way or multi-way tables withe each row and column becoming a point ion a multidimensional graphical map called biplot. It can be used to extract essential dimensions allowing simplification of the data matrix. This study aims to measure the association between the location and the leading causes of death in the United States of America and to determine the location where a particular disease is highly associated. The research data consists of two variables with 510 data points. Results show that there is a significant association between the location ad leading cause of death in the United States, and 61% of the variance in the model are explained by the first two dimensions.


title: Transfer Learning from Synthetic to Real-Noise Denoising with Adaptive Instance Normalization

abstract: Real-noise denoising is a challenging task because the statistics of real-noise do not follow the normal distribution, and they are also spatially and temporally changing. In order to cope with various and complex real-noise, we propose a well-generalized denoising architecture and a transfer learning scheme. Specifically, we adopt an adaptive instance normalization to build a denoiser, which can regularize the feature map and prevent the network from overfitting to the training set. We also introduce a transfer learning scheme that transfers knowledge learned from synthetic-noise data to the real-noise denoiser. From the proposed transfer learning, the synthetic-noise denoiser can learn general features from various synthetic-noise data, and the real-noise denoiser can learn the real-noise characteristics from real data. From the experiments, we find that the proposed denoising method has great generalization ability, such that our network trained with synthetic-noise achieves the best performance for Darmstadt Noise Dataset (DND) among the methods from published papers. We can also see that the proposed transfer learning scheme robustly works for real-noise images through the learning with a very small number of labeled data.


title: Linear-Time Variational Integrators in Maximal Coordinates

abstract: Most dynamic simulation tools parameterize the configuration of multi-body robotic systems using minimal coordinates, also called generalized or joint coordinates. However, maximal-coordinate approaches have several advantages over minimal-coordinate parameterizations, including native handling of closed kinematic loops and nonholonomic constraints. This paper describes a linear-time variational integrator that is formulated in maximal coordinates. Due to its variational formulation, the algorithm does not suffer from constraint drift and has favorable energy and momentum conservation properties. A sparse matrix factorization technique allows the dynamics of a loop-free articulated mechanism with $n$ links to be computed in $O(n)$ (linear) time. Additional constraints that introduce loops can also be handled by the algorithm without incurring much computational overhead. Experimental results show that our approach offers speed competitive with state-of-the-art minimal-coordinate algorithms while outperforming them in several scenarios, especially when dealing with closed loops and configuration singularities.


title: Supervised Categorical Metric Learning with Schatten p-Norms

abstract: Metric learning has been successful in learning new metrics adapted to numerical datasets. However, its development on categorical data still needs further exploration. In this paper, we propose a method, called CPML for \emph{categorical projected metric learning}, that tries to efficiently~(i.e. less computational time and better prediction accuracy) address the problem of metric learning in categorical data. We make use of the Value Distance Metric to represent our data and propose new distances based on this representation. We then show how to efficiently learn new metrics. We also generalize several previous regularizers through the Schatten $p$-norm and provides a generalization bound for it that complements the standard generalization bound for metric learning. Experimental results show that our method provides


title: Collision Avoidance Based on Robust Lexicographic Task Assignment

abstract: Traditional task assignment approaches for multi-agent motion control do not take the possibility of collisions into account. This can lead to challenging requirements for path planning. We derive an assignment method that not only minimises the largest distance between an agent and its assigned destination but also provides local constraints for guaranteed collision avoidance. To this end, we introduce a sequential bottleneck optimisation problem and define a notion of robustness of an optimising assignment to changes of individual assignment costs. Conditioned on a sufficient level of robustness in relation to the size of the agents, we construct time-varying position bounds for every individual agent. These local constraints are a direct byproduct of the assignment procedure and only depend on the initial agent positions, the destinations that are to be visited, and a timing parameter. We prove that no agent that is assigned to move to one of the target locations collides with any other agent if all agents satisfy their local position constraints. We demonstrate the method in a illustrative case study.


title: Super-Resolving Commercial Satellite Imagery Using Realistic Training Data

abstract: In machine learning based single image super-resolution, the degradation model is embedded in training data generation. However, most existing satellite image super-resolution methods use a simple down-sampling model with a fixed kernel to create training images. These methods work fine on synthetic data, but do not perform well on real satellite images. We propose a realistic training data generation model for commercial satellite imagery products, which includes not only the imaging process on satellites but also the post-process on the ground. We also propose a convolutional neural network optimized for satellite images. Experiments show that the proposed training data generation model is able to improve super-resolution performance on real satellite images.


title: Concatenated Convolution-Polar codes over Rayleigh Channels degraded as Erasure Channels

abstract: Polar codes have been shown to approach capacity of symmetric binary erasure channels and also have low encoding and decoding complexity. Wireless channels are bursty in nature and a method to apply polar codes over wireless channels is proposed. This paper investigates the performance of polar codes over typical urban fading wireless channels. The wireless channel is modeled as a degraded erasure channel by using an inner convolution code and interleaving. When the convolution code fails to correct any errors, the corresponding bits are discarded and treated as erasures. We present the performance of polar codes over such degraded wireless channels using different code rates for the inner convolution code.


title: Dataset of raw and pre-processed speech signals, Mel Frequency Cepstral Coefficients of Speech and Heart Rate measurements

abstract: Heart rate is an important vital sign used in the diagnosis of many medical conditions. Conventionally, heart rate is measured using a medical device such as pulse oxymeter. Physiological parameters such as heart rate bear a correlation to speech characteristics of an individual. Hence, there is a possibility to measure heart rate from speech signals using machine learning and deep learning, which would also allow non-invasive, non contact based and remote monitoring of patients. However, to design such a scheme and verify its accuracy, it is necessary to collect speech recordings along with heart rates measured using a medical device, simultaneously during the recording. This article provides a dataset as well as the procedure used to create the dataset which could be used to facilitate research in developing techniques to estimate heart rate accurately by observing speech signal.


title: Back to the Future: Joint Aware Temporal Deep Learning 3D Human Pose Estimation

abstract: We propose a new deep learning network that introduces a deeper CNN channel filter and constraints as losses to reduce joint position and motion errors for 3D video human body pose estimation. Our model outperforms the previous best result from the literature based on mean per-joint position error, velocity error, and acceleration errors on the Human 3.6M benchmark corresponding to a new state-of-the-art mean error reduction in all protocols and motion metrics. Mean per joint error is reduced by 1%, velocity error by 7% and acceleration by 13% compared to the best results from the literature. Our contribution increasing positional accuracy and motion smoothness in video can be integrated with future end to end networks without increasing network complexity. Our model and code are available at https://vnmr.github.io/ Keywords: 3D, human, image, pose, action, detection, object, video, visual, supervised, joint, kinematic


title: AutoEmb: Automated Embedding Dimensionality Search in Streaming Recommendations

abstract: Deep learning based recommender systems (DLRSs) often have embedding layers, which are utilized to lessen the dimensionality of categorical variables (e.g. user/item identifiers) and meaningfully transform them in the low-dimensional space. The majority of existing DLRSs empirically pre-define a fixed and unified dimension for all user/item embeddings. It is evident from recent researches that different embedding sizes are highly desired for different users/items according to their popularity. However, manually selecting embedding sizes in recommender systems can be very challenging due to the large number of users/items and the dynamic nature of their popularity. Thus, in this paper, we propose an AutoML based end-to-end framework (AutoEmb), which can enable various embedding dimensions according to the popularity in an automated and dynamic manner. To be specific, we first enhance a typical DLRS to allow various embedding dimensions; then we propose an end-to-end differentiable framework that can automatically select different embedding dimensions according to user/item popularity; finally we propose an AutoML based optimization algorithm in a streaming recommendation setting. The experimental results based on widely used benchmark datasets demonstrate the effectiveness of the AutoEmb framework.


title: Reconstructing particle number distributions with convoluting volume fluctuations

abstract: We propose methods to reconstruct particle distributions with and without considering initial volume fluctuations. This approach enables us to correct for detector efficiencies and initial volume fluctuations simultaneously. Our study suggests such a tool could investigate the possible bimodal structure of net-proton distribution in Au+Au collisions at $\sqrt{s_{\rm NN}}=$7.7 GeV a signature of first-order phase transition and critical point [arXiv:1804.04463,arXiv:1811.04456].


title: Star order automorphisms on the poset of type 1 operators

abstract: Let $\mathcal{H}$ be a complex infinite dimensional Hilbert space and $\mathcal{B}(\mathcal{H})$ the algebra of all bounded linear operators on $\mathcal H$. The star partial order is defined by $A\overset{*}{\leq}B$ if and only if $A^*A=A^*B$ and $AA^*=AB^*$ for any $A$ and $B$ in $\mathcal B(\mathcal H)$. We give a type decomposition of operators with respect to star order. For any $A\in\mathcal B(\mathcal H)$, there are unique type 1 operator $A_1 $ which is $0$ or the supremum of those rank 1 operators less than $A_1$ and type 2 operator $A_2$ which is not greater than any rank 1 operator in star order such that $A_i\overset{*}{\leq}A$$(i=1,2)$ and $A=A_1+A_2$. Moreover, we determine all automorphisms on the poset of type 1 operators. As a consequence, we characterize continuous automorphisms on $\mathcal B(\mathcal H)$.


title: An Optimal Statistical and Computational Framework for Generalized Tensor Estimation

abstract: This paper describes a flexible framework for generalized low-rank tensor estimation problems that includes many important instances arising from applications in computational imaging, genomics, and network analysis. The proposed estimator consists of finding a low-rank tensor fit to the data under generalized parametric models. To overcome the difficulty of non-convexity in these problems, we introduce a unified approach of projected gradient descent that adapts to the underlying low-rank structure. Under mild conditions on the loss function, we establish both an upper bound on statistical error and the linear rate of computational convergence through a general deterministic analysis. Then we further consider a suite of generalized tensor estimation problems, including sub-Gaussian tensor PCA, tensor regression, and Poisson and binomial tensor PCA. We prove that the proposed algorithm achieves the minimax optimal rate of convergence in estimation error. Finally, we demonstrate the superiority of the proposed framework via extensive experiments on both simulated and real data.


title: Incorporating Expert Prior Knowledge into Experimental Design via Posterior Sampling

abstract: Scientific experiments are usually expensive due to complex experimental preparation and processing. Experimental design is therefore involved with the task of finding the optimal experimental input that results in the desirable output by using as few experiments as possible. Experimenters can often acquire the knowledge about the location of the global optimum. However, they do not know how to exploit this knowledge to accelerate experimental design. In this paper, we adopt the technique of Bayesian optimization for experimental design since Bayesian optimization has established itself as an efficient tool for optimizing expensive black-box functions. Again, it is unknown how to incorporate the expert prior knowledge about the global optimum into Bayesian optimization process. To address it, we represent the expert knowledge about the global optimum via placing a prior distribution on it and we then derive its posterior distribution. An efficient Bayesian optimization approach has been proposed via posterior sampling on the posterior distribution of the global optimum. We theoretically analyze the convergence of the proposed algorithm and discuss the robustness of incorporating expert prior. We evaluate the efficiency of our algorithm by optimizing synthetic functions and tuning hyperparameters of classifiers along with a real-world experiment on the synthesis of short polymer fiber. The results clearly demonstrate the advantages of our proposed method.


title: Photoinduced $\eta$ pairing in the Kondo lattice model

abstract: The previous theoretical study has shown that pulse irradiation to the Mott insulating state in the Hubbard model can induce the enhancement of superconducting correlation due to the generation of $\eta$ pairs. Here, we show that the same mechanism can be applied to the Kondo lattice model, an effective model for heavy electron systems, by demonstrating that the pulse irradiation indeed enhances the $\eta$-pairing correlation. As in the case of the Hubbard model, the non-linear optical process is essential to increase the number of photoinduced $\eta$ pairs and thus the enhancement of the superconducting correlation. We also find the diffusive behavior of the spin dynamics after the pulse irradiation, suggesting that the increase of the number of $\eta$ pairs leads to the decoupling between the conduction band and the localized spins in the Kondo lattice model, which is inseparably related to the photodoping effect.


title: Fast Lower and Upper Estimates for the Price of Constrained Multiple Exercise American Options by Single Pass Lookahead Search and Nearest-Neighbor Martingale

abstract: This article presents fast lower and upper estimates for a large class of options: the class of constrained multiple exercise American options. Typical options in this class are swing options with volume and timing constraints, and passport options with multiple lookback rights. The lower estimate algorithm uses the artificial intelligence method of lookahead search. The upper estimate algorithm uses the dual approach to option pricing on a nearest-neighbor basis for the martingale space. Probabilistic convergence guarantees are provided. Several numerical examples illustrate the approaches including a swing option with four constraints, and a passport option with 16 constraints.


title: Dimensional Analysis in Statistical Modelling

abstract: Building on recent work in statistical science, the paper presents a theory for modelling natural phenomena that unifies physical and statistical paradigms based on the underlying principle that a model must be nondimensionalizable. After all, such phenomena cannot depend on how the experimenter chooses to assess them. Yet the model itself must be comprised of quantities that can be determined theoretically or empirically. Hence, the underlying principle requires that the model represents these natural processes correctly no matter what scales and units of measurement are selected. This goal was realized for physical modelling through the celebrated theories of Buckingham and Bridgman and for statistical modellers through the invariance principle of Hunt and Stein. Building on recent research in statistical science, the paper shows how the latter can embrace and extend the former. The invariance principle is extended to encompass the Bayesian paradigm, thereby enabling an assessment of model uncertainty. The paper covers topics not ordinarily seen in statistical science regarding dimensions, scales, and units of quantities in statistical modelling. It shows the special difficulties that can arise when models involve transcendental functions, such as the logarithm which is used e.g. in likelihood analysis and is a singularity in the family of Box-Cox family of transformations. Further, it demonstrates the importance of the scale of measurement, in particular how differently modellers must handle ratio- and interval-scales


title: Heat flow reversal in a trapped-ion simulator

abstract: We propose a trapped-ion platform to simulate a reconfigurable spin-spin Hamiltonian related to quantum thermodynamic processes. Starting from an experimental model describing two trapped-ions under slightly off-resonant first sideband driving with individually controlled driving phases, we follow an operational quantum optics approach to show that it produces an effective model appearing in recent quantum thermodynamics proposals. We show that projection into the vibrational ground state manifold allows full analytic treatment. As a practical example, we take experimental data from a $^{187}$Yb$^{+}$ trap and numerically simulate the reversal of heat flow between two thermal spins controlled by their quantum correlations.


title: Multi-Attribute Guided Painting Generation

abstract: Controllable painting generation plays a pivotal role in image stylization. Currently, the control way of style transfer is subject to exemplar-based reference or a random one-hot vector guidance. Few works focus on decoupling the intrinsic properties of painting as control conditions, e.g., artist, genre and period. Under this circumstance, we propose a novel framework adopting multiple attributes from the painting to control the stylized results. An asymmetrical cycle structure is equipped to preserve the fidelity, associating with style preserving and attribute regression loss to keep the unique distinction of colors and textures between domains. Several qualitative and quantitative results demonstrate the effect of the combinations of multiple attributes and achieve satisfactory performance.


title: DLSpec: A Deep Learning Task Exchange Specification

abstract: Deep Learning (DL) innovations are being introduced at a rapid pace. However, the current lack of standard specification of DL tasks makes sharing, running, reproducing, and comparing these innovations difficult. To address this problem, we propose DLSpec, a model-, dataset-, software-, and hardware-agnostic DL specification that captures the different aspects of DL tasks. DLSpec has been tested by specifying and running hundreds of DL tasks.


title: Learning Light Field Angular Super-Resolution via a Geometry-Aware Network

abstract: The acquisition of light field images with high angular resolution is costly. Although many methods have been proposed to improve the angular resolution of a sparsely-sampled light field, they always focus on the light field with a small baseline, which is captured by a consumer light field camera. By making full use of the intrinsic \textit{geometry} information of light fields, in this paper we propose an end-to-end learning-based approach aiming at angularly super-resolving a sparsely-sampled light field with a large baseline. Our model consists of two learnable modules and a physically-based module. Specifically, it includes a depth estimation module for explicitly modeling the scene geometry, a physically-based warping for novel views synthesis, and a light field blending module specifically designed for light field reconstruction. Moreover, we introduce a novel loss function to promote the preservation of the light field parallax structure. Experimental results over various light field datasets including large baseline light field images demonstrate the significant superiority of our method when compared with state-of-the-art ones, i.e., our method improves the PSNR of the second best method up to 2 dB in average, while saves the execution time 48$\times$. In addition, our method preserves the light field parallax structure better.


title: Update on the IsoDAR 60 MeV/amu Cyclotron Acceleration Simulations

abstract: In this technical report, the work done by members of the DAE$\delta$ALUS collaboration on simulating the acceleration of 5 mA of H2+ in the IsoDAR compact isochronous cyclotron is summarized. The findings are that sufficient turn separation between the (N-1)th and Nth turn can be achieved with either initially matched but also with mismatched beams by careful placement of collimators to scrape away halo particles before the beam energy has reached 2 MeV/amu. Here we define "sufficient" as < 200 W beam loss on the septum. In the acceleration simulations, beam losses on the collimators vary from 10% to 25% depending on the initial beam conditions. The collimator placement in the central region, halo formation, and turn separation at the final turn are discussed. Finally, a scheme to further mitigate the risk of beam loss on the extraction septum by placing a narrow stripping foil in front of it is investigated. We also present the findings of a separate study subcontracted to the company AIMA, who investigated injection of H2+ ions through a spiral inflector and accelerating them for four turns. Here the reported beam loss before 2 MeV/amu was 58%, but bunches were not accelerated further, so we are not reporting this as a final result.


title: Exceptional band touching for strongly correlated systems in equilibrium

abstract: Quasi-particles described by Green's functions of equilibrium systems exhibit non-Hermitian topological phenomena because of their finite lifetime. This non-Hermitian perspective on equilibrium systems provides new insights into correlated systems and attracts much interest because of its potential to solve open questions in correlated compounds. In this paper, we provide a concise review of the non-Hermitian topological band structures for quantum many-body systems in equilibrium as well as their classification.

categories: cond-mat.str-el cond-mat.dis-nn cond-mat.mes-hall quant-ph


title: Wide-Sense 2-Frameproof Codes

abstract: Various kinds of fingerprinting codes and their related combinatorial structures are extensively studied for protecting copyrighted materials. This paper concentrates on one specialised fingerprinting code named wide-sense frameproof codes in order to prevent innocent users from being framed. Let $Q$ be a finite alphabet of size $q$. Given a $t$-subset $X=\{x ^1,\ldots, x ^t\}\subseteq Q^n$, a position $i$ is called undetectable for $X$ if the values of the words of $X$ match in their $i$th position: $x_i^1=\cdots=x_i^t$. The wide-sense descendant set of $X$ is defined by $\wdesc(X)=\{y\in Q^n:y_i=x_i^1,i\in {U}(X)\},$ where ${U}(X)$ is the set of undetectable positions for $X$. A code ${\cal C}\subseteq Q^n$ is called a wide-sense $t$-frameproof code if $\wdesc(X) \cap{\cal C} = X$ for all $X \subseteq {\cal C}$ with $|X| \le t$. The paper improves the upper bounds on the sizes of wide-sense $2$-frameproof codes by applying techniques on non $2$-covering Sperner families and intersecting families in extremal set theory.


title: Ross-Macdonald Models: Which one should we use?

abstract: Ross-Macdonald models are the building blocks of most vector-borne disease models. Even for the same disease, different authors use different model formulations, but a study of the dynamical consequences of assuming different hypotheses is missing. In this work we present different formulations of the basic Ross-Macdonald model together with a careful discussion of the assumptions behind each model. The most general model presented is an agent based model for which arbitrary distributions for latency and infectious periods for both, host and vectors, is considered. At population level we also developed a deterministic Volterra integral equations model for which also arbitrary distributions in the waiting times are included. We compare the model solutions using different distributions for the infectious and latency periods using statistics, like the epidemic peak, or epidemic final size, to characterize the epidemic curves. The basic reproduction number ($R_0$) for each formulation is computed and compared with empirical estimations obtained with the agent based models. The importance of considering realistic distributions for the latent and infectious periods is highlighted and discussed. We also show that seasonality is a key driver of vector-borne disease dynamics shaping the epidemic curve and its duration.


title: A Density Ratio Approach to Language Model Fusion in End-To-End Automatic Speech Recognition

abstract: This article describes a density ratio approach to integrating external Language Models (LMs) into end-to-end models for Automatic Speech Recognition (ASR). Applied to a Recurrent Neural Network Transducer (RNN-T) ASR model trained on a given domain, a matched in-domain RNN-LM, and a target domain RNN-LM, the proposed method uses Bayes' Rule to define RNN-T posteriors for the target domain, in a manner directly analogous to the classic hybrid model for ASR based on Deep Neural Networks (DNNs) or LSTMs in the Hidden Markov Model (HMM) framework (Bourlard & Morgan, 1994). The proposed approach is evaluated in cross-domain and limited-data scenarios, for which a significant amount of target domain text data is used for LM training, but only limited (or no) {audio, transcript} training data pairs are used to train the RNN-T. Specifically, an RNN-T model trained on paired audio & transcript data from YouTube is evaluated for its ability to generalize to Voice Search data. The Density Ratio method was found to consistently outperform the dominant approach to LM and end-to-end ASR integration, Shallow Fusion.


title: Is the OWASP Top 10 list comprehensive enough for writing secure code?

abstract: The OWASP Top 10 is a list that is published by the Open Web Application Security Project (OWASP). The general purpose is to serve as a watchlist for bugs to avoid while writing code. This paper compares how many of those weakness as described in the top ten list are actually reported in vulnerabilities listed in the National Vulnerability Database (NVD). That way it makes it possible to empirically show whether the OWASP Top 10 list is comprehensive enough or not, for code weaknesses that have been found in the past decade.


title: DNN-Chip Predictor: An Analytical Performance Predictor for DNN Accelerators with Various Dataflows and Hardware Architectures

abstract: The recent breakthroughs in deep neural networks (DNNs) have spurred a tremendously increased demand for DNN accelerators. However, designing DNN accelerators is non-trivial as it often takes months/years and requires cross-disciplinary knowledge. To enable fast and effective DNN accelerator development, we propose DNN-Chip Predictor, an analytical performance predictor which can accurately predict DNN accelerators' energy, throughput, and latency prior to their actual implementation. Our Predictor features two highlights: (1) its analytical performance formulation of DNN ASIC/FPGA accelerators facilitates fast design space exploration and optimization; and (2) it supports DNN accelerators with different algorithm-to-hardware mapping methods (i.e., dataflows) and hardware architectures. Experiment results based on 2 DNN models and 3 different ASIC/FPGA implementations show that our DNN-Chip Predictor's predicted performance differs from those of chip measurements of FPGA/ASIC implementation by no more than 17.66% when using different DNN models, hardware architectures, and dataflows. We will release code upon acceptance.


title: Comparing copy-number profiles under multi-copy amplifications and deletions

abstract: During cancer progression, malignant cells accumulate somatic mutations that can lead to genetic aberrations. In particular, evolutionary events akin to segmental duplications or deletions can alter the copy-number profile (CNP) of a set of genes in a genome. Our aim is to compute the evolutionary distance between two cells for which only CNPs are known. This asks for the minimum number of segmental amplifications and deletions to turn one CNP into another. This was recently formalized into a model where each event is assumed to alter a copy-number by $1$ or $-1$, even though these events can affect large portions of a chromosome. We propose a general cost framework where an event can modify the copy-number of a gene by larger amounts. We show that any cost scheme that allows segmental deletions of arbitrary length makes computing the distance strongly NP-hard. We then devise a factor $2$ approximation algorithm for the problem when copy-numbers are non-zero and provide an implementation called \textsf{cnp2cnp}. We evaluate our approach experimentally by reconstructing simulated cancer phylogenies from the pairwise distances inferred by \textsf{cnp2cnp} and compare it against two other alternatives, namely the \textsf{MEDICC} distance and the Euclidean distance. The experimental results show that our distance yields more accurate phylogenies on average than these alternatives if the given CNPs are error-free, but that the \textsf{MEDICC} distance is slightly more robust against error in the data. In all cases, our experiments show that either our approach or the \textsf{MEDICC} approach should preferred over the Euclidean distance.


title: A class of three-dimensional gyroviscous magnetohydrodynamic models

abstract: A Hamiltonian and action principle formalism for deriving three-dimensional gyroviscous magnetohydrodynamic models is presented. The uniqueness of the approach in constructing the gyroviscous tensor from first principles and its ability to explain the origin of the gyromap and the gyroviscous terms are highlighted. The procedure allows for the specification of free functions, which can be used to generate a wide range of gyroviscous models. Through the process of reduction, the noncanonical Hamiltonian bracket is obtained and briefly analysed.


title: A Systematic Survey of General Sparse Matrix-Matrix Multiplication

abstract: General Sparse Matrix-Matrix Multiplication (SpGEMM) has attracted much attention from researchers in graph analyzing, scientific computing, and deep learning. Many optimization techniques have been developed for different applications and computing architectures over the past decades. The objective of this paper is to provide a structured and comprehensive overview of the researches on SpGEMM. Existing researches have been grouped into different categories based on target architectures and design choices. Covered topics include typical applications, compression formats, general formulations, key problems and techniques, architecture-oriented optimizations, and programming models. The rationales of different algorithms are analyzed and summarized. This survey sufficiently reveals the latest progress of SpGEMM research to 2021. Moreover, a thorough performance comparison of existing implementations is presented. Based on our findings, we highlight future research directions, which encourage better design and implementations in later studies.


title: Structures, stellar population properties, AGN fractions, and environments of massive compact galaxies at $1 < z < 2$ in 3D--{\it HST}/CANDELS

abstract: We present a study on structures and physical properties of massive ($M_* >10^{10} M_{\sun} $) compact galaxies at $1.0<z<2.0$ in five 3D--{\it HST}/CANDELS fields. Compared with the extended star-forming galaxies (eSFGs), compact star-forming galaxies (cSFGs) are found to have the lower level of star formation, and mainly distribute in the quiescent region of the {\it UVJ} diagram. The distributions of dust attenuation and S{\'e}rsic index support that the progenitors of cQGs are cSFGs, and cSFGs are at a transitional phase between eSFGs and cQGs. The prevalence of X-ray selected AGNs ($\sim 28\%$) is confirmed in the cSFGs at $1<z<2$ which indicates that the violent gas-rich processes such as merger and disk instability could drive the structure to be more compact, and trigger both star formation and black hole growth in the central regions. Our results support the "two-step" scenario that the cSFGs at $1<z<2$ are the intermediate population after compaction but before a quick quenching. Our analysis of parametric and nonparametric morphologies shows that cQGs (eQGs) are more concentrated and have less substructures than cSFGs (eSFGs), and quenching and compactness should be associated with each other. The cSFGs at $1.5<z<2$ ($1<z<1.5$) prefer to be in higher (lower) density environment, similar as cQGs (eSFGs). It suggests that merger or strong interaction might be the main driving mechanism of compaction at higher redshifts, whereas the disk instability of individual galaxies might play a more important role on the formation of cSFGs at lower redshifts.


title: Adversarial Monte Carlo Meta-Learning of Optimal Prediction Procedures

abstract: We frame the meta-learning of prediction procedures as a search for an optimal strategy in a two-player game. In this game, Nature selects a prior over distributions that generate labeled data consisting of features and an associated outcome, and the Predictor observes data sampled from a distribution drawn from this prior. The Predictor's objective is to learn a function that maps from a new feature to an estimate of the associated outcome. We establish that, under reasonable conditions, the Predictor has an optimal strategy that is equivariant to shifts and rescalings of the outcome and is invariant to permutations of the observations and to shifts, rescalings, and permutations of the features. We introduce a neural network architecture that satisfies these properties. The proposed strategy performs favorably compared to standard practice in both parametric and nonparametric experiments.


title: A Balancing Weight Framework for Estimating the Causal Effect of General Treatments

abstract: In observational studies, weighting methods that directly optimize the balance between treatment and covariates have received much attention lately; however these have mainly focused on binary treatments. Inspired by domain adaptation, we show that such methods can be actually reformulated as specific implementations of a discrepancy minimization problem aimed at tackling a shift of distribution from observational to interventional data. More precisely, we introduce a new framework, Covariate Balance via Discrepancy Minimization (CBDM), that provably encompasses most of the existing balancing weight methods and formally extends them to treatments of arbitrary types (e.g., continuous or multivariate). We establish theoretical guarantees for our framework that both offer generalizations of properties known when the treatment is binary, and give a better grasp on what hyperparameters to choose in non-binary settings. Based on such insights, we propose a particular implementation of CBDM for estimating dose-response curves and demonstrate through experiments its competitive performance relative to other existing approaches for continuous treatments.


title: Learning Product Graphs Underlying Smooth Graph Signals

abstract: Real-world data is often times associated with irregular structures that can analytically be represented as graphs. Having access to this graph, which is sometimes trivially evident from domain knowledge, provides a better representation of the data and facilitates various information processing tasks. However, in cases where the underlying graph is unavailable, it needs to be learned from the data itself for data representation, data processing and inference purposes. Existing literature on learning graphs from data has mostly considered arbitrary graphs, whereas the graphs generating real-world data tend to have additional structure that can be incorporated in the graph learning procedure. Structure-aware graph learning methods require learning fewer parameters and have the potential to reduce computational, memory and sample complexities. In light of this, the focus of this paper is to devise a method to learn structured graphs from data that are given in the form of product graphs. Product graphs arise naturally in many real-world datasets and provide an efficient and compact representation of large-scale graphs through several smaller factor graphs. To this end, first the graph learning problem is posed as a linear program, which (on average) outperforms the state-of-the-art graph learning algorithms. This formulation is of independent interest itself as it shows that graph learning is possible through a simple linear program. Afterwards, an alternating minimization-based algorithm aimed at learning various types of product graphs is proposed, and local convergence guarantees to the true solution are established for this algorithm. Finally the performance gains, reduced sample complexity, and inference capabilities of the proposed algorithm over existing methods are also validated through numerical simulations on synthetic and real datasets.


title: Merging Black Holes in the Low-mass and High-mass Gaps from 2+2 Quadruple Systems

abstract: The origin of the black hole (BH) binary mergers observed by LIGO-Virgo is still uncertain, as are the boundaries of the stellar BH mass function. Stellar evolution models predict a dearth of BHs both at masses $\gtrsim 50$ M$_\odot$ and $\lesssim 5$ M$_\odot$, thus leaving low- and high-mass gaps in the BH mass function. A natural way to form BHs of these masses is through mergers of neutron stars (NSs; for the low-mass gap) or lower-mass BHs (for the high-mass gap); the low- or high-mass-gap BH produced as a merger product can then be detected by LIGO-Virgo if it merges again with a new companion. We show that the evolution of a 2+2 quadruple system can naturally lead to BH mergers with component masses in the low- or high-mass gaps. In our scenario, the BH in the mass gap originates from the merger of two NSs, or two BHs, in one of the two binaries and the merger product is imparted a recoil kick (from anisotropic gravitational wave emission), which triggers its interaction with the other binary component of the quadruple system. The outcome of this three-body interaction is usually a new eccentric compact binary containing the BH in the mass gap, which can then merge again. The merger rate is $\sim 10^{-7} - 10^{-2}$ Gpc$^{-3}$ yr$^{-1}$ and $\sim 10^{-3} - 10^{-2}$ Gpc$^{-3}$ yr$^{-1}$ for BHs in the low-mass and high-mass gap, respectively. As the sensitivity of gravitational wave detectors improves, tighter constraints will soon be placed on the stellar BH mass function.


title: $D \to P(\pi,K)$ helicity form factors within light-cone sum rule approach

abstract: In this paper, the $D\to P(\pi, K)$ helicity form factors (HFFs) are studied by applying the QCD light-cone sum rule (LCSR) approach. The calculation accuracy is up to next-to-leading order (NLO) gluon radiation correction of twist-(2,3) distribution amplitude. The resultant HFFs at large recoil point are ${\cal P}_{t,0}^\pi(0) = 0.688^{+0.020}_{-0.024}$, ${\cal P}_{t,0}^K(0)=0.780^{+0.024}_{-0.029}$. In which, the contributions from three particles of the leading order (LO) are so small that can be safely neglected, and the maximal contribution of the NLO gluon radiation correction for ${\cal P}_{t,0}^{\pi,K}(0)$ is less than $3\%$. After extrapolating the LCSR predictions for these HFFs to whole $q^2$-region, we obtain the decay widths for semileptonic decay processes $D\to P\ell\nu_\ell$, which are consistent with BES-III collaboration predictions within errors. After considering the $D^{+}/D^{0}$-meson lifetime, we give the branching fractions of $D\to P\ell\nu_\ell$ with $\ell = e, \mu$, our predictions also agree with BES-III collaboration within errors, especially for $D\to \pi \ell\nu_\ell$ decay process. Finally, we present the forward-backward asymmetry ${\cal A}_{\rm FB}^\ell(q^2)$ and lepton convexity parameter ${\cal C}_F^\ell(q^2)$, and further calculate the mean value of these two observations $\langle{\cal A}_{\rm FB}^\ell\rangle$ and $\langle{\cal C}_F^\ell\rangle$, which may provide a way to test those HFFs in future experiments.


title: Matem\'aticas B\'asicas (con Software Libre)

abstract: This is a book on elementary mathematics suitable to use in courses for first-year students of science and engineering. Students often find the formal approach to math courses rather discouraging, and demand 'rel-life applications' in order to grasp the need for these courses. The text has plenty of such applications, treated with the free software Maxima (a CAS) and Geogebra (a DGS). Its level is more or less that of a 'Pre-Calculus' course, but the contents and style of presentation are different in that they are determined by the applications in each chapter.


title: Generalized Product Quantization Network for Semi-supervised Image Retrieval

abstract: Image retrieval methods that employ hashing or vector quantization have achieved great success by taking advantage of deep learning. However, these approaches do not meet expectations unless expensive label information is sufficient. To resolve this issue, we propose the first quantization-based semi-supervised image retrieval scheme: Generalized Product Quantization (GPQ) network. We design a novel metric learning strategy that preserves semantic similarity between labeled data, and employ entropy regularization term to fully exploit inherent potentials of unlabeled data. Our solution increases the generalization capacity of the quantization network, which allows overcoming previous limitations in the retrieval community. Extensive experimental results demonstrate that GPQ yields state-of-the-art performance on large-scale real image benchmark datasets.


title: Transient Convective Spin-up Dynamics

abstract: We study the formation, longevity and breakdown of convective rings during impulsive spin-up in square and cylindrical containers using direct numerical simulations. The rings, which are axisymmetric alternating regions of up- and down-welling flow that can last for O (100) rotation times, were first demonstrated experimentally and arise due to a balance between Coriolis and viscous effects. We study the formation of these rings in the context of the Greenspan-Howard spin-up process, the disruption of which modifies ring formation and evolution. We show that, unless imprinted by boundary geometry, convective rings can only form when the surface providing buoyancy forcing is a free-slip surface, thereby explaining an apparent disagreement between experimental results in the literature. For Prandtl numbers from 1--5 we find that the longest-lived rings occur for intermediate Prandtl numbers, with a Rossby number dependence. Finally, we find that the constant evaporative heat-flux conditions imposed in the experiments are essential in sustaining the rings and in maintaining the vortices that form in consequence of the ring breakdown.


title: Age-optimal Service and Decision Scheduling in Internet of Things

abstract: We consider an Internet of Things (IoT) system in which a sensor observes a phenomena of interest with exponentially distributed intervals and delivers the updates to a monitor with the First-come-First-served (FCFS) policy. At the monitor, the received updates are used to make decisions with deterministic or random intervals. For this system, we investigate the freshness of the updates at these decision epochs using the age upon decisions (AuD) metric. Theoretical results show that 1) when the decisions are made with exponentially distributed intervals, the average AuD of the system is smaller if the service time (e.g., transmission time) is uniformly distributed than when it is exponentially distributed, and would be the smallest if it is deterministic; 2)when the decisions are made periodically, the average AuD of the system is larger than, and decreases with decision rate to, the average AuD of the corresponding system with Poisson decision intervals; 3)the probability of missing to use a received update for any decisions is decreasing with the decision rate, and is the smallest if the service time is deterministic. For IoT monitoring systems, therefore, it is suggested to use deterministic monitoring schemes, deterministic transmitting schemes, and Poisson decision schemes, so that the received updates are as fresh as possible at the time they are used to make decisions.


title: Achieving Single-Sensor Complex Activity Recognition from Multi-Sensor Training Data

abstract: In this study, we propose a method for single sensor-based activity recognition, trained with data from multiple sensors. There is no doubt that the performance of complex activity recognition systems increases when we use enough sensors with sufficient quality, however using such rich sensors may not be feasible in real-life situations for various reasons such as user comfort, privacy, battery-preservation, and/or costs. In many cases, only one device such as a smartphone is available, and it is challenging to achieve high accuracy with a single sensor, more so for complex activities. Our method combines representation learning with feature mapping to leverage multiple sensor information made available during training while using a single sensor during testing or in real usage. Our results show that the proposed approach can improve the F1-score of the complex activity recognition by up to 17\% compared to that in training while utilizing the same sensor data in a new user scenario.


title: Automorphic pairs of distributions on $\mathbb{R}$, and Maass forms of real weights

abstract: We give a correspondence between automorphic pairs of distributions on $\mathbb{R}$ and Dirichlet series satisfying functional equations and some additional analytic conditions. Moreover, we show that the notion of automorphic pairs of distributions on $\mathbb{R}$ can be regarded as a generalization of automorphic distributions on smooth principal series representations of the universal covering group of $SL(2,\mathbb{R})$. As an application, we prove Weil type converse theorems for automorphic distributions and Maass forms of real weights.


title: Steady-states of thin film droplets on chemically heterogeneous substrates

abstract: We study steady-state thin films on a chemically heterogeneous substrates of finite size, subject to no-flux boundary conditions. Based on the structure of the bifurcation diagram, we classify the one-dimensional steady-state solutions that exist on such substrates into six different branches and develop asymptotic estimates for the steady-states on each branch. We show using perturbation expansions, that leading order solutions provide good predictions of the steady-state thin films on stepwise-patterned substrates. The analysis in one dimension can be extended to axisymmetric solutions. We also examine the influence of the wettability contrast on linear stability and dynamics. Results are also applied to describe two-dimensional droplets on hydrophilic square patches and striped regions used in microfluidic applications.


title: Stressor-Layer-Induced Elastic Strain Sharing in SrTiO3 Complex Oxide Sheets

abstract: A precisely selected elastic strain can be introduced in submicron-thick single-crystal SrTiO3 sheets using a silicon nitride stressor layer. A conformal stressor layer deposited using plasma-enhanced chemical vapor deposition produces an elastic strain in the sheet consistent with the magnitude of the nitride residual stress. Synchrotron x-ray nanodiffraction reveals that the strain introduced in the SrTiO3 sheets is on the order of 10 4, matching the predictions of an elastic model. This approach to elastic strain sharing in complex oxides allows the strain to be selected within a wide and continuous range of values, an effect not achievable in heteroepitaxy on rigid substrates.


title: Peripheral-free Device Pairing by Randomly Switching Power

abstract: The popularity of Internet-of-Things (IoT) comes with security concerns. Attacks against wireless communication venues of IoT (e.g., Man-in-the-Middle attacks) have grown at an alarming rate over the past decade. Pairing, which allows the establishment of the secure communicating channels for IoT devices without a prior relationship, is thus a paramount capability. Existing secure pairing protocols require auxiliary equipment/peripheral (e.g., displays, speakers and sensors) to achieve authentication, which is unacceptable for low-priced devices such as smart lamps. This paper studies how to design a peripheral-free secure pairing protocol. Concretely, we design the protocol, termed SwitchPairing, via out-of-box power supplying chargers and on-board clocks, achieving security and economics at the same time. When a user wants to pair two or more devices, he/she connects the pairing devices to the same power source, and presses/releases the switch on/off button several times. Then, the press and release timing can be used to derive symmetric keys. We implement a prototype via two CC2640R2F development boards from Texas Instruments (TI) due to its prevalence. Extensive experiments and user studies are also conducted to benchmark our protocol in terms of efficiency and security.


title: LORAX: Loss-Aware Approximations for Energy-Efficient Silicon Photonic Networks-on-Chip

abstract: The approximate computing paradigm advocates for relaxing accuracy goals in applications to improve energy-efficiency and performance. Recently, this paradigm has been explored to improve the energy efficiency of silicon photonic networks-on-chip (PNoCs). In this paper, we propose a novel framework (LORAX) to enable more aggressive approximation during communication over silicon photonic links in PNoCs. Given that silicon photonic interconnects have significant power dissipation due to the laser sources that generate the wavelengths for photonic communication, our framework attempts to reduce laser power overheads while intelligently approximating communication such that application output quality is not distorted beyond an acceptable limit. To the best of our knowledge, this is the first work that considers loss-aware laser power management and multilevel signaling to enable effective data approximation and energy-efficiency in PNoCs. Simulation results show that our framework can achieve up to 31.4% lower laser power consumption and up to 12.2% better energy efficiency than the best known prior work on approximate communication with silicon photonic interconnects, for the same application output quality


title: Dynamical Scattering in Coherent Hard X-Ray Nanobeam Bragg Diffraction

abstract: Unique intensity features arising from dynamical diffraction arise in coherent x-ray nanobeam diffraction patterns of crystals having thicknesses larger than the x-ray extinction depth or exhibiting combinations of nanoscale and mesoscale features. We demonstrate that dynamical scattering effects can be accurately predicted using an optical model combined with the Darwin theory of dynamical x-ray diffraction. The model includes the highly divergent coherent x-ray nanobeams produced by Fresnel zone plate focusing optics and accounts for primary extinction, multiple scattering, and absorption. The simulation accurately reproduces the dynamical scattering features of experimental diffraction patterns acquired from a GaAs/AlGaAs epitaxial heterostructure on a GaAs (001) substrate.


title: Fitting the trajectories of particles in the equatorial plane of a magnetic dipole with epicycloids

abstract: In this paper we discuss epicycloid approximation of the trajectories of charged particles in axisymmetric magnetic fields. Epicycloid trajectories are natural in the Guiding Center approximation and we study in detail the errors arising in this approach. We also discuss how using exact results for particle motion in the equatorial plane of a magnetic dipole the accuracy of this approximation can be significantly extended. We also show that the epicycloids approximate the trajectory of a charged particle more accurately than the position of the particle along the trajectory.


title: Cocycles of $G$-Alexander biquandles and $G$-Alexander multiple conjugation biquandles

abstract: Biquandles and multiple conjugation biquandles are algebras which are related to links and handlebody-links in $3$-space. Cocycles of them can be used to construct state-sum type invariants of links and handlebody-links. In this paper we discuss cocycles of a certain class of biquandles and multiple conjugation biquandles, which we call $G$-Alexander biquandles and $G$-Alexander multiple conjugation biquandles, with a relationship with group cocycles. We give a method to obtain a (biquandle or multiple conjugation biquandle) cocycle of them from a group cocycle.


title: Adversarial Ranking Attack and Defense

abstract: Deep Neural Network (DNN) classifiers are vulnerable to adversarial attack, where an imperceptible perturbation could result in misclassification. However, the vulnerability of DNN-based image ranking systems remains under-explored. In this paper, we propose two attacks against deep ranking systems, i.e., Candidate Attack and Query Attack, that can raise or lower the rank of chosen candidates by adversarial perturbations. Specifically, the expected ranking order is first represented as a set of inequalities, and then a triplet-like objective function is designed to obtain the optimal perturbation. Conversely, a defense method is also proposed to improve the ranking system robustness, which can mitigate all the proposed attacks simultaneously. Our adversarial ranking attacks and defense are evaluated on datasets including MNIST, Fashion-MNIST, and Stanford-Online-Products. Experimental results demonstrate that a typical deep ranking system can be effectively compromised by our attacks. Meanwhile, the system robustness can be moderately improved with our defense. Furthermore, the transferable and universal properties of our adversary illustrate the possibility of realistic black-box attack.


title: Geometry of varieties for graded maximal Cohen--Macaulay modules

abstract: We study a variety for graded maximal Cohen--Macaulay modules, which was introduced by Dao and Shipman. The main result of this paper asserts that there are only a finite number of isomorphism classes of graded maximal Cohen--Macaulay modules with fixed Hilbert series over Cohen--Macaulay algebras of graded countable representation type.


title: Machine Learning based prediction of noncentrosymmetric crystal materials

abstract: Noncentrosymmetric materials play a critical role in many important applications such as laser technology, communication systems,quantum computing, cybersecurity, and etc. However, the experimental discovery of new noncentrosymmetric materials is extremely difficult. Here we present a machine learning model that could predict whether the composition of a potential crystalline structure would be centrosymmetric or not. By evaluating a diverse set of composition features calculated using matminer featurizer package coupled with different machine learning algorithms, we find that Random Forest Classifiers give the best performance for noncentrosymmetric material prediction, reaching an accuracy of 84.8% when evaluated with 10 fold cross-validation on the dataset with 82,506 samples extracted from Materials Project. A random forest model trained with materials with only 3 elements gives even higher accuracy of 86.9%. We apply our ML model to screen potential noncentrosymmetric materials from 2,000,000 hypothetical materials generated by our inverse design engine and report the top 20 candidate noncentrosymmetric materials with 2 to 4 elements and top 20 borate candidates


title: Sparse Sinkhorn Attention

abstract: We propose Sparse Sinkhorn Attention, a new efficient and sparse method for learning to attend. Our method is based on differentiable sorting of internal representations. Concretely, we introduce a meta sorting network that learns to generate latent permutations over sequences. Given sorted sequences, we are then able to compute quasi-global attention with only local windows, improving the memory efficiency of the attention module. To this end, we propose new algorithmic innovations such as Causal Sinkhorn Balancing and SortCut, a dynamic sequence truncation method for tailoring Sinkhorn Attention for encoding and/or decoding purposes. Via extensive experiments on algorithmic seq2seq sorting, language modeling, pixel-wise image generation, document classification and natural language inference, we demonstrate that our memory efficient Sinkhorn Attention method is competitive with vanilla attention and consistently outperforms recently proposed efficient Transformer models such as Sparse Transformers.


title: Generalized ODIN: Detecting Out-of-distribution Image without Learning from Out-of-distribution Data

abstract: Deep neural networks have attained remarkable performance when applied to data that comes from the same distribution as that of the training set, but can significantly degrade otherwise. Therefore, detecting whether an example is out-of-distribution (OoD) is crucial to enable a system that can reject such samples or alert users. Recent works have made significant progress on OoD benchmarks consisting of small image datasets. However, many recent methods based on neural networks rely on training or tuning with both in-distribution and out-of-distribution data. The latter is generally hard to define a-priori, and its selection can easily bias the learning. We base our work on a popular method ODIN, proposing two strategies for freeing it from the needs of tuning with OoD data, while improving its OoD detection performance. We specifically propose to decompose confidence scoring as well as a modified input pre-processing method. We show that both of these significantly help in detection performance. Our further analysis on a larger scale image dataset shows that the two types of distribution shifts, specifically semantic shift and non-semantic shift, present a significant difference in the difficulty of the problem, providing an analysis of when ODIN-like strategies do or do not work.


title: Theory for Self-Bound States of Dipolar Bose-Einstein Condensates

abstract: We investigate the self-bound states of dipolar Dy condensates with the Gaussian-state ansatz which improves the conventional coherent-state ansatz with multimode squeezed coherent states. We show that the self-bound states consist of the experimentally observed self-bound liquid phase and the unobserved self-bound gas phase. The numerically obtained gas-liquid boundary is in good agreement with experimental data. Our theory also allows one to extract the real part of the three-body coupling constant of the Dy atoms from the particle number distribution of the condensates. In particular, we results show that the self-bound states are stabilized by the short-range three-body repulsion. Our study shed a different light to understand the self-bound droplets of Bose-Einstein condensates.


title: Embeddings of $4$--manifolds in $\CP^3$

abstract: In this article we show that every closed orientable smooth $4$--manifold admits a smooth embedding in the complex projective $3$--space.


title: Self-supervised Image Enhancement Network: Training with Low Light Images Only

abstract: This paper proposes a self-supervised low light image enhancement method based on deep learning. Inspired by information entropy theory and Retinex model, we proposed a maximum entropy based Retinex model. With this model, a very simple network can separate the illumination and reflectance, and the network can be trained with low light images only. We introduce a constraint that the maximum channel of the reflectance conforms to the maximum channel of the low light image and its entropy should be largest in our model to achieve self-supervised learning. Our model is very simple and does not rely on any well-designed data set (even one low light image can complete the training). The network only needs minute-level training to achieve image enhancement. It can be proved through experiments that the proposed method has reached the state-of-the-art in terms of processing speed and effect.


title: Quantifying Convective Aggregation using the Tropical Moist Margin's Length

abstract: On small scales, the tropical atmosphere tends to be either moist or very dry. This defines two states that, on large scales, are separated by a sharp margin, well-identified by the anti-mode of the bimodal tropical column water vapor distribution. Despite recent progress in understanding physical processes governing the spatio-temporal variability of tropical water vapor, the behavior of this margin remains elusive, and we lack a simple framework to understand the bimodality of tropical water vapor in observations. Motivated by the success of coarsening theory in explaining bimodal distributions, we leverage its methodology to relate the moisture field's spatial organization to its time-evolution. This results in a new diagnostic framework for the bimodality of tropical water vapor, from which we argue that the length of the margin separating moist from dry regions should evolve towards a minimum in equilibrium. As the spatial organization of moisture is closely related to the organization of tropical convection, we hereby introduce a new organization index (BLW) measuring the ratio of the margin's length to the circumference of a well-defined equilibrium shape. Using BLW, we assess the evolution of self-aggregation in idealized cloud-resolving simulations of radiative-convective equilibrium and contrast it to the time-evolution of the Atlantic Inter-Tropical Convergence Zone (ITCZ) in the ERA5 meteorological re-analysis product. We find that BLW successfully captures aspects of convective organization ignored by more traditional metrics, while offering a new perpective on the seasonal cycle of convective organization in the Atlantic ITCZ.


title: Bandwidth-Optimized Parallel Algorithms for Sparse Matrix-Matrix Multiplication using Propagation Blocking

abstract: Sparse matrix-matrix multiplication (SpGEMM) is a widely used kernel in various graph, scientific computing and machine learning algorithms. It is well known that SpGEMM is a memory-bound operation, and its peak performance is expected to be bound by the memory bandwidth. Yet, existing algorithms fail to saturate the memory bandwidth, resulting in suboptimal performance under the Roofline model. In this paper we characterize existing SpGEMM algorithms based on their memory access patterns and develop practical lower and upper bounds for SpGEMM performance. We then develop an SpGEMM algorithm based on outer product matrix multiplication. The newly developed algorithm called PB-SpGEMM saturates memory bandwidth by using the propagation blocking technique and by performing in-cache sorting and merging. For many practical matrices, PB-SpGEMM runs 20%-50% faster than the state-of-the-art heap and hash SpGEMM algorithms on modern multicore processors. Most importantly, PB-SpGEMM attains performance predicted by the Roofline model, and its performance remains stable with respect to matrix size and sparsity.