title: Epitaxial binding and strain effects of monolayer stanene on the Al$_{2}$O$_{3}$(0001) surface

abstract: Stanene, the two-dimensional monolayer form of tin, has been predicted to be a 2D topological insulator due to its large spin--orbit interaction. However, a clear experimental demonstration of stanene's topologically nontrivial properties has eluded observation, in part because of the difficulty of choosing a substrate on which stanene will remain topologically nontrivial. In this paper, we present first-principles density functional theory (DFT) calculations of epitaxial monolayer stanene grown on the (0001) surface of alumina, Al$_{2}$O$_{3}$, as well as free-standing decorated stanene under strain. By describing the energetics and nature of how monolayer stanene binds to alumina, we show a strong energetic drive for the monolayer to be coherently strained and epitaxial to the substrate. By analyzing the electronic structure of strained stanene, we find it to be a quantum spin Hall insulator on Al$_{2}$O$_{3}$. We also describe the effect of \emph{in situ} fluorine decoration on the bound stanene monolayer, including on its potential for mechanical exfoliation.


title: Distributed Cooperation Under Uncertainty in Drone-Based Wireless Networks: A Bayesian Coalitional Game

abstract: We study the resource sharing problem in a drone-based wireless network. We consider a distributed control setting under uncertainty (i.e. unavailability of full information). In particular, the drones cooperate in serving the users while pooling their spectrum and energy resources in the absence of prior knowledge about different system characteristics such as the amount of available power at the other drones. We cast the aforementioned problem as a Bayesian cooperative game in which the agents (drones) engage in a coalition formation process, where the goal is to maximize the overall transmission rate of the network. The drones update their beliefs using a novel technique that combines the maximum likelihood estimation with Kullback-Leibler divergence. We propose a decision-making strategy for repeated coalition formation that converges to a stable coalition structure. We analyze the performance of the proposed approach by both theoretical analysis and simulations.


title: A revision of the vdB 130 cluster stellar content based on GAIA DR2 Data. Interstellar extinction toward the Cyg OB1 supershell

abstract: Two star-forming regions are studied: the young embedded open cluster vdB 130 and the protocluster neighbourhood observed in the head and tail of the cometary molecular cloud located in the wall of the expanding supershell surrounding the Cyg OB1 association. The GAIA DR2 catalogue is employed to verify the stellar composition of the vdB 130 cluster whose members were earlier selected using the UCAC4 catalogue. The new sample of vdB 130 members contains 68 stars with close proper motions (within 1 mas yr$^{-1}$) and close trigonometric parallaxes (ranging from 0.50 to 0.70 mas). The relative parallax error is shown to increase with distance to objects and depend on their magnitude. At a distance of 1.5-2 kpc it is of about 3-7 per cent and 20-30 per cent for bright and faint stars, respectively. The cluster is not older than ~10 Myr. New spectroscopic and photometric observations carried out on Russian telescopes are combined with GAIA DR2 to search for optical components in the protocluster region - a new starburst. An analysis of 20 stars in the vicinity of the protocluster revealed no concentration of either proper motions or parallaxes. According to spectroscopic, photometric, and trigonometric estimates, the distances to these stars range from 0.4 to 2.5 kpc, and colour excess is shown to increase with a distance D (kpc) in accordance with the law: $E(B-V)\simeq 0.6\times D$ mag.


title: The ExoMolOP Database: Cross-sections and k-tables for Molecules of Interest in High-Temperature Exoplanet Atmospheres

abstract: A publicly available database of opacities for molecules of astrophysical interest, ExoMolOP, has been compiled for over 80 species, based on the latest line list data from the ExoMol, HITEMP and MoLLIST databases. These data are generally suitable for characterising high temperature exoplanet or cool stellar/substellar atmospheres, and have been computed at a variety of pressures and temperatures, with a few molecules included at room-temperature only from the HITRAN database. The data are formatted in different ways for four different exoplanet atmosphere retrieval codes; ARCiS, TauREx, NEMESIS and petitRADTRANS, and include both cross-sections (at R~=~$\frac{\lambda}{\Delta \lambda}$~=~15,000) and k-tables (at R~=~$\frac{\lambda}{\Delta \lambda}$~=~1000) for the 0.3~-~50$\mu$m wavelength region. Opacity files can be downloaded and used directly for these codes. Atomic data for alkali metals Na and K are also included, using data from the NIST database and the latest line shapes for the resonance lines. Broadening parameters have been taken from the literature where available, or from those for a known molecule with similar molecular properties where no broadening data are available. The data are available from www.exomol.com.


title: Graph neural network for 3D classification of ambiguities and optical crosstalk in scintillator-based neutrino detectors

abstract: Deep learning tools are being used extensively in high energy physics and are becoming central in the reconstruction of neutrino interactions in particle detectors. In this work, we report on the performance of a graph neural network in assisting with particle flow event reconstruction. The three-dimensional reconstruction of particle tracks produced in neutrino interactions can be subject to ambiguities due to high multiplicity signatures in the detector or leakage of signal between neighboring active detector volumes. Graph neural networks potentially have the capability of identifying all these features to boost the reconstruction performance. As an example case study, we tested a graph neural network, inspired by the GraphSAGE algorithm, on a novel 3D-granular plastic-scintillator detector, that will be used to upgrade the near detector of the T2K experiment. The developed neural network has been trained and tested on diverse neutrino interaction samples, showing very promising results: the classification of particle track voxels produced in the detector can be done with efficiencies and purities of 94-96% per event and most of the ambiguities can be identified and rejected, while being robust against systematic effects.


title: Multi-pulse Photoacoustic Microscopy Based on k-space Pseudospectral Method: Simulation Study

abstract: Photoacoustic imaging has emerged in the past decades. Compared with the traditional medical imaging mode, it has better imaging performance and has great development potential in the field of biological imaging. In traditional photoacoustic microscopy, a single laser pulse is generally used to irradiate the sample to produce photoacoustic signal. And signal-to-noise ratio (SNR) is a very important indicator for photoacoustic imaging. In order to obtain the image with high SNR, multiple acquisition or increasing laser pulse energy is usually adopted. The former will lead to slower imaging speed, and the latter will lead to photobleaching or phototoxicity. Here, we proposes multi-pulse photoacoustic microscopy, the photoacoustic signals were stimulated sequentially using multiple laser pulses in each A line data acquisition. In order to verify the feasibility of this method, a multi-pulse photoacoustic imaging simulation platform is established using k-Wave toolbox. The performance of multi-pulse photoacoustic imaging is verified through the three scanning modes of photoacoustic microscopy A-scan, B -scan, and C-scan. The results indicate that the SNR is proportion to the number of laser pulses used, high SNR can be achieved by low-energy laser pulse. This work will help to expand the application of photoacoustic imaging.


title: Improved Bilevel Model: Fast and Optimal Algorithm with Theoretical Guarantee

abstract: Due to the hierarchical structure of many machine learning problems, bilevel programming is becoming more and more important recently, however, the complicated correlation between the inner and outer problem makes it extremely challenging to solve. Although several intuitive algorithms based on the automatic differentiation have been proposed and obtained success in some applications, not much attention has been paid to finding the optimal formulation of the bilevel model. Whether there exists a better formulation is still an open problem. In this paper, we propose an improved bilevel model which converges faster and better compared to the current formulation. We provide theoretical guarantee and evaluation results over two tasks: Data Hyper-Cleaning and Hyper Representation Learning. The empirical results show that our model outperforms the current bilevel model with a great margin. \emph{This is a concurrent work with \citet{liu2020generic} and we submitted to ICML 2020. Now we put it on the arxiv for record.}


title: Bad places for the approximation property for finite groups

abstract: Given a number field $k$ and a finite $k$-group $G$, the Tame Approximation Problem for $G$ asks whether the restriction map $H^1(k,G)\to\prod_{v\in\Sigma}H^1(k_v,G)$ is surjective for every finite set of places $\Sigma\subseteq\Omega_k$ disjoint from $\text{Bad}_G$, where $\text{Bad}_G$ is the finite set of places that either divides the order of $G$ or ramifies in the minimal extension splitting $G$. In this paper we prove that the set $\text{Bad}_G$ is "sharp". To achieve this we prove that there are finite abelian $k$-groups $A$ where the map $H^1(k,A)\to\prod_{v\in\Sigma_0}H^1(k_v,A)$ is not surjective in a set $\Sigma_0\subseteq\text{Bad}_A$ with particular properties, namely $\Sigma_0$ is the set of places that do not divide the order of $A$ and ramify in the minimal extension splitting $A$.


title: Amplification of brightness variability by active-region nesting in solar-like stars

abstract: Kepler observations revealed that hundreds of stars with near-solar fundamental parameters and rotation periods have much stronger and more regular brightness variations than the Sun. Here we identify one possible reason for the peculiar behaviour of these stars. Inspired by solar nests of activity, we assume that the degree of inhomogeneity of active-region (AR) emergence on such stars is higher than on the Sun. To test our hypothesis, we model stellar light curves by injecting ARs consisting of spots and faculae on stellar surfaces at various rates and nesting patterns, using solar AR properties and differential rotation. We show that a moderate increase of the emergence frequency from the solar value combined with the increase of the degree of nesting can explain the full range of observed amplitudes of variability of Sun-like stars with nearly the solar rotation period. Furthermore, nesting in the form of active longitudes, in which ARs tend to emerge in the vicinity of two longitudes separated by $180^\circ$, leads to highly regular, almost sine-like variability patterns, rather similar to those observed in a number of solar-like stars.


title: Most Generalized Petersen graphs of girth 8 have cop number 4

