title: A 100 au-Wide Bipolar Rotating Shell Emanating From The HH 212Protostellar Disk: A Disk Wind?

abstract: HH 212 is a Class 0 protostellar system found to host a "hamburger"-shaped dusty disk with a rotating disk atmosphere and a collimated SiO jet at a distance of ~ 400 pc. Recently, a compact rotating outflow has been detected in SO and SO2 toward the center along the jet axis at ~ 52 au (0.13") resolution. Here we resolve the compact outflow into a small-scale wide-opening rotating outflow shell and a collimated jet, with the observations in the same S-bearing molecules at ~ 16 au (0.04") resolution. The collimated jet is aligned with the SiO jet, tracing the shock interactions in the jet. The wide-opening outflow shell is seen extending out from the inner disk around the SiO jet and has a width of ~ 100 au. It is not only expanding away from the center, but also rotating around the jet axis. The specific angular momentum of the outflow shell is ~ 40 au km/s. Simple modeling of the observed kinematics suggests that the rotating outflow shell can trace either a disk wind or disk material pushed away by an unseen wind from the inner disk or protostar. We also resolve the disk atmosphere in the same S-bearing molecules, confirming the Keplerian rotation there.


title: The Power Allocation Game on A Network: A Paradox

abstract: The well-known Braess paradox in congestion games states that adding an additional road to a transportation network may increase the total travel time, and consequently decrease the overall efficiency. Motivated by this, this paper presents a paradox in a similar spirit emerging from another distributed resource allocation game on networks, namely the power allocation game between countries developed in \cite{allocation}. The paradox is that by having additional friends may actually decrease a country's total welfare in equilibrium. Conditions for this paradox to occur as well as some price of anarchy results are also derived.


title: Non-Axisymmetric Line Driven Disc Winds II - Full Velocity Gradient

abstract: We study non-axisymetric features of 3D line driven winds in the Sobolev approximation, where the optical depth is calculated using the full velocity gradient. We find that non-axisymmetric density features, so called clumps, form primarily at the base of the wind on super-Sobolev length scales. The density of clumps differs by a factor of $\sim 3$ from the azimuthal average, the magnitude of their velocity dispersion is comparable to the flow velocity and they produce $\sim 20\%$ variations in the column density. Clumps may be observable because differences in density produce enhancements in emission and absorption profiles or through their velocity dispersion which enhances line broadening.


title: Faster Distributed Shortest Path Approximations via Shortcuts

abstract: A long series of recent results and breakthroughs have led to faster and better distributed approximation algorithms for single source shortest paths (SSSP) and related problems in the CONGEST model. The runtime of all these algorithms, however, is $\tilde{\Omega}(\sqrt{n})$, regardless of the network topology, even on nice networks with a (poly)logarithmic network diameter $D$. While this is known to be necessary for some pathological networks, most topologies of interest are arguably not of this type. We give the first distributed approximation algorithms for shortest paths problems that adjust to the topology they are run on, thus achieving significantly faster running times on many topologies of interest. The running time of our algorithms depends on and is close to $Q$, where $Q$ is the quality of the best shortcut that exists for the given topology. While $Q = \tilde{\Theta}(\sqrt{n} + D)$ for pathological worst-case topologies, many topologies of interest have $Q = \tilde{\Theta}(D)$, which results in near instance optimal running times for our algorithm, given the trivial $\Omega(D)$ lower bound. The problems we consider are as follows: (1) an approximate shortest path tree and SSSP distances, (2) a polylogarithmic size distance label for every node such that from the labels of any two nodes alone one can determine their distance (approximately), and (3) an (approximately) optimal flow for the transshipment problem. Our algorithms have a tunable tradeoff between running time and approximation ratio. Our fastest algorithms have an arbitrarily good polynomial approximation guarantee and an essentially optimal $\tilde{O}(Q)$ running time. On the other end of the spectrum, we achieve polylogarithmic approximations in $\tilde{O}(Q \cdot n^{\epsilon})$ rounds for any $\epsilon > 0$.


title: On the Scalar-Vector-Tensor Gravity: Black Hole, Thermodynamics and Geometrothermodynamics

abstract: Recently, a new class of modified gravity theories formulated via an additional scalar and vector field on top of the standard tensor field has been proposed. The direct implications of these theories are expected to be relevant for cosmology and astrophysics. In the present work, we revisit the modified framework of the scalar-vector-tensor theories of gravity. Surprisingly, we discover novel metric function for the black hole solutions. We also investigate the semi-classical thermodynamics of the black holes and study the thermodynamic properties of the obtained solutions. Moreover, we quantify the entropy and the temperature of the new black hole and also calculate the heat capacity. Finally, we also apply the formalism of the geometrothermodynamics to examine thermodynamic properties of the new black hole. This formalism yields results consistent with those obtained from the usual thermodynamic implementation.


title: Two-dimensional Iron Monocarbide with Planar Hypercoordinate Iron and Carbon

abstract: We report on the theoretical discovery of Iron monocarbide binary sheets stabilized at two-dimensional confined space, which we call tetragonal-FeC (t-FeC) and orthorhombic-FeC (o-FeC), respectively. From the energy viewpoint, the proposed t-FeC is the global minimum configuration in the 2D space, and each carbon atom is four-coordinated with ambient four Iron atoms. Strikingly, the o-FeC monolayer is an orthorhombic phase with planar pentacoordinate carbon moiety and planar seven-coordinate Fe moiety. To our knowledge, this monolayer is the first example of a simultaneously pentacoordinate carbon and planar seven-coordinate Fe-containing material. State-of-the-art theoretical calculations confirm that all these monolayers have significantly dynamic, mechanical, and thermal stabilities. Among these two monolayers, t-FeC monolayer shows a higher theoretical capacity (395 mAh g-1 ), and can stably adsorb Li up to t-FeCLi3 . Low migration energy barrier is predicted as small as 0.26 eV for Li, which result in the fast diffusion of Li atom on this monolayer. Moreover, electron-phonon calculations coupled with Bardeen-Cooper-Schrieffer arguments suggest t-FeC can be potential two-dimensional superconductors with 6.77 K superconducting transition temperature.


title: Compressive Spectrum Sensing for Cognitive Radio Networks

abstract: A cognitive radio system has the ability to observe and learn from the environment, adapt to the environmental conditions, and use the radio spectrum more efficiently. It allows secondary users (SUs) to use the primary users (PUs) channels when they are not being utilized. Cognitive radio involves three main processes: spectrum sensing, deciding, and acting. In the spectrum sensing process, the channel occupancy is measured with spectrum sensing techniques in order to detect unused channels. In the deciding process, sensing results are analyzed and decisions are made based on these results. In the acting process, actions are made by adjusting the transmission parameters to enhance the cognitive radio performance. One of the main challenges of cognitive radio is the wideband spectrum sensing. Existing spectrum sensing techniques are based on a set of observations sampled by an ADC at the Nyquist rate. However, those techniques can sense only one channel at a time because of the hardware limitations on the sampling rate. In addition, in order to sense a wideband spectrum, the wideband is divided into narrow bands or multiple frequency bands. SUs have to sense each band using multiple RF frontends simultaneously, which can result in a very high processing time, hardware cost, and computational complexity. In order to overcome this problem, the signal sampling should be as fast as possible even with high dimensional signals. Compressive sensing has been proposed as a low-cost solution to reduce the processing time and accelerate the scanning process. It allows reducing the number of samples required for high dimensional signal acquisition while keeping the essential information.


title: Understanding Convolutional Networks with APPLE : Automatic Patch Pattern Labeling for Explanation

abstract: With the success of deep learning, recent efforts have been focused on analyzing how learned networks make their classifications. We are interested in analyzing the network output based on the network structure and information flow through the network layers. We contribute an algorithm for 1) analyzing a deep network to find neurons that are 'important' in terms of the network classification outcome, and 2)automatically labeling the patches of the input image that activate these important neurons. We propose several measures of importance for neurons and demonstrate that our technique can be used to gain insight into, and explain how a network decomposes an image to make its final classification.


title: Differentiable Dynamic Programming for Structured Prediction and Attention

abstract: Dynamic programming (DP) solves a variety of structured combinatorial problems by iteratively breaking them down into smaller subproblems. In spite of their versatility, DP algorithms are usually non-differentiable, which hampers their use as a layer in neural networks trained by backpropagation. To address this issue, we propose to smooth the max operator in the dynamic programming recursion, using a strongly convex regularizer. This allows to relax both the optimal value and solution of the original combinatorial problem, and turns a broad class of DP algorithms into differentiable operators. Theoretically, we provide a new probabilistic perspective on backpropagating through these DP operators, and relate them to inference in graphical models. We derive two particular instantiations of our framework, a smoothed Viterbi algorithm for sequence prediction and a smoothed DTW algorithm for time-series alignment. We showcase these instantiations on two structured prediction tasks and on structured and sparse attention for neural machine translation.


title: The possibility of leptonic CP-violation measurement with JUNO

abstract: The existence of CP-violation in the leptonic sector is one of the most important issues for modern science. Neutrino physics is a key to the solution of this problem. JUNO (under construction) is the near future of neutrino physics. However CP-violation is not a priority for the current scientific program. We estimate the capability of $\delta_{\rm CP}$ measurement, assuming a combination of the JUNO detector and a superconductive cyclotron as the antineutrino source. This method of measuring CP-violation is an alternative to conventional beam experiments. A significance level of 3$\sigma$ can be reached for 22% of the $\delta_{\rm CP}$ range. The accuracy of measurement lies between 8$^{\rm o}$ and 22$^{\rm o}$. It is shown that the dominant influence on the result is the uncertainty in the mixing angle $\Theta_{23}$.


title: Boson-boson pure-dephasing model with non-Markovian properties

abstract: Understanding decoherence processes is crucial in the study of open quantum systems. In this paper, we discuss the mechanism of pure-dephasing process with a newly proposed boson-boson model, namely, a bosonic field coupled to another bosonic bath in thermal equilibrium. Our model is fully solvable and can reproduce the pure-dephasing process which is usually described by the well-known spin-boson model, therefore offering a new perspective to understanding decoherence processes in open quantum systems of high dimension. We also show that this model admits a generically non-Markovian dynamics with respect to various different non-Markovian measures.


title: Searching for the light Higgsinos in MSSM at the future e-p colliders

