title: Theory of diffusive fluctuations

abstract: The recently developed effective field theory of fluctuations around thermal equilibrium is used to compute late-time correlation functions of conserved densities. Specializing to systems with a single conservation law, we find that the diffusive pole is shifted in the presence of non-linear hydrodynamic self-interactions, and that the density-density Green's function acquires a branch point half way to the diffusive pole, at frequency $\omega= -\frac{i}{2}Dk^2$. We discuss the relevance of diffusive fluctuations for strongly correlated transport in condensed matter and cold atomic systems.


title: A Rprop-Neural-Network-Based PV Maximum Power Point Tracking Algorithm with Short-Circuit Current Limitation

abstract: This paper proposes a resilient-backpropagation-neural-network-(Rprop-NN) based algorithm for Photovoltaic (PV) maximum power point tracking (MPPT). A supervision mechanism is proposed to calibrate the Rprop-NN-MPPT reference and limit short-circuit current caused by incorrect prediction. Conventional MPPT algorithms (e.g., perturb and observe (P&O), hill climbing, and incremental conductance (Inc-Cond) etc.) are trial-and-error-based, which may result in steady-state oscillations and loss of tracking direction under fast-changing ambient environment. In addition, partial shading is also a challenge due to the difficulty of finding the global maximum power point on a multi-peak characteristic curve. As an attempt to address the aforementioned issues, a novel Rprop-NN MPPT algorithm is developed and elaborated in this work. Multiple case studies are carried out to verify the effectiveness of the proposed algorithm.


title: Blue-Noise Sampling on Graphs

abstract: In the area of graph signal processing, a graph is a set of nodes arbitrarily connected by weighted links; a graph signal is a set of scalar values associated with each node; and sampling is the problem of selecting an optimal subset of nodes from which a graph signal can be reconstructed. This paper proposes the use of spatial dithering on the vertex domain of the graph, as a way to conveniently find statistically good sampling sets. This is done establishing that there is a family of good sampling sets characterized on the vertex domain by a maximization of the distance between sampling nodes; in the Fourier domain, these are characterized by spectrums that are dominated by high frequencies referred to as blue-noise. The theoretical connection between blue-noise sampling on graphs and previous results in graph signal processing is also established, explaining the advantages of the proposed approach. Restricting our analysis to undirected and connected graphs, numerical tests are performed in order to compare the effectiveness of blue-noise sampling against other approaches.


title: Topology-Aware Surface Reconstruction for Point Clouds

abstract: We present an approach to inform the reconstruction of a surface from a point scan through topological priors. The reconstruction is based on basis functions which are optimized to provide a good fit to the point scan while satisfying predefined topological constraints. We optimize the parameters of a model to obtain likelihood function over the reconstruction domain. The topological constraints are captured by persistence diagrams which are incorporated in the optimization algorithm promote the correct topology. The result is a novel topology-aware technique which can: 1.) weed out topological noise from point scans, and 2.) capture certain nuanced properties of the underlying shape which could otherwise be lost while performing surface reconstruction. We showcase results reconstructing shapes with multiple potential topologies, compare to other classical surface construction techniques, and show the completion of real scan data.


title: Edwards curve point counting method and supersingular Edwards and Montgomery curves

abstract: We consider algebraic affine and projective curves of Edwards \cite{E, SkOdProj} over a finite field $\text{F}_{p^n}$. Most cryptosystems of the modern cryptography \cite{SkBlock} can be naturally transform into elliptic curves \cite{Kob}. We research Edwards algebraic curves over a finite field, which at the present time is one of the most promising supports of sets of points that are used for fast group operations \cite{Bir}. New method of counting Edwards curve order over finite field was constructed. It can be applied to order of elliptic curve due to birational equivalence between elliptic curve and Edwards curve. We find not only a specific set of coefficients with corresponding field characteristics, for which these curves are supersingular but also a general formula by which one can determine whether a curve $E_d[\mathbb{F}_{p^n}]$ is supersingular over this field or not. The embedding degree of the supersingular curve of Edwards over $\mathbb{F}_{p^n}$ in a finite field is investigated, the field characteristic, where this degree is minimal, was found. The criterion of supersungularity of the Edwards curves is found over $\mathbb{F}_{p^n}$. Also the generator of crypto stable sequence on an elliptic curve with a deterministic lower estimate of its period is proposed.


title: Berry-Esseen type estimate and return sequence for parabolic iteration in the upper half-plane

abstract: Two different aspects of parabolic iteration in the complex upper half-plane are considered here. First, from a noncommutative probability perspective, a Berry-Esseen type estimate for the convergence speed of the monotone central limit theorem is proved. Secondly, if the underlying measure in this central limit process is singular to the Lebesgue measure on the real line, then the iteration is shown to be an infinite-measure preserving dynamical system that has a regularly varying return sequence of index 1/2.


title: Lightweight and Efficient Image Super-Resolution with Block State-based Recursive Network

abstract: Recently, several deep learning-based image super-resolution methods have been developed by stacking massive numbers of layers. However, this leads too large model sizes and high computational complexities, thus some recursive parameter-sharing methods have been also proposed. Nevertheless, their designs do not properly utilize the potential of the recursive operation. In this paper, we propose a novel, lightweight, and efficient super-resolution method to maximize the usefulness of the recursive architecture, by introducing block state-based recursive network. By taking advantage of utilizing the block state, the recursive part of our model can easily track the status of the current image features. We show the benefits of the proposed method in terms of model size, speed, and efficiency. In addition, we show that our method outperforms the other state-of-the-art methods.


title: The inverse Voronoi problem in graphs

abstract: We introduce the inverse Voronoi diagram problem in graphs: given a graph $G$ with positive edge-lengths and a collection $\mathbb{U}$ of subsets of vertices of $V(G)$, decide whether $\mathbb{U}$ is a Voronoi diagram in $G$ with respect to the shortest-path metric. We show that the problem is NP-hard, even for planar graphs where all the edges have unit length. We also study the parameterized complexity of the problem and show that the problem is W[1]-hard when parameterized by the number of Voronoi cells or by the pathwidth of the graph. For trees we show that the problem can be solved in $O(N+n \log^2 n)$ time, where $n$ is the number of vertices in the tree and $N=n+\sum_{U\in \mathbb{U}}|U|$ is the size of the description of the input. We also provide a lower bound of $\Omega(n \log n)$ time for trees with $n$ vertices.


title: Improved bounds for Hadwiger's covering problem via thin-shell estimates

abstract: A central problem in discrete geometry, known as Hadwiger's covering problem, asks what the smallest natural number $N\left(n\right)$ is such that every convex body in ${\mathbb R}^{n}$ can be covered by a union of the interiors of at most $N\left(n\right)$ of its translates. Despite continuous efforts, the best general upper bound known for this number remains as it was more than sixty years ago, of the order of ${2n \choose n}n\ln n$. In this note, we improve this bound by a sub-exponential factor. That is, we prove a bound of the order of ${2n \choose n}e^{-c\sqrt{n}}$ for some universal constant $c>0$. Our approach combines ideas from previous work by Artstein-Avidan and the second named author with tools from Asymptotic Geometric Analysis. One of the key steps is proving a new lower bound for the maximum volume of the intersection of a convex body $K$ with a translate of $-K$; in fact, we get the same lower bound for the volume of the intersection of $K$ and $-K$ when they both have barycenter at the origin. To do so, we make use of measure concentration, and in particular of thin-shell estimates for isotropic log-concave measures. Using the same ideas, we establish an exponentially better bound for $N\left(n\right)$ when restricting our attention to convex bodies that are $\psi_{2}$. By a slightly different approach, an exponential improvement is established also for classes of convex bodies with positive modulus of convexity.


title: Holographic Complexity for Defects Distinguishes Action from Volume

abstract: We explore the two holographic complexity proposals for the case of a 2d boundary CFT with a conformal defect. We focus on a Randall-Sundrum type model of a thin AdS$_2$ brane embedded in AdS$_3$. We find that, using the "complexity=volume" proposal, the presence of the defect generates a logarithmic divergence in the complexity of the full boundary state with a coefficient which is related to the central charge and to the boundary entropy. For the "complexity=action" proposal we find that the complexity is not influenced by the presence of the defect. This is the first case in which the results of the two holographic proposals differ so dramatically. We consider also the complexity of the reduced density matrix for subregions enclosing the defect. We explore two bosonic field theory models which include two defects on opposite sides of a periodic domain. We point out that for a compact boson, current free field theory definitions of the complexity would have to be generalized to account for the effect of zero-modes.


title: Site and orbital selective correlations in $\beta$-Pu

abstract: We investigate the electronic structure of the highly anisotropic $\beta$ phase of metallic plutonium, within the combination of density functional theory (DFT) and dynamical mean field theory (DMFT). Its crystal structure gives rise to site and orbital selective electronic correlations, with coherent Pu-5$f_{5/2}$ states and very incoherent Pu-5$f_{7/2}$ states. The Hund's coupling is essential for determining the level of correlations of electrons in Pu-5$f$ states, and for the quasiparticle multiplets features in the Pu-5$f$ spectral function.


title: Identification of kinetic order parameters for non-equilibrium dynamics

abstract: A popular approach to analyze the dynamics of high-dimensional many-body systems, such as macromolecules, is to project the trajectories onto a space of slowly-varying collective variables, where subsequent analyses are made, such as clustering or estimation of free energy profiles or Markov state models (MSMs). However, existing "dynamical" dimension reduction methods, such as the time-lagged independent component analysis (TICA) are only valid if the dynamics obeys detailed balance (microscopic reversibility) and typically require long, equilibrated simulation trajectories. Here we develop a dimension reduction method for non-equilibrium dynamics based on the recently developed Variational Approach for Markov Processes (VAMP) by Wu and No\'e. VAMP is illustrated by obtaining a low-dimensional description of a single file ion diffusion model and by identifying long-lived states from molecular dynamics simulations of the KcsA channel protein in an external electrochemical potential. This analysis provides detailed insights into the coupling of conformational dynamics, the configuration of the selectivity filter, and the conductance of the channel. We recommend VAMP as a replacement for the less general TICA method.