abstract: A generalized Petersen graph $GP(n,k)$ is a regular cubic graph on $2n$ vertices (the parameter $k$ is used to define some of the edges). It was previously shown (Ball et al., 2015) that the cop number of $GP(n,k)$ is at most $4$, for all permissible values of $n$ and $k$. In this paper we prove that the cop number of "most" generalized Petersen graphs is exactly $4$. More precisely, we show that unless $n$ and $k$ fall into certain specified categories, then the cop number of $GP(n,k)$ is $4$. The graphs to which our result applies all have girth $8$. In fact, our argument is slightly more general: we show that in any cubic graph of girth at least $8$, unless there exist two cycles of length $8$ whose intersection is a path of length $2$, then the cop number of the graph is at least $4$. Even more generally, in a graph of girth at least $9$ and minimum valency $\delta$, the cop number is at least $\delta+1$.


title: Automatic Assignment of Radiology Examination Protocols Using Pre-trained Language Models with Knowledge Distillation

abstract: Selecting radiology examination protocol is a repetitive, and time-consuming process. In this paper, we present a deep learning approach to automatically assign protocols to computer tomography examinations, by pre-training a domain-specific BERT model ($BERT_{rad}$). To handle the high data imbalance across exam protocols, we used a knowledge distillation approach that up-sampled the minority classes through data augmentation. We compared classification performance of the described approach with the statistical n-gram models using Support Vector Machine (SVM), Gradient Boosting Machine (GBM), and Random Forest (RF) classifiers, as well as the Google's $BERT_{base}$ model. SVM, GBM and RF achieved macro-averaged F1 scores of 0.45, 0.45, and 0.6 while $BERT_{base}$ and $BERT_{rad}$ achieved 0.61 and 0.63. Knowledge distillation improved overall performance on the minority classes, achieving a F1 score of 0.66.


title: Closure under Coupling of Cellular-DEVS for the Optimization of Memory Resource: Wildfire Spread Case Study

abstract: The present work aims to show one of the advantages of using the property of closure under coupling in the DEVS specification. The advantage concerned in this paper attempts to address the need for memory resources during the simulation of systems by cellular-DEVS. This improvement of performance is based on the usage of the property closure under coupling in the DEVS formalism. With this property and taking account of the iterative behavior of each cellular-DEVS atomic model, we provide simulation of many models simultaneously. The method starts with the specification of the cellular-DEVS coupled model which is then converted into its equivalent DEVS atomic model. Thus, the goal of this conversion is to transform large quantities of atomic models coupled together, which require huge computational resources, into one DEVS atomic model. A case study is presented at the end of the work on modeling and simulation of forest fire propagation using DEVS and cellular-DEVS. A specification by cellular-DEVS of the forest fire model and its non-modular equivalent DEVS atomic model are presented. Finally, a comparison between both methods is presented in term of consumption of resources.


title: Conley Index Theory and the Attractor-Repeller Decomposition for Differential Inclusions

abstract: The Conley index theory is a powerful topological tool for describing the basic structure of dynamical systems. One important feature of this theory is the attractor-repeller decomposition of isolated invariant sets. In this decomposition, all points in the invariant set belong to the attractor, its associated dual repeller, or a connecting region. In this connecting region, points tend towards the attractor in forwards time and the repeller in backwards time. This decomposition is also, in a certain topological sense, stable under perturbation. Conley theory is well-developed for flows and homomorphisms, and has also been extended to some more abstract settings such as semiflows and relations. In this paper we aim to extend the attractor-repeller decomposition, including its stability under perturbation, to continuous time set-valued dynamical systems. The most common of these systems are differential inclusions such as Filippov systems.


title: Simpler and Stronger Approaches for Non-Uniform Hypergraph Matching and the F\"uredi, Kahn, and Seymour Conjecture

abstract: A well-known conjecture of F\"uredi, Kahn, and Seymour (1993) on non-uniform hypergraph matching states that for any hypergraph with edge weights $w$, there exists a matching $M$ such that the inequality $\sum_{e\in M} g(e) w(e) \geq \mathrm{OPT}_{\mathrm{LP}}$ holds with $g(e)=|e|-1+\frac{1}{|e|}$, where $\mathrm{OPT}_{\mathrm{LP}}$ denotes the optimal value of the canonical LP relaxation. While the conjecture remains open, the strongest result towards it was very recently obtained by Brubach, Sankararaman, Srinivasan, and Xu (2020)---building on and strengthening prior work by Bansal, Gupta, Li, Mestre, Nagarajan, and Rudra (2012)---showing that the aforementioned inequality holds with $g(e)=|e|+O(|e|\exp(-|e|))$. Actually, their method works in a more general sampling setting, where, given a point $x$ of the canonical LP relaxation, the task is to efficiently sample a matching $M$ containing each edge $e$ with probability at least $\frac{x(e)}{g(e)}$. We present simpler and easy-to-analyze procedures leading to improved results. More precisely, for any solution $x$ to the canonical LP, we introduce a simple algorithm based on exponential clocks for Brubach et al.'s sampling setting achieving $g(e)=|e|-(|e|-1)x(e)$. Apart from the slight improvement in $g$, our technique may open up new ways to attack the original conjecture. Moreover, we provide a short and arguably elegant analysis showing that a natural greedy approach for the original setting of the conjecture shows the inequality for the same $g(e)=|e|-(|e|-1)x(e)$ even for the more general hypergraph $b$-matching problem.


title: The Bramson delay in a Fisher-KPP equation with log-singular non-linearity

abstract: We consider a class of reaction-diffusion equations of Fisher-KPP type in which the nonlinearity (reaction term) $f$ is merely $C^1$ at $u=0$ due to a logarithmic competition term. We first derive the asymptotic behavior of (minimal speed) traveling wave solutions that is, we obtain precise estimates on the decay to zero of the traveling wave profile at infinity. We then use this to characterize the Bramson shift between the traveling wave solutions and solutions of the Cauchy problem with localized initial data. We find a phase transition depending on how singular $f$ is near $u=0$ with quite different behavior for more singular $f$. This is in contrast to the smooth case, that is, when $f \in C^{1,\delta}$, where these behaviors are completely determined by $f'(0)$. In the singular case, several scales appear and require new techniques to understand.


title: On Generalised Petersen Graphs of Girth 7 that have Cop Number 4

abstract: We show that if $n=7k/i$ with $i \in \{1,2,3\}$ then the cop number of the generalised Petersen graph $GP(n,k)$ is $4$, with some small previously-known exceptions. It was previously proved by Ball et al. (2015) that the cop number of any generalised Petersen graph is at most $4$. The results in this paper explain all of the known generalised Petersen graphs that actually have cop number $4$ but were not previously explained by Morris et al. in a recent preprint, and places them in the context of infinite families. (More precisely, the preprint by Morris et al. explains all known generalised Petersen graphs with cop number $4$ and girth $8$, while this paper explains those that have girth $7$.)


title: Multimodal Inductive Transfer Learning for Detection of Alzheimer's Dementia and its Severity

abstract: Alzheimer's disease is estimated to affect around 50 million people worldwide and is rising rapidly, with a global economic burden of nearly a trillion dollars. This calls for scalable, cost-effective, and robust methods for detection of Alzheimer's dementia (AD). We present a novel architecture that leverages acoustic, cognitive, and linguistic features to form a multimodal ensemble system. It uses specialized artificial neural networks with temporal characteristics to detect AD and its severity, which is reflected through Mini-Mental State Exam (MMSE) scores. We first evaluate it on the ADReSS challenge dataset, which is a subject-independent and balanced dataset matched for age and gender to mitigate biases, and is available through DementiaBank. Our system achieves state-of-the-art test accuracy, precision, recall, and F1-score of 83.3% each for AD classification, and state-of-the-art test root mean squared error (RMSE) of 4.60 for MMSE score regression. To the best of our knowledge, the system further achieves state-of-the-art AD classification accuracy of 88.0% when evaluated on the full benchmark DementiaBank Pitt database. Our work highlights the applicability and transferability of spontaneous speech to produce a robust inductive transfer learning model, and demonstrates generalizability through a task-agnostic feature-space. The source code is available at https://github.com/wazeerzulfikar/alzheimers-dementia


title: A new electromechanical analogy approach based on electrostatic coupling for vertical dynamic analysis of planar vehicle models

abstract: Analogies between mechanical and electrical systems have been developed and applied for almost a century, and they have proved their usefulness in the study of mechanical and electrical systems. The development of new elements such as the inerter or the memristor is a clear example. However, new applications and possibilities of using these analogues still remain to be explored. In this work, the electrical analogues of different vehicle models are presented. A new and not previously reported analogy between inertial coupling and electrostatic capacitive coupling is found and described. Several examples are provided to highlight the benefits of this analogy. Well-known mechanical systems like the half-car or three three-axle vehicle models are discussed and some numerical results are presented. To the best of the author's knowledge, such systems were never dealt with by using a full electromechanical analogy. The mechanical equations are also derived and compared with those of the electrical domain for harmonic steady state analysis.


title: Unsupervised Single-Image Reflection Separation Using Perceptual Deep Image Priors

abstract: Reflections often degrade the quality of the image by obstructing the background scene. This is not desirable for everyday users, and it negatively impacts the performance of multimedia applications that process images with reflections. Most current methods for removing reflections utilize supervised-learning models. However, these models require an extensive number of image pairs to perform well, especially on natural images with reflection, which is difficult to achieve in practice. In this paper, we propose a novel unsupervised framework for single-image reflection separation. Instead of learning from a large dataset, we optimize the parameters of two cross-coupled deep convolutional networks on a target image to generate two exclusive background and reflection layers. In particular, we design a new architecture of the network to embed semantic features extracted from a pre-trained deep classification network, which gives more meaningful separation similar to human perception. Quantitative and qualitative results on commonly used datasets in the literature show that our method's performance is at least on par with the state-of-the-art supervised methods and, occasionally, better without requiring large training datasets. Our results also show that our method significantly outperforms the closest unsupervised method in the literature for removing reflections from single images.


title: On tomography velocity uncertainty in relation with structural imaging

