title: Representation, Exploration and Recommendation of Music Playlists

abstract: Playlists have become a significant part of our listening experience because of the digital cloud-based services such as Spotify, Pandora, Apple Music. Owing to the meteoric rise in the usage of playlists, recommending playlists is crucial to music services today. Although there has been a lot of work done in playlist prediction, the area of playlist representation hasn't received that level of attention. Over the last few years, sequence-to-sequence models, especially in the field of natural language processing, have shown the effectiveness of learned embeddings in capturing the semantic characteristics of sequences. We can apply similar concepts to music to learn fixed length representations for playlists and use those representations for downstream tasks such as playlist discovery, browsing, and recommendation. In this work, we formulate the problem of learning a fixed-length playlist representation in an unsupervised manner, using Sequence-to-sequence (Seq2seq) models, interpreting playlists as sentences and songs as words. We compare our model with two other encoding architectures for baseline comparison. We evaluate our work using the suite of tasks commonly used for assessing sentence embeddings, along with a few additional tasks pertaining to music, and a recommendation task to study the traits captured by the playlist embeddings and their effectiveness for the purpose of music recommendation.


title: Predicting Treatment Initiation from Clinical Time Series Data via Graph-Augmented Time-Sensitive Model

abstract: Many computational models were proposed to extract temporal patterns from clinical time series for each patient and among patient group for predictive healthcare. However, the common relations among patients (e.g., share the same doctor) were rarely considered. In this paper, we represent patients and clinicians relations by bipartite graphs addressing for example from whom a patient get a diagnosis. We then solve for the top eigenvectors of the graph Laplacian, and include the eigenvectors as latent representations of the similarity between patient-clinician pairs into a time-sensitive prediction model. We conducted experiments using real-world data to predict the initiation of first-line treatment for Chronic Lymphocytic Leukemia (CLL) patients. Results show that relational similarity can improve prediction over multiple baselines, for example a 5% incremental over long-short term memory baseline in terms of area under precision-recall curve.


title: A comment on "Can observations inside the Solar System reveal the gravitational properties of the quantum vacuum?" by D.S. Hajdukovic

abstract: The modified gravitational theory by Hajdukovic [arXiv:1210.7421], based on the idea that quantum vacuum contains virtual gravitational dipoles, predicts, among other things, anomalous secular precessions of the planets of the Solar System as large as $\simeq 700-6,000$ milliarceconds per century. We demonstrate that they are ruled out by several orders of magnitude by the existing bounds on any anomalous orbital secular rates obtained with the EPM and INPOP ephemerides.


title: A Simulation Study of Social-Networking-Driven Smart Recommendations for Internet of Vehicles

abstract: Social aspects of connectivity and information dispersion are often ignored while weighing the potential of Internet of Things (IoT). In the specialized domain of Internet of Vehicles (IoV), Social IoV (SIoV) is introduced realization its importance. Assuming a more commonly acceptable standardization of Big Data generated by IoV, the social dimensions enabling its fruitful usage remains a challenge. In this paper, an agent-based model of information sharing between vehicles for context-aware recommendations is presented. The model adheres to social dimensions as that of human society. Some important hypotheses are tested under reasonable connectivity and data constraints. The simulation results reveal that closure of social ties and its timing impacts dispersion of novel information (necessary for a recommender system) substantially. It was also observed that as the network evolves as a result of incremental interactions, recommendations guaranteeing a fair distribution of vehicles across equally good competitors is not possible.


title: Nature Inspired Dimensional Reduction Technique for Fast and Invariant Visual Feature Extraction

abstract: Fast and invariant feature extraction is crucial in certain computer vision applications where the computation time is constrained in both training and testing phases of the classifier. In this paper, we propose a nature-inspired dimensionality reduction technique for fast and invariant visual feature extraction. The human brain can exchange the spatial and spectral resolution to reconstruct missing colors in visual perception. The phenomenon is widely used in the printing industry to reduce the number of colors used to print, through a technique, called color dithering. In this work, we adopt a fast error-diffusion color dithering algorithm to reduce the spectral resolution and extract salient features by employing novel Hessian matrix analysis technique, which is then described by a spatial-chromatic histogram. The computation time, descriptor dimensionality and classification performance of the proposed feature are assessed under drastic variances in orientation, viewing angle and illumination of objects comparing with several different state-of-the-art handcrafted and deep-learned features. Extensive experiments on two publicly available object datasets, coil-100 and ALOI carried on both a desktop PC and a Raspberry Pi device show multiple advantages of using the proposed approach, such as the lower computation time, high robustness, and comparable classification accuracy under weakly supervised environment. Further, it showed the capability of operating solely inside a conventional SoC device utilizing a small fraction of the available hardware resources.


title: Molecular activity prediction using graph convolutional deep neural network considering distance on a molecular graph

abstract: Machine learning is often used in virtual screening to find compounds that are pharmacologically active on a target protein. The weave module is a type of graph convolutional deep neural network that uses not only features focusing on atoms alone (atom features) but also features focusing on atom pairs (pair features); thus, it can consider information of nonadjacent atoms. However, the correlation between the distance on the graph and the three-dimensional coordinate distance is uncertain. In this paper, we propose three improvements for modifying the weave module. First, the distances between ring atoms on the graph were modified to bring the distances on the graph closer to the coordinate distance. Second, different weight matrices were used depending on the distance on the graph in the convolution layers of the pair features. Finally, a weighted sum, by distance, was used when converting pair features to atom features. The experimental results show that the performance of the proposed method is slightly better than that of the weave module, and the improvement in the distance representation might be useful for compound activity prediction.


title: Isolation Kernel: The X Factor in Efficient and Effective Large Scale Online Kernel Learning

abstract: Large scale online kernel learning aims to build an efficient and scalable kernel-based predictive model incrementally from a sequence of potentially infinite data points. A current key approach focuses on ways to produce an approximate finite-dimensional feature map, assuming that the kernel used has a feature map with intractable dimensionality---an assumption traditionally held in kernel-based methods. While this approach can deal with large scale datasets efficiently, this outcome is achieved by compromising predictive accuracy because of the approximation. We offer an alternative approach which overrides the assumption and puts the kernel used at the heart of the approach. It focuses on creating an exact, sparse and finite-dimensional feature map of a kernel called Isolation Kernel. Using this new approach, to achieve the above aim of large scale online kernel learning becomes extremely simple---simply use Isolation Kernel instead of a kernel having a feature map with intractable dimensionality. We show that, using Isolation Kernel, large scale online kernel learning can be achieved efficiently without sacrificing accuracy.


title: Energy conservative SBP discretizations of the acoustic wave equation in covariant form on staggered curvilinear grids

abstract: We develop a numerical method for solving the acoustic wave equation in covariant form on staggered curvilinear grids in an energy conserving manner. The use of a covariant basis decomposition leads to a rotationally invariant scheme that outperforms a Cartesian basis decomposition on rotated grids. The discretization is based on high order Summation-By-Parts (SBP) operators and preserves both symmetry and positive definiteness of the contravariant metric tensor. To improve accuracy and decrease computational cost, we also derive a modified discretization of the metric tensor that leads to a conditionally stable discretization. Bounds are derived that yield a point-wise condition that can be evaluated to check for stability of the modified discretization. This condition shows that the interpolation operators should be constructed such that their norm is close to one.


title: An efficient method to solve the mathematical model of HIV infection for CD8+ T-cells

abstract: In this paper, the mathematical model of HIV infection for CD8+ T-cells is illustrated. The homotopy analysis method and the Laplace transformations are combined for solving this model. Also, the convergence theorem is proved to demonstrate the abilities of presented method for solving non-linear mathematical models. The numerical results for N = 5, 10 are presented. Several h-curves are plotted to show the convergence regions of solutions. The plots of residual error functions indicate the precision of presented method.


title: The fourth moment of quadratic Dirichlet $L$-functions

abstract: We study the fourth moment of quadratic Dirichlet $L$-functions at $s= \frac{1}{2}$. We show an asymptotic formula under the generalized Riemann hypothesis, and obtain a precise lower bound unconditionally. The proofs of these results follow closely arguments of Soundararajan and Young [19] and Soundararajan [17].


title: Language2Pose: Natural Language Grounded Pose Forecasting

abstract: Generating animations from natural language sentences finds its applications in a a number of domains such as movie script visualization, virtual human animation and, robot motion planning. These sentences can describe different kinds of actions, speeds and direction of these actions, and possibly a target destination. The core modeling challenge in this language-to-pose application is how to map linguistic concepts to motion animations. In this paper, we address this multimodal problem by introducing a neural architecture called Joint Language to Pose (or JL2P), which learns a joint embedding of language and pose. This joint embedding space is learned end-to-end using a curriculum learning approach which emphasizes shorter and easier sequences first before moving to longer and harder ones. We evaluate our proposed model on a publicly available corpus of 3D pose data and human-annotated sentences. Both objective metrics and human judgment evaluation confirm that our proposed approach is able to generate more accurate animations and are deemed visually more representative by humans than other data driven approaches.


title: Stochastic Lagrangians for Statistical Dynamics

abstract: The concept of stochastic Lagrangian and its use in statistical dynamics is illustrated theoretically, and with some examples. Dynamical variables undergoing stochastic differential equations are stochastic processes themselves, and their realization probability functional within a given time interval arises from the interplay between the deterministic parts of dynamics and noise statistics. The stochastic Lagrangian is a tool to formulate realization probabilities via functional integrals, once the statistics of noises involved in the stochastic dynamical equations is known. In principle, it allows to highlight the invariance properties of the statistical dynamics of the system. In this work, after a review of the stochastic Lagrangian formalism, some applications of it to physically relevant cases are illustrated.


title: Robust analogs to the Coefficient of Variation