title: Double bit in-plane magnetic skyrmions on a track

abstract: A magnetic skyrmion isusually refers to a twisted spin texture surrounded by uniformly aligned out-of-plane spinsin the background of a uniformly magnetized state. The invariance of the magnetic skyrmion conserves its topological charge under any continuous transformations of the spin textures, leads to which represents the robustness of a magnetic skyrmion their texture against external perturbations, making it ideal to use skyrmions as . Such a behaviour is required for an ideal information carriers. To date, most magnetic skyrmion studies have been performedfocused on in perpendicularly magnetized systems, where the skyrmion topological number is determined by the relative orientation between the core and /outer perpendicular magnetization directions of the skyrmion is either up/down or down/up. Here we show that there also exists a new type of magnetic skyrmion with surrounding spins to be uniformly aligned to the in-plane direction. By continuous transformation and relaxation of the spin textures of out-of-plane skyrmions, we showed that an in-plane skyrmion, where the background magnetization is in the in-plane direction, is also possible. Different from the conventional perpendicular magnetic skyrmionsContrary to the skyrmions in a perpendicularly magnetized state, the in-plane magnetic skyrmions with opposite signs of topological charges can inherentlycharges can inherently coexist in the in-plane magnetization system coexist. Moreover, the iIn-plane skyrmions of opposite charge can move together by an electric current and exhibit with an opposite spin Hall effect. These findings showdemonstrate the inherent possibility of a double-bit transfer in a single magnetic wire that is not possible in a perpendicularly magnetized system.


title: Development of a Multi-Agent System for Optimal Sizing of a Commercial Complex Microgrid

abstract: In this paper, a novel intelligent method based on a multi-agent system (MAS) is applied to the problem of optimal sizing in a stand-alone office complex microgrid such that the electricity demand of the office building and the charging demand of the plug-in hybrid electric vehicle (PHEV) charging station are met. The proposed MAS-based architecture consists of five different agents, namely generation agent (GA), electrical load agent (LA), PHEV charging station agent (SA), control agent (CA), and design agent (DA) which are organized in three levels. In the proposed MAS, control agent coordinates the interactions between the generation agent and electrical load and charging station agents and following a request from design agent, sends the information on the operation of the microgrid with determined sizes to the design agent. According to the received information, design agent finds the optimal sizes of the system's components such that the electricity and charging demands are met considering two reliability indices to investment decisions. This study is performed for Kish Island in south Iran.


title: Fast Algorithms for Knapsack via Convolution and Prediction

abstract: The \Problem{knapsack} problem is a fundamental problem in combinatorial optimization. It has been studied extensively from theoretical as well as practical perspectives as it is one of the most well-known NP-hard problems. The goal is to pack a knapsack of size $t$ with the maximum value from a collection of $n$ items with given sizes and values. Recent evidence suggests that a classic $O(nt)$ dynamic-programming solution for the \Problem{knapsack} problem might be the fastest in the worst case. In fact, solving the \Problem{knapsack} problem was shown to be computationally equivalent to the \Problem{$(\min, +)$ convolution} problem, which is thought to be facing a quadratic-time barrier. This hardness is in contrast to the more famous \Problem{$(+, \cdot)$ convolution} (generally known as \Problem{polynomial multiplication}), that has an $O(n\log n)$-time solution via Fast Fourier Transform. Our main results are algorithms with near-linear running times (in terms of the size of the knapsack and the number of items) for the \Problem{knapsack} problem, if either the values or sizes of items are small integers. More specifically, if item sizes are integers bounded by $\smax$, the running time of our algorithm is $\tilde O((n+t)\smax)$. If the item values are integers bounded by $\vmax$, our algorithm runs in time $\tilde O(n+t\vmax)$. Best previously known running times were $O(nt)$, $O(n^2\smax)$ and $O(n\smax\vmax)$ (Pisinger, J. of Alg., 1999).


title: Ensemble Bayesian Decision Making with Redundant Deep Perceptual Control Policies

abstract: This work presents a novel ensemble of Bayesian Neural Networks (BNNs) for control of safety-critical systems. Decision making for safety-critical systems is challenging due to performance requirements with significant consequences in the event of failure. In practice, failure of such systems can be avoided by introducing redundancies of control. Neural Networks (NNs) are generally not used for safety-critical systems as they can behave in unexpected ways in response to novel inputs. In addition, there may not be any indication as to when they will fail. BNNs have been recognized for their ability to produce not only viable outputs but also provide a measure of uncertainty in these outputs. This work combines the knowledge of prediction uncertainty obtained from BNNs and ensemble control for a redundant control methodology. Our technique is applied to an agile autonomous driving task. Multiple BNNs are trained to control a vehicle in an end-to-end fashion on different sensor inputs provided by the system. We show that an individual network is successful in maneuvering around the track but crashes in the presence of unforeseen input noise. Our proposed ensemble of BNNs shows successful task performance even in the event of multiple sensor failures.


title: How to Organize your Deep Reinforcement Learning Agents: The Importance of Communication Topology

abstract: In this empirical paper, we investigate how learning agents can be arranged in more efficient communication topologies for improved learning. This is an important problem because a common technique to improve speed and robustness of learning in deep reinforcement learning and many other machine learning algorithms is to run multiple learning agents in parallel. The standard communication architecture typically involves all agents intermittently communicating with each other (fully connected topology) or with a centralized server (star topology). Unfortunately, optimizing the topology of communication over the space of all possible graphs is a hard problem, so we borrow results from the networked optimization and collective intelligence literatures which suggest that certain families of network topologies can lead to strong improvements over fully-connected networks. We start by introducing alternative network topologies to DRL benchmark tasks under the Evolution Strategies paradigm which we call Network Evolution Strategies. We explore the relative performance of the four main graph families and observe that one such family (Erdos-Renyi random graphs) empirically outperforms all other families, including the de facto fully-connected communication topologies. Additionally, the use of alternative network topologies has a multiplicative performance effect: we observe that when 1000 learning agents are arranged in a carefully designed communication topology, they can compete with 3000 agents arranged in the de facto fully-connected topology. Overall, our work suggests that distributed machine learning algorithms would learn more efficiently if the communication topology between learning agents was optimized.


title: Deep Multi-Agent Reinforcement Learning with Relevance Graphs

abstract: Over recent years, deep reinforcement learning has shown strong successes in complex single-agent tasks, and more recently this approach has also been applied to multi-agent domains. In this paper, we propose a novel approach, called MAGnet, to multi-agent reinforcement learning (MARL) that utilizes a relevance graph representation of the environment obtained by a self-attention mechanism, and a message-generation technique inspired by the NerveNet architecture. We applied our MAGnet approach to the Pommerman game and the results show that it significantly outperforms state-of-the-art MARL solutions, including DQN, MADDPG, and MCTS.


title: Super-Keplerian Equatorial Outflows in SS 433. Centrifugal Ejection of the Circumbinary Disk

abstract: The microquasar SS 433 is the only known steady supercritical accretor in the Galaxy. It is well-known for its relativistic baryonic jets, but the system also drives equatorial outflows. These have been routinely detected in radio images, and components associated with a circumbinary disk have also been suggested in optical emission lines. We wish to spatially resolve the regions producing the stationary emission lines of SS 433 to shed light on its circumbinary structure and outflows. With an estimated binary orbit size $\lesssim 0.1 \text{ mas}$, this requires optical interferometry. We use the optical interferometer VLTI+GRAVITY to spatially resolve SS 433 in the near-infrared K band at high spectral resolution ($R\approx 4000$) on three nights in July 2017. The stationary Br$\gamma$ line is clearly dominated by an extended $\sim 1 \text{ mas} \sim 5 \text{ AU}$ circumbinary structure perpendicular to the jets and with a strong rotation component. The rotation direction is retrograde relative to the jet precession, in accordance with the slaved disk precession model. The structure has a very high specific angular momentum and is too extended to be a stable circumbinary disk in Keplerian rotation; interpreting it as such leads to a very high enclosed mass $M \gtrsim 400 M_{\odot}$. We instead interpret it as the centrifugal ejection of the circumbinary disk, with the implication that there must be an efficient transfer of specific angular momentum from the binary to the disk. We suggest that the equatorial outflows sometimes seen in radio images result from similar episodes of circumbinary disk centrifugal ejection. In addition to the equatorial structure, we find a very extended $\sim 6 \text{ mas} \sim 30 \text{ AU}$ spherical wind component to the Br$\gamma$ line: the entire binary is engulfed in an optically thin spherical line emission envelope.


title: An a.e. lower bound for Hausdorff dimension under vertical projections in the Heisenberg group

abstract: An improved a.e. lower bound is given for Hausdorff dimension under vertical projections in the first Heisenberg group.


title: An Introduction to Deep Reinforcement Learning

abstract: Deep reinforcement learning is the combination of reinforcement learning (RL) and deep learning. This field of research has been able to solve a wide range of complex decision-making tasks that were previously out of reach for a machine. Thus, deep RL opens up many new applications in domains such as healthcare, robotics, smart grids, finance, and many more. This manuscript provides an introduction to deep reinforcement learning models, algorithms and techniques. Particular focus is on the aspects related to generalization and how deep RL can be used for practical applications. We assume the reader is familiar with basic machine learning concepts.


title: The correlation between isotropic energy and duration of gamma-ray bursts

abstract: In this paper, we study the correlation between isotropic energy and duration of gamma-ray bursts (GRBs) for the first time. The correlation is found to be $T_d \propto {E_{iso}}^{0.34\pm 0.03}$ from the {\em Swift} GRB sample. After comparing with solar flares from {\em RHESSI} and stellar superflares from {\em Kepler} satellite, we find that the correlation of GRBs shows similar exponent with those of solar flares and stellar superflares. Inspired by the physical mechanism of solar flares and stellar superflares which is magnetic reconnection, we interpret the correlation using magnetic reconnection theory. This similarity hints that magnetic reconnection may dominate energy releasing process of GRBs.