abstract: The search of the light Higgsino {in the Minimal Supersymmetric Standard Model (MSSM)} is a crucial test for the criteria of the naturalness in supersymmetry. On the other hand, the direct production of light Higgsino is also known as one of the most challenging SUSY searches at the current CERN Large Hadron collider (LHC). The lack of visible leptons due to the compressed spectrum and their small production rates limits their discovery potential in both mono-jet plus MET as well as the weak boson fusion (WBF) production. The signal $S/B$ ratio prediction is usually within the background systematic uncertainties at the LHC. Without color exchange between the beams, the $e-p$ colliders are well-known to have the WBF feature. Therefore, we study the search of the light Higgsinos at two future $e-p$ colliders at CERN, LHeC and FCC-eh.~The light Higgsinos will be produced in pair through weak boson fusion with controlled background at these colliders. We find the Higgsino of 95/145~GeV can be reached at 2$\sigma$ level at the future LHeC/FCC-eh with a luminosity $3 ~\text{ab}^{-1}$ respectively.


title: RoadTracer: Automatic Extraction of Road Networks from Aerial Images

abstract: Mapping road networks is currently both expensive and labor-intensive. High-resolution aerial imagery provides a promising avenue to automatically infer a road network. Prior work uses convolutional neural networks (CNNs) to detect which pixels belong to a road (segmentation), and then uses complex post-processing heuristics to infer graph connectivity. We show that these segmentation methods have high error rates because noisy CNN outputs are difficult to correct. We propose RoadTracer, a new method to automatically construct accurate road network maps from aerial images. RoadTracer uses an iterative search process guided by a CNN-based decision function to derive the road network graph directly from the output of the CNN. We compare our approach with a segmentation method on fifteen cities, and find that at a 5% error rate, RoadTracer correctly captures 45% more junctions across these cities.


title: On the boundary of the zero set of super-Brownian motion and its local time

abstract: If $X(t,x)$ is the density of one-dimensional super-Brownian motion, we prove that $\text{dim}(\partial\{x:X(t,x)>0\})=2-2\lambda_0\in(0,1)$ a.s. on $\{X_t\neq 0\}$, where $-\lambda_0\in(-1,-1/2)$ is the lead eigenvalue of a killed Ornstein-Uhlenbeck process. This confirms a conjecture of Mueller, Mytnik and Perkins who proved the above with positive probability. To establish this result we derive some new basic properties of a recently introduced boundary local time and analyze the behaviour of $X(t,\cdot)$ near the upper edge of its support. Numerical estimates of $\lambda_0$ suggest that the above Hausdorff dimension is approximately $.224$.


title: Non-Equilibrium Strongly Hyperuniform Fluids of Circle Active Particles with Large Local Density Fluctuations

abstract: Disordered hyperuniform structures are an exotic state of matter having vanishing long-wavelength density fluctuations similar to perfect crystals but without long-range order. Although its importance in materials science has been brought to the fore in past decades, the rational design of experimentally realizable disordered strongly hyperuniform microstructures remains challenging. Here we find a new type of non-equilibrium fluid with strong hyperuniformity in two-dimensional systems of chiral active particles, where particles perform independent circular motions of the radius R with the same handedness. This new hyperuniform fluid features a special length scale, i.e., the diameter of the circular trajectory of particles, below which large density fluctuations are observed. By developing a dynamic mean-field theory, we show that the large local density fluctuations can be explained as a motility-induced microphase separation, while the Fickian diffusion at large length scales and local center-of-mass-conserved noises are responsible for the global hyperuniformity.


title: Abelian quotients of the categories of short exact sequences

abstract: We mainly investigate abelian quotients of the categories of short exact sequences. The natural framework to consider the question is via identifying quotients of morphism categories as modules categories. These ideas not only can be used to recover the abelian quotients produced by cluster-tilting subcategories of both exact categories and triangulated categories, but also can be used to reach our goal. Let $(\mathcal{C},\mathcal{E})$ be an exact category. We denote by $\mathcal{E}(\mathcal{C})$ the category of bounded complexes whose objects are given by short exact sequences in $\mathcal{E}$ and by $S\mathcal{E}(\mathcal{C})$ the full subcategory formed by split short exact sequences. In general, $\mathcal{E}(\mathcal{C})$ is just an exact category, but the quotient $\mathcal{E}(\mathcal{C})/[S\mathcal{E}(\mathcal{C})]$ turns out to be abelian. In particular, if $(\mathcal{C},\mathcal{E})$ is Frobenius, we present three equivalent abelian quotients of $\mathcal{E}(\mathcal{C})$ and point out that the equivalences are actually given by left and right rotations. The abelian quotient $\mathcal{E}(\mathcal{C})/[S\mathcal{E}(\mathcal{C})]$ admits some nice properties. We explicitly describe the abelian structure, projective objects, injective objects and simple objects, which provide a new viewpoint to understanding Hilton-Rees Theorem and Auslander-Reiten theory. Furthermore, we present some analogous results both for $n$-exact versions and for triangulated versions.


title: Ball Prolate Spheroidal Wave Functions In Arbitrary Dimensions

abstract: In this paper, we introduce the prolate spheroidal wave functions (PSWFs) of real order $\alpha>-1$ on the unit ball in arbitrary dimension, termed as ball PSWFs. They are eigenfunctions of both a weighted concentration integral operator, and a Sturm-Liouville differential operator. Different from existing works on multi-dimensional PSWFs, the ball PSWFs are defined as a generalisation of orthogonal {\em ball polynomials} in primitive variables with a tuning parameter $c>0$, through a "perturbation" of the Sturm-Liouville equation of the ball polynomials. From this perspective, we can explore some interesting intrinsic connections between the ball PSWFs and the finite Fourier and Hankel transforms. We provide an efficient and accurate algorithm for computing the ball PSWFs and the associated eigenvalues, and present various numerical results to illustrate the efficiency of the method. Under this uniform framework, we can recover the existing PSWFs by suitable variable substitutions.


title: Learning a SAT Solver from Single-Bit Supervision

abstract: We present NeuroSAT, a message passing neural network that learns to solve SAT problems after only being trained as a classifier to predict satisfiability. Although it is not competitive with state-of-the-art SAT solvers, NeuroSAT can solve problems that are substantially larger and more difficult than it ever saw during training by simply running for more iterations. Moreover, NeuroSAT generalizes to novel distributions; after training only on random SAT problems, at test time it can solve SAT problems encoding graph coloring, clique detection, dominating set, and vertex cover problems, all on a range of distributions over small random graphs.


title: Physical renormalization condition for de Sitter QED

abstract: We considered a new renormalization condition for the vacuum expectation values of the scalar and spinor currents induced by a homogeneous and constant electric field background in de Sitter spacetime. Following a semiclassical argument, the condition named maximal subtraction imposes the exponential suppression on the massive charged particle limit of the renormalized currents. The maximal subtraction changes the behaviors of the induced currents previously obtained by the conventional minimal subtraction scheme. The maximal subtraction is favored for a couple of physically decent predictions including the identical asymptotic behavior of the scalar and spinor currents, the removal of the infrared (IR) hyperconductivity from the scalar current, and the finite current for the massless fermion.


title: Computation of Transmission Eigenvalues for Elastic Waves

abstract: The goal of this paper is to develop numerical methods computing a few smallest elasticity transmission eigenvalues, which are of practical importance in inverse scattering theory. The problem is challenging since it is nonlinear, non-self-adjoint, and of fourth order. We construct a nonlinear function whose values are generalized eigenvalues of a series of self-adjoint fourth order problems. The roots of the function are the transmission eigenvalues. Using an $H^2$-conforming finite element for the self-adjoint fourth order eigenvalue problems, we employ a secant method to compute the roots of the nonlinear function. The convergence of the proposed method is proved. In addition, a mixed finite element method is developed for the purpose of verification. Numerical examples are presented to verify the theory and demonstrate the effectiveness of the two methods.


title: On the Rates of Convergence from Surrogate Risk Minimizers to the Bayes Optimal Classifier

abstract: We study the rates of convergence from empirical surrogate risk minimizers to the Bayes optimal classifier. Specifically, we introduce the notion of \emph{consistency intensity} to characterize a surrogate loss function and exploit this notion to obtain the rate of convergence from an empirical surrogate risk minimizer to the Bayes optimal classifier, enabling fair comparisons of the excess risks of different surrogate risk minimizers. The main result of the paper has practical implications including (1) showing that hinge loss is superior to logistic and exponential loss in the sense that its empirical minimizer converges faster to the Bayes optimal classifier and (2) guiding to modify surrogate loss functions to accelerate the convergence to the Bayes optimal classifier.


title: Dual Control Memory Augmented Neural Networks for Treatment Recommendations

abstract: Machine-assisted treatment recommendations hold a promise to reduce physician time and decision errors. We formulate the task as a sequence-to-sequence prediction model that takes the entire time-ordered medical history as input, and predicts a sequence of future clinical procedures and medications. It is built on the premise that an effective treatment plan may have long-term dependencies from previous medical history. We approach the problem by using a memory-augmented neural network, in particular, by leveraging the recent differentiable neural computer that consists of a neural controller and an external memory module. But differing from the original model, we use dual controllers, one for encoding the history followed by another for decoding the treatment sequences. In the encoding phase, the memory is updated as new input is read; at the end of this phase, the memory holds not only the medical history but also the information about the current illness. During the decoding phase, the memory is write-protected. The decoding controller generates a treatment sequence, one treatment option at a time. The resulting dual controller write-protected memory-augmented neural network is demonstrated on the MIMIC-III dataset on two tasks: procedure prediction and medication prescription. The results show improved performance over both traditional bag-of-words and sequence-to-sequence methods.


title: On the Generalization of Equivariance and Convolution in Neural Networks to the Action of Compact Groups

abstract: Convolutional neural networks have been extremely successful in the image recognition domain because they ensure equivariance to translations. There have been many recent attempts to generalize this framework to other domains, including graphs and data lying on manifolds. In this paper we give a rigorous, theoretical treatment of convolution and equivariance in neural networks with respect to not just translations, but the action of any compact group. Our main result is to prove that (given some natural constraints) convolutional structure is not just a sufficient, but also a necessary condition for equivariance to the action of a compact group. Our exposition makes use of concepts from representation theory and noncommutative harmonic analysis and derives new generalized convolution formulae.


title: Tree-to-tree Neural Networks for Program Translation