abstract: The coefficient of variation (CV) is commonly used to measure relative dispersion. However, since it is based on the sample mean and standard deviation, outliers can adversely affect the CV. Additionally, for skewed distributions the mean and standard deviation do not have natural interpretations and, consequently, neither does the CV. Here we investigate the extent to which quantile-based measures of relative dispersion can provide appropriate summary information as an alternative to the CV. In particular, we investigate two measures, the first being the interquartile range (in lieu of the standard deviation), divided by the median (in lieu of the mean), and the second being the median absolute deviation (MAD), divided by the median, as robust estimators of relative dispersion. In addition to comparing the influence functions of the competing estimators and their asymptotic biases and variances, we compare interval estimators using simulation studies to assess coverage.


title: How trapped particles interact with and sample superfluid vortex excitations

abstract: Particles have been used for more than a decade to visualize and study the dynamics of quantum vortices in superfluid helium. In this work we study how the dynamics of a collection of particles set inside a vortex reflects the motion of the vortex. We use a self-consistent model based on the Gross-Pitaevskii equation coupled with classical particle dynamics. We find that each particle oscillates with a natural frequency proportional to the number of vortices attached to it. We then study the dynamics of an array of particles trapped in a quantum vortex and use particle trajectories to measure the frequency spectrum of the vortex excitations. Surprisingly, due to the discreetness of the array, the vortex excitations measured by the particles exhibits bands, gaps and Brillouin zones, analogous to the ones of electrons moving in crystals. We then establish a mathematical analogy where vortex excitations play the role of electrons and particles that of the potential barriers constituting the crystal. We find that the height of the effective potential barriers is proportional to the particle mass and the frequency of the incoming waves. We conclude that large-scale vortex excitations could be in principle directly measured by particles and novel physics could emerge from particle-vortex interaction.


title: On the Optimal Refresh Power Allocation for Energy-Efficient Memories

abstract: Refresh is an important operation to prevent loss of data in dynamic random-access memory (DRAM). However, frequent refresh operations incur considerable power consumption and degrade system performance. Refresh power cost is especially significant in high-capacity memory devices and battery-powered edge/mobile applications. In this paper, we propose a principled approach to optimizing the refresh power allocation. Given a model for the bit error rate dependence on power, we formulate a convex optimization problem to minimize the word mean squared error for a refresh power constraint; hence we can guarantee the optimality of the obtained refresh power allocations. In addition, we provide an integer programming problem to optimize the discrete refresh interval assignments. For an 8-bit accessed word, numerical results show that the optimized nonuniform refresh intervals reduce the refresh power by 29% at a peak signal-to-noise ratio of 50dB compared to the uniform assignment.


title: Robust Tensor Completion Using Transformed Tensor SVD

abstract: In this paper, we study robust tensor completion by using transformed tensor singular value decomposition (SVD), which employs unitary transform matrices instead of discrete Fourier transform matrix that is used in the traditional tensor SVD. The main motivation is that a lower tubal rank tensor can be obtained by using other unitary transform matrices than that by using discrete Fourier transform matrix. This would be more effective for robust tensor completion. Experimental results for hyperspectral, video and face datasets have shown that the recovery performance for the robust tensor completion problem by using transformed tensor SVD is better in PSNR than that by using Fourier transform and other robust tensor completion methods.


title: A Minimizing Valuation is Quasi-monomial

abstract: We prove a version of Jonsson-Musta\c{t}\v{a}'s Conjecture, which says for any graded sequence of ideals, there exists a quasi-monomial valuation computing its log canonical threshold. As a corollary, we confirm Chi Li's conjecture that a minimizer of the normalized volume function is always quasi-monomial. Applying our techniques to a family of klt singularities, we show that the volume of klt singularities is a constructible function. As a corollary, we prove that in a family of klt log Fano pairs, the K-semistable ones form a Zariski open set. Together with [Jia17], we conclude that all K-semistable klt Fano varieties with a fixed dimension and volume are parametrized by an Artin stack of finite type, which then admits a separated good moduli space by [BX18, ABHLX19], whose geometric points parametrize K-polystable klt Fano varieties.


title: Neural Semantic Parsing with Anonymization for Command Understanding in General-Purpose Service Robots

abstract: Service robots are envisioned to undertake a wide range of tasks at the request of users. Semantic parsing is one way to convert natural language commands given to these robots into executable representations. Methods for creating semantic parsers, however, rely either on large amounts of data or on engineered lexical features and parsing rules, which has limited their application in robotics. To address this challenge, we propose an approach that leverages neural semantic parsing methods in combination with contextual word embeddings to enable the training of a semantic parser with little data and without domain specific parser engineering. Key to our approach is the use of an anonymized target representation which is more easily learned by the parser. In most cases, this simplified representation can trivially be transformed into an executable format, and in others the parse can be completed through further interaction with the user. We evaluate this approach in the context of the RoboCup@Home General Purpose Service Robot task, where we have collected a corpus of paraphrased versions of commands from the standardized command generator. Our results show that neural semantic parsers can predict the logical form of unseen commands with 89% accuracy. We release our data and the details of our models to encourage further development from the RoboCup and service robotics communities.


title: Counting bounded elements of a number field

abstract: We estimate, in a number field, the number of elements and the maximal number of linearly independent elements, with prescribed bounds on their valuations. As a by-product, we obtain new bounds for the successive minima of ideal lattices. Our arguments combine group theory, ramification theory, and the geometry of numbers.


title: Exploring Feasible Design Spaces for Heterogeneous Constraints

abstract: We demonstrate an approach of exploring design spaces to simultaneously satisfy kinematics- and physics-based requirements. We present a classification of constraints and solvers to enable postponing optimization as far down the design workflow as possible. The solvers are organized into two broad classes of design space 'pruning' and 'exploration' by considering the types of constraints they can satisfy. We show that pointwise constraints define feasible design subspaces that can be represented and computed as first-class entities by their maximal feasible elements. The design space is pruned upfront by intersecting maximal elements, without premature optimization. To solve for other constraints, we apply topology optimization (TO), starting from the pruned feasible space. The optimization is steered by a topological sensitivity field (TSF) that measures the global changes in violation of constraints with respect to local topological punctures. The TSF for global objective functions is augmented with TSF for global constraints, and penalized/filtered to incorporate local constraints, including set constraints converted to differentiable (in)equality constraints. We demonstrate application of the proposed workflow to nontrivial examples in design and manufacturing. Among other examples, we show how to explore pruned design spaces via TO to simultaneously satisfy physics-based constraints (e.g., minimize compliance and mass) as well as kinematics-based constraints (e.g., maximize accessibility for machining).


title: Neural Machine Reading Comprehension: Methods and Trends

abstract: Machine reading comprehension (MRC), which requires a machine to answer questions based on a given context, has attracted increasing attention with the incorporation of various deep-learning techniques over the past few years. Although research on MRC based on deep learning is flourishing, there remains a lack of a comprehensive survey summarizing existing approaches and recent trends, which motivated the work presented in this article. Specifically, we give a thorough review of this research field, covering different aspects including (1) typical MRC tasks: their definitions, differences, and representative datasets; (2) the general architecture of neural MRC: the main modules and prevalent approaches to each; and (3) new trends: some emerging areas in neural MRC as well as the corresponding challenges. Finally, considering what has been achieved so far, the survey also envisages what the future may hold by discussing the open issues left to be addressed.


title: Comparative analysis of layered structures in empirical investor networks and cellphone communication networks

abstract: Empirical investor networks (EIN) proposed by \cite{Ozsoylev-Walden-Yavuz-Bildik-2014-RFS} are assumed to capture the information spreading path among investors. Here, we perform a comparative analysis between the EIN and the cellphone communication networks (CN) to test whether EIN is an information exchanging network from the perspective of the layer structures of ego networks. We employ two clustering algorithms ($k$-means algorithm and $H/T$ break algorithm) to detect the layer structures for each node in both networks. We find that the nodes in both networks can be clustered into two groups, one that has a layer structure similar to the theoretical Dunbar Circle corresponding to that the alters in ego networks exhibit a four-layer hierarchical structure with the cumulative number of 5, 15, 50 and 150 from the inner layer to the outer layer, and the other one having an additional inner layer with about 2 alters compared with the Dunbar Circle. We also find that the scale ratios, which are estimated based on the unique parameters in the theoretical model of layer structures \citep{Tamarit-Cuesta-Dunbar-Sanchez-2018-PNAS}, conform to a log-normal distribution for both networks. Our results not only deepen our understanding on the topological structures of EIN, but also provide empirical evidence of the channels of information diffusion among investors.


title: Topological Superconductivity in Doped Symmetry Protected Topological State

abstract: We propose an exotic scenario that topological superconductivity can emerge by doping strongly interacting fermionic systems whose spin degrees of freedom form bosonic symmetry protected topological (SPT) state. Specifically, we study a 1-dimensional (1D) example where the spin degrees of freedom form a spin-1 Haldane phase. {Before doping, the charge and spin degrees of freedom are both gapped.} Upon doping, {the charge channel becomes gapless and is described by a $c=1$ compactified bosonic conformal field theory (CFT), while the spin channel remains gapped and still form a bosonic SPT state. Interestingly,} an instability toward $p$-wave topological superconductivity is induced coexisting with the symmetry protected spin edge modes that are inherited from the Haldane phase. This scenario is confirmed by density-matrix renormalization group simulation of a concrete lattice model, where we find that topological superconductivity is robust against interactions. We further show that by stacking doped Haldane phases an exotic 2D anisotropic superconductor can be realized, {where the boundaries transverse to the {chain-}direction are either gapless or spontaneously symmetry-broken due to the Lieb-Schultz-Mattis (LSM) anomaly.}


title: An Iteratively Re-weighted Method for Problems with Sparsity-Inducing Norms

abstract: This work aims at solving the problems with intractable sparsity-inducing norms that are often encountered in various machine learning tasks, such as multi-task learning, subspace clustering, feature selection, robust principal component analysis, and so on. Specifically, an Iteratively Re-Weighted method (IRW) with solid convergence guarantee is provided. We investigate its convergence speed via numerous experiments on real data. Furthermore, in order to validate the practicality of IRW, we use it to solve a concrete robust feature selection model with complicated objective function. The experimental results show that the model coupled with proposed optimization method outperforms alternative methods significantly.