title: On least squares problems with certain Vandermonde--Khatri--Rao structure with applications to DMD

abstract: This paper proposes a new computational method for solving structured least squares problems that arise in the process of identification of coherent structures in fluid flows. It is deployed in combination with dynamic mode decomposition (DMD) which provides a non-orthogonal set of modes --- corresponding to particular temporal frequencies --- a subset of which is used to represent time snapshots of the underlying dynamics. The coefficients of the representation are determined from a solution of a structured linear least squares problem with the matrix that involves the Khatri--Rao product of a triangular and a Vandermonde matrix. Such a structure allows a very efficient normal equation based least squares solution, which is used in state of the art CFD tools such as the sparsity promoting DMD (DMDSP). A new numerical analysis of the normal equations approach provides insights about its applicability and its limitations. Relevant condition numbers that determine numerical robustness are identified and discussed. Further, the paper offers a corrected semi-normal solution and QR factorization based algorithms. It is shown how to use the Vandermonde--Khatri--Rao structure to efficiently compute the QR factorization of the least squares coefficient matrix, thus providing a new computational tool for the ill-conditioned cases where the normal equations may fail to compute a sufficiently accurate solution. Altogether, the presented material provides a firm numerical linear algebra framework for a class of structured least squares problems arising in a variety of applications.


title: Deep Multimodal Learning: An Effective Method for Video Classification

abstract: Videos have become ubiquitous on the Internet. And video analysis can provide lots of information for detecting and recognizing objects as well as help people understand human actions and interactions with the real world. However, facing data as huge as TB level, effective methods should be applied. Recurrent neural network (RNN) architecture has wildly been used on many sequential learning problems such as Language Model, Time-Series Analysis, etc. In this paper, we propose some variations of RNN such as stacked bidirectional LSTM/GRU network with attention mechanism to categorize large-scale video data. We also explore different multimodal fusion methods. Our model combines both visual and audio information on both video and frame level and received great result. Ensemble methods are also applied. Because of its multimodal characteristics, we decide to call this method Deep Multimodal Learning(DML). Our DML-based model was trained on Google Cloud and our own server and was tested in a well-known video classification competition on Kaggle held by Google.


title: The Collimated Radiation in SS 433. Constraints from Spatially Resolved Optical Jets and $\texttt{Cloudy}$ Modeling of the Optical Bullets

abstract: The microquasar SS 433 is well-known for its precessing, relativistic baryonic jets. Depending on their heating mechanism, the optical jet bullets may serve as a probe of the collimated radiation coming from the inner region close to the compact object. The optical interferometer VLTI/GRAVITY has allowed to spatially resolved the optical jets in SS 433 for the first time. We present here the second such observation taken over three nights in July 2017. In addition, we use the multi-wavelength XSHOOTER spectrograph at VLT to study the optical bullets in SS 433 in detail. GRAVITY reveals elongated exponential-like spatial profiles for the optical jets, suggestive of a heating mechanism acting throughout a long portion of the jet and naturally explained by photoionization by the collimated radiation. We also spatially resolve the movement of the optical bullets for the first time, detecting extended jet components corresponding to previous ejections. \texttt{Cloudy} photoionization models explain both the spatial intensity profiles measured with GRAVITY and the line ratios from XSHOOTER, and constrain the properties of the optical bullets and the ionizing radiation. We find that the latter must peak in the UV with an isotropic luminosity (as inferred by a face on observer) $\approx 10^{41}$ erg/s. Provided that the X-ray SED is sufficiently hard, the collimated X-ray luminosity could still be high enough so that the face on observer would see SS 433 as an ULX ($L_X \lesssim 10^{40}$ erg/s) and it would still be compatible with the H/He/He+ ionization balance of the optical bullets. The kinetic power in the optical jets is constrained to $3-20 \times 10^{38}$ erg/s, and the extinction in the optical jets to $A_V = 6.7 \pm 0.1$. We suggest there may be substantial $A_V \gtrsim 1$ and structured circumstellar extinction in SS 433, likely arising from dust formed in equatorial outflows.


title: Eigenvalue Corrected Noisy Natural Gradient

abstract: Variational Bayesian neural networks combine the flexibility of deep learning with Bayesian uncertainty estimation. However, inference procedures for flexible variational posteriors are computationally expensive. A recently proposed method, noisy natural gradient, is a surprisingly simple method to fit expressive posteriors by adding weight noise to regular natural gradient updates. Noisy K-FAC is an instance of noisy natural gradient that fits a matrix-variate Gaussian posterior with minor changes to ordinary K-FAC. Nevertheless, a matrix-variate Gaussian posterior does not capture an accurate diagonal variance. In this work, we extend on noisy K-FAC to obtain a more flexible posterior distribution called eigenvalue corrected matrix-variate Gaussian. The proposed method computes the full diagonal re-scaling factor in Kronecker-factored eigenbasis. Empirically, our approach consistently outperforms existing algorithms (e.g., noisy K-FAC) on regression and classification tasks.


title: On a class of anharmonic oscillators

abstract: In this work we study a class of anharmonic oscillators within the framework of the Weyl-H\"ormander calculus. The anharmonic oscillators arise from several applications in mathematical physics as natural extensions of the harmonic oscillator. A prototype is an operator on $\mathbb{R}^n$ of the form $(-\Delta)^{\ell}+|x|^{2k}$ for $k,\ell$ integers $\geq 1$. The simplest case corresponds to Hamiltonians of the form $|\xi|^2+|x|^{2k}$. Here by associating a H\"ormander metric $g$ to a given anharmonic oscillator we investigate several properties of the anharmonic oscillators. We obtain spectral properties in terms of Schatten-von Neumann classes for their negative powers. We also study some examples of anharmonic oscillators arising from the analysis on Lie groups


title: Stochastic Dynamics for Earthquake Ruptures