abstract: Program translation is an important tool to migrate legacy code in one language into an ecosystem built in a different language. In this work, we are the first to employ deep neural networks toward tackling this problem. We observe that program translation is a modular procedure, in which a sub-tree of the source tree is translated into the corresponding target sub-tree at each step. To capture this intuition, we design a tree-to-tree neural network to translate a source tree into a target one. Meanwhile, we develop an attention mechanism for the tree-to-tree model, so that when the decoder expands one non-terminal in the target tree, the attention mechanism locates the corresponding sub-tree in the source tree to guide the expansion of the decoder. We evaluate the program translation capability of our tree-to-tree model against several state-of-the-art approaches. Compared against other neural translation models, we observe that our approach is consistently better than the baselines with a margin of up to 15 points. Further, our approach can improve the previous state-of-the-art program translation approaches by a margin of 20 points on the translation of real-world projects.


title: Nearly Optimal Adaptive Procedure with Change Detection for Piecewise-Stationary Bandit

abstract: Multi-armed bandit (MAB) is a class of online learning problems where a learning agent aims to maximize its expected cumulative reward while repeatedly selecting to pull arms with unknown reward distributions. We consider a scenario where the reward distributions may change in a piecewise-stationary fashion at unknown time steps. We show that by incorporating a simple change-detection component with classic UCB algorithms to detect and adapt to changes, our so-called M-UCB algorithm can achieve nearly optimal regret bound on the order of $O(\sqrt{MKT\log T})$, where $T$ is the number of time steps, $K$ is the number of arms, and $M$ is the number of stationary segments. Comparison with the best available lower bound shows that our M-UCB is nearly optimal in $T$ up to a logarithmic factor. We also compare M-UCB with the state-of-the-art algorithms in numerical experiments using a public Yahoo! dataset to demonstrate its superior performance.


title: Stable Cosmic Time Crystals

abstract: Cosmological time crystals are created when a scalar field moves periodically through phase space in a spatially flat Friedmann-Robertson-Walker spacetime due to the presence of a limit cycle. All such cosmological time crystals in the literature suffer from gradient instabilities occurring at Null Energy Condition violating phases where the square sound speed for cosmological perturbations becomes negative. Here we present stable cosmological time crystals. Our analysis suggests this new form of scalar matter--cosmic time crystals--may be considered as a physically viable cosmological matter source.


title: The Spectrum of the Universe

abstract: The cosmic background (CB) radiation, encompassing the sum of emission from all sources outside our own Milky Way galaxy across the entire electromagnetic spectrum, is a fundamental phenomenon in observational cosmology. Many experiments have been conceived to measure it (or its constituents) since the extragalactic Universe was first discovered; in addition to estimating the bulk (cosmic monopole) spectrum, directional variations have also been detected over a wide range of wavelengths. Here we gather the most recent of these measurements and discuss the current status of our understanding of the CB from radio to $\gamma$-ray energies. Using available data in the literature we piece together the sky-averaged intensity spectrum, and discuss the emission processes responsible for what is observed. We examine the effect of perturbations to the continuum spectrum from atomic and molecular line processes and comment on the detectability of these signals. We also discuss how one could in principle obtain a complete census of the CB by measuring the full spectrum of each spherical harmonic expansion coefficient. This set of spectra of multipole moments effectively encodes the entire statistical history of nuclear, atomic and molecular processes in the Universe.


title: The lamplighter group of rank two generated by a bireversible automaton

abstract: We construct a 4-state 2-letter bireversible automaton generating the lamplighter group $(\mathbb Z_2^2)\wr\mathbb Z$ of rank two. The action of the generators on the boundary of the tree can be induced by the affine transformations on the ring $\mathbb Z_2[[t]]$ of formal power series over $\mathbb Z_2$.


title: A novel numerical scheme for nonlinear electron-plasma oscillations

abstract: In this work, we suggest an easy-to-code higher-order finite volume semi-discrete scheme to analyze the nonlinear behavior of the electron-plasma oscillations by solving electron fluid equations numerically. The present method employs a fourth-order accurate centrally weighted essentially nonoscillatory reconstruction (CWENO4) polynomial for estimating the numerical flux at the grid-cell interfaces, and a fourth-order Runge-Kutta method for the time integration. The numerical implementation is validated by reproducing earlier results for both non-dissipative and dissipative cold plasmas. The stability of the present scheme is illustrated by evolving the nonlinear electron plasma oscillations in a cold non-dissipative plasma for hundred plasma periods, which also display a negligible numerical dissipation. The fourth-order accuracy of the existing approach is also confirmed by evaluating the convergence of errors for the nonlinear electron plasma oscillations in a cold non-dissipative plasma.


title: Holographic Magnetic Susceptibility

abstract: The (2+1)-dimensional static magnetic susceptibility in strong-coupling is studied via a Reissner-Nordstr\"{o}m-AdS geometry. The analyticity of the susceptibility on the complex momentum $\mathfrak{q}$-plane in relation to the Friedel-like oscillation in coordinate space is explored. In contrast to the branch-cuts crossing the real momentum-axis for a Fermi liquid, we prove that the holographic magnetic susceptibility remains an analytic function of the complex momentum around the real axis in the limit of zero temperature, At zero temperature, we located analytically two pairs of branch-cuts that are parallel to the imaginary momentum-axis for large $|\text{Im}\ \mathfrak{q}|$ but become warped with the end-points keeping away from the real and imaginary momentum-axes. We conclude that these branch-cuts give rise to the exponential decay behaviour of Friedel-like oscillation of magnetic susceptibility in coordinate space. We also derived the analytical forms of the susceptibility in large and small-momentum, respectively.


title: A Smooth Curve as a Fractal Under the Third Definition

abstract: It is commonly believed in the literature that smooth curves, such as circles, are not fractal, and only non-smooth curves, such as coastlines, are fractal. However, this paper demonstrates that a smooth curve can be fractal, under the new, relaxed, third definition of fractal - a set or pattern is fractal if the scaling of far more small things than large ones recurs at least twice. The scaling can be rephrased as a hierarchy, consisting of numerous smallest, a very few largest, and some in between the smallest and the largest. The logarithmic spiral, as a smooth curve, is apparently fractal because it bears the self-similar property, or the scaling of far more small squares than large ones recurs multiple times, or the scaling of far more small bends than large ones recurs multiple times. A half-circle or half-ellipse and the UK coastline (before or after smooth processing) are fractal, if the scaling of far more small bends than large ones recurs at least twice. Keywords: Third definition of fractal, head/tail breaks, bends, ht-index, scaling hierarchy


title: PCA-Based Missing Information Imputation for Real-Time Crash Likelihood Prediction Under Imbalanced Data

abstract: The real-time crash likelihood prediction has been an important research topic. Various classifiers, such as support vector machine (SVM) and tree-based boosting algorithms, have been proposed in traffic safety studies. However, few research focuses on the missing data imputation in real-time crash likelihood prediction, although missing values are commonly observed due to breakdown of sensors or external interference. Besides, classifying imbalanced data is also a difficult problem in real-time crash likelihood prediction, since it is hard to distinguish crash-prone cases from non-crash cases which compose the majority of the observed samples. In this paper, principal component analysis (PCA) based approaches, including LS-PCA, PPCA, and VBPCA, are employed for imputing missing values, while two kinds of solutions are developed to solve the problem in imbalanced data. The results show that PPCA and VBPCA not only outperform LS-PCA and other imputation methods (including mean imputation and k-means clustering imputation), in terms of the root mean square error (RMSE), but also help the classifiers achieve better predictive performance. The two solutions, i.e., cost-sensitive learning and synthetic minority oversampling technique (SMOTE), help improve the sensitivity by adjusting the classifiers to pay more attention to the minority class.


title: Stochastic Non-preemptive Co-flow Scheduling with Time-Indexed Relaxation

abstract: Co-flows model a modern scheduling setting that is commonly found in a variety of applications in distributed and cloud computing. A stochastic co-flow task contains a set of parallel flows with randomly distributed sizes. Further, many applications require non-preemptive scheduling of co-flow tasks. This paper gives an approximation algorithm for stochastic non-preemptive co-flow scheduling. The proposed approach uses a time-indexed linear relaxation, and uses its solution to come up with a feasible schedule. This algorithm is shown to achieve a competitive ratio of $(2\log{m}+1)(1+\sqrt{m}\Delta)(1+m{\Delta}){(3+\Delta)}/{2}$ for zero-release times, and $(2\log{m}+1)(1+\sqrt{m}\Delta)(1+m\Delta)(2+\Delta)$ for general release times, where $\Delta$ represents the upper bound of squared coefficient of variation of processing times, and $m$ is the number of servers.


title: Formal Ontology Learning from English IS-A Sentences

abstract: Ontology learning (OL) is the process of automatically generating an ontological knowledge base from a plain text document. In this paper, we propose a new ontology learning approach and tool, called DLOL, which generates a knowledge base in the description logic (DL) SHOQ(D) from a collection of factual non-negative IS-A sentences in English. We provide extensive experimental results on the accuracy of DLOL, giving experimental comparisons to three state-of-the-art existing OL tools, namely Text2Onto, FRED, and LExO. Here, we use the standard OL accuracy measure, called lexical accuracy, and a novel OL accuracy measure, called instance-based inference model. In our experimental results, DLOL turns out to be about 21% and 46%, respectively, better than the best of the other three approaches.


title: Disentangling Entanglements in Biopolymer Solutions

abstract: Reptation theory has been highly successful in explaining the unusual material properties of entangled polymer solutions. It reduces the complex many-body dynamics to a single-polymer description where each polymer is envisaged to be confined to a tube through which it moves in a snake-like fashion. For flexible polymers, reptation theory has been amply confirmed by both experiments and simulations. In contrast, for semiflexible polymers experimental and numerical tests are either limited to the onset of reptation, or were performed for tracer polymers in a fixed, static matrix. Here we report Brownian dynamics simulations of entangled solutions of semiflexible polymers, which show that curvilinear motion along a tube (reptation) is no longer the dominant mode of dynamics. Instead, we find that polymers disentangle due to correlated constraint release which leads to equilibration of internal bending modes before polymers diffuse the full tube length. The physical mechanism underlying terminal stress relaxation is rotational diffusion mediated by disentanglement rather than curvilinear motion along a tube.


title: Stochastic Spectral and Conjugate Descent Methods