title: Lepton scattering from $^{40}$Ar and $^{48}$Ti in the quasielastic peak region

abstract: Neutron and proton spectral functions of $^{40}$Ar, $^{40}$Ca, and $^{48}$Ti isotopes are computed using the ab initio self-consistent Green's function approach. The resulting radii and charge distributions are in good agreement with available experimental data. The spectral functions of Ar and Ti are then utilized to calculate inclusive ($e$,$e$') cross sections within a factorization scheme and are found to correctly reproduce the recent Jefferson Lab measurements. Based on these successful agreements, the weak charged and neutral current double-differential cross sections for neutrino-$^{40}$Ar scattering are predicted in the quasielastic region. Results obtained by replacing the (experimentally inaccessible) neutron spectral distribution of $^{40}$Ar with the (experimentally accessible) proton distribution of $^{48}$Ti are compared and the accuracy of this approximation is assessed.


title: Growth Process of Hexagonal Boron Nitride in the Diffusion and Precipitation Method Studied by X-ray Photoelectron Spectroscopy

abstract: Submonolayer h-BN was grown on Ni foil in ultra-high vacuum by the diffusion and precipitation method and the growth process was studied by x-ray photoelectron spectroscopy. Formation of h-BN started to be observed at 600 C. All through the process, the surface was always slightly B-rich, which is consistent with the fact that B which is soluble in Ni at a high temperature can diffuse in Ni by the conventional bulk diffusion and insoluble N cannot. Moreover, both formation and decomposition of h-BN were found to occur at elevated temperatures possibly depending on provision of N atoms to the surface. On the Ni surface, decomposition of h-BN was observed at a relatively low temperature of 800 C.


title: Rotational evolution of the Vela pulsar during the 2016 glitch

abstract: The 2016 Vela glitch observed by the Mt Pleasant radio telescope provides the first opportunity to study pulse-to-pulse dynamics of a pulsar glitch, opening up new possibilities to study the neutron star's interior. We fit models of the star's rotation frequency to the pulsar data, and present three new results. First, we constrain the glitch rise time to less than 12.6s with 90% confidence, almost three times shorter than the previous best constraint. Second, we find definitive evidence for a rotational-frequency overshoot and fast relaxation following the glitch. Third, we find evidence for a slow-down of the star's rotation immediately prior to the glitch. The overshoot is predicted theoretically by some models; we discuss implications of the glitch rise and overshoot decay times on internal neutron-star physics. The slow down preceding the glitch is unexpected; we propose the slow-down may trigger the glitch by causing a critical lag between crustal superfluid and the crust.


title: Fitting functions on the cheap: the relative nonlinear matter power spectrum

abstract: We propose an alternative approach to the construction of fitting functions to the nonlinear matter power spectrum extracted from $N$-body simulations based on the relative matter power spectrum $\delta(k,a)$, defined as the fractional deviation in the absolute matter power spectrum produced by a target cosmology away from a reference $\Lambda$CDM prediction. From the computational perspective, $\delta(k,a)$ is fairly insensitive to the specifics of the simulation settings, and numerical convergence at the 1%-level can be readily achieved without the need for huge computing capacity. Furthermore, $\delta(k,a)$ exhibits several interesting properties that enable a piece-wise construction of the full fitting function, whereby component fitting functions are sought for single-parameter variations and then multiplied together to form the final product. Then, to obtain 1%-accurate absolute power spectrum predictions for any target cosmology only requires that the community as a whole invests in producing one single ultra-precise reference $\Lambda$CDM absolute power spectrum, to be combined with the fitting function to produce the desired result. To illustrate the power of this approach, we have constructed the fitting function RelFit using only five relatively inexpensive $w$CDM simulations (box length $L=256 h^{-1}$Mpc, $N=1024^3$ particles, initialised at $z_i=49$). In a 6-parameter space spanning $\{\omega_m,A_s,n_s,w,\omega_b,h\}$, the output relative power spectra of RelFit are consistent with the predictions of the CosmicEmu emulator to 1% or better for a wide range of cosmologies up to $k\simeq 10$/Mpc. Thus, our approach could provide an inexpensive and democratically accessible route to fulfilling the 1%-level accuracy demands of the upcoming generation of large-scale structure probes, especially in the exploration of "non-standard" or "exotic" cosmologies on nonlinear scales.


title: Nonlinear stability of explicit self-similar solutions for the timelike extremal hypersurfaces in R^{1+3}

abstract: This paper is devoted to the study of the singularity phenomenon of timelike extremal hypersurfaces in Minkowski spacetime $\mathbb{R}^{1+3}$. We find that there are two explicit lightlike self-similar solutions to a graph representation of timelike extremal hypersurfaces in Minkowski spacetime $\mathbb{R}^{1+3}$, the geometry of them are two spheres. The linear mode unstable of those lightlike self-similar solutions for the radially symmetric membranes equation is given. After that, we show those self-similar solutions of the radially symmetric membranes equation are nonlinearly stable inside a strictly proper subset of the backward lightcone. This means that the dynamical behavior of those two spheres is as attractors. Meanwhile, we overcome the double roots case (the theorem of Poincar\'{e} can't be used) in solving the difference equation by construction of a Newton's polygon when we carry out the analysis of spectrum for the linear operator.


title: Convergence Rates of Smooth Message Passing with Rounding in Entropy-Regularized MAP Inference

abstract: Maximum a posteriori (MAP) inference is a fundamental computational paradigm for statistical inference. In the setting of graphical models, MAP inference entails solving a combinatorial optimization problem to find the most likely configuration of the discrete-valued model. Linear programming (LP) relaxations in the Sherali-Adams hierarchy are widely used to attempt to solve this problem, and smooth message passing algorithms have been proposed to solve regularized versions of these LPs with great success. This paper leverages recent work in entropy-regularized LPs to analyze convergence rates of a class of edge-based smooth message passing algorithms to $\epsilon$-optimality in the relaxation. With an appropriately chosen regularization constant, we present a theoretical guarantee on the number of iterations sufficient to recover the true integral MAP solution when the LP is tight and the solution is unique.


title: A class of large solution of 2D tropical climate model without thermal diffusion

abstract: In this paper, we consider the Cauchy problem of 2D tropical climate model without thermal diffusion and construct global smooth solutions by choosing a class of special initial data whose $L^{\infty}$ norm can be arbitrarily large.


title: New Results for the Complexity of Resilience for Binary Conjunctive Queries with Self-Joins

abstract: The resilience of a Boolean query is the minimum number of tuples that need to be deleted from the input tables in order to make the query false. A solution to this problem immediately translates into a solution for the more widely known problem of deletion propagation with source-side effects. In this paper, we give several novel results on the hardness of the resilience problem for $\textit{binary conjunctive queries with self-joins}$ (i.e. conjunctive queries with relations of maximal arity 2) with one repeated relation. Unlike in the self-join free case, the concept of triad is not enough to fully characterize the complexity of resilience. We identify new structural properties, namely chains, confluences and permutations, which lead to various $NP$-hardness results. We also give novel involved reductions to network flow to show certain cases are in $P$. Overall, we give a dichotomy result for the restricted setting when one relation is repeated at most 2 times, and we cover many of the cases for 3. Although restricted, our results provide important insights into the problem of self-joins that we hope can help solve the general case of all conjunctive queries with self-joins in the future.


title: Nodal Line Spin-gapless Semimetals and High-quality Candidate Materials

abstract: Spin-gapless semimetals (SGSMs), which generate 100\% spin polarization, are viewed as promising semi-half-metals in spintronics with high speed and low consumption. We propose and characterize a new $\mathbb{Z_{\mathrm{2}}}$ class of topological nodal line (TNL) in SGSMs. The proposed TNLSGSMs are protected by space-time inversion symmetry or glide mirror symmetry with two-dimensional (2D) fully spin-polarized nearly flat surface states. Based on first-principles calculations and effective model analysis, a series of high-quality materials with $\textit{R}\overline{3}\textit{c}$ and $\textit{R}{3}\textit{c}$ space groups are predicted to realize such TNLSGSMs (chainlike). The 2D fully spin-polarized nearly flat surface states may provide a route to achieving equal spin pairing topological superconductivity as well as topological catalysts.


title: Learnable Gated Temporal Shift Module for Deep Video Inpainting

abstract: How to efficiently utilize temporal information to recover videos in a consistent way is the main issue for video inpainting problems. Conventional 2D CNNs have achieved good performance on image inpainting but often lead to temporally inconsistent results where frames will flicker when applied to videos (see https://www.youtube.com/watch?v=87Vh1HDBjD0&list=PLPoVtv-xp_dL5uckIzz1PKwNjg1yI0I94&index=1); 3D CNNs can capture temporal information but are computationally intensive and hard to train. In this paper, we present a novel component termed Learnable Gated Temporal Shift Module (LGTSM) for video inpainting models that could effectively tackle arbitrary video masks without additional parameters from 3D convolutions. LGTSM is designed to let 2D convolutions make use of neighboring frames more efficiently, which is crucial for video inpainting. Specifically, in each layer, LGTSM learns to shift some channels to its temporal neighbors so that 2D convolutions could be enhanced to handle temporal information. Meanwhile, a gated convolution is applied to the layer to identify the masked areas that are poisoning for conventional convolutions. On the FaceForensics and Free-form Video Inpainting (FVI) dataset, our model achieves state-of-the-art results with simply 33% of parameters and inference time.


title: Astraea: Self-balancing Federated Learning for Improving Classification Accuracy of Mobile Deep Learning Applications

abstract: Federated learning (FL) is a distributed deep learning method which enables multiple participants, such as mobile phones and IoT devices, to contribute a neural network model while their private training data remains in local devices. This distributed approach is promising in the edge computing system where have a large corpus of decentralized data and require high privacy. However, unlike the common training dataset, the data distribution of the edge computing system is imbalanced which will introduce biases in the model training and cause a decrease in accuracy of federated learning applications. In this paper, we demonstrate that the imbalanced distributed training data will cause accuracy degradation in FL. To counter this problem, we build a self-balancing federated learning framework call Astraea, which alleviates the imbalances by 1) Global data distribution based data augmentation, and 2) Mediator based multi-client rescheduling. The proposed framework relieves global imbalance by runtime data augmentation, and for averaging the local imbalance, it creates the mediator to reschedule the training of clients based on Kullback-Leibler divergence (KLD) of their data distribution. Compared with FedAvg, the state-of-the-art FL algorithm, Astraea shows +5.59% and +5.89% improvement of top-1 accuracy on the imbalanced EMNIST and imbalanced CINIC-10 datasets, respectively. Meanwhile, the communication traffic of Astraea can be 82% lower than that of FedAvg.