abstract: In this paper, we propose a stochastic dynamic model for earthquake rupture and suggest that the Langevin equation of frictions may be used for interpreting the slip distributions of rupture processes in earthquakes. The steady-state solution of the derived Langevin equation analytically attains the truncated exponential (TEX) distribution that is empirically characterized in many rupture models of earthquake events worldwide, as demonstrated by Thingbaijam and Mai (2016, https://doi.org/10.1785/0120150291). Our proposed stochastic dynamic faulting model for earthquake rupture intrinsically includes fluctuations and uncertainties in the heterogeneity of faulting planes as random variables. Specifically, we related the characteristic parameter u_c in TEX functions to the ratio of diffusion and friction coefficients D and \gamma of the Langevin equation.


title: Parallelizing greedy for submodular set function maximization in matroids and beyond

abstract: We consider parallel, or low adaptivity, algorithms for submodular function maximization. This line of work was recently initiated by Balkanski and Singer and has already led to several interesting results on the cardinality constraint and explicit packing constraints. An important open problem is the classical setting of matroid constraint, which has been instrumental for developments in submodular function maximization. In this paper we develop a general strategy to parallelize the well-studied greedy algorithm and use it to obtain a randomized $\left(\frac{1}{2} - \epsilon\right)$-approximation in $\operatorname{O}\left( \frac{\log^2 n}{\epsilon^2} \right)$ rounds of adaptivity. We rely on this algorithm, and an elegant amplification approach due to Badanidiyuru and Vondr\'ak to obtain a fractional solution that yields a near-optimal randomized $\left( 1 - 1/e - \epsilon \right)$-approximation in $O\left( {\frac{\log^2 n}{\epsilon^3}} \right) $ rounds of adaptivity. For non-negative functions we obtain a $\left( {3-2\sqrt{2}}\right)$-approximation and a fractional solution that yields a $\left( {\frac{1}{e} - \epsilon}\right)$-approximation. Our approach for parallelizing greedy yields approximations for intersections of matroids and matchoids, and the approximation ratios are comparable to those known for sequential greedy.


title: Are All Training Examples Created Equal? An Empirical Study

abstract: Modern computer vision algorithms often rely on very large training datasets. However, it is conceivable that a carefully selected subsample of the dataset is sufficient for training. In this paper, we propose a gradient-based importance measure that we use to empirically analyze relative importance of training images in four datasets of varying complexity. We find that in some cases, a small subsample is indeed sufficient for training. For other datasets, however, the relative differences in importance are negligible. These results have important implications for active learning on deep networks. Additionally, our analysis method can be used as a general tool to better understand diversity of training examples in datasets.


title: Light vector mesons ($\omega$, $\rho$ and $\phi$) in strong magnetic fields: A QCD sum rule approach

abstract: The mass modifications of the light vector mesons ($\omega$, $\rho$ and $\phi$) are investigated in asymmetric nuclear matter in the presence of strong magnetic fields, using a QCD sum rule approach. These are computed from the medium modifications of the non-strange and strange light quark condensates as well as scalar gluon condensate. The quark and gluon condensates are calculated from the medium changes of the scalar fields (non-strange and strange) and a scalar dilaton field in the magnetized nuclear matter, within a chiral SU(3) model. The scalar dilaton field within the model breaks the scale invariance of QCD and simulates the gluon condensate. The anomalous magnetic moments for the nucleons are taken into account in the present study.


title: Testing S-duality with non-orientable surfaces

abstract: Kapustin and Witten showed that a twisted version of N=4 gauge theory in four dimensions compactifies to a two-dimensional sigma-model whose target space is the Hitchin moduli space. In this talk, I consider the reduction of the gauge theory on a four dimensional orientable spacetime manifold which is not a global product of two surfaces but contains embedded non-orientable surfaces. The low energy theory is a sigma-model on a two dimensional worldsheet whose boundary components end on branes constructed from the Hitchin moduli space associated to a non-orientable surface. I will also compare the discrete topological fluxes in four and two dimensional theories and verify the mirror symmetry on branes as predicted by the S-duality in gauge theory. This provides another non-trivial test of $S$-duality using reduction along possibly non-orientable surfaces. Finally, I consider the quantisation of the Hitchin moduli space from a non-orientable surface as an example of quantisation via branes and mirror symmetry.


title: Decoupling between thermodynamics and dynamics during rejuvenation in colloidal glasses

abstract: We rejuvenate well-aged quasi-2D binary colloidal glasses by thermal cycling, and systematically measure both the statistical responses and particle-level structural evolutions during rejuvenation. While the moduli and boson peak are continuously rejuvenated with increasing number of cycles, the mean square displacement (MSD) fluctuates significantly between different groups of thermal cycles. The decoupling between the thermodynamical and dynamical evolutions suggests different microscopic origins for different bulk properties of glasses. We find that a small fraction of structural rearrangements triggered by thermal cycling could alter the whole elastic continuum and lead to the significant thermodynamic rejuvenation, while localized defects could be activated and deactivated at the positions close to the rearrangements with significantly high mobility change and hence result in the fluctuated dynamics even with only about 10% of particles as fast regions. Our results offer a comprehensive picture for the microscopic mechanisms underlying bulk glass rejuvenation, which could be readily used to refine glass properties or to formulate further statistical theories in glassy systems.

categories: cond-mat.soft cond-mat.dis-nn cond-mat.mtrl-sci cond-mat.stat-mech


title: ContextServ: Towards Model-Driven Development of Context-AwareWeb Services

abstract: In the era of Web of Things and Services, Context-aware Web Services (CASs) are emerging as an important technology for building innovative context-aware applications. CASs enable the information integration from both the physical and virtual world, which affects human living. However, it is challenging to build CASs, due to the lack of context provisioning management approach and limited generic approach for formalizing the development process. We therefore propose ContextServ, a platform that uses a model-driven approach to support the full life cycle of CASs development, hence offering significant design and management flexibility. ContextServ implements a proposed UML-based modelling language ContextUML to support multiple modelling languages. It also supports dynamic adaptation of WS-BPEL based context-aware composite services by weaving context-aware rules into the process. Extensive experimental evaluations on ContextServ and its components showcase that ContextServ can support effective development and efficient execution of context-aware Web services.


title: The role of protection zone on species spreading governed by a reaction-diffusion model with strong Allee effect

abstract: It is known that a species dies out in the long run for small initial data if its evolution obeys a reaction of bistable nonlinearity. Such a phenomenon, which is termed as the strong Allee effect, is well supported by numerous evidence from ecosystems, mainly due to the environmental pollution as well as unregulated harvesting and hunting. To save an endangered species, in this paper we introduce a protection zone that is governed by a Fisher-KPP nonlinearity, and examine the dynamics of a reaction-diffusion model with strong Allee effect and protection zone. We show the existence of two critical values $0<L_*\leq L^*$, and prove that a vanishing-transition-spreading trichotomy result holds when the length of protection zone is smaller than $L_*$; a transition-spreading dichotomy result holds when the length of protection zone is between $L_*$ and $L^*$; only spreading happens when the length of protection zone is larger than $L^*$. This suggests that the protection zone works when its length is larger than the critical value $L_*$. Furthermore, we compare two types of protection zone with the same length: a connected one and a separate one, and our results reveal that the former is better for species spreading than the latter.


title: Constant gap between conventional strategies and those based on C*-dynamics for self-embezzlement

abstract: We consider a bipartite transformation that we call self-embezzlement and use it to prove a constant gap between the capabilities of two models of quantum information: the conventional model, where bipartite systems are represented by tensor products of Hilbert spaces; and a natural model of quantum information processing for abstract states on C*-algebras, where joint systems are represented by tensor products of C*-algebras. We call this the C*-circuit model and show that it is a special case of the commuting-operator model (in that it can be translated into such a model). For the conventional model, we show that there exists a constant $\epsilon_0 > 0$ such that self-embezzlement cannot be achieved with precision parameter less than $\epsilon_0$ (i.e., the fidelity cannot be greater than $1 - \epsilon_0$); whereas, in the C*-circuit model -- as well as in a commuting-operator model -- the precision can be $0$ (i.e., fidelity~$1$).


title: Parametric Timed Broadcast Protocols

abstract: In this paper we consider state reachability in networks composed of many identical processes running a parametric timed broadcast protocol (PTBP). PTBP are a new model extending both broadcast protocols and parametric timed automata. This work is, up to our knowledge, the first to consider the combination of both a parametric network size and timing parameters in clock guard constraints. Since the communication topology is of utmost importance in broadcast protocols, we investigate reachability problems in both clique semantics where every message reaches every processes, and in reconfigurable semantics where the set of receivers is chosen non-deterministically. In addition, we investigate the decidability status depending on whether the timing parameters in guards appear only as upper bounds in guards, as lower bounds or when the set of parameters is partitioned in lower-bound and upper-bound parameters.


title: Arc closures and the local isomorphism problem

abstract: We give an answer in the "geometric" setting to a question of de Fernex, Ein, and Ishii, asking when local isomorphisms of $k$-schemes can be detected on the associated maps of local arc or jet schemes. In particular, we show that their ideal-closure operation $\mathfrak a\mapsto \mathfrak a^{\mathrm{ac}}$ (the arc-closure) on a local $k$-algebra $(R,\mathfrak m,L)$ is trivial when $R$ is Noetherian and $k\hookrightarrow L$ is separable, and thus that such a germ $\mathrm{Spec}\ R$ has the (embedded) local isomorphism property.


title: Using ATLAS@Home to exploit extra CPU from busy grid sites

abstract: Grid computing typically provides most of the data processing resources for large High Energy Physics experiments. However typical grid sites are not fully utilized by regular workloads. In order to increase the CPU utilization of these grid sites, the ATLAS@Home volunteer computing framework can be used as a backfilling mechanism. Results show an extra 15% to 42% of CPU cycles can be exploited by backfilling grid sites running regular workloads while the overall CPU utilization can remain over 90%. Backfilling has no impact on the failure rate of the grid jobs, and the impact on the CPU efficiency of grid jobs varies from 1% to 11% depending on the configuration of the site. In addition the throughput of backfill jobs in terms of CPU time per simulated event is the same as for resources dedicated to ATLAS@Home. This approach is sufficiently generic that it can easily be extended to other clusters.


title: Inverse obstacle scattering problem for elastic waves with phased or phaseless far-field data

abstract: This paper concerns an inverse elastic scattering problem which is to determine the location and the shape of a rigid obstacle from the phased or phaseless far-field data for a single incident plane wave. By introducing the Helmholtz decomposition, the model problem is reduced to a coupled boundary value problem of the Helmholtz equations. The relation is established between the compressional or shear far-field pattern for the elastic wave equation and the corresponding far-field pattern for the coupled Helmholtz equations. An efficient and accurate Nystr\"{o}m type discretization for the boundary integral equation is developed to solve the coupled system. The translation invariance of the phaseless compressional and shear far-field patterns are proved. A system of nonlinear integral equations is proposed and two iterative reconstruction methods are developed for the inverse problem. In particular, for the phaseless data, a reference ball technique is introduced to the scattering system in order to break the translation invariance. Numerical experiments are presented to demonstrate the effectiveness and robustness of the proposed method.


title: varRhoTurbVOF: a new set of volume of fluid solvers for turbulent isothermal multiphase flows in OpenFOAM

abstract: The volume of fluid (VOF) method is a popular approach for multiphase flow modeling. The open-source computational fluid dynamics (CFD) software, OpenFOAM, implements a variety of VOF-based solvers and provides users a wide range of turbulence models. Since isothermal multiphase flows under the VOF framework belong to the variable-density incompressible flow category, the isothermal VOF-based solvers in OpenFOAM fail to use the correct turbulence models. varRhoTurbVOF is designed to solve this issue and with the hope to replace all the corresponding existing solvers in the future. With the object-oriented paradigm, varRhoTurbVOF guarantees the usability, reusability and maintainability of the codes. Aside from turbulence modeling, all other features in the original solvers are preserved in varRhoTurbVOF.


title: Knotting statistics for polygons in lattice tubes

abstract: We study several related models of self-avoiding polygons in a tubular subgraph of the simple cubic lattice, with a particular interest in the asymptotics of the knotting statistics. Polygons in a tube can be characterised by a finite transfer matrix, and this allows for the derivation of pattern theorems, calculation of growth rates and exact enumeration. We also develop a static Monte Carlo method which allows us to sample polygons of a given size directly from a chosen Boltzmann distribution. Using these methods we accurately estimate the growth rates of unknotted polygons in the $2\times1\times\infty$ and $3\times1\times\infty$ tubes, and confirm that these are the same for any fixed knot-type $K$. We also confirm that the entropic exponent for unknots is the same as that of all polygons, and that the exponent for fixed knot-type $K$ depends only on the number of prime factors in the knot decomposition of $K$. For the simplest knot-types, this leads to a good approximation for the polygon size at which the probability of the given knot-type is maximized, and in some cases we are able to sample sufficiently long polygons to observe this numerically.

categories: cond-mat.stat-mech math-ph math.CO math.GN math.MP


title: A Pseudospectral Approach to High Index DAE Optimal Control Problems

abstract: Historically, solving optimal control problems with high index differential algebraic equations (DAEs) has been considered extremely hard. Computational experience with Runge-Kutta (RK) methods confirms the difficulties. High index DAE problems occur quite naturally in many practical engineering applications. Over the last two decades, a vast number of real-world problems have been solved routinely using pseudospectral (PS) optimal control techniques. In view of this, we solve a "provably hard," index-three problem using the PS method implemented in DIDO, a state-of-the-art MATLAB optimal control toolbox. In contrast to RK-type solution techniques, no laborious index-reduction process was used to generate the PS solution. The PS solution is independently verified and validated using standard industry practices. It turns out that proper PS methods can indeed be used to "directly" solve high index DAE optimal control problems. In view of this, it is proposed that a new theory of difficulty for DAEs be put forth.


title: Rethinking clinical prediction: Why machine learning must consider year of care and feature aggregation

abstract: Machine learning for healthcare often trains models on de-identified datasets with randomly-shifted calendar dates, ignoring the fact that data were generated under hospital operation practices that change over time. These changing practices induce definitive changes in observed data which confound evaluations which do not account for dates and limit the generalisability of date-agnostic models. In this work, we establish the magnitude of this problem on MIMIC, a public hospital dataset, and showcase a simple solution. We augment MIMIC with the year in which care was provided and show that a model trained using standard feature representations will significantly degrade in quality over time. We find a deterioration of 0.3 AUC when evaluating mortality prediction on data from 10 years later. We find a similar deterioration of 0.15 AUC for length-of-stay. In contrast, we demonstrate that clinically-oriented aggregates of raw features significantly mitigate future deterioration. Our suggested aggregated representations, when retrained yearly, have prediction quality comparable to year-agnostic models.


title: Blowing up extremal Poincar\'e type manifolds

abstract: We prove a version of the Arezzo-Pacard-Singer blow-up theorem in the setting of Poincar\'e type metrics. We apply this to give new examples of extremal Poincar\'e type metrics. A key feature is an additional obstruction which has no analogue in the compact case. This condition is conjecturally related to ensuring the metrics remain of Poincar\'e type.


title: On the Letelier spacetime with quintessence: solution, thermodynamics and Hawking radiation

abstract: We obtain the solution corresponding to a static and spherically symmetric black hole with a cloud of strings (Letelier spacetime) immersed in a quintessential fluid. We discuss some aspects of its thermodynamics. We also present a discussion about Hawking radiation of particles, in the background under consideration.


title: A Model Problem for Nematic-Isotropic Transitions with Highly Disparate Elastic Constants

abstract: We analyze a model problem based on highly disparate elastic constants that we propose in order to understand corners and cusps that form on the boundary between the nematic and isotropic phases in a liquid crystal. For a bounded planar domain $\Omega$ we investigate the $\varepsilon \to 0$ asymptotics of the variational problem \[\inf \frac{1}{2}\int_\Omega \left( \frac{1}{\varepsilon} W(u)+\varepsilon |\nabla u|^2 + L_\varepsilon(\mathrm{div}\, u)^2 \right) \,dx\] within various parameter regimes for $L_\varepsilon > 0.$ Here $u:\Omega\to\mathbb{R}^2$ and $W$ is a potential vanishing on the unit circle and at the origin. When $\varepsilon\ll L_\varepsilon\to 0$, we show that these functionals $\Gamma-$converge to a constant multiple of the perimeter of the phase boundary and the divergence penalty is not felt. However, when $L_\varepsilon \equiv L > 0$, we find that a tangency requirement along the phase boundary for competitors in the conjectured $\Gamma$-limit becomes a mechanism for development of singularities. We establish criticality conditions for this limit and under a non-degeneracy assumption on the potential we prove compactness of energy bounded sequences in $L^2$. The role played by this tangency condition on the formation of interfacial singularities is investigated through several examples: each of these examples involves analytically rigorous reasoning motivated by numerical experiments. We argue that generically, "wall" singularities between $\mathbb{S}^1$-valued states are expected near the defects along the phase boundary.


title: Restricted Boltzmann Machine with Multivalued Hidden Variables: a model suppressing over-fitting

abstract: Generalization is one of the most important issues in machine learning problems. In this study, we consider generalization in restricted Boltzmann machines (RBMs). We propose an RBM with multivalued hidden variables, which is a simple extension of conventional RBMs. We demonstrate that the proposed model is better than the conventional model via numerical experiments for contrastive divergence learning with artificial data and a classification problem with MNIST.


title: Two-term relative cluster tilting subcategories, $\tau-$tilting modules and silting subcategories

abstract: Let $\mathcal{C}$ be a triangulated category with shift functor $[1]$ and $\mathcal{R}$ a rigid subcategory of $\mathcal{C}$. We introduce the notions of two-term $\mathcal{R}[1]$-rigid subcategories, two-term (weak) $\mathcal{R}[1]$-cluster tilting subcategories and two-term maximal $\mathcal{R}[1]$-rigid subcategories, and discuss relationship between them. Our main result shows that there exists a bijection between the set of two-term $\mathcal{R}[1]$-rigid subcategories of $\mathcal{C}$ and the set of $\tau$-rigid subcategories of $\mod\mathcal{R}$, which induces a one-to-one correspondence between the set of two-term weak $\mathcal{R}[1]$-cluster tilting subcategories of $\mathcal{C}$ and the set of support $\tau$-tilting subcategories of $\mod\mathcal{R}$. This generalizes the main results in \cite{YZZ} where $\mathcal{R}$ is a cluster tilting subcategory. When $\mathcal{R}$ is a silting subcategory, we prove that the two-term weak $\mathcal{R}[1]$-cluster tilting subcategories are precisely two-term silting subcategories in \cite{IJY}. Thus the bijection above induces the bijection given by Iyama-J{\o}rgensen-Yang in \cite{IJY}


title: Time Aggregation and Model Interpretation for Deep Multivariate Longitudinal Patient Outcome Forecasting Systems in Chronic Ambulatory Care

abstract: Clinical data for ambulatory care, which accounts for 90% of the nations healthcare spending, is characterized by relatively small sample sizes of longitudinal data, unequal spacing between visits for each patient, with unequal numbers of data points collected across patients. While deep learning has become state-of-the-art for sequence modeling, it is unknown which methods of time aggregation may be best suited for these challenging temporal use cases. Additionally, deep models are often considered uninterpretable by physicians which may prevent the clinical adoption, even of well performing models. We show that time-distributed-dense layers combined with GRUs produce the most generalizable models. Furthermore, we provide a framework for the clinical interpretation of the models.


title: Least Squares Rational Approximation

abstract: Rational approximation appears in many contexts throughout science and engineering, playing a central role in linear systems theory, special function approximation, and many others. There are many existing methods for solving the rational approximation problem, from fixed point methods like the Sanathanan-Koerner iteration and Vector Fitting, to partial interpolation methods like Adaptive Anderson Antoulas (AAA). While these methods can often find rational approximations with a small residual norm, they are unable to find optimizers with respect to a weighted l2 norm with a square dense weighting matrix. Here we develop a nonlinear least squares approach constructing rational approximations with respect to this norm. We explore this approach using two parameterizations of rational functions: a ratio of two polynomials and a partial fraction expansion. In both cases, we show how we can use Variable Projection (VARPRO) to reduce the dimension of the optimization problem. As many applications seek a real rational approximation that can be described as a ratio of two real polynomials, we show how this constraint can be enforced in both parameterizations. Although this nonlinear least squares approach often converge to suboptimal local minimizers, we find this can be largely mitigated by initializing the algorithm using the poles of the AAA algorithm applied to the same data. This combination of initialization and nonlinear least squares enables us to construct rational approximants using dense and potentially ill-conditioned weight matrices such as those that appear as a step in new H2 model reduction algorithm recently developed by the authors.


title: Active Learning in Recommendation Systems with Multi-level User Preferences

abstract: While recommendation systems generally observe user behavior passively, there has been an increased interest in directly querying users to learn their specific preferences. In such settings, considering queries at different levels of granularity to optimize user information acquisition is crucial to efficiently providing a good user experience. In this work, we study the active learning problem with multi-level user preferences within the collective matrix factorization (CMF) framework. CMF jointly captures multi-level user preferences with respect to items and relations between items (e.g., book genre, cuisine type), generally resulting in improved predictions. Motivated by finite-sample analysis of the CMF model, we propose a theoretically optimal active learning strategy based on the Fisher information matrix and use this to derive a realizable approximation algorithm for practical recommendations. Experiments are conducted using both the Yelp dataset directly and an illustrative synthetic dataset in the three settings of personalized active learning, cold-start recommendations, and noisy data -- demonstrating strong improvements over several widely used active learning methods.


title: Noise-resolution uncertainty principle in classical and quantum systems

abstract: It is proved that the width of a function and the width of the distribution of its values cannot be made arbitrarily small simultaneously. In the case of ergodic stochastic processes, an ensuing uncertainty relationship is demonstrated for the product of correlation length and variance. A closely related uncertainty principle is also established for the average degree of fourth-order coherence and the spatial width of modes of bosonic quantum fields. However, it is shown that, in the case of stochastic and quantum observables, certain non-classical states with sub-Poissonian statistics, such as for example photon number squeezed states in quantum optics, can overcome the "classical" noise-resolution uncertainty limit. The described uncertainty relationship, which is fundamentally different from the Heisenberg and related uncertainty principles, can define an upper limit for the information capacity of communication and imaging systems. It is expected to be useful in a variety of problems in classical and quantum physics.


title: Two-sample Test of Community Memberships of Weighted Stochastic Block Models

abstract: Suppose two networks are observed for the same set of nodes, where each network is assumed to be generated from a weighted stochastic block model. This paper considers the problem of testing whether the community memberships of the two networks are the same. A test statistic based on singular subspace distance is developed. Under the weighted stochastic block models with dense graphs, the limiting distribution of the proposed test statistic is developed. Simulation results show that the test has correct empirical type 1 errors under the dense graphs. The test also behaves as expected in empirical power, showing gradual changes when the intra-block and inter-block distributions are close and achieving 1 when the two distributions are not so close, where the closeness of the two distributions is characterized by Renyi divergence of order 1/2. The Enron email networks are used to demonstrate the proposed test.


title: Homogeneous Einstein metrics on Euclidean spaces are Einstein solvmanifolds

abstract: We show that homogeneous Einstein metrics on Euclidean spaces are Einstein solvmanifolds, using that they admit periodic, integrally minimal foliations by homogeneous hypersurfaces. For the geometric flow induced by the orbit-Einstein condition, we construct a Lyapunov function based on curvature estimates which come from real GIT.


title: Experimental Investigation of the Transition Energy $\gamma_t$ in the Isochronous Mode of the HIRFL-CSRe

abstract: The Isochronous Mass Spectrometry (IMS) based on storage rings is a powerful technique for mass measurement of short-lived exotic nuclei. The transition energy $\gamma_t$ of the storage ring is a vital parameter of the IMS technique. It is difficult to measure the $\gamma_t$ and its relation to momentum spread or circulating length, especially to monitor the variation of $\gamma_t$ during experiments. An experimental investigation on the $\gamma_t$ has been performed for the IMS experiment at the Cooler Storage Ring of the Heavy Ion Research Facility in Lanzhou (HIRFL-CSRe). With the velocity measured by two time-of-flight (TOF) detectors, the $\gamma_t$ as a function of orbital length can be determined. The influences of higher order magnetic field components on the $\gamma_t$ function were inferred for isochronous correction. This paper introduces and investigates the influence of dipole magnetic fields, quadrupole magnetic fields and sextupole magnetic fields on the $\gamma_t$ function. With the quadrupole magnets and sextupole magnets corrections, a mass resolution of 171332 (FWHM) and $\sigma(T)/T=1.34\times10^{-6}$ were reached, which shall be compared with 31319 (FWHM) and $\sigma(T)/T=7.35\times10^{-6}$ obtained without correction.


title: Parsing R-CNN for Instance-Level Human Analysis

abstract: Instance-level human analysis is common in real-life scenarios and has multiple manifestations, such as human part segmentation, dense pose estimation, human-object interactions, etc. Models need to distinguish different human instances in the image panel and learn rich features to represent the details of each instance. In this paper, we present an end-to-end pipeline for solving the instance-level human analysis, named Parsing R-CNN. It processes a set of human instances simultaneously through comprehensive considering the characteristics of region-based approach and the appearance of a human, thus allowing representing the details of instances. Parsing R-CNN is very flexible and efficient, which is applicable to many issues in human instance analysis. Our approach outperforms all state-of-the-art methods on CIHP (Crowd Instance-level Human Parsing), MHP v2.0 (Multi-Human Parsing) and DensePose-COCO datasets. Based on the proposed Parsing R-CNN, we reach the 1st place in the COCO 2018 Challenge DensePose Estimation task. Code and models are public available.


title: Constructing Trivariate B-splines with Positive Jacobian by Pillow Operation and Geometric Iterative Fitting

abstract: The advent of isogeometric analysis has prompted a need for methods to generate Trivariate B-spline Solids (TBS) with positive Jacobian. However, it is difficult to guarantee a positive Jacobian of a TBS since the geometric pre-condition for ensuring the positive Jacobian is very complicated. In this paper, we propose a method for generating TBSs with guaranteed positive Jacobian. For the study, we used a tetrahedral (tet) mesh model and segmented it into sub-volumes using the pillow operation. Then, to reduce the difficulty in ensuring a positive Jacobian, we separately fitted the boundary curves and surfaces and the sub-volumes using a geometric iterative fitting algorithm. Finally, the smoothness between adjacent TBSs is improved. The experimental examples presented in this paper demonstrate the effectiveness and efficiency of the developed algorithm.


title: X*: Anytime Multi-Agent Path Finding for Sparse Domains using Window-Based Iterative Repairs

abstract: Real-world multi-agent systems such as warehouse robots operate under significant time constraints -- in such settings, rather than spending significant amounts of time solving for optimal paths, it is instead preferable to find valid collision-free paths quickly, even if suboptimal, and given additional time, to iteratively refine such paths to improve their cost. In such domains, we observe that agent-agent collisions are sparse -- they involve small local subsets of agents, and are geographically contained within a small region of the overall space. Leveraging this insight, we can first plan paths for each agent individually, and in the cases of collisions between agents, perform small local repairs limited to local subspace windows. As time permits, these windows can be successively grown and the repairs within them refined, thereby improving the path quality, and eventually converging to the global joint optimal solution. Using these insights, we present two algorithmic contributions: 1) the Windowed Anytime Multiagent Planning Framework (WAMPF) for a class of anytime planners that quickly generate valid paths with suboptimality estimates and generate optimal paths given sufficient time, and 2) X*, an efficient WAMPF-based planner. X* is able to efficiently find successive valid solutions by employing re-use techniques during the repair growth step of WAMPF. Experimentally, we demonstrate that in sparse domains: 1) X* outperforms state-of-the-art anytime or optimal MAPF solvers in time to valid path, 2) X* is competitive with state-of-the-art anytime or optimal MAPF solvers in time to optimal path, 3) X* quickly converges to very tight suboptimality bounds, and 4) X* is competitive with state-of-the-art suboptimal MAPF solvers in time to valid path for small numbers of agents while providing much higher quality paths.