abstract: Evaluating structural uncertainties associated with seismic imaging and target horizons can be of critical importance for decision-making related to oil and gas exploration and production. An important breakthrough for industrial applications has been made with the development of industrial approaches to velocity model building. We propose an extension of these approaches, sampling an equi-probable contour of the tomography posterior probability density function (pdf) rather than the full pdf, and using non-linear slope tomography (rather than standard tomographic migration velocity analysis as in previous publications). Our approach allows to assess the quality of uncertainty-related assumptions (linearity and Gaussian hypothesis within the Bayesian theory) and estimate volumetric migration positioning uncertainties (a generalization of horizon uncertainties), in addition to the advantages in terms of efficiency. We derive the theoretical concepts underlying this approach and unify our derivations with those of previous publications. As the method works in the full model space rather than in a preconditioned model space, we split the analysis into the resolved and unresolved tomography spaces. We argue that the resolved space uncertainties are to be used in further steps leading to decision-making and can be related to the output of methods that work in a preconditioned model space. The unresolved space uncertainties represent a qualitative byproduct specific to our method, strongly highlighting the most uncertain gross areas, thus useful for QCs. These concepts are demonstrated on a synthetic data. Complementarily, the industrial viability of the method is illustrated on two different 3D field datasets.


title: Superconvergent Interpolatory HDG methods for reaction diffusion equations II: HHO-inspired methods

abstract: In J. Sci. Comput., 81: 2188-2212, 2019, we considered a superconvergent hybridizable discontinuous Galerkin (HDG) method, defined on simplicial meshes, for scalar reaction diffusion equations and showed how to define an interpolatory version which maintained its convergence properties. The interpolatory approach uses a locally postprocessed approximate solution to evaluate the nonlinear term, and assembles all HDG matrices once before the time intergration leading to a reduction in computational cost. The resulting method displays a superconvergent rate for the solution for polynomial degree $k\ge 1$. In this work, we take advantage of the link found between the HDG and the hybrid high-order (HHO) methods, in ESAIM Math. Model. Numer. Anal., 50: 635-650, 2016 and extend this idea to the new, HHO-inspired HDG methods, defined on meshes made of general polyhedral elements, uncovered therein}. We prove that the resulting interpolatory HDG methods converge at the same rate as for the linear elliptic problems. Hence, we obtain superconvergent methods for $k\ge 0$ by some methods. We present numerical results to illustrate the convergence theory.


title: Efficient p-multigrid spectral element model for water waves and marine offshore structures

abstract: In marine offshore engineering, cost-efficient simulation of unsteady water waves and their nonlinear interaction with bodies are important to address a broad range of engineering applications at increasing fidelity and scale. We consider a fully nonlinear potential flow (FNPF) model discretized using a Galerkin spectral element method to serve as a basis for handling both wave propagation and wave-body interaction with high computational efficiency within a single modellingapproach. We design and propose an efficientO(n)-scalable computational procedure based on geometric p-multigrid for solving the Laplace problem in the numerical scheme. The fluid volume and the geometric features of complex bodies is represented accurately using high-order polynomial basis functions and unstructured meshes with curvilinear prism elements. The new p-multigrid spectralelement model can take advantage of the high-order polynomial basis and thereby avoid generating a hierarchy of geometric meshes with changing number of elements as required in geometric h-multigrid approaches. We provide numerical benchmarks for the algorithmic and numerical efficiency of the iterative geometric p-multigrid solver. Results of numerical experiments are presented for wave propagation and for wave-body interaction in an advanced case for focusing design waves interacting with a FPSO. Our study shows, that the use of iterative geometric p-multigrid methods for theLaplace problem can significantly improve run-time efficiency of FNPF simulators.


title: Industrial scale large eddy simulations (LES) with adaptive octree meshes using immersogeometric analysis

abstract: We present a variant of the immersed boundary method integrated with octree meshes for highly efficient and accurate Large-Eddy Simulations (LES) of flows around complex geometries. We demonstrate the scalability of the proposed method up to $\mathcal{O}(32K)$ processors. This is achieved by (a) rapid in-out tests; (b) adaptive quadrature for an accurate evaluation of forces; (c) tensorized evaluation during matrix assembly. We showcase this method on two non-trivial applications: accurately computing the drag coefficient of a sphere across Reynolds numbers $1-10^6$ encompassing the drag crisis regime; simulating flow features across a semi-truck for investigating the effect of platooning on efficiency.


title: Pursuing a Prospective Perspective

abstract: Retrospective testing of predictive models does not consider the real-world context in which models are deployed. Prospective validation, on the other hand, enables meaningful comparisons between data generation processes by incorporating trained models and considering the subjective decisions that affect reproducibility. Prospective experiments are essential for consistent progress in modeling.


title: The convexity of optimal transport-based waveform inversion for certain structured velocity models

abstract: Full--waveform inversion (FWI) is a method used to determine properties of the Earth from information on the surface. We use the squared Wasserstein distance (squared $W_2$ distance) as an objective function to invert for the velocity of seismic waves as a function of position in the Earth, and we discuss its convexity with respect to the velocity parameter. In one dimension, we consider constant, piecewise increasing, and linearly increasing velocity models as a function of position, and we show the convexity of the squared $W_2$ distance with respect to the velocity parameter on the interval from zero to the true value of the velocity parameter when the source function is a probability measure. Furthermore, we consider a two--dimensional model where velocity is linearly increasing as a function of depth and prove the convexity of the squared $W_2$ distance in the velocity parameter on large regions containing the true value. We discuss the convexity of the squared $W_2$ distance compared with the convexity of the squared $L^2$ norm, and we discuss the relationship between frequency and convexity of these respective distances. We also discuss multiple approaches to optimal transport for non--probability measures by first converting the wave data into probability measures.


title: Equation of state of atomic solid hydrogen by stochastic many-body wave function methods

abstract: We report a numerical study of the equation of state of crystalline body-centered-cubic (BCC) hydrogen, tackled with a variety of complementary many-body wave function methods. These include continuum stochastic techniques of fixed-node diffusion and variational quantum Monte Carlo, and the Hilbert space stochastic method of full configuration-interaction quantum Monte Carlo. In addition, periodic coupled-cluster methods were also employed. Each of these methods is underpinned with different strengths and approximations, but their combination in order to perform reliable extrapolation to complete basis set and supercell size limits gives confidence in the final results. The methods were found to be in good agreement for equilibrium cell volumes for the system in the BCC phase, with a lattice parameter of 3.307 Bohr.

categories: cond-mat.mtrl-sci cond-mat.other physics.chem-ph physics.comp-ph


title: Invariant conservative difference schemes for shallow water equations in Eulerian and Lagrangian coordinates

abstract: The one-dimensional shallow water equations in Eulerian coordinates are considered. Relations between symmetries and conservation laws for the potential form of the equations, and symmetries and conservation laws in Eulerian coordinates are shown. An invariant difference scheme for equations in Eulerian coordinates with arbitrary bottom topography is constructed. It possesses all the finite-difference analogues of the conservation laws. Some bottom topographies require moving meshes in Eulerian coordinates, which are stationary meshes in mass Lagrangian coordinates. The developed invariant conservative difference schemes are verified numerically using examples of flow with various bottom topographies.

categories: math.NA cs.NA math-ph math.MP physics.flu-dyn


title: Uniformly bounded Lebesgue constants for scaled cardinal interpolation with Mat\'{e}rn kernels

abstract: For $h>0$ and positive integers $m$, $d$, such that $m>d/2$, we study non-stationary interpolation at the points of the scaled grid $h\mathbb{Z}^d$ via the Mat\'{e}rn kernel $\Phi_{m,d}$---the fundamental solution of $(1-\Delta)^m$ in $\mathbb{R}^d$. We prove that the Lebesgue constants of the corresponding interpolation operators are uniformly bounded as $h\to0$ and deduce the convergence rate $O(h^{2m})$ for the scaled interpolation scheme. We also provide convergence results for approximation with Mat\'{e}rn and related compactly supported polyharmonic kernels.


title: Short-term Traffic Prediction with Deep Neural Networks: A Survey

abstract: In modern transportation systems, an enormous amount of traffic data is generated every day. This has led to rapid progress in short-term traffic prediction (STTP), in which deep learning methods have recently been applied. In traffic networks with complex spatiotemporal relationships, deep neural networks (DNNs) often perform well because they are capable of automatically extracting the most important features and patterns. In this study, we survey recent STTP studies applying deep networks from four perspectives. 1) We summarize input data representation methods according to the number and type of spatial and temporal dependencies involved. 2) We briefly explain a wide range of DNN techniques from the earliest networks, including Restricted Boltzmann Machines, to the most recent, including graph-based and meta-learning networks. 3) We summarize previous STTP studies in terms of the type of DNN techniques, application area, dataset and code availability, and the type of the represented spatiotemporal dependencies. 4) We compile public traffic datasets that are popular and can be used as the standard benchmarks. Finally, we suggest challenging issues and possible future research directions in STTP.


title: WaveGrad: Estimating Gradients for Waveform Generation

abstract: This paper introduces WaveGrad, a conditional model for waveform generation which estimates gradients of the data density. The model is built on prior work on score matching and diffusion probabilistic models. It starts from a Gaussian white noise signal and iteratively refines the signal via a gradient-based sampler conditioned on the mel-spectrogram. WaveGrad offers a natural way to trade inference speed for sample quality by adjusting the number of refinement steps, and bridges the gap between non-autoregressive and autoregressive models in terms of audio quality. We find that it can generate high fidelity audio samples using as few as six iterations. Experiments reveal WaveGrad to generate high fidelity audio, outperforming adversarial non-autoregressive baselines and matching a strong likelihood-based autoregressive baseline using fewer sequential operations. Audio samples are available at https://wavegrad.github.io/.


title: The Dirichlet isospectral problem for trapezoids

abstract: We show that non-obtuse trapezoids are uniquely determined by their Dirichlet Laplace spectrum. This extends our previous result, which was only concerned with the Neumann Laplace spectrum.

categories: math.AP math-ph math.DG math.DS math.MP math.SP


title: A Survey on Recent Hardware Data Prefetching Approaches with An Emphasis on Servers