abstract: The state-of-the-art methods for solving optimization problems in big dimensions are variants of randomized coordinate descent (RCD). In this paper we introduce a fundamentally new type of acceleration strategy for RCD based on the augmentation of the set of coordinate directions by a few spectral or conjugate directions. As we increase the number of extra directions to be sampled from, the rate of the method improves, and interpolates between the linear rate of RCD and a linear rate independent of the condition number. We develop and analyze also inexact variants of these methods where the spectral and conjugate directions are allowed to be approximate only. We motivate the above development by proving several negative results which highlight the limitations of RCD with importance sampling.


title: TARM: A Turbo-type Algorithm for Affine Rank Minimization

abstract: The affine rank minimization (ARM) problem arises in many real-world applications. The goal is to recover a low-rank matrix from a small amount of noisy affine measurements. The original problem is NP-hard, and so directly solving the problem is computationally prohibitive. Approximate low-complexity solutions for ARM have recently attracted much research interest. In this paper, we design an iterative algorithm for ARM based on message passing principles. The proposed algorithm is termed turbo-type ARM (TARM), as inspired by the recently developed turbo compressed sensing algorithm for sparse signal recovery. We show that, when the linear operator for measurement is right-orthogonally invariant (ROIL), a scalar function called state evolution can be established to accurately predict the behaviour of the TARM algorithm. We also show that TARM converges much faster than the counterpart algorithms for low-rank matrix recovery. We further extend the TARM algorithm for matrix completion, where the measurement operator corresponds to a random selection matrix. We show that, although the state evolution is not accurate for matrix completion, the TARM algorithm with carefully tuned parameters still significantly outperforms its counterparts.


title: The Gaussian wave packets transform for the semi-classical Schr\"odinger equation with vector potentials

abstract: In this paper, we reformulate the semi-classical Schr\"odinger equation in the presence of electromagnetic field by the Gaussian wave packets transform. With this approach, the highly oscillatory Schr\"odinger equation is equivalently transformed into another Schr\"odinger type wave equation, the $w$ equation, which is essentially not oscillatory and thus requires much less computational effort. We propose two numerical methods to solve the $w$ equation, where the Hamiltonian is either divided into the kinetic, the potential and the convection part, or into the kinetic and the potential-convection part. The convection, or the potential-convection part is treated by a semi-Lagrangian method, while the kinetic part is solved by the Fourier spectral method. The numerical methods are proved to be unconditionally stable, spectrally accurate in space and second order accurate in time, and in principle they can be extended to higher order schemes in time. Various one dimensional and multidimensional numerical tests are provided to justify the properties of the proposed methods.


title: FDM-Structured Preamble Optimization for Channel Estimation in MIMO-OQAM/FBMC Systems

abstract: In this paper, we consider the problem of preamble design in multiple-input multiple-output (MIMO) systems employing offset quadrature amplitude modulation based filter bank multicarrier (OQAM/FBMC) and propose a preamble optimization method for the frequency division multiplexing (FDM)-structured preamble. Specifically, we formulate an optimization problem to determine the frequency division multiplexed preambles, where the objective is to minimize the mean square error (MSE) of the channel estimation, subject to the constraint on the transmit energy. For two transmit antennas, we find the relationship between preambles and the intrinsic interference from neighboring symbols to achieve the minimum channel estimation MSE, and derive the optimal closed-form solution. For more than two transmit antennas, the constrained preamble optimization problem is nonconvex quadratic. Therefore, we convert the original optimization problem into a quadratically constrained quadratic program (QCQP) and obtain the suboptimal solution by relaxing the nonconvex constraint. Simulation results demonstrate that, in terms of MSE and bit error rate (BER) performances, the proposed method outperforms the conventional FDM preamble design method at all signal-to-noise ratio (SNR) regimes and outperforms the interference approximation method-complex (IAM-C) preamble design method at low to medium SNR regimes with lower preamble overhead.


title: The Need for Speed of AI Applications: Performance Comparison of Native vs. Browser-based Algorithm Implementations

abstract: AI applications pose increasing demands on performance, so it is not surprising that the era of client-side distributed software is becoming important. On top of many AI applications already using mobile hardware, and even browsers for computationally demanding AI applications, we are already witnessing the emergence of client-side (federated) machine learning algorithms, driven by the interests of large corporations and startups alike. Apart from mathematical and algorithmic concerns, this trend especially demands new levels of computational efficiency from client environments. Consequently, this paper deals with the question of state-of-the-art performance by presenting a comparison study between native code and different browser-based implementations: JavaScript, ASM.js as well as WebAssembly on a representative mix of algorithms. Our results show that current efforts in runtime optimization push the boundaries well towards (and even beyond) native binary performance. We analyze the results obtained and speculate on the reasons behind some surprises, rounding the paper off by outlining future possibilities as well as some of our own research efforts.


title: A Time-Varying Network for Cryptocurrencies

abstract: Cryptocurrencies return cross-predictability and technological similarity yield information on risk propagation and market segmentation. To investigate these effects, we build a time-varying network for cryptocurrencies, based on the evolution of return cross-predictability and technological similarities. We develop a dynamic covariate-assisted spectral clustering method to consistently estimate the latent community structure of cryptocurrencies network that accounts for both sets of information. We demonstrate that investors can achieve better risk diversification by investing in cryptocurrencies from different communities. A cross-sectional portfolio that implements an inter-crypto momentum trading strategy earns a 1.08% daily return. By dissecting the portfolio returns on behavioral factors, we confirm that our results are not driven by behavioral mechanisms.

categories: stat.ME econ.EM q-fin.PM q-fin.RM stat.AP


title: Considerations for a Multi-beam Multi-purpose Survey with FAST

abstract: Having achieved 'first-light' right before the opening ceremony on September 25, 2016, the Five-hundred-meter Aperture Spherical radio Telescope (FAST) is being busily commissioned. Its innovative design requires ~1000 points to be measured and driven instead of just the two axes of motion, e.g. Azimuth and Elevation for most of the conventional antennae, to realize pointing and tracking. We have devised a survey plan to utilized the full sensitivity of FAST, while minimizing the complexities in operation the system. The 19-beam L band focal plan array will be rotated to specific angles and taking continuous data streams while the surface shape and the focal cabin stay fixed. Such a survey will cover the northern sky in about 220 full days. Our aim is to obtain data for pulsar search, HI (neutral hydrogen) galaxies, HI imaging, and radio transients, simultaneously, through multiple backends. These data sets could be a significant contribution to all related fields in radio astronomy and remain relevant for decades.


title: Simultaneous 13 cm/3 cm single-pulse observations of PSR B0329+54

abstract: We have investigated the mode-changing properties of PSR B0329+54 using 31 epochs of simultaneous 13 cm/3 cm single-pulse observations obtained with Shanghai Tian Ma 65 m telescope. The pulsar was found in the abnormal emission mode 17 times, accounting for ~13% of the 41.6 hours total observation time. Single pulse analyses indicate that mode changes took place simultaneously at 13 cm/3 cm within a few rotational periods. We detected occasional bright and narrow pulses whose peak flux densities were 10 times higher than that of the integrated profile in both bands. At 3 cm, about 0.66% and 0.27% of single pulses were bright in the normal mode and abnormal mode respectively, but at 13 cm the occurrence rate was only about 0.007%. We divided the pulsar radiation window into three components (C1, C2 and C3) corresponding to the main peaks of the integrated profile. The bright pulses preferentially occurred at pulse phases corresponding to the peaks of C2 and C3. Fluctuation spectra showed that C2 had excess red noise in the normal mode, but broad quasi-periodic features with central frequencies around 0.12 cycles/period in the abnormal mode. At 3 cm, C3 had a stronger quasi-periodic modulation centered around 0.06 cycles/period in the abnormal mode. Although there were some asymmetries in the two-dimensional fluctuation spectra, we found no clear evidence for systematic subpulse drifting. Consistent with previous low-frequency observations, we found a very low nulling probability for B0329+54 with upper limits of 0.13% and 1.68% at 13 cm/3 cm respectively.


title: Kazhdan-Lusztig polynomials of fan matroids, wheel matroids and whirl matroids

abstract: The Kazhdan-Lusztig polynomial of a matroid was introduced by Elias, Proudfoot and Wakefield, whose properties need to be further explored. In this paper we prove that the Kazhdan-Lusztig polynomials of fan matroids coincide with Motzkin polynomials, which was recently conjectured by Gedeon. As a byproduct, we determine the Kazhdan-Lusztig polynomials of graphic matroids of squares of paths. We further obtain explicit formulas of the Kazhdan-Lusztig polynomials of wheel matroids and whirl matroids. We prove the real-rootedness of the Kazhdan-Lusztig polynomials of these matroids, which provides positive evidence for a conjecture due to Gedeon, Proudfoot and Young. Based on the results on the Kazhdan-Lusztig polynomials, we also determine the $Z$-polynomials of fan matroids, wheel matroids and whirl matroids, and prove their real-rootedness, which provides further evidence in support of a conjecture of Proudfoot, Xu, and Young.


title: Syntax and Semantics of Italian Poetry in the First Half of the 20th Century

abstract: In this paper we study, analyse and comment rhetorical figures present in some of most interesting poetry of the first half of the twentieth century. These figures are at first traced back to some famous poet of the past and then compared to classical Latin prose. Linguistic theory is then called in to show how they can be represented in syntactic structures and classified as noncanonical structures, by positioning discontinuous or displaced linguistic elements in Spec XP projections at various levels of constituency. Then we introduce LFG (Lexical Functional Grammar) as the theory that allows us to connect syntactic noncanonical structures with informational structure and psycholinguistic theories for complexity evaluation. We end up with two computational linguistics experiments and then evaluate the results. The first one uses best online parsers of Italian to parse poetic structures; the second one uses Getarun, the system created at Ca Foscari Computational Linguistics Laboratory. As will be shown, the first approach is unable to cope with these structures due to the use of only statistical probabilistic information. On the contrary, the second one, being a symbolic rule based system, is by far superior and allows also to complete both semantic an pragmatic analysis.


title: $\mathcal{G}$-SGD: Optimizing ReLU Neural Networks in its Positively Scale-Invariant Space