title: Gregory Solid Construction for Polyhedral Volume Parameterization by Sparse Optimization

abstract: In isogeometric analysis, it is frequently required to handle the geometric models enclosed by four-sided or non-four-sided boundary patches, such as trimmed surfaces. In this paper, we develop a Gregory solid based method to parameterize those models. First, we extend the Gregory patch representation to the trivariate Gregory solid representation. Second, the trivariate Gregory solid representation is employed to interpolate the boundary patches of a geometric model, thus generating the polyhedral volume parametrization. To improve the regularity of the polyhedral volume parametrization, we formulate the construction of the trivariate Gregory solid as a sparse optimization problem, where the optimization objective function is a linear combination of some terms, including a sparse term aiming to reduce the negative Jacobian area of the Gregory solid. Then, the alternating direction method of multipliers (ADMM) is used to solve the sparse optimization problem. Lots of experimental examples illustrated in this paper demonstrate the effectiveness and efficiency of the developed method.


title: Half-metallicity of Mn2VAl ferrimagnet revealed by resonant inelastic soft x-ray scattering under magnetic field

abstract: Detailed information on the electronic states of both V and Mn 3d electrons in the ferrimagnet Mn2VAl is obtained by the bulk sensitive resonant inelastic soft x-ray scattering (SX-RIXS) excited with the circularly polarized light under an external magnetic field for the first time. The results under the V L-edge excitation have revealed the negligible partial density of states (PDOS) of the V 3d states around the Fermi energy as well as their rather localized character. Under the Mn L-edge excitation, on the other hand, the spectra are dominated by fluorescence with clear magnetic circular dichroism with noticeable excitation photon energy dependence. Compared with the theoretical prediction of the RIXS spectra based on the density-functional-theory band structure calculation, an itinerant, spin-dependent character of the Mn 3d states and decays of the Mn 2p core states are confirmed in consistence with the half-metallicity of the Mn 3d states.