abstract: Data prefetching, i.e., the act of predicting application's future memory accesses and fetching those that are not in the on-chip caches, is a well-known and widely-used approach to hide the long latency of memory accesses. The fruitfulness of data prefetching is evident to both industry and academy: nowadays, almost every high-performance processor incorporates a few data prefetchers for capturing various access patterns of applications; besides, there is a myriad of proposals for data prefetching in the research literature, where each proposal enhances the efficiency of prefetching in a specific way. In this survey, we discuss the fundamental concepts in data prefetching and study state-of-the-art hardware data prefetching approaches. Additional Key Words and Phrases: Data Prefetching, Scale-Out Workloads, Server Processors, and Spatio-Temporal Correlation.


title: MAKE: a Matrix Action Key Exchange

abstract: We offer a public key exchange protocol based on a semidirect product of two cyclic (semi)groups of matrices over Z_p. One of the (semi)groups is additive, the other one multiplicative. This allows us to take advantage of both operations on matrices to diffuse information. We note that in our protocol, no power of any matrix or of any element of Z_p is ever exposed, so all standard attacks on Diffie-Hellman-like protocols (including Shor's quantum algorithm attack) are not applicable.


title: Subtask Analysis of Process Data Through a Predictive Model

abstract: Response process data collected from human-computer interactive items contain rich information about respondents' behavioral patterns and cognitive processes. Their irregular formats as well as their large sizes make standard statistical tools difficult to apply. This paper develops a computationally efficient method for exploratory analysis of such process data. The new approach segments a lengthy individual process into a sequence of short subprocesses to achieve complexity reduction, easy clustering and meaningful interpretation. Each subprocess is considered a subtask. The segmentation is based on sequential action predictability using a parsimonious predictive model combined with the Shannon entropy. Simulation studies are conducted to assess performance of the new methods. We use the process data from PIAAC 2012 to demonstrate how exploratory analysis of process data can be done with the new approach.


title: When the Differences in Frequency Domain are Compensated: Understanding and Defeating Modulated Replay Attacks on Automatic Speech Recognition

abstract: Automatic speech recognition (ASR) systems have been widely deployed in modern smart devices to provide convenient and diverse voice-controlled services. Since ASR systems are vulnerable to audio replay attacks that can spoof and mislead ASR systems, a number of defense systems have been proposed to identify replayed audio signals based on the speakers' unique acoustic features in the frequency domain. In this paper, we uncover a new type of replay attack called modulated replay attack, which can bypass the existing frequency domain based defense systems. The basic idea is to compensate for the frequency distortion of a given electronic speaker using an inverse filter that is customized to the speaker's transform characteristics. Our experiments on real smart devices confirm the modulated replay attacks can successfully escape the existing detection mechanisms that rely on identifying suspicious features in the frequency domain. To defeat modulated replay attacks, we design and implement a countermeasure named DualGuard. We discover and formally prove that no matter how the replay audio signals could be modulated, the replay attacks will either leave ringing artifacts in the time domain or cause spectrum distortion in the frequency domain. Therefore, by jointly checking suspicious features in both frequency and time domains, DualGuard can successfully detect various replay attacks including the modulated replay attacks. We implement a prototype of DualGuard on a popular voice interactive platform, ReSpeaker Core v2. The experimental results show DualGuard can achieve 98% accuracy on detecting modulated replay attacks.


title: Excited state R\'enyi entropy and subsystem distance in two-dimensional non-compact bosonic theory -- I. Single-particle states

abstract: We investigate the R\'enyi entropy of the excited states produced by the current and its derivatives in the two-dimensional free massless non-compact bosonic theory, which is a two-dimensional conformal field theory. We also study the subsystem Schatten distance between these states. The two-dimensional free massless non-compact bosonic theory is the continuum limit of the finite periodic gapless harmonic chains with the local interactions. We identify the excited states produced by current and its derivatives in the massless bosonic theory as the single-particle excited states in the gapless harmonic chain. We calculate analytically the second R\'enyi entropy and the second Schatten distance in the massless bosonic theory. We then use the wave functions of the excited states and calculate the second R\'enyi entropy and the second Schatten distance in the gapless limit of the harmonic chain, which match perfectly with the analytical results in the massless bosonic theory. We verify that in the large momentum limit the single-particle state R\'enyi entropy takes a universal form. We also show that in the limit of large momenta and large momentum difference the subsystem Schatten distance takes a universal form but it is replaced by a corrected form when the momentum difference is small. Finally we also comment on the mutual R\'enyi entropy of two disjoint intervals in the excited states of the two-dimensional free non-compact bosonic theory.


title: Locally Homogeneous Non-gradient Quasi-Einstein 3-Manifolds

abstract: In this paper, we classify the compact locally homogeneous non-gradient $m$-quasi Einstein 3-manifolds. Along the way, we prove that given a compact quotient of a Lie group of any dimension that is $m$-quasi Einstein, the potential vector field $X$ must be left invariant and Killing. We also classify the nontrivial $m$-quasi Einstein metrics that are a compact quotient of be the product of two Einstein metrics. We also show that $S^1$ is the only compact manifold of any dimension which admits a metric which is nontrivially $m$-quasi Einstein and Einstein.


title: Simultaneous Observation of Carrier-Specific Redistribution and Coherent Lattice Dynamics in 2H-MoTe$_{2}$ with Femtosecond Core-Level Spectroscopy

abstract: We employ few-femtosecond extreme ultraviolet (XUV) transient absorption spectroscopy to reveal simultaneously the intra- and interband carrier relaxation and the light-induced structural dynamics in nanoscale thin films of layered 2H-MoTe$_{2}$ semiconductor. By interrogating the valence electronic structure via localized Te 4$\textit{d}$ (39-46 eV) and Mo 4$\textit{p}$ (35-38 eV) core levels, the relaxation of the photoexcited hole distribution is directly observed in real time. We obtain hole thermalization and cooling times of 15$\pm$5 fs and 380$\pm$90 fs, respectively, and an electron-hole recombination time of 1.5$\pm$0.1 ps. Furthermore, excitations of coherent out-of-plane A$_{1g}$ (5.1 THz) and in-plane E$_{1g}$ (3.7 THz) lattice vibrations are visualized through oscillations in the XUV absorption spectra. By comparison to Bethe-Salpeter equation simulations, the spectral changes are mapped to real-space excited-state displacements of the lattice along the dominant A$_{1g}$ coordinate. By directly and simultaneously probing the excited carrier distribution dynamics and accompanying femtosecond lattice displacement in 2H-MoTe$_{2}$ within a single experiment, our work provides a benchmark for understanding the interplay between electronic and structural dynamics in photoexcited nanomaterials.


title: Encodable: Configurable Grammar for Visualization Components

abstract: There are so many libraries of visualization components nowadays with their APIs often different from one another. Could these components be more similar, both in terms of the APIs and common functionalities? For someone who is developing a new visualization component, how should the API look like? This work drew inspiration from visualization grammar, decoupled the grammar from its rendering engine and adapted it into a configurable grammar for individual components called Encodable. Encodable helps component authors define grammar for their components, and parse encoding specifications from users into utility functions for the implementation. This paper explains the grammar design and demonstrates how to build components with it.


title: Chemical Guidance in the Search for Past and Extant Life on Mars

abstract: NASA should design missions to Mars for the purpose of generating "Aha!" discoveries to jolt scientists contemplating the molecular origins of life. These missions should be designed with an understanding of the privileged chemistry that likely created RNA prebiotically on Earth, and universal chemical principles that constrain the structure of Darwinian molecules generally.


title: On Throughput Improvement using Immediate Re-transmission in Grant-Free Random Access with Massive MIMO

abstract: To support machine-type communication (MTC), massive multiple-input multiple-output (MIMO) has been considered for grant-free random access. In general, the performance of grant-free random access with massive MIMO is limited by the number of preambles and the number of active devices. In particular, when there are a number of active devices transmitting data packets simultaneously, the signal-to-interference-plus-noise ratio (SINR) cannot be high enough for successful decoding. In this paper, in order to improve performance, we consider immediate re-transmissions for an active device that has a low SINR although it does not experience preamble collision to exploit re-transmission diversity (RTD) gain. To see the performance of the proposed approach, we perform throughput analysis with certain approximations and assumption. Since the proposed approach can be unstable due to immediate re-transmissions, conditions for stable systems are also studied. Simulations are carried out and it is shown that analysis results reasonably match simulation results.


title: Conditional Constrained Graph Variational Autoencoders for Molecule Design

abstract: In recent years, deep generative models for graphs have been used to generate new molecules. These models have produced good results, leading to several proposals in the literature. However, these models may have troubles learning some of the complex laws governing the chemical world. In this work, we explore the usage of the histogram of atom valences to drive the generation of molecules in such models. We present Conditional Constrained Graph Variational Autoencoder (CCGVAE), a model that implements this key-idea in a state-of-the-art model, and shows improved results on several evaluation metrics on two commonly adopted datasets for molecule generation.


title: SPAN: Spatial Pyramid Attention Network forImage Manipulation Localization

abstract: We present a novel framework, Spatial Pyramid Attention Network (SPAN) for detection and localization of multiple types of image manipulations. The proposed architecture efficiently and effectively models the relationship between image patches at multiple scales by constructing a pyramid of local self-attention blocks. The design includes a novel position projection to encode the spatial positions of the patches. SPAN is trained on a generic, synthetic dataset but can also be fine tuned for specific datasets; The proposed method shows significant gains in performance on standard datasets over previous state-of-the-art methods.


title: Performance Analysis and Non-Quadratic Lyapunov Functions for Linear Time-Varying Systems

abstract: Performance analysis for linear time-invariant (LTI) systems has been closely tied to quadratic Lyapunov functions ever since it was shown that LTI system stability is equivalent to the existence of such a Lyapunov function. Some metrics for LTI systems, however, have resisted treatment via means of quadratic Lyapunov functions. Among these, point-wise-in-time metrics, such as peak norms, are not captured accurately using these techniques, and this shortcoming has prevented the development of tools to analyze system behavior by means other than e.g. time-domain simulations. This work demonstrates how the more general class of homogeneous polynomial Lyapunov functions can be used to approximate point-wise-in-time behavior for LTI systems with greater accuracy, and we extend this to the case of linear time-varying (LTV) systems as well. Our findings rely on the recent observation that the search for homogeneous polynomial Lyapunov functions for LTV systems can be recast as a search for quadratic Lyapunov functions for a related hierarchy of time-varying Lyapunov differential equations; thus, performance guarantees for LTV systems are attainable without heavy computation. Numerous examples are provided to demonstrate the findings of this work.


title: The cosmological constant and the use of cutoffs

abstract: Of the contributions to the cosmological constant, zero-point energy and self energy contributions scale as $\Lambda^4$ where $\Lambda$ is an ultraviolet cutoff used to regulate the calculations. I show that such contributions vanish when calculated in perturbation theory. This demonstration uses a little-known modification to perturbation theory found by Honerkamp and Meetz and by Gerstein, Jackiw, Lee and Weinberg which comes into play when using cutoffs and interactions with multiple derivatives, as found in chiral theories and gravity. In a path integral treatment, the new interaction arises from the path integral measure. This reduces the sensitivity of the cosmological constant to the high energy cutoff significantly, although it does not resolve the cosmological constant problem. The feature removes one of the common motivations for supersymmetry. It also calls into question some of the results of the Asymptotic Safety program. Covariance and quadratic cutoff dependence are also briefly discussed.


title: Evaluating Catchment Models as Multiple Working Hypotheses: on the Role of Error Metrics, Parameter Sampling, Model Structure, and Data Information Content

abstract: To evaluate models as hypotheses, we developed the method of Flux Mapping to construct a hypothesis space based on dominant runoff generating mechanisms. Acceptable model runs, defined as total simulated flow with similar (and minimal) model error, are mapped to the hypothesis space given their simulated runoff components. In each modeling case, the hypothesis space is the result of an interplay of factors: model structure and parameterization, chosen error metric, and data information content. The aim of this study is to disentangle the role of each factor in model evaluation. We used two model structures (SACRAMENTO and SIMHYD), two parameter sampling approaches (Latin Hypercube Sampling of the parameter space and guided-search of the solution space), three widely used error metrics (Nash-Sutcliffe Efficiency - NSE, Kling-Gupta Efficiency skill score - KGEss, and Willmott refined Index of Agreement - WIA), and hydrological data from a large sample of Australian catchments. First, we characterized how the three error metrics behave under different error types and magnitudes independent of any modeling. We then conducted a series of controlled experiments to unpack the role of each factor in runoff generation hypotheses. We show that KGEss is a more reliable metric compared to NSE and WIA for model evaluation. We further demonstrate that only changing the error metric -- while other factors remain constant -- can change the model solution space and hence vary model performance, parameter sampling sufficiency, and or the flux map. We show how unreliable error metrics and insufficient parameter sampling impair model-based inferences, particularly runoff generation hypotheses.

categories: stat.ME physics.data-an physics.geo-ph stat.AP


title: On the first $\delta$ Sct--roAp hybrid pulsator and the stability of p and g modes in chemically peculiar A/F stars

abstract: Strong magnetic fields in chemically peculiar A-type (Ap) stars typically suppress low-overtone pressure modes (p modes) but allow high-overtone p modes to be driven. KIC 11296437 is the first star to show both. We obtained and analysed a Subaru spectrum, from which we show that KIC 11296437 has abundances similar to other magnetic Ap stars, and we estimate a mean magnetic field modulus of $2.8\pm0.5$ kG. The same spectrum rules out a double-lined spectroscopic binary, and we use other techniques to rule out binarity over a wide parameter space, so the two pulsation types originate in one $\delta$ Sct--roAp hybrid pulsator. We construct stellar models depleted in helium and demonstrate that helium settling is second to magnetic damping in suppressing low-overtone p modes in Ap stars. We compute the magnetic damping effect for selected p and g modes, and find that modes with frequencies similar to the fundamental mode are driven for polar field strengths $\lesssim4$ kG, while other low-overtone p modes are driven for polar field strengths up to $\sim$1.5 kG. We find that the high-order g modes commonly observed in $\gamma$ Dor stars are heavily damped by polar fields stronger than 1--4 kG, with the damping being stronger for higher radial orders. We therefore explain the observation that no magnetic Ap stars have been observed as $\gamma$ Dor stars. We use our helium-depleted models to calculate the $\delta$ Sct instability strip for metallic-lined A (Am) stars, and find that driving from a Rosseland mean opacity bump at $\sim$$5\times10^4$ K caused by the discontinuous H-ionization edge in bound-free opacity explains the observation of $\delta$ Sct pulsations in Am stars.


title: On the decay in $W^{1,\infty}$ for the 1D semilinear damped wave equation on a bounded domain

abstract: In this paper we study a semilinear wave equation with nonlinear, time-dependent damping in one space dimension. For this problem, we prove a well-posedness result in $W^{1,\infty}$ in the space-time domain $(0,1)\times [0,+\infty)$. Then we address the problem of the time-asymptotic stability of the zero solution and show that, under appropriate conditions, the solution decays to zero at an exponential rate in the space $W^{1,\infty}$. The proofs are based on the analysis of the corresponding semilinear system for the first order derivatives, for which we show a contractive property of the invariant domain.


title: On Stars in Caterpillars and Lobsters

abstract: The family of all $k$-independent sets of a graph containing a fixed vertex $v$ is called a {star} and $v$ is called its center. Stars are interesting for their relation to Erd\"{o}s-Ko-Rado graphs. Hurlbert and Kamat conjectured that in trees the largest stars are centered in leafs. This conjecture was disproven independently by Baber, Borg, and Feghali, Johnson, and Thomas. In this paper we introduce a tool to bound the size of stars centered at certain vertices by stars centered at leafs. We use this tool to show that caterpillars and sunlet graphs satisfy Hurlbert and Kamat's conjecture, and to show that the centers of the largest stars in lobsters are either leafs or spinal vertices of degree 2.


title: Cosmicflows-4: The Catalog of ~10000 Tully-Fisher Distances

abstract: We present the distances of 9792 spiral galaxies lying within 15,000 km/s using the relation between luminosity and rotation rate of spiral galaxies. The sample is dominantly, but not exclusively, drawn from galaxies detected in the course of the ALFALFA HI survey with the Arecibo Telescope. Relations between \hi line widths and luminosity are calibrated at SDSS u, g, r, i, z bands and WISE W1 and W2 bands. By exploiting secondary parameters, particularly color indices, we address discrepancies between measured distances at different wave bands with unprecedented detail. We provide a catalog that includes reduced kinematic, photometric, and inclination parameters. We also describe a machine learning algorithm, based on the random forest technique that predicts the dust attenuation in spirals lacking infrared photometry. We determine a Hubble Constant value of H0 = 75.1+-0.2 (stat.), with potential systematics up to +-3 km/s/Mpc.


title: Disaggregating Customer-level Behind-the-Meter PV Generation Using Smart Meter Data and Solar Exemplars

abstract: Customer-level rooftop photovoltaic (PV) has been widely integrated into distribution systems. In most cases, PVs are installed behind-the-meter (BTM), and only the net demand is recorded. Therefore, the native demand and PV generation are unknown to utilities. Separating native demand and solar generation from net demand is critical for improving grid-edge observability. In this paper, a novel approach is proposed for disaggregating customer-level BTM PV generation using low-resolution but widely available hourly smart meter data. The proposed approach exploits the strong correlation between monthly nocturnal and diurnal native demands and the high similarity among PV generation profiles. First, a joint probability density function (PDF) of monthly nocturnal and diurnal native demands is constructed for customers without PVs, using Gaussian mixture modeling (GMM). Deviation from the constructed PDF is utilized to probabilistically assess the monthly solar generation of customers with PVs. Then, to identify hourly BTM solar generation for these customers, their estimated monthly solar generation is decomposed into an hourly timescale; to do this, we have proposed a maximum likelihood estimation (MLE)-based technique that utilizes hourly typical solar exemplars. Leveraging the strong monthly native demand correlation and high PV generation similarity enhances our approach's robustness against the volatility of customers' hourly load and enables highly accurate disaggregation. The proposed approach has been verified using real native demand and PV generation data.


title: On Open and Strong-Scaling Tools for Atom Probe Crystallography: High-Throughput Methods for Indexing Crystal Structure and Orientation

abstract: Volumetric crystal structure indexing and orientation mapping are key data processing steps for virtually any quantitative study of spatial correlations between the local chemistry and the microstructure of a material. For electron and X-ray diffraction methods it is possible to develop indexing tools which compare measured and analytically computed patterns to decode the structure and relative orientation within local regions of interest. Consequently, a number of numerically efficient and automated software tools exist to solve the above characterisation tasks. For atom probe tomography (APT) experiments, however, the strategy of making comparisons between measured and analytically computed patterns is less robust because many APT datasets may contain substantial noise. Given that general enough predictive models for such noise remain elusive, crystallography tools for APT face several limitations: Their robustness to noise, and therefore, their capability to identify and distinguish different crystal structures and orientation is limited. In addition, the tools are sequential and demand substantial manual interaction. In combination, this makes robust uncertainty quantifying with automated high-throughput studies of the latent crystallographic information a difficult task with APT data. To improve the situation, we review the existent methods and discuss how they link to those in the diffraction communities. With this we modify some of the APT methods to yield more robust descriptors of the atomic arrangement. We report how this enables the development of an open-source software tool for strong-scaling and automated identifying of crystal structure and mapping crystal orientation in nanocrystalline APT datasets with multiple phases.

categories: physics.comp-ph cond-mat.mtrl-sci cs.AI cs.CV


title: Error estimates for a pointwise tracking optimal control problem of a semilinear elliptic equation

abstract: We consider a pointwise tracking optimal control problem for a semilinear elliptic partial differential equation. We derive the existence of optimal solutions and analyze first and, necessary and sufficient, second order optimality conditions. We devise two strategies of discretization to approximate a solution of the optimal control problem: a semidiscrete scheme where the control variable is not discretized -- the so-called variational discretization approach -- and a fully discrete scheme where the control variable is discretized with piecewise constant functions. For both solution techniques, we analyze convergence properties of discretizations and derive error estimates.


title: Fermions in the presence of topological structures under geometric constrictions

abstract: In this work we study modifications of the spectrum of fermions interacting with kinklike structures in two-dimensional spacetime. We consider the Yukawa coupling between fermions and scalar fields that engender nontrivial internal structure and investigate how the fermion spectra change in terms of the parameters that control the kinklike configuration and the Yukawa coupling. We consider models that allow the internal structure of the kinklike solution to respond to the presence of a geometrical constriction, and show the fermion spectra may also appear directly affected by the constriction. The main results are of current interest and may be used to propose the construction of electronic devices capable of engendering new effects at the nanometric scale.


title: A Deontic Logic Analysis of Autonomous Systems' Safety

abstract: We consider the pressing question of how to model, verify, and ensure that autonomous systems meet certain \textit{obligations} (like the obligation to respect traffic laws), and refrain from impermissible behavior (like recklessly changing lanes). Temporal logics are heavily used in autonomous system design; however, as we illustrate here, temporal (alethic) logics alone are inappropriate for reasoning about obligations of autonomous systems. This paper proposes the use of Dominance Act Utilitarianism (DAU), a deontic logic of agency, to encode and reason about obligations of autonomous systems. We use DAU to analyze Intel's Responsibility-Sensitive Safety (RSS) proposal as a real-world case study. We demonstrate that DAU can express well-posed RSS rules, formally derive undesirable consequences of these rules, illustrate how DAU could help design systems that have specific obligations, and how to model-check DAU obligations.


title: Non-asymptotic Identification of Linear Dynamical Systems Using Multiple Trajectories

abstract: This paper considers the problem of linear time-invariant (LTI) system identification using input/output data. Recent work has provided non-asymptotic results on partially observed LTI system identification using a single trajectory but is only suitable for stable systems. We provide finite-time analysis for learning Markov parameters based on the ordinary least-squares (OLS) estimator using multiple trajectories, which covers both stable and unstable systems. For unstable systems, our results suggest that the Markov parameters are harder to estimate in the presence of process noise. Without process noise, our upper bound on the estimation error is independent of the spectral radius of system dynamics with high probability. These two features are different from fully observed LTI systems for which recent work has shown that unstable systems with a bigger spectral radius are easier to estimate. Extensive numerical experiments demonstrate the performance of our OLS estimator.


title: Control Framework for a Hybrid-steel Bridge Inspection Robot

abstract: Autonomous navigation of steel bridge inspection robots is essential for proper maintenance. The majority of existing robotic solutions for bridge inspection require human intervention to assist in the control and navigation. In this paper, a control system framework has been proposed for a previously designed ARA robot [1], which facilitates autonomous real-time navigation and minimizes human involvement. The mechanical design and control framework of ARA robot enables two different configurations, namely the mobile and inch-worm transformation. In addition, a switching control was developed with 3D point clouds of steel surfaces as the input which allows the robot to switch between mobile and inch-worm transformation. The surface availability algorithm (considers plane, area, and height) of the switching control enables the robot to perform inch-worm jumps autonomously. Themobiletransformationallows the robot to move on continuous steel surfaces and perform visual inspection of steel bridge structures. Practical experiments on actual steel bridge structures highlight the effective performance of ARA robot with the proposed control framework for autonomous navigation during a visual inspection of steel bridges.


title: Radius, Girth and Minimum Degree

abstract: Given a connected graph $G$ on $n$ vertices, with minimum degree $\delta\geq 2$ and girth at least $g \geq 4$, what is the maximum radius $r$ this graph can have? Erd\H{o}s, Pach, Pollack and Tuza established in the triangle-free case ($g=4$) that $r \leq \frac{n-2}{\delta}+12$, and noted that up to the value of the additive constant, this is tight. We determine the exact value for the triangle-free case. For higher $g$ little is known. We settle the order of $r$ for $g=6,8,12$ and prove an upper bound to the order for general even $g$. Finally, we show that proving the corresponding lower bound for general even $g$ is equivalent to the Erd\H{o}s girth conjecture.


title: Non-triviality in a totally asymmetric one-dimensional Boolean percolation model on a half-line

abstract: It is well known that there are two regimes in a standard one-dimensional Boolean percolation model: either the entire space is covered a.s., or the covered volume fraction is strictly less than one. The aim of this work is to demonstrate that there is a third possibility in a Boolean model with totally asymmetric grains on a half-line: a.s. there is no unbounded component, but the covered volume fraction is one. An explicit condition is given characterizing the existence of an unbounded occupied component.


title: Bidirectional Attention Network for Monocular Depth Estimation

abstract: In this paper, we propose a Bidirectional Attention Network (BANet), an end-to-end framework for monocular depth estimation (MDE) that addresses the limitation of effectively integrating local and global information in convolutional neural networks. The structure of this mechanism derives from a strong conceptual foundation of neural machine translation, and presents a light-weight mechanism for adaptive control of computation similar to the dynamic nature of recurrent neural networks. We introduce bidirectional attention modules that utilize the feed-forward feature maps and incorporate the global context to filter out ambiguity. Extensive experiments reveal the high degree of capability of this bidirectional attention model over feed-forward baselines and other state-of-the-art methods for monocular depth estimation on two challenging datasets -- KITTI and DIODE. We show that our proposed approach either outperforms or performs at least on a par with the state-of-the-art monocular depth estimation methods with less memory and computational complexity.


title: Computational screening of repurposed drugs and natural products against SARS-Cov-2 main protease (Mpro) as potential COVID-19 therapies

abstract: There remains an urgent need to identify existing drugs that might be suitable for treating patients suffering from COVID-19 infection. Drugs rarely act at a single molecular target, with off target effects often being responsible for undesirable side effects and sometimes, beneficial synergy between targets for a specific illness. Off target activities have also led to blockbuster drugs in some cases, e.g. Viagra for erectile dysfunction and Minoxidil for male pattern hair loss. Drugs already in use or in clinical trials plus approved natural products constitute a rich resource for discovery of therapeutic agents that can be repurposed for existing and new conditions, based on the rationale that they have already been assessed for safety in man. A key question then is how to rapidly and efficiently screen such compounds for activity against new pandemic pathogens such as COVID-19. Here we show how a fast and robust computational process can be used to screen large libraries of drugs and natural compounds to identify those that may inhibit the main protease of SARS-Cov-2 (3CL pro, Mpro). We show how the resulting shortlist of candidates with strongest binding affinities is highly enriched in compounds that have been independently identified as potential antivirals against COVID-19. The top candidates also include a substantial number of drugs and natural products not previously identified as having potential COVID-19 activity, thereby providing additional targets for experimental validation. This in silico screening pipeline may also be useful for repurposing of existing drugs and discovery of new drug candidates against other medically important pathogens and for use in future pandemics.


title: The LISE package: solvers for static and time-dependent superfluid local density approximation equations in three dimensions

abstract: Nuclear implementation of the density functional theory (DFT) is at present the only microscopic framework applicable to the whole nuclear landscape. The extension of DFT to superfluid systems in the spirit of the Kohn-Sham approach, the superfluid local density approximation (SLDA) and its extension to time-dependent situations, time-dependent superfluid local density approximation (TDSLDA), have been extensively used to describe various static and dynamical problems in nuclear physics, neutron star crust, and cold atom systems. In this paper, we present the codes that solve the static and time-dependent SLDA equations in three-dimensional coordinate space without any symmetry restriction. These codes are fully parallelized with the message passing interface (MPI) library and take advantage of graphic processing units (GPU) for accelerating execution. The dynamic codes have checkpoint/restart capabilities and for initial conditions one can use any generalized Slater determinant type of wave function. The code can describe a large number of physical problems: nuclear fission, collisions of heavy ions, the interaction of quantized vortices with nuclei in the nuclear star crust, excitation of superfluid fermion systems by time dependent external fields, quantum shock waves, domain wall generation and propagation, the dynamics of the Anderson-Bogoliubov-Higgs mode, dynamics of fragmented condensates, vortex rings dynamics, generation and dynamics of quantized vortices, their crossing and recombinations and the incipient phases of quantum turbulence.


title: Single-photon switch controlled by a qubit embedded in an engineered electromagnetic environment

abstract: A single-photon switch is an important element for the building of scalable quantum networks. In this paper, we propose a feasible scheme for efficient single-photon switching. The proposed switch is controlled by a state of a qubit formed by the pair of the lowest levels of a three-level system (qutrit) coupled to a resonator. This resonator-qutrit system comprises a switching unit of the considered setup. For suppression of the Purcell relaxation of the control qubit, the switching unit is embedded into a coupled-resonator array serving as an engineered electromagnetic environment with a band gap on a qubit transition frequency. We discuss the possible implementation of the considered single-photon switch on the microwave circuit QED architecture. We demonstrate that high switching contrasts can be attained for the parameters achievable for the state-of-the-art superconducting circuit QED setups.


title: SM antideuteron background to indirect dark matter signals in galactic cosmic rays

abstract: Antideuteron production cross-sections estimated using EPOS-LHC with a coalescence afterburner, tuned to reproduce published experimental data over a wide range of energy were used here as input to the galactic propagator code GALPROP, validated with comparing to existing proton, helium fluxes as well as boron-to-carbon ratio data. The resulting near-Earth antideuteron flux, including solar modulation, is compared to previous estimates. An overall factor of two increments in the antideuteron flux is predicted, the origin of which is also discussed. However, this standard model source of antideuteron background still lies well below the AMS-02, and the expected GAPS, sensitivities, as well as the fluxes predicted by several dark matter models.


title: TensorDash: Exploiting Sparsity to Accelerate Deep Neural Network Training and Inference

abstract: TensorDash is a hardware level technique for enabling data-parallel MAC units to take advantage of sparsity in their input operand streams. When used to compose a hardware accelerator for deep learning, TensorDash can speedup the training process while also increasing energy efficiency. TensorDash combines a low-cost, sparse input operand interconnect comprising an 8-input multiplexer per multiplier input, with an area-efficient hardware scheduler. While the interconnect allows a very limited set of movements per operand, the scheduler can effectively extract sparsity when it is present in the activations, weights or gradients of neural networks. Over a wide set of models covering various applications, TensorDash accelerates the training process by $1.95{\times}$ while being $1.89\times$ more energy-efficient, $1.6\times$ more energy efficient when taking on-chip and off-chip memory accesses into account. While TensorDash works with any datatype, we demonstrate it with both single-precision floating-point units and bfloat16.


title: Iris Liveness Detection Competition (LivDet-Iris) -- The 2020 Edition

abstract: Launched in 2013, LivDet-Iris is an international competition series open to academia and industry with the aim to assess and report advances in iris Presentation Attack Detection (PAD). This paper presents results from the fourth competition of the series: LivDet-Iris 2020. This year's competition introduced several novel elements: (a) incorporated new types of attacks (samples displayed on a screen, cadaver eyes and prosthetic eyes), (b) initiated LivDet-Iris as an on-going effort, with a testing protocol available now to everyone via the Biometrics Evaluation and Testing (BEAT)(https://www.idiap.ch/software/beat/) open-source platform to facilitate reproducibility and benchmarking of new algorithms continuously, and (c) performance comparison of the submitted entries with three baseline methods (offered by the University of Notre Dame and Michigan State University), and three open-source iris PAD methods available in the public domain. The best performing entry to the competition reported a weighted average APCER of 59.10\% and a BPCER of 0.46\% over all five attack types. This paper serves as the latest evaluation of iris PAD on a large spectrum of presentation attack instruments.


title: Visualizing information on watch faces: A survey with smartwatch users

abstract: People increasingly wear smartwatches that can track a wide variety of data. However, it is currently unknown which data people consume and how it is visualized. To better ground research on smartwatch visualization, it is important to understand the current use of these representation types on smartwatches, and to identify missed visualization opportunities. We present the findings of a survey with 237 smartwatch wearers, and assess the types of data and representations commonly displayed on watch faces. We found a predominant display of health & fitness data, with icons accompanied by text being the most frequent representation type. Combining these results with a further analysis of online searches of watch faces and the data tracked on smartwatches that are not commonly visualized, we discuss opportunities for visualization research.


title: Text Modular Networks: Learning to Decompose Tasks in the Language of Existing Models

abstract: We propose a general framework called Text Modular Networks(TMNs) for building interpretable systems that learn to solve complex tasks by decomposing them into simpler ones solvable by existing models. To ensure solvability of simpler tasks, TMNs learn the textual input-output behavior (i.e., language) of existing models through their datasets. This differs from prior decomposition-based approaches which, besides being designed specifically for each complex task, produce decompositions independent of existing sub-models. Specifically, we focus on Question Answering (QA) and show how to train a next-question generator to sequentially produce sub-questions targeting appropriate sub-models, without additional human annotation. These sub-questions and answers provide a faithful natural language explanation of the model's reasoning. We use this framework to build ModularQA, a system that can answer multi-hop reasoning questions by decomposing them into sub-questions answerable by a neural factoid single-span QA model and a symbolic calculator. Our experiments show that ModularQA is more versatile than existing explainable systems for DROP and HotpotQA datasets, is more robust than state-of-the-art blackbox (uninterpretable) systems, and generates more understandable and trustworthy explanations compared to prior work.


title: Continuous slowing of a gadolinium atomic beam

abstract: The article presents the development of a new and innovative experimental method to fully characterize a solenoidal "spin-flip" Zeeman slower (ZS) using a Quartz Crystal $\mu$-balance (QCM) as a kinetic energy sensor. In this experiment, we focus a 447.1 nm laser into a counter-propagating beam of gadolinium (Gd) atoms in order to drive the dipole transition between ground $^{9}$D$^{0}_{2}$ state and $^{9}$D$_{3}$ excited state. The changes in the velocity of the beam were measured using a QCM during this process, as a novel and alternative method to characterize the efficiency of a 1 m-long spin-flip Zeeman slower. The QCM, normally used in solid-state physics, is continuously and carefully monitored to determine the change in its natural frequency of oscillation. These changes reveal a direct relation with changes in the deposition rate and the momentum exchanged between the QCM and Gd atoms. Hence, in terms of ultracold atom physics, it might be used to study the time-evolution of the velocity distribution of the atoms during the cooling process. By this method, we obtain a maximum atom average velocity reduction of (43.5 $\pm$ 6.4)$\%$ produced by our apparatus. Moreover, we estimate an experimental lifetime of $\tau_{e}$ = 8.2 ns for the used electronic transition, and then we compared it with the reported lifetime for 443.06 nm and 451.96 nm electronic transitions of Gd. These results confirm that the QCM offers an accessible and simple solution to take into account for laser cooling experiments. Therefore, a novel and innovative technique can be available for future experiments.


title: Covariant propagator and chiral power counting

abstract: Some one-loop diagrams with one and two external baryons/nucleons are revisited using covariant baryon propagators in chiral effective theory. We showed that it is enough to separate and subtract all the local terms that violate chiral power counting to recover chiral power counting, no need to introduce extra operations. The structures of leading chiral corrections and IR enhancement or threshold effects are 'stable' or persist as long as covariant propagators are employed for all particles.


title: Constraining Einstein's Equivalence Principle With Multi-Wavelength Polarized astrophysical Sources

abstract: The observed time delays between photons with different circular polarizations from an astrophysical object provide a new, interesting way of testing the Einstein Equivalence Principle (EEP). In this paper, we constrain the EEP by considering both Shapiro time delay and Faraday rotation effects. We continue to search for astronomical sources that are suitable for testing the EEP accuracy, and obtain 60 extragalactic radio sources with multi-wavelength polarization angles in three different radio bands (20, 8.6, and 4.8 GHz) and 29 brightest stars within our own Milky Way galaxy with multi-colour linear polarimetric data in five optical bands ($UBVRI$). We apply the Metropolis-Hastings Markov Chain to simulate the fit parameters. The final results show that the values of the parameterized post-Newtonian parameter $\gamma$ discrepancy ($\Delta \gamma_{p}$) are constrained to be in the range of $10^{-26}-10^{-23}$ for 60 radio sources and in the range of $10^{-23}-10^{-20}$ for 29 optical polarization stars. Compared to previous EEP tests that based on the single polarization measurement in the gamma-ray band, our results have profound superiority that nearly a few tens of astrophysical sources with multi-wavelength polarization observations commonly in the optical and radio bands are available. It ensures that these sources can give more significantly robust bounds on the EEP. Although the presented method is straightforward, the resulting constraints on the EEP should be taken as upper limits as other more complex astrophysical effects affecting a polarization rotation are hardly considered.


title: Turning Machines: a simple algorithmic model for molecular robotics

abstract: Molecular robotics is challenging, so it seems best to keep it simple. We consider an abstract molecular robotics model based on simple folding instructions that execute asynchronously. Turning Machines are a simple 1D to 2D folding model, also easily generalisable to 2D to 3D folding. A Turning Machine starts out as a line of connected monomers in the discrete plane, each with an associated turning number. A monomer turns relative to its neighbours, executing a unit-distance translation that drags other monomers along with it, and through collective motion the initial set of monomers eventually folds into a programmed shape. We provide a suite of tools for reasoning about Turning Machines by fully characterising their ability to execute line rotations: executing an almost-full line rotation of $5\pi/3$ radians is possible, yet a full $2\pi$ rotation is impossible. Furthermore, line rotations up to $5\pi/3$ are executed efficiently, in $O(\log n)$ expected time in our continuous time Markov chain time model. We then show that such line-rotations represent a fundamental primitive in the model, by using them to efficiently and asynchronously fold shapes. In particular, arbitrarily large zig-zag-rastered squares and zig-zag paths are foldable, as are $y$-monotone shapes albeit with error (bounded by perimeter length). Finally, we give shapes that despite having paths that traverse all their points, are in fact impossible to fold, as well as techniques for folding certain classes of (scaled) shapes without error. Our approach relies on careful geometric-based analyses of the feats possible and impossible by a very simple robotic system, and pushes conceptional hardness towards mathematical analysis and away from molecular implementation.


title: Point estimates, Simpson's paradox and nonergodicity in biological sciences

abstract: Modern biomedical, behavioral and psychological inference about cause-effect relationships respects an ergodic assumption, that is, that mean response of representative samples allow predictions about individual members of those samples. Recent empirical evidence in all of the same fields indicates systematic violations of the ergodic assumption. Indeed, violation of ergodicity in biomedical, behavioral and psychological causes is precisely the inspiration behind our research inquiry. Here, we review the long term costs to scientific progress in these domains and a practical way forward. Specifically, we advocate the use of statistical measures that can themselves encode the degree and type of non-ergodicity in measurements. Taking such steps will lead to a paradigm shift, allowing researchers to investigate the nonstationary, far-from-equilibrium processes that characterize the creativity and emergence of biological and psychological behavior.


title: Properties of f-divergences and f-GAN training

abstract: In this technical report we describe some properties of f-divergences and f-GAN training. We present an elementary derivation of the f-divergence lower bounds which form the basis of f-GAN training. We derive informative but perhaps underappreciated properties of f-divergences and f-GAN training, including a gradient matching property and the fact that all f-divergences agree up to an overall scale factor on the divergence between nearby distributions. We provide detailed expressions for computing various common f-divergences and their variational lower bounds. Finally, based on our reformulation, we slightly generalize f-GAN training in a way that may improve its stability.


title: Fully analog memristive circuits for optimization tasks: a comparison

abstract: We introduce a Lyapunov function for the dynamics of memristive circuits, and compare the effectiveness of memristors in minimizing the function to widely used optimization software. We study in particular three classes of problems which can be directly embedded in a circuit topology, and show that memristors effectively attempt at (quickly) extremizing these functionals.


title: The mean square radius of the neutron distribution and the skin thickness derived from electron scattering

abstract: The second-order moment of the nuclear charge density($R^2_c$) is dominated by the mean square radius(msr) of the point proton distribution($R_p^2$), while the fourth-order moment($Q^4_c$) depends on the msr of the point neutron one($R_n^2$) also. Moreover, $R^2_n$ is strongly correlated to $R^2_c$ in nuclear models. According to these facts, the linear relationship between various moments in the nuclear mean field models are investigated with use of the least squares method for $^{40}$Ca, $^{48}$Ca and $^{208}$Pb. From the intersection points of the obtained straight lines with those of the experimental values for $R^2_c$ and $Q^4_c$ determined through electron scattering, the values of $R_p$ and $R_n$ are estimated. Since relativistic and non-relativistic models provide different lines, the obtained values of $R_n$ and the skin thickness($R_n-R_p$) differ from each other in the two frameworks.


title: Distribution of peak heights modulo $k$ and double descents on $k$-Dyck paths

abstract: We show that the distribution of the number of peaks at height $i$ modulo $k$ in $k$-Dyck paths of a given length is independent of $i\in[0,k-1]$ and is the reversal of the distribution of the total number of peaks. Moreover, these statistics, together with the number of double descents, are jointly equidistributed with any of their permutations. We also generalize this result to generalized Motzkin paths and generalized ballot paths.


title: A Survey of Singular Value Decomposition Methods for Distributed Tall/Skinny Data

abstract: The Singular Value Decomposition (SVD) is one of the most important matrix factorizations, enjoying a wide variety of applications across numerous application domains. In statistics and data analysis, the common applications of SVD such as Principal Components Analysis (PCA) and linear regression. Usually these applications arise on data that has far more rows than columns, so-called "tall/skinny" matrices. In the big data analytics context, this may take the form of hundreds of millions to billions of rows with only a few hundred columns. There is a need, therefore, for fast, accurate, and scalable tall/skinny SVD implementations which can fully utilize modern computing resources. To that end, we present a survey of three different algorithms for computing the SVD for these kinds of tall/skinny data layouts using MPI for communication. We contextualize these with common big data analytics techniques, principally PCA. Finally, we present both CPU and GPU timing results from the Summit supercomputer, and discuss possible alternative approaches.


title: Multiparameter Quantum Estimation Theory in Quantum Gaussian states

abstract: Multiparameter quantum estimation theory aims to determine simultaneously the ultimate precision of all parameters contained in the state of a given quantum system. Determining this ultimate precision depends on the quantum Fisher information matrix (QFIM) which is essential to obtaining the quantum Cram\'er-Rao bound. This is the main motivation of this work which concerns the computation of the analytical expression of the QFIM. Inspired by the results reported in J. Phys. A 52, 035304 (2019), the general formalism of the multiparameter quantum estimation theory of quantum Gaussian states in terms of their first and second moments are given. We give the analytical formulas of right logarithmic derivative (RLD) and symmetric logarithmic derivative (SLD) operators. Then we derive the general expressions of the corresponding quantum Fisher information matrices. We also derive an explicit expression of the condition which ensures the saturation of the quantum Cram\'er-Rao bound in estimating several parameters. Finally, we examine some examples to clarify the use of our results


title: Depth Range Reduction for 3D Range Geometry Compression

abstract: Three-dimensional (3D) shape measurement devices and techniques are being rapidly adopted within a variety of industries and applications. As acquiring 3D range data becomes faster and more accurate it becomes more challenging to efficiently store, transmit, or stream this data. One prevailing approach to compressing 3D range data is to encode it within the color channels of regular 2D images. This paper presents a novel method for reducing the depth range of a 3D geometry such that it can be stored within a 2D image using lower encoding frequencies (or a fewer number of encoding periods). This allows for smaller compressed file sizes to be achieved without a proportional increase in reconstruction errors. Further, as the proposed method occurs prior to encoding, it is readily compatible with a variety of existing image-based 3D range geometry compression methods.


title: Monocular 3D Detection with Geometric Constraints Embedding and Semi-supervised Training

abstract: In this work, we propose a novel single-shot and keypoints-based framework for monocular 3D objects detection using only RGB images, called KM3D-Net. We design a fully convolutional model to predict object keypoints, dimension, and orientation, and then combine these estimations with perspective geometry constraints to compute position attribute. Further, we reformulate the geometric constraints as a differentiable version and embed it into the network to reduce running time while maintaining the consistency of model outputs in an end-to-end fashion. Benefiting from this simple structure, we then propose an effective semi-supervised training strategy for the setting where labeled training data is scarce. In this strategy, we enforce a consensus prediction of two shared-weights KM3D-Net for the same unlabeled image under different input augmentation conditions and network regularization. In particular, we unify the coordinate-dependent augmentations as the affine transformation for the differential recovering position of objects and propose a keypoints-dropout module for the network regularization. Our model only requires RGB images without synthetic data, instance segmentation, CAD model, or depth generator. Nevertheless, extensive experiments on the popular KITTI 3D detection dataset indicate that the KM3D-Net surpasses all previous state-of-the-art methods in both efficiency and accuracy by a large margin. And also, to the best of our knowledge, this is the first time that semi-supervised learning is applied in monocular 3D objects detection. We even surpass most of the previous fully supervised methods with only 13\% labeled data on KITTI.


title: Shell structure of $^{43}$S and collapse of the $N=28$ shell closure

abstract: The single-particle structure of the $N=27$ isotones provides insights into the shell evolution of neutron-rich nuclei from the doubly-magic $^{48}$Ca toward the drip line. $^{43}$S was studied employing the one-neutron knockout reaction from a radioactive $^{44}$S beam. Using a combination of prompt and delayed $\gamma$-ray spectroscopy the level structure of $^{43}$S was clarified. Momentum distributions were analyzed and allowed for spin and parity assignments. The deduced spectroscopic factors show that the $^{44}$S ground-state configuration has a strong intruder component. The results were confronted with shell model calculations using two effective interactions. General agreement was found between the calculations, but strong population of states originating from the removal of neutrons from the $2p_{3/2}$ orbital in the experiment indicates that the breakdown of the $N=28$ magic number is more rapid than the theoretical calculations suggest.


title: A short note on reconstruction variables in shock capturing schemes for magnetohydrodynamics

abstract: We propose a set of quick and easy approximate characteristic variables for higher-order reconstructions of shock capturing schemes for magnetohydrodynamics (MHD). Numerical experiments suggest that the reconstructions using the approximate characteristic variables are more robust than those using the conservative or primitive variables, while their computational efficiencies are comparable. The approximate characteristic variables are simple compared to the full characteristic variables for MHD, and can be a practical choice of reconstruction variables.


title: Summation of certain locally bilinear forms and its applications to the Fast Multipole Method

abstract: The Fast Multipole Method (FMM) reduces the computation of pairwise two-body interactions among $N$-particles to order $N$, whose computation cost should be of order $N^2$ by brute force. However, its implementation is somewhat complicated and requires a considerable amount of time to write the code. In this paper, I show a method that enables us to implement and write FMM algorithm code simply and briefly. FMM algorithm is composed of several steps. The main steps are Upward Pass and Downward Pass. Both the Upward Pass and Downward Pass include shift processes by which we move the centers of local expansions and multipole expansions. In this paper, I show a method that enables us to get rid of these processes.As a result of this simplification, the coding of FMM becomes much easier, and we can save considerable computation time. I compared the accuracy and time required to calculate potential fields with that of the existing FMM code.


title: Speaker Representation Learning using Global Context Guided Channel and Time-Frequency Transformations

abstract: In this study, we propose the global context guided channel and time-frequency transformations to model the long-range, non-local time-frequency dependencies and channel variances in speaker representations. We use the global context information to enhance important channels and recalibrate salient time-frequency locations by computing the similarity between the global context and local features. The proposed modules, together with a popular ResNet based model, are evaluated on the VoxCeleb1 dataset, which is a large scale speaker verification corpus collected in the wild. This lightweight block can be easily incorporated into a CNN model with little additional computational costs and effectively improves the speaker verification performance compared to the baseline ResNet-LDE model and the Squeeze&Excitation block by a large margin. Detailed ablation studies are also performed to analyze various factors that may impact the performance of the proposed modules. We find that by employing the proposed L2-tf-GTFC transformation block, the Equal Error Rate decreases from 4.56% to 3.07%, a relative 32.68% reduction, and a relative 27.28% improvement in terms of the DCF score. The results indicate that our proposed global context guided transformation modules can efficiently improve the learned speaker representations by achieving time-frequency and channel-wise feature recalibration.


title: An Explicit Upper Bound for $|\zeta(1+it)|$

abstract: In this paper we provide an explicit bound for $|\zeta(1+it)|$ in the form of $|\zeta(1+it)|\leq \min\left(\log t, \frac{1}{2}\log t+1.93, \frac{1}{5}\log t+44.02 \right)$. This improves on the current best-known explicit bound of $|\zeta(1+it)|\leq 62.6(\log t)^{2/3}$ up until $t$ of the magnitude $10^{10^7}$.


title: Strong interlayer coupling in two-dimensional PbSe with high thermoelectric performance

abstract: It was generally assumed that weak van der Waals interactions exist between neighboring layers in the two-dimensional group-IV chalcogenides. Using PbSe as a prototypal example, however, we find additional strong coupling between the Pb-Pb layers, as evidenced by detailed analysis of the differential charge density. The coupling resembles covalent-like bond and exhibits strong harmonicity around the equilibrium distance, which can be fine tuned to obviously reduce the phonon thermal conductivity but slightly change the electronic transport of PbSe. As a consequence, a maximum ZT value of 2.5 can be realized at 900 K for the p-type system. Our work offers an effective and feasible design strategy to enhance the thermoelectric performance of similar layered structures.