abstract: It is well known that neural networks with rectified linear units (ReLU) activation functions are positively scale-invariant. Conventional algorithms like stochastic gradient descent optimize the neural networks in the vector space of weights, which is, however, not positively scale-invariant. This mismatch may lead to problems during the optimization process. Then, a natural question is: \emph{can we construct a new vector space that is positively scale-invariant and sufficient to represent ReLU neural networks so as to better facilitate the optimization process }? In this paper, we provide our positive answer to this question. First, we conduct a formal study on the positive scaling operators which forms a transformation group, denoted as $\mathcal{G}$. We show that the value of a path (i.e. the product of the weights along the path) in the neural network is invariant to positive scaling and prove that the value vector of all the paths is sufficient to represent the neural networks under mild conditions. Second, we show that one can identify some basis paths out of all the paths and prove that the linear span of their value vectors (denoted as $\mathcal{G}$-space) is an invariant space with lower dimension under the positive scaling group. Finally, we design stochastic gradient descent algorithm in $\mathcal{G}$-space (abbreviated as $\mathcal{G}$-SGD) to optimize the value vector of the basis paths of neural networks with little extra cost by leveraging back-propagation. Our experiments show that $\mathcal{G}$-SGD significantly outperforms the conventional SGD algorithm in optimizing ReLU networks on benchmark datasets.


title: Lightweight Classification of IoT Malware based on Image Recognition

abstract: The Internet of Things (IoT) is an extension of the traditional Internet, which allows a very large number of smart devices, such as home appliances, network cameras, sensors and controllers to connect to one another to share information and improve user experiences. Current IoT devices are typically micro-computers for domain-specific computations rather than traditional functionspecific embedded devices. Therefore, many existing attacks, targeted at traditional computers connected to the Internet, may also be directed at IoT devices. For example, DDoS attacks have become very common in IoT environments, as these environments currently lack basic security monitoring and protection mechanisms, as shown by the recent Mirai and Brickerbot IoT botnets. In this paper, we propose a novel light-weight approach for detecting DDos malware in IoT environments.We firstly extract one-channel gray-scale images converted from binaries, and then utilize a lightweight convolutional neural network for classifying IoT malware families. The experimental results show that the proposed system can achieve 94.0% accuracy for the classification of goodware and DDoS malware, and 81.8% accuracy for the classification of goodware and two main malware families.


title: Defect-Free Carbon Nanotube Coils

abstract: Carbon nanotubes are promising building blocks for various nanoelectronic components. A highly desirable geometry for such applications is a coil. However, coiled nanotube structures reported so far were inherently defective or had no free ends accessible for contacting. Here we demonstrate the spontaneous self-coiling of single-wall carbon nanotubes into defect-free coils of up to more than 70 turns with identical diameter and chirality, and free ends. We characterize the structure, formation mechanism, and electrical properties of these coils by different microscopies, molecular dynamics simulations, Raman spectroscopy, and electrical and magnetic measurements. The coils are highly conductive, as expected for defect-free carbon nanotubes, but adjacent nanotube segments in the coil are more highly coupled than in regular bundles of single-wall carbon nanotubes, owing to their perfect crystal momentum matching, which enables tunneling between the turns. Although this behavior does not yet enable the performance of these nanotube coils as inductive devices, it does point a clear path for their realization. Hence, this study represents a major step toward the production of many different nanotube coil devices, including inductors, electromagnets, transformers, and dynamos.


title: Bimodal logics with contingency and accident

abstract: Contingency and accident are two important notions in philosophy and philosophical logic. Their meanings are so close that they are mixed sometimes, in both everyday discourse and academic research. This indicates that it is necessary to study them in a unified framework. However, there has been no logical research on them together. In this paper, we propose a language of a bimodal logic with these two concepts, investigate its model-theoretical properties such as expressivity and frame definability. We axiomatize this logic over various classes of frames, whose completeness proofs are shown with the help of a crucial schema. The interactions between contingency and accident can sharpen our understanding of both notions. Then we extend the logic to a dynamic case: public announcements. By finding the required reduction axioms, we obtain a complete axiomatization, which gives us a good application to Moore sentences.


title: Localized peaking regimes for quasilinear parabolic equations

abstract: This paper deals with the asymptotic behavior as $t\rightarrow T<\infty$ of all weak (energy) solutions of a class of equations with the following model representative: \begin{equation*} (|u|^{p-1}u)_t-\Delta_p(u)+b(t,x)|u|^{\lambda-1}u=0 \quad (t,x)\in(0,T)\times\Omega,\,\Omega\in{R}^n,\,n>1, \end{equation*} with prescribed global energy function \begin{equation*} E(t):=\int_{\Omega}|u(t,x)|^{p+1}dx+ \int_0^t\int_{\Omega}|\nabla_xu(\tau,x)|^{p+1}dxd\tau \rightarrow\infty\ \text{ as }t\rightarrow T. \end{equation*} Here $\Delta_p(u)=\sum_{i=1}^n\left(|\nabla_xu|^{p-1}u_{x_i}\right)_{x_i}$, $p>0$, $\lambda>p$, $\Omega$ is a bounded smooth domain, $b(t,x)\geq0$. Particularly, in the case \begin{equation*} E(t)\leq F_\mu(t)=\exp\left(\omega(T-t)^{-\frac1{p+\mu}}\right)\quad\forall\,t<T,\,\mu>0,\,\omega>0, \end{equation*} it is proved that solution $u$ remains uniformly bounded as $t\rightarrow T$ in an arbitrary subdomain $\Omega_0\subset\Omega:\overline{\Omega}_0\subset\Omega$ and the sharp upper estimate of $u(t,x)$ when $t\rightarrow T$ has been obtained depending on $\mu>0$ and $s=dist(x,\partial\Omega)$. In the case $b(t,x)>0$ $\forall\,(t,x)\in(0,T)\times\Omega$ sharp sufficient conditions on degeneration of $b(t,x)$ near $t=T$ that guarantee mentioned above boundedness for arbitrary (even large) solution have been found and the sharp upper estimate of a final profile of solution when $t\rightarrow T$ has been obtained.


title: Parametrization of quantum states and the quantum state discrimination problem

abstract: A discrimination problem consists of $N$ linearly independent pure quantum states $\Phi=\{\ket{\phi_i}\}$ and the corresponding occurrence probabilities $\eta=\{\eta_i\}$. To any such problem we associate, up to a permutation over the probabilities $\{\eta_i\}$, a unique pair of density matrices $\boldsymbol{\rho_{_{T}}}$ and $\boldsymbol{\eta_{{p}}}$ defined on the $N$-dimensional Hilbert space $\mathcal{H}_N$. The first one, $\boldsymbol{\rho_{_{T}}}$, provides a new parametrization of a generic full-rank density matrix in terms of the parameters of the discrimination problem, i.e. the mutual overlaps $\gamma_{ij}=\bra{\phi_i}\phi_j\rangle$ and the occurrence probabilities $\{\eta_i\}$. The second one is defined as a diagonal density matrix $\boldsymbol{\eta_p}$ with the diagonal entries given by the probabilities $\{\eta_i\}$ with the ordering induced by the permutation $p$ of the probabilities. $\boldsymbol{\rho_{_{T}}}$ and $\boldsymbol{\eta_{{p}}}$ capture information about the quantum and classical versions of the discrimination problem, respectively. In this sense, when the set $\Phi$ can be discriminated unambiguously with probability one, i.e. when the states to be discriminated are mutually orthogonal and can be distinguished by a classical observer, then $\boldsymbol{\rho_{_{T}}}\rightarrow \boldsymbol{\eta_{{p}}}$. Moreover, if the set lacks its independency and cannot be discriminated anymore the distinguishability of the pair, measured by the fidelity $F(\boldsymbol{\rho_{_{T}}}, \boldsymbol{\eta_{{p}}})$, becomes minimum. This enables one to associate to each discrimination problem a measure of discriminability defined by the fidelity $F(\boldsymbol{\rho_{_{T}}}, \boldsymbol{\eta_{{p}}})$. This quantity, has the advantage of being easy to calculate and in this respect it can find useful applications in estimating the extent to which the set is discriminable.


title: Encoding and avoiding 2-connected patterns in polygon dissections and outerplanar graphs

abstract: Let $\Delta =\{ \delta_1,\delta_2,...,\delta_m \} $ be a finite set of 2-connected patterns, i.e. graphs up to vertex relabelling. We study the generating function $D_{\Delta }(z,u_1,u_2,...,u_m),$ which counts polygon dissections and marks subgraph copies of $\delta_i$ with the variable $u_i$. We prove that this is always algebraic, through an explicit combinatorial decomposition depending on $\Delta $. The decomposition also gives a defining system for $D_{\Delta }(z,\mathbf{0})$, which encodes polygon dissections that restrict these patterns as subgraphs. In this way, we are able to extract normal limit laws for the patterns when they are encoded, and perform asymptotic enumeration of the resulting classes when they are avoided. The results can be directly transferred in the case of labelled outerplanar graphs. We give examples and compute the relevant constants when the patterns are small cycles or dissections.


title: Double Minimum Variance Beamforming Method to Enhance Photoacoustic Imaging

abstract: One of the common algorithms used to reconstruct photoacoustic (PA) images is the non-adaptive Delay-and-Sum (DAS) beamformer. However, the quality of the reconstructed PA images obtained by DAS is not satisfying due to its high level of sidelobes and wide mainlobe. In contrast, adaptive beamformers, such as minimum variance (MV), result in an improved image compared to DAS. In this paper, a novel beamforming method, called Double MV (D-MV) is proposed to enhance the image quality compared to the MV. It is shown that there is a summation procedure between the weighted subarrays in the output of the MV beamformer. This summation can be interpreted as the non-adaptive DAS beamformer. It is proposed to replace the existing DAS with the MV algorithm to reduce the contribution of the off-axis signals caused by the DAS beamformer between the weighted subarrays. The numerical results show that the proposed technique improves the full-width-half-maximum (FWHM) and signal-to-noise ratio (SNR) for about 28.83 \mu m and 4.8 dB in average, respectively, compared to MV beamformer. Also, quantitative evaluation of the experimental results indicates that the proposed D-MV leads to 0.15 mm and 1.96 dB improvement in FWHM and SNR, in comparison with MV beamformer.


title: Lack of superconductivity in the phase diagram of single-crystalline Eu(Fe1-xCox)2As2 grown by transition metal arsenide flux