title: Adversarial Examples as an Input-Fault Tolerance Problem

abstract: We analyze the adversarial examples problem in terms of a model's fault tolerance with respect to its input. Whereas previous work focuses on arbitrarily strict threat models, i.e., $\epsilon$-perturbations, we consider arbitrary valid inputs and propose an information-based characteristic for evaluating tolerance to diverse input faults.


title: Local inversion-free estimation of spatial Gaussian processes

abstract: Maximizing the likelihood has been widely used for estimating the unknown covariance parameters of spatial Gaussian processes. However, evaluating and optimizing the likelihood function can be computationally intractable, particularly for large number of (possibly) irregularly spaced observations, due to the need to handle the inverse of ill-conditioned and large covariance matrices. Extending the "inversion-free" method of Anitescu, Chen and Stein \cite{anitescu2017inversion}, we investigate a broad class of covariance parameter estimators based on inversion-free surrogate losses and block diagonal approximation schemes of the covariance structure. This class of estimators yields a spectrum for negotiating the trade-off between statistical accuracy and computational cost. We present fixed-domain asymptotic properties of our proposed method, establishing $\sqrt{n}$-consistency and asymptotic normality results for isotropic Matern Gaussian processes observed on a multi-dimensional and irregular lattice. Simulation studies are also presented for assessing the scalability and statistical efficiency of the proposed algorithm for large data sets.


title: Hamel Spaces and Distal Expansions

abstract: In this note, we construct a distal expansion for the structure $(\mathbb{R}; +,<,H)$, where $H\subseteq \mathbb{R}$ is a dense $\mathbb{Q}$-vector space basis of $\mathbb{R}$ (a so-called Hamel basis). Our construction is also an expansion of the dense pair $(\mathbb{R}; +,<,\mathbb{Q})$ and has full quantifier elimination in a natural language.


title: Metric Based Quadrilateral Mesh Generation

abstract: This work proposes a novel metric based algorithm for quadrilateral mesh generating. Each quad-mesh induces a Riemannian metric satisfying special conditions: the metric is a flat metric with cone signualrites conformal to the original metric, the total curvature satisfies the Gauss-Bonnet condition, the holonomy group is a subgroup of the rotation group $\{e^{ik\pi/2}\}$, furthermore there is cross field obtained by parallel translation which is aligned with the boundaries, and its streamlines are finite geodesics. Inversely, such kind of metric induces a quad-mesh. Based on discrete Ricci flow and conformal structure deformation, one can obtain a metric satisfying all the conditions and obtain the desired quad-mesh. This method is rigorous, simple and automatic. Our experimental results demonstrate the efficiency and efficacy of the algorithm.


title: Minimizing Age-of-Information with Throughput Requirements in Multi-Path Network Communication

abstract: We consider the scenario where a sender periodically sends a batch of data to a receiver over a multi-hop network, possibly using multiple paths. Our objective is to minimize peak/average Age-of-Information (AoI) subject to throughput requirements. The consideration of batch generation and multi-path communication differentiates our AoI study from existing ones. We first show that our AoI minimization problems are NP-hard, but only in the weak sense, as we develop an optimal algorithm with a pseudo-polynomial time complexity. We then prove that minimizing AoI and minimizing maximum delay are "roughly" equivalent, in the sense that any optimal solution of the latter is an approximate solution of the former with bounded optimality loss. We leverage this understanding to design a general approximation framework for our problems. It can build upon any $\alpha$-approximation algorithm of the maximum delay minimization problem, to construct an $(\alpha+c)$-approximate solution for minimizing AoI. Here $c$ is a constant depending on the throughput requirements. Simulations over various network topologies validate the effectiveness of our approach.


title: Millimeter Wave Systems for Wireless Cellular Communications