title: A Local Perspective on the Edge Removal Problem

abstract: The edge removal problem studies the loss in network coding rates that results when a network communication edge is removed from a given network. It is known, for example, that in networks restricted to linear coding schemes and networks restricted to Abelian group codes, removing an edge e* with capacity Re* reduces the achievable rate on each source by no more than Re*. In this work, we seek to uncover larger families of encoding functions for which the edge removal statement holds. We take a local perspective: instead of requiring that all network encoding functions satisfy certain restrictions (e.g., linearity), we limit only the function carried on the removed edge e*. Our central results give sufficient conditions on the function carried by edge e* in the code used to achieve a particular rate vector under which we can demonstrate the achievability of a related rate vector once e* is removed.


title: Comparing introductory physics and astronomy students' attitudes and approaches to problem solving

abstract: We examined how introductory physics students' attitudes and approaches to problem solving compare to those of introductory astronomy students, using a previously validated survey, the Attitudes and Approaches to Problem Solving (AAPS) survey. In addition, we compared the performance of the introductory physics and astronomy students on the factors which were identified in a factor analysis in the original validation study. We found that introductory astronomy students' overall average AAPS score was significantly more favorable than that of introductory physics students (p < 0.01), and the effect size was large (Cohen's d = 0.81). We also found that introductory astronomy students' scores were more favorable in all clusters of questions except for one factor involving drawing diagrams and writing scratchwork while solving problems. We also found that introductory physics and astronomy students were equally capable of solving two isomorphic problems posed to them, and that the majority of introductory physics and introductory astronomy students reported that the problem posed in the astronomy context was more interesting to them. Interviews suggest that the context of astronomy in problem solving may be more interesting for students and could be one possible explanation for the more favorable AAPS scores amongst introductory astronomy students compared to introductory physics students. Instructors of introductory physics courses should heed these findings which indicate that it may be beneficial for instructors of introductory physics courses to incorporate problems into their instruction which contain real-world contexts, which may serve to increase student interest-level, and which could help create more favorable attitudes and approaches towards problem solving.


title: On strong exceptional collections of line bundles of maximal length on Fano toric Deligne-Mumford stacks

abstract: We study strong exceptional collections of line bundles on Fano toric Deligne-Mumford stacks $\mathbb{P}_{\mathbf{\Sigma}}$ with rank of Picard group at most two. We prove that any strong exceptional collection of line bundles generates the derived category of $\mathbb{P}_{\mathbf{\Sigma}}$, as long as the number of elements in the collection equals the rank of the (Grothendieck) $K$-theory group of $\mathbb{P}_{\mathbf{\Sigma}}$.


title: Finding Outliers in Gaussian Model-Based Clustering

abstract: Clustering, or unsupervised classification, is a task often plagued by outliers. Yet there is a paucity of work on handling outliers in clustering. Outlier identification algorithms tend to fall into three broad categories: outlier inclusion, outlier trimming, and post hoc outlier identification methods, with the former two often requiring pre-specification of the number of outliers. The fact that sample squared Mahalanobis distance is beta-distributed is used to derive an approximate distribution for the log-likelihoods of subset finite Gaussian mixture models. An algorithm is then proposed that removes the least plausible points according to the subset log-likelihoods, which are deemed outliers, until the subset log-likelihoods adhere to the reference distribution. This results in a trimming method, called OCLUST, that inherently estimates the number of outliers.


title: Asteroseismology of the DOV star PG 1159-035

abstract: Grids of DOV star models are evolved by \texttt{WDEC} with fixed atmospheric constituent to the spectral values of $X_{C}/X_{He}/X_{O}$ = 50/33/17. The core compositions are from white dwarf models at highest $T_{eff}$ evolved by \texttt{MESA}. The eigenfrequencies are calculated and used to fit the observed modes. Based on 264.1 hours of photometric observations on PG 1159-035, Winget et al. identified 125 individual frequencies. Costa et al. identified 198 pulsation modes for PG 1159-035 according to the WET photometric data from 1983, 1985, 1989, and 2002. Both of them derived frequency splitting values of $\delta\sigma_{l=1}$ $\sim$ 4.2\,$\mu$Hz and $\delta\sigma_{l=2}$ $\sim$ 6.9\,$\mu$Hz. According to the values of $\delta\sigma_{l=1}$ and $\delta\sigma_{l=2}$, 20 triplets and 9 quintuplets are selected and used to constrain the fitting models. Our optimal model has $T_{eff}$ = 129000\,K, $M_{*}$ = 0.63\,$M_{\odot}$, log$g$ = 7.59, log($M_{\rm env}/M_{*}$) = -5.0, and $\sigma_{RMS}$ = 1.97\,s. The values of $T_{eff}$ and log$g$ are consistent with that values of C\'orsico et al.. The calculated modes of minimum rate of period change correspond to modes with maximum kinetic energy distributed in the envelope. The observed rates of period change with positive and negative values can also be partially reproduced. In particular, there are negative rates of period change for the calculated modes from our optimal model, which is not found in previous work.


title: Multimessenger Asteroseismology of Core-Collapse Supernovae

abstract: We investigate correlated gravitational wave and neutrino signals from rotating core-collapse supernovae with simulations. Using an improved mode identification procedure based on mode function matching, we show that a linear quadrupolar mode of the core produces a dual imprint on gravitational waves and neutrinos in the early post-bounce phase of the supernova. The angular harmonics of the neutrino emission are consistent with the mode energy around the neutrinospheres, which points to a mechanism for the imprint on neutrinos. Thus, neutrinos carry information about the mode amplitude in the outer region of the core, whereas gravitational waves probe deeper in. We also find that the best-fit mode function has a frequency bounded above by $\sim 420$ Hz, and yet the mode's frequency in our simulations is $\sim 15\%$ higher, due to the use of Newtonian hydrodynamics and a widely used pseudo-Newtonian gravity approximation. This overestimation is particularly important for the analysis of gravitational wave detectability and asteroseismology, pointing to limitations of pseudo-Newtonian approaches for these purposes, possibly even resulting in excitation of incorrect modes. In addition, mode frequency matching (as opposed to mode function matching) could be resulting in mode misidentification in recent work. Lastly, we evaluate the prospects of a multimessenger detection of the mode using current technology. The detection of the imprint on neutrinos is most challenging, with a maximum detection distance of $\sim\!1$ kpc using the IceCube Neutrino Observatory. The maximum distance for detecting the complementary gravitational wave imprint is $\sim\!5$ kpc using Advanced LIGO at design sensitivity.


title: A posteriori error analysis for Schwarz overlapping domain decomposition methods

abstract: Domain decomposition methods are widely used for the numerical solution of partial differential equations on high performance computers. We develop an adjoint-based a posteriori error analysis for both multiplicative and additive overlapping Schwarz domain decomposition methods. The numerical error in a user-specified functional of the solution (quantity of interest) is decomposed into contributions that arise as a result of the finite iteration between the subdomains and from the spatial discretization. The spatial discretization contribution is further decomposed into contributions arising from each subdomain. This decomposition of the numerical error is used to construct a two stage solution strategy that efficiently reduces the error in the quantity of interest by adjusting the relative contributions to the error.


title: Three-terminal interface as a thermoelectric generator beyond Seebeck effect

abstract: We investigate thermoelectric transport through interfaces with inelastic scatterings by developing a quantum theory, which has been extensively validated by existing theories. We find that under temperature bias, while a two-terminal conductor-insulator interface behaves only as a thermal resistor, a three-terminal conductor-insulator-conductor interface can function as an electricity generator caused by phonon-mediated electron scatterings with heat-charge current separation. Unlike conventional thermoelectrics which is a property of a bulk caused by the Seebeck effect, this thermoelectric behavior is a property of an interface driven by electron-phonon scatterings.


title: High-speed Railway Fastener Detection and Localization Method based on convolutional neural network

abstract: Railway transportation is the artery of China's national economy and plays an important role in the development of today's society. Due to the late start of China's railway security inspection technology, the current railway security inspection tasks mainly rely on manual inspection, but the manual inspection efficiency is low, and a lot of manpower and material resources are needed. In this paper, we establish a steel rail fastener detection image dataset, which contains 4,000 rail fastener pictures about 4 types. We use the regional suggestion network to generate the region of interest, extracts the features using the convolutional neural network, and fuses the classifier into the detection network. With online hard sample mining to improve the accuracy of the model, we optimize the Faster RCNN detection framework by reducing the number of regions of interest. Finally, the model accuracy reaches 99% and the speed reaches 35FPS in the deployment environment of TITAN X GPU.


title: Fast Algorithms for Surface Reconstruction from Point Cloud

abstract: We consider constructing a surface from a given set of point cloud data. We explore two fast algorithms to minimize the weighted minimum surface energy in [Zhao, Osher, Merriman and Kang, Comp.Vision and Image Under., 80(3):295-319, 2000]. An approach using Semi-Implicit Method (SIM) improves the computational efficiency through relaxation on the time-step constraint. An approach based on Augmented Lagrangian Method (ALM) reduces the run-time via an Alternating Direction Method of Multipliers-type algorithm, where each sub-problem is solved efficiently. We analyze the effects of the parameters on the level-set evolution and explore the connection between these two approaches. We present numerical examples to validate our algorithms in terms of their accuracy and efficiency.


title: Bootstrapping Inflationary Correlators in Mellin Space

abstract: We develop a Mellin space approach to boundary correlation functions in anti-de Sitter (AdS) and de Sitter (dS) spaces. Using the Mellin-Barnes representation of correlators in Fourier space, we show that the analytic continuation between AdS$_{d+1}$ and dS$_{d+1}$ is encoded in a collection of simple relative phases. This allows us to determine the late-time tree-level three-point correlators of spinning fields in dS$_{d+1}$ from known results for Witten diagrams in AdS$_{d+1}$ by multiplication with a simple trigonometric factor. At four point level, we show that Conformal symmetry fixes exchange four-point functions both in AdS$_{d+1}$ and dS$_{d+1}$ in terms of the dual Conformal Partial Wave (which in Fourier space is a product of boundary three-point correlators) up to a factor which is determined by the boundary conditions. In this work we focus on late-time four-point correlators with external scalars and an exchanged field of integer spin-$\ell$. The Mellin-Barnes representation makes manifest the analytic structure of boundary correlation functions, providing an analytic expression for the exchange four-point function which is valid for general $d$ and generic scaling dimensions, in particular massive, light and (partially-)massless fields. When $d=3$ we reproduce existing explicit results available in the literature for external conformally coupled and massless scalars. From these results, assuming the weak breaking of the de Sitter isometries, we extract the corresponding correction to the inflationary three-point function of general external scalars induced by a general spin-$\ell$ field at leading order in slow roll. These results provide a step towards a more systematic understanding of de Sitter observables at tree level and beyond using Mellin space methods.


title: Disentangled Makeup Transfer with Generative Adversarial Network

abstract: Facial makeup transfer is a widely-used technology that aims to transfer the makeup style from a reference face image to a non-makeup face. Existing literature leverage the adversarial loss so that the generated faces are of high quality and realistic as real ones, but are only able to produce fixed outputs. Inspired by recent advances in disentangled representation, in this paper we propose DMT (Disentangled Makeup Transfer), a unified generative adversarial network to achieve different scenarios of makeup transfer. Our model contains an identity encoder as well as a makeup encoder to disentangle the personal identity and the makeup style for arbitrary face images. Based on the outputs of the two encoders, a decoder is employed to reconstruct the original faces. We also apply a discriminator to distinguish real faces from fake ones. As a result, our model can not only transfer the makeup styles from one or more reference face images to a non-makeup face with controllable strength, but also produce various outputs with styles sampled from a prior distribution. Extensive experiments demonstrate that our model is superior to existing literature by generating high-quality results for different scenarios of makeup transfer.


title: The generalized orthogonal Procrustes problem in the high noise regime

abstract: We consider the problem of estimating a cloud of points from numerous noisy observations of that cloud after unknown rotations, and possibly reflections. This is an instance of the general problem of estimation under group action, originally inspired by applications in 3-D imaging and computer vision. We focus on a regime where the noise level is larger than the magnitude of the signal, so much so that the rotations cannot be estimated reliably. We propose a simple and efficient procedure based on invariant polynomials (effectively: the Gram matrices) to recover the signal, and we assess it against fundamental limits of the problem that we derive. We show our approach adapts to the noise level and is statistically optimal (up to constants) for both the low and high noise regimes. In studying the variance of our estimator, we encounter the question of the sensivity of a type of thin Cholesky factorization, for which we provide an improved bound which may be of independent interest.


title: Generalized hyperbolicity for linear operators

abstract: It is introduced an open class of linear operators on Banach and Hilbert spaces such that their non-wandering set is an infinite dimensional topologically mixing subspace. In certain cases, the non-wandering set coincides with the whole space.


title: Localization of Fr\'echet frames and expansion of generalized functions

abstract: Matrix type operators with the off-diagonal decay of polynomial or sub-exponential types are revisited with weaker assumptions concerning row or column estimates, still giving the continuity results for the frame type operators. Such results are extended from Banach to Fr\'{e}chet spaces. Moreover, the localization of Fr\'{e}chet frames is used for the frame expansions of tempered distributions and a class of Beurling ultradistributions.


title: The Impact of Time Delay in a Tumor Model

abstract: In this paper we consider a free boundary tumor growth model with a time delay in cell proliferation and study how time delay affects the stability and the size of the tumor. The model is a coupled system of an elliptic equation, a parabolic equation and an ordinary differential equation. It incorporates the cell location under the presence of time delay, with the tumor boundary as a free boundary. A parameter $\mu$ in the model is proportional to the "aggressiveness" of the tumor. It is proved that there exists a unique classical radially symmetric stationary solution $(\sigma_*, p_*, R_*)$ which is stable for any $\mu > 0$ with respect to all radially symmetric perturbations (c.f. \cite{delay1}). However, under non-radially symmetric perturbations, we prove that there exists a critical value $\mu_*$, such that if $\mu<\mu_*$ then the stationary solution $(\sigma_*, p_*, R_*)$ is linearly stable; whereas if $\mu>\mu_*$ the stationary solution is unstable. It is actually unrealistic to expect the problem to be stable for large tumor aggressiveness parameter, therefore our result is more reasonable. Furthermore, we established that adding the time delay in the model would result in a larger stationary tumor, and if the tumor aggressiveness parameter is larger, then the time delay would have a greater impact on the size of the tumor.


title: Hyperspectral Super-Resolution via Global-Local Low-Rank Matrix Estimation

abstract: Hyperspectral super-resolution (HSR) is a problem that aims to estimate an image of high spectral and spatial resolutions from a pair of co-registered multispectral (MS) and hyperspectral (HS) images, which have coarser spectral and spatial resolutions, respectively. In this paper we pursue a low-rank matrix estimation approach for HSR. We assume that the spectral-spatial matrices associated with the whole image and the local areas of the image have low-rank structures. The local low-rank assumption, in particular, has the aim of providing a more flexible model for accounting for local variation effects due to endmember variability. We formulate the HSR problem as a global-local rank-regularized least-squares problem. By leveraging on the recent advances in non-convex large-scale optimization, namely, the smooth Schatten-p approximation and the accelerated majorization-minimization method, we develop an efficient algorithm for the global-local low-rank problem. Numerical experiments on synthetic, semi-real and real data show that the proposed algorithm outperforms a number of benchmark algorithms in terms of recovery performance.


title: Multi-scale Template Matching with Scalable Diversity Similarity in an Unconstrained Environment

abstract: We propose a novel multi-scale template matching method which is robust against both scaling and rotation in unconstrained environments. The key component behind is a similarity measure referred to as scalable diversity similarity (SDS). Specifically, SDS exploits bidirectional diversity of the nearest neighbor (NN) matches between two sets of points. To address the scale-robustness of the similarity measure, local appearance and rank information are jointly used for the NN search. Furthermore, by introducing penalty term on the scale change, and polar radius term into the similarity measure, SDS is shown to be a well-performing similarity measure against overall size and rotation changes, as well as non-rigid geometric deformations, background clutter, and occlusions. The properties of SDS are statistically justified, and experiments on both synthetic and real-world data show that SDS can significantly outperform state-of-the-art methods.


title: Polarization difference between hyperons and anti-hyperons induced by external magnetic field

abstract: We investigate the quantum correlated $\Lambda \bar{\Lambda} $ production in the reaction $e^{+}e^{-} \to J/\psi \to \Lambda \bar{\Lambda}$. Since the $\Lambda$ or $\bar{\Lambda}$ has a nonzero magnetic moment, its spin will undergo a Larmor precession in the magnetic field of the detector, such as the BESIII experiment. Because of the spin precession, the angular distribution of the $\Lambda$ and $\bar{\Lambda}$ is slightly modified. Therefore, we obtain the corresponding term of the modified angular distribution due to the effect of the Larmor precession. We also estimate its potential effect on the measurements of $CP$ violation, as well as the decay asymmetry parameter and polarization of $\Lambda$. The polarization of the $\Lambda$ or $\bar{\Lambda}$ at the production vertex will rotate around the $B$-field axis, over an angle depending on the flight length in, but it still could be measured by fit to the corrected angular distribution. Of important, We conclude that a nonzero $CP$ asymmetry of order $10^{-4}$ will be caused once neglecting spin precession of the $\Lambda$ and $\bar{\Lambda}$ in the $e^{+}e^{-} \to J/\psi \to \Lambda \bar{\Lambda}$ process. The size of this $CP$ asymmetry is several times that of predicted within Standard Model in the hyperon decay. Although this effect is small, it will play an important role in future high precision experiments, such as the super-tau-charm factory.


title: Covariant momentum map for non-Abelian topological BF field theory

abstract: We analyze the inherent symmetries associated to the non-Abelian topological BF theory from the geometric and covariant perspectives of the Lagrangian and the multisymplectic formalisms. At the Lagrangian level, we classify the symmetries of the theory as natural and Noether symmetries and construct the associated Noether currents, while at the multisymplectic level the symmetries of the theory arise as covariant canonical transformations. These transformations allowed us to build within the multisymplectic approach, in a complete covariant way, the momentum maps which are analogous to the conserved Noether currents. The covariant momentum maps are fundamental to recover, after the space plus time decomposition of the background manifold, not only the extended Hamiltonian of the BF theory but also the generators of the gauge transformations which arise in the instantaneous Dirac-Hamiltonian analysis of the first-class constraint structure that characterizes the BF model under study. To the best of our knowledge, this is the first non-trivial physical model associated to General Relativity for which both natural and Noether symmetries have been analyzed at the multisymplectic level. Our study shed some light on the understanding of the manner in which the generators of gauge transformations may be recovered from the multisymplectic formalism for field theory.


title: Second-order topological phases protected by chiral symmetry

abstract: We study second-order topological insulators and semimetals characterized by chiral symmetry. We investigate topological phase transitions of a model for construction of the two-dimensional second-order topological insulators protected only by chiral symmetry. By the theory of the phase transitions, we propose a second-order topological semimetal and insulators with flat hinge bands in chiral-symmetric three-dimensional systems. The three-dimensional second-order topological phases can be obtained from the stacked two-dimensional second-order topological insulators with chiral symmetry. Moreover, we show that broken chiral symmetry in the three-dimensional second-order topological phase allows a second-order topological insulator with chiral hinge states. We also demonstrate the second-order topological phases by using a lattice model.


title: Adaptive Music Composition for Games

abstract: The generation of music that adapts dynamically to content and actions has an important role in building more immersive, memorable and emotive game experiences. To date, the development of adaptive music systems for video games is limited by both the nature of algorithms used for real-time music generation and the limited modelling of player action, game world context and emotion in current games. We propose that these issues must be addressed in tandem for the quality and flexibility of adaptive game music to significantly improve. Cognitive models of knowledge organisation and emotional affect are integrated with multi-modal, multi-agent composition techniques to produce a novel Adaptive Music System (AMS). The system is integrated into two stylistically distinct games. Gamers reported an overall higher immersion and correlation of music with game-world concepts with the AMS than with the original game soundtracks in both games.


title: Isometric copies of $\ell_\infty^n$ and $\ell_1^n$ in transportation cost spaces on finite metric spaces

abstract: Main results: (a) If a metric space contains $2n$ elements, the transportation cost space on it contains a $1$-complemented isometric copy of $\ell_1^n$. (b) An example of a finite metric space whose transportation cost space contains an isometric copy of $\ell_\infty^4$. Transportation cost spaces are also known as Arens-Eells, Lipschitz-free, or Wasserstein $1$ spaces.


title: Distinguishing open symplectic mapping tori via their wrapped Fukaya categories

abstract: In this paper, we present partial results towards a classification of symplectic mapping tori using dynamical properties of wrapped Fukaya categories. More precisely, we construct a symplectic manifold $T_\phi$ associated to a Weinstein domain $M$, and an exact, compactly supported symplectomorphism $\phi$. $T_\phi$ is another Weinstein domain and its contact boundary is independent of $\phi$. In this paper, we distinguish $\phi$ from $T_{1_M}$, under certain assumptions (Theorem 1.1). As an application, we obtain pairs of diffeomorphic Weinstein domains with the same contact boundary and whose symplectic cohomology groups are the same, as vector spaces, but that are different as Liouville domains. To our knowledge, this is the first example of such pairs that can be distinguished by their wrapped Fukaya category. Previously, we have suggested a categorical model $M_\phi$ for the wrapped Fukaya category $\mathcal{W}(T_\phi)$, and we have distinguished $M_\phi$ from the mapping torus category of the identity. In this paper, we prove $\mathcal{W}(T_\phi)$ and $M_\phi$ are derived equivalent (Theorem 1.9); hence, deducing the promised Theorem 1.1. Theorem 1.9 is of independent interest as it preludes an algebraic description of wrapped Fukaya categories of locally trivial symplectic fibrations as twisted tensor products.


title: Some $q$-exponential formulas involving the double lowering operator $\psi$ for a tridiagonal pair

abstract: Let $\mathbb{K}$ denote an algebraically closed field and let $V$ denote a vector space over $\mathbb{K}$ with finite positive dimension. Let $A,A^*$ denote a tridiagonal pair on $V$. We assume that $A,A^*$ belongs to a family of tridiagonal pairs said to have $q$-Racah type. Let $\{U_i\}_{i=0}^d$ and $\{U_i^\Downarrow\}_{i=0}^{d}$ denote the first and second split decompositions of $V$. In an earlier paper we introduced a double lowering operator $\psi:V\to V$ with the notable feature that both $\psi U_i\subseteq U_{i-1}$ and $\psi U_i^\Downarrow\subseteq U_{i-1}^\Downarrow$ for $0\leq i\leq d$, where $U_{-1}=0$ and $U_{-1}^\Downarrow=0$. In the same paper, we showed that there exists a unique linear transformation $\Delta:V\to V$ such that $\Delta(U_i)\subseteq U_i^{\Downarrow}$ and $(\Delta -I)U_i\subseteq U_0+U_1+\cdots +U_{i-1}$ for $0\leq i \leq d$. In the present paper, we show that $\Delta$ can be expressed as a product of two linear transformations; one is a $q$-exponential in $\psi$ and the other is a $q^{-1}$-exponential in $\psi$. We view $\Delta$ as a transition matrix from the first split decomposition of $V$ to the second. Consequently, we view the $q^{-1}$-exponential in $\psi$ as a transition matrix from the first split decomposition to a decomposition of $V$ which we interpret as a kind of halfway point. This halfway point turns out to be the eigenspace decomposition of a certain linear transformation $\mathcal{M}$. We discuss the eigenspace decomposition of $\mathcal{M}$ and give the actions of various operators on this decomposition.


title: $^{16}O(p,\alpha)^{13}N$ makes explosive oxygen burning sensitive to the metallicity of the progenitors of type Ia supernovae

abstract: Even though the main nucleosynthetic products of type Ia supernovae belong to the iron-group, intermediate-mass alpha-nuclei (silicon, sulfur, argon, and calcium) stand out in their spectra up to several weeks past maximum brightness. Recent measurements of the abundances of calcium, argon, and sulfur in type Ia supernova remnants have been interpreted in terms of metallicity-dependent oxygen burning, in accordance with previous theoretical predictions. It is known that $\alpha$-rich oxygen burning results from $^{16}$O$\rightarrow^{12}$C followed by efficient $^{12}$C+$^{12}$C fusion reaction, as compared to oxygen consumption by $^{16}$O fusion reactions, but the precise mechanism of dependence on the progenitor metallicity has remained unidentified so far. I show that the chain $^{16}$O(p,$\alpha$)$^{13}$N($\gamma$,p)$^{12}$C boosts $\alpha$-rich oxygen burning when the proton abundance is large, increasing the synthesis of argon and calcium with respect to sulfur and silicon. For high-metallicity progenitors, the presence of free neutrons leads to a drop in the proton abundance and the above chain is not efficient. Although the rate of $^{16}$O(p,$\alpha$)$^{13}$N can be found in astrophysical reaction rate libraries, its uncertainty is unconstrained. Assuming that all reaction rates other than $^{16}$O(p,$\alpha$)$^{13}$N retain their standard values, an increase by a factor of approximately seven of the $^{16}$O(p,$\alpha$)$^{13}$N rate at temperatures in the order $3-5\times10^9$ K is enough to explain the whole range of calcium-to-sulfur mass ratios measured in Milky Way and LMC supernova remnants. These same measurements provide a lower limit to the $^{16}$O(p,$\alpha$)$^{13}$N rate in the mentioned temperature range, on the order of a factor of 0.5 with respect to the rate reported in widely used literature tabulations.


title: Bundled Causal History Interaction

abstract: Complex system arises as a result of the nonlinear interactions between components. In particular, the evolutionary dynamics of a multivariate system encodes the ways in which different variables interact with each other individually or in groups. One fundamental question that remains unanswered is: how do two non-overlapping multivariate subsets of variables interact to causally determine the outcome of a specific variable? Here we provide an information based approach to address this problem. We delineate the temporal interactions between the bundles in a probabilistic graphical model. The strength of the interactions, captured by partial information decomposition, then exposes complex behavior of dependencies and memory within the system. The proposed approach successfully illustrates complex dependence between cations and anions as determinants of \textit{pH} in an observed stream chemistry system. This example demonstrates the potentially broad applicability of the approach, establishing the foundation to study the interaction between groups of variables in a range of complex systems.


title: WHAM!: Extending Speech Separation to Noisy Environments

abstract: Recent progress in separating the speech signals from multiple overlapping speakers using a single audio channel has brought us closer to solving the cocktail party problem. However, most studies in this area use a constrained problem setup, comparing performance when speakers overlap almost completely, at artificially low sampling rates, and with no external background noise. In this paper, we strive to move the field towards more realistic and challenging scenarios. To that end, we created the WSJ0 Hipster Ambient Mixtures (WHAM!) dataset, consisting of two speaker mixtures from the wsj0-2mix dataset combined with real ambient noise samples. The samples were collected in coffee shops, restaurants, and bars in the San Francisco Bay Area, and are made publicly available. We benchmark various speech separation architectures and objective functions to evaluate their robustness to noise. While separation performance decreases as a result of noise, we still observe substantial gains relative to the noisy signals for most approaches.


title: Heteroclinic Orbits and Nonintegrability in Two-Degree-of-Freedom Hamiltonian Systems with Saddle-Centers

abstract: We consider a class of two-degree-of-freedom Hamiltonian systems with saddle-centers connected by heteroclinic orbits and discuss some relationships between the existence of transverse heteroclinic orbits and nonintegrability. By the Lyapunov center theorem there is a family of periodic orbits near each of the saddle-centers, and the Hessian matrices of the Hamiltonian at the two saddle-centers are assumed to have the same number of positive eigenvalues. We show that if the associated Jacobian matrices have the same pair of purely imaginary eigenvalues, then the stable and unstable manifolds of the periodic orbits intersect transversely on the same Hamiltonian energy surface when sufficient conditions obtained in previous work for real-meromorphic nonintegrability of the Hamiltonian systems hold; if not, then these manifolds intersect transversely on the same energy surface, have quadratic tangencies or do not intersect whether the sufficient conditions hold or not. Our theory is illustrated for a system with quartic single-well potential and some numerical results are given to support the theoretical results.


title: Sample Adaptive Multiple Kernel Learning for Failure Prediction of Railway Points

abstract: Railway points are among the key components of railway infrastructure. As a part of signal equipment, points control the routes of trains at railway junctions, having a significant impact on the reliability, capacity, and punctuality of rail transport. Traditionally, maintenance of points is based on a fixed time interval or raised after the equipment failures. Instead, it would be of great value if we could forecast points' failures and take action beforehand, minimising any negative effect. To date, most of the existing prediction methods are either lab-based or relying on specially installed sensors which makes them infeasible for large-scale implementation. Besides, they often use data from only one source. We, therefore, explore a new way that integrates multi-source data which are ready to hand to fulfil this task. We conducted our case study based on Sydney Trains rail network which is an extensive network of passenger and freight railways. Unfortunately, the real-world data are usually incomplete due to various reasons, e.g., faults in the database, operational errors or transmission faults. Besides, railway points differ in their locations, types and some other properties, which means it is hard to use a unified model to predict their failures. Aiming at this challenging task, we firstly constructed a dataset from multiple sources and selected key features with the help of domain experts. In this paper, we formulate our prediction task as a multiple kernel learning problem with missing kernels. We present a robust multiple kernel learning algorithm for predicting points failures. Our model takes into account the missing pattern of data as well as the inherent variance on different sets of railway points. Extensive experiments demonstrate the superiority of our algorithm compared with other state-of-the-art methods.


title: Analysis on hadron spectra in heavy-ion collisions with a new non-extensive approach

abstract: The transverse momentum spectra of identified charged hadrons stemming from high energy collisions at different beam energies are described by a new non-extensive distribution, the Kaniadakis $\kappa$-distribution, with respect to the constraints in non-extensive quantum statistics. All fittings are also compared with the Tsallis distributions as well as the usual Boltzmann-Gibbs one. $\chi^2/ndf$ is also used to test the fitting goodness of all functions. Our results show that these different non-extensive approaches can be well applied in high energy collisions rather than the classical one. The Kaniadakis statistics is typically better applied into such systems with both positive and negative particles considered. This provides an alternative non-extensive view to study high energy physics. Analysis on the fitting parameters are present as well. The similar relationships of all functions remind us of the further understanding of the non-extensivity.


title: Learning to Traverse Latent Spaces for Musical Score Inpainting

abstract: Music Inpainting is the task of filling in missing or lost information in a piece of music. We investigate this task from an interactive music creation perspective. To this end, a novel deep learning-based approach for musical score inpainting is proposed. The designed model takes both past and future musical context into account and is capable of suggesting ways to connect them in a musically meaningful manner. To achieve this, we leverage the representational power of the latent space of a Variational Auto-Encoder and train a Recurrent Neural Network which learns to traverse this latent space conditioned on the past and future musical contexts. Consequently, the designed model is capable of generating several measures of music to connect two musical excerpts. The capabilities and performance of the model are showcased by comparison with competitive baselines using several objective and subjective evaluation methods. The results show that the model generates meaningful inpaintings and can be used in interactive music creation applications. Overall, the method demonstrates the merit of learning complex trajectories in the latent spaces of deep generative models.


title: An iterative scheme for the generalized Peierls-Nabarro model based on the inverse Hilbert transform

abstract: A new semi-analytical iterative scheme is proposed in this work for solving the generalized Peierls-Nabarro model. The numerical method developed here exploits certain basic properties of the Hilbert transform to achieve the desired reduction of the non-local and non-linear equations characterizing the generalized Peierls-Nabarro model to a local fixed point iteration scheme. The method is validated with simple examples involving the 1D Peierls-Nabarro model corresponding to a sinusoidal stacking fault energy, and with calculations of the core structure of both edge and screw dislocations on the close-packed $\{111\}$ planes in Aluminium. An approximate technique to incorporate external stresses within the framework of the proposed iterative scheme is also discussed with applications to the equilibration of a dislocation dipole. Finally, the advantages, limitations and avenues for future extension of the proposed method are discussed.


title: Multimodal Transformer Networks for End-to-End Video-Grounded Dialogue Systems

abstract: Developing Video-Grounded Dialogue Systems (VGDS), where a dialogue is conducted based on visual and audio aspects of a given video, is significantly more challenging than traditional image or text-grounded dialogue systems because (1) feature space of videos span across multiple picture frames, making it difficult to obtain semantic information; and (2) a dialogue agent must perceive and process information from different modalities (audio, video, caption, etc.) to obtain a comprehensive understanding. Most existing work is based on RNNs and sequence-to-sequence architectures, which are not very effective for capturing complex long-term dependencies (like in videos). To overcome this, we propose Multimodal Transformer Networks (MTN) to encode videos and incorporate information from different modalities. We also propose query-aware attention through an auto-encoder to extract query-aware features from non-text modalities. We develop a training procedure to simulate token-level decoding to improve the quality of generated responses during inference. We get state of the art performance on Dialogue System Technology Challenge 7 (DSTC7). Our model also generalizes to another multimodal visual-grounded dialogue task, and obtains promising performance. We implemented our models using PyTorch and the code is released at https://github.com/henryhungle/MTN.


title: A Tandem Learning Rule for Effective Training and Rapid Inference of Deep Spiking Neural Networks

abstract: Spiking neural networks (SNNs) represent the most prominent biologically inspired computing model for neuromorphic computing (NC) architectures. However, due to the non-differentiable nature of spiking neuronal functions, the standard error back-propagation algorithm is not directly applicable to SNNs. In this work, we propose a tandem learning framework, that consists of an SNN and an Artificial Neural Network (ANN) coupled through weight sharing. The ANN is an auxiliary structure that facilitates the error back-propagation for the training of the SNN at the spike-train level. To this end, we consider the spike count as the discrete neural representation in the SNN, and design ANN neuronal activation function that can effectively approximate the spike count of the coupled SNN. The proposed tandem learning rule demonstrates competitive pattern recognition and regression capabilities on both the conventional frame-based and event-based vision datasets, with at least an order of magnitude reduced inference time and total synaptic operations over other state-of-the-art SNN implementations. Therefore, the proposed tandem learning rule offers a novel solution to training efficient, low latency, and high accuracy deep SNNs with low computing resources.


title: Kleene Theorems for Free Choice Nets Labelled with Distributed Alphabets

abstract: We provided (PNSE'2014) expressions for free choice nets having "distributed choice property" which makes the nets "direct product" representable. In a recent work (PNSE'2016), we gave equivalent syntax for a larger class of free choice nets obtained by dropping distributed choice property. In both these works, the classes of free choice nets were restricted by a "product condition" on the set of final markings. In this paper we do away with this restriction and give expressions for the resultant classes of nets which correspond to "free choice synchronous products and Zielonka automata". For free choice nets with distributed choice property, we give an alternative characterization using properties checkable in polynomial time. Free choice nets we consider are 1-bounded, S-coverable, and are labelled with distributed alphabets, where S-components of the associated S-cover respect the given alphabet distribution.


title: Can a Robot Hear the Shape and Dimensions of a Room?

abstract: Knowing the geometry of a space is desirable for many applications, e.g. sound source localization, sound field reproduction or auralization. In circumstances where only acoustic signals can be obtained, estimating the geometry of a room is a challenging proposition. Existing methods have been proposed to reconstruct a room from the room impulse responses (RIRs). However, the sound source and microphones must be deployed in a feasible region of the room for it to work, which is impractical when the room is unknown. This work propose to employ a robot equipped with a sound source and four acoustic sensors, to follow a proposed path planning strategy to moves around the room to collect first image sources for room geometry estimation. The strategy can effectively drives the robot from a random initial location through the room so that the room geometry is guaranteed to be revealed. Effectiveness of the proposed approach is extensively validated in a synthetic environment, where the results obtained are highly promising.


title: Bayesian Analysis of High-dimensional Discrete Graphical Models

abstract: This work introduces a Bayesian methodology for fitting large discrete graphical models with spike-and-slab priors to encode sparsity. We consider a quasi-likelihood approach that enables node-wise parallel computation resulting in reduced computational complexity. We introduce a scalable Langevin MCMC algorithm for sampling from the quasi-posterior distribution which enables variable selection and estimation simultaneously. We present extensive simulation results to demonstrate scalability and accuracy of the method. We also analyze the 16 Personality Factors (PF) dataset to illustrate performance of the method.


title: Perforating freestanding molybdenum disulfide monolayers with highly charged ions

abstract: Porous single layer molybdenum disulfide (MoS$_2$) is a promising material for applications such as DNA sequencing and water desalination. In this work, we introduce irradiation with highly charged ions (HCIs) as a new technique to fabricate well-defined pores in MoS$_2$. Surprisingly, we find a linear increase of the pore creation efficiency over a broad range of potential energies. Comparison to atomistic simulations reveals the critical role of energy deposition from the ion to the material through electronic excitation in the defect creation process, and suggests an enrichment in molybdenum in the vicinity of the pore edges at least for ions with low potential energies. Analysis of the irradiated samples with atomic resolution scanning transmission electron microscopy reveals a clear dependence of the pore size on the potential energy of the projectiles, establishing irradiation with highly charged ions as an effective method to create pores with narrow size distributions and radii between ca. 0.3 and 3 nm.


title: Procedure Planning in Instructional Videos

abstract: In this paper, we study the problem of procedure planning in instructional videos, which can be seen as a step towards enabling autonomous agents to plan for complex tasks in everyday settings such as cooking. Given the current visual observation of the world and a visual goal, we ask the question "What actions need to be taken in order to achieve the goal?". The key technical challenge is to learn structured and plannable state and action spaces directly from unstructured videos. We address this challenge by proposing Dual Dynamics Networks (DDN), a framework that explicitly leverages the structured priors imposed by the conjugate relationships between states and actions in a learned plannable latent space. We evaluate our method on real-world instructional videos. Our experiments show that DDN learns plannable representations that lead to better planning performance compared to existing planning approaches and neural network policies.


title: Numerical modelling of shock-bubble interactions using a pressure-based algorithm without Riemann solvers

abstract: The interaction of a shock wave with a bubble features in many engineering and emerging technological applications, and has been used widely to test new numerical methods for compressible interfacial flows. Recently, density-based algorithms with pressure-correction methods as well as fully-coupled pressure-based algorithms have been established as promising alternatives to classical density-based algorithms based on Riemann solvers. The current paper investigates the predictive accuracy of fully-coupled pressure-based algorithms without Riemann solvers in modelling the interaction of shock waves with one-dimensional and two-dimensional bubbles in gas-gas and liquid-gas flows. For a gas bubble suspended in another gas, the mesh resolution and the applied advection schemes are found to only have a minor influence on the bubble shape and position, as well as the behaviour of the dominant shock waves and rarefaction fans. For a gas bubble suspended in a liquid, however, the mesh resolution has a critical influence on the shape, the position and the post-shock evolution of the bubble, as well as the pressure and temperature distribution.


title: Multi-TeV flaring in nearby High Energy Blazars: A photohadronic scenario

abstract: Blazars are a subclass of AGN and flaring in multi-TeV gamma-ray seems to be the major activity in high energy blazars a subgrup of blazars. Flaring is also unpredictable and switches between quiescent and active states involving different time scales and fluxes. While in some high energy blazars a strong temporal correlation between X-ray and multi-TeV gamma-ray has been observed, outbursts in some other have no low energy counterparts and explanation of such extreme activity needs to be addressed through different mechanisms as it is not understood well. The extragalactic background light (EBL) plays an important role in the observation of these high energy gamma-rays as it attenuates through pair production of electron-positron and also changes the spectral shape of the high energy photons. In the context of the photohadronic model and taking EBL correction into account, flaring can be explained very well. In a series of papers we have developed this model to explain multi-TeV flaring events form many blazars. Here in this review, the photohadronic model is discussed and applied to explain the multi-TeV flaring from nearby high energy blazars: Markarian 421, Markarian 501 and 1ES1959+650.


title: Volatility Analysis with Realized GARCH-Ito Models

abstract: This paper introduces a unified approach for modeling high-frequency financial data that can accommodate both the continuous-time jump-diffusion and discrete-time realized GARCH model by embedding the discrete realized GARCH structure in the continuous instantaneous volatility process. The key feature of the proposed model is that the corresponding conditional daily integrated volatility adopts an autoregressive structure where both integrated volatility and jump variation serve as innovations. We name it as the realized GARCH-Ito model. Given the autoregressive structure in the conditional daily integrated volatility, we propose a quasi-likelihood function for parameter estimation and establish its asymptotic properties. To improve the parameter estimation, we propose a joint quasi-likelihood function that is built on the marriage of daily integrated volatility estimated by high-frequency data and nonparametric volatility estimator obtained from option data. We conduct a simulation study to check the finite sample performance of the proposed methodologies and an empirical study with the S&P500 stock index and option data.


title: Multi-Cue Vehicle Detection for Semantic Video Compression In Georegistered Aerial Videos

abstract: Detection of moving objects such as vehicles in videos acquired from an airborne camera is very useful for video analytics applications. Using fast low power algorithms for onboard moving object detection would also provide region of interest-based semantic information for scene content aware image compression. This would enable more efficient and flexible communication link utilization in lowbandwidth airborne cloud computing networks. Despite recent advances in both UAV or drone platforms and imaging sensor technologies, vehicle detection from aerial video remains challenging due to small object sizes, platform motion and camera jitter, obscurations, scene complexity and degraded imaging conditions. This paper proposes an efficient moving vehicle detection pipeline which synergistically fuses both appearance and motion-based detections in a complementary manner using deep learning combined with flux tensor spatio-temporal filtering. Our proposed multi-cue pipeline is able to detect moving vehicles with high precision and recall, while filtering out false positives such as parked vehicles, through intelligent fusion. Experimental results show that incorporating contextual information of moving vehicles enables high semantic compression ratios of over 100:1 with high image fidelity, for better utilization of limited bandwidth air-to-ground network links.


title: The quantum trace as a quantum non-abelianization map

abstract: We prove that the balanced Chekhov-Fock algebra of a punctured triangulated surface is isomorphic to a skein algebra which is a deformation of the algebra of regular functions of some abelian character variety. We first deduce from this observation a classification of the irreducible representations of the balanced Chekhov-Fock algebra at odd roots of unity, which generalizes to open surfaces the classification of Bonahon, Liu and Wong. We re-interpret Bonahon and Wong's quantum trace map as a non-commutative deformation of some regular morphism between this abelian character variety and the SL2-character variety. This algebraic morphism shares many resemblance with the non-abelianization map of Gaiotto, Moore, Hollands and Neitzke. When the punctured surface is closed, we prove that this algebraic non-abelianization map induces a birational morphism between a smooth torus and the relative SL2 character variety.


title: Trees and flowers on a billiard table

abstract: In this work we completely describe the dynamics of triangle tiling billiards. In the first part of this work, we propose a geometric approach of dynamics by introducing natural foliations associated to it. In the second part, we exploit the relationship between triangle tiling billiards and a family of fully flipped $3$-interval exchange transformations on the circle. We give a combinatorial approach of dynamics via renormalization. By uniting the two approaches, we prove several conjectures on the dynamics of triangle tiling billiards. First, we prove the Tree Conjecture and the 4n+2 Conjecture, both concerning the symbolic dynamics of periodic trajectories, and both stated by Baird-Smith, Davis, Fromm and Iyer. Second, we study a family of exceptional trajectories which are closely related to the orbits of minimal Arnoux-Rauzy maps. We prove that all of these exceptional trajectories pass by all tiles, which confirms our own conjecture with P. Hubert on their non-linear escape. Moreover, we use tiling billiards to prove the convergence, up to rescaling, of arithmetic orbits of the Arnoux-Yoccoz map to the Rauzy fractal, conjectured by Hooper and Weiss. All of these conjectures have been stated in print in the last three years.


title: On the reliability of magnetic energy and helicity computations based on nonlinear force-free coronal magnetic field models

abstract: We demonstrate the sensitivity of magnetic energy and helicity computations regarding the quality of the underlying coronal magnetic field model. We apply the method of Wiegelmann & Inhester (2010) to a series of SDO/HMI vector magnetograms, and discuss nonlinear force-free (NLFF) solutions based on two different sets of the free model parameters. The two time series differ from each other concerning their force-free and solenoidal quality. Both force- and divergence-freeness are required for a consistent NLFF solution. Full satisfaction of the solenoidal property is inherent in the definition of relative magnetic helicity in order to insure gauge-independence. We apply two different magnetic helicity computation methods (Thalmann et al. 2011; Valori et al. 2012) to both NLFF time series and find that the output is highly dependent on the level to which the NLFF magnetic fields satisfy the divergence-free condition, with the computed magnetic energy being less sensitive than the relative helicity. Proxies for the non-potentiality and eruptivity derived from both quantities are also shown to depend strongly on the solenoidal property of the NLFF fields. As a reference for future applications, we provide quantitative thresholds for the force- and divergence-freeness, for the assurance of reliable computation of magnetic energy and helicity, and of their related eruptivity proxies.


title: Conservative Q-Improvement: Reinforcement Learning for an Interpretable Decision-Tree Policy

abstract: There is a growing desire in the field of reinforcement learning (and machine learning in general) to move from black-box models toward more "interpretable AI." We improve interpretability of reinforcement learning by increasing the utility of decision tree policies learned via reinforcement learning. These policies consist of a decision tree over the state space, which requires fewer parameters to express than traditional policy representations. Existing methods for creating decision tree policies via reinforcement learning focus on accurately representing an action-value function during training, but this leads to much larger trees than would otherwise be required. To address this shortcoming, we propose a novel algorithm which only increases tree size when the estimated discounted future reward of the overall policy would increase by a sufficient amount. Through evaluation in a simulated environment, we show that its performance is comparable or superior to traditional tree-based approaches and that it yields a more succinct policy. Additionally, we discuss tuning parameters to control the tradeoff between optimizing for smaller tree size or for overall reward.


title: Adaptive Partitioning Design and Analysis for Emulation of a Complex Computer Code

abstract: Computer models are used as replacements for physical experiments in a large variety of applications. Nevertheless, direct use of the computer model for the ultimate scientific objective is often limited by the complexity and cost of the model. Historically, Gaussian process regression has proven to be the almost ubiquitous choice for a fast statistical emulator for such a computer model, due to its flexible form and analytical expressions for measures of predictive uncertainty. However, even this statistical emulator can be computationally intractable for large designs, due to computing time increasing with the cube of the design size. Multiple methods have been proposed for addressing this problem. We discuss several of them, and compare their predictive and computational performance in several scenarios. We then propose solving this problem using an adaptive partitioning emulator (APE). The new approach is motivated by the idea that most computer models are only complex in particular regions of the input space. By taking a data-adaptive approach to the development of a design, and choosing to partition the space in the regions of highest variability, we obtain a higher density of points in these regions and hence accurate prediction.


title: Nonlinear spectrums of Finsler manifolds

abstract: In this paper we investigate the spectral problem in Finsler geometry. Due to the nonlinearity of the Finsler-Laplacian operator, we introduce \textit{faithful dimension pairs} by means of which the spectrum of a compact reversible Finsler metric measure manifold is defined. Various upper and lower bounds of such eigenvalues are provided in the spirit of Cheng, Buser and Gromov, which extend in several aspects the results of Hassannezhad, Kokarev and Polterovich. Moreover, we construct several faithful dimension pairs based on Lusternik-Schnirelmann category, Krasnoselskii genus and essential dimension, respectively; however, we also show that the Lebesgue covering dimension pair is not faithful. As an application, we show that the Bakry-\'Emery spectrum of a closed weighted Riemannian manifold can be characterized by the faithful Lusternik-Schnirelmann dimension pair.