title: Uniformly rotating, axisymmetric and triaxial quark stars in general relativity

abstract: Quasi-equilibrium models of uniformly rotating axisymmetric and triaxial quark stars are computed in general relativistic gravity scenario. The Isenberg-Wilson-Mathews (IWM) formulation is employed and the Compact Object CALculator (COCAL) code is extended to treat rotating stars with finite surface density and new equations of state (EOSs). Besides the MIT bag model for quark matter which is composed of de-confined quarks, we examine a new EOS proposed by Lai and Xu that is based on quark clustering and results in a stiff EOS that can support masses up to $3.3M_\odot$ in the case we considered. We perform convergence tests for our new code to evaluate the effect of finite surface density in the accuracy of our solutions and construct sequences of solutions for both small and high compactness. The onset of secular instability due to viscous dissipation is identified and possible implications are discussed. An estimate of the gravitational wave amplitude and luminosity based on quadrupole formulas is presented and comparison with neutron stars is discussed.


title: PoseCNN: A Convolutional Neural Network for 6D Object Pose Estimation in Cluttered Scenes

abstract: Estimating the 6D pose of known objects is important for robots to interact with the real world. The problem is challenging due to the variety of objects as well as the complexity of a scene caused by clutter and occlusions between objects. In this work, we introduce PoseCNN, a new Convolutional Neural Network for 6D object pose estimation. PoseCNN estimates the 3D translation of an object by localizing its center in the image and predicting its distance from the camera. The 3D rotation of the object is estimated by regressing to a quaternion representation. We also introduce a novel loss function that enables PoseCNN to handle symmetric objects. In addition, we contribute a large scale video dataset for 6D object pose estimation named the YCB-Video dataset. Our dataset provides accurate 6D poses of 21 objects from the YCB dataset observed in 92 videos with 133,827 frames. We conduct extensive experiments on our YCB-Video dataset and the OccludedLINEMOD dataset to show that PoseCNN is highly robust to occlusions, can handle symmetric objects, and provide accurate pose estimation using only color images as input. When using depth data to further refine the poses, our approach achieves state-of-the-art results on the challenging OccludedLINEMOD dataset. Our code and dataset are available at https://rse-lab.cs.washington.edu/projects/posecnn/.


title: A singular local minimizer for the volume constrained minimal surface problem in a nonconvex domain

abstract: It has recently been established byWang and Xia [WX] that local minimizers of perimeter within a ball subject to a volume constraint must be spherical caps or planes through the origin. This verifies a conjecture of the authors and is in contrast to the situation of area minimizing surfaces with prescribed boundary where singularities can be present in high dimensions. This result lends support to the more general conjecture that volume-constrained minimizers in arbitrary convex sets may enjoy better regularity properties than their boundary-prescribed cousins. Here, we show the importance of the convexity condition by exhibiting a simple example, given by the Simons cone, of a singular volume-constrained locally area-minimizing surface within a nonconvex domain that is arbitrarily close to the unit ball.


title: Credimus

abstract: We believe that economic design and computational complexity---while already important to each other---should become even more important to each other with each passing year. But for that to happen, experts in on the one hand such areas as social choice, economics, and political science and on the other hand computational complexity will have to better understand each other's worldviews. This article, written by two complexity theorists who also work in computational social choice theory, focuses on one direction of that process by presenting a brief overview of how most computational complexity theorists view the world. Although our immediate motivation is to make the lens through which complexity theorists see the world be better understood by those in the social sciences, we also feel that even within computer science it is very important for nontheoreticians to understand how theoreticians think, just as it is equally important within computer science for theoreticians to understand how nontheoreticians think.


title: Interplay between morphological and shielding effects in field emission via Schwarz-Christoffel transformation

abstract: It is well known that sufficiently strong electrostatic fields are able to change the morphology of Large Area Field Emitters (LAFEs). This phenomenon affects the electrostatic interactions between adjacent sites on a LAFE during field emission and may lead to several consequences, such as: the emitter's degradation, diffusion of absorbed particles on the emitter's surface, deflection due to electrostatic forces and mechanical stress. These consequences are undesirable for technological applications, since they may significantly affect the macroscopic current density on the LAFE. Despite the technological importance, these processes are not completely understood yet. Moreover, the electrostatic effects due to the proximity between emitters on a LAFE may compete with the morphological ones. The balance between these effects may lead to a non trivial behavior in the apex-Field Enhancement Factor (FEF). The present work intends to study the interplay between proximity and morphological effects by studying a model amenable for an analytical treatment. In order to do that, a conducting system under an external electrostatic field, with a profile limited by two mirror-reflected triangular protrusions on an infinite line, is considered. The FEF near the apex of each emitter is obtained as a function of their shape and the distance between them via a Schwarz-Christoffel transformation. Our results suggest that a tradeoff between morphological and proximity effects on a LAFE may provide an explanation for the observed reduction of the local FEF and its variation at small distances between the emitter sites.


title: A unified method for maximal truncated Calder\'on-Zygmund operators in general function spaces by sparse domination

abstract: In this note we give simple proofs of several results involving maximal truncated Calde\'on-Zygmund operators in the general setting of rearrangement invariant quasi-Banach function spaces by sparse domination. Our techniques allow us to track the dependence of the constants in weighted norm inequalities; additionally, our results hold in $\mathbb{R}^n$ as well as in many spaces of homogeneous type.


title: Negative differential thermal conductance and heat amplification in a nonequilibrium triangle-coupled spin-boson system at strong coupling

abstract: We investigate the nonequilibrium quantum heat transfer in a triangle-coupled spin-boson system within a three-terminal setup. By including the nonequilibrium noninteracting blip approximation approach combined with the full counting statistics, we analytically obtain the steady state populations and heat currents. The negative differential thermal conductance and giant heat amplification factor are clearly observed at strong qubit-bath coupling. %and the heat amplification is dramatically suppressed in the moderate coupling regime. Moreover, the strong interaction between the gating qubit and gating thermal bath is unraveled to be compulsory to exhibit these far-from equilibrium features.


title: Towards Effective Low-bitwidth Convolutional Neural Networks

abstract: This paper tackles the problem of training a deep convolutional neural network with both low-precision weights and low-bitwidth activations. Optimizing a low-precision network is very challenging since the training process can easily get trapped in a poor local minima, which results in substantial accuracy loss. To mitigate this problem, we propose three simple-yet-effective approaches to improve the network training. First, we propose to use a two-stage optimization strategy to progressively find good local minima. Specifically, we propose to first optimize a net with quantized weights and then quantized activations. This is in contrast to the traditional methods which optimize them simultaneously. Second, following a similar spirit of the first method, we propose another progressive optimization approach which progressively decreases the bit-width from high-precision to low-precision during the course of training. Third, we adopt a novel learning scheme to jointly train a full-precision model alongside the low-precision one. By doing so, the full-precision model provides hints to guide the low-precision model training. Extensive experiments on various datasets ( i.e., CIFAR-100 and ImageNet) show the effectiveness of the proposed methods. To highlight, using our methods to train a 4-bit precision network leads to no performance decrease in comparison with its full-precision counterpart with standard network architectures ( i.e., AlexNet and ResNet-50).


title: Steady state current fluctuations and dynamical control in a nonequilibrium single-site Bose-Hubbard system

abstract: We investigate nonequilibrium energy transfer in a single-site Bose-Hubbard model coupled to two thermal baths. By including a quantum kinetic equation combined with full counting statistics, we investigate the steady state energy flux and noise power. The influence of the nonlinear Bose-Hubbard interaction on the transfer behaviors is analyzed, and the nonmonotonic features are clearly exhibited. Particularly, in the strong on-site repulsion limit, the results become identical with the nonequilibrium spin-boson model. We also extend the quantum kinetic equation to study the geometric-phase-induced energy pump. An interesting reversal behavior is unraveled by enhancing the Bose-Hubbard repulsion strength.


title: Learning deep features for source color laser printer identification based on cascaded learning

abstract: Color laser printers have fast printing speed and high resolution, and forgeries using color laser printers can cause significant harm to society. A source printer identification technique can be employed as a countermeasure to those forgeries. This paper presents a color laser printer identification method based on cascaded learning of deep neural networks. The refiner network is trained by adversarial training to refine the synthetic dataset for halftone color decomposition. The halftone color decomposing ConvNet is trained with the refined dataset, and the trained knowledge is transferred to the printer identifying ConvNet to enhance the accuracy. The robustness about rotation and scaling is considered in training process, which is not considered in existing methods. Experiments are performed on eight color laser printers, and the performance is compared with several existing methods. The experimental results clearly show that the proposed method outperforms existing source color laser printer identification methods.


title: Privacy Preserving and Collusion Resistant Energy Sharing

abstract: Energy has been increasingly generated or collected by different entities on the power grid (e.g., universities, hospitals and householdes) via solar panels, wind turbines or local generators in the past decade. With local energy, such electricity consumers can be considered as "microgrids" which can simulataneously generate and consume energy. Some microgrids may have excessive energy that can be shared to other power consumers on the grid. To this end, all the entities have to share their local private information (e.g., their local demand, local supply and power quality data) to each other or a third-party to find and implement the optimal energy sharing solution. However, such process is constrained by privacy concerns raised by the microgrids. In this paper, we propose a privacy preserving scheme for all the microgrids which can securely implement their energy sharing against both semi-honest and colluding adversaries. The proposed approach includes two secure communication protocols that can ensure quantified privacy leakage and handle collusions.


title: Coherence and entanglement dynamics of vibrating qubits

abstract: We investigate the dynamics of coherence and entanglement of vibrating qubits. Firstly, we consider a single trapped ion qubit inside a perfect cavity and successively we use it to construct a bipartite system made of two of such subsystems, taken identical and noninteracting. As a general result, we find that qubit vibration can lead to prolonging initial coherence in both single-qubit and two-qubit system. However, despite of this coherence preservation, we show that the decay of the entanglement between the two qubits is sped up by the vibrational motion of the qubits. Furthermore, we highlight how the dynamics of photon-phonon correlations between cavity mode and vibrational mode, which may serve as a further useful resource stored in the single-qubit system, is strongly affected by the initial state of the qubit. These results provide new insights about the ability of systems made of moving qubits in maintaining quantum resources compared to systems of stationary qubits.