abstract: This thesis considers channel estimation and multiuser (MU) data transmission for massive MIMO systems with fully digital/hybrid structures in mmWave channels. It contains three main contributions. In this thesis, we first propose a tone-based linear search algorithm to facilitate the estimation of angle-of-arrivals of the strongest components as well as scattering components of the users at the base station (BS) with fully digital structure. Our results show that the proposed maximum-ratio transmission (MRT) based on the strongest components can achieve a higher data rate than that of the conventional MRT, under the same mean squared errors (MSE). Second, we develop a low-complexity channel estimation and beamformer/precoder design scheme for hybrid mmWave systems. In addition, the proposed scheme applies to both non-sparse and sparse mmWave channel environments. We then leverage the proposed scheme to investigate the downlink achievable rate performance. The results show that the proposed scheme obtains a considerable achievable rate of fully digital systems. Taking into account the effect of various types of errors, we investigate the achievable rate performance degradation of the considered scheme. Third, we extend our proposed scheme to a multi-cell MU mmWave MIMO network. We derive the closed-form approximation of the normalized MSE of channel estimation under pilot contamination over Rician fading channels. Furthermore, we derive a tight closed-form approximation and the scaling law of the average achievable rate. Our results unveil that channel estimation errors, the intra-cell interference, and the inter-cell interference caused by pilot contamination over Rician fading channels can be efficiently mitigated by simply increasing the number of antennas equipped at the desired BS.


title: Learning Dynamics from Kinematics: Estimating 2D Foot Pressure Maps from Video Frames

abstract: Pose stability analysis is the key to understanding locomotion and control of body equilibrium, with applications in numerous fields such as kinesiology, medicine, and robotics. In biomechanics, Center of Pressure (CoP) is used in studies of human postural control and gait. We propose and validate a novel approach to learn CoP from pose of a human body to aid stability analysis. More specifically, we propose an end-to-end deep learning architecture to regress foot pressure heatmaps, and hence the CoP locations, from 2D human pose derived from video. We have collected a set of long (5min +) choreographed Taiji (Tai Chi) sequences of multiple subjects with synchronized foot pressure and video data. The derived human pose data and corresponding foot pressure maps are used jointly in training a convolutional neural network with residual architecture, named PressNET. Cross-subject validation results show promising performance of PressNET, significantly outperforming the baseline method of K-Nearest Neighbors. Furthermore, we demonstrate that our computation of center of pressure (CoP) from PressNET is not only significantly more accurate than those obtained from the baseline approach but also meets the expectations of corresponding lab-based measurements of stability studies in kinesiology.


title: DeepFlux for Skeletons in the Wild

abstract: Computing object skeletons in natural images is challenging, owing to large variations in object appearance and scale, and the complexity of handling background clutter. Many recent methods frame object skeleton detection as a binary pixel classification problem, which is similar in spirit to learning-based edge detection, as well as to semantic segmentation methods. In the present article, we depart from this strategy by training a CNN to predict a two-dimensional vector field, which maps each scene point to a candidate skeleton pixel, in the spirit of flux-based skeletonization algorithms. This "image context flux" representation has two major advantages over previous approaches. First, it explicitly encodes the relative position of skeletal pixels to semantically meaningful entities, such as the image points in their spatial context, and hence also the implied object boundaries. Second, since the skeleton detection context is a region-based vector field, it is better able to cope with object parts of large width. We evaluate the proposed method on three benchmark datasets for skeleton detection and two for symmetry detection, achieving consistently superior performance over state-of-the-art methods.


title: Invited Article: 4D Printing as a New Paradigm for Manufacturing with Minimum Energy Consumption

abstract: 4D printing is a new manufacturing paradigm that combines stimuli-responsive materials, mathematics, and multi-material additive manufacturing to yield encoded multi-material structures with intelligent behavior over time. This emerging field has received growing interests from various disciplines such as space exploration, renewable energy, bioengineering, textile industry, infrastructures, soft robotics, and so on. Here, as a first attempt, we consider the energy aspect of 4D printing. By a thermodynamic analysis, we obtain the theoretical limit of energy consumption in 4D printing and prove that 4D printing can be the most energy-efficient manufacturing process. Before that, we clearly underpin 4D printing as a new manufacturing process and identify its unique attributes.


title: An Energy-Efficient Transaction Model for the Blockchain-enabled Internet of Vehicles (IoV)

abstract: The blockchain is a safe, reliable and innovative mechanism for managing numerous vehicles seeking connectivity. However, following the principles of the blockchain, the number of transactions required to update ledgers pose serious issues for vehicles as these may consume the maximum available energy. To resolve this, an efficient model is presented in this letter which is capable of handling the energy demands of the blockchain-enabled Internet of Vehicles (IoV) by optimally controlling the number of transactions through distributed clustering. Numerical results suggest that the proposed approach is 40.16% better in terms of energy conservation and 82.06% better in terms of the number of transactions required to share the entire blockchain-data compared with the traditional blockchain.


title: Virtual Class Enhanced Discriminative Embedding Learning

abstract: Recently, learning discriminative features to improve the recognition performances gradually becomes the primary goal of deep learning, and numerous remarkable works have emerged. In this paper, we propose a novel yet extremely simple method \textbf{Virtual Softmax} to enhance the discriminative property of learned features by injecting a dynamic virtual negative class into the original softmax. Injecting virtual class aims to enlarge inter-class margin and compress intra-class distribution by strengthening the decision boundary constraint. Although it seems weird to optimize with this additional virtual class, we show that our method derives from an intuitive and clear motivation, and it indeed encourages the features to be more compact and separable. This paper empirically and experimentally demonstrates the superiority of Virtual Softmax, improving the performances on a variety of object classification and face verification tasks.


title: An analytical approach to the field amplification and particle production by parametric resonance during inflation and reheating

abstract: Field amplification and particle production due to parametric resonance are highly nontrivial predictions of quantum fields that couple to an oscillating source during inflation and reheating. Understanding this two effects is crucial for the connection between the resonance phenomenon and precise observational data. In this paper, we give a general and analytic analysis of parametric resonance of relevant field modes evolving during inflation and reheating by using the uniform asymptotic approximation. This analysis can provide a clear and quantitative explanation for the field amplification and particle production during the resonance. The potential applications of our results to several examples, including sound resonance during inflation, particle productions during reheating, and parametric resonance due to self-resonance potentials, have also been explored. The formalism developed in this paper is also applicable to parametric resonance in a broad areas of modern science.


title: Selective transport of atomic excitations in a driven chiral-coupled atomic chain

abstract: We theoretically investigate the flow of the atomic excitations in a driven chiral-coupled atomic chain with nonreciprocal decay channels. This one-dimensional system allows infinite-range dipole-dipole interaction, and enables directional guided modes of scattered light. Under a weakly driven condition, we are able to simulate the transport properties of atomic excitations between the left and right parts of the chain. In the steady states, the transport is highly dependent on the equidistant positions of the ordered array, the excitation field detunings, and the directionality of such chiral-coupled system. We discuss the parameter regimes which are resilient or sensitive to position fluctuations, and provide insights to precise and genuine state preparations. Furthermore, we study the effect of position fluctuations on the transport of excitations. Our results can shed light on deterministic state preparations and pave the way for many-body state manipulations in the driven and dissipative chiral-coupled systems.


title: X-ray Diffraction Analysis of Cu2+ Doped Zn1-xCuxFe2O4 Spinel Nanoparticles using Williamson-Hall Plot Method

abstract: The nanoparticles (NPs) of Zn1-xCuxFe2O4 (ZCFO) spinels with x = 0, 0.2, 0.4, 0.6 and 0.8 were synthesized by a sol-gel combustion method using acetate precursor. The NPs of ZCFO were prepared by following calcination process at 600C for 8hrs. The synthesized NPs of ZCFO were characterized by X-ray diffraction (XRD) analysis using Rietveld refinement. The Rietveld refinement of the XRD patterns revealed that the ZCFO spinels crystallize into single diamond cubic structure with Fd-3m space group. The lattice constant and unit cell volume for ZCFO NPs shrink with enhancing doping concentration of Cu2+ ion. The crystalline growth in the NPs of ZCFO was examined by peak broadening present in the XRD pattern. The Williamson-Hall (W-H) plot method were used to study the individual involvements of crystallite sizes and lattice strain on the peak broadening of the NPs of ZCFO spinels. Whereas, particle size of the ZCFO sample with x = 0.40 was estimated by high-resolution scanning electron microscopy micrographs


title: An Interpretable Model with Globally Consistent Explanations for Credit Risk

abstract: We propose a possible solution to a public challenge posed by the Fair Isaac Corporation (FICO), which is to provide an explainable model for credit risk assessment. Rather than present a black box model and explain it afterwards, we provide a globally interpretable model that is as accurate as other neural networks. Our "two-layer additive risk model" is decomposable into subscales, where each node in the second layer represents a meaningful subscale, and all of the nonlinearities are transparent. We provide three types of explanations that are simpler than, but consistent with, the global model. One of these explanation methods involves solving a minimum set cover problem to find high-support globally-consistent explanations. We present a new online visualization tool to allow users to explore the global model and its explanations.


title: Probability of radiation of twisted photons by axially symmetric bunches of particles

abstract: The effect of a finite width of a particle bunch on radiation of twisted photons is studied. The general formulas connecting the radiation probability distribution of twisted photons produced by bunches of identical particles with the radiation probability distribution of twisted photons generated by one particle are obtained for axially symmetric bunches. The bunch is called axially symmetric if it is axially symmetric with respect to the detector axis at some instant of time and all the particles in the bunch move along parallel trajectories. The general sum rules for the probability of radiation of twisted photons by axially symmetric bunches are established. In particular, we prove that the projection of the average total angular momentum of radiated twisted photons per particle in the bunch does not depend on the radial profile of the bunch. The uniform, Gaussian, and exponential radial bunch profiles are considered in detail. The radiation of axially symmetric bunches in ordinary and crystalline undulators is investigated. The selection rules for radiation of twisted photons by one particle in undulators are violated when the finite width of the particle bunch is taken into account. We find the condition when this violation is marginal. The form of the radiation probability distribution of twisted photons becomes universal for wide incoherent axially symmetric particle bunches. We completely describe these universal distributions.