abstract: The interplay of magnetism and superconductivity (SC) has been a focus of interest in condensed matter physics for decades. EuFe2As2 has been identified as a potential platform to investigate interactions between structural, magnetic, electronic effects as well as coexistence of magnetism and SC with similar transition temperatures. However, there are obvious inconsistencies in the reported phase diagrams of Eu(Fe1-xCox)2As2 crystals grown by different methods. For transition metal arsenide (TMA)-flux-grown crystals, even the existence of SC is open for dispute. Here we re-examine the phase diagram of single-crystalline Eu(Fe1-xCox)2As2 grown by TMA flux. We found that the lattice parameter c shrinks linearly with Co doping, almost twice as fast as that of the tin-flux-grown crystals. With Co doping, the spin-density-wave (SDW) order of Fe sublattice is quickly suppressed, being detected only up to x = 0.08. The magnetic ordering temperature of the Eu2+ sublattice (TEu) shows a systematic evolution with Co doping, first going down and reaching a minimum at x = 0.08, then increasing continuously up to x = 0.24. Over the whole composition range investigated, no signature of SC is observed.


title: Influence of pulsatile blood flow on allometry of aortic wall shear stress

abstract: Shear stress plays an important role in the creation and evolution of atherosclerosis. An key element for in-vivo measurements and extrapolations is the dependence of shear stress on body mass. In the case of a Poiseuille modeling of the blood flow, P. Weinberg and C. Ethier have shown that shear stress on the aortic endothelium varies like body mass to the power $-\frac{3}{8}$, and is therefore 20-fold higher in mice than in men. However, by considering a more physiological oscillating Poiseuille + Womersley combinated flow in the aorta, we show that results differ notably: at larger masses ($M>10 \ kg$) shear stress varies as body mass to the power $-\frac{1}{8}$ and modifies the man to mouse ratio to 1:8. The allometry and values of temporal gradient of shear stress also change: $\partial\tau/\partial t$ varies as $M^{-3/8}$ instead of $M^{-5/8}$ at larger masses, and the 1:150 ratio from man to mouse becomes 1:61. Lastly, we show that the unsteady component of blood flow does not influence the constant allometry of peak velocity on body mass: $u_{max} \propto M^{0}$. This work extends our knowledge on the dependence of hemodynamic parameters on body mass and paves the way for a more precise extrapolation of in-vivo measurements to humans and bigger mammals.


title: Asymptotic uniformity of the quantization error for Moran measures on $\mathbb{R}^1$

abstract: Let $E$ be a Moran set on $\mathbb{R}^1$ associated with a closed interval $J$ and two sequences $(n_k)_{k=1}^\infty$ and $(\mathcal{C}_k=(c_{k,j})_{j=1}^{n_k})_{k\geq1}$. Let $\mu$ be the infinite product measure (Moran measure) on $E$ associated with a sequence $(\mathcal{P}_k)_{k\geq1}$ of positive probability vectors with $\mathcal{P}_k=(p_{k,j})_{j=1}^{n_k},k\geq 1$. We assume that \[ \inf_{k\geq1}\min_{1\leq j\leq n_k}c_{k,j}>0,\;\inf_{k\geq1}\min_{1\leq j\leq n_k}p_{k,j}>0. \] For every $n\geq 1$, let $\alpha_n$ be an $n$ optimal set in the quantization for $\mu$ of order $r\in(0,\infty)$ and $\{P_a(\alpha_n)\}_{a\in\alpha_n}$ an arbitrary Voronoi partition with respect to $\alpha_n$. For every $a\in\alpha_n$, we write $I_a(\alpha,\mu):=\int_{P_a(\alpha_n)}d(x,\alpha_n)^rd\mu(x)$ and \[ \underline{J}(\alpha_n,\mu):=\min_{a\in\alpha_n}I_a(\alpha,\mu),\; \overline{J}(\alpha_n,\mu):=\max_{a\in\alpha_n}I_a(\alpha,\mu). \] We show that $\underline{J}(\alpha_n,\mu),\overline{J}(\alpha_n,\mu)$ and $e^r_{n,r}(\mu)-e^r_{n+1,r}(\mu)$ are of the same order as $\frac{1}{n}e^r_{n,r}(\mu)$, where $e^r_{n,r}(\mu):=\int d(x,\alpha_n)^rd\mu(x)$ is the $n$th quantization error for $\mu$ of order $r$. In particular, for the class of Moran measures on $\mathbb{R}^1$, our result shows that a weaker version of Gersho's conjecture holds.


title: Photoacoustic Image Formation Based on Sparse Regularization of Minimum Variance Beamformer

abstract: Delay-and-Sum (DAS) is the most common algorithm used in photoacoustic (PA) image formation. However, this algorithm results in a reconstructed image with a wide mainlobe and high level of sidelobes. Minimum variance (MV), as an adaptive beamformer, overcomes these limitations and improves the image resolution and contrast. In this paper, a novel algorithm, named modified-sparse-MV (MS-MV) is proposed in which a L1-norm constraint is added to the MV minimization problem after some modifications, in order to suppress the sidelobes more efficiently, compared to MV. The added constraint can be interpreted as the sparsity of the output of the MV beamformed signals. Since the final minimization problem is convex, it can be solved efficiently using a simple iterative algorithm. The numerical results show that the proposed method, MS-MV beamformer, improves the signal-to-noise (SNR) about 19.48 dB, in average, compared to MV. Also, the experimental results, using a wire-target phantom, show that MS-MV leads to SNR improvement of about 2.64 dB in comparison with the MV.


title: A Generative Model for Dynamic Networks with Applications

abstract: Networks observed in real world like social networks, collaboration networks etc., exhibit temporal dynamics, i.e. nodes and edges appear and/or disappear over time. In this paper, we propose a generative, latent space based, statistical model for such networks (called dynamic networks). We consider the case where the number of nodes is fixed, but the presence of edges can vary over time. Our model allows the number of communities in the network to be different at different time steps. We use a neural network based methodology to perform approximate inference in the proposed model and its simplified version. Experiments done on synthetic and real world networks for the task of community detection and link prediction demonstrate the utility and effectiveness of our model as compared to other similar existing approaches.


title: Event Structures for Petri nets with Persistence

abstract: Event structures are a well-accepted model of concurrency. In a seminal paper by Nielsen, Plotkin and Winskel, they are used to establish a bridge between the theory of domains and the approach to concurrency proposed by Petri. A basic role is played by an unfolding construction that maps (safe) Petri nets into a subclass of event structures, called prime event structures, where each event has a uniquely determined set of causes. Prime event structures, in turn, can be identified with their domain of configurations. At a categorical level, this is nicely formalised by Winskel as a chain of coreflections. Contrary to prime event structures, general event structures allow for the presence of disjunctive causes, i.e., events can be enabled by distinct minimal sets of events. In this paper, we extend the connection between Petri nets and event structures in order to include disjunctive causes. In particular, we show that, at the level of nets, disjunctive causes are well accounted for by persistent places. These are places where tokens, once generated, can be used several times without being consumed and where multiple tokens are interpreted collectively, i.e., their histories are inessential. Generalising the work on ordinary nets, Petri nets with persistence are related to a new subclass of general event structures, called locally connected, by means of a chain of coreflections relying on an unfolding construction.


title: Separation choosability and dense bipartite induced subgraphs

abstract: We study a restricted form of list colouring, for which every pair of lists that correspond to adjacent vertices may not share more than one colour. The optimal list size such that a proper list colouring is always possible given this restriction, we call separation choosability. We show for bipartite graphs that separation choosability increases with (the logarithm of) the minimum degree. This strengthens results of Molloy and Thron and, partially, of Alon. One attempt to drop the bipartiteness assumption precipitates a natural class of Ramsey-type questions, of independent interest. For example, does every triangle-free graph of minimum degree $d$ contain a bipartite induced subgraph of minimum degree $\Omega(\log d)$ as $d\to\infty$?


title: Strongly suppressed superconducting proximity effect and ferromagnetism in trilayers of $\rm Bi_2Se_3$ / $\rm SrRuO_3$ / underdoped $\rm YBa_2Cu_3O_x$: A possible new platform for Majorana nano-electronics

abstract: We report properties of topological insulator - ferromagnet - superconductor trilayers comprised of thin films of 20 nm thick $\rm Bi_2Se_3$ on 10 nm $\rm SrRuO_3$ on 30 nm $\rm YBa_2Cu_3O_x$. As deposited trilayers are underdoped and have a superconductive transition with $\rm T_c$ onset at 75 K, zero resistance at 65 K, $\rm T_{Cueri}$ at 150 K and $\rm T^*$ of about 200 K. Further reannealing under vacuum yields the 60 K phase of $\rm YBa_2Cu_3O_x$ which still has zero resistance below about 40 K. Only when $10\times 100$ micro-bridges were patterned in the trilayer, some of the bridges showed resistive behavior all the way down to low temperatures. Magnetoresistance versus temperature of the superconductive ones showed the typical peak due to flux flow against pinning below $\rm T_c$, while the resistive ones showed only the broad leading edge of such a peak. All this indicates clearly weak-link superconductivity in the resistive bridges between superconductive $\rm YBa_2Cu_3O_x$ grains via the topological and ferromagnetic cap layers. Comparing our results to those of a reference trilayer with the topological $\rm Bi_2Se_3$ layer substituted by a non-superconducting highly overdoped $\rm La_{1.65}Sr_{0.35}CuO_4$, indicates that the superconductive proximity effect as well as ferromagnetism in the topological trilayer are actually strongly suppressed compared to the non-topological reference trilayer. This strong suppression is likely to originate in strong proximity induced edge currents in the SRO/YBCO layer that can lead to Majorana bound states, a possible signature of which is observed in the present study as zero bias conductance peaks.


title: The three point gauge algebra $\mathcal V\ltimes \mathfrak{sl}(2, \mathcal R) \oplus\left(\Omega_{\mathcal R}/d{\mathcal R}\right)$ and an action on a Fock space

abstract: The three point current algebra $\mathfrak{sl}(2,\mathcal R)$ where $\mathcal R=\mathbb C[t,t^{-1},u\,|\,u^2=t^2+4t ]$ and three-point Virasoro algebra both act on a previously constructed Fock space. In this paper we prove that the semi-direct product, i.e. the gauge algebra acts on the Fock space as well.


title: Rederiving the Upper Bound for Halving Edges using Cardano's Formula