title: On the complete weight enumerators of some linear codes with a few weights

abstract: Linear codes with a few weights have important applications in authentication codes, secret sharing, consumer electronics, etc.. The determination of the parameters such as Hamming weight distributions and complete weight enumerators of linear codes are important research topics. In this paper, we consider some classes of linear codes with a few weights and determine the complete weight enumerators from which the corresponding Hamming weight distributions are derived with help of some sums involving Legendre symbol.


title: Stability of the simplex bound for packings by equal spherical caps determined by simplicial regular polytopes

abstract: It is well known that the vertices of any Euclidean simplicial regular polytope determine an optimal packing of equal spherical balls. We prove a stability version of optimal order of this result.


title: Large induced acyclic and outerplanar subgraphs of 2-outerplanar graph

abstract: Albertson and Berman conjectured that every planar graph has an induced forest on half of its vertices. The best known lower bound, due to Borodin, is that every planar graph has an induced forest on two fifths of its vertices. In a related result, Chartran and Kronk, proved that the vertices of every planar graph can be partitioned into three sets, each of which induce a forest. We show tighter results for 2-outerplanar graphs. We show that every 2-outerplanar graph has an induced forest on at least half the vertices by showing that its vertices can be partitioned into two sets, each of which induces a forest. We also show that every 2-outerplanar graph has an induced outerplanar graph on at least two-thirds of its vertices.


title: Closed Form Solutions of Combinatorial Graph Laplacian Estimation under Acyclic Topology Constraints

abstract: How to obtain a graph from data samples is an important problem in graph signal processing. One way to formulate this graph learning problem is based on Gaussian maximum likelihood estimation, possibly under particular topology constraints. To solve this problem, we typically require iterative convex optimization solvers. In this paper, we show that when the target graph topology does not contain any cycle, then the solution has a closed form in terms of the empirical covariance matrix. This enables us to efficiently construct a tree graph from data, even if there is only a single data sample available. We also provide an error bound of the objective function when we use the same solution to approximate a cyclic graph. As an example, we consider an image denoising problem, in which for each input image we construct a graph based on the theoretical result. We then apply low-pass graph filters based on this graph. Experimental results show that the weights given by the graph learning solution lead to better denoising results than the bilateral weights under some conditions.


title: A Coalition Formation Approach to Coordinated Task Allocation in Heterogeneous UAV Networks

abstract: The problem of adversary target detection and the subsequent task completion using a heterogeneous network of resource-constrained UAVs is considered. No prior knowledge about locations and required resources to identify these targets is available to the UAVs. In the proposed leader-follower coalition formation model, the UAV that first locates a target serves as the coalition leader and selects a group of follower UAVs to complete the task associated with the identified target. The goal of the coalition formation is to complete the designated tasks with minimal resource utilization. Another role of coalition members is to make the ground station aware of the detected adversary target by forwarding its signal to the station via a distributed cooperative relaying scheme. We also propose a reputation-based mechanism for coalition formation to monitor the cooperative behavior of the UAVs over the course of time and exclude potentially untrustworthy UAVs. Simulation results show the efficiency of the proposed method in forming optimal coalitions compared to alternative methods.


title: Minimum Energy Quantized Neural Networks

abstract: This work targets the automated minimum-energy optimization of Quantized Neural Networks (QNNs) - networks using low precision weights and activations. These networks are trained from scratch at an arbitrary fixed point precision. At iso-accuracy, QNNs using fewer bits require deeper and wider network architectures than networks using higher precision operators, while they require less complex arithmetic and less bits per weights. This fundamental trade-off is analyzed and quantified to find the minimum energy QNN for any benchmark and hence optimize energy-efficiency. To this end, the energy consumption of inference is modeled for a generic hardware platform. This allows drawing several conclusions across different benchmarks. First, energy consumption varies orders of magnitude at iso-accuracy depending on the number of bits used in the QNN. Second, in a typical system, BinaryNets or int4 implementations lead to the minimum energy solution, outperforming int8 networks up to 2-10x at iso-accuracy. All code used for QNN training is available from https://github.com/BertMoons.


title: Observation of Hopping and Blockade of Bosons in a Trapped Ion Spin Chain

abstract: The local phonon modes in a Coulomb crystal of trapped ions can represent a Hubbard system of coupled bosons. We selectively prepare single excitations at each site and observe free hopping of a boson between sites, mediated by the long-range Coulomb interaction between ions. We then implement phonon blockades on targeted sites by driving a Jaynes-Cummings interaction on individually addressed ions to couple their internal spin to the local phonon mode. The resulting dressed states have energy splittings that can be tuned to suppress phonon hopping into the site. This new experimental approach opens up the possibility of realizing large-scale Hubbard systems from the bottom up with tunable interactions at the single-site level.


title: Limiting Laws for Divergent Spiked Eigenvalues and Largest Non-spiked Eigenvalue of Sample Covariance Matrices

abstract: We study the asymptotic distributions of the spiked eigenvalues and the largest nonspiked eigenvalue of the sample covariance matrix under a general covariance matrix model with divergent spiked eigenvalues, while the other eigenvalues are bounded but otherwise arbitrary. The limiting normal distribution for the spiked sample eigenvalues is established. It has distinct features that the asymptotic mean relies on not only the population spikes but also the nonspikes and that the asymptotic variance in general depends on the population eigenvectors. In addition, the limiting Tracy-Widom law for the largest nonspiked sample eigenvalue is obtained. Estimation of the number of spikes and the convergence of the leading eigenvectors are also considered. The results hold even when the number of the spikes diverges. As a key technical tool, we develop a Central Limit Theorem for a type of random quadratic forms where the random vectors and random matrices involved are dependent. This result can be of independent interest.


title: Spatio-Temporal Reference Frames as Geographic Objects

abstract: It is often desirable to analyse trajectory data in local coordinates relative to a reference location. Similarly, temporal data also needs to be transformed to be relative to an event. Together, temporal and spatial contextualisation permits comparative analysis of similar trajectories taken across multiple reference locations. To the GIS professional, the procedures to establish a reference frame at a location and reproject the data into local coordinates are well known, albeit tedious. However, GIS tools are now often used by subject matter experts who may not have the deep knowledge of coordinate frames and projections required to use these techniques effectively. We introduce a novel method for representing spatio-temporal reference frames using ordinary geographic objects available in GIS tools. We argue that our method both reduces the number of manual steps required to reproject data to a local reference frame, in addition to reducing the number of concepts a novice user would need to learn.


title: Cumulants, Spreadability and the Campbell-Baker-Hausdorff Series

abstract: We define spreadability systems as a generalization of exchangeability systems in order to unify various notions of independence and cumulants known in noncommutative probability. In particular, our theory covers monotone independence and monotone cumulants which do not satisfy exchangeability. To this end we study generalized zeta and M\"obius functions in the context of the incidence algebra of the semilattice of ordered set partitions and prove an appropriate variant of Faa di Bruno's theorem. With the aid of this machinery we show that our cumulants cover most of the previously known cumulants. Due to noncommutativity of independence the behaviour of these cumulants with respect to independent random variables is more complicated than in the exchangeable case and the appearance of Goldberg coefficients exhibits the role of the Campbell-Baker-Hausdorff series in this context. In a final section we exhibit an interpretation of the Campbell-Baker-Hausdorff series as a sum of cumulants in a particular spreadability system, thus providing a new derivation of the Goldberg coefficients.


title: The Hardness of Synthesizing Elementary Net Systems from Highly Restricted Inputs

abstract: Elementary net systems (ENS) are the most fundamental class of Petri nets. Their synthesis problem has important applications in the design of digital hardware and commercial processes. Given a labeled transition system (TS) $A$, feasibility is the NP-complete decision problem whether $A$ can be equivalently synthesized into an ENS. It is well known that $A$ is feasible if and only if it has the event state separation property (ESSP) and the state separation property (SSP). Recently, these properties have also been studied individually for their practical implications. A fast ESSP algorithm, for instance, would allow applications to at least validate the language equivalence of $A$ and a synthesized ENS. Being able to efficiently decide SSP, on the other hand, could serve as a quick-fail preprocessing mechanism for synthesis. Although a few tractable subclasses have been found, this paper destroys much of the hope that many practically meaningful input restrictions make feasibility or at least one of ESSP and SSP efficient. We show that all three problems remain NP-complete even if the input is restricted to linear TSs where every event occurs at most three times or if the input is restricted to TSs where each event occurs at most twice and each state has at most two successor and two predecessor states.


title: Stochastic Variational Inference for Bayesian Sparse Gaussian Process Regression

abstract: This paper presents a novel variational inference framework for deriving a family of Bayesian sparse Gaussian process regression (SGPR) models whose approximations are variationally optimal with respect to the full-rank GPR model enriched with various corresponding correlation structures of the observation noises. Our variational Bayesian SGPR (VBSGPR) models jointly treat both the distributions of the inducing variables and hyperparameters as variational parameters, which enables the decomposability of the variational lower bound that in turn can be exploited for stochastic optimization. Such a stochastic optimization involves iteratively following the stochastic gradient of the variational lower bound to improve its estimates of the optimal variational distributions of the inducing variables and hyperparameters (and hence the predictive distribution) of our VBSGPR models and is guaranteed to achieve asymptotic convergence to them. We show that the stochastic gradient is an unbiased estimator of the exact gradient and can be computed in constant time per iteration, hence achieving scalability to big data. We empirically evaluate the performance of our proposed framework on two real-world, massive datasets.


title: Latitude Distribution of Sunspots: Analysis Using Sunspot Data and A Dynamo Model