title: Imaging black holes through AdS/CFT

abstract: Thermal states in some quantum field theories (QFTs) correspond to black holes in asymptotically AdS spacetime in the AdS/CFT correspondence. We propose a direct procedure to construct holographic images of the black hole in the bulk from a given response function of the QFT on the boundary. The response function with respect to an external source corresponds to the asymptotic data of the bulk field generated by the source on the AdS boundary. According to the wave optics, we can obtain the images from the bulk field propagating in the bulk spacetime. For a thermal state on two-dimensional sphere dual to Schwarzschild-AdS$_4$ black hole, we demonstrate that the holographic images gravitationally lensed by the black hole can be constructed from the response function. In particular, the Einstein rings on the image can be clearly observed and their radius depends on the total energy of the QFT thermal state. These results are consistent with the size of the photon sphere of the black hole calculated in geometrical optics. This implies that, if there exists a dual gravitational picture for a given quantum system, we would be able to probe existence of the dual black hole by the Einstein rings constructed from observables of the quantum system.

categories: hep-th astro-ph.HE cond-mat.str-el gr-qc


title: Structure and Dynamics of the Quasi-Liquid Layer at the Surface of Ice from Molecular Simulations

abstract: We characterized the structural and dynamical properties of the quasi-liquid layer (QLL) at the surface of ice by molecular dynamics simulations with a thermodynamically consistent water model. Our simulations show that for three low-index ice surfaces only the outermost molecular layer presents short-range and mid-range disorder and is diffusive. The onset temperature for normal diffusion is much higher than the glass temperature of supercooled water, although the diffusivity of the QLL is higher than that of bulk water at the corresponding temperature. The underlying subsurface layers impose an ordered template, which produces a regular patterning of the ice/water interface at any temperature, and is responsible for the major differences between QLL and bulk water, especially for what concern the dynamics and the mid-range structure of the hydrogen-bonded network. Our work highlights the need of a holistic approach to the characterization of QLL, as a single experimental technique may probe only one specific feature, missing part of the complexity of this fascinating system.


title: Dirichlet and Neumann problems for elliptic equations with singular drifts on Lipschitz domains

abstract: We consider the Dirichlet and Neumann problems for second-order linear elliptic equations: \[ -\triangle u +\mathrm{div}(u\mathbf{b}) =f \quad\text{ and }\quad -\triangle v -\mathbf{b} \cdot \nabla v =g \] in a bounded Lipschitz domain $\Omega$ in $\mathbb{R}^n$ $(n\geq 3)$, where $\mathbf{b}:\Omega \rightarrow \mathbb{R}^n$ is a given vector field. Under the assumption that $\mathbf{b} \in L^{n}(\Omega)^n$, we first establish existence and uniqueness of solutions in $L_{\alpha}^{p}(\Omega)$ for the Dirichlet and Neumann problems. Here $L_{\alpha}^{p}(\Omega)$ denotes the Sobolev space (or Bessel potential space) with the pair $(\alpha,p)$ satisfying certain conditions. These results extend the classical works of Jerison-Kenig [17] and Fabes-Mendez-Mitrea [12] for the Poisson equation. We also prove existence and uniqueness of solutions of the Dirichlet problem with boundary data in $L^{2}(\partial\Omega)$. Our results for the Dirichlet problems hold even for the case $n=2$.


title: Change Point Models for Real-time V2I Cyber Attack Detection in a Connected Vehicle Environment

abstract: Connected vehicle (CV) systems are cognizant of potential cyber attacks because of increasing connectivity between its different components such as vehicles, roadside infrastructure and traffic management centers. However, it is a challenge to detect security threats in real-time and develop appropriate/effective countermeasures for a CV system because of the dynamic behavior of such attacks, high computational power requirement and a historical data requirement for training detection models. To address these challenges, statistical models, especially change point models, have potentials for real-time anomaly detections. Thus, the objective of this study is to investigate the efficacy of two change point models, Expectation Maximization (EM) and Cumulative Sum (CUSUM), for real-time V2I cyber attack detection in a CV Environment. To prove the efficacy of these models, we evaluated these two models for three different type of cyber attack, denial of service (DOS), impersonation, and false information, using basic safety messages (BSMs) generated from CVs through simulation. Results from numerical analysis revealed that EM and CUSUM could detect these cyber attacks, DOS, impersonation, and false information, with an accuracy of 99\%, 100\%, and 98\%, and 100\%, 100\% and 98\%, respectively.


title: A Core Ontology for Privacy Requirements Engineering

abstract: Nowadays, most companies need to collect, store, and manage personal information in order to deliver their services. Accordingly, privacy has emerged as a key concern for these companies since they need to comply with privacy laws and regulations. To deal with them properly, such privacy concerns should be considered since the early phases of system design. Ontologies have proven to be a key factor for elaborating high-quality requirements models. However, most existing work deals with privacy as a special case of security requirements, thereby missing essential traits of this family of requirements. In this paper, we introduce COPri, a Core Ontology for Privacy requirements engineering that adopts and extends our previous work on privacy requirements engineering ontology that has been mined through a systematic literature review. Additionally, we implement, validate and then evaluate our ontology.


title: Thermal Transport Across Graphene Step Junctions

abstract: Step junctions are often present in layered materials, i.e. where single-layer regions meet multi-layer regions, yet their effect on thermal transport is not understood to date. Here, we measure heat flow across graphene junctions (GJs) from monolayer to bilayer graphene, as well as bilayer to four-layer graphene for the first time, in both heat flow directions. The thermal conductance of the monolayer-bilayer GJ device ranges from ~0.5 to 9.1x10^8 Wm-2K-1 between 50 K to 300 K. Atomistic simulations of such GJ device reveal that graphene layers are relatively decoupled, and the low thermal conductance of the device is determined by the resistance between the two dis-tinct graphene layers. In these conditions the junction plays a negligible effect. To prove that the decoupling between layers controls thermal transport in the junction, the heat flow in both directions was measured, showing no evidence of thermal asymmetry or rectification (within experimental error bars). For large-area graphene applications, this signifies that small bilayer (or multilayer) islands have little or no contribution to overall thermal transport.


title: Reservoir Computing with Random Skyrmion Textures

abstract: The Reservoir Computing (RC) paradigm posits that sufficiently complex physical systems can be used to massively simplify pattern recognition tasks and nonlinear signal prediction. This work demonstrates how random topological magnetic textures present sufficiently complex resistance responses for the implementation of RC as applied to A/C current pulses. In doing so, we stress how the applicability of this paradigm hinges on very general dynamical properties which are satisfied by a large class of physical systems where complexity can be put to computational use. By harnessing the complex resistance response exhibited by random magnetic skyrmion textures and using it to demonstrate pattern recognition, we explain how spintronically accessible magnetic systems offer an advantage in the search for an ideal reservoir computer. The low-power properties of compact skyrmion fabrics, coupled with their CMOS integrability operating on similar length and timescales, open the door for their RC employment geared towards industrial application.


title: Modality-based Factorization for Multimodal Fusion

abstract: We propose a novel method, Modality-based Redundancy Reduction Fusion (MRRF), for understanding and modulating the relative contribution of each modality in multimodal inference tasks. This is achieved by obtaining an $(M+1)$-way tensor to consider the high-order relationships between $M$ modalities and the output layer of a neural network model. Applying a modality-based tensor factorization method, which adopts different factors for different modalities, results in removing information present in a modality that can be compensated by other modalities, with respect to model outputs. This helps to understand the relative utility of information in each modality. In addition it leads to a less complicated model with less parameters and therefore could be applied as a regularizer avoiding overfitting. We have applied this method to three different multimodal datasets in sentiment analysis, personality trait recognition, and emotion recognition. We are able to recognize relationships and relative importance of different modalities in these tasks and achieves a 1\% to 4\% improvement on several evaluation measures compared to the state-of-the-art for all three tasks.


title: Performance analysis of Schottky-junction-based near-field thermophotovoltaic system

abstract: Numerous studies have reported performance enhancement of a thermophotovoltaic (TPV) system when an emitter is separated by nanoscale gaps from a TPV cell. Although a $p$-$n$-junction-based TPV cell has been widely used for the near-field TPV system, a Schottky-junction-based near-field TPV system has drawn attention recently with the advantage of the easy fabrication. However, existing studies mostly focused on the generated photocurrent only in the metal side due to the fact that required energy for the metal-side photocurrent (i.e., Schottky barrier height) is smaller than the bandgap energy. Here, we suggest the precise performance analysis model for the Schottky-junction-based near-field TPV system, including photocurrent generation on the semiconductor side by considering the transport of minority carriers within the semiconductor. It is found that most of the total photocurrent in the Schottky-junction-based near-field TPV system is generated in the semiconductor side. We also demonstrate that further enhancement in the photocurrent generation can be achieved by re-absorbing the usable photon energy in the metal with the help of a backside reflector. The present work will provide a design guideline for the Schottky-junction-based near-field TPV system taking into account three types of photocurrents.


title: Dynamical precession of spin in the two-dimensional spin-orbit coupled systems

abstract: We investigate the spin dynamics in the two-dimensional spin-orbit coupled system subject to an in-plane ($x$-$y$ plane) constant electric field, which is assumed to be turned on at the moment $t=0$. The equation of spin precession in linear response to the switch-on of the electric field is derived in terms of Heisenberg's equation by the perturbation method up to the first order of the electric field. The dissipative effect, which is responsible for bringing the dynamical response to an asymptotic result, is phenomenologically implemented \`{a} la the Landau-Lifshitz-Gilbert equation by introducing damping terms upon the equation of spin dynamics. Mediated by the dissipative effect, the resulting spin dynamics asymptotes to a stationary state, where the spin and the momentum-dependent effective magnetic field are aligned again and have nonzero components in the out-of-plane ($z$) direction. In the linear response regime, the asymptotic response obtained by the dynamical treatment is in full agreement with the stationary response as calculated in the Kubo formula, which is a time-independent approach treating the applied electric field as completely time-independent. Our method provides a new perspective on the connection between the dynamical and stationary responses.