abstract: In this paper we rederive an old upper bound on the number of halving edges present in the halving graph of an arbitrary set of $n$ points in 2-dimensions which are placed in general position. We provide a different analysis of an identity discovered by Andrejak et al, to rederive this upper bound of $O(n^{4/3})$. In the original paper of Andrejak et al. the proof is based on a naive analysis whereas in this paper we obtain the same upper bound by tightening the analysis thereby opening a new door to derive these upper bounds using the identity. Our analysis is based on a result of Cardano for finding the roots of a cubic equation. We believe that our technique has the potential to derive improved bounds on the number of halving edges.


title: Robust Private Information Retrieval from Coded Systems with Byzantine and Colluding Servers

abstract: A private information retrieval (PIR) scheme on coded storage systems with colluding, byzantine, and non-responsive servers is presented. Furthermore, the scheme can also be used for symmetric PIR in the same setting. An explicit scheme using an $[n,k]$ generalized Reed-Solomon storage code is designed, protecting against $t$-collusion and handling up to $b$ byzantine and $r$ non-responsive servers, when $n\geq n'= (\nu +1) k+t+2b+r-1$, for some integer $\nu \geq 1$. This scheme achieves a PIR rate of $1-\frac{k+2b+t+r-1}{n'}$. In the case where the capacity is known, namely when $k=1$, it is asymptotically capacity achieving as the number of files grows.


title: Nef anti-canonical divisors and rationally connected fibrations

abstract: We study the Iitaka-Kodaira dimension of nef relative anti-canonical divisors. As a consequence, we prove that given a complex projective variety with klt singularities, if the anti-canonical divisor is nef, then the dimension of a general fibre of the maximal rationally connected fibration is at least the Iitaka-Kodaira dimension of the anti-canonical divisor.


title: Zero limit of entropic relaxation time for the Shliomis model of ferrofluids

abstract: We construct solutions for the Shilomis model of ferrofluids in a critical space, uniformly in the entropic relaxation time $ \tau \in\left(0, \tau_0\right) $. This allows us to study the convergence when $ \tau\to 0 $ for such solutions.


title: The Use of Presence Data in Modelling Demand for Transportation

abstract: We consider the applicability of the data from operators of cellular systems to modelling demand for transportation. While individual-level data may contain precise paths of movement, stringent privacy rules prohibit their use without consent. Presence data aggregate the individual-level data to information on the numbers of transactions at each base transceiver station (BTS) per each time period. Our work is aimed at demonstrating value of such aggregate data for mobility management while maintaining privacy of users. In particular, given mobile subscriber activity aggregated to short time intervals for a zone, a convex optimisation problem estimates most likely transitions between zones. We demonstrate the method on presence data from Warsaw, Poland, and compare with official demand estimates obtained with classical econometric methods.


title: Structural Estimation of Behavioral Heterogeneity

abstract: We develop a behavioral asset pricing model in which agents trade in a market with information friction. Profit-maximizing agents switch between trading strategies in response to dynamic market conditions. Due to noisy private information about the fundamental value, the agents form different evaluations about heterogeneous strategies. We exploit a thin set---a small sub-population---to pointly identify this nonlinear model, and estimate the structural parameters using extended method of moments. Based on the estimated parameters, the model produces return time series that emulate the moments of the real data. These results are robust across different sample periods and estimation methods.


title: On a Riemannian manifold with a circulant structure whose third power is the identity

abstract: It is studied a 3-dimensional Riemannian manifold equipped with a tensor structure of type (1,1), whose third power is the identity. This structure has a circulant matrix with respect to some basis, i.e. the structure is circulant. On such a manifold a fundamental tensor by the metric and by the covariant derivative of the circulant structure is defined. An important characteristic identity for this tensor is obtained. It is established that the image of the fundamental tensor with respect to the usual conformal transformation satisfies the same identity. A Lie group as a manifold of the considered type is constructed and some of its geometrical characteristics are found.


title: Probing shear-induced rearrangements in Fourier Space. I. Dynamic Light Scattering

abstract: Understanding the microscopic origin of the rheological behavior of soft matter is a long-lasting endeavour. While early efforts concentrated mainly on the relationship between rheology and structure, current research focuses on the role of microscopic dynamics. We present in two companion papers a thorough discussion of how Fourier space-based methods may be coupled to rheology to shed light on the relationship between the microscopic dynamics and the mechanical response of soft systems. In this first companion paper, we report a theoretical, numerical and experimental investigation of dynamic light scattering coupled to rheology. While in ideal solids and simple viscous fluids the displacement field under a shear deformation is purely affine, additional non-affine displacements arise in many situations of great interest, for example in elastically heterogeneous materials or due to plastic rearrangements. We show how affine and non-affine displacements can be separately resolved by dynamic light scattering, and discuss in detail the effect of several non-idealities in typical experiments.


title: Efficient Machine Learning Representations of Surface Code with Boundaries, Defects, Domain Walls and Twists

abstract: Machine learning representations of many-body quantum states have recently been introduced as an ansatz to describe the ground states and unitary evolutions of many-body quantum systems. We explore one of the most important representations, restricted Boltzmann machine (RBM) representation, in stabilizer formalism. We give the general method of constructing RBM representation for stabilizer code states and find the exact RBM representation for several types of stabilizer groups with the number of hidden neurons equal or less than the number of visible neurons, which indicates that the representation is extremely efficient. Then we analyze the surface code with boundaries, defects, domain walls and twists in full detail and find that all the models can be efficiently represented via RBM ansatz states. Besides, the case for Kitaev's $D(\Zb_d)$ model, which is a generalized model of surface code, is also investigated.


title: Dark matter: an efficient catalyst for intermediate-mass-ratio-inspiral events

abstract: Gravitational waves (GWs) can be produced if a stellar compact object, such as a black hole (BH) or neutron star, inspirals into an intermediate-massive black hole (IMBH) of $(10^3 \sim 10^5)\,M_\odot$. Such a system may be produced in the center of a globular cluster (GC) or a nuclear star cluster (NSC), and is known as an intermediate- or extreme-mass-ratio inspiral (IMRI or EMRI). Motivated by the recent suggestions that dark matter minispikes could form around IMBHs, we study the effect of dynamical friction against DM on the merger rate of IMRIs/EMRIs. We find that the merger timescale of IMBHs with BHs and NSs would be shortened by two to three orders of magnitude. As a result, the event rate of IMRIs/EMRIs are enhanced by orders of magnitude relative to that in the case of no DM minispikes. In the most extreme case where IMBHs are small and the DM minispikes have a steep density profile, all the BH in GCs and NSCs might be exhausted so that the mergers with NSs would dominate the current IMRIs/EMRIs. Our results suggest that the mass function of the IMBHs below $10^4 \,M_\odot$ would bear imprints of the distribution of DM minispikes because these low-mass IMBHs can grow efficiently in the presence of DM minispikes by merging with BHs and NSs. Future space-based GW detectors, like LISA, Taiji, and Tianqin, can measure the IMRI/EMRI rate and hence constrain the distribution of DM around IMBHs.


title: Metamaterial analogue of Ising model

abstract: The interaction between microscopic particles has always been a fascinating and intriguing area of science. Direct interrogation of such interactions is often difficult or impossible. Structured electromagnetic systems offer a rich toolkit for mimicking and reproducing the key dynamics that governs the microscopic interactions, and thus provide an avenue to explore and interpret the microscopic phenomena. In particular, metamaterials offer the freedom to artificially tailor light-matter coupling and to control the interaction between unit cells in the metamaterial array. Here we demonstrate a terahertz metamaterial that mimics spin-related interactions of microscopic particles in a 2D lattice via complex electromagnetic multipole interactions within the metamaterial array. Fano resonances featured by distinct mode properties due to strong nearest-neighbor interactions are discussed that draw parallels with the 2D Ising model. Interestingly, a hyperfine Fano splitting spectrum is observed by manipulating the 2D interactions without applying external magnetic or electric fields, which provides a passive multispectral platform for applications in super-resolution imaging, biosensing, and selective thermal emission. The dynamic approach to reproduce the static interaction between microscopic particles would enable more profound significance in exploring the unknown physical world by the macroscopic analogues.


title: Power-law/exponential transport of electromagnetic field in one-dimensional metallic nanoparticle arrays

abstract: Based on the coupled-dipole analysis and finite-difference time-domain simulation, we have investigated the surface plasmon propagation in one-dimensional metallic nanoparticle (NP) chains. Our systematic studies reveal that the interplay between the localized plasmon excitation and the lattice collective behavior leads to two phases (I and II) of different electromagnetic (EM) field transport properties. In phase I, the EM field decays follow the power-law. In phase II, the EM field shows the exponential decay in the short distance regime and the power-law decay in the long distance regime. Moreover, universal power-law exponents have been found in the long propagation distance. The two different EM field propagation behaviors can be transformed to each other by tuning the parameters of the excitation fields and/or those of the NP chains. The EM field transport mechanisms we have found are very useful in the design of plasmonic waveguide with both strong field confinement and efficient field/energy transfer, which has important applications in integrated nanophotonic circuits.


title: On a linearization trick

abstract: In several situations, mainly involving a self-adjoint set of unitary generators of a $C^*$-algebra, we show that any matrix polynomial in the generators and the unit that is in the open unit ball can be written as a product of matrix polynomials of degree 1 also in the open unit ball.


title: Morrey meets Muckenhoupt: A note on Nakai's generalized Morrey spaces and applications

abstract: The goal of this paper is to extend Nakai's generalized Morrey spaces to a wider function class, the one-sided Muckenhoupt weighted case. Morrey matching Muckenhoupt enables us to study the weak and strong type boundedness of one-sided sublinear operators satisfying certain size conditions on the one-sided weighted Morrey spaces. We also establish one-sided Fefferman-Stein inequalities on one-sided weighted Morrey spaces in this paper. Meanwhile, the boundedness and compactness of Riemann-Liouville integral operators on locally one-sided weighted Morrey space are considered. As applications, we establish the existence and uniqueness of solutions to a Cauchy type problem associated with fractional differential equations.


title: Denotational semantics for guarded dependent type theory