abstract: In this paper, we explore the evolution of sunspot latitude distribution and explore its relations with the cycle strength. With the progress of the solar cycle, the distributions in two hemispheres from mid-latitudes propagate toward the equator and then (before the usual solar minimum) these two distributions touch each other. By visualizing the evolution of the distributions in two hemispheres, we separate the solar cycles by excluding this hemispheric overlap. From these isolated solar cycles in two hemispheres, we generate latitude distributions for each cycle, starting from cycle 8 to cycle 23. We find that the parameters of these distributions, namely, the central latitude ($C$), width ($\delta$) and height ($H$) evolve with the cycle number and they show some hemispheric asymmetries. Although the asymmetries in these parameters persist for a few successive cycles, they get corrected within a few cycles and the new asymmetries appear again. In agreement with the previous study, we find that distribution parameters are correlated with the strengths of the cycles, although these correlations are significantly different in two hemispheres. The general trend that stronger cycles begin sunspot eruptions at relatively higher latitudes and have wider bands of sunspot emergence latitudes are confirmed when combining the data from two hemispheres. We explore these features using a flux transport dynamo model with stochastic fluctuations. We find that these features are correctly reproduced in this model. The solar cycle evolution of the distribution center is also in good agreement with observations. Possible explanations of the observed features based on this dynamo model are presented.


title: Life-Span of Semilinear Wave Equations with Scale-invariant Damping: Critical Strauss Exponent Case

abstract: The blow up problem of the semilinear scale-invariant damping wave equation with critical Strauss type exponent is investigated. The life span is shown to be: $T(\varepsilon)\leq C\exp(\varepsilon^{-2p(p-1)})$ when $p=p_S(n+\mu)$ for $0<\mu<\frac{n^2+n+2}{n+2}$. This result completes our previous study \cite{Tu-Lin} on the sub-Strauss type exponent $p<p_S(n+\mu)$. Our novelty is to construct the suitable test function from the modified Bessel function. This approach might be also applied to the other type damping wave equations.


title: The MOSDEF survey: a stellar mass-SFR-metallicity relation exists at $z\sim2.3$

abstract: We investigate the nature of the relation among stellar mass, star-formation rate, and gas-phase metallicity (the M$_*$-SFR-Z relation) at high redshifts using a sample of 260 star-forming galaxies at $z\sim2.3$ from the MOSDEF survey. We present an analysis of the high-redshift M$_*$-SFR-Z relation based on several emission-line ratios for the first time. We show that a M$_*$-SFR-Z relation clearly exists at $z\sim2.3$. The strength of this relation is similar to predictions from cosmological hydrodynamical simulations. By performing a direct comparison of stacks of $z\sim0$ and $z\sim2.3$ galaxies, we find that $z\sim2.3$ galaxies have $\sim0.1$ dex lower metallicity at fixed M$_*$ and SFR. In the context of chemical evolution models, this evolution of the M$_*$-SFR-Z relation suggests an increase with redshift of the mass-loading factor at fixed M$_*$, as well as a decrease in the metallicity of infalling gas that is likely due to a lower importance of gas recycling relative to accretion from the intergalactic medium at high redshifts. Performing this analysis simultaneously with multiple metallicity-sensitive line ratios allows us to rule out the evolution in physical conditions (e.g., N/O ratio, ionization parameter, and hardness of the ionizing spectrum) at fixed metallicity as the source of the observed trends with redshift and with SFR at fixed M$_*$ at $z\sim2.3$. While this study highlights the promise of performing high-order tests of chemical evolution models at high redshifts, detailed quantitative comparisons ultimately await a full understanding of the evolution of metallicity calibrations with redshift.


title: The multiset dimension of graphs

abstract: We introduce a variation of metric dimension, called the multiset dimension. The representation multiset of a vertex $v$ with respect to $W$ (which is a subset of the vertex set of a graph $G$), $r_m (v|W)$, is defined as a multiset of distances between $v$ and the vertices in $W$. If $r_m (u |W) \neq r_m(v|W)$ for every pair of distinct vertices $u$ and $v$, then $W$ is called an m-resolving set of $G$. If $G$ has an m-resolving set, then the cardinality of a smallest m-resolving set is called the multiset dimension of $G$, denoted by $md(G)$. If $G$ does not contain an m-resolving set, we write $md(G) = \infty$. In this paper we present basic results on the multiset dimension. We obtain some (sharp) bounds for multiset dimension of arbitrary graphs in term of its metric dimension, order, or diameter. We provide some necessary conditions for a graph to have finite multiset dimension, with an example of an infinite family of graphs where those necessary conditions are also sufficient. We also show that the multiset dimension of any graph other than a path is at least $3$ and finally we provide two families of graphs having the multiset dimension $3$.


title: Diffractive Scattering: Problems in Theory and Praxis

abstract: A concise survey of problems in modern strong interaction physics at high energies is given.


title: Vertex-Context Sampling for Weighted Network Embedding

abstract: In recent years, network embedding methods have garnered increasing attention because of their effectiveness in various information retrieval tasks. The goal is to learn low-dimensional representations of vertexes in an information network and simultaneously capture and preserve the network structure. Critical to the performance of a network embedding method is how the edges/vertexes of the network is sampled for the learning process. Many existing methods adopt a uniform sampling method to reduce learning complexity, but when the network is non-uniform (i.e. a weighted network) such uniform sampling incurs information loss. The goal of this paper is to present a generalized vertex sampling framework that works seamlessly with most existing network embedding methods to support weighted instead of uniform vertex/edge sampling. For efficiency, we propose a delicate sequential vertex-to-context graph data structure, such that sampling a training pair for learning takes only constant time. For scalability and memory efficiency, we design the graph data structure in a way that keeps space consumption low without requiring additional space. In addition to implementing existing network embedding methods, the proposed framework can be used to implement extensions that feature high-order proximity modeling and weighted relation modeling. Experiments conducted on three datasets, including a commercial large-scale one, verify the effectiveness and efficiency of the proposed weighted network embedding methods on a variety of tasks, including word similarity search, multi-label classification, and item recommendation.


title: Solid hulls and cores of weighted $H^\infty$-spaces

abstract: We determine the solid hull and solid core of weighted Banach spaces $H_v^\infty$ of analytic functions $f$ such that $v|f|$ is bounded, both in the case of the holomorphic functions on the disc and on the whole complex plane, for a very general class of radial weights $v$. Precise results are presented for concrete weights on the disc that could not be treated before. It is also shown that if $H_v^\infty$ is solid, then the monomials are an (unconditional) basis of the closure of the polynomials in $H_v^\infty$. As a consequence $H_v^\infty$ does not coincide with its solid hull and core in the case of the disc. An example shows that this does not hold for weighted spaces of entire functions.


title: Reducing Model Complexity for DNN Based Large-Scale Audio Classification

abstract: Audio classification is the task of identifying the sound categories that are associated with a given audio signal. This paper presents an investigation on large-scale audio classification based on the recently released AudioSet database. AudioSet comprises 2 millions of audio samples from YouTube, which are human-annotated with 527 sound category labels. Audio classification experiments with the balanced training set and the evaluation set of AudioSet are carried out by applying different types of neural network models. The classification performance and the model complexity of these models are compared and analyzed. While the CNN models show better performance than MLP and RNN, its model complexity is relatively high and undesirable for practical use. We propose two different strategies that aim at constructing low-dimensional embedding feature extractors and hence reducing the number of model parameters. It is shown that the simplified CNN model has only 1/22 model parameters of the original model, with only a slight degradation of performance.


title: On the $\Gamma$-equivalence of binary quadratic forms

abstract: For a congruence subgroup $\Gamma$, we define the notion of $\Gamma$-equivalence on binary quadratic forms which is the same as proper equivalence if $\Gamma = \mathrm{SL}_2(\mathbb Z)$. We develop a theory on $\Gamma$-equivalence such as the finiteness of $\Gamma$-reduced forms, the isomorphism between $\Gamma_0(N)$-form class group and the ideal class group, $N$-representation of integers, and $N$-genus of binary quadratic forms. As an application, we deal with representations of integers by binary quadratic forms under certain congruence condition on variables.


title: Dynamic Load Balancing Strategies for Graph Applications on GPUs

abstract: Acceleration of graph applications on GPUs has found large interest due to the ubiquitous use of graph processing in various domains. The inherent \textit{irregularity} in graph applications leads to several challenges for parallelization. A key challenge, which we address in this paper, is that of load-imbalance. If the work-assignment to threads uses node-based graph partitioning, it can result in skewed task-distribution, leading to poor load-balance. In contrast, if the work-assignment uses edge-based graph partitioning, the load-balancing is better, but the memory requirement is relatively higher. This makes it unsuitable for large graphs. In this work, we propose three techniques for improved load-balancing of graph applications on GPUs. Each technique brings in unique advantages, and a user may have to employ a specific technique based on the requirement. Using Breadth First Search and Single Source Shortest Paths as our processing kernels, we illustrate the effectiveness of each of the proposed techniques in comparison to the existing node-based and edge-based mechanisms.


title: Re-DPoctor: Real-time health data releasing with w-day differential privacy

abstract: Wearable devices enable users to collect health data and share them with healthcare providers for improved health service. Since health data contain privacy-sensitive information, unprotected data release system may result in privacy leakage problem. Most of the existing work use differential privacy for private data release. However, they have limitations in healthcare scenarios because they do not consider the unique features of health data being collected from wearables, such as continuous real-time collection and pattern preservation. In this paper, we propose Re-DPoctor, a real-time health data releasing scheme with $w$-day differential privacy where the privacy of health data collected from any consecutive $w$ days is preserved. We improve utility by using a specially-designed partition algorithm to protect the health data patterns. Meanwhile, we improve privacy preservation by applying newly proposed adaptive sampling technique and budget allocation method. We prove that Re-DPoctor satisfies $w$-day differential privacy. Experiments on real health data demonstrate that our method achieves better utility with strong privacy guarantee than existing state-of-the-art methods.


title: Super unitary representations revisited

abstract: With the usual definition of a super Hilbert space and a super unitary representation, it is easy to show that there are lots of super Lie groups for which the left-regular representation is not super unitary. I will argue that weakening the definition of a super Hilbert space (by allowing the super scalar product to be non-homogeneous, not just even) will allow the left-regular representation of all (connected) super Lie groups to be super unitary (with an adapted definition). Along the way I will introduce a (super) metric on a supermanifold that will allow me to define super and non-super scalar products on function spaces and I will show that the former are intimately related to the Hodge-star operation and the Fermionic Fourier transform. The latter also allows me to decompose certain super unitary representations as a direct integral over odd parameters of a family of super unitary representations depending on these odd parameters.


title: A Study of Fermi-LAT GeV gamma-ray Emission towards the Magnetar-harboring Supernova Remnant Kesteven 73 and Its Molecular Environment

abstract: We report our independent GeV gamma-ray study of the young shell-type supernova remnant (SNR) Kes 73 which harbors a central magnetar, and CO-line millimeter observations toward the SNR. Using 7.6 years of Fermi-LAT observation data, we detected an extended gamma-ray source ("source A") with the centroid on the west of the SNR, with a significance of 21.6 sigma in 0.1-300 GeV and an error circle of 5.4 arcminute in angular radius. The gamma-ray spectrum cannot be reproduced by a pure leptonic emission or a pure emission from the magnetar, and thus a hadronic emission component is needed. The CO-line observations reveal a molecular cloud (MC) at V_LSR~90 km/s, which demonstrates morphological correspondence with the western boundary of the SNR brightened in multiwavelength. The 12CO (J=2-1)/12CO (J=1-0) ratio in the left (blue) wing 85-88 km/s is prominently elevated to ~1.1 along the northwestern boundary, providing kinematic evidence of the SNR-MC interaction. This SNR-MC association yields a kinematic distance 9 kpc to Kes 73. The MC is shown to be capable of accounting for the hadronic gamma-ray emission component. The gamma-ray spectrum can be interpreted with a pure hadronic emission or a magnetar+hadronic hybrid emission. In the case of pure hadronic emission, the spectral index of the protons is 2.4, very similar to that of the radio-emitting electrons, essentially consistent with the diffusive shock acceleration theory. In the case of magnetar+hadronic hybrid emission, a magnetic field decay rate >= 10^36 erg/s is needed to power the magnetar's curvature radiation.


title: Group-Average and Convex Clustering for Partially Heterogeneous Linear Regression

abstract: In this paper, a subgroup least squares and a convex clustering are introduced for inferring a partially heterogenous linear regression that has potential application in the areas of precision marketing and precision medicine. The homogenous parameter and the subgroup-average of the heterogenous parameters can be consistently estimated by the subgroup least squares, without need of the sparsity assumption on the heterogenous parameters. The heterogenous parameters can be consistently clustered via the convex clustering. Unlike the existing methods for regression clustering, our clustering procedure is a standard mean clustering, although the model under study is a type of regression, and the corresponding algorithm only involves low dimensional parameters. Thus, it is simple and stable even if the sample size is large. The advantage of the method is further illustrated via simulation studies and the analysis of car sales data.


title: Spectroscopic properties of a two-dimensional time-dependent Cepheid model II. Determination of stellar parameters and abundances

abstract: Standard spectroscopic analyses of variable stars are based on hydrostatic one-dimensional model atmospheres. This quasi-static approach has theoretically not been validated. We aim at investigating the validity of the quasi-static approximation for Cepheid variables. We focus on the spectroscopic determination of the effective temperature $T_\mathrm{eff}$, surface gravity $\log \,g$, microturbulent velocity $\xi_\mathrm{t}$, and a generic metal abundance $\log\,A$ -- here taken as iron. We calculate a grid of 1D hydrostatic plane-parallel models covering the ranges in effective temperature and gravity encountered during the evolution of a two-dimensional time-dependent envelope model of a Cepheid computed with the radiation-hydrodynamics code CO5BOLD. We perform 1D spectral syntheses for artificial iron lines in local thermodynamic equilibrium varying the microturbulent velocity and abundance. We fit the resulting equivalent widths to corresponding values obtained from our dynamical model. For the four-parametric case, the stellar parameters are typically underestimated exhibiting a bias in the iron abundance of $\approx-0.2\,\mbox{dex}$. To avoid biases of this kind it is favourable to restrict the spectroscopic analysis to photometric phases $\phi_\mathrm{ph}\approx0.3\ldots 0.65$ using additional information to fix effective temperature and surface gravity. Hydrostatic 1D model atmospheres can provide unbiased estimates of stellar parameters and abundances of Cepheid variables for particular phases of their pulsations. We identified convective inhomogeneities as the main driver behind potential biases. For obtaining a complete view on the effects when determining stellar parameters with 1D models, multi-dimensional Cepheid atmosphere models are necessary for variables of longer period than investigated here.


title: Spin current noise of the spin Seebeck effect and spin pumping

abstract: We theoretically investigate the fluctuation of a pure spin current induced by the spin Seebeck effect and spin pumping in a normal metal (NM)/ferromagnet (FM) bilayer system. Starting with a simple FI--NM interface model with both spin-conserving and spin-non-conserving processes, we derive general expressions of the spin current and the spin-current noise at the interface within second-order perturbation of the FI--NM coupling strength, and estimate them for an yttrium iron garnet (YIG) --platinum interface. We show that the spin-current noise can be used to determine the effective spin carried by a magnon, modified by the spin-non-conserving process at the interface. In addition, we show that it provides information on the effective spin of a magnon, heating at the interface under spin pumping, and spin Hall angle of the NM.


title: 3D-SSD: Learning Hierarchical Features from RGB-D Images for Amodal 3D Object Detection

abstract: This paper aims at developing a faster and a more accurate solution to the amodal 3D object detection problem for indoor scenes. It is achieved through a novel neural network that takes a pair of RGB-D images as the input and delivers oriented 3D bounding boxes as the output. The network, named 3D-SSD, composed of two parts: hierarchical feature fusion and multi-layer prediction. The hierarchical feature fusion combines appearance and geometric features from RGB-D images while the multi-layer prediction utilizes multi-scale features for object detection. As a result, the network can exploit 2.5D representations in a synergetic way to improve the accuracy and efficiency. The issue of object sizes is addressed by attaching a set of 3D anchor boxes with varying sizes to every location of the prediction layers. At the end stage, the category scores for 3D anchor boxes are generated with adjusted positions, sizes and orientations respectively, leading to the final detections using non-maximum suppression. In the training phase, the positive samples are identified with the aid of 2D ground truth to avoid the noisy estimation of depth from raw data, which guide to a better converged model. Experiments performed on the challenging SUN RGB-D dataset show that our algorithm outperforms the state-of-the-art Deep Sliding Shape by 10.2% mAP and 88x faster. Further, experiments also suggest our approach achieves comparable accuracy and is 386x faster than the state-of-art method on the NYUv2 dataset even with a smaller input image size.


title: Secure Classification With Augmented Features

abstract: With the evolution of data collection ways, it is possible to produce abundant data described by multiple feature sets. Previous studies show that including more features does not necessarily bring positive effect. How to prevent the augmented features worsening classification performance is crucial but rarely studied. In this paper, we study this challenging problem by proposing a secure classification approach, whose accuracy is never degenerated when exploiting augmented features. We propose two ways to achieve the security of our method named as SEcure Classification (SEC). Firstly, to leverage augmented features, we learn various types of classifiers and adapt them by employing a specially designed robust loss. It provides various candidate classifiers to meet the following assumption of security operation. Secondly, we integrate all candidate classifiers by approximately maximizing the performance improvement. Under a mild assumption, the integrated classifier has theoretical security guarantee. Several new optimization methods have been developed to accommodate the problems with proved convergence. Besides evaluating SEC on 16 data sets, we also apply SEC in the application of diagnostic classification of schizophrenia since it has vast application potentiality. Experimental results demonstrate the effectiveness of SEC in both tackling security problem and discriminating schizophrenic patients from healthy controls.


title: Atmospheric gas dynamics in the Perseus cluster observed with Hitomi

abstract: Extending the earlier measurements reported in Hitomi collaboration (2016, Nature, 535, 117), we examine the atmospheric gas motions within the central 100~kpc of the Perseus cluster using observations obtained with the Hitomi satellite. After correcting for the point spread function of the telescope and using optically thin emission lines, we find that the line-of-sight velocity dispersion of the hot gas is remarkably low and mostly uniform. The velocity dispersion reaches maxima of approximately 200~km~s$^{-1}$ toward the central active galactic nucleus (AGN) and toward the AGN inflated north-western `ghost' bubble. Elsewhere within the observed region, the velocity dispersion appears constant around 100~km~s$^{-1}$. We also detect a velocity gradient with a 100~km~s$^{-1}$ amplitude across the cluster core, consistent with large-scale sloshing of the core gas. If the observed gas motions are isotropic, the kinetic pressure support is less than 10\% of the thermal pressure support in the cluster core. The well-resolved optically thin emission lines have Gaussian shapes, indicating that the turbulent driving scale is likely below 100~kpc, which is consistent with the size of the AGN jet inflated bubbles. We also report the first measurement of the ion temperature in the intracluster medium, which we find to be consistent with the electron temperature. In addition, we present a new measurement of the redshift to the brightest cluster galaxy NGC~1275.


title: Extending the Propagation Distance of a Silver Nanowire Plasmonic Waveguide with a Dielectric Multilayer Substrate

abstract: Chemical synthesized silver nanowires have been proved to be the efficient architecture for Plasmonic waveguides, but the high propagation loss prevents their widely applications. Here, we demonstrate that the propagation distance of the plasmons along the Ag NW can be extended if the Ag NW was placed on a dielectric multilayer substrate containing a photonic band gap, but not placed on a commonly used glass substrate. The propagation distance at 630 nm wavelength can reach 16 um even that the Ag NW is as thin as 90 nm in diameter. Experimental and simulation results further show that the polarization of this propagating plasmon mode was nearly parallel to the surface of the dielectric multilayer, so it was excited by a transverse-electric polarized Bloch surface wave propagating along a polymer nanowire with diameter at only about 170 nm on the same dielectric multilayer. Numerical simulations were also carried out and consistent with the experiment results. Our work provides a platform to extend the propagation distance of plasmonic waveguide and also for the integration between photonic and plasmonic waveguides on the nanometre scale.


title: On an inverse elastic wave imaging scheme for nearly incompressible materials

abstract: This paper is devoted to the algorithmic development of inverse elastic scattering problems. We focus on reconstructing the locations and shapes of elastic scatterers with known dictionary data for the nearly incompressible materials. The scatterers include non-penetrable rigid obstacles and penetrable mediums, and we use time-harmonic elastic point signals as the incident input waves. The scattered waves are collected in a relatively small backscattering aperture on a bounded surface. A two-stage algorithm is proposed for the reconstruction and only two incident waves of different wavenumbers are required. The unknown scatterer is first approximately located by using the measured data at a small wavenumber, and then the shape of the scatterer is determined by the computed location of the scatterer together with the measured data at a regular wavenumber. The corresponding mathematical principle with rigorous analysis is presented. Numerical tests illustrate the effectiveness and efficiency of the proposed method.


title: Scattered Short Gamma-Ray Bursts as Electromagnetic Counterparts to Gravitational Waves and Implications of GW170817 and GRB 170817A

abstract: In the faint short gamma-ray burst sGRB 170817A followed by the gravitational waves (GWs) from a merger of two neutron stars (NSs) GW170817, the spectral peak energy is too high to explain only by canonical off-axis emission. We investigate off-axis appearance of an sGRB prompt emission scattered by a cocoon, which is produced through the jet-merger-ejecta interaction, with either sub-relativistic or mildly-relativistic velocities. We show that the observed properties of sGRB 170817A, in particular the high peak energy, can be consistently explained by the Thomson-scattered emission with a typical sGRB jet, together by its canonical off-axis emission, supporting that an NS-NS merger is the origin of sGRBs. The scattering occurs at $\lesssim 10^{10}$--$10^{12}\,{\rm cm}$ not far from the central engine, implying the photospheric or internal shock origin of the sGRB prompt emission. The boundary between the jet and cocoon is sharp, which could be probed by future observations of off-axis afterglows. The scattering model predicts a distribution of the spectral peak energy that is similar to the observed one but with a cutoff around $\sim$ MeV energy, and its correlations with the luminosity, duration, and time lag from GWs, providing a way to distinguish it from alternative models.


title: Efficient Inferencing of Compressed Deep Neural Networks

abstract: Large number of weights in deep neural networks makes the models difficult to be deployed in low memory environments such as, mobile phones, IOT edge devices as well as "inferencing as a service" environments on cloud. Prior work has considered reduction in the size of the models, through compression techniques like pruning, quantization, Huffman encoding etc. However, efficient inferencing using the compressed models has received little attention, specially with the Huffman encoding in place. In this paper, we propose efficient parallel algorithms for inferencing of single image and batches, under various memory constraints. Our experimental results show that our approach of using variable batch size for inferencing achieves 15-25\% performance improvement in the inference throughput for AlexNet, while maintaining memory and latency constraints.


title: Star-Forming Environments Throughout the M101 Group

abstract: We present a multiwavelength study of star formation within the nearby M101 Group, including new deep H$\alpha$ imaging of M101 and its two companions. We perform a statistical analysis of the H$\alpha$ to FUV flux ratios in HII regions located in three different environments: M101's inner disk, M101's outer disk, and M101's lower mass companion galaxy NGC~5474. We find that, once bulk radial trends in extinction are taken into account, both the median and scatter in F$_{\rm H\alpha}/$F$_{\rm FUV}$ in HII regions are invariant across all of these environments. Also, using Starburst99 models, we are able to qualitatively reproduce the distributions of F$_{\rm H\alpha}/$F$_{\rm FUV}$ throughout these different environments using a standard Kroupa IMF, hence we find no need to invoke truncations in the upper mass end of the IMF to explain the young star-forming regions in the M101 Group even at extremely low surface density. This implies that star formation in low density environments differs from star formation in high density environments only by intensity and not by cloud-to-cloud physics.


title: Output Consensus of Networked Hammerstein and Wiener Systems

abstract: In this paper we consider the output consensus problem of networked Hammerstein and Wiener systems in a noisy environment. The Hammerstein or Wiener system is assumed to be open-loop stable, and its static nonlinearity is allowed to grow up but not faster than a polynomial. A control algorithm based on the distributed stochastic approximation algorithm with expanding truncations is designed and it is shown that under the designed control the output consensus is achieved. The numerical simulation given in the paper justifies the theoretical assertions.


title: Improved Text Language Identification for the South African Languages

abstract: Virtual assistants and text chatbots have recently been gaining popularity. Given the short message nature of text-based chat interactions, the language identification systems of these bots might only have 15 or 20 characters to make a prediction. However, accurate text language identification is important, especially in the early stages of many multilingual natural language processing pipelines. This paper investigates the use of a naive Bayes classifier, to accurately predict the language family that a piece of text belongs to, combined with a lexicon based classifier to distinguish the specific South African language that the text is written in. This approach leads to a 31% reduction in the language detection error. In the spirit of reproducible research the training and testing datasets as well as the code are published on github. Hopefully it will be useful to create a text language identification shared task for South African languages.


title: Query-free Clothing Retrieval via Implicit Relevance Feedback

abstract: Image-based clothing retrieval is receiving increasing interest with the growth of online shopping. In practice, users may often have a desired piece of clothing in mind (e.g., either having seen it before on the street or requiring certain specific clothing attributes) but may be unable to supply an image as a query. We model this problem as a new type of image retrieval task in which the target image resides only in the user's mind (called "mental image retrieval" hereafter). Because of the absence of an explicit query image, we propose to solve this problem through relevance feedback. Specifically, a new Bayesian formulation is proposed that simultaneously models the retrieval target and its high-level representation in the mind of the user (called the "user metric" hereafter) as posterior distributions of pre-fetched shop images and heterogeneous features extracted from multiple clothing attributes, respectively. Requiring only clicks as user feedback, the proposed algorithm is able to account for the variability in human decision-making. Experiments with real users demonstrate the effectiveness of the proposed algorithm.


title: Atomistic potential for graphene and other sp$^2$ carbon systems

abstract: We introduce a torsional force field for sp$^2$ carbon to augment an in-plane atomistic potential of a previous work (Kalosakas et al, J. Appl. Phys. {\bf 113}, 134307 (2013)) so that it is applicable to out-of-plane deformations of graphene and related carbon materials. The introduced force field is fit to reproduce DFT calculation data of appropriately chosen structures. The aim is to create a force field that is as simple as possible so it can be efficient for large scale atomistic simulations of various sp$^2$ carbon structures without significant loss of accuracy. We show that the complete proposed potential reproduces characteristic properties of fullerenes and carbon nanotubes. In addition, it reproduces very accurately the out-of-plane ZA and ZO modes of graphene's phonon dispersion as well as all phonons with frequencies up to 1000~cm$^{-1}$.


title: Structural Bifurcation Analysis in Chemical Reaction Networks

abstract: In living cells, chemical reactions form a complex network. Complicated dynamics arising from such networks are the origins of biological functions. We propose a novel mathematical method to analyze bifurcation behaviors of a reaction system from the network structure alone. The whole network is decomposed into subnetworks based on "buffering structures". For each subnetwork, the bifurcation condition is studied independently, and the parameters that can induce bifurcations and the chemicals that can exhibit bifurcations are determined. We demonstrate our theory using hypothetical and real networks.


title: Unifying different interpretations of the nonlinear response in glass-forming liquids

abstract: This work aims at reconsidering several interpretations coexisting in the recent literature concerning non-linear susceptibilities in supercooled liquids. We present experimental results on glycerol and propylene carbonate showing that the three independent cubic susceptibilities have very similar frequency and temperature dependences, both for their amplitudes and phases. This strongly suggests a unique physical mechanism responsible for the growth of these non-linear susceptibilities. We show that the framework proposed by two of us [BB, Phys. Rev. B 72, 064204 (2005)], where the growth of non-linear susceptibilities is intimately related to the growth of "glassy domains", accounts for all the salient experimental features. We then review several complementary and/or alternative models, and show that the notion of cooperatively rearranging glassy domains is a key (implicit or explicit) ingredient to all of them. This paves the way for future experiments which should deepen our understanding of glasses.


title: Pressure-Induced Superconductivity from Doping-Induced Antiferromagnetic Phase of 112-type Ca$_{1-x}$La$_{x}$FeAs$_{2}$

abstract: The effects of pressure on antiferromagnetic (AFM) and superconducting phase transitions of 112-type Ca$_{1-x}$La$_{x}$FeAs$_{2}$ were studied, and the in-plane electrical resistivity $\rho_{ab}$ was measured with an indenter-type pressure cell. The AFM phase transition temperatures of $T_{\rm N}$ = 47, 63, and 63 K at ambient pressure for $x$ = 0.18, 0.21, and 0.26 was suppressed by applying pressure $P$, with superconductivity emerging at critical pressures of $P_{\rm c}$ $\simeq$ 0, 1.5, and 3.4 GPa, respectively. Correspondingly, the slope of $T_{\rm N}$ against $P$ decreased as $dT_{\rm N}/P$ $\simeq$ $-$15 and $-$2 K/GPa for $x$ = 0.21 and 0.26, respectively. Thus, although the AFM phase was stabilized with La doping $x$, the AFM phase was suppressed by pressure, and superconductivity eventually emerged.


title: Adversarial Learning of Structure-Aware Fully Convolutional Networks for Landmark Localization

abstract: Landmark/pose estimation in single monocular images have received much effort in computer vision due to its important applications. It remains a challenging task when input images severe occlusions caused by, e.g., adverse camera views. Under such circumstances, biologically implausible pose predictions may be produced. In contrast, human vision is able to predict poses by exploiting geometric constraints of landmark point inter-connectivity. To address the problem, by incorporating priors about the structure of pose components, we propose a novel structure-aware fully convolutional network to implicitly take such priors into account during training of the deep network. Explicit learning of such constraints is typically challenging. Instead, inspired by how human identifies implausible poses, we design discriminators to distinguish the real poses from the fake ones (such as biologically implausible ones). If the pose generator G generates results that the discriminator fails to distinguish from real ones, the network successfully learns the priors. Training of the network follows the strategy of conditional Generative Adversarial Networks (GANs). The effectiveness of the proposed network is evaluated on three pose-related tasks: 2D single human pose estimation, 2D facial landmark estimation and 3D single human pose estimation. The proposed approach significantly outperforms the state-of-the-art methods and almost always generates plausible pose predictions, demonstrating the usefulness of implicit learning of structures using GANs.


title: On anomalous localized resonance and plasmonic cloaking beyond the quasistatic limit

abstract: In this paper, we give the mathematical construction of novel core-shell plasmonic structures that can induce anomalous localized resonance and invisibility cloaking at certain finite frequencies beyond the quasistatic limit. The crucial ingredient in our study is that the plasmon constant and the loss parameter are constructed in a delicate way that are correlated and depend on the source and the size of the plasmonic structure. As a significant byproduct of this study, we also derive the complete spectrum of the Neumann-Poinc\'are operator associated to the Helmholtz equation with finite frequencies in the radial geometry. The spectral result is the first one in its type and is of significant mathematical interest for its own sake.


title: An extended hybrid numerical simulation of isotropic compressible turbulence

abstract: This paper presents an extension of the hybrid scheme proposed by Wang et al. (J. Comput. Phys. 229 (2010) 169-180) for numerical simulation of compressible isotropic turbulence to flows with higher turbulent Mach numbers. The scheme still utilizes an 8th-order compact scheme with built-in hyperviscosity for smooth regions and a 7th-order WENO scheme for highly compressive regions, but now both in their conservation formulations and for the latter with the Roe type characteristic-wise reconstruction. To enhance the robustness of the WENO scheme without compromising its high-resolution and accuracy, the recursive-order-reduction procedure is adopted, where a new type of reconstruction-failure-detection criterion is constructed. To capture the upwind direction properly in extreme conditions, the global Lax-Friedrichs numerical flux is used. In addition, a new form of cooling function is proposed, which is proved to be positivity-preserving. With these techniques, the new scheme not only inherits the good properties of the original one but also extends largely the computable range of turbulent Mach number, which has been further confirmed by numerical results.


title: Sensitivity of Coronal Loop Sausage Mode Frequencies and Decay Rates to Radial and Longitudinal Density Inhomogeneities: A Spectral Approach

abstract: Fast sausage modes in solar magnetic coronal loops are only fully contained in unrealistically short dense loops. Otherwise they are leaky, losing energy to their surrounds as outgoing waves. This causes any oscillation to decay exponentially in time. Simultaneous observations of both period and decay rate therefore reveal the eigenfrequency of the observed mode, and potentially insight into the tubes' nonuniform internal structure. In this article, a global spectral description of the oscillations is presented that results in an implicit matrix eigenvalue equation where the eigenvalues are associated predominantly with the diagonal terms of the matrix. The off-diagonal terms vanish identically if the tube is uniform. A linearized perturbation approach, applied with respect to a uniform reference model, is developed that makes the eigenvalues explicit. The implicit eigenvalue problem is easily solved numerically though, and it is shown that knowledge of the real and imaginary parts of the eigenfrequency is sufficient to determine the width and density contrast of a boundary layer over which the tubes' enhanced internal densities drop to ambient values. Linearized density kernels are developed that show sensitivity only to the extreme outside of the loops for radial fundamental modes, especially for small density enhancements, with no sensitivity to the core. Higher radial harmonics do show some internal sensitivity, but these will be more difficult to observe. Only kink modes are sensitive to the tube centres. {Variation in internal and external Alfv\'en speed along the loop is shown to have little effect on the fundamental dimensionless eigenfrequency, though the associated eigenfunction becomes more compact at the loop apex as stratification increases, or may even displace from the apex.


title: Bluetooth 5: a concrete step forward towards the IoT

abstract: Six years after the adoption of the standard 4.0, the Bluetooth Special Interest Group (SIG), a non-profit association that deals with the study and the development of technology standards including those of Bluetooth, has officially released the main features of Bluetooth 5.0. It is one of the significant developments in short-range wireless communication technology. As stated by the SIG, the new standard will forever change the way people approach the Internet of Things (IoT), turning it into something that takes place around them in an almost natural and transparent way. In this article, the future IoT scenarios and use cases that justify the push for Bluetooth 5 are introduced. A set of new technical features that are included in Bluetooth 5 are presented, and their advantages and drawbacks are described.


title: Smooth Neighbors on Teacher Graphs for Semi-supervised Learning

abstract: The recently proposed self-ensembling methods have achieved promising results in deep semi-supervised learning, which penalize inconsistent predictions of unlabeled data under different perturbations. However, they only consider adding perturbations to each single data point, while ignoring the connections between data samples. In this paper, we propose a novel method, called Smooth Neighbors on Teacher Graphs (SNTG). In SNTG, a graph is constructed based on the predictions of the teacher model, i.e., the implicit self-ensemble of models. Then the graph serves as a similarity measure with respect to which the representations of "similar" neighboring points are learned to be smooth on the low-dimensional manifold. We achieve state-of-the-art results on semi-supervised learning benchmarks. The error rates are 9.89%, 3.99% for CIFAR-10 with 4000 labels, SVHN with 500 labels, respectively. In particular, the improvements are significant when the labels are fewer. For the non-augmented MNIST with only 20 labels, the error rate is reduced from previous 4.81% to 1.36%. Our method also shows robustness to noisy labels.


title: Rarity of extremal edges in random surfaces and other theoretical applications of cluster algorithms

abstract: Motivated by questions on the delocalization of random surfaces, we prove that random surfaces satisfying a Lipschitz constraint rarely develop extremal gradients. Previous proofs of this fact relied on reflection positivity and were thus limited to random surfaces defined on highly symmetric graphs, whereas our argument applies to general graphs. Our proof makes use of a cluster algorithm and reflection transformation for random surfaces of the type introduced by Swendsen-Wang, Wolff and Evertz et al. We discuss the general framework for such cluster algorithms, reviewing several particular cases with emphasis on their use in obtaining theoretical results. Two additional applications are presented: A reflection principle for random surfaces and a proof that pair correlations in the spin $O(n)$ model have monotone densities, strengthening Griffiths' first inequality for such correlations.


title: The scaling and skewness of optimally transported meshes on the sphere

abstract: In the context of numerical solution of PDEs, dynamic mesh redistribution methods (r-adaptive methods) are an important procedure for increasing the resolution in regions of interest, without modifying the connectivity of the mesh. Key to the success of these methods is that the mesh should be sufficiently refined (locally) and flexible in order to resolve evolving solution features, but at the same time not introduce errors through skewness and lack of regularity. Some state-of-the-art methods are bottom-up in that they attempt to prescribe both the local cell size and the alignment to features of the solution. However, the resulting problem is overdetermined, necessitating a compromise between these conflicting requirements. An alternative approach, described in this paper, is to prescribe only the local cell size and augment this an optimal transport condition to provide global regularity. This leads to a robust and flexible algorithm for generating meshes fitted to an evolving solution, with minimal need for tuning parameters. Of particular interest for geophysical modelling are meshes constructed on the surface of the sphere. The purpose of this paper is to demonstrate that meshes generated on the sphere using this optimal transport approach have good a-priori regularity and that the meshes produced are naturally aligned to various simple features. It is further shown that the sphere's intrinsic curvature leads to more regular meshes than the plane. In addition to these general results, we provide a wide range of examples relevant to practical applications, to showcase the behaviour of optimally transported meshes on the sphere. These range from axisymmetric cases that can be solved analytically to more general examples that are tackled numerically. Evaluation of the singular values and singular vectors of the mesh transformation provides a quantitative measure of the mesh aniso...


title: Transient Behavior of Redox Flow Battery Connected to Circuit Based on Global Phase Structure

abstract: A Redox Flow Battery (RFB) is one of the promising energy storage systems in power grid. An RFB has many advantages such as a quick response, a large capacity, and a scalability. Due to these advantages, an RFB can operate in mixed time scale. Actually, it has been demonstrated that an RFB can be used for load leveling, compensating sag, and smoothing the output of the renewable sources. An analysis on transient behaviors of an RFB is a key issue for these applications. An RFB is governed by electrical, chemical, and fluid dynamics. The hybrid structure makes the analysis difficult. To analyze transient behaviors of an RFB, the exact model is necessary. In this paper, we focus on a change in a concentration of ions in the electrolyte, and simulate the change with a model which is mainly based on chemical kinetics. The simulation results introduces transient behaviors of an RFB in a response to a load variation. There are found three kinds of typical transient behaviors including oscillations. As results, it is clarified that the complex transient behaviors, due to slow and fast dynamics in the system, arise by the quick response to load.


title: Non-expansive bijections to the unit ball of $\ell_1$-sum of strictly convex Banach spaces

abstract: Extending recent results by Cascales, Kadets, Orihuela and Wingler (2016), Kadets and Zavarzina (2017), and Zavarzina (2017) we demonstrate that for every Banach space $X$ and every collection $Z_i, i\in I$ of strictly convex Banach spaces every non-expansive bijection from the unit ball of $X$ to the unit ball of sum of $Z_i$ by $\ell_1$ is an isometry.


title: Tunable topological phases with fermionic atoms in a one-dimensional flux lattice

abstract: We present a simple scheme for implementing a one-dimensional (1D) magnetic-flux lattice of ultracold fermionic spin-$1/2$ atoms. The resulting tight-binding model supports gapped and gapless topological phases, and chiral currents for Meissner and vortex phases. Its single-particle spectra exhibit topological flat bands at small flux, and the flatness sensitively depends on hopping strength. An effective $p$-wave interaction arises in a $s$-wave paired superfluid. Treating atomic internal states as forming a synthetic dimension and balancing the interplay of magnetic flux and Zeeman field, our model describes a tunable topological Fermi superfluid, which paves the way towards experimental explorations of non-Abelian topological matter in 1D atomic quantum gases.


title: Role of sufficient statistics in stochastic thermodynamics and its implication to sensory adaptation

abstract: A sufficient statistic is a significant concept in statistics, which means a probability variable that has sufficient information required for an inference task. We investigate the roles of sufficient statistics and related quantities in stochastic thermodynamics. Specifically, we prove that for general continuous-time bipartite networks, the existence of a sufficient statistic implies that an informational quantity called the sensory capacity takes the maximum. Since the maximal sensory capacity imposes a constraint that the energetic efficiency cannot exceed one-half, our result implies that the existence of a sufficient statistic is inevitably accompanied by energetic dissipation. We also show that, in a particular parameter region of linear Langevin systems, there exists the optimal noise intensity, at which the sensory capacity, the information-thermodynamic efficiency, and the total entropy production are optimized at the same time. We apply our general result to a model of sensory adaptation of E. Coli, and find that the sensory capacity is nearly maximal with experimentally realistic parameters, suggesting that E. Coli approximately realizes a sufficient statistic in signal transduction at the cost of energetic dissipation.


title: $\eta_c'$ Hadroproduction at Next-to-Leading Order and its Relevance to $\psi'$ Production

abstract: We proceed to the first study of $\eta_c'$ prompt hadroproduction at next-to-leading order in $\alpha_s$. Based on heavy-quark-spin symmetry, which is systematically used in quarkonium-production phenomenology, we demonstrate that prompt $\eta_c'$ can be studied at the LHC with the existing data. We emphasize its relevance to constrain $\psi'$ production, in the same way as the first prompt $\eta_c'$ data at the LHC lately strongly impacted the phenomenology of $J/\psi$ studies.


title: High Energy Transients: the Millisecond Domain

abstract: Search for high energy transients in the millisecond domain has come to the focus in recent times due to the detection of Gravitational Wave events and the identification of Fast Radio Bursts as cosmological sources. I will highlight the sensitivity limitations in the currently operating hard X-ray telescopes and give some details of the search for millisecond events in the AstroSat CZT Imager data.


title: Acquiring Target Stacking Skills by Goal-Parameterized Deep Reinforcement Learning

abstract: Understanding physical phenomena is a key component of human intelligence and enables physical interaction with previously unseen environments. In this paper, we study how an artificial agent can autonomously acquire this intuition through interaction with the environment. We created a synthetic block stacking environment with physics simulation in which the agent can learn a policy end-to-end through trial and error. Thereby, we bypass to explicitly model physical knowledge within the policy. We are specifically interested in tasks that require the agent to reach a given goal state that may be different for every new trial. To this end, we propose a deep reinforcement learning framework that learns policies which are parametrized by a goal. We validated the model on a toy example navigating in a grid world with different target positions and in a block stacking task with different target structures of the final tower. In contrast to prior work, our policies show better generalization across different goals.


title: Axiomatic quantum mechanics: Necessity and benefits for the physics studies

abstract: The ongoing progress in quantum theory emphasizes the crucial role of the very basic principles of quantum theory. However, this is not properly followed in teaching quantum mechanics on the graduate and undergraduate levels of physics studies. The existing textbooks typically avoid the axiomatic presentation of the theory. We emphasize usefulness of the systematic, axiomatic approach to the basics of quantum theory as well as its importance in the light of the modern scientific-research context.


title: $p$-adic Dynamics of Hecke Operators on Modular Curves

abstract: In this paper we study the $p$-adic dynamics of prime-to-$p$ Hecke operators on the set of points of modular curves in both cases of good ordinary and supersingular reduction. We pay special attention to the dynamics on the set of CM points. In the case of ordinary reduction we employ the Serre-Tate coordinates, while in the case of supersingular reduction, we use a parameter on the deformation space of the unique formal group of height $2$ over $\overline{\mathbb{F}}_p$, and take advantage of the Gross-Hopkins period map.


title: Determination of Checkpointing Intervals for Malleable Applications

abstract: Selecting optimal intervals of checkpointing an application is important for minimizing the run time of the application in the presence of system failures. Most of the existing efforts on checkpointing interval selection were developed for sequential applications while few efforts deal with parallel applications where the applications are executed on the same number of processors for the entire duration of execution. Some checkpointing systems support parallel applications where the number of processors on which the applications execute can be changed during the execution. We refer to these kinds of parallel applications as {\em malleable} applications. In this paper, we develop a performance model for malleable parallel applications that estimates the amount of useful work performed in unit time (UWT) by a malleable application in the presence of failures as a function of checkpointing interval. We use this performance model function with different intervals and select the interval that maximizes the UWT value. By conducting a large number of simulations with the traces obtained on real supercomputing systems, we show that the checkpointing intervals determined by our model can lead to high efficiency of applications in the presence of failures.


title: On hp-Streamline Diffusion and Nitsche schemes for the Relativistic Vlasov-Maxwell System

abstract: We study stability and convergence of $hp$-streamline diffusion (SD) finite element, and Nitsche's schemes for the three dimensional, relativistic (3 spatial dimension and 3 velocities), time dependent Vlasov-Maxwell system and Maxwell's equations, respectively. For the $hp$ scheme for the Vlasov-Maxwell system, assuming that the exact solution is in the Sobolev space $H^{s+1}(\Omega)$, we derive global {\sl a priori} error bound of order ${\mathcal O}(h/p)^{s+1/2}$, where $h (= \max_K h_K)$ is the mesh parameter and $p (= \max_K p_K)$ is the spectral order. This estimate is based on the local version with $h_K=\mbox{ diam } K$ being the diameter of the {\sl phase-space-time} element $K$ and $p_K$ is the spectral order (the degree of approximating finite element polynomial) for $K$. As for the Nitsche's scheme, by a simple calculus of the field equations, first we convert the Maxwell's system to an {\sl elliptic type} equation. Then, combining the Nitsche's method for the spatial discretization with a second order time scheme, we obtain optimal convergence of ${\mathcal O}(h^2+k^2)$, where $h$ is the spatial mesh size and $k$ is the time step. Here, as in the classical literature, the second order time scheme requires higher order regularity assumptions. Numerical justification of the results, in lower dimensions, is presented and is also the subject of a forthcoming computational work [20].


title: Secondary fields and partial wave expansion. Self consistency conditions in a conformal model

abstract: A nontrivial conformally invariant model is obtained via generalization the method of obtaining conformally invariant models in $2D$ Euclidean space to the Euclidean space with dimension $D>2$. This method was previously developed by E.S. Fradkin and M.Ya. Palchik (see [7] and reference therein). The partial wave expansion of a four-point function $\langle j_{\mu}(x_{1})j_{\nu}(x_{2})\varphi(x_{3})\varphi(x_{4})\rangle$ containing two conserved vector fields $j_{\mu}$ and two scalars $\varphi$ of dimension $ d $ in a $D$ -dimensional Euclidean space is considered. The requirement of the absence the vector operator of the dimension $d + 1$ in this expansion allows us to find the relationship between all the coupling constants in such a model.


title: Saturation of a turbulent mixing layer over a cavity: response to harmonic forcing around mean flows

abstract: Turbulent mixing layers over cavities can couple with acoustic waves and lead to undesired oscillations. To understand the nonlinear aspects of this phenomenon, a turbulent mixing layer over a deep cavity at Reynolds number 150 000 is considered and its response to harmonic forcing is analysed with large-eddy simulations (LES) and linearised Navier-Stokes equations (LNSE). As a model of incoming acoustic perturbations, spatially uniform time-harmonic forcing is applied at the cavity end, with amplitudes in the wide range 0.045-8.9% of the bulk velocity. Compressible LES provide reference nonlinear responses of the shear layer, and the associated mean flows. Linear responses are calculated with the incompressible LNSE around the LES mean flows; they predict well the amplification (both measured with kinetic energy and with a proxy for vortex sound production) and capture the nonlinear saturation observed as the forcing amplitude increases and the mixing layer thickens. Perhaps surprisingly, LNSE calculations based on a monochromatic (single frequency) assumption yield a good agreement even though higher harmonics and their nonlinear interaction (Reynolds stresses) are not negligible. However, the leading Reynolds stresses do not force the mixing layer efficiently, as shown by a comparison with the optimal volume forcing obtained in a resolvent analysis. Thus, they cannot fully benefit from the potential for amplification available in the flow. Finally, the sensitivity of the optimal harmonic forcing at the cavity end is computed with an adjoint method. The sensitivities to mean flow modification and to a localised feedback (structural sensitivity) both identify the upstream cavity corner as the region where a small-amplitude modification has the strongest effect. This can guide in a systematic way the design of strategies for the control of amplification and saturation mechanisms.


title: Flux large deviations of weakly interacting jump processes via well-posedness of an associated Hamilton-Jacobi equation

abstract: We establish uniqueness for a class of first-order Hamilton-Jacobi equations with Hamiltonians that arise from the large deviations of the empirical measure and empirical flux pair of weakly interacting Markov jump processes. As a corollary we obtain a large deviation principle for the trajectory of the empirical measure and empirical flux pair of such processes. As a second corollary we get same result in the setting where the jump-rates are time-periodic, with period-length that decreases to 0.


title: Fast Dynamic Arrays

abstract: We present a highly optimized implementation of tiered vectors, a data structure for maintaining a sequence of $n$ elements supporting access in time $O(1)$ and insertion and deletion in time $O(n^\epsilon)$ for $\epsilon > 0$ while using $o(n)$ extra space. We consider several different implementation optimizations in C++ and compare their performance to that of vector and multiset from the standard library on sequences with up to $10^8$ elements. Our fastest implementation uses much less space than multiset while providing speedups of $40\times$ for access operations compared to multiset and speedups of $10.000\times$ compared to vector for insertion and deletion operations while being competitive with both data structures for all other operations.


title: Modelling of Vortex-Induced Aviation Turbulence

abstract: Aviation turbulence is modelled as an interaction between an aircraft and a vortex tube. The vortex tube can have an arbitrary orientation/offset with respect to the aircraft. We compare modelling the aircraft (i) as a point and (ii) having a finite area (wing and fuselage). We consider both vertical and horizontal acceleration experienced by the aircraft. The baseline vortex tube has an area which is of the order of the aircraft area.


title: Approximating the nonlinear Schr\"odinger equation by a two level linearly implicit finite element method

abstract: We consider the study of a numerical scheme for an initial- and Dirichlet boundary- value problem for a nonlinear Schr\"odinger equation. We approximate the solution using a, local (non-uniform) two level scheme in time (see C. Besse [6] and [7]) combined with, an optimal, finite element strategy for the discretization in the spatial variable based on studies outlined as, e.g. in [2] and [10]. For the proposed fully discrete scheme, we show convergence both in $L_2 $ and $H^1$ norms.


title: High-Temperature Superconductivity in Th-H System at Pressure Conditions

abstract: New stable thorium decahydride $Fm\bar{3}m$-$ThH_{10}$, a record high-temperature superconductor with $T_C$ up to 241 K ($-32^{\circ}C$), critical field $H_C$ up to 71T and superconducting gap ${\Delta_0}$ = 52 meV at 80-100 GPa was predicted by evolutionary algorithm USPEX. Another phase ${P2_1c}$-$ThH_7$ was found to be a high-temperature superconductor with $T_C$ of 65 K. Analysis of superconducting state was performed within Eliashberg formalism and the dependencies of $H_C(T)$, ${\Delta_0}(T)$, $T_C(P)$ together with jump in the specific heat at critical temperature were calculated. Several other new thorium hydrides were predicted to be stable under pressure including $ThH_3$, $Th_3H_{10}$, $ThH_4$, $ThH_6$. Thorium (which has $s^2d^2$ electronic configuration) forms high-$T_C$ polyhydrides similar to those formed by $s^2d^1$ metals (Y-La-Ac). Thorium - is the next member in Mg-Ca-Sc-Y-La-Ac family of elements forming high-$T_C$ superconducting hydrides.


title: Paraphrase Generation with Deep Reinforcement Learning

abstract: Automatic generation of paraphrases from a given sentence is an important yet challenging task in natural language processing (NLP), and plays a key role in a number of applications such as question answering, search, and dialogue. In this paper, we present a deep reinforcement learning approach to paraphrase generation. Specifically, we propose a new framework for the task, which consists of a \textit{generator} and an \textit{evaluator}, both of which are learned from data. The generator, built as a sequence-to-sequence learning model, can produce paraphrases given a sentence. The evaluator, constructed as a deep matching model, can judge whether two sentences are paraphrases of each other. The generator is first trained by deep learning and then further fine-tuned by reinforcement learning in which the reward is given by the evaluator. For the learning of the evaluator, we propose two methods based on supervised learning and inverse reinforcement learning respectively, depending on the type of available training data. Empirical study shows that the learned evaluator can guide the generator to produce more accurate paraphrases. Experimental results demonstrate the proposed models (the generators) outperform the state-of-the-art methods in paraphrase generation in both automatic evaluation and human evaluation.


title: Stable and unstable vortex knots in a trapped Bose-Einstein condensate

abstract: The dynamics of a quantum vortex torus knot ${\cal T}_{P,Q}$ and similar knots in an atomic Bose-Einstein condensate at zero temperature in the Thomas-Fermi regime has been considered in the hydrodynamic approximation. The condensate has a spatially nonuniform equilibrium density profile $\rho(z,r)$ due to an external axisymmetric potential. It is assumed that $z_*=0$, $r_*=1$ is a maximum point for function $r\rho(z,r)$, with $\delta (r\rho)\approx-(\alpha-\epsilon) z^2/2 -(\alpha+\epsilon) (\delta r)^2/2$ at small $z$ and $\delta r$. Configuration of knot in the cylindrical coordinates is specified by a complex $2\pi P$-periodic function $A(\varphi,t)=Z(\varphi,t)+i [R(\varphi,t)-1]$. In the case $|A|\ll 1$ the system is described by relatively simple approximate equations for re-scaled functions $W_n(\varphi)\propto A(2\pi n+\varphi)$, where $n=0,\dots,P-1$, and $iW_{n,t}=-(W_{n,\varphi\varphi}+\alpha W_n -\epsilon W_n^*)/2-\sum_{j\neq n}1/(W_n^*-W_j^*)$. At $\epsilon=0$, numerical examples of stable solutions as $W_n=\theta_n(\varphi-\gamma t)\exp(-i\omega t)$ with non-trivial topology have been found for $P=3$. Besides that, dynamics of various non-stationary knots with $P=3$ was simulated, and in some cases a tendency towards a finite-time singularity has been detected. For $P=2$ at small $\epsilon\neq 0$, rotating around $z$ axis configurations of the form $(W_0-W_1)\approx B_0\exp(i\zeta)+\epsilon C(B_0,\alpha)\exp(-i\zeta) + \epsilon D(B_0,\alpha)\exp(3i\zeta)$ have been investigated, where $B_0>0$ is an arbitrary constant, $\zeta=k_0\varphi -\Omega_0 t+\zeta_0$, $k_0=Q/2$, $\Omega_0=(k_0^2-\alpha)/2-2/B_0^2$. In the parameter space $(\alpha, B_0)$, wide stability regions for such solutions have been found. In unstable bands, a recurrence of the vortex knot to a weakly excited state has been noted to be possible.


title: Hadronic light-by-light contribution to $(g-2)_\mu$: a dispersive approach

abstract: After a brief introduction on ongoing experimental and theoretical activities on $(g-2)_\mu$, we report on recent progress in approaching the calculation of the hadronic light-by-light contribution with dispersive methods. General properties of the four-point function of the electromagnetic current in QCD, its Lorentz decomposition and dispersive representation are discussed. On this basis a numerical estimate for the pion box contribution and its rescattering corrections is obtained. We conclude with an outlook for this approach to the calculation of hadronic light-by-light.


title: Inapproximability of the independent set polynomial in the complex plane

abstract: We study the complexity of approximating the independent set polynomial $Z_G(\lambda)$ of a graph $G$ with maximum degree $\Delta$ when the activity $\lambda$ is a complex number. This problem is already well understood when $\lambda$ is real using connections to the $\Delta$-regular tree $T$. The key concept in that case is the "occupation ratio" of the tree $T$. This ratio is the contribution to $Z_T(\lambda)$ from independent sets containing the root of the tree, divided by $Z_T(\lambda)$ itself. If $\lambda$ is such that the occupation ratio converges to a limit, as the height of $T$ grows, then there is an FPTAS for approximating $Z_G(\lambda)$ on a graph $G$ with maximum degree $\Delta$. Otherwise, the approximation problem is NP-hard. Unsurprisingly, the case where $\lambda$ is complex is more challenging. Peters and Regts identified the complex values of $\lambda$ for which the occupation ratio of the $\Delta$-regular tree converges. These values carve a cardioid-shaped region $\Lambda_\Delta$ in the complex plane. Motivated by the picture in the real case, they asked whether $\Lambda_\Delta$ marks the true approximability threshold for general complex values $\lambda$. Our main result shows that for every $\lambda$ outside of $\Lambda_\Delta$, the problem of approximating $Z_G(\lambda)$ on graphs $G$ with maximum degree at most $\Delta$ is indeed NP-hard. In fact, when $\lambda$ is outside of $\Lambda_\Delta$ and is not a positive real number, we give the stronger result that approximating $Z_G(\lambda)$ is actually #P-hard. If $\lambda$ is a negative real number outside of $\Lambda_\Delta$, we show that it is #P-hard to even decide whether $Z_G(\lambda)>0$, resolving in the affirmative a conjecture of Harvey, Srivastava and Vondrak. Our proof techniques are based around tools from complex analysis - specifically the study of iterative multivariate rational maps.


title: Center-of-Mass Equations of Motion and Conserved Integrals of Compact Binary Systems at the Fourth Post-Newtonian Order

abstract: The dynamics of compact binary systems at the fourth post-Newtonian (4PN) approximation of general relativity has been recently completed in a self-consistent way. In this paper, we compute the ten Poincar\'e constants of the motion and present the equations of motion in the frame of the center of mass (CM), together with the corresponding CM Lagrangian, conserved energy and conserved angular momentum. Next, we investigate the reduction of the CM dynamics to the case of quasi-circular orbits. The non local (in time) tail effect at the 4PN order is consistently included, as well as the relevant radiation-reaction dissipative contributions to the energy and angular momentum.


title: Improved Approximation Schemes for the Restricted Shortest Path Problem

abstract: The Restricted Shortest Path (RSP) problem, also known as the Delay-Constrained Least-Cost (DCLC) problem, is an NP-hard bicriteria optimization problem on graphs with $n$ vertices and $m$ edges. In a graph where each edge is assigned a cost and a delay, the goal is to find a min-cost path which does not exceed a delay bound. In this paper, we present improved approximation schemes for RSP on several graph classes. For planar graphs, undirected graphs with positive integer resource (= delay) values, and graphs with $m \in \Omega(n \log n)$, we obtain $(1 + \varepsilon)$-approximations in time $O(mn/\varepsilon)$. For general graphs and directed acyclic graphs, we match the results by Xue et al. (2008, [10]) and Ergun et al. (2002, [1]), respectively, but with arguably simpler algorithms.


title: Personalized Schedules for Surveillance of Low Risk Prostate Cancer Patients

abstract: Low risk prostate cancer patients enrolled in active surveillance (AS) programs commonly undergo biopsies on a frequent basis for examination of cancer progression. AS programs employ a fixed schedule of biopsies for all patients. Such fixed and frequent schedules, may schedule unnecessary biopsies for the patients. Since biopsies have an associated risk of complications, patients do not always comply with the schedule, which increases the risk of delayed detection of cancer progression. Motivated by the world's largest AS program, Prostate Cancer Research International Active Surveillance (PRIAS), in this paper we present personalized schedules for biopsies to counter these problems. Using joint models for time to event and longitudinal data, our methods combine information from historical prostate-specific antigen (PSA) levels and repeat biopsy results of a patient, to schedule the next biopsy. We also present methods to compare personalized schedules with existing biopsy schedules.


title: Bounded $H^\infty$-calculus for a Degenerate Elliptic Boundary Value Problem

abstract: On a manifold $X$ with boundary and bounded geometry we consider a strongly elliptic second order operator $A$ together with a degenerate boundary operator $T$ of the form $T=\varphi_0\gamma_0 + \varphi_1\gamma_1$. Here $\gamma_0$ and $\gamma_1$ denote the evaluation of a function and its exterior normal derivative, respectively, at the boundary. We assume that $\varphi_0,\varphi_1\in C^{\infty}_b(\partial X)$, $\varphi_0,\varphi_1\ge 0$, and $\varphi_0+\varphi_1\geq c$, for some $c>0$. We also assume that the highest order coefficients of $A$ belong to $C^\tau(X)$ for some $\tau>0$ and the lower order coefficients are in $L_\infty(X)$. We show that the $L_p(X)$-realization of $A$ which respect to the boundary operator $T$ has a bounded $H^\infty$-calculus.