abstract: We present a new model of Guarded Dependent Type Theory (GDTT), a type theory with guarded recursion and multiple clocks in which one can program with, and reason about coinductive types. Productivity of recursively defined coinductive programs and proofs is encoded in types using guarded recursion, and can therefore be checked modularly, unlike the syntactic checks implemented in modern proof assistants. The model is based on a category of covariant presheaves over a category of time objects, and quantification over clocks is modelled using a presheaf of clocks. To model the clock irrelevance axiom, crucial for programming with coinductive types, types must be interpreted as presheaves orthogonal to the object of clocks. In the case of dependent types, this translates to a lifting condition similar to the one found in homotopy theoretic models of type theory, but here with an additional requirement of uniqueness of lifts. Since the universes defined by the standard Hofmann-Streicher construction in this model do not satisfy this property, the universes in GDTT must be indexed by contexts of clock variables. We show how to model these universes in such a way that inclusions of clock contexts give rise to inclusions of universes commuting with type operations on the nose.


title: Geant4 simulation of the ELIMED transport and dosimetry beam line for high-energy laser-driven ion beam multidisciplinary applications

abstract: The ELIMED (MEDical and multidisciplinary application at ELI Beamlines) beam line is being developed at INFN-LNS with the aim of transporting and selecting in energy proton and ion beams accelerated by laser-matter interaction at ELI Beamlines in Prague. It will be a section of the ELIMAIA (ELI Multidisciplinary Applications of laser-Ions Acceleration) beam line, dedicated to applications, including the medical one, of laser-accelerated ion beams [1,2]. A Monte Carlo model has been developed to support the design of the beam line in terms of particle transport efficiency, to optimize the transport parameters at the irradiation point in air and, furthermore, to predict beam parameters in order to deliver dose distributions of clinical relevance. The application has been developed using the Geant4 [3] Monte Carlo toolkit and has been designed in a modular way in order to easily switch on/off geometrical components according to different experimental setups and users requirements, as reported in [4], describing the early-stage code and simulations. The application has been delivered to ELI Beamlines and will be available for future ELIMAIA's users as ready-to-use tool useful during experiment preparation and analysis. The final version of the developed application will be described in detail in this contribution, together with the final results, in terms of energy spectra and transmission efficiency along the in-vacuum beam line, obtained by performing end-to-end simulations.


title: Conserving symmetries in Bose-Einstein condensate dynamics requires many-body theory

abstract: We explain from first principles why satisfying conservation laws in Bose Einstein condensate dynamics requires many-body theory. For the Gross-Pitaevskii mean-field we show analytically and numerically that conservation laws are violated. We provide examples for angular momentum and linear momentum conservation. Arbitrarily large violations occur despite negligible depletion and interaction energy. For the case of angular momentum we show through extensive many-body simulations how the conservation law can be gradually restored on the many-body level. Implications are discussed.


title: Bounds for the first non-zero Steklov eigenvalue

abstract: Let $\Omega$ be a star-shaped bounded domain in $(\mathbb{S}^{n}, ds^{2})$ with smooth boundary. In this article, we give a sharp lower bound for the first non-zero eigenvalue of the Steklov eigenvalue problem in $\Omega.$ This result is the generalization of a result given by Kuttler and Sigillito for a star-shaped bounded domain in $\mathbb{R}^2.$ Further we also obtain a two sided bound for the first non-zero eigenvalue of the Steklov problem on the ball in $\mathbb{R}^n$ with rotationally invariant metric and with bounded radial curvature.


title: Binary Pebbling Algorithms for In-Place Reversal of One-Way Hash Chains

abstract: We present optimal binary pebbling algorithms for in-place reversal (backward traversal) of one-way hash chains. For a hash chain of length $2^k$, the number of hashes performed in each output round does not exceed $\lceil k/2 \rceil$, whereas the number of hash values stored (pebbles) throughout is at most $k$. We introduce a framework for rigorous comparison of explicit binary pebbling algorithms, including simple speed-1 binary pebbling, Jakobsson's speed-2 binary pebbling, and our optimal binary pebbling algorithm. Explicit schedules describe for each pebble exactly how many hashes need to be performed in each round. The optimal schedule turns out to be essentially unique and exhibits a nice recursive structure, which allows for fully optimized implementations that can readily be deployed. In particular, we develop the first in-place implementations with minimal storage overhead (essentially, storing only hash values), and fast implementations with minimal computational overhead.


title: Locality Optimized Unstructured Mesh Algorithms on GPUs

abstract: Unstructured-mesh based numerical algorithms such as finite volume and finite element algorithms form an important class of applications for many scientific and engineering domains. The key difficulty in achieving higher performance from these applications is the indirect accesses that lead to data-races when parallelized. Current methods for handling such data-races lead to reduced parallelism and suboptimal performance. Particularly on modern many-core architectures, such as GPUs, that has increasing core/thread counts, reducing data movement and exploiting memory locality is vital for gaining good performance. In this work we present novel locality-exploiting optimizations for the efficient execution of unstructured-mesh algorithms on GPUs. Building on a two-layered coloring strategy for handling data races, we introduce novel reordering and partitioning techniques to further improve efficient execution. The new optimizations are then applied to several well established unstructured-mesh applications, investigating their performance on NVIDIA's latest P100 and V100 GPUs. We demonstrate significant speedups ($1.1\text{--}1.75\times$) compared to the state-of-the-art. A range of performance metrics are benchmarked including runtime, memory transactions, achieved bandwidth performance, GPU occupancy and data reuse factors and are used to understand and explain the key factors impacting performance. The optimized algorithms are implemented as an open-source software library and we illustrate its use for improving performance of existing or new unstructured-mesh applications.


title: FD-MobileNet: Improved MobileNet with a Fast Downsampling Strategy

abstract: We present Fast-Downsampling MobileNet (FD-MobileNet), an efficient and accurate network for very limited computational budgets (e.g., 10-140 MFLOPs). Our key idea is applying an aggressive downsampling strategy to MobileNet framework. In FD-MobileNet, we perform 32$\times$ downsampling within 12 layers, only half the layers in the original MobileNet. This design brings three advantages: (i) It remarkably reduces the computational cost. (ii) It increases the information capacity and achieves significant performance improvements. (iii) It is engineering-friendly and provides fast actual inference speed. Experiments on ILSVRC 2012 and PASCAL VOC 2007 datasets demonstrate that FD-MobileNet consistently outperforms MobileNet and achieves comparable results with ShuffleNet under different computational budgets, for instance, surpassing MobileNet by 5.5% on the ILSVRC 2012 top-1 accuracy and 3.6% on the VOC 2007 mAP under a complexity of 12 MFLOPs. On an ARM-based device, FD-MobileNet achieves 1.11$\times$ inference speedup over MobileNet and 1.82$\times$ over ShuffleNet under the same complexity.


title: Artificial intelligence meets minority game: toward optimal resource allocation

abstract: Resource allocation systems provide the fundamental support for the normal functioning and well being of the modern society, and can be modeled as minority games. A ubiquitous dynamical phenomenon is the emergence of herding, where a vast majority of the users concentrate on a small number of resources, leading to a low efficiency in resource allocation. To devise strategies to prevent herding is thus of high interest. Previous works focused on control strategies that rely on external interventions, such as pinning control where a fraction of users are forced to choose a certain action. Is it possible to eliminate herding without any external control? The main point of this paper is to provide an affirmative answer through exploiting artificial intelligence (AI). In particular, we demonstrate that, when agents are empowered with reinforced learning in that they get familiar with the unknown game environment gradually and attempt to deliver the optimal actions to maximize the payoff, herding can effectively be eliminated. Computations reveal the striking phenomenon that, regardless of the initial state, the system evolves persistently and relentlessly toward the optimal state in which all resources are used efficiently. However, the evolution process is not without interruptions: there are large fluctuations that occur but only intermittently in time. The statistical distribution of the time between two successive fluctuating events is found to depend on the parity of the evolution, i.e., whether the number of time steps in between is odd or even. We develop a physical analysis and derive mean-field equations to gain an understanding of these phenomena. As minority game dynamics and the phenomenon of herding are common in social, economic, and political systems, and since AI is becoming increasingly widespread, we expect our AI empowered minority game system to have broad applications.


title: Supervised classification of Dermatological diseases by Deep learning

abstract: This paper introduces a deep-learning based efficient classifier for common dermatological conditions, aimed at people without easy access to skin specialists. We report approximately 80% accuracy, in a situation where primary care doctors have attained 57% success rate, according to recent literature. The rationale of its design is centered on deploying and updating it on handheld devices in near future. Dermatological diseases are common in every population and have a wide spectrum in severity. With a shortage of dermatological expertise being observed in several countries, machine learning solutions can augment medical services and advise regarding existence of common diseases. The paper implements supervised classification of nine distinct conditions which have high occurrence in East Asian countries. Our current attempt establishes that deep learning based techniques are viable avenues for preliminary information to aid patients.


title: Sample Efficient Deep Reinforcement Learning for Dialogue Systems with Large Action Spaces

abstract: In spoken dialogue systems, we aim to deploy artificial intelligence to build automated dialogue agents that can converse with humans. A part of this effort is the policy optimisation task, which attempts to find a policy describing how to respond to humans, in the form of a function taking the current state of the dialogue and returning the response of the system. In this paper, we investigate deep reinforcement learning approaches to solve this problem. Particular attention is given to actor-critic methods, off-policy reinforcement learning with experience replay, and various methods aimed at reducing the bias and variance of estimators. When combined, these methods result in the previously proposed ACER algorithm that gave competitive results in gaming environments. These environments however are fully observable and have a relatively small action set so in this paper we examine the application of ACER to dialogue policy optimisation. We show that this method beats the current state-of-the-art in deep learning approaches for spoken dialogue systems. This not only leads to a more sample efficient algorithm that can train faster, but also allows us to apply the algorithm in more difficult environments than before. We thus experiment with learning in a very large action space, which has two orders of magnitude more actions than previously considered. We find that ACER trains significantly faster than the current state-of-the-art.


title: Partial immunization of trees

abstract: For a graph $G$ and an integer-valued function $\tau$ on its vertex set, a dynamic monopoly is a set of vertices of $G$ such that iteratively adding to it vertices $u$ of $G$ that have at least $\tau(u)$ neighbors in it eventually yields the vertex set of $G$. We study the problem of maximizing the minimum order of a dynamic monopoly by increasing the threshold values of individual vertices subject to vertex-dependent lower and upper bounds, and fixing the total increase. We solve this problem efficiently for trees, which extends a result of Khoshkhah and Zaker (On the largest dynamic monopolies of graphs with a given average threshold, Canadian Mathematical Bulletin 58 (2015) 306